triggering the collapse of the presolar cloud with an interstellar shock wave propagating away from a site of stellar nucleosynthesis is a favored explanation for the widespread evidence of short - lived radioisotopes ( slri ) in chondritic refractory inclusions ( lee et al .
1976 ; cameron and truran 1977 ; macpherson , davis and zinner 1995 ) and , much more rarely , in chondrules ( russell et al .
1996 ) found in primitive meteorites .
the goal of this paper is to continue the theoretical exploration of the trigggering and injection scenario for slris , in the context of shock waves striking a dense molecular cloud core that could have collapsed to form the solar system .
the dozen or so confirmed ( @xmath4be , @xmath5ca , @xmath6al , @xmath7fe , @xmath8mn , @xmath9pd , @xmath10hf , @xmath11i , and @xmath12pu ) or suspected ( @xmath13tc , @xmath14cl , @xmath15pb , and @xmath16nb ) short - lived radioisotopes may require a fairly involved history for their complete explanation ( goswami and vanhala 2000 ; meyer and clayton 2000 ; mckeegan & davis 2003 ; wadhwa et al .
a stellar nucleosynthetic source ( a supernova or an agb star ) has been the leading explanation for most of these nuclei ( cameron 1993 , 2001 ; wasserburg et al . 1994 , 1995 , 1998 ;
trigo - rodrguez et al .
2009 ; huss et al .
2009 ) , though there are other possibilities as well , in particular local production ( i.e. , in the solar nebula ) of short - lived radioisotopes produced during spallation reactions involving energetic particles emanating from protosolar flares ( shu et al .
such local irradiation models appear to have a problem with being able to match the observed abundance ratio of @xmath6al to @xmath5ca ( srinivasan et al .
1996 ; sahijpal et al . 1998 ; lee et al .
the observed abundance of @xmath6al thus seems to require its production by stellar nucleosynthesis ( mckeegan et al .
however , this problem can be avoided by assuming that the refractory inclusions are shielded by a less refractory mantle during the irradiation , with the mantle being lost later on during heating of the inclusions , thereby yielding the approximate abundance ratio of @xmath6al to @xmath5ca observed in certain meteorites ( gounelle et al .
2001 ) . on the other hand , producing @xmath6al by multiple episodes of local irradiation would negate its use as a precise chronometer for the early solar system ( bizzarro et al . 2004 ; halliday 2004 ; krot et al . 2005 ; thrane et al .
2006 ) , which seems to be ruled out by the agreement of @xmath6al ages with those derived from the pb - pb dating system ( connelly et al .
2008 ) .
evidence has also appeared for the presence of the short - lived isotope @xmath4be in an allende inclusion ( mckeegan et al .
2000 ) . because @xmath4be is thought to be produced only by nuclear spallation reactions , its existence has been used to argue strongly in favour of local irradiation ( mckeegan et al .
2000 ; gounelle et al .
sahijpal & gupta ( 2009 ) have calculated that even if all of the @xmath4be was produced by local irradiation , then the amount of @xmath6al also produced by local irradiation was about 10% of the total amount of @xmath6al , so that the bulk of the @xmath6al was probably synthesized by a massive star . however ,
if irradiation is responsible for the @xmath4be , it is unclear if the irradiation occurred in the solar nebula , or in an earlier phase of evolution .
arnould et al .
( 2000 ) point out that spallation can occur in the winds ejected from h - depleted wolf - rayet ( wr ) stars .
desch et al .
( 2004 ) showed that the @xmath4be might well have originated from @xmath4be galactic cosmic rays that were stopped in the presolar cloud .
evidence has also been advanced for the presence of live @xmath17be in ca , al - rich , refractory inclusions ( cais ) , which are believed to represent the earliest solids formed in the solar nebula that have survived relatively unaltered ( chaussidon et al .
because of the extremely short half - life of @xmath17be of 53 days , this evidence , if correct , would require formation of @xmath17be in the solar nebula .
desch & ouellette ( 2006 ) have disputed the @xmath17be claim , which has not been confirmed by other groups to date .
the short - lived isotope @xmath7fe ( goswami and vanhala 2000 ) can not be produced in the appropriate amounts by spallation , and requires a stellar nucleosynthetic source ( tachibana & huss 2003 ) , as does the bulk of the @xmath6al observed to be polluting the interstellar medium .
the half - life of @xmath7fe is 2.6 million years ( rugel et al .
2009 ) , roughly four times that of @xmath6al , so in any case , the evidence for live short - lived isotopes in refractory inclusions seems to require that no more than about 1 million years elapsed between the nucleosynthesis of some of the short - lived isotopes in a star and the formation of refractory inclusions in the solar nebula .
solar - type stars are believed to form from the collapse of dense molecular cloud cores , which are supported against collapse primarily by magnetic fields ( e.g. , mouschovias , tassis , & kunz 2006 ) , though turbulence also plays a role ( e.g. , kudoh & basu 2008 ) .
collapse of a quiescent cloud core begins once the magnetic field support decreases sufficiently through the process of ambipolar diffusion .
recently kunz & mouschovias ( 2009 ) have shown that ambipolar diffusion leading to collapse and fragmentation is able to reproduce the observed distribution of molecular cloud core masses , i.e. , the initial core mass function , suggesting the importance of magnetic fields for star formation in general .
ambipolar diffusion is estimated to require of order 10 myr to lead to collapse ( mouschovias , tassis , & kunz 2006 ) , a period considerably longer than the half - life of @xmath6al of 0.73 myr .
if the solar system s @xmath6al was produced in a massive star , it may have been injected promptly into the protosolar cloud , which must have then collapsed and formed cm - sized solids , all within @xmath18 myr .
this constraint assumes that the same stellar outflow that carried the @xmath6al may have triggered the collapse of the protosolar cloud and injected other newly - synthesized elements , including other short - lived isotopes ( cameron & truran 1977 ; boss 1995 ) .
abundant observational support exists for the triggering of star formation by expanding supernova shells in upper scorpius ( preibisch & zinnecker 1999 ; preibisch et al .
2002 ) and the cygnus loop ( patnaude et al . 2002 ) , by superbubbles in ob associations ( oey et al . 2005 ; lee & chen 2009 ) , by ionization fronts associated with hii regions ( leppanen , liljestrom , & diamond 1998 ; healey , hester , & claussen 2004 ; hester & desch 2005 ; snider et al .
2009 ) , by generic external shocks ( tachihara et al .
2002 ) , and by protostellar outflows ( barsony et al . 1998 ; sandell & knee 2001 ; yokogawa et al . 2003 ) . here
we consider generic shocks , with a special emphasis on supernovae and agb stars as shock sources .      boss ( 1995 ) showed that shock fronts from a nearby agb star or a relatively distant supernova could trigger the collapse of a 3d dense cloud core and inject shock front material into the collapsing cloud .
foster & boss ( 1996 , 1997 ) studied this process in greater detail for axisymmetric , 2d clouds , and pointed out the crucial role of the assumed isothermal shock front for achieving both goals of triggered collapse and injection .
a supernova shock passes through three phases : ejecta - dominated , sedov blast wave , and radiative ( chevalier 1974 ) .
the latter phase occurs at distances of about 10 pc , after which the shock front sweeps up a cool shell of gas and dust as it propagates .
several recent amr studies ( nakamura et al .
2006 ; melioli et al .
2006 ) have confirmed the results of boss ( 1995 ) and foster & boss ( 1996 , 1997 ) that shock - triggered star formation is likely to occur when the supernova shock front has evolved into a radiative shock , i.e. , the shock front is able to cool so rapidly by radiation that the thin shock front gas is essentially at the same temperature as the ambient gas , which is the same situation as in the isothermal shocks considered by boss ( 1995 ) and foster & boss ( 1996 , 1997 ) .
rayleigh - taylor ( r - t ) fingers were identified as the physical mechanism for achieving injection of dust grains and gas into the collapsing presolar cloud ( foster & boss 1997 ; vanhala & boss 2000 , 2002 ) .
because the r - t fingers strike the outermost layers of the presolar cloud , inducing collapse , the r - t fingers do not reach the central regions until shortly after the central protosun and the early solar nebula have formed , possibly explaining the lack of @xmath6al in certain ( fraction unknown nuclear = fun ) refractory inclusions ( sahijpal & goswami 1998 ) that may have formed before the r - t fingers arrived .
boss ( 2007 ) modeled the r - t finger injection process in the context of this scenario simply by imagining spraying the @xmath6al onto the surface of an existing solar nebula .
a related but alternative scenario involves having a nearby ( @xmath19 0.1 pc ) supernova inject @xmath6al directly into the solar nebula ( rather than into the presolar cloud ) , as studied by ouellette , desch & hester ( 2007 ) .
they found that the gas from a shock front could not be injected efficiently into a protoplanetary disk because of the disk s much higher density compared to the presolar cloud .
chevalier ( 2000 ) had found the same result and attributed it to the hot shock gas not having enough time to cool down . as a result ,
ouellette et al . (
2007 ) suggested that the @xmath6al resided primarily in dust grains that shot through the stalled shock - front gas and thereby penetrated into the disk .
dust grains smaller than 0.1 micron would be deflected , but micron - sized and larger dust grains would be injected with a high efficiency ( ouellette et al .
2009 ) .
supernovae are known to produce large amounts of dust grains , but theoretical models suggest that the newly condensed grains are essentially all smaller than 0.1 micron , and are sputtered to even smaller sizes in the reverse shock front driven into the expanding supernova remnant ( snr ) by the interstellar medium that the snr encounters ( bianchi & schneider 2007 ) .
the models are in accord with dust extinction estimates for an observed reddened qso .
an insignificant fraction of the total dust grain mass is contained in grains larger than 0.1 micron , presenting a problem for scenarios that rely on large dust grains for injection .
however , nittler ( 2007 ) has argued that a sub - class of presolar grains appears to have been formed in a single supernova , conceivably the same supernova that produced many slris .
these presolar grains have sizes of 0.1 to 10 micron , large enough to be injected into the disk , or the presolar cloud core .
what fraction of the mass of the initial dust grain population these relatively large grains represent is difficult to determine , given the processing associated with detecting presolar grains in meterorites and the typical limitation to the study of grains larger than about 0.1 micron ( e.g. , amari , lewis , & anders 1994 ) .
recently it has been suggested that @xmath7fe was injected by a supernova directly into the solar nebula roughly 1 myr after solar system formation ( bizzarro et al .
2007 ) , as in the ouellette et al .
( 2007 ) models . in this scenario ,
the @xmath6al was derived from the wind from a wolf - rayet ( wr ) star , a massive star that would later become a supernova and inject the @xmath7fe .
wr winds do indeed carry large amounts of @xmath6al ( arnould , goriely , & meynet 2006 ) , with winds that are comparable in speed ( @xmath19 1500 km / sec ) to supernova shocks ( marchenko et al .
2006 ) , meaning that such winds would need to be slowed down by sweeping up interstellar gas and dust to speeds less than @xmath20 km / sec if they are to trigger cloud collapse rather than simply shred clouds to pieces ( foster & boss 1996 , 1997 ) .
however , other groups have not been able to replicate the ni isotope data that forms the basis for the @xmath7fe scenario ( dauphas et al . 2008 ; regelous et al .
nevertheless , wr stars should be considered as a possible source of @xmath6al in addition to that obtained from a supernova ( e.g. , gaidos et al .
2009 ) .
williams & gaidos ( 2007 ) estimated that the likelihood of a protoplanetary disk being struck by a supernova shock was less than 1% , but considering that we do not know if the solar system s inventory of short - lived isotopes is rare or not , such an argument can not be considered decisive .
in fact , it has been argued recently that significant @xmath6al is necessary for the development of technological civilizations ( gilmour & middleton 2009 ) .
however , gounelle & meibom ( 2007 ) and krot et al . (
2008 ) argued that injection could not have occurred directly onto a relatively late - phase , low mass solar nebula , as @xmath6al from a massive star supernova would have been accompanied by sufficient oxygen to lead to an oxygen isotope distribution in the solar nebula that would be distinct from that inferred for the sun based on mass - independent fractionation of carbonaceous chondrites and from that recently measured by the genesis mission ( mckeegan et al .
krot et al .
( 2008 ) therefore argued that injection must have instead occurred into the presolar cloud , so that the sun and the solar nebula shared a common reservoir of oxygen isotopes that could then undergo fractionation .
ellinger , young , & desch ( 2009 ) , however , pointed out that supernova explosions are not spherically symmetric , and so some degree of anisotropy in the ejecta is to be expected , perhaps enough to permit @xmath6al injection into the presolar nebula without disturbing the oxygen isotope ratios .
while isothermal shock fronts are capable of simultaneous triggering and injection ( boss 1995 ; foster & boss 1996 , 1997 ; vanhala & boss 2000 , 2002 ) , it has been less clear what happens when detailed heating and cooling processes in the shock front are considered .
vanhala and cameron ( 1998 ) found that when they allowed nonisothermal shocks in their models , they could not find a combination of target cloud and shock wave parameters that permitted both triggered collapse and injection to occur : they could trigger cloud collapse , or they could inject particles , but not both in the same simulation .
such an outcome would be fatal to the triggering and injection hypothesis if definitive .
however , vanhala & cameron s ( 1998 ) models employed a smoothed - particle hydrodynamics ( sph ) code , which has since been shown to be poor at resolving dynamical instabilities such as the rayleigh - taylor or kelvin - helmholtz instabilities ( agertz et al .
2007 ) . furthermore , vanhala and cameron s ( 1998 ) thermodynamical routines led to post - shock thermal profiles that were quite different from those of kaufman & neufeld ( 1996 ) , who found that for a 40 km / sec shock , the post - shock gas cooled down from a maximum of over 3000 k to less than 100 k within a distance of 0.001 pc . with a 25 km / sec shock in the vanhala & cameron ( 1996 ) sph code
, however , the post - shock temperature rose to 3000 k and showed no signs of decreasing over a distance on the order of 0.5 pc .
kaufmann & neufeld ( 1996 ) studied c - type shock fronts , which result in the most successful models of the shock emission from the kleinmann - low nebula in the orion molecular cloud .
kaufman & neufeld ( 1996 ) assumed preshock magnetic field strengths in their c - shock models that are consistent with magnetic field strengths measured by zeeman splitting in molecular clouds ( crutcher 1999 ) .
such nondissociative , magnetohydrodynamic c - type shocks appear to be the correct analogue for the shock speeds ( @xmath19 5 to 40 km s@xmath21 ) that we expect will be necessary to simultaneously achieve triggered collapse and injection . the relatively low postshock temperatures in c - type shocks are crucial for this scenario .
postshock cooling depends sensitively on the detailed microphysics , e.g. , on the emission from rotational states of h@xmath3o , h@xmath3 , co , and oh , and hence on quantities such as the ratio of ortho- to para - hydrogen molecules ( neufeld & kaufmann 1993 ; kaufmann & neufeld 1996 ) .
atomic species are also important coolants , as are dust grains .
more recently morris , desch , & ciesla ( 2009 ) have reexamined the question of cooling shocked gas by h@xmath3o line emission , finding that the cooling rates in the optically thin limit are at least as high as those calculated by neufeld & kaufmann ( 1993 ) .
boss et al .
( 2008 ) found that simultaneous triggering and injection was possible for a 20 km / sec shock striking a 1 @xmath2 cloud , provided that neufeld & kaufmann s ( 1993 ) molecular cooling functions were employed .
gounelle et al .
( 2009 ) interpreted this result as meaning that only a very narrow range of shock speeds was consistent with the triggering and injection hypothesis . in this paper
we examine what happens for a wide range of shock speeds ( and different mass clouds ) compared to the single shock speed considered by boss et al .
( 2008 ) , in order to determine the robustness of the shock wave trigger hypothesis .
achieving adequate spatial resolution of shock - compressed regions of presolar clouds hit with supernova shock fronts , while minimizing the overall computational burden , is an insurmountable problem for a fixed grid hydro code of the type previously used on these problems ( boss 1995 ; foster & boss 1996 , 1997 ; vanhala & boss 2000 , 2002 ) .
clearly the demands of the shock - triggered collapse and injection problem requires adoption of the adaptive mesh refinement ( amr ; e.g. , truelove et al .
1997 ; poludnenko et al .
2002 ) technique , which was designed to handle just this type of situation .
amr techniques automatically insert new grid points in regions of strong physical gradients , and remove them in regions without strong gradients , in order to maximize the spatial resolution in the crucial regions while minimizing the computational burden .
flash employs a block - structured adaptive grid approach using the paramesh package .
advection is handled by the piecewise parabolic method ( ppm ) , which features a riemann solver at cell boundaries that handles shock fronts exceptionally well . in flash
, ppm is incorporated in a form that is second - order accurate in space and time .
we have tested the flash2.5 code s ability to reproduce the results of several different test cases that are relevant to the problem of triggering cloud collapse , namely the sod shock tube problem and the collapse of a pressureless sphere . while the performance of flash on the pressureless sphere collapse is not as accurate as with codes designed to study collapse problems ( e.g. , boss & myhill 1992 ) , flash does a superb job of handling the sod shock tube problem on a cartesian grid , when the shock flows parallel to one axis or at a 45 degree angle .
while the standard flash test cases were run on cartesian grids , we have also reproduced the correct results for the sod shock tube and for pressureless cloud collapse on the cylindrical coordinate ( @xmath22 ) grid that is used in the present calculations .    in the absence of cooling or an isothermal constraint
, flash produces an adiabatic evolution with an effective @xmath23 .
the flash equation of state routines were taken to be those for a simple perfect gas with a mean molecular weight of @xmath24 .
we have adapted these flash routines to simulate isothermal shock - cloud interactions , where the entire computational grid is forced to remain isothermal , typically at 10 k , in order to compare with the results of foster & boss ( 1996 ) . in these models , @xmath25 is set equal to 1.01 , as a value of @xmath26 is prohibited by the riemann solver .
in addition to isothermal models , we also present models that employed the same compressional heating and radiative cooling that was studied by boss et al .
( 2008 ) . our model for radiative cooling is based on the results of neufeld & kaufman ( 1993 ) , who calculated the radiative cooling caused by rotational and vibrational transitions of optically thin , warm molecular gas composed of h@xmath3o , co , and h@xmath3 and found h@xmath3o to be the dominant cooling agent .
neufeld & kaufman s ( 1993 ) figure 3 shows that over the range of temperatures from 100 k to 4000 k , the total cooling rate coefficient @xmath27 can be approximated as @xmath28 erg @xmath29 s@xmath21 . the cooling rate @xmath30 , where @xmath31 is the number density of hydrogen molecules and @xmath32 is the number density of the molecular species under consideration . assuming that @xmath33 ( neufeld & kaufman 1993 ) , we take @xmath34 , leading to a radiative cooling rate of @xmath35 erg @xmath36 s@xmath21 , where @xmath37 is the gas density in g @xmath36 .
boss et al .
( 2008 ) found that @xmath38 could be increased or decreased by factors of two without having a major effect on the outcome of shock triggering and injection , so the precise value of @xmath38 does not appear to be critical to the results .
kaufman & neufeld ( 1996 ) found that the peak temperatures in their mhd shocks were typically of order 1000 k for shock speeds in the range from 5 km / sec to 45 km / sec . as a result , the temperatures in the present models as well as those of boss et al .
( 2008 ) have been restricted to values between 10 k and 1000 k. kaufman & neufeld ( 1996 ) also found a typical shock thickness of order 0.001 pc = @xmath39 cm for a 40 km / sec shock propagating in a magnetized gas with a preshock density similar to that of dense cloud cores .
a similar shock thickness occurs in the models of boss et al .
( 2008 ) , using the neufeld & kaufman ( 1993 ) radiative cooling rate @xmath38 , though with a 20 km / sec shock speed and an unmagnetized cloud : evidently the higher shock speed in kaufman & neufeld ( 1996 ) is roughly compensated for by the presence of magnetic fields , compared to the boss et al .
( 2008 ) results .    as in boss et al .
( 2008 ) , in the present models we used the two dimensional , cylindrical coordinate ( @xmath40 ) version of flash2.5 , with axisymmetry about the rotational axis ( @xmath41 ) .
multipole gravity was used , including legendre polynomials up to @xmath42 .
the cylindrical grid is typically 0.197 pc long in @xmath43 and 0.063 pc wide in @xmath44 , though in the higher shock speed models , the grid was extended to be 0.320 pc long in order to better follow the interaction .
we set the number of blocks in @xmath44 ( @xmath45 ) to be 5 in all cases , while the number of blocks in @xmath43 ( @xmath46 ) has been varied from 5 to 20 .
most models have had 15 blocks in @xmath43 , leading to an approximately uniform grid spacing in @xmath44 and @xmath43 . with each block consisting of @xmath47 grid points ,
this is equivalent to an initial grid of @xmath48 for most models .
the number of levels of refinement ( @xmath49 ) has been varied from five to six . with five levels of refinement employed ,
flash is able to follow small - scale structures with the effective resolution of a grid 16 times finer in scale , or effectively @xmath50 , somewhat better than the highest resolution of @xmath51 used by vanhala & boss ( 2000 ) . with six levels ,
the resolution is increased by another factor of two in each direction .
our target dense cloud cores consist of bonner - ebert ( be ) spheres ( bonnor 1956 ) , which are the equilibrium structures for self - gravitating , isothermal spheres of gas .
be spheres are excellent models for the structure of pre - collapse dense molecular cloud cores seen in star - forming regions ( e.g. , shirley et al . 2005 ) .
for the models that compare results to the isothermal models of foster & boss ( 1996 , 1997 ) , the target dense cloud core is a be sphere with a mass of 1.1 @xmath2 , a radius of 0.058 pc , a temperature of @xmath52 k , and a maximum density of @xmath53 g @xmath36 , located at rest near the top of the cylindrical grid .
the be sphere is embedded in an intercloud medium with a density of @xmath54 g @xmath36 and a temperature of 10 k. the shock wave begins at the top of the grid and propagates downward at a specified speed toward the be sphere .
the shock wave has a thickness of 0.003 pc with a uniform density of @xmath55 g @xmath36 , a mass of 0.016 @xmath2 , and a temperature of 10 k. for the models where cooling is included , the shock wave begins with a temperature of 1000 k and is followed by a post - shock wind with a density of @xmath54 g @xmath36 and temperature of 1000 k , also moving downward at the same speed as the shock wave .
the assumed shock structure is the same as that used in the standard case of foster & boss ( 1996 , 1997 ) and investigated by boss et al .
the shock structure was chosen to resemble the expected conditions in a planetary nebula wind ( e.g. , plait & soker 1990 ; see discussion in foster & boss 1996 ) .
however , it is also consistent with a supernova shock that has swept - up considerable material and slowed down as a result .
chevalier ( 1974 ) considered a supernova shock propagating into a cold ( t = 10 k ) medium with a number density of 1 @xmath36 . by 0.25 myr ,
the shock has slowed to a top speed of @xmath19 60 km / sec and has travelled a distance of @xmath19 2.5 pc .
the amount of swept - up mass contained in the shock front that is incident on the target clouds used in the present calculations is 0.015 @xmath2 , quite close to the value of 0.016 @xmath2 in the standard case .
the shock wave material is represented by a color field , initially defined to be equal to 1 inside the shock wave and 0 elsewhere , which allows the shock wave material to be tracked in time ( foster & boss 1997 ) .
the slri are assumed to be contained primarily in dust grains of sub - micron size ( e.g. , bianchi & schneider 2007 ) , small enough for the grains to remain coupled to the gas .
grains larger than this size could shoot through the shock front as it strikes the target cloud and increase the injection efficiency ( e.g. , ouellette et al .
2007 , 2009 ) , so injection efficiencies derived solely from the color field approach should be considered as lower bounds on the true injection efficiencies .
we present results for several related studies with the flash code , namely a set of comparisons with the results on the standard case of foster & boss ( 1996 , 1997 ) , the stability of be - like spheres with higher central densities ( and higher masses ) in the absence of shock waves , and , finally , the effects of varying the shock wave speed across a wide range of values for a be - like target sphere with a mass of @xmath56 .
we first used flash2.5 to reproduce the standard case of triggered isothermal collapse of foster & boss ( 1996 ) , and verified that flash2.5 was able to produce simultaneous triggered collapse and injection of shock wave material in the case of a 20 km / sec shock .
table 1 lists four models that duplicate the standard case , with varied spatial resolution in flash2.5 .
table 1 also gives the results of the models : the maximum density obtained ( @xmath57 ) in g @xmath36 , the fraction of the incident color field that is injected ( @xmath58 ) , and the final time of the model ( @xmath59 ) in seconds .
once maximum densities of @xmath60 g @xmath36 are reached , even with five or six levels of refinement , flash is unable to follow the collapse to even higher densities , either because it does not have sufficient spatial resolution on the scale of the density maximum , which typically occupies only a few grid cells , or else because the poisson solver is unable to properly represent the gravitational potential of what has become in essence a point mass . as a result , the final times listed in table 1 do not correspond to the time when the density maximum is reached , as flash simply continues to do its best to evolve the entire cloud - shock system past that instant of time .
once maximum densities of @xmath60 g @xmath36 are reached , it is appropriate to terminate these calculations anyway , because by such densities the collapsing cores have become optically thick and can no longer cool by our assumed molecular cooling law , which assumes optically thin clouds . a full radiative transfer treatment ( e.g. , boss & myhill 1992 ) is required to treat optically thick regions , a capability that does not exist in flash2.5 . finally , as in boss et al .
( 2008 ) , @xmath58 is defined to be that fraction of the initial color field that is incident on the initial target cloud and that is injected into regions of the collapsing cloud core with density greater than @xmath61 g @xmath36 .
that is to say , if the total amount of color that is incident on the target cloud is 1 in dimensionless units , @xmath62 means that 0.001 is the amount of color that was injected .
it is clear from table 1 that in all four models , the be sphere was triggered into self - gravitational collapse , given that maximum densities of @xmath60 g @xmath36 or higher were achieved in each case .
this maximum density is well over a factor of @xmath63 times higher than the initial maximum density in the target be sphere of @xmath53 g @xmath36 , indicating that dynamical collapse has been induced in these cloud cores .
it is also clear that as the spatial resolution in flash is increased , the collapsing cloud is able to reach slightly higher densities , implying that the highest resolution calculations are approaching the continuum limit .
note also that the injection efficiency obtained for model fbd , 0.002 , is only slightly lower than the efficiency of 0.003 obtained for the 20 km / sec shock model with heating and cooling presented by boss et al .
( 2008 ) , which had the same spatial resolution as model fbd , showing that when nonisothermal processes are considered , the injection process becomes only somewhat less efficient than when isothermality is assumed .
the evolution of these four isothermal models is very similar to that of the models to be presented below with compressional heating and radiative cooling , and so separate figures for the isothermal models are not displayed .
we have also used flash to verify the long - term stability of the target clouds in the absence of a triggering shock front , as clouds that collapse on their own on time scales similar to the shock wave passage can not be considered to have been triggered into collapse .
shocked - triggered collapse for the standard case occurs within a time span of about @xmath64 yrs ( foster & boss 1996 , 1997 ; boss et al .
when the target be sphere used in the models in table 1 is evolved isothermally in flash2.5 without being struck by a shock wave , it does not collapse , but instead oscillates around its initial equilibrium structure , over a time period of at least @xmath63 yrs , sufficiently long to validate the claim that triggering has occurred in the foster & boss ( 1996 ) comparison models .
in addition , we have tested the stability of higher mass cloud cores , obtained simply by multipying the initial densities of the standard case be sphere by factors of 2 , 3 , or 4 .
formally speaking , these clouds are not true be spheres in equilibrium , but any deleterious effects of this simple computational convenience can be discerned by seeing if the clouds do indeed collapse on their own without being subjected to shock triggering . in the case of the two clouds where
the densities had been increased by factors of 3 or 4 , the clouds did indeed begin to undergo dynamic collapse to densities of @xmath60 g @xmath36 or higher within times of @xmath65 yrs , showing their gravitational instability and unsuitability for the shock - triggering models .
however , the cloud with twice the be sphere density , i.e. , an initial central density of @xmath66 g @xmath36 and a mass of 2.2 @xmath2 , remained stable in flash for at least @xmath67 yrs , proving its suitability for the present studies .
the next section then presents the results for when a 2.2 @xmath2 cloud core is struck by shocks with varying speeds .
we now turn to the main focus of this paper , a consideration of what happens when shock waves with a wide variety of speeds are allowed to strike a dense cloud core , while including a fully nonisothermal treatment with compressional heating and radiative cooling .
table 2 lists the shock speeds employed for these models of a 2.2 @xmath2 target cloud core , as well as the results of the calculations , as in table 1 .
several variations on the spatial resolution were also calculated : models v4 - 4 and v5 - 4 had four levels of adaptive mesh refinement , rather than five , while model v20 - 6 had six .
the higher speed models ( v75l , v80l , v90l , and v100l ) had longer grids in the @xmath43 ( vertical ) direction , in order to better follow the evolution of the cloud as it is swept downward by the shock front .
figure 1 depicts the initial density distribution for all of the nonisothermal models with varied shock speeds , showing the target cloud and the incoming shock front , where the color field is taken to be initially uniform with a value of unity .
figures 1 through 3 show the evolution of model v20 , which typifies the results for models resulting in simultaneous collapse and injection .
model v20 is also identical to model c of boss et al .
( 2008 ) , except for having a cloud mass of 2.2 @xmath2 instead of 1 @xmath2 , and behaves in a very similar manner , implying that the process works equally well for target clouds in this mass range .
figure 2 shows that the model v20 shock front is able to compress the top edge of the target cloud , while rayleigh - taylor ( r - t ) fingers and kelvin - helmholtz ( k - h ) vortices form around the shock - cloud interface .
the r - t fingers drive into the target cloud , while the k - h vortices tend to ablate material off the cloud edge and force it to join the downstream flow . the contours for the color field in figure 2
show that most of the color field impinging on the target cloud is diverted from the cloud by these vortices and disappears downstream .
figure 4 shows a close - up of the shock - cloud interface of model v20 at the same time as figure 2 , providing a better look at the structure of the r - t and k - h features , as well as of the color field .
figure 3 shows that model v20 is able to form a dense clump along the symmetry axis less than 0.1 myr after the evolution began .
figure 5 shows a close - up of the dense clump at the same time as figure 3 , this time with the temperature field contoured in black , instead of the color field .
temperatures rise above 100 k only in regions immediately adjacent to the shock - cloud interface , i.e. , in the region with the largest density gradients and hence the strongest compressional heating . except for this highly shocked region , then , the bulk of the cloud is able to remain nearly isothermal at the assumed background temperature of 10 k as a result of the cooling by molecular species such as h@xmath3o , co , and h@xmath3 .
the strength of this molecular cooling explains why the results of the nonisothermal model presented by boss et al .
( 2008 ) are so similar to those of the isothermal standard test case of foster & boss ( 1996 , 1997 ) .
figure 6 shows an even closer - in view of the density maximum of model v20 after a time of 0.1 myr .
the velocity vectors indicate that the region in the vicinity of the density maximum is trying to collapse onto the clump with speeds as high as several km / sec , which is highly supersonic , considering that the sound speed in 10 k gas is 0.2 km / sec .
evidently , a protostar with a maximum density of @xmath68 g @xmath36 has formed and is growing by the accretion of gas from the target cloud s envelope .
meanwhile , the protostar has been accelerated by the shock front to a speed of order 1 km / sec , and is moving downward as a result .
figure 7 plots the color field over the same region as figure 6 , showing that significant shock front material has been injected into the collapsing protostar and the infalling cloud envelope .
multiple waves of shock front material should be accreted by the protostar , considering the velocity field evident in figure 7 , including the regions with even higher color densities than that within the density maximum , which appear likely to collapse onto the protostar within a few thousand more years .
figure 8 , 9 , and 10 show the evolution of model v4 , where dynamic collapse was triggered , but no significant injection occurred . comparing figure 8 for model v4 with figure 2 for model v20 , at comparable phases of shock - cloud interaction , it is apparent that the much slower speed shock front in model v4 is unable to compress the target cloud s edge to the extent achieved by the shock in model v20 .
in fact , by a time of 0.105 myr ( figure 9 ) , the cloud in model v4 has been clearly triggered into collapse by the shock front , yet the color field lags behind and seems unlikely to achieve a significant injection efficiency .
figure 10 depicts a close - up of the density maximum for model v4 at the same time as in figure 9 , showing that the cloud has been triggered into roughly spherically symmetric collapse , judging from the density distribution and the velocity vectors . at the same time
, the color field within the region plotted in figure 10 is essentially zero , implying that if any shock front material is to be accreted by the growing protostar , it must occur at some later phase in its evolution .
finally , figures 11 , 12 , and 13 show the evolution for model v80l , where injection occurred , but dynamical collapse did not ensue . in this case
, the strong shock front tends to shred the target cloud into streamers ( figure 11 ) , an excellent situation for injecting shock - front material into the same region ( figure 12 ) , but not well - suited for inducing sustained dynamic collapse .
model v80l achieved a maximum density of @xmath69 g @xmath36 by the time shown in figure 11 of 0.056 myr , and by 0.1 myr ( figure 13 ) , the maximum density has dropped to @xmath70 g @xmath36 .
figure 12 also demonstrates that even with the more vigorous shock compressional heating in model v80l , the molecular cooling is able to limit the shock - heated regions to the close vicinity of the shock - cloud boundary .
in summary , figure 14 shows how the critical outcomes of achieving injection and sustained collapse depend on the assumed shock speed , for all the models listed in table 2 .
low speed shocks can induce collapse , but not injection , while high speed shocks result in significant injection , but not in collapse .
shocks falling in the range of about 5 km / sec to 70 km / sec appear to be able to simultaneously induce collapse and achieve injection of significant amounts of shock wave material .
this result is consistent with the prediction by foster & boss ( 1996 ) that shock speeds of @xmath71 km / sec or higher would shred the target clouds and prevent the formation of a collapsing protostar .
table 2 also shows that these basic results are relatively insensitive to the amount of spatial resolution employed , specifically to the number of amr grid levels allowed .
model v5 with the standard 5 levels of amr and model v5 - 4 with 4 levels resulted in quite similar outcomes , as was the case for model v20 with 5 levels and model v20 - 6 with 6 levels .
boss et al .
( 2008 ) found that for a 1.0 @xmath2 target cloud and a 20 km / sec shock , the injection efficiency @xmath58 was 0.003 .
for the comparable 2.2 @xmath2 target cloud in model v20 , @xmath58 was 0.001 , a factor of three times lower for the higher mass cloud .
this difference may be attributed to the fact that while the shock front was identical in both models , the larger mass of the target cloud in model v20 made the task of the shock front more difficult : evidently somewhat lower mass clouds are easier to trigger into collapse and pollute with shock front material than somewhat higher mass clouds , at least to the phase studied by these models ( i.e. , maximum densities less than @xmath60 g @xmath36 ) .
estimates of the injection efficiency @xmath58 have dropped steadily as the spatial resolution and physical modeling have improved .
boss ( 1995 ) found that about half of the impinging shock material ( @xmath72 ) entered the collapsing cloud in his coarsely - gridded , 3d isothermal models .
foster & boss ( 1997 ) found @xmath58 between 0.1 and 0.2 in their relatively coarsely - gridded , 2d isothermal models , while vanhala & boss ( 2000 , 2002 ) found @xmath73 in their increasingly higher spatial resolution , 2d isothermal models .
broadly speaking , a typical value of @xmath74 characterizes all of the successful triggering and injection models listed in table 2 , a value considerably lower than those previously found .
this difference may be attributed to a number of factors , principally the improved spatial resolution of the current models and the inclusion of nonisothermal heating effects , both of which appear to have the effect of reducing @xmath58 compared to lower resolution , isothermal calculations .
the superior shock - handling ability of the ppm hydrodynamics method that flash is based upon undoubtedly also plays a role .
finally , there is the question of how @xmath58 is defined , and when it is evaluated : in vanhala & boss ( 2002 ) , e.g. , @xmath58 was typically evaluated at earlier times than in the present models , and the region over which the color field was considered to have been injected was liberally interpreted to extend quite some distance from the collapsing protostar ( e.g. , figure 3 in vanhala & boss 2002 ) . when the region expected to be accreted by the protostar is defined in a similar manner to that used in the present models , an estimate of @xmath75 results from the vanhala & boss ( 2002 ) models , an estimate more in line with the current values .
nevertheless , such a low injection efficiency may still be in accord with a supernova as the source of the shock wave .
based on the estimates of cameron et al .
( 1995 ) , foster & boss ( 1997 ) noted that the @xmath6al - containing gas and dust in a supernova shock wave would have to be diluted by a factor of @xmath76 in order to explain the inferred initial abundance of @xmath6al in the solar nebula , i.e. , @xmath77 of supernova shock - wave material should be injected into a 1 @xmath2 presolar cloud .
more recently , takigawa et al . (
2008 ) used detailed models of a faint supernova with mixing and fallback to attempt to match the inferred initial abundances of the slris @xmath6al , @xmath5ca , @xmath8mn , and @xmath7fe in the solar nebula , based on nucleosynthetic yield calculations ( e.g. , rauscher et al .
they found that a dilution factor of @xmath78 and a time interval of 1 myr between the supernova explosion and the formation of the first refractory solids in the solar nebula was able to do a good job of matching all four initial abundances .
the dilution factor @xmath79 is the ratio of the amount of mass derived from the supernova that ends up in the solar nebula to the amount of mass in the solar nebula that did not derive from the supernova , i.e. , the mass derived from the target cloud in the present models .
takigawa et al .
( 2008 ) found that the best estimates for @xmath79 depended on the assumed mass of the pre - supernova star , ranging from @xmath80 for a 25 @xmath2 star to @xmath81 for a 20 @xmath2 star ; stars with masses of 30 and 40 @xmath2 led to intermediate values of @xmath79 .
gaidos et al . ( 2009 ) , on the other hand , suggest a value of @xmath79 of at least @xmath82 for a 25 @xmath2 progenitor star , significantly larger than the estimate by takigawa et al .
( 2008 ) .
trigo - rodrguez et al .
( 2009 ) have shown that a 6.5 @xmath2 agb star could have produced the inferred initial abundances of the slris @xmath6al , @xmath5ca , @xmath7fe , and @xmath9pd in the solar nebula , with a similar dilution factor of @xmath83 .
the planetary nebulae formed by agb stars typically have slow wind speeds of 10 km / sec , while the fast winds caused by agb flashes can overtake the slow winds and produce swept - up shells with speeds of 30 km / sec ( frank & mellema 1994 , their figures 1 and 2 ) .
in fact , the swept - up shell in the frank & mellema ( 1994 ) model closely resembles the structure assumed for the shock front here and in the standard cases of foster & boss ( 1996 , 1997 ) : a shock front with a number density of @xmath76 @xmath36 and a thickness of @xmath84 cm , leading to a total shock front mass of 0.016 @xmath2 impacting the target cloud .
based on the present models , then , agb - derived shocks moving at speeds of @xmath19 30 km / sec should be able to trigger collapse and injection in the same manner as a supernova shock with similar properties .    in the present models , the mass of the shock wave that is incident on the target cloud is 0.016 @xmath2 , so values of @xmath74 imply that about @xmath85 of shock front material is injected into the collapsing protostar . if the final result is a well - mixed @xmath86 protostar and protoplanetary disk , then the dilution factor @xmath79 is @xmath87 , considerably lower than required by either takigawa et al .
( 2008 ) or gaidos et al .
in fact , the mismatch is even worse than this , because in the case of a supernova shock front , the material ejected by the supernova must be diluted by the swept - up , intervening interstellar medium that is necessary to slow down the shock front to speeds capable of achieving triggering and injection .
a supernova shock launched at 1000 km / sec must snowplow at least 15 times more mass in order to slow below 70 km / sec .
hence the dilution factor for a supernova shock must be decreased by this same factor , to @xmath88 .
clearly there is a need to learn if these crudely estimated dilution factors can be increased to values closer to those advocated by takigawa et al .
( 2008 ) or gaidos et al .
not all of the 2.2 @xmath2 target cloud will be accreted by the growing protostar , so the injection efficiency will be larger by a proportionate amount .
in addition , if the bulk of the shock front material infalls somewhat later than the earliest arrivals ( as suggested by the color waves in figure 7 ) , then the shock front material may end up preferentially in the protoplanetary disk , rather than in the star , thereby increasing proportionally the dilution factor in the disk . since the present models do not include rotation of the target cloud or shock front
, there is no possibility for a rotationally supported disk to form , and so the present models can not fully answer the question of the dilution factor appropriate for the solar nebula , as opposed to the presolar cloud as a whole .
calculations are currently underway on the dc101 cluster at dtm that include rotation for 2d target clouds , in order to address this key question .
given the lower injection efficiency for a 2.2 @xmath2 cloud compared to a 1.0 @xmath2 cloud , higher cloud densities evidently result in lower injection efficiencies , as might be expected .
target clouds with lower initial densities should then have higher injection efficiencies , and their larger radii ( for a given mass cloud ) will also go in the direction of increasing the total amount of injected shock wave material .
furthermore , allowing injection of shock wave material from behind the leading edge of the shock front ( here considered to be only 0.003 pc thick ) will also increase the amount injected .
increasing the assumed density and thickness of the shock front should also lead to higher injection efficiencies .
future work will include the study of shock fronts with different densities and thicknesses compared to the standard case of foster & boss ( 1996 , 1997 ) and employed in the present models , as well as different density and radii target clouds , in order to better ascertain the suitability of a wider range of possible shock fronts and target clouds for triggering and injection .    coupled with these concerns over the low injection efficiency and the dilution factor
are the implications for the star - formation environment where simultaneous triggering and injection might have occurred .
looney , tobin , & fields ( 2006 ) showed that for the target cloud considered here and an injection efficiency of @xmath89 , due to geometric dilution alone , a supernova would have to occur within about 0.06 pc to 1.2 pc from the target cloud core in order to inject the desired amount of slris .
for the much lower values of @xmath58 found here , @xmath74 , these distance estimates decrease by factors of 10 , to 0.006 pc to 0.12 pc .
such distances are appropriate for the proplyds in the orion nebula that are being photoevaporated by the orion trapezium s four o stars ( e.g. , williams , andrews , & wilner 2005 ) .
however , the proplyds have already collapsed to form protostars , and the density of the surrounding hii region is too low to slow down a supernova shock wave by snowplowing to the required speeds .
thus it remains to be seen if a combination of target cloud and shock front parameters can be found that will increase the injection efficiencies sufficiently to produce a scenario that is consistent with observations of regions of high mass star formation ( e.g. , hester & desch 2005 ) .
when cooling by appropriate molecular species is included , shocks with speeds in the range from 5 km / sec to 70 km / sec are able to trigger the gravitational collapse of otherwise stable , dense cloud cores , as well as to inject shock wave material into the collapsing cloud cores .
this injected material consists of shock wave gas as well as dust grains small enough to remain coupled to the gas , i.e. , sub - micron - sized grains , which are expected to characterize supernova shock waves ( bianchi & schneider 2007 ) and to carry the slri whose decay products have been found in refractory inclusions of chondritic meteorites .
evidently a radiative - phase supernova shock wave ( chevalier 1974 ) is able to cool sufficiently rapidly to behave in much the same way as a shock wave that is assumed to remain isothermal with the target cloud ( e.g. , boss 1995 ) . given that wolf - rayet star winds and supernova shocks both are launched with shock speeds on the order of @xmath90 km / sec , these shock waves can only trigger collapse after they have travelled some distance ( typically about 10 pc ) and been slowed down to 5 km / sec to 70 km / sec by the snowplowing of intervening interstellar cloud gas and dust .
the distance a fast shock must travel in order to slow down to speeds consistent with simultaneous triggering and injection is inversely proportional to the mean density of the intervening material ( assuming this material to be moving much less than the shock speed ) ; typical interstellar medium densities lead to distances of a few pc , depending on the desired shock speed .
an agb star wind with a typical speed of 10 to 30 km / sec could also have triggered collapse and injection without the need for snowplowing .
the low injection efficiencies of the present models , however , point to the need to consider shock fronts and target clouds with different parameters , in order to learn if the injection efficiencies can be increased to the levels thought to be necessary to explain the observed abundances of fossil slris in meteorites .
we are currently running three - dimensional ( 3d ) models on the xenia cluster at dtm .
the need for a 3d treatment of the shock triggering process is evident from previous 3d studies of the r - t instability ( stone & gardiner 2007 ) and shock fronts ( stone & norman 1992 ; whalen & norman 2008 ) .
the r - t `` sheets '' that form in axisymmetric 2d models become true r - t fingers in 3d , allowing better penetration into the target dense cloud cores .
we will present the results of these ongoing 3d models as well as of models with varied shock fronts and target clouds in future papers .
we thank the referee , steven desch , for several constructive reports that have led to a number of significant improvements in this paper .
the calculations were performed primarily on the dc101 cluster at dtm .
this research was supported in part by nasa origins of solar systems grant nng05gi10 g and nasa planetary geology and geophysics grant nnx07ap46 g , and is contributed in part to nasa astrobiology institute grant ncc2 - 1056 .
the software used in this work was in part developed by the doe - supported asc / alliances center for astrophysical thermonuclear flashes at the university of chicago .
v1 & 1.0 & 5 & 15 & 5 & 5.e-12 & 0.0 & 3.8e12 + v2 & 2.0 & 5 & 15 & 5 & 1.e-12 & 0.0 & 8.5e12 + v2.5 & 2.5 & 5 & 15 & 5 & 2.e-12 & 0.0 & 4.8e12 + v4 - 4 & 4.0 & 5 & 15 & 4 & 1.e-12 & 0.0 & 1.2e12 + v4 & 4.0 & 5 & 15 & 5 & 5.e-12 & 0.0 & 3.6e12 + v5 - 4 & 5.0 & 5 & 15 & 4 & 2.e-12 & 2.e-4 & 7.5e12 + v5 & 5.0 & 5 & 15 & 5 & 1.e-12 & 3.e-4 & 8.5e12 + v7 & 7.0 & 5 & 15 & 5 & 2.e-12 & 6.e-4 & 8.5e12 + v9 & 9.0 & 5 & 15 & 5 & 1.e-12 & 2.e-4 & 1.2e12 + v10 & 10.0 & 5 & 15 & 5 & 5.e-12 & 2.e-3 & 6.5e12 + v20 & 20.0 & 5 & 15 & 5 & 1.e-12 & 1.e-3 & 5.3e12 + v20 - 6 & 20.0 & 5 & 15 & 6 & 1.e-12 & 4.e-4 & 5.2e12 + v30 & 30.0 & 5 & 15 & 5 & 4.e-12 & 3.e-3 & 4.0e12 + v40 & 40.0 & 5 & 15 & 5 &
2.e-12 & 1.e-3 & 1.2e13 + v50 & 50.0 & 5 & 15 & 5 & 1.e-12 & 4.e-4 & 4.7e12 + v60 & 60.0 & 5 & 15 & 5 & 1.e-12 & 4.e-4 & 6.5e12 + v70 & 70.0 & 5 & 15 & 5 & 1.e-12 & 3.e-4 & 6.5e12 + v75 & 75.0 & 5 & 20 & 5 & 3.e-13 & 3.e-4 & 8.5e12 + v75l & 75.0 & 5 & 20 & 5 & 1.e-15 & 6.e-4 & 4.5e12 + v80l & 80.0 & 5 & 20 & 5 & 1.e-15 & 5.e-4 & 3.7e12 + v90l & 90.0 & 5 & 20 & 5 & 1.e-15 & 4.e-4 & 3.4e12 + v100l & 100.0 & 5 & 20 & 5 & 1.e-15 & 6.e-4 & 3.0e12 +