the coalescence of binary black holes is one the most promising sources of gravitational waves for interferometric gravitational - wave detectors , such as ligo , virgo and geo600  @xcite .
the first - generation ligo detectors have achieved their design sensitivity and recorded over one year of coincident data  @xcite .
this data , together with data from the virgo detector , are currently being searched for gravitational waves from compact binary coalescence  @xcite .
upgrades to improve the sensitivity of these detectors by a factor of two , and ultimately 10 , are underway .
optimal searches using the enhanced detectors in 2009 will be sensitive to black - hole coalescence out to hundreds of megaparsecs  @xcite .
the advanced detectors , operational next decade , could detect black - hole binaries at distances of over  @xcite .
optimal searches for gravitational waves use matched filtering , which requires accurate knowledge of the waveform  @xcite .
previous searches in ligo data have used post - newtonian and phenomenological templates to search for the coalescence of black - hole binaries  @xcite . over the last several years numerical relativity
has made remarkable breakthroughs in simulating the late inspiral , merger and ringdown of black - hole binaries .
the computational cost of these simulations is high , however , making it impractical to use them directly as template waveforms for use in a matched - filter search .
it has been shown that there is good agreement between the waveforms generated by numerical relativity with analytic post - newtonian waveforms to within just a few orbits of merger  @xcite .
this paper uses the high - accuracy caltech
cornell numerical - relativity waveforms to suggest improvements to the analytic waveforms currently used in gravitational - wave searches by ligo and virgo .
a similar study has been performed by pan  _ et al._using numerical data from pretorius and the goddard groups  @xcite .
our main results are in agreement with their conclusion that a simple extension of the existing stationary - phase approximation to the adiabatic post - newtonian waveforms ( called _ taylorf2 _ in ref .
@xcite ) yields high overlaps with numerical waveforms .    in sec .
[ sec : searches ] , we review the current techniques used for searching for gravitational waves in gravitational - wave detector data .
we discuss the construction of the waveform  a pn
nr hybrid  in sec .
[ sec : pnnrhybridwaveform ] . in sec .
[ sec : efficiency ] we employ the hybrid waveform in a comparison of the detection efficiency of gravitational - wave templates that may be used in upcoming searches of ligo and virgo data . finally , in sec .
[ sec : recommendations ] , we discuss improvements that may be made to the current data - analysis techniques to optimize overlaps .    throughout this paper , we use only the @xmath5 component of the waveform @xmath6 ( as defined , e.g. , in  @xcite ) . for convenience
, we drop the superscript .
whenever possible , we use dimensionless quantities , like @xmath7 , where @xmath8 is the areal radius of the observation sphere , and @xmath9 is the total apparent - horizon mass of the holes in the initial data .
however , for any calculation involving the ligo noise curve , we have a physical scale , and thus use standard mks units .
current searches for gravitational waves from binary black - hole coalescence use matched filtering to search for a waveform buried in noise .
the matched filter is the optimal filter for detecting a signal in stationary gaussian noise .
suppose that @xmath10 is a stationary gaussian noise process with one - sided power spectral density @xmath11 given by @xmath12 . for long integration times ,
the data stream @xmath13 output by the detector will always be dominated by the noise .
thus , we can simply approximate @xmath14 to calculate @xmath15 .    using this power spectral density ( psd )
, we can define the inner product between two real - valued signals  the data stream @xmath16 and the filter template @xmath17by @xmath18 then , given data @xmath16 which may contain either noise @xmath19 or noise and a gravitational wave signal @xmath17 , @xmath20 the matched - filter signal - to - noise ratio ( snr ) is defined as @xmath21    the snr can then be used to construct a detection statistic ( directly or in combination with other statistics ) .
it is therefore important to ensure that the templates used in searches accurately model the expected waveforms to avoid reduction in the value of @xmath22 .
the _ overlap _ between two templates @xmath17 and @xmath23 is defined as @xmath24 the overlap encodes the fractional loss in snr that results from using the template @xmath23 rather than the true waveform @xmath17 . in a search that uses @xmath22 as a detection statistic this corresponds to the fractional loss in distance to which the search is sensitive .
the filter template includes arbitrary time and phase offsets , encoded by the arrival time and phase , @xmath25 and @xmath26 . under a change of these quantities ,
the fourier transform behaves as @xmath27 we maximize over these two variables by calculating the inner product as    @xmath28    note that this integral is just the ( inverse ) fourier transform of the quantity @xmath29 evaluated at @xmath25 . thus finding the maximum over @xmath25 involves taking the fourier
transform and selecting the largest element of the finite set that results from discretization .
in this paper we are concerned with overlaps between pn waveforms and nr signals ; to weight the inner product we use the following psds for initial and advanced ligo : for initial ligo we use an analytic approximation to the ligo design psd given by @xmath30\end{aligned}\ ] ] all integrals start from 40 hz . as shown in fig .
[ fig : stildesandinitialpsd ] , at this frequency the noise is an order of magnitude higher than its lowest value , and below this frequency it rises rapidly as @xmath31 .
the region below 40 hz therefore contributes very little signal power to the snr  @xcite .
the psd for enhanced ligo , which will begin operation in mid 2009 , has a similar shape to that for initial ligo although it has a factor of @xmath32 increase in strain sensitivity .
our results using the initial - ligo psd are therefore valid for enhanced ligo , as the sensitivity factor cancels in eq .
( [ eq : overlapdefinition ] ) ; overlaps depend on the _ shape _ of the psd .    for advanced ligo
we use the gwinc program  @xcite to generate the psd .
bench reports the psd in increments of 0.0124 hz . when calculating discrete integrals against signals sampled at other frequencies we obtain values for the psd by linearly interpolating between the values provided by bench .
we start integrals at 10 hz as that is the point where the noise has increased by two orders of magnitude above its minimum , as also shown in fig .
[ fig : stildesandinitialpsd ] .      searches for gravitational waves in ligo and virgo use a post - newtonian waveform known as__taylorf2 _
_  @xcite .
this is a frequency - domain waveform obtained via the stationary - phase approximation  @xcite , which assumes that the frequency - domain amplitude is simply proportional to @xmath33 ( the lowest - order behavior ) , while its phasing is given by the phase of the time - domain waveform , as a function of frequency . for a binary consisting of masses
@xmath34 and @xmath35 , located at an `` effective '' distance @xmath36 , we have @xmath37{mpc}}{{\ensuremath{d_{\mathrm{eff}}}}}\right)\ ,    { \mathcal{a}}_{\unit[1]{mpc}}(m,\eta)\ , f^{-7/6}\ ,    \e^{i\varphi(f;m,\eta)}\ ,    \label{e : waveform}\ ] ] where @xmath38{mpc}}(m,\eta ) & = &    \left(\frac{5\pi}{24}\right)^{1/2 }    \left(\frac{{g}{\ensuremath{m_{\odot}}}/{c}^2}{\unit[1]{mpc}}\right ) \times \nonumber \\    & \quad & \times \left(\frac{\pi { g}{\ensuremath{m_{\odot}}}}{{c}^3}\right)^{-1/6 }    \left(\frac{\eta}{{\ensuremath{m_{\odot}}}}\right)^{1/2 }    \left(\frac{m}{{\ensuremath{m_{\odot}}}}\right)^{1/3}\ , \end{aligned}\ ] ] and the phasing @xmath39 of the frequency - domain waveform is given to 3.5pn accuracy by the formula  @xcite @xmath40 \left[1 + 3      \ln\left(\frac{v}{v_0}\right ) \right ] \nonumber \\    & & + \left[\frac{11\,583\,231\,236\,531}{4\,694\,215\,680 }      - \frac{640}{3 } \pi^2 - \frac{6\,848}{21 } \gamma \right ] v \nonumber \\    & & + \bigg[\left(-\frac{15\,335\,597\,827}{3\,048\,192 } +      \frac{2\,255}{12 } \pi^2 - \frac{47\,324.0}{63.0 } -      \frac{794\,8}{9}\right ) \eta \nonumber \\    & & \qquad + \frac{76\,055}{1\,728 } \eta^2    - \frac{127\,825}{1\,296 } \eta^3 \bigg ] v \nonumber \\     & & + \pi    \left [ \frac{77\,096\,675}{254\,016 } + \frac{378\,515}{1\,512 }      \eta - \frac{74\,045}{756 } \eta^2 \right ] v^2\bigg]\ .\end{aligned}\ ] ] where @xmath41 , @xmath42 and @xmath43 .
the overall frequency scale is set by the total mass @xmath9 , as can be seen by observing that each occurrence of @xmath44 is accompanied by a factor of @xmath9 .
might be rewritten as @xmath45 .
this term accounts for a time offset altering the phase of the fourier transform . ]
thus , going to a higher - mass system shifts the waveform to lower frequencies . on the other hand ,
to first order , the timescale for the rate of change of the frequency is given by the _ chirp mass _ : @xmath46 clearly , the total mass and chirp mass give us two very different handles on the behavior of the waveform .
these two handles will be important when we try to match the template to our waveform in regions where the post - newtonian and stationary - phase approximations are poor .
this is typically the case for high - mass systems , which only enter the detector band late in the inspiral . in this case
, we can still obtain a high match , at the cost of using templates with the wrong values of @xmath9 and @xmath1 .
we also note that physical binary systems are restricted to @xmath47 .
however , for higher values of @xmath1 , the formulas shown above still give plausible waveforms ; in fact , in some cases these templates match the true waveform better than any physical template .
we will explore the implications of allowing unphysical values for @xmath1 in searches over the templates in sec .
[ sec : unrestrictedeta ] .
note the heaviside function in eq .
( [ e : waveform ] ) .
this contains a cutoff frequency @xmath48 which is used to ensure that the template does not extend to frequencies much greater than the frequencies contained in the expected signal .
this is essentially a third parameter for the template waveform , and will be searched over .
see sec .
[ sec : effectofupperfreqcutoff ] for a discussion of strategies for optimizing detection by changing this cutoff .
frequencies which are often used to characterize coalescing binary black holes are : ( i ) the frequency at the innermost stable circular orbit ( isco ) , @xmath49 , around a single schwarzschild black hole with the total mass of the binary system ; ( ii ) the frequency at the light ring , @xmath50 , around a single schwarzschild black hole with the total mass of the binary system ; ( iii ) the ringdown frequency of the final black hole , which depends on both its spin and mass ; and ( iv ) an `` effective ringdown , '' @xmath51 defined in  @xcite .
current searches use 2pn stationary - phase approximation ( spa ) _ taylorf2
_ templates  @xcite .
it has previously been shown that such waveforms provide acceptable detection templates for binary neutron stars and sub - solar mass black holes  @xcite , but not necessarily for higher - mass black holes .
this is an issue we will investigate below by testing 2pn and 3.5pn templates .
in order to perform our comparison we need to construct a `` true '' black - hole binary waveform , which we might expect to observe with detectors .
a numerical simulation will provide the data for the crucial nonlinear merger phase .
we carefully extract the data and extrapolate it to large radius , and investigate the effects of numerical error on the final result .
because this waveform is very computationally expensive to produce , it covers only about 32 cycles , which is not sufficient for a thorough investigation of the possibility of detecting it in searches of data from gravitational - wave detectors .
thus , we match the numerical waveform to a post - newtonian waveform , producing a hybrid which extends for many thousands of cycles , covering the entire band of interest .
the numerical simulation is the same as that described in refs .
@xcite : an equal - mass , non - spinning , black - hole binary with reduced eccentricity  @xcite , beginning roughly 16 orbits before merger , continuing through merger and ringdown  @xcite .
it is performed with the caltech
cornell pseudospectral code , using boundary conditions designed to prevent constraint violations and gravitational radiation from entering the domain  @xcite .
data is extracted from the simulation in the form of the newman
penrose scalar @xmath52 where @xmath53 and the complex vector @xmath54 are constructed with reference to the coordinate basis . along the positive @xmath55 axis
, we have @xmath56    here , @xmath57 is the timelike unit normal to the spatial hypersurface , and @xmath58 is the unit vector in the positive @xmath55 direction . the vectors @xmath59 and @xmath60 are the standard coordinate vectors , which are not normalized .
@xmath61 is extracted as a function of time , at various radii along the positive @xmath55 axis .
this is then extrapolated to large radii , as described in ref .
@xcite , and in greater detail in ref .
@xcite .
the measured ( instantaneous ) frequency at the beginning of the simulation is @xmath62{hz}\ ,    \frac{{\ensuremath{m_{\odot}}}}{m}\ .\ ] ] the measured ringdown frequency is @xmath63{hz}\ ,    \frac{{\ensuremath{m_{\odot}}}}{m}\ .\ ] ] the measured _ christodoulou _ mass and spin of the final black hole are @xmath64 using this value for the spin , a quasi - analytic formula due to echeverria  @xcite predicts a value of @xmath65{hz}\ , \frac{{\ensuremath{m_{\odot}}}}{m}$ ] , for the ringdown frequency , in close agreement with the measured frequency .      the numerical waveform will be the standard against which we will judge the _ taylorf2 _ waveforms used in ligo data analysis . to understand
how precisely we should trust our final results , we need to understand the accuracy of the waveform itself
. the most obvious measure of the error in this fiducial waveform is its convergence with increasing numerical resolution .
[ f : accuracy ] shows the overlap ( eq .
( [ eq : overlapdefinition ] ) ) between waveforms computed at different resolutions
. the data used here are the extrapolated @xmath66 waveforms , integrated in time twice .        because of the short extent of the numerical waveforms , we need to be careful when using their fourier transforms . the signal can be corrupted easily by the non - periodicity of the waveforms , and the discontinuous jumps that result . for fig .
[ f : accuracy ] we mitigate this problem by increasing the sampling frequency of the input data , and restricting the fourier transform to frequencies corresponding to instantaneous frequencies contained in the data .
the input data can easily be upsampled in the time domain by interpolating the phase and amplitude of the complex data to a finer time grid .
we then perform the transform , and explicitly set the data to zero at frequencies below @xmath67 and above @xmath68 , as given in eqs .
[ eq : measuredstartingfreq ] and  [ eq : measuredringdownfreq ] .
while the results do depend on whether or not we impose these cutoffs , they do not depend sensitively on the actual cutoff frequencies .
the overlap between the lowest- and highest - resolution simulations ( dashed lines ) actually passes through zero , as shown in the upper panel .
presumably , this is because of loss of phase accuracy over the course of the simulation .
all three simulations begin with the same initial data , so the waveforms are most similar at the beginning .
masses for which this is the most important segment ( the lowest masses ) will naturally have the highest overlap between resolutions .
as the simulation progresses , numerical error accumulates  notably in the phase  so the overlap decreases with masses for which later segments dominate the overlap ( higher masses ) .
when the overlap is optimized over arrival time and phase , we can see that the overlap becomes much better , as shown in the lower panel , indicating sufficient accuracy within any frequency band for which phase coherence is required . in either case , the medium and highest resolutions are much more nearly the same . without optimization ,
their overlap is within a few tenths of a percent of 1 ; after optimization , the overlap is within @xmath69 of 1 .
in the rest of our analysis we use the highest - resolution waveform . because we always optimize over arrival time and phase , the lower panel of fig .
[ f : accuracy ] is the most relevant , and shows that the waveform has converged to very high accuracy .
the overlaps we quote below will only be given to three decimal places at most , because this is roughly the accuracy of the single - precision numerical methods used in the rest of the paper .
this accuracy is also sufficient for searches of gravitational - wave data .
thus , the truncation error of the simulated waveform is irrelevant for those purposes .
other sources of error include residual eccentricity and spin , the influence of the outer boundary of the simulation , extrapolation errors , and coordinate effects , as discussed in ref .
the eccentricity had a disproportionately large effect on the error quoted in that paper because of the matching technique , which is not used here . restricting attention to the other effects of eccentricity
, the uncertainty falls below that due to numerical error .
similarly , using the techniques of ref .
@xcite , the initial spins of the black holes have been measured more reliably , and found to be more than an order of magnitude smaller than previously determined , allowing us to reduce the estimate for that error to less than the numerical truncation error .
the various coordinate effects were all estimated to be of roughly the same magnitude as the numerical error .    with the numerical error being many times more accurate than needed for this analysis , and the other sources of uncertainty being of roughly the same size , these considerations indicate that the overall error in our fiducial waveform is substantially less than the precision needed for this analysis
.      numerical simulations can not simulate a very large portion of the inspiral of a black - hole binary system .
indeed , the longest such simulation currently in the literature is the one used here  which extends over just 32 gravitational wave cycles before merger .
fortunately , this is the only stage in which simulations are needed .
it has been shown previously  @xcite that the _ taylort4 _ waveform with 3.5-pn phase and 3.0-pn amplitude matches the early part of this simulation to very high accuracy .
we generate a _
taylort4 _ waveform of over 8000 gravitational wave cycles ( @xmath70 , starting at @xmath71 ) , and transition between the two to create a hybrid .
this long waveform is sufficient to ensure that  even for the lowest - mass systems we will consider  the waveform begins well before it enters the frequency band of interest to ligo .
we begin with @xmath66 data , which will later be integrated to obtain @xmath17 . following ref .
@xcite , we match the numerical waveform to the pn waveform by adjusting the time and phase offsets of the pn waveform to minimize the quantity @xmath72^{2}\ , d t \ .\ ] ] here , we choose @xmath73 and @xmath74 , which is closer to the beginning of the waveform than in the previous paper .
this particular interval is chosen to begin and end at troughs of the small oscillations due to the residual eccentricity @xmath75 in our numerical waveform .
taking a range from trough to trough or peak to peak  rather than node to node , for example  of the eccentricity effects minimizes their influence on the matching .
the eccentricity oscillations can be seen more easily after low - pass filtering the waveform , though we find filtering to be unnecessary for this paper .
the junk radiation apparent in the waveform as shown here has no effect on the resulting match  as we have verified by filtering , and redoing the match . because the final waveform will incorporate no numerical data before @xmath76 and very little immediately thereafter ( as explained below )
, the junk radiation will have no effect on any of our results  as we have also explicitly verified .
in particular , by integrating @xmath61 to obtain @xmath17 , we will effectively smooth the junk radiation .
, that are blended to create the hybrid waveform .
the vertical lines at @xmath77 and @xmath78 denote the region over which matching and hybridization occur .
note that the agreement is well within the numerical accuracy of the simulation , represented by the horizontal bands , throughout the matching region .
also note that the phase difference is fairly flat for a significant period of time after the matching range , which indicates that the match is not sensitive to the particular range chosen for matching . ]    in fig .
[ fig : matchingphasecomparison ] we compare the phase of the numerical and pn waveforms .
the quantities plotted are @xmath79 shown over the interval on which both data sets exist .
the vertical bars denote the matching region .
note that the phase difference is well within the accuracy of the simulation ( about 0.01 radians , represented by the horizontal band ) over a range extending later than the matching region .
also , the difference between the two is fairly flat , which implies that the match is not very sensitive to the region chosen for matching .
because of this , we expect that the phase coherence between the early pn data and the late nr data will be physically accurate to high precision .
the hybrid waveform is then constructed by blending the two matched waveforms together according to @xmath80\ ,    a_{{\mathrm{pn}}}(t)\ , \\
\phi_{{\mathrm{hyb}}}(t ) & = & \tau(t)\ , \phi_{{\mathrm{nr } } } + \left [ 1 - \tau(t )    \right]\ , \phi_{{\mathrm{pn}}}(t)\ .\end{aligned}\ ] ] the blending function @xmath81 is defined by    @xmath82    the values of @xmath76 and @xmath83 are the same as those used for the matching .
the amplitude discrepancy between the pn waveform and the nr waveform over this interval is within numerical uncertainty  roughly @xmath84 . as with the matching technique ( eq .  ( [ eq : matchingchisquared ] ) )
, this method is similar to that of ref .
@xcite , but distinct , in that we blend the phase and amplitude , rather than the real and imaginary parts .
this leads to a smoothly blended waveform , shown in fig  [ fig : waveformsnapshot ] .     of the hybrid waveform used in this analysis : the @xmath85 waveform seen by an observer on the positive @xmath55 axis .
the vertical lines denote the matching and hybridization region .
the @xmath86 on the time axis corresponds to the beginning of data from the numerical simulation . ]    up to this point , the waveform has been @xmath61 data . with the long waveform in hand , we numerically integrate twice to obtain @xmath17 , and set the four integration constants so that the final waveform has zero average and first moment  @xcite .
because of the very long duration of the waveforms , this gives a reasonable result , which is only incorrect at very low frequencies  lower than any frequency of interest to us .
we have also checked that our results do not change when we effectively integrate in the frequency domain by taking @xmath87 which is the frequency - domain analog of the equation @xmath88 .
we now compare the signal described in the previous section to restricted , stationary phase _
taylorf2 _ post - newtonian templates with terms up to order 2.0 , order 3.5 , and a `` pseudo-4.0 pn - order '' term recommended in ref .
overlaps are calculated using the techniques of sec .
[ sec : matchedfiltering ] , with the signal @xmath16 being the hybrid waveform described in sec .
[ sec : pnnrhybridwaveform ] scaled to a range of masses .
we consider both the initial- and advanced - ligo noise curves .
plots of the hybrid waveforms in comparison to the initial - ligo noise curve are shown in fig .
[ fig : stildesandinitialpsd ] .     use the initial - ligo noise curve , with sources at a distance of 100mpc . ]
the masses are chosen so that various frequencies of interest ( the final stitching frequency , the isco , and the ringdown ) occur at the `` seismic wall '' for initial ligo : .
the waveforms @xmath89 are scaled to depict the detectability of the signal , typically quantified by the snr introduced in  ( [ eq : innerproductsnr ] ) , which may be written as @xmath90 in the final expression , the numerator and denominator have the same units , and are directly comparable . because the square root of the denominator is familiar
, we plot that along with the square root of the numerator .
plotting these two quantities together gives a graphical impression of the detectability of the waveform , and the relative importance of each part of the waveform , by its height above the noise curve . in ref .
@xcite , brady and creighton define a slightly different quantity , the characteristic strain @xmath91 the relative factor of @xmath92 they use is present so that they can plot @xmath93 against @xmath94 .
cutler and thorne  @xcite define still another quantity , the signal strength @xmath95 , which is related to the fourier transform by @xmath96 the factor of @xmath97 comes from averaging over the orientation of the binary , which we do not do .
@xmath98 is the ratio of the threshold to the rms noise at the endpoint of signal processing .
, with randomized initial conditions .
note the logarithmic scale on the vertical axis .
the majority of instances produced a lower overlap than the optimum .
we interpret this as pointing to the existence of a broad local maximum which did not coincide with the global maximum . ]    for each template family we initially optimize over signal mass @xmath9 , symmetric mass ratio @xmath99 , and upper cutoff frequency @xmath100 .
the optimization is performed using a nelder
mead ( `` amoeba '' ) algorithm  @xcite .
the amoeba starts with a simplex in the parameter space , and proceeds through a series of steps , each of which will improve the value of the function at at least one vertex .
the algorithm terminates when all vertices have converged to the same point to within a specified tolerance .
this process is deterministic , and amounts to an enhanced steepest - ascent algorithm .
it is therefore only guaranteed to find a local maximum , and indeed we find that an amoeba instance started at a random point in the parameter space is most likely to converge to a point that does not give the highest possible overlap .
we interpret this as being due to a large region in parameter space containing a local maximum and a relatively smaller region containing the global maximum .
we therefore supplement the basic amoeba by running 300 instances with random starting values , and taking the best match obtained over all instances . in repeated
runs the same optimal parameters were found by at least some of the amoebas , which supports the claim that this is the true maximum .
the results of optimizing over all of @xmath101 and @xmath100 for selected masses for initial ligo are given in table  [ tab : threeparamoverlapdetailinitial ] and summarized in fig .
[ fig : threeparamoverlapsummaries ] . for initial ligo , in the range covered by the current compact binary coalescence ( cbc )
low - mass search
@xmath102  @xcite , the pseudo-4.0 pn _ taylorf2 _
waveforms achieve the highest overlaps , exceeding those obtained with 3.5 pn waveforms by @xmath103 . above @xmath104
the 3.5 pn waveforms produce overlaps as much as 4% greater than those obtained with pseudo-4.0 pn waveforms over a range from @xmath105@xmath106 . with the advanced - ligo noise curve , in the cbc low - mass range ,
the 3.5 pn and pseudo-4.0 pn waveforms produce overlaps within 2% of each other , with 3.5 pn producing higher overlaps below 20 @xmath107 and pseudo-4.0 pn producing higher overlaps in the range @xmath108@xmath104 .
pseudo-4.0 pn continues to give the highest overlaps up to @xmath109 , producing overlaps as much as 4% greater than those obtained with 3.5 pn waveforms . above @xmath109 3.5 pn waveforms again
yield the best overlaps , by as much as 6% around 90 @xmath107 .    [ cols="<,^,^,^,^",options="header " , ]     a significant feature of tables  [ tab : threeparamoverlapdetailinitial ] and  [ tab : threeparamoverlapdetail ] is the size of the error bars on the cutoff frequencies . for @xmath110 the cutoff frequency can vary as much as 128% above and 28% below the optimal value while losing no more than 1% of overlap
this leads us to consider the range of possible template parameters which may give high overlaps . in the next section we consider the reduction in overlap as the parameters @xmath48 and @xmath1 are independently varied from the optimal value .     and @xmath1 , which the cutoff frequency @xmath48 is prescribed by the weighted average described below .
the mass ratio @xmath1 is allowed to range over unphysical values .
the best - fit values found for the pseudo-4.0 pn templates are always physical in this case . see sec .
[ sec : unrestrictedeta ] .
right : the same , for the advanced - ligo psd , title="fig : " ]   and @xmath1 , which the cutoff frequency @xmath48 is prescribed by the weighted average described below .
the mass ratio @xmath1 is allowed to range over unphysical values .
the best - fit values found for the pseudo-4.0 pn templates are always physical in this case . see sec .
[ sec : unrestrictedeta ] .
right : the same , for the advanced - ligo psd , title="fig : " ]      as shown in fig .
[ fig : stildesandinitialpsd ] the amplitude of the nr waveforms drops sharply at around the lightring frequency , which depends on the total mass of the binary .
the _ taylorf2 _
waveforms do not model the late inspiral , merger or ringdown and hence will continue to evolve as @xmath33 at all frequencies , increasingly deviating from the nr waveform .
this suggests that the upper frequency cutoff of the _ taylorf2 _
waveform should be chosen to be below the frequency at which the two diverge .
however , the effect of the divergence is mitigated by the psd .
the denominator of the overlap , eq .
( [ eq : overlapdefinition ] ) , depends on @xmath111 which is a constant , and @xmath112 which would increase without limit if not for the psd .
[ fig : filterintegrand ] shows @xmath113the integrand of @xmath112for the initial - ligo noise curve for an example _ taylorf2 _ waveform for an equal - mass @xmath114 binary .
we see that above about 450 hz there is very little contribution to the integrand , and so extending the cutoff frequency above this will not impact the overlap .    ) for a _
taylorf2 _ , 3.5 pn waveform with @xmath115 and @xmath116 , at a distance of 100mpc , using the initial - ligo noise curve .
note that the shape of this curve does not change as we change @xmath9 and @xmath1 ; only the vertical scale changes .
right : overlap between caltech
cornell waveform scaled to @xmath117 and restricted _
taylorf2 _ , 3.5 pn waveform using the best - match values for @xmath9 and @xmath1 , as a function of the cutoff frequency @xmath100 , with the initial - ligo noise curve .
the vertical bars are meant to delineate 1% loss .
note that the upper bound extends to higher frequencies indefinitely .
, title="fig : " ] ) for a _ taylorf2 _ , 3.5 pn waveform with @xmath115 and @xmath116 , at a distance of 100mpc , using the initial - ligo noise curve .
note that the shape of this curve does not change as we change @xmath9 and @xmath1 ; only the vertical scale changes .
right : overlap between caltech
cornell waveform scaled to @xmath117 and restricted _
taylorf2 _ , 3.5 pn waveform using the best - match values for @xmath9 and @xmath1 , as a function of the cutoff frequency @xmath100 , with the initial - ligo noise curve .
the vertical bars are meant to delineate 1% loss .
note that the upper bound extends to higher frequencies indefinitely .
, title="fig : " ]    the numerator of the overlap , @xmath118 , can only increase as the cutoff frequency is raised , however frequencies above the lightring where the waveforms have diverged will contribute very little .
the effect of including higher frequencies on the overlap is therefore determined by the @xmath112 term in the denominator . for systems with ringdown frequencies well above the peak of the integrand in fig .
[ fig : filterintegrand ] , this term will not significantly reduce the overlap .
for example , binaries of total mass roughly @xmath119 have ringdown frequencies at roughly 450hz .
only a small fraction of the snr comes from higher frequencies .
thus , we expect that systems with lower masses should not suffer great loss in overlap if the cutoff frequency is higher than ringdown .
however for higher - mass systems the overlap can be significantly reduced if the upper frequency cutoff is too large .
this is indeed what we find , as shown by a representative example on the right in fig .
[ fig : filterintegrand ] . for this @xmath119 system , using the initial - ligo noise curve , the optimal cutoff frequency is around 450hz  roughly the ringdown frequency . decreasing
the cutoff quickly decreases the overlap .
the cutoff may be increased almost indefinitely , however , with only 0.5% loss in overlap .
this , of course , changes when using the advanced - ligo noise curve .
we revisit this issue in sec .
[ sec : recommendations ] .
the physical symmetric mass ratio is restricted to the range @xmath120 , values above this imply complex - valued masses .
however the pn waveforms are well - behaved for @xmath121 , and as seen from tables  [ tab : threeparamoverlapdetailinitial ] and  [ tab : threeparamoverlapdetail ] , the highest overlaps are often obtained at unphysical values of @xmath1 . in fig .
[ fig : physicaleta ] we show the effect of limiting the optimization to physical @xmath1 . at high masses
, the limitation reduces the optimal overlap by up to 12% .
_ taylorf2 _ waveforms with @xmath122 would not be expected to accurately model the late - inspiral and merger part of the waveform , as non - newtonian effects are increasingly significant in this region .
we find that allowing unphysical @xmath1 broadens the space of waveforms covered by the _
approximation sufficiently to capture more of the late - inspiral and merger .     to range over unphysical values , compared to those obtained by restricting the range of @xmath1 .
these overlaps are generated using 3.5 pn taylorf2 templates , searching over values of the total mass and mass ratio .
extending to unphysical values of @xmath1 improves the match by up to 11% . ]
based on the analysis of the previous sections we propose a series of adjustments to searches using _
taylorf2 _ template waveforms to enhance the efficiency of those searches .
first , as seen in fig .
[ fig : threeparamoverlapsummaries ] for initial ligo , adding terms up to 3.5 pn order produces overlaps as large or larger than the current 2.0 pn templates over most of the mass range , while the pseudo-4.0 pn templates recommended in ref .
@xcite produce slightly larger overlaps at masses near @xmath123 .
thus , we recommend pseudo-4.0 pn templates for the low mass range , @xmath124 , and 3.5 pn templates for higher masses .
the improvement due to 3.5 pn templates over 2.0 pn generally holds for advanced ligo as well .
the 3.5 pn templates produce larger overlaps than 2.0 pn templates above @xmath125 without a significant loss ( within 1% ) at lower masses .
however , there is a large region for which the pseudo-4.0 pn term does significantly better .
when using an advanced - ligo noise curve , we recommend 3.5 pn templates generally , 2.0 pn templates in the range @xmath126@xmath127 and pseudo-4.0 pn templates for masses in the range @xmath128@xmath4 .    as a second improvement
, we note from fig .
[ fig : physicaleta ] that allowing @xmath1 to range over unphysical values significantly improves matches with 3.5 pn templates above @xmath129 . in preliminary studies
we have found that extending to @xmath130 roughly doubles the size of the template bank , and the advantages must therefore be weighed against the increase in false alarm rate .
values for 3.5 pn templates with initial ligo .
the dark gray band contains cutoff frequencies with matches within 1% of the value at which the best overlap was obtained .
the light gray band contains frequencies with matches within 3% .
right : candidate @xmath100 values for 3.5 pn templates with advanced ligo .
the dark gray band contains cutoff frequencies with matches within 1% of the value at which the best overlap was obtained .
the light gray band contains frequencies with matches within 3% .
note that the weighted - average cutoff extends past the 1% error bars for @xmath131 .
however , in that same region , the 3.5 pn templates do poorly overall , and we recommend pseudo-4.0 pn templates .
the optimal cutoff frequency for pseudo-4.0 pn templates is much closer to the weighted - average cutoff in this mass range.,title="fig : " ]   values for 3.5 pn templates with initial ligo .
the dark gray band contains cutoff frequencies with matches within 1% of the value at which the best overlap was obtained .
the light gray band contains frequencies with matches within 3% .
right : candidate @xmath100 values for 3.5 pn templates with advanced ligo .
the dark gray band contains cutoff frequencies with matches within 1% of the value at which the best overlap was obtained .
the light gray band contains frequencies with matches within 3% .
note that the weighted - average cutoff extends past the 1% error bars for @xmath131 .
however , in that same region , the 3.5 pn templates do poorly overall , and we recommend pseudo-4.0 pn templates .
the optimal cutoff frequency for pseudo-4.0 pn templates is much closer to the weighted - average cutoff in this mass range.,title="fig : " ]    our third recommendation involves the cutoff frequency used for the template waveform .
optimization over the cutoff frequency is too computationally intensive to be done in searches .
currently , the cutoff frequency is typically taken to be the schwarzschild isco frequency . to examine the effect of this choice we vary @xmath100 while keeping the mass and @xmath1 at their optimal values , for each of the signal masses in our range .
the result of one such variation is shown in fig .
[ fig : filterintegrand ] ( right ) . figs .
[ fig : fcrecomendations ] shows the variations for all masses , highlighting the regions within which the overlap drops by less than 1% ( dark gray ) and 3% ( light gray ) of the optimal value .
this figure also shows the isco and erd frequencies , neither of which stays within the 1% band for both initial and advanced ligo . in particular
, the isco is a poor choice for both initial and advanced ligo except at very low masses , where the precise value of the cutoff is almost irrelevant .
the isco is often pointed to  somewhat arbitrarily  as a good estimate of the breakdown of post - newtonian approximations  @xcite .
so , for instance , if we were to match a pn template to a physical waveform , beginning at some point in the distant past , we might expect them to separate quite badly near the isco .
of course , for realistic black - hole binaries , the gravitational waves will only enter the ligo band late in the inspiral  just before the isco for low - mass systems , or after the isco for high - mass systems .
we can see from fig .
[ fig : stildesandinitialpsd ] that , for masses below about @xmath129 , the isco is high enough that lower - frequency parts of the waveform contribute the most to the snr . for very high masses
, however , this basically cuts the waveform down to nothing . in initial ligo
, the isco is completely buried in seismic noise for masses above about @xmath132 .
thus , we must move the cutoff frequency up .
we can not push the cutoff far above ringdown , because the physical waveform simply ceases to exist ( see fig .
[ fig : stildesandinitialpsd ] ) .
it has been suggested that an `` effective ringdown '' ( erd ) frequency @xmath51 is a useful upper limit  @xcite . for intermediate masses
, we would like to interpolate somehow between these two extremes of isco and erd .
we suggest setting the cutoff frequency to a weighted average of the two , where the weights are the contributions to the snr below the given frequency .
if we assume coherent phasing between the template and the physical waveform , we can simply take the amplitudes of the two waveforms .
also , note that the restricted spa approximation for the amplitude is reasonable .
thus , define @xmath133 we have already dropped constant factors in the expressions for @xmath22 that will cancel out .
note that these expressions only depend on the total mass by way of the limits of integrations  which are very simple , known functions of the mass  so these integrals could be done just once for a given noise curve , storing the intermediate values .
when the cutoff needs to be calculated , the cumulative integral could be evaluated at the given isco and ringdown frequencies .
hence , this would be a fast way of calculating the cutoff , with no need to do the integrals each time the cutoff is needed .
we can test this recommended frequency by comparing it to the optimal cutoff frequency found by the amoeba search described in sec .
[ sec : efficiency ] . for 3.5 pn templates in initial ligo
, we find that it is an excellent match to the optimal frequency .
[ fig : fcrecomendations ] shows these two values , along with dark and light bands showing the regions in which changing @xmath48 results in a loss of overlap of 1% and 3% , respectively .
of course , the same figure shows that using the erd recommendation would stay within the 1% error bounds .
nonetheless , the close match between this recommendation and the true optimum suggests that it is sound .
thus , our final recommendation is to use the weighted - average frequency cutoff throughout the entire mass range .
while our analysis has been restricted to equal - mass systems , the cutoff frequency we have defined here could be applied to unequal - mass systems as well .
it will be interesting to see how this cutoff fares in those situations .
similar results hold for advanced ligo , when using our recommended template for each mass .
that is , in regions where 3.5 pn templates do poorly ( see fig .  [
fig : threeparamoverlapsummaries ] ) , the weighted average is a poor predictor of the optimal cutoff frequency using those templates , as shown in fig .  [
fig : fcrecomendations ] .
however , in those same regions  where pseudo-4.0 pn templates do
well  the weighted average is a good predictor of the optimal cutoff frequency for 4.0 pn templates .
thus , again , we recommend using the weighted - average frequency cutoff throughout the entire mass range with advanced ligo .    by prescribing a cutoff frequency , the search does not need to extend over that parameter .
similarly , by prescribing a post - newtonian order , we need use only one template for a given total mass . on the other hand , if these recommendations decrease the overlap found by too much when using them compared to the overlap found by an unconstrained search , it may be better to search the larger parameter space .
we can evaluate the loss in overlap by comparing the results found using our recommendations to the results found when searching over the set of all three template families , and all masses , mass ratios , and cutoff frequencies .
we have determined that this loss in overlap when using our recommendations is always less than 0.0025 for initial ligo , and less than 0.007 for advanced ligo .
we have compared high - accuracy nr waveforms for equal - mass binary black holes from the caltech  cornell group to stationary phase post - newtonian waveforms .
we examined a number of factors that influence the matches between the two , with the goal of optimizing the matches and hence improving the efficiency of templated searches in initial and advanced ligo .
we first considered the effect of the post - newtonian order to which the phase evolution is taken , and found that adding terms up to 3.5 pn or pseudo-4.0 pn to the currently - used 2.0 pn templates significantly improves the matches over a large range of masses , as shown in fig .
[ fig : threeparamoverlapsummaries ] .
we then studied the effect of varying the upper cutoff frequency of the templates .
the frequency that achieves the optimal match is a function of mass , and we find this function is well - approximated by an average between isco and erd , weighted by contribution to the snr , as shown in fig .
[ fig : fcrecomendations ] . finally , we allow the symmetric mass ratio @xmath1 to range over unphysical values up to @xmath134 , and find that this dramatically improves matches , as shown in fig .
[ fig : physicaleta ] . based on the results we recommend adjusting the searches using
_ taylorf2 _ template waveforms by going up to 3.5 pn or 4.0 pn over most of the mass range , integrating up to our recommended cutoff , and allowing allowing @xmath1 to extend up to 1 . for initial ligo , the overlaps obtained using these parameters is always within 0.0025 of overlaps achievable by optimizing over all three parameters .    in future work
we plan to extend this analysis to unequal - mass and spinning black - hole systems .
we have found that allowing unphysical values of @xmath1 roughly doubles the size of the template bank , and we also plan to study the impact of this on the false alarm rate .
we thank luisa buchman , yanbei chen , curt cutler , lisa goggin , larry kidder , luis lehner , lee lindblom , harald pfeiffer , mark scheel and kip thorne for useful discussions .
this work was supported in part by grants from the sherman fairchild foundation and by nsf grants phy-0601459 , phy-0652995 , and dms-0553302 to caltech .
10 url # 1#1urlprefix[2][]#2 thorne k  s 1987 _ three hundred years of gravitation _ ed hawking s  w and israel w ( cambridge : cambridge university press )                                                                      ajith p , babak s , chen y , hewitson m , krishnan b , sintes a  m , whelan j  t , brgmann b , diener p , dorband n , gonzalez j , hannam m , husa s , pollney d , rezzolla l , santamara l , sperhake u and thornburg j 2008 _ phys . rev .
d _ * 77 * 104017


the phenomenon of shear localization during high strain - rate deformations of metals @xcite is a striking instance of material instability in mechanics , that has attracted considerable attention in the mechanics literature ( _ e.g. _ @xcite ) as well as a number of mathematical studies ( @xcite and references therein ) .
material instability is typically associated with ill - posedness of an underlying initial value problem , what has coined the term _ hadamard instability _ for its description in the mechanics literature .
it should however be emphasized that while hadamard instability indicates the catastrophic growth of oscillations around a mean state , what is observed in localization is the orderly albeit extremely fast development of coherent structures : the shear bands .
the latter is a nonlinear phenomenon and a linearized analysis can only serve as a precursor to such behavior .
we work with the simplest model capturing the mechanism of shear band formation , a simple shear motion of a non - newtonian fluid with temperature decreasing viscosity @xmath0 , @xmath1 where @xmath2 describes the strain - rate sensitivity which typically in these problems satisfies @xmath3 . mathematical studies of , for power law viscosities @xmath4 , @xmath5 , show that the uniform shear is asymptotically stable for @xmath6 @xcite , but indicate unstable response in the complementary region @xmath7 , @xcite .
the mathematical analysis has not so far characterized the precise behavior in the instability region ; however , detailed and resolved numerical studies ( _ e.g. _ @xcite ) show that instability is followed by formation of shear bands .
the difficulty is that in order to analyze the unstable regime one needs to account for the combined effect of parabolic regularizations in conjunction to ill - posed initial value problems .
an asymptotic theory for calculating an effective equation for the compound effect is proposed in @xcite .
it is based on the theory of relaxation systems and the chapman - enskog expansion and produces an effective asymptotic equation . for a power law the effective equation changes type from forward to backward parabolic , across the threshold @xmath8 of the parameter space .
the criteria of @xcite are based on asymptotic analysis and as such are formal in nature .
our objective , here , is to develop a mathematical theory for the onset of localisation and the emergence of a nonlinear localized state , thus providing a rigorous justification to the behavior conjectured by the effective equation in @xcite .
we employ the constitutive hypothesis of an _ exponential law _ for the viscosity coefficient , @xmath9 the model for the exponential law admits a special class of solutions describing uniform shearing @xmath10 the graph @xmath11 describing the average stress vs. average strain response is monotonically decreasing , and thus the deformation of uniform shear is characterized by strain softening response .
we set the following goals :    1 .   to study the linearized problem around the uniform shearing solution and to assess the effects of the various parameters of the problem 2 .   to inquire whether there are special solutions of the nonlinear problem that indicate the emergence of localized deformations and formation of shear bands .
since the base solution around which we linearize is time - dependent the linearized analysis necessarily involves analysis of non - autonomous systems and also perturbations have to be compared against the base solution .
we perform a linearization using ideas of relative perturbations @xmath12 leading to linearized system @xmath13 with @xmath14 computed by @xmath15 .
this problem turns out to have the following behavior :    1 .
[ en : introa]for the case without thermal diffusion @xmath16 : + when @xmath17 the high frequency modes grow exponentially fast with the exponent of the order of the frequency what is characteristic of ill - posedness and _ hadamard instability_. + when @xmath2 the modes are still unstable but the rate of growth is bounded independently of the frequency .
in fact , the behavior in that range is characteristic of _ turing instability_. 2 .
the effect of thermal diffusion ( @xmath18 ) for the non - autonomous linearized model can also be assessed : perturbations do grow initially in the early stages of deformation , but over time the effect of the diffusion intensifies and eventually the system gets stabilized . for a quantitative description of the above qualitative statement
we refer to lemmas [ lemdecay ] and [ lembound ] .
in addition , this behavior can also occur for the full solution of the nonlinear problem subject to simple shear boundary conditions as indicated in numerical simulations presented in figure [ mtstb ] .
this type of metastable response was already conjectured ( for a somewhat different model ) in numerical experiments of @xcite .
the linearized analysis presented here offers a theoretical explanation of this interesting behavior .    in the second part of this article we consider the adiabatic variant of the model with exponential constitutive law .
the system is expressed for @xmath19 in the form @xmath20 and is now considered on the whole space @xmath21 . for the system , we construct a class of localizing solutions @xmath22 \sigma ( x , t ) & = \sigma_s ( t ) \frac{1}{\phi_\lambda ( t ) }   \ ,   \sigma_\lambda \big ( \sqrt{\lambda } x \phi_\lambda ( t ) \big ) , \\[2pt ] \theta(x , t ) & = \theta_s(t ) + \lambda \tfrac{n+1}{\alpha } ( \theta_s(t ) - \theta_0 ) + \theta_\lambda \big ( \sqrt{\lambda } x \phi_\lambda ( t ) \big ) , \end{aligned}\ ] ] where @xmath23 depends on the parameter @xmath24 , and is determined via the solution @xmath25 of the singular initial value problem , with @xmath26 .
the existence of the solution @xmath25 is established in theorem [ connection ] where most importantly its behavior at the origin and at infinity is made precise .
the emerging solution depends on three parameters : @xmath27 linked to the uniform shear in , a parameter @xmath28 which can be viewed as a length scale of initial data , and a parameter @xmath29 determining the size of the initial profile .    to our knowledge
this is the first instance that the compound effect of hadamard instability and parabolic regularizations is mathematically analyzed in a nonlinear context .
the reader should contrast the nonlinear response captured by the exact solution to the response of the associated linearized problem .
the solutions indicate that the competition among hadamard instability and momentum diffusion in a nonlinear context can lead to a nonlinear coherent structure that localizes in a narrow band ( see theorem [ shearband ] and remark [ rmkbeh ] ) .
by contrast , for a linear model , the combined effect of hadamard instability and parabolicity leads to unstable oscillatory modes growing in amplitude ( see ) .
this justifies the scenario that nonlinearity can kill the oscillations caused by hadamard instability so that the non - uniformity of deformation produced by the instability merges into a single localized zone
. it would be interesting to explore this type of response for other models related to material instability in solid mechanics .
the exposition is organized as follows : section [ secdescription ] presents a description of the problem and the formulation via relative perturbations .
the linearized analysis is presented in section [ stabar ] , first the part that can be done via spectral analysis for autonomous systems , and then the part concerning the effect of thermal diffusion , which is based on energy estimates for a non - autonomous linear system .
detailed asymptotics for the eigenvalues and the connection with hadamard and turing instabilities are presented there .
the presentation of the localized solutions is split in sections [ seclocal]-[seclocalii ] : section [ seclocal ] presents the invariance properties , the _ ansatz _ of localized solutions , and also introduces the relevant auxiliary problem , .
the construction of solutions to the singular initial - value problem , and the proof of their properties is the centerpiece of our analysis and forms the objective of section [ secss ] .
the construction is based on a series of steps , involving desingularization , a nonlinear transformation , and eventually the existence of a heteroclinic orbit via the poincar - bendixson theorem .
the construction of the localized solutions is then effected in section [ seclocalii ] .
in this work we study the onset of localization and formation of shear bands for the model describing shear motions of a non - newtonian fluid with viscosity @xmath30 decreasing with temperature . for concreteness
we take a viscosity satisfying the exponential law and are interested in the range where the rate sensitivity @xmath31 is positive and small .
the model then becomes @xmath32 and describes shear motions between two parallel plates located at @xmath33 and @xmath34 with @xmath35 the velocity in the shear direction , @xmath36 the shear stress , @xmath37 the temperature , and @xmath38 depicting the strain rate .
the boundary conditions are taken @xmath39 equation reflects prescribed velocities at the boundaries , while manifests that the bounding plates are adiabatic .
note that implies that the stress satisfies @xmath40 and that @xmath41 the initial conditions are @xmath42 throughout the study we will consider two cases : when @xmath43 the process is adiabatic and the resulting model will be the main object of study , as it is the simplest model proposed for studying the phenomenon of formation of shear bands .
the case @xmath18 will also be studied in section [ stabar ] as a paradigm to assess the effect of thermal diffusion .
the problem , , admits a class of special solutions describing uniform shear @xmath44 for the special case of the exponential law , @xmath45 and @xmath46 read @xmath47 these are universal solutions , _
i.e. _ they hold for all values of the parameters @xmath31 , @xmath48 , @xmath49 in the model including the limiting elliptic initial value problem when @xmath50 , @xmath43 . the graph @xmath51 is viewed as describing the stress vs. average strain response . as the graph is decreasing ,
it means that there is always softening response in the course of the deformation .
the system has served as a paradigm for the mathematical study of the phenomenon of shear bands ( _ e.g. _ @xcite ) occurring during the high strain - rate deformations of metals .
formation of shear bands is a well known material instability that usually leads to rupture , and has received considerable attention in the mechanics literature ( _ e.g. _ @xcite ) .
the reader is referred to @xcite for a description of the problem intended for a mathematically oriented reader that is briefly also outlined below .
it was recognized by zener and hollomon @xcite ( see also @xcite ) that the effect of the deformation speed is twofold : ( a ) to change the deformation conditions from isothermal to nearly adiabatic ; under nearly adiabatic conditions the combined effect of thermal softening and strain hardening tends to produce net softening response .
( b ) strain rate has an effect _ per se _ and induces momentum diffusion that needs to be included in the constitutive modeling . a theory of thermoviscoplasticity is commonly used to model shear bands and we refer to @xcite for various accounts of the modeling aspects .
the basic mechanism for localization is encoded in the behavior of the following simple system of differential equations @xmath52 ( which is the adiabatic case @xmath43 of ) .
this model is appropriate for a non - newtonian fluid with temperature dependent viscosity , and is also a simplification of the models proposed in @xcite to capture the mechanism for shear band formation in high strain - rate plastic deformations of metals .
the argument - adapted to the language of a temperature dependent fluid - goes as follows : as the fluid gets sheared the dissipation produces heat and elevates the temperature of the fluid
. the viscosity decreases with the elevation of the temperature and the fluid becomes softer and easier to shear .
one eventuality is that the shear deformation distributes uniformly across the material in the same manner as it does for the usual newtonian fluids .
this arrangement is described by the uniform shearing solution .
a second mode of response is suggested when comparing the behavior of the fluid in two spots : one hot and one cold .
since the viscosity is weaker in the hot spot and stronger in the cold spots , the amount of shear generated in the hot spot will be larger than the shear generated in the cold spot and in turn the temperature difference will be intensified .
it is then conceivable that the deformation localizes into a narrow band where the material is considerably hotter ( and weaker ) than the surrounding environment .
the two possibilities are depicted in fig .
[ shearflow ] .
the outcome of the competition of thermal softening and momentum diffusion is captured by the behavior of solutions to .
a series of works analyze the system mostly for power law viscosities @xmath53 these studies show stable response @xcite when @xmath54 , and indicate unstable response @xcite but without characterizing the behavior in the complementary region @xmath55 .
more precisely , in the unstable regime tzavaras @xcite shows that development of shear bands can be induced by stress boundary conditions leading to finite - time blow - up , and that their occurrence is associated with a collapse of stress diffusion across the band . on the other hand , bertsch , peletier and verduyn - lunel
@xcite show that if the problem has global solutions and there is a sufficiently large initial temperature perturbation , then a shear band appears as the asymptotic in time state . careful and detailed numerical investigations @xcite
indicate that instability is widespread in the parameter regime @xmath55 and is typically followed by the development of shear bands .
the behavior in the  instability domain " is at present poorly understood , especially regarding the precise mechanism of formation of the shear band .
the reason is the response in this regime results from the competition between an ill - posed problem and a parabolic regularization , and little mathematical theory is available to study such effects .
nevertheless , an asymptotic theory similar to the chapman - enskog expansion is used by katsaounis - tzavaras @xcite to calculate an effective equation for the compound effect of thermal softening and strain - rate sensitivity .
the effective equation changes type from forward to backward parabolic across the instability threshold , and thus provides an asymptotic criterion for the onset of localization . in the present work
, we will develop a mathematical framework for studying the competition and for examining the conjecture of @xcite .
our analysis is facilitated by the use of the exponential law inducing special scaling properties on the resulting system of partial differential equations .      in the mechanics literature shear localization
is typically associated with hadamard - instability .
this is easily put on display for , by considering the limiting case @xmath50 and @xmath43 . in the parlance of mechanics
, this corresponds to rate - insensitive deformations and absence of heat conduction , and leads to the system of conservation laws @xmath56 for @xmath57 the wave speeds are both imaginary , @xmath58 , and the linearized initial value problem for is ill - posed .
two questions then emerge :    1 .   to assess the combined effect of the competition between hadamard instability and momentum diffusion and/or heat conduction ( as manifested by small @xmath2 and @xmath18 respectively ) .
while hadamard instability indicates the catastrophic growth of oscillations around a mean state , it does not by itself explain the formation of coherent structures typically observed in localization .
the latter is an outcome of the nonlinear effects and has to be assessed at the level of the nonlinear problem .
the first objective is to develop quantitative criteria for the stability of the uniform shearing solutions .
since the basic solutions are time dependent an issue arises how to define their stability .
following @xcite , we call the uniform shearing solution _ asymptotically stable _ if perturbations of the solution decay faster than the basic solution , and we call it _ unstable _ if they grow faster than the uniform shearing solution
.    it is expedient to introduce a relative perturbation formulation of the problem : for the _ exponential law _ , we introduce the state @xmath59 describing the relative perturbation and defined through the transformation @xmath60 where @xmath61 is given by and satisfies in particular that @xmath62 using and we derive the system satisfied by the relative perturbation @xmath59 in the form @xmath63 \theta_t - \kappa \theta_{xx } & = \sigma_s ( t )    \ , \big ( \sigma \ , u -1   \big ) , \\[2pt ] \sigma & = e^{-\alpha \theta } u^n .
\end{aligned}\ ] ] note that the uniform shearing solution is mapped to the state @xmath64 , @xmath65 and @xmath66 , and the latter is an equilibrium for the system .
the problem then becomes to study the stability of the equilibrium @xmath67 for the non - autonomous system .
in this section we present an analysis of the linearized system . since the uniform shearing solution is time - dependent the linearized analysis will seek to classify the precise growth in - time of perturbations .
we compute the linearized problem , using the ansatz @xmath68 where @xmath69 is a small parameter measuring the size of the perturbation and @xmath70 .
expansion of the equation @xmath71 gives to the leading order of @xmath69 @xmath72 while linearizing the other two equations gives @xmath73 collecting the equations together and neglecting the terms of @xmath74 we obtain the linearized system : @xmath75 { \bar \theta}_t - \kappa { \bar \theta}_{xx } & = \sigma_s   ( t ) \ ,   \big ( ( n+1 ) { \bar u}- \alpha { \bar \theta}\big ) , \\[2pt ] { \bar \sigma}&= n { \bar u}- \alpha { \bar \theta}. \end{aligned } \right .\ ] ] the boundary conditions and and the transformations and , imply that the induced boundary conditions for the linearized system are @xmath76    in view of the transformations and , the relation between the original solution @xmath77 and the solution of the linear problem is as follows @xmath78      our analysis of the linearized equation is facilitated by introducing the transformation of variables @xmath79 where @xmath80 is a rescaling of time , defined by @xmath81 the map @xmath80 is invertible and its inverse is given by @xmath82 one computes that @xmath83 satisfy the linear system @xmath84 \frac{{\partial}{\hat \theta}}{{\partial}\tau }    - \kappa c_0 e^{\alpha \tau }   { \hat \theta}_{xx } & =   ( n+1 ) { \hat u}- \alpha { \hat \theta } , \\[2pt ] { \hat \sigma}&= n { \hat u}- \alpha { \hat \theta } , \end{aligned } \right .\ ] ] with boundary conditions @xmath85    for further reference , we record that the original variables @xmath86 and the solution @xmath87 of the rescaled problem - are related through the formulas @xmath88      we next consider the system @xmath89 \frac{{\partial}{\hat \theta}}{{\partial}\tau }    - k { \hat \theta}_{xx } & =   ( n+1 ) { \hat u}- \alpha { \hat \theta } , \end{aligned } \right .\ ] ] with boundary conditions .
note that in the thermal diffusion coefficient @xmath90 is time dependent , while in this coefficient is frozen in time .
we analyze the solutions of via an eigenvalue analysis .
this provides the following information :    1 .   for @xmath91 it provides the stability or instability properties of the linearized system .
2 .   for @xmath92
the eigenvalue analysis will only offer an indication of the effect of thermal diffusion on the linearized analysis .
this will be complemented by an energy estimate analysis of the exact system in the following section .
the solutions of - are expressed as a cosine fourier series @xmath93 where the fourier coefficients @xmath94 satisfy the ordinary differential equations @xmath95    the eigenvalues of the matrix @xmath96 are the roots of the binomial equation @xmath97    * the mode @xmath98*. for @xmath99 the equation becomes @xmath100 and has two eigenvalues @xmath101 and @xmath102 .
this corresponds to marginal stability for the linearized system .
however , any instability associated with the eigenvalue @xmath101 is eliminated by , which imply for the perturbation @xmath103 in the condition @xmath104 .
indeed , the system for @xmath105 takes the form @xmath106 this implies , using , , and the boundary condition , that @xmath107 and clearly @xmath108 decays .
thus , the @xmath109-th mode decays .    * the modes @xmath110*. for these modes
the eigenvalues are @xmath111 where the discriminant @xmath112 since @xmath113 we conclude that @xmath114 in summary :    1 .
if @xmath115 then all modes @xmath116 are unstable .
2 .   if @xmath92 the first few modes might be unstable , that is the modes @xmath117 for which @xmath118 , but the high frequency modes are stable .
3 .   the higher @xmath119 the fewer the number of unstable modes
. for @xmath119 larger than a certain threshold ( here equal to @xmath120 ) all modes are asymptotically stable .
we analyze the behavior of the linearized modes as a function of @xmath117 and discuss the nature of the instability in the various special cases of interest .
the analysis reveals the role of the parameters @xmath121 regarding the nature of the instability .      in this regime @xmath122
one eigenvalue is negative and the other strictly positive and increasing with @xmath117 . using the taylor series development @xmath123
we obtain the following asymptotic developments for the eigenvalues @xmath124    referring to the formulas , and , we see that the @xmath117-th mode of the perturbation grows to the leading order like @xmath125 clearly the perturbation grows much faster than the uniform shearing solution @xmath126 given by where the growth is logarithmic , and the behavior is characteristic of _ hadamard instability _ , that is the high frequency oscillations grow at a catastrophic rate .      in this case
the eigenvalues of the @xmath117-th mode are given by the formulas @xmath127 again we have one negative and one positive eigenvalue and thus unstable response in this regime as well .
using once again the taylor expansion we obtain the asymptotic expansions @xmath128 for the positive and negative eigenvalues respectively .
the positive eigenvalue @xmath129 satisfies the following properties : @xmath130    set @xmath131 and recall that @xmath132 satisfies the identity @xmath133 differentiating , we derive @xmath134 next , using the elementary bound @xmath135 for @xmath136 we obtain @xmath137    _ + ( x ) & = - + ( n x + ) < < .    @xmath137 this shows that @xmath138 and provides and .
then follows from .
the following remarks are in order : observe first that , by virtue of , and , the @xmath117-th mode of the perturbation grows like @xmath139 clearly the perturbation grows much faster than the uniform shearing solution @xmath126 , but the rate of growth is bounded in @xmath117 . on the other hand , the growth becomes very fast as the rate sensitivity @xmath140 .
the behavior in this regime is characterized by _
turing instability_. recall , that turing instability @xcite corresponds to the following scenario that occurs in some instances of dynamical systems .
there are two dynamical systems @xmath141 where the spectrum of both matrices @xmath142 and @xmath143 lies in the left - hand plane , but the trotter product of the two systems @xmath144 exhibits unstable response .
indeed , the dynamical system with @xmath115 may be visualized as the trotter product of the two systems @xmath145 and @xmath146 each of and has spectrum in the left half space but the eigenvalues of have unstable response .      in this regime
the eigenvalues are given by and they are real and have the properties @xmath147 therefore , the high modes are stable , but the low modes might be unstable if @xmath119 is sufficiently small .
the  most unstable " mode is the first mode .
one can again carry out the large @xmath117 asymptotics of the eigenvalues .
the result , in the case that @xmath148 , reads @xmath149 a similar formula holds when @xmath150 .
consider now the linearized system with boundary condition @xmath151 .
since we are interested in perturbations of the uniform shearing solution we will consider initial data that satisfy @xmath152 and thus focus on solutions with @xmath153 the objective is now to analyze the effect of the time - dependent diffusion on the dynamics of the system .
the linearized analysis of the system with frozen coefficients suggests to expect that diffusion stabilizes the system and its effect is intensified with the passage of time .
indeed , here we provide a rigorous proof of this statement via energy estimates .
[ lemdecay ] there exist @xmath154 and @xmath155 such that @xmath156 and @xmath157    we multiply @xmath158 by @xmath159 and @xmath160 by @xmath161 and obtain the identities @xmath162 using cauchy - schwarz inequality , gives @xmath163 in view of , we may use the poincare inequality @xmath164 combining the above , we note that for @xmath142 a positive constant we have @xmath165 next we select @xmath142 such that @xmath166 , and deduce that for @xmath142 thus selected there is a choice @xmath167 such that @xmath168    select now @xmath169 such that @xmath170 for @xmath171 ( recall from that @xmath172 when @xmath173 ) . on the interval @xmath171
, we have the differential inequality @xmath174 the latter implies and yields , upon using gronwall s inequality and , @xmath175 and concludes the proof    the next lemma provides control of the intermediate times .
[ lembound ] there exist @xmath176 and a constant @xmath177 such that @xmath178    using again and we have for any @xmath179 @xmath180 select now @xmath143 large so that @xmath181 for all @xmath182 , and then we obtain @xmath183 from where follows via gronwall s inequality .
we conclude now the analysis of stability .
suppose that for some @xmath184 the initial data satisfy @xmath185 and are controlled in @xmath186 by @xmath187 then , imply @xmath188 together they imply that the state @xmath189 is stable in @xmath186 , and in fact on account of it is asymptotically stable .
we state the result in the following :    for the initial - boundary value problem , with @xmath18 , if the initial perturbation @xmath190 satisfies @xmath191 then the equilibrium @xmath189 is asymptotically stable in @xmath186 .
the proof of the above lemmas and the analysis of the linearized problem with frozen coefficients indicate that the perturbation of the base solution may drift far from the origin until eventually it is recaptured as the effect of the diffusion intensifies with time .
this phenomenon of _ metastability _ also appears in numerical experiments for the nonlinear problem ; an instance of these experiments appears in figure [ mtstb ] .
we solve numerically the rescaled system for @xmath192 .
the initial conditions are perturbations of the uniform shearing solutions and they are : @xmath193 and @xmath194 , where @xmath195 is a small gaussian perturbation centered at the middle of the interval . the numerical method is based on adaptivity in space and time . an adaptive finite element method is used for the spatial discretisation while a runge - kutta method with variable time step is used for the discretisation in time .
further details concerning the numerical scheme and the adaptivity criteria can be found in @xcite . in figure [ mtstb ] the evolution in time ( in natural logarithmic scale ) is presented for the original variables @xmath35 and @xmath37 in ; @xmath37 is sketched in natural logarithmic scale .
we notice that , after a transient period of localization , thermal diffusion dominates the process , inhomogeneities of the state variables dissipate , and eventually the solution approaches a uniform shearing solution .
( left ) , @xmath196 ( right).,title="fig:",width=302,height=226 ]   ( left ) , @xmath196 ( right).,title="fig:",width=302,height=226 ]
we consider next the problem on the whole real line and put aside the boundary conditions . the objective is to construct special solutions that capture a process of localization .
this is done over three sections .
the present one is preliminary in nature and discusses various properties that make the construction feasible and the setting of the problem .
the main step is performed in section [ secss ] where a special class of self - similar solutions are constructed .
the localized solutions are then presented in section [ seclocalii ] .
we consider on @xmath197 and recall the notation @xmath198 .
the uniform shearing solutions are defined by solving .
introduce the transformation @xmath199 which involves a passage to relative perturbations ( as in section [ secrelpert ] ) and a change of time scale from the original to a rescaled time @xmath80 .
the scaling @xmath80 is precisely the one used in the linearized analysis of sections [ secchatime ] and [ seclinfroz ] ; it is defined by @xmath200 it satisfies @xmath201 where @xmath202 is as in , @xmath203 .
the scaling @xmath204 satisfies @xmath205 as @xmath206 , and @xmath80 is invertible with inverse map given by @xmath82 the functions @xmath207 satisfy the system of partial differential equations @xmath208
{ \hat \theta}_\tau - \kappa c_0 e^{\alpha \tau }   { \hat \theta}_{xx } & =   \big ( { \hat \sigma}\ , { \hat u}-1   \big ) , \\[2pt ] { \hat \sigma}&= e^{-\alpha { \hat \theta } } { \hat u}^n .
\end{aligned}\ ] ] we emphasize that    * for @xmath18 the system is non - autonomous , which is natural in view of the nature of the transformation . * notably , for @xmath43 the system becomes autonomous .
next , we consider the adiabatic case @xmath43 where takes the form @xmath209 { \hat \theta}_\tau   & =   \big ( { \hat \sigma}\ , { \hat u}-1   \big ) , \\[2pt ] { \hat \sigma}&= e^{-\alpha { \hat \theta } } { \hat u}^n \ , .
\end{aligned}\ ] ]    to understand the long time response of we perform a formal asymptotic analysis following ideas from @xcite .
fix an observational time scale @xmath169 and introduce the change of variables @xmath210 one easily checks that @xmath211 satisfies the system @xmath212 \tfrac{1}{t }   \ , { \tilde \theta}_s   & =   \big ( { \tilde \sigma}\ , { \tilde u}-1   \big ) , \\[2pt ] { \tilde \sigma}&= e^{-\alpha { \tilde \theta } } { \tilde u}^n .
\end{aligned}\ ] ] we are interested in calculating the equation describing the effective response of for @xmath169 sufficiently large .
this asymptotic analysis involves a parabolic equation @xmath158 , which can be thought as a moment equation , and equations @xmath160 and @xmath71 , which may be thought as describing a relaxation process towards the equilibrium curve described by @xmath213 note that the equilibria are a one - dimensional curve parametrized by @xmath214 .    to calculate the effective equation
we perform a procedure analogous to the chapman - enskog expansion of kinetic theory .
we refer to ( * ? ? ?
* sec 5 ) for the details of the asymptotic analysis .
the procedure produces the following effective equation @xmath215 which captures the effective response up to order @xmath216 .
the leading term is backward parabolic and the first order correction is fourth order .
one checks that the effect of the fourth order term is to stabilize the equation , in the sense that the linearized equation around the equilibrium @xmath217 is stable ( see @xcite for the details ) .
it was further shown in @xcite that admits special solutions that exhibit localization , which is not surprising as the leading response in is backward parabolic .
localized solutions of a similar nature will be constructed here directly for the system .
this is a more  unexpected " behavior , because is at least _ pro - forma _ a  hyperbolic - parabolic " system , and the explanation for this behavior lies precisely in the asymptotic relation between and .
it is easy to check that for @xmath43 the system is invariant under the one - parameter scaling transformation @xmath218 for any @xmath219 .
it is also easy to check that this is the only scaling transformation under which the problem is invariant .
this transformation has the distinctive feature that the scaling of space is independent of the scaling of time .
finally , the system is not scaling invariant when the thermal diffusion @xmath49 is positive .
the scaling invariance motivates a change of variables and a related _ ansatz _ of solutions .
introduce a variable in time scaling @xmath220 , the variable @xmath221 , and the change of variables @xmath222 we denote by @xmath223 and by @xmath224 . introducing to
, we obtain that the function @xmath225 satisfies the system of partial differential equations @xmath226 { \bar \theta}_\tau - \frac{\dot r}{r } \left ( \tfrac{n+1}{\alpha } + \xi { \bar \theta}_{\xi } \right ) , -   \kappa c_0 \frac{e^{\alpha \tau}}{r^2 ( \tau ) }    { \bar \theta}_{\xi \xi } & =   \big ( { \bar \sigma}\ , { \bar u}-1   \big ) , \\[2pt ] { \bar \sigma}&= e^{-\alpha { \bar \theta } } { \bar u}^n .
\end{aligned}\ ] ]    the resulting system depends on the choice of the scaling function @xmath227 .
we will be interested in solutions that are steady states in the new variables .
of course such solutions are dynamically evolving in the original variables .
there are two cases that produce simplified results of interest :    1 .
consider the case @xmath18 , and select @xmath228 .
then the system is consistent with the ansatz of steady state solutions @xmath229 and the latter are constructed by solving the system of ordinary differential equations @xmath230 -\frac{\alpha}{2 } \left ( \tfrac{n+1}{\alpha } + \xi { \bar \theta } ' \right )   -   \kappa c_0    { \bar \theta } ''   & =   { \bar \sigma}\ , { \bar u}-1 , \\[2pt ] { \bar \sigma}&= e^{-\alpha { \bar \theta } } { \bar u}^n .
\end{aligned}\ ] ] + since @xmath231 as @xmath232 increases , solutions of will produce defocusing solutions for .
a study of is expected to give information about the convergence to the uniform shearing state in the non adiabatic setting .
this is not the focus of the present work and we will not deal further with the system .
2 .   instead we consider the case @xmath43 .
we may now select any exponential function for @xmath227 but are most interested in @xmath233 with @xmath24 as this produces focusing solutions . with this selection ,
@xmath229 is a steady solution for with @xmath43 provided it satisfies @xmath234 \lambda \left ( \tfrac{n+1}{\alpha } + \xi { \bar \theta } ' \right )     & =   { \bar \sigma}\ , { \bar u}-1 , \\[2pt ] { \bar \sigma}&= e^{-\alpha { \bar \theta } } { \bar u}^n , \end{aligned}\ ] ] where @xmath24 is a parameter .
note that :    * the form of the transformation is motivated by the scaling invariance . *
the scaling connects the two variables and produces interesting systems even in the case @xmath49 positive when the system is not scale invariant .
we are ready to state and prove the main result of the paper .
this concerns the adiabatic case @xmath378 of system , taken now on the whole real line @xmath379 and expressed in terms of the variable @xmath19 in the form : @xmath380 our analysis is based on the transformations of section [ seclocal ] and the self - similar solution constructed in section [ secss ] .
[ shearband ] let @xmath381 and @xmath247 be given .
for any @xmath24 , the system admits a special class of solutions describing localization of the form : @xmath382 \label{locsigma } \sigma ( x , t ) & = \sigma_s ( t )   \left ( \tfrac{\alpha}{c_0 }   t + 1 \right ) ^ { - \tfrac{\lambda}{\alpha } } \sigma_\lambda \left ( \sqrt{\lambda }   x   \left ( \tfrac{\alpha}{c_0 } t + 1 \right ) ^{\tfrac{\lambda}{\alpha } }    \right ) , \\[2pt ] \label{loctheta } \theta(x , t ) & = \left ( 1 + \lambda \tfrac{n+1}{\alpha } \right ) \theta_s ( t ) - \lambda \tfrac{n+1}{\alpha } \theta_0 + \theta_\lambda \left ( \sqrt{\lambda }   x   \left ( \tfrac{\alpha}{c_0 } t + 1 \right ) ^{\tfrac{\lambda}{\alpha } }    \right ) , \end{aligned}\ ] ] where @xmath46 , @xmath45 are the uniform shear solutions in , @xmath203 , while @xmath383 is an even function defined on @xmath384 , which solves the initial value problem , with @xmath385 for @xmath386 .
[ rmkbeh ]    1 .    2 .
because @xmath387 goes to infinity in time , the level curves of @xmath388 accumulate along the @xmath389 axis in the @xmath390 plane . the functions , , thus describe a localizing solution where the flow tends to concentrate around @xmath33 as the time proceeds .
the solution , , is defined on @xmath391 and emanates from initial data @xmath392 \sigma ( x,0 ) & = \sigma_s ( 0 )    \sigma_\lambda   \big ( \sqrt{\lambda }   x     \big ) , \\[2pt ] \theta(x,0 ) & =   \theta_s ( 0 )   + \theta_\lambda   \big   ( \sqrt{\lambda }   x      \big )   \ , . \end{aligned}\ ] ] the initial data depend on two parameters : the parameter @xmath28 which can be thought as a length scale in initial data , and @xmath29 in which describes the amplitude of the initial perturbation .
this parametric family of solutions is valid for any @xmath24 .
if @xmath28 is large then the length scale of the initial perturbation is short while the growth in time and the rate of localization for the solution become fast .
we introduce the changes of variables and the stationary variant of , which once combined together reads @xmath393 we select @xmath394 with @xmath24 . according to the analysis in section [ secinvariance ] the function @xmath395 is sought as a solution of defined on the interval @xmath396 .
the desired @xmath395 is constructed as follows .
first , observe the relation between the system and the system .
let @xmath29 be fixed , and let @xmath397 be the solution of with @xmath26 and , defined for @xmath398 .
this solution is well defined by theorem [ connection ] . using a direct computation , the smoothness properties in ( iii ) of theorem [ connection ] , and
we conclude that    1 .
the function @xmath399 defined by @xmath400 satisfies on @xmath401 and the condition @xmath402 .
if we extend @xmath399 on @xmath403 by setting @xmath404 for @xmath405 , the resulting function @xmath406 , @xmath407 is a solution of defined on the whole real line that is even .    noting that @xmath408 we combine with and to arrive at , , .    in figure [ uqsol ] the profiles of @xmath409 and @xmath410 in and , respectively , are drawn at various instances of time from @xmath411 to @xmath412 . to achieve these profiles ,
we construct numerically the heteroclinic orbit for whose existence is justified in proposition [ thode ] . the construction is simple , if one exploits the special properties of system , and proceeds as follows : initial conditions @xmath413 are selected close to the equilibrium @xmath303 inside the region @xmath311 and in the direction of the stable manifold of @xmath284 for the linearized problem . then is solved backward in time , using an ode solver , and since @xmath414 is an unstable node ( for the forward problem ) the resulting orbit is a good approximation of the heteroclinic . for the integration of we use an explicit runge  kutta predictor - corrector method of order ( 4,5 ) with adaptive time stepping .
then @xmath29 is selected by choosing an appropriate parametrisation @xmath344 following the construction in .    , and @xmath410 , .,title="fig:",width=302,height=226 ] , and @xmath410 , .,title="fig:",width=302,height=226 ]    for the numerical runs the constitutive parameters are @xmath415 , @xmath416 .
the parameters related to the solution are @xmath417 , @xmath418 while for the selected reparametrization of the heteroclinic orbit the value of @xmath419 .
once the heteroclinic has been constructed , the profiles @xmath420 are computed via the changes of variables and . finally , @xmath409 and @xmath410 are computed by , respectively with @xmath417 , and are presented in fig [ uqsol ] .
in this section we consider the system of ordinary differential equations @xmath235 \nu \left ( \tfrac{n+1}{\alpha } + \xi \theta ' \right )     & =   \sigma \ , u -1 , \\[2pt ] \sigma & = e^{-\alpha \theta } u^n \ , , \end{aligned}\ ] ] on the domain @xmath236 with @xmath237 , subject to the initial condition @xmath238 this is an auxiliary problem that will be studied in detail in this section .
we emphasize that it does not coincide with the integrated form of , however it will be used together with additional properties in order to construct solutions of and localizing solutions for in the following section .    in this section ,
we are interested in constructing solutions @xmath239 for the singular initial value problem , for any @xmath237 .
the following elementary remarks are in order :    1 .   when @xmath240 , simplifies to @xmath241 \sigma ( 0 ) & = \sigma_0   \end{aligned }   \right .
\left \ {   \begin{aligned }    \sigma \ , u   & = 1   \\[2pt ] \theta   & = \frac{n+1}{\alpha } \log u   \end{aligned } \right . .\ ] ] the solution of is computed explicitly by @xmath242 2 .   when @xmath243 , the system is invariant under the family of scaling transformations @xmath244 for @xmath245 .
there is a special solution of that is self - similar with respect to the scaling invariance and reads @xmath246 this solution does not however satisfy the initial condition .
we next construct a function @xmath239 which satisfies , for all values of @xmath237 and @xmath247 .
we will show that the constructed solution will have the inner behavior of and the outer behavior of .
[ connection ] let @xmath247 be fixed .
let us note @xmath248 .
given @xmath237 there exists a solution @xmath249 of , defined for @xmath250 and satisfying the properties :    1 .
it achieves the initial data @xmath251 where @xmath29 as in and @xmath252 , @xmath253 are defined via @xmath254 2 .   as @xmath255
it has the limiting behavior @xmath256 \theta ( \xi ) & =    - \tfrac{n+1}{\alpha } \log \xi + o(1 ) .
\end{aligned}\ ] ] 3 .
the solution has the regularity @xmath257 , @xmath258 , @xmath259 .
it satisfies @xmath260 and has the taylor expansion as @xmath261 , @xmath262 where as usual @xmath263 ( but the rate may degenerate as @xmath264 tends to @xmath109 ) .
_ step 1 . _ motivated by the special solution we introduce a change of variables in two steps : @xmath265 \theta ( \xi ) & = - \tfrac{n+1}{\alpha } \log \xi   + { \bar \theta } ( \xi ) & & =   - \tfrac{n+1}{\alpha } \log \xi   + \theta ( \log \xi ) , \nonumber\end{aligned}\ ] ] indeed , we easily check that @xmath266 satisfies the system of equations @xmath267
\nu   \xi \ , { \bar \theta } '   & =   { \bar \sigma } \ , { \bar u } -1 , \\[2pt ] { \bar \sigma } & = e^{-\alpha { \bar \theta } } { \bar u}^n \ , .
\end{aligned}\ ] ] in the second step , we introduce in the change of independent variable @xmath268 and derive that @xmath269 satisfy @xmath270 \nu   \frac{d\theta}{d\eta }    & =   \sigma \ , u   -1 , \\[2pt ] \sigma   & = e^{-\alpha \theta } u^n \ , .
\end{aligned}\ ] ] the system is far simpler than the original in that it is autonomous and non singular .
we next consider and perform a further simplification , via the change of dependent variables @xmath271 note that @xmath272 .
we now compute @xmath273 and @xmath274    in summary , we conclude that @xmath275 satisfies the autonomous system of ordinary differential equations @xmath276 where @xmath277    there are two advantages in relative to :    1 .
is expressed via only two independent variables @xmath278 .
the equilibrium at infinity in is pulled at the axis for via the transformation .
we perform a phase space analysis of and establish the existence of a heteroclinic connection .
we show :    [ thode ] consider the autonomous planar system of differential equations with @xmath279 defined by . for any @xmath237
the system has the following properties :    * @xmath280 has two equilibria @xmath281 , @xmath282 in the first quadrant .
both @xmath283 and @xmath284 are hyperbolic points , @xmath283 is a repelling node while @xmath284 is a saddle point .
* there exists a heteroclinic orbit connecting @xmath283 to @xmath284 .
_ proof of proposition [ thode ] . _    1 .
the system has the special solution @xmath285 this solution splits the phase portrait into the left and right plane and acts as a barrier from crossing from the one side to the other .
the equilibria of are computed by solving @xmath286 and they are @xmath287 and @xmath288 and @xmath289 . focusing on the behavior in the first quadrant we neglect the equilibrium @xmath290 .
a computation shows @xmath291 -2 \tfrac{n+1}{n } a &    \tfrac{\alpha}{n\nu } c_\nu \end{pmatrix }   .\ ] ] at the equilibrium @xmath283 we have @xmath292 and we have eigenvalues @xmath293 with eigenvector @xmath294 , and @xmath295 with eigenvector @xmath296 . in view of
, we have @xmath297 .
the local behavior of solutions of around the equilibrium @xmath283 is that of a repelling node .
standard phase plane theory for second order systems ( _ e.g. _ @xcite ) states that there is a small neighbourhood of the equilibrium @xmath283 such that for any given solution @xmath298 of there are constants @xmath299 , @xmath300 such that @xmath301 as @xmath302 .
4 .   at the point
@xmath303 equation gives @xmath304 the eigenvalues @xmath305 , computed by @xmath306 are @xmath307 with eigenvector @xmath308 and @xmath309 with eigenvector @xmath310 .
@xmath284 is a saddle point .
5 .   consider the region @xmath311 bounded by the curves * @xmath312 is the @xmath313-axis defined by the equation @xmath314 , * @xmath315 is the horizontal line @xmath316 , * @xmath317 is the parabola @xmath318 .
+ the region @xmath311 is expressed as @xmath319 see figure [ fig : phase2 ] .
+    + we will show that @xmath311 is negatively invariant .
+ indeed , the flow of is parallel to @xmath312 as indicated in ( b ) .
on the horizontal line @xmath320 the system gives @xmath321 on the curve @xmath322 the system gives @xmath323 therefore @xmath311 is negatively invariant . 6 .   for completeness
we also show that the vector @xmath324 associated to the eigenvalue @xmath325 at @xmath284 points towards the interior of the region @xmath311 .
recall that @xmath284 is the upper right corner of the region @xmath311 and that * the tangent vector to @xmath315 is @xmath326 , * the tangent vector to @xmath317 at @xmath303 is @xmath327 , * the eigenvector @xmath328 .
+ we will show that @xmath329 this then implies that the vector @xmath330 when placed at the corner @xmath284 points strictly towards the interior of @xmath311 . to show recall that @xmath307 and use and to rewrite @xmath331 where @xmath332 this shows @xmath333 and concludes .
it is easy to check that there can be no closed limit cycle in @xmath311 because in that region @xmath334 the poincar - bendixson theorem implies that there exists a heteroclinic orbit joining @xmath335 and @xmath336 which emanates along the direction of the stable manifold at @xmath336 and progresses backwards in time towards the node at @xmath335 .
we turn now to defining the solution of , .
let @xmath337 stand for a parametrization of the heteroclinic orbit .
in view of and the fact that @xmath338 , there exists a constant @xmath299 such that @xmath339 due to the fact that the vertical axis is a trajectory of its own and thus distinct from the trajectory we are considering , and due to the form of the phase portrait we can assert that @xmath340 .    for any @xmath341
the function @xmath342 offers an alternative parametrization of the same heteroclinic orbit . applying to this parametrization yields @xmath343 choose @xmath344 such that @xmath345 and let @xmath298 denote the corresponding parametrization of the heteroclinic : @xmath346 the selected parametrization satisfies as @xmath302 @xmath347 where as usual @xmath348 as @xmath302 . recalling the changes of variables and we define @xmath349 the function clearly satisfies the system and , since @xmath350 as @xmath351
, it satisfies @xmath352 \theta ( \xi ) & =    - \tfrac{n+1}{\alpha } \log \xi + o(1 ) , \end{aligned}\ ] ] as @xmath353 goes to infinity , which shows .    in the last step ,
we study the behavior near the ( apparent ) singular point at @xmath354 and prove that the solution is in fact regular there and satisfies , and .
clearly and @xmath160 imply that @xmath355 as @xmath261 .
let @xmath252 and @xmath253 be defined by . since @xmath356 , gives @xmath357 now , , , imply @xmath358 and together with and provide .
observe now that @xmath359 using , and , we obtain @xmath360 \\ & = \tfrac{1}{n } u^3 b - u \tfrac{\alpha}{\nu n }   \frac{c_\nu}{b } \big ( b - \frac{1}{c_\nu } \big )   \\ & = \tfrac{1}{n } \tfrac{u_0 ^ 3}{c_\nu } e^{3 \eta}(1 + o(1 ) ) -
\tfrac{\alpha}{\nu n }   u_0 o(1 ) e^{2\eta }   \qquad \mbox{as $ \eta \to -\infty$ } \end{aligned}\ ] ] and thus @xmath361 since @xmath362 \cap c^1(0,1)$ ] the mean value theorem implies @xmath363 and implies @xmath364 exists and @xmath365 .
we conclude that @xmath366 has the taylor expansion @xmath158 . in turn ,
@xmath367 , @xmath368 , @xmath369 and @xmath370 achieves the taylor expansion @xmath160 .    turning to @xmath371 ,
observe that by and @xmath71 @xmath372 while @xmath373 in addition , @xmath374 combined with the formula @xmath375 gives @xmath376 this shows @xmath377 has the taylor expansion @xmath71 at the origin .
we considered the system as a paradigm to study the onset of localization and formation of shear bands and to examine the role of thermal softening , strain - rate hardening and thermal diffusion in their development .
the stability of the uniform shear solution leads , in general , to the analysis of a non - autonomous system .
the present model has the remarkable property that the system of relative perturbations is autonomous in the adiabatic case ( @xmath421 ) .
a study of linearized stability for the case with no thermal diffusion @xmath378 shows that : ( i ) for @xmath17 the high frequency modes grow exponentially fast , with the exponent proportional to the frequency , what is characteristic of ill - posedness and hadamard instability .
( ii ) when @xmath2 the modes are still unstable but the rate of growth is bounded independently of the frequency .
the behavior in that range is that characteristic of turing instability .
the effect of thermal diffusion ( @xmath18 ) for the non - autonomous linearized model was also assessed ; it was shown that perturbations grow in the early stages of deformation , but over time the compound effect of diffusion stabilizes the response of the system .
then we turned to the nonlinear system and constructed a class of self - similar solutions on the real line , leading to the explicit solutions , , .
their construction is based on the self - similar profiles of , , in turn determined by a suitable heteroclinic connection .
they depend on three parameters : @xmath27 linked to the uniform shear in , a parameter @xmath28 which can be viewed as a length scale of initial data , and @xmath29 determining the size of the initial profile .
the solutions presented in fig [ uqsol ] evidently exhibit localisation : as time increases a coherent structure forms with the deformation localizing in a narrow band . in this self - similar solution
there is no trace of the oscillations .
a comparison with the modes of the linearized problem indicates that the nonlinearity suppresses the oscillations and a seamless coherent structure emerges as the net outcome of hadamard instability , rate dependence , and the nonlinear structure of the problem .
, _ localization and shear bands in high strain - rate plasticity_. in nonlinear conservation laws and applications , a. bressan , g .- q .
chen , m. lewiscka and d. wang , eds ; i m a vol . math .
appl . , 153 , springer , new york , 2011 , pp 365 - 378 .