mathematical modeling is a prominent tool used to better understand complex mechanical or biological systems @xcite .
a common problem that arises when developing a model of a biological or mechanical system is that some of its parameters are unknown .
this is especially important when those parameters have special meaning but can not be directly measured .
thus a natural question arises : can all , or at least some , of the model s parameters be estimated indirectly and _ uniquely _ from observations of the system s input and output ?
this is the question of _ structural identifiability_. sometimes the uniqueness holds only within a certain range . in this case , we say that a system is only _ locally _ structurally identifiable .
there are numerous reasons why one would be interested in establishing identifiability .
structural identifiability is a necessary condition for the practical or numerical identifiability problem , which involves parameter estimation with real , and often noisy , data .
the unobservable biologically meaningful parameters of a model can only be determined ( or approximated ) if the model is structurally identifiable .
moreover , optimization schemes can not be employed reliably since they will find difficulties when trying to estimate unidentifiable parameters @xcite .
the concept of structural identifiability was introduced for the first time in the work of bellman and strm @xcite . since then
, numerous techniques have been developed to analyze the identifiability of linear and nonlinear systems with and without controls @xcite ; see also @xcite for a review of different approaches .
viscoelastic mechanical models that utilize springs and dashpots in various configurations have been widely used in numerous areas of research including material sciences @xcite , computer graphics @xcite , and biomedical engineering to describe mechanical properties of biological systems @xcite . to achieve a desirable response , networks with different numbers of springs and dashpots in various configurations have been constructed .
for example , it is well - known that the simplest models of viscoelastic materials such as voigt ( spring and dashpot in parallel ) or maxwell ( spring and dashpot in series ) do not offer satisfactory representation of the nature of real materials @xcite .
thus more complicated configurations are usually constructed and analyzed @xcite .    in this paper
we investigate the identifiability problem of viscoelastic models represented by an arbitrarily complex spring - dashpot network .
although there exist numerous methods that can determine the type of identifiability of a system of ordinary differential equations , generally they are difficult to apply
. our results will show in a remarkably simple way how to verify whether the studied model is ( locally or globally ) structurally identifiable . in case
it is unidentifiable , our method provides an explanation why this is the case and how to reformulate the problem .
moreover , the existing methods usually allow to establish the identifiability only _ a posteriori _ , i.e.  after concrete systems have been established
thus , we also introduce `` identifiability tables '' , which allow not only to check but also to construct an arbitrarily complex _
identifiable _ spring - dashpot network .
a particular motivation for this work comes from cardiovascular modeling @xcite , although the results of this paper can be applied to any viscoelastic modeling approach .    _
arterial wall .
_ changing blood pressure causes periodic expansion and contraction of the arterial wall ( see fig .  [
fig : cv ] ) .
it is well - known that the stress - strain curves of the artery walls exhibit hysteresis , which is understood to be a consequence of the fact that the wall is viscoelastic .
another manifestation of the viscoelasticity of the arterial tissue is the stress relaxation experiments under constant stretch ( strain ) .
spring - dashpot ( s - d ) networks are often used in order to describe the biomechanical properties of the arterial tissue @xcite .
identifiable networks can be determined using the results of this paper ( see theorem [ thm : parameqcoeff ] ) .
_ neural activity .
_ it is common to use the spring - dashpot network to describe the neural firing of various sensors ( e.g.  muscle spindle , baroreceptors ) , see @xcite .
typically one assumes that the firing activity is proportional to the strain sensed by a spring connected in series with a spring - dashpot network , which represents a local integration of the nerve endings to the arterial wall ( see fig .
[ fig : cv ] ) . then the arterial wall and neural activity models are combined . although separately each model is structurally identifiable , there is no guarantee that the resulting viscoelastic structure is identifiable .
thus , using our results given in theorem [ thm : localmt ] , we can establish whether the combined viscoelastic model is identifiable , and if not , what needs to be modified .
after reviewing basic concepts of viscoelasticity of systems , we present and discuss our main results related to local and global structural identifiability of such systems .
finally , we illustrate our results with a number of examples from the literature .
the ideal linear elastic material follows hooke s law @xmath0 , where @xmath1 is a young s modulus ( or a spring constant ) , which describes the relationship between the stress @xmath2 and the strain @xmath3 .
analogously , the relation @xmath4 describes the viscous material , where @xmath5 and @xmath6 is a viscous constant @xcite . in the basic linear viscoelasticity theory ,
the elastic and viscous elements are combined . in this work ,
we shall be concerned with the problem of identifiability of networks of springs and dashpots that are essentially one - dimensional .
the elements can be combined either in series or in parallel . in order to obtain the relationship between the total stress ( force ) @xmath2 and the total strain ( extension ) @xmath3 for a given spring - dashpot network , we use two fundamental rules . for two viscoelastic elements connected in series ,
the stress is the same in both elements , but the total strain is the sum of individual strains on each element . on the other hand , for elements connected in parallel
, the strain is the same for both elements , but the total stress is the sum of individual stresses on each element .
now we consider concrete viscoelastic networks , starting with the simplest configurations .
[ ex : m ] the series combination of a spring , denoted by its constant @xmath1 , and a dashpot , denoted by its constant @xmath6 , is known as a maxwell element ( see fig .  [
fig : simple](a ) ) . since the elements are connected in series , the stress @xmath2 is the same on both elements and the total strain @xmath3 is the sum of strains @xmath7 and @xmath8 corresponding to the spring and dashpot , respectively .
now , the relationship between the total strain and stress for this system is @xmath9    [ ex : v ] another simple example is the voigt element ( also known as kelvin or kelvin - voigt ) given in fig .
[ fig : simple](b ) . following the steps outlined in the previous example ,
we obtain the @xmath10 relationship @xmath11    [ ex : burgers ] in our third example we consider a particularly popular four - element model , represented by a maxwell element combined in series with a voigt element , and known as the burgers model ( fig .
[ fig : simple](c ) ) .
denote by subscript @xmath12 and @xmath13 the spring and viscous constants of the maxwell and voigt elements , respectively .
note that the stress @xmath2 is the same on all three elements connected in series ( voigt , spring and dashpot ) .
eliminating the corresponding local strains , we obtain the following relationship @xmath14\dot\sigma+\frac{e_m e_v}{\eta_m\eta_v}\sigma.\ ] ]      first note ( cf .
examples [ ex : m ] , [ ex : v ] , and [ ex : burgers ] ) that for any configuration of springs @xmath15 and dashpots @xmath16 , the total strain  stress relationship can always be written as the following @xmath17-st order linear ordinary differential equation @xmath18 where the coefficients @xmath19 and @xmath20 are functions of the spring and dashpot constants .
the precise value of @xmath21 and the forms of @xmath22 and @xmath23 will depend on the particular structure of the spring - dashpot model .
equation [ eq : ode ] is known as the _ constitutive equation_. in the context of spring - dashpot networks , identifiability concerns whether or not it is possible to recover the unknown parameters ( @xmath24 and @xmath25 ) of the system from the governing equation of the model , given only the total stress @xmath2 and total strain @xmath3 . in other words
, we assume that we know the stress and the strain at the bounding nodes only and ask if it is possible to determine the unknown parameters ( @xmath24 and @xmath25 ) . in order to uniquely fix the coefficients of the constitutive equation
, we require that be _
normalized _ so that the leading term ( in @xmath2 or @xmath3 , depending on the situation ) is monic . thus , letting the @xmath26 non - monic coefficients of be represented by the vector @xmath27 , we have the following formal definition of identifiability .
[ defn : id ] let @xmath28 be a function @xmath29 , where @xmath30 is the parameter space .
the model is _ globally identifiable _ from @xmath28 if and only if the map @xmath28 is one - to - one .
the model is _ locally identifiable _ from @xmath28 if and only if the map @xmath28 is finite - to - one .
the model is _ unidentifiable _ from @xmath28 if and only if the map @xmath28 is infinite - to - one .
note that local identifiability is equivalent to saying that around each point in parameter space there exists a neighborhood on which the function @xmath28 is one - to - one .
for example , for the burgers model considered in example [ ex : burgers ] , the coefficient function @xmath31 is defined as @xmath32 technically speaking , in this paper we will consider the slightly weaker notion of _ generic global identifiability _ ( or generic local identifiability , or generic unidentifiability ) , where _ generic _ means that the property holds almost everywhere .
we will omit the use of the term generic when speaking of identifiability .    definition [ defn : id ] implies that if there are more parameters than non - monic coefficients , then the system must be unidentifiable .
thus , a necessary condition for structural identifiability is that the number of parameters @xmath33 ( elements of the network ) is less than or equal to the number of non - monic coefficients in the constitutive equation .
we will soon show that the number of non - monic coefficients is bounded by the number of parameters in spring - dashpot networks .
thus , in this case , a necessary condition for structural identifiability is that the number of parameters and non - monic coefficients are equal .
we will prove that , remarkably , in the case of viscoelastic models represented by a spring - dashpot network , the converse to this statement holds as well .
[ local identifiability ] [ thm : parameqcoeff ] a viscoelastic model represented by a spring - dashpot network is locally identifiable if and only if the number of non - monic coefficients of the corresponding constitutive equation equals the total number of its moduli @xmath34 and viscosity parameters @xmath35 .
note that although the constitutive equation is a linear differential equation , its coefficients considered as functions of spring and viscous constants are not linear functions of the parameters ( see ) .
thus , theorem [ thm : parameqcoeff ] allows to reduce the difficult problem of checking one - to - one or finite - to - one behavior of nonlinear functions to simply counting the number of parameters ( springs and dashpots ) and non - monic coefficients of the constitutive equation and asking whether the two numbers are equal .
the positive answer implies local identifiability , whereas a negative answer implies unidentifiability .
consider , for example , the maxwell and voigt elements , and the burgers model .
we note that the constitutive equations , , and for all three models are already in the normalized form .
now , simply by counting the number of parameters and the non - monic coefficients of the constitutive equations , we see that the two are equal for each model .
thus , by the above theorem , all three models are locally structurally identifiable .
now we examine when combining two identifiable models results also in an identifiable model .
this will allow us to construct arbitrarily complex and identifiable spring - dashpot networks .
we start with an observation , which we prove in the following section , related to the possible form of any differential equation that describes a spring - dashpot network .    [
prop : types ] every spring - dashpot network , given by equation , has one of the four possible types @xmath36    recall that @xmath21 is the highest derivative of the stress component @xmath37 , which appears in the total strain - stress equation .
now we illustrate the different types of networks defined in the above proposition by considering the simplest elements .    for a spring , given by @xmath38 , we have @xmath39 ( only @xmath2 appears in the constitutive equation , but none of its derivatives ) , @xmath40 , and @xmath41
. therefore a spring is of type a. note that for a dashpot , which is given by @xmath42 , we also have @xmath39 but @xmath43 and @xmath44 .
thus , according to notation given in proposition [ prop : types ] , a dashpot is of type b. for the voigt element , given by , we have @xmath39 as well as @xmath45 and @xmath46 ( that is @xmath47 ) .
we conclude that it is of type c. finally , a maxwell element is given by .
note that here @xmath48 ( since @xmath49 appears in ) , @xmath50 , @xmath51 , and @xmath52 .
thus a maxwell element is of type d.    once the constitutive equation has been determined for a given spring - dashpot network , it is very easy to establish the type that it belongs to .
unfortunately , @xmath21 does not always have a physical significance .
the value of @xmath21 is determined by the specific network and can not be easily related to the number of springs and dashpots as we will illustrate later on .
[ thm : localmt ] consider two locally identifiable spring - dashpot systems @xmath53 and @xmath54 of one of the four types @xmath55 , @xmath56 , @xmath57 , @xmath58 .
then the new model resulting in joining @xmath53 and @xmath54 either in parallel or in series is of the type indicated by the identifiability tables ( table [ thetable ] ) . the letter _ u _ indicates that the network is unidentifiable , otherwise it is identifiable of the given type .
there are several ways one could use the above theorem .
one way is to establish the local identifiability of a given spring - dashpot network .
contrary to our similar result given in theorem [ thm : parameqcoeff ] , this can be done without actually calculating the constitutive equation .
we will show how to apply theorem [ thm : localmt ] to establish structural identifiability after first introducing some notation .
given any two spring - dashpot models @xmath59 and @xmath60 , we use the following notation @xmath61 and @xmath62 to denote respectively the parallel and series combination of @xmath59 and @xmath60 .
let @xmath63 denote the function that takes a spring and dashpot model @xmath59 and outputs its type ( @xmath64 ) if it is locally identifiable , and @xmath65 if it is unidentifiable . to apply @xmath63 to a complicated model built up from springs and dashpots using series and parallel connections
, we replace any springs and dashpots with their respective types @xmath55 and @xmath56 as well as the operations @xmath66 and @xmath67 with @xmath68 and @xmath69 , respectively .
then we apply the operations in the identifiability tables ( see table [ thetable ] ) .
when connecting two identifiable spring - dashpot networks of one of the types @xmath55 , @xmath56 , @xmath57 , @xmath58 , or an unidentifiable @xmath65 either in series or in parallel , the above tables establish the type of the resulting identifiable system .
if the resulting structure is unidentifiable it is indicated by @xmath65 .
for example , a parallel connection of two networks of types @xmath55 and @xmath58 gives rise to an identifiable network of type @xmath55 ( see ( a ) ) , but the series connection results in an unidentifiable structure ( see ( b ) ) .
[ thetable ]    note that the maxwell model shown in fig .
[ fig : simple](a ) can be symbolically written as @xmath70 in this formula , we simply replace the spring and the dashpot with @xmath55 and @xmath56 , respectively , as well as the operations @xmath66 and @xmath67 with @xmath68 and @xmath69 , respectively , to obtain @xmath71 thus we conclude that the maxwell model is locally identifiable and is of type @xmath58 .
[ ex : burgers ] similarly , the burgers model shown in fig .
[ fig : simple](c ) can be symbolically written as @xmath72 to check the local identifiability , we find @xmath73 and use table [ thetable ] to obtain @xmath74 we conclude that the burgers model is locally identifiable and of type @xmath58 .    in the next example we show how we can easily establish local structural identifiability of a more complicated network .
[ ex : turski ] consider a viscoelastic material studied in @xcite and represented by a spring - dashpot network shown in fig .
[ fig : turski](a )
. it can be symbolically represented by @xmath75\land e_4.\ ] ] again , we can verify the local identifiability of the above model using table  [ thetable ] and obtain @xmath76\odot a\\ & = \big[(((\overbrace{(c\odot a)}^a\odot b)\oplus b)\odot a)\oplus b\big]\odot a\\ & = \big[((\overbrace{(a\odot b)}^d\oplus b)\odot a)\oplus b\big]\odot a=\ldots = d\\ \end{aligned}\ ] ] this simple computation confirms that the model is locally structurally identifiable .    our method can also verify if a network is unidentifiable , providing the reason for the lack of its identifiability . consider the following example .
consider a viscoelastic model used in @xcite and shown in fig .
[ fig : turski](b ) . using the notation previously introduced ,
it can symbolically be written as @xmath77\land e_5\land \eta_5.\ ] ] now applying table [ thetable ] , we obtain @xmath78\odot \overbrace{a\odot b}^d\\ & = \big[(\overbrace{(d\oplus d)}^d\odot d)\oplus d\big]\odot d\\ & = \big[\overbrace{(d\odot d)}^{u}\oplus d\big]\odot d=\emph{u}. \end{aligned}\ ] ] whenever table [ thetable ] indicates @xmath65 ( i.e.  the corresponding substructure is unidentifiable ) , this inevitably leads to the whole model being unidentifiable .
moreover , our method can also explain what is the reason for the lack of identifiability .
in this example the situation is simple : joining in series a maxwell element ( type d ) with a generalized maxwell model leads to an unidentifiable network .
so far we have considered only _ local _ identifiability of mechanical systems .
now we complete the presentation and discussion of our results by introducing a criterium , which establishes when a given network is _ globally _ structurally identifiable .
[ thm : global ] a viscoelastic model represented by a spring - dashpot network is globally identifiable if and only if it is locally identifiable and the network is constructed by adding either in parallel or in series at the bounding nodes exactly one basic element ( spring or dashpot ) at a time .    note that the network given in fig .
[ fig : turski](a ) and considered in example [ ex : turski ] was deemed locally structurally identifiable .
we note that it can be constructed by adding just one element at a time and therefore it is _ globally _ structurally identifiable .
similarly , all the simple models shown in fig .
[ fig : simple ] can also by constructed adding only one element at a time , and since they are locally identifiable , we conclude that they are also globally structurally identifiable .
now consider a model which is locally , but not globally , structurally identifiable .
consider a generalized kelvin - voigt model shown fig .
[ fig : turski](c ) and used in @xcite ) in the context of cardiovascular modeling .
it can be symbolically represented by @xmath79 thus the local identifiability can be checked by computing @xmath80 we immediately conclude that the network is locally identifiable . in order to verify whether it is also globally identifiable , note that this network _ can not _ be constructed by adding only one element at a time .
thus the system is only locally , but not globally , identifiable .
however , in this case the non - global identifiability arises from merely permuting the parameters among the three voigt elements .
in this section , we prove the main results from the previous section . to do this
requires a careful analysis of the structure of the constitutive equation after combining a pair of systems in series or in parallel .
let @xmath53 and @xmath54 be spring - dashpot models whose respective constitutive equations are @xmath81 and @xmath82 , where @xmath83 represent linear differential operators .
we can write the differential operators ( in general form ) as :    @xmath84    table [ shapestable ] shows that there are restrictions on the values of the @xmath85 and @xmath86 , e.g.  the differential order of the lowest order term in @xmath2 is always zero and the differential order of the lowest order term in @xmath3 is zero or one , but we leave the operators in general form for simplicity .
we now show the form of the resulting constitutive equation after combining these systems in series or in parallel , in terms of these differential operators . in
what follows , we will treat the differential operators @xmath83 as polynomial functions in the variable @xmath87 . for example , @xmath88 can be thought of as a polynomial @xmath89 .
suppose that @xmath90 is a series connection of models @xmath53 and @xmath54 , whose constitutive equations are @xmath91 and @xmath92 , respectively .
then the stresses ( @xmath2 ) are the same for the two systems while the strains ( @xmath3 ) are added .
if @xmath88 and @xmath93 are relatively prime , then the constitutive equation of @xmath59 is : @xmath94 we assume that @xmath95 , so that the constitutive equation is monic .
if @xmath88 and @xmath93 have a common factor , then the constitutive equation of @xmath59 is obtained by dividing ( [ eq : series ] ) by the greatest common divisor of @xmath88 and @xmath93 .
suppose that @xmath96 is a parallel connection of models @xmath53 and @xmath54 , whose constitutive equations are @xmath91 and @xmath92 , respectively .
then the strains ( @xmath3 ) are the same for the two systems while the stresses ( @xmath2 ) are added .
if @xmath97 and @xmath98 are relatively prime , then the constitutive equation is : @xmath99 we assume that @xmath100 , so that the constitutive equation is monic .
if @xmath97 and @xmath98 have a common factor , then the constitutive equation is obtained by dividing ( [ eq : parallel ] ) by the greatest common divisor of @xmath97 and @xmath98 .
now we prove proposition [ prop : types ] , that is , we show that every spring - dashpot network , given by equation , has one of the four possible types displayed in table [ shapestable ] , which are defined by the _ shapes _ of the linear operators @xmath83 acting on @xmath3 and @xmath2 .
we make this notion precise :    the _ shape _ of a linear operator @xmath83 is a pair of numbers , written @xmath101 $ ] , where @xmath85 is the highest differential order and @xmath86 is the lowest different order .    .
* possible types of constitutive equations * [ cols="<,<,<,<,<",options="header " , ]     two systems of types a , b , c , or d are combined in parallel , where in the first system @xmath102 and in the second system @xmath103 .
[ paralleltable ]    we note that if a type b or d is combined in series with a type b or d , then @xmath88 and @xmath93 have a common factor ( since both lacked a constant term ) , so the equation @xmath104 is divided by @xmath105 to arrive at the shapes listed in the table .
in addition to the type of equation that results after combining two equations of types @xmath106 , we have in tables [ seriestable ] and [ paralleltable ] the resulting identifiability properties of each equation , which we will obtain in the next section .
note that definition [ defn : id ] implies that if there are more parameters than non - monic coefficients , then the system must be unidentifiable .
the tables show that the number of non - monic coefficients is bounded by the number of parameters , thus a necessary condition for identifiability is that the number of parameters equals the number of non - monic coefficients in the constitutive equation . in the next section ,
we show that this is also a sufficient condition .
consider a spring - dashpot system @xmath59 whose final step connection is a series connection of two systems @xmath53 and @xmath54 , i.e.  @xmath90 . since the number of non - monic coefficients in any spring - dashpot model is always less than or equal to the number of parameters in that model , we know that a necessary condition for this system to be locally identifiable is that @xmath53 and @xmath54 are both locally identifiable .
let @xmath91 be the constitutive equation for @xmath53 and @xmath92 be the constitutive equation for @xmath54 .
each of the operators @xmath107 , @xmath108 , @xmath109 , and @xmath110 will have a fixed shape determined by the structure of @xmath53 and @xmath54 .
assuming that @xmath53 and @xmath54 are locally identifiable , we can choose parameters in each of the models @xmath53 and @xmath54 so that the coefficients of these constitutive equations are arbitrary numbers .
thus , deciding identifiability of this system amounts to determining whether the map that takes the pair of equations @xmath111 to the constitutive equation @xmath112 , where @xmath113 , @xmath114 , @xmath115 , and @xmath116 ( cf . ) , for the system @xmath59 is finite - to - one or not .
the same reasoning works _
mutatis mutandis _ for parallel connections , where we now concern ourselves with the map from the pair of equations @xmath111 with generic coefficients to the constitutive equation for @xmath117 given in .    to make the above , intuitive , statements precise we introduce the following definition .
the _ shape factorization problem _ for a quadruple of shapes @xmath118 , [ n_{2 } , m_{2 } ] , [ n_{3 } , m_{3 } ] , [ n_{4 } , m_{4}])\ ] ] is the following problem : for a generic pair of polynomials @xmath119 with @xmath120 monic such that the @xmath121 $ ] and @xmath122 $ ] , do there exist finitely many quadruples of polynomials @xmath123 with @xmath124 $ ] , @xmath107 and @xmath109 are monic , and such that @xmath125 and @xmath126 ?
a quadruple of shapes @xmath127 is said to be _ good _ if the shape factorization problem for that quadruple has a positive solution .    since the above definition introduces one of the key concepts of the paper , in the following example we shall further illustrate the meaning of the shape factorization problem .
suppose that our quadruple @xmath128 , [ n_{2 } , m_{2 } ] , [ n_{3 } , m_{3 } ] , [ n_{4 } , m_{4 } ] ) =   ( [ 2,0 ] , [ 2,0 ] , [ 3,0],[2,0])\ ] ] which is a special case of joining models of types @xmath55 and @xmath57 in series . the shape factorization problem in this case asks the following question :    let @xmath119 be a generic pair of polynomials where @xmath120 and @xmath129 are degree @xmath130 polynomials with nonzero constant term and @xmath120 is monic : @xmath131 do there exist finitely many polynomials @xmath132 such that @xmath125 and @xmath133 ? or to say it another way , for generic values of @xmath134 and @xmath135 , does the system of @xmath136 equations in @xmath136 unknowns : @xmath137 have only finitely many solutions ?    the language of shape factorization problems and the remarks in the preceding paragraphs allow us to reduce the local identifiability problem for a spring - dashpot system to determining whether a certain quadruple is a good quadruple .
let @xmath90 be a spring - dashpot model joined in series from @xmath53 and @xmath54 , where @xmath53 has constitutive equation @xmath138 of shapes @xmath139 $ ] and @xmath140 $ ] , respectively , and @xmath54 has constitutive equation @xmath141 of shapes @xmath142 $ ] and @xmath143 $ ] , respectively .
then the model @xmath59 is locally identifiable if and only if    1 .
@xmath53 and @xmath54 are locally identifiable , and 2 .
@xmath144 , [ n_{2 } , m_{2 } ] ,   [ n_{3 } , m_{3 } ] , [ n_{4 } , m_{4}])$ ] is a good quadruple .    similarly , if @xmath96 is a spring - dashpot model joined in parallel from @xmath53 and @xmath54 , then @xmath59 is locally identifiable if and only if    1 .
@xmath53 and @xmath54 are locally identifiable , and 2 .
@xmath145 , [ n_{1 } , m_{1 } ] ,     [ n_{4 } , m_{4 } ] , [ n_{3 } , m_{3}])$ ] is a good quadruple .
so what remains to show is that , for the shapes that arise in spring - dashpot models , whether a quadruple of shapes is a good quadruple only depends on the types ( @xmath146 or @xmath58 ) of the systems being combined .
the proof of this statement will occupy the rest of this section .
let @xmath120 and @xmath129 be two polynomials .
note that for given fixed shapes , @xmath139 $ ] and @xmath142 $ ] , there are at most finitely many factorizations @xmath125 , where @xmath107 has shape @xmath139 $ ] and @xmath109 has shape @xmath142 $ ] and both are monic .
this is because there are at most finitely many ways to factorize a monic polynomial into monic factors .
once we fix one of these finitely many choices for @xmath107 and @xmath109 , the equation @xmath147 is a linear system in the ( unknown ) coefficients of @xmath108 and @xmath110 .    for a polynomial @xmath148 of shape
@xmath149 $ ] , we can write the coefficients of @xmath120 as a vector , which we denote @xmath150 : = \begin{pmatrix } f_{n } \\ \vdots \\
f_{m } \end{pmatrix}.\ ] ]    let @xmath83 have shape @xmath101 $ ] , as defined in equation .
the vector of coefficients of @xmath151 can be written as the result of a matrix vector product as : @xmath152= \begin{pmatrix } a_{n_1 } & 0 & \cdots & 0 \\ \vdots & a_{n_1 } & \cdots & 0 \\ a_{m_1 } & \vdots & \cdots & \vdots \\   0 & a_{m_1 } & \cdots & 0 \\ \vdots & 0 & \cdots & a_{n_1 } \\ \vdots & \vdots & \cdots & \vdots \\ 0 & 0 & \cdots & a_{m_1 }    \end{pmatrix } \begin{pmatrix } e_{n_4 } \\ \vdots \\
e_{m_4 } \end{pmatrix}.\ ] ] we will refer to this product as @xmath153 $ ] , where @xmath154 is a @xmath155 by @xmath156 matrix .
likewise , the coefficients of @xmath157 can be written as the result of a matrix vector product as : @xmath158= \begin{pmatrix } c_{n_3 } & 0 & \cdots & 0 \\ \vdots & c_{n_3 } & \cdots & 0 \\ c_{m_3 } & \vdots & \cdots & \vdots \\   0 & c_{m_3 } & \cdots & 0 \\ \vdots & 0 & \cdots & c_{n_3 } \\ \vdots & \vdots & \cdots & \vdots \\ 0 & 0 & \cdots & c_{m_3 }    \end{pmatrix } \begin{pmatrix } b_{n_2 } \\ \vdots \\
b_{m_2 } \end{pmatrix}.\ ] ] we will refer to this product as @xmath159 $ ] , where @xmath160 is a @xmath161 by @xmath162 matrix . then we call the _ matrix factored form _ of @xmath163 $ ] the expression : @xmath164+h[l_4],\ ] ] where the matrices @xmath165 and @xmath166 are the matrices @xmath160 and @xmath154 padded with rows of zeros so that coefficients corresponding to monomials of the same degree appear in the same row
. this makes @xmath167 a @xmath168 by @xmath169 matrix .
we can now state a criteria for determining if the shape factorization problem has finitely many solutions :    [ prop : abinvert ] the quadruple @xmath144 , [ n_{2 } , m_{2 } ] ,   [ n_{3 } , m_{3 } ] , [ n_{4 } , m_{4}])$ ] is a good quadruple if and only if the matrix @xmath167 is generically invertible .    we can write the shape factorization problem of type @xmath144 , [ n_{2 } , m_{2 } ] ,   [ n_{3 } , m_{3 } ] , [ n_{4 } , m_{4}])$ ] in matrix factored form as @xmath170+h[l_4]=[g]$ ] ( see ) , so that @xmath171.\ ] ] this system has a unique solution if and only if @xmath167 is generically invertible , i.e.  invertible for a generic choice of parameter values .
[ ex : quadruple ] suppose that our quadruple @xmath144 , [ n_{2 } , m_{2 } ] ,   [ n_{3 } , m_{3 } ] , [ n_{4 } , m_{4}])$ ] is @xmath172 , [ 2,0 ] , [ 3,0],[2,0])$ ] , which is a special case of joining models of types @xmath55 and @xmath57 in series .
the resulting matrix @xmath167 is the matrix @xmath173    we now determine when this matrix @xmath167 is generically invertible , i.e.  square and full rank .
sylvester matrix _ associated to two polynomials @xmath174 and @xmath175 is the @xmath176 by @xmath176 matrix that has the coefficients of @xmath177 repeated @xmath21 times as columns and the coefficients of @xmath178 repeated @xmath12 times as columns in the following way :    @xmath179    @xmath180    the determinant of the sylvester matrix of the two polynomials @xmath181 and @xmath182 is the _ resultant _ , which is zero if and only if the two polynomials have a common root . in particular , for generic polynomials @xmath181 and @xmath182 , the sylvester matrix is invertible ( * ? ? ? * chapter 3 ) .
we will use the sylvester matrix in the following way .
we will show that there are submatrices of @xmath167 that correspond to the sylvester matrix associated to @xmath88 and @xmath93 .
[ prop : squarefullrank ] if the matrix @xmath167 is square , then it is generically invertible
.    we claim that the columns of @xmath167 can be ordered so that the resulting matrix has the shape @xmath183 where @xmath184 is the sylvester matrix associated to the nonzero coefficients of @xmath107 and @xmath109 .
note that this means that we might shift the coefficients down if necessary so there are no extraneous zero terms of low degree ( i.e.  if the shape is @xmath185 $ ] with @xmath186 ) .
the matrix @xmath187 is a square lower triangular matrix with nonzero entries on the diagonal , and @xmath188 is a square upper triangular matrix with nonzero entries on the diagonal .
this will prove that @xmath167 is invertible , since its determinant will be the product of the determinants of @xmath184 , @xmath187 and @xmath188 , all of which are nonzero .
to prove that claim requires a careful case analysis .
the number of columns of @xmath167 is @xmath189 and the number of rows is @xmath190 . without loss of generality
, we can assume that the maximum is attained by @xmath191 .
we need to distinguish between the two cases where the minimum is attained by @xmath192 and by @xmath193 .
* case 1 : @xmath194 = m_{1 } + m_{4}$ ] . *
since @xmath167 is a square matrix , this implies that @xmath195 . in this case
we group the columns of @xmath167 in the following order .    1 .   the first @xmath196 columns of @xmath165 2 .
then the next @xmath197 columns of @xmath165 3
.   then all @xmath198 columns of @xmath166 4
.   then the remaining @xmath199 columns of @xmath165 .
this choice has the property that the middle two blocks of columns together have the desired form , since we have chosen to start including columns from @xmath165 and @xmath166 precisely when they both have nonzero entries in the same rows , and stopping the formation of these when they stop having nonzero entries in the same rows , which has the correct form .
note we have used all columns of @xmath165 since @xmath200    * case 2 : @xmath194 = m_{2 } + m_{3}$ ] . * note that since @xmath167 is square , this implies that @xmath201 . in this case , we do not need to reorder the columns to obtain the desired form .
we mention how to block the columns to obtain the desired form .    1 .   the first @xmath196 columns of @xmath165 2
.   then the next @xmath197 columns of @xmath165 3 .
then the first @xmath202 columns of @xmath166 4
.   then the remaining @xmath203 columns of @xmath166 .
note that we have the desired number of columns from the second and third blocks , and we have chosen them so that that those columns have nonzero entries at exactly the same rows .
furthermore , we have used all columns of @xmath165 since @xmath204 and all columns of @xmath166 since @xmath205    we can rewrite the matrix in example [ ex : quadruple ] as @xmath206 here the the @xmath207 matrix in the lower righthand corner is the sylvester matrix , the matrix @xmath187 is the @xmath208 matrix in the upper lefthand corner , and the matrix @xmath188 is the empty matrix .
[ proof of theorem [ thm : parameqcoeff ] ] we will show that if the number of parameters equals the number of non - monic coefficients , then the matrix @xmath167 is square . by propositions [ prop : abinvert ] and [ prop :
squarefullrank ] , this will imply that the model is locally identifiable .
let @xmath90 be a spring - dashpot model joined in series from @xmath53 and @xmath54 , where @xmath53 has constitutive equation @xmath138 of shapes @xmath139 $ ] and @xmath140 $ ] , respectively , and @xmath54 has constitutive equation @xmath141 of shapes @xmath142 $ ] and @xmath143 $ ] , respectively . by induction , we can assume that the number of parameters equals the number of non - monic coefficients for the systems @xmath53 and @xmath54 , i.e.  there are @xmath209 parameters in the first and @xmath210 in the second
. assume the number of parameters equals the number of non - monic coefficients in this full system , i.e.   @xmath211 subtracting @xmath212 from both sides , we get that @xmath213 from the definition of @xmath167 , this means the number of rows equals the number of columns , so that @xmath167 is square .    the argument for the parallel extension is identical and is omitted .
[ proof of theorem [ thm : localmt ] ] theorem [ thm : parameqcoeff ] shows that the model is locally identifiable if and only if the number of parameters equals the number of non - monic coefficients .
thus the identifiability properties of the @xmath214 cases in tables [ seriestable ] and [ paralleltable ] are determined by checking if the numbers in the columns corresponding to the number of parameters and the number of non - monic coefficients are equal .
we now determine necessary and sufficient conditions for global identifiability .
let @xmath90 be a spring - dashpot model joined in series from @xmath53 and @xmath54 , where @xmath53 has constitutive equation @xmath138 of shapes @xmath139 $ ] and @xmath140 $ ] , respectively , and @xmath54 has constitutive equation @xmath141 of shapes @xmath142 $ ] and @xmath143 $ ] , respectively .
then the model @xmath59 is globally identifiable if and only if    1 .
@xmath53 and @xmath54 are globally identifiable , 2 .
the shape factorization problem for the quadruple @xmath144 , [ n_{2 } , m_{2 } ] ,   [ n_{3 } , m_{3 } ] , [ n_{4 } , m_{4}])$ ] generically has a unique solution .    similarly ,
if @xmath96 is a spring - dashpot model joined in parallel from @xmath53 and @xmath54 , then @xmath59 is globally identifiable if and only if    1 .
@xmath53 and @xmath54 are globally identifiable , and 2 .   the shape factorization problem for the quadruple @xmath145 , [ n_{1 } , m_{1 } ]
,     [ n_{4 } , m_{4 } ] , [ n_{3 } , m_{3}])$ ] generically has a unique solution .
we handle the case of series extensions , parallel extensions being identical .
let @xmath90 . clearly , @xmath53 and @xmath54 must be globally identifiable otherwise we could give two sets of parameters yielding the same constitutive equation for @xmath53 , which could then be combined with parameters for @xmath54 to get two sets of parameters for @xmath59 yielding the same constitutive equation .
now if the shape factorization problem has a unique solution , there is a unique way to take the constitutive equation for @xmath59 and solve for the constitutive equations for @xmath53 and @xmath54 , since @xmath53 and @xmath54 are globally identifiable , there is a unique way to solve for parameters of those models giving a unique solution for parameters for @xmath59 .
conversely , if there were multiple solutions to the shape factorization problem , then by global identifiability of @xmath53 and @xmath54 , we could solve all the way back to get multiple parameter choices for the same parameter choice for @xmath59 .
note that in our analysis of the shape factorization problem in the previous section , we saw that once @xmath107 and @xmath109 are chosen among all their finitely many values , when the model is locally identifiable there is a unique way to then construct @xmath108 and @xmath110 .
hence , the shape factorization problem has a unique solution when there is a unique way to factor @xmath125 .
this happens if and only if either @xmath215 or @xmath216 , otherwise , generically , we can exchange roots of @xmath107 and @xmath109 giving multiple solutions .
[ cor : globidconstzero ] suppose that @xmath90 is globally identifiable .
then either @xmath53 or @xmath54 must have been one of a spring , a dashpot , or a maxwell model .
suppose that @xmath96 is globally identifiable .
then either @xmath53 or @xmath54 must have been one of a spring , a dashpot , or a voigt model .
the four models given by the spring , dashpot , voigt , and maxwell elements are the only four locally identifiable models that have the property that at least one of the differential operators in its constitutive equation has exactly one term . this can be seen by analyzing the four types ( a , b , c , d ) and looking at all possibilities that arise on combining two equations .
once both operators do not have a single term , no model combined from such a model can have an operator with a single term .
the three choices for the series connection ( spring , a dashpot , or a maxwell model ) are the three of four models that put a differential operator with a single term in the correct place so there could be a unique solution to the shape factorization probelm .
similarly for the parallel connection .
[ proof of theorem [ thm : global ] ] clearly a globally identifiable model is locally identifiable . by corollary [ cor : globidconstzero ] , we must be able to construct such a globally identifiable model by adding at each step either a spring , dashpot , maxwell or voigt element at each step , but when adding a maxwell element it must be used in series and when using a voigt element it must have been added in parallel . however , adding a maxwell element in series can be achieved by adding a spring and then a dashpot both in series .
similarly , adding a voigt element in parallel can be achieved by adding a spring and then a dashpot both in parallel .
hence , we can work only adding springs or dashpots at each step .
adam mahdi was partially supported by the vpr project under nih - nigms grant # 1p50gm094503 - 01a0 sub - award to ncsu .
nicolette meshkat was partially supported by the david and lucille packard foundation .
seth sullivant was partially supported by the david and lucille packard foundation and the us national science foundation ( dms 0954865 ) .
nichols ww ( 2011 ) mcdonald s blood flow in arteries : theoretical , experimental and clinical principles .
london : hodder arnold , 6th ed .
/ wilmer w. nichols , michael f. orourke , charalambos vlachopoulos . edition .


one of the key open questions in the pursuit of a quantum theory of gravity is that of understanding its behaviour in time - dependent scenarios .  in the context of string theory
this has recently led to intensive investigation of asymptotically de sitter spacetimes to see what the prospects are for developing a de sitter / cft correspondence @xcite .
another line of approach has been to construct simple time - dependent solutions that provide ( at least to leading order ) consistent backgrounds for string theory .
this has been carried out in asymptotically flat spacetimes , and recently extended to include the asymptotically anti de sitter ( ads ) case .
this latter situation is of interest since the ads / cft correspondence conjecture could be employed to relate the time - dependence to the behaviour of the non - perturbative field theory dual . in so doing one might be able to avoid ( or at least ameliorate ) difficulties that arise in understanding the behaviour of string theory on time - dependent backgrounds .    in this paper
we extend the search for time - dependent backgrounds to asymptotically locally flat / ads spacetimes .
such spacetimes have a form of magnetic mass referred to as a nut charge .
their euclideanized solutions have been a subject of study in recent years @xcite , in part because in the asymptotically locally ads case they are relevant to 2 + 1 dimensional  `` exotic '' conformal field theories that live on the world volume of m2-branes ( and closely related theories ) , after placing them on squashed three spheres .
we consider here the multiple analytic continuation of such spacetimes to form time - dependent bubble spacetimes with nut charge .
such `` nutty bubbles '' form interesting time - dependent backgrounds that are asymptotically locally flat / ads , extending the variety of time - dependent cases previously studied @xcite .
the stability of these time - dependent bubble solutions has also recently been investigated @xcite .
we consider in this paper the different ( even - dimensional ) taub - nut - ads spacetimes and study their different analytic continuations . after a brief overview of bubble solutions obtained from spherically symmetric spacetimes in section 2
, we go on in section 3 to find three bubble solutions in the four dimensional case , all time - dependent .
two of them are obtained from analytic continuation of taub - nut - ads spacetime and the last one is inspired from the first two bubble solutions . as all of these ` nutty ' bubbles arefour dimensional ,
no issue of classical gravitational instability arises in contrast to the situation considered in  ref.@xcite .
we then compute the boundary stress - energies associated with such solutions in section 4 .  in section 5
we consider double analytic continuations of four dimensional topological nut - charged spacetimes .
we find here some new static bubble - type spacetimes , and an interesting new time - dependent solution .
we find that the presence of a nut charge induces an ergoregion , in which timelike curves must have non - zero momenta in a compact direction . in section 6
we consider double analytic continuations of four dimensional nut - charged spacetimes with a flat subspace and find a new static bubble - type spacetime that is a nut - charged generalization of the ads soliton , as well as an interesting new time - dependent solution : it is asymptotically ads , yet the bubble evolution is milne - like throughout .
similar to topological nut spacetimes , the presence of a nut charge induces an ergoregion in the static case . in section 7
we search for well - behaved bubble solutions with nut charge in more than four dimensions .
our search proves fruitless : we find in more than four dimensions that double ( or multiple ) analytic continuations give us spacetimes with at least two timelike coordinates .
we then study the four dimensional kerr - nut spacetimes @xcite in section 8 and find that their asymptotic behaviour at late times is similar to that of the zero - nut charge case @xcite , and so such spacetimes can be candidate time - dependent backgrounds for string theory .
however the presence of a nut charge can induce quantitatively different ( and sometimes even asymmetric ) particle creation and temporal evolution of the bubble .  in section 9
we consider the ads generalizations of such spacetimes and find that they are beset with the same quantum instability problems as their non - rotating counterparts .
we conclude with some brief comments about generalizations to de sitter spacetimes .
the general form of a four - dimensional static , spherically symmetric metric can be written as @xmath3 where @xmath4 and @xmath5 are functions that are determined by the einstein equations .
a bubble solution is obtained by a double analytic continuation in the time coordinate ( @xmath6 ) and in some other combination of coordinates on the @xmath7-section , effectively changing it into a two - dimensional de sitter spacetime .
three distinct possibilities arise from this procedure   @xmath8 that are characterized by the topology of the two - dimensional de sitter spacetime at every fixed @xmath9 .
locally these are all equivalent , being related to one another by the coordinate transformations @xmath10 globally they are all equivalent provided the coordinates @xmath11 are all unwrapped .
however if any or all of them are identified , then the de sitter spacetime corresponds three distinct evolutions for the circle geometry they describe . for the metric ( [ bubds1 ] ) the de sitter evolution
is that of a circle with exponentially large initial radius as @xmath12 that exponentially shrinks with proper time to a minimal value and then expands again   the worldsheet is like a candlestick . for ( [ bubds2 ] )
the circle begins with zero radius at @xmath13 and then expands exponentially like a bowl . for ( [ bubds3 ] )
the circle begins with infinitesimally small initial radius at @xmath14 and exponentially expands with increasing proper time like a trumpet .
the bubble is located at the largest positive root @xmath15 of @xmath4 , assumed to be at least as great as the largest positive root of @xmath5 ( otherwise there is a wormhole geometry ) .
the coordinate @xmath16 can either be unwrapped or identified with period @xmath17 unless @xmath18 , in which @xmath16 must have period @xmath19 .
the radial variable is then restricted to the range @xmath20 and the solutions ( [ bubds1]-[bubds3 ] ) have a @xmath21 symmetry .
the spacetime is geodesically complete for each .  at any fixed time
, @xmath22 corresponds to the circle of minimal circumference , referred to as the bubble .
this circumference exponentially contracts to @xmath23 after which time it exponentially expands for the solution ( [ bubds1 ] ) .  for ( [ bubds2 ] )
this circumference has zero size at @xmath13 after which it expands exponentially .
this spacetime could matched smoothly at @xmath13 onto any spacetime with a hypersurface of matching extrinsic curvature and 3-metric ( including its time - reversed counterpart ) , and so can be regarded as describing the nucleation of a bubble ( or its collapse if @xmath24 ) .  for ( [ bubds3 ] )
the circumference exponentially expands from infinitesimally small values beginning at @xmath25 .
it is not time - reversal invariant , and so has a counterpart in which the bubble undergoes exponential contraction .
it is straightforward to show that any null curve in spacetimes for which the bubble is expanding has @xmath26 at late times , where @xmath27 is any pair @xmath28 ,  and the overdot refers to a derivative with respect to proper time .
hence observers at different values of @xmath29 will eventually lose causal contact since at late times any null geodesic can only traverse exponentially vanishing changes in @xmath29 .
null geodesics at fixed @xmath29 obey the equation @xmath30 where @xmath31 is the conserved momentum in the @xmath16-direction and @xmath32 is the conserved energy .
( [ nullfixphi ] ) is an equation for a particle in an effective potential @xmath33  @xmath34 , with all allowed motions having @xmath35 .
if @xmath18 and the spacetime is asymptotically flat , then null geodesics oscillate between some minimal and maximal values of @xmath36 , and so by appropriate choices of  @xmath37 observers at any two differing values of @xmath36 can be causally connected .
however if @xmath38 then the effective potential diverges at @xmath22 , and so there will always be some region @xmath39 which is not in causal contact with observers in regions @xmath40 .  if the spacetime is not asymptotically flat then the effective potential may become positive for some @xmath41 , in which case observers in this region are causally disconnected from the region @xmath42 . since @xmath43 it is always possible for observers at different values of  @xmath44 to remain in causal contact at fixed @xmath36 , and so any two observers in causal contact at differing values of @xmath36 will also be in causal contact if they are also at differing values of @xmath16 .
the above results are generalizable to higher - dimensional spacetimes @xcite .  the same classical stability and quantum instability properties hold for all of these bubbles as for the ones analyzed in ref .
the ( 3 + 1)-dimensional taub - nut - ads solution is @xcite : @xmath45 where @xmath46 where @xmath47 is the nonvanishing nut charge , @xmath48 is related to the cosmological constant @xmath49 by the relation @xmath50 .
its euclidean version is found by setting @xmath51 : @xmath52 where @xmath53 the coordinate @xmath16 parametrizes a circle , which is fibered over the sphere parametrized by @xmath54 the geometry of a constant-@xmath36 surface is that of a hopf fibration of @xmath55 over @xmath7 . to avoid singularities the periodicity of @xmath16 must simultaneously be @xmath56 and @xmath57 , which forces one of two conditions : @xmath58 or @xmath59 depending on the value of the largest root @xmath15 of @xmath60 .
the solution ( [ nutmass ] ) is referred to as the nut solution  the fixed point of the @xmath61 isometry is 0-dimensional .
for the special case @xmath62 and @xmath63 , it is pure ads .
solution ( [ boltmass ] ) is referred to as the bolt solution  the fixed point of the @xmath64 isometry is 2-dimensional , as is the case for a black hole .
we can obtain bubble spacetimes by forming several other analytic continuations under the constraint that there be only one timelike direction .
these metrics are all generalizations of ( [ bubds1]-[bubds3 ] ) : @xmath65 @xmath66 @xmath67 and all are exact solutions to the einstein equations .
each is nonsingular and geodesically complete .
for large @xmath36 they all can easily be shown to approach ads@xmath68 if @xmath69 . the geometry of sections of fixed @xmath36 at each constant time slice is that of a  twisted torus @xmath70  using the notation of the previous section , at each constant time slice @xmath71 , the twisted torus has one fundamental circle respectively parametrized by a coordinate @xmath72 the other twisted circle is respectively parameterized by a shifted coordinate @xmath73 .
the torus becomes exponentially more twisted with time . if @xmath74 the geometry of fixed @xmath36 sections is @xmath75ds@xmath76    the quartic function @xmath4 in eq .
( [ ftnlo ] ) has two real roots , one at @xmath77 and one at @xmath78 respectively .
without loss of generality we can locate the bubble at the positive root of @xmath4 ,  since bubble solutions exist for any values of the parameter @xmath79 , positive , zero or negative ( unlike the spherical bubbles of the previous section ) .
the penrose diagrams of the asymptotically flat bubble spacetimes obtained from ( [ tnbub1])-([tnbub3 ] ) are given in figures ( [ fig1])-([fig3 ] ) .
the penrose diagrams of the asymptotically ads bubble spacetimes ( [ tnbub1])-([tnbub3 ] ) are similar to ( [ fig1])-([fig3 ] ) , by changing the diagonal and off - diagonal lines in these figures to a cylinder which encloses the bubble .
the bubble is located on the surface of the shaded region at @xmath80 , where @xmath15  is the root of @xmath81    in the asymptotically flat case ( @xmath82 ) we obtain from ( [ ftnlo ] ) @xmath83 for the location of the bubble .
if @xmath84 the function @xmath4 monotonically approaches unity as @xmath85 , whereas if @xmath86 it has a local maximum
.  in the asymptotically ads case @xmath4 has a local maximum and minimum outside of the bubble if @xmath87 ; otherwise it monotonically increases .
the only restriction on the periodicity of @xmath16 is that it has the value@xmath88 which equals @xmath89 in the asymptotically flat case .     using the notation of the previous section , it is straightforward to check that observers at different values of @xmath29 eventually fall out of causal contact , and that observers at fixed @xmath29 but differing values of @xmath90 can remain in causal contact in both the asymptotically flat and ads cases .
consider next the double analytic continuation @xmath91 applied to the metric ( [ tnlo ] ) : @xmath92 this metric is simply a local patch of one of the solutions ( [ bubds1]-[bubds3 ] ) ; for example under the coordinate transformation @xmath93
\nonumber\end{aligned}\ ] ] the metric ( [ tnbub1 ] ) is recovered .
the lapse function is @xmath94 and so is non - vanishing everywhere except @xmath95 .     in the coordinates ( [ tnb1 ] )
the metric has a killing vector @xmath96 and so describes a static bubble , not evolving in time .
the norm of this killing vector is @xmath97\sin ^{2}\theta \label{timecoeff}\ ] ] and becomes spacelike unless @xmath98 where @xmath99   \label{thetalimits}\ ] ] taking @xmath100 .
these limits define an ` ergocone ' for the spacetime , in which @xmath101 ( @xmath102 ) for all timelike curves in the region @xmath103 ( @xmath104 ) . for @xmath105 ( or @xmath106 )
the angle ( [ thetalimits ] ) becomes vanishingly small at large @xmath36 , and the ergocone shrinks to a pair of lines on the @xmath107 axis where the lapse function vanishes .
observers using the coordinates ( [ tnb1 ] ) can , on a surface fixed @xmath36 ( with respective spacelike / timelike unit normals @xmath108 and @xmath109 ) , associate a quasilocal mass density with this spacetime    @xmath110    where we have used the counterterm method @xcite to compute the mass density .
we therefore obtain a total quasilocal mass @xmath111 of @xmath112 where the integration over @xmath113 is from  @xmath114  to @xmath115  and then from @xmath116  to  @xmath117 *  * we note that if @xmath118 the total mass can be either positive or negative , a situation reminiscent of one that can occur for topological black holes @xcite .     if @xmath119 is finite then this surface penetrates the ergocone , where @xmath120 is restricted to the interval @xmath121 at large @xmath36 .
since the lapse function is positive in this region , we can perform a similar calculation ; the mass density is again given by ( [ massdens ] ) and we obtain    @xmath122    where @xmath15 is the largest positive root of @xmath4 .
somewhat counterintuitively , there is no angular momentum associated with the killing vector @xmath123 in either case , despite the presence of an ergoregion .
it is straightforward to show that both the extrinsic curvature and counterterm contributions to the angular momentum density are identically zero .
interpretation of the preceding results must be taken with some care . since at fixed @xmath124 the metric ( [ tnb1 ] ) describes a torus ,
observers at @xmath125 are separated from observers at @xmath126 by a line of vanishing lapse , and so their surface does not actually enclose the bubble .
however observers on one side of this boundary ( e.g. @xmath127 ) could infer the above results by making use of the @xmath128 symmetry of their metric and assuming that a corresponding set of observers existed on the other side that were making similar calculations / measurements .
the large-@xmath36 boundary stress - energy tensor for the bubble solution ( [ tnbub1 ] ) can be calculated by the well known counterterm method , @xmath129 as one can see , @xmath130 which means that the bubble solution has a negative mass and @xmath131 which means there exists a negative pressure in the @xmath16 direction .
we see that the presence of a nut charge enhances these respective values in the next - to - leading - order term .
the boundary energy momentum tensor components for the bubble solution ( [ tnbub2 ] ) also are given by the same relations ( [ tnb3stress ] ) changing coordinates to @xmath132 and @xmath133 in the last equation of ( [ tnb3stress ] ) .
similarly , the components for the bubble solution ( [ tnbub3 ] ) also are given in ( [ tnb3stress ] ) changing coordinates @xmath134 and sinh@xmath135  in the last equation .
we can compare the above boundary stress - energy tensor components of the bubble solution with the components of boundary energy momentum tensor for the taub - nut - ads spacetime ( [ tnlo ] ) .
the non - vanishing components of the latter are given by : @xmath136 as a consistency check , we see that in special case of @xmath74 for which ( [ tnlo ] ) reduces to the 4-dimensional schwarzschild - ads black hole , the covariant energy - momentum tensor components of  taub - nut - ads spacetime ( [ tnlo ] ) reduce to the same as obtained in @xcite .
the equivalence of the non - vanishing boundary energy - momentum tensor of the bubble solution with the components of boundary energy momentum tensor for the taub - nut - ads spacetime ( [ tnlo ] ) is a direct consequence of analytic continuation .
so the boundary cft quantities related to the energy - momentum tensor ( like the central charge or casimir energy ) are the same for the two nutty bubble spacetimes and the taub - nut - ads spacetime .
such an equivalence was observed for the boundary energy - momentum tensors of  bubble solutions derived from analytic continuation of the five - dimensional schwarzschild - ads metric and schwarzschild - ads metric @xcite .
although the boundary theory is not known for the bubble metrics ( [ tnbub1])-([tnbub3 ] ) ( in contrast to the 5d - schwarzschild bubble @xcite , for which the boundary theory is @xmath137 sym living on @xmath75ds@xmath138 ) , since the boundary is conformally flat we can use the following relation for the stress tensor @xcite : @xmath139 where @xmath140and @xmath141 are conserved quantities constructed from the curvatures via @xmath142 and @xmath143 is the traceless part .
for the boundary metric of the bubble ( [ tnbub1 ] ) , the functions @xmath144and @xmath141 are : @xmath145 @xmath146 by fixing the constants @xmath147 and @xmath148 from the boundary theory , the result of ( [ s ] ) is in agreement with ( [ tnb3stress ] ) in the limit @xmath149 , since the two constants must vanish due to absence of the conformal anomaly in the odd dimensional boundary theory . consequently the above comparison of  the stress tensor ( [ tnb3stress ] ) to ( [ s ] ) does not result in any non - trivial connection between them .
in other words , in the limit @xmath85  , all of the stress tensor components ( [ tnb3stress ] ) vanish and since the boundary conformal anomaly vanishes in odd dimensions , so the boundary stress tensor ( [ s ] ) also vanishes .
these bubbles have in four - dimensional euclidean space the general form @xmath150 where @xmath151 there are three distinct analytic continuations : @xmath152 @xmath153 @xmath154 and these all satisfy the einstein equations .
the @xmath113-coordinate is no longer periodic and takes on all real values .
these metrics can , after a coordinate transformation be written in the form @xmath155 where @xmath156 for the respective cases ( [ top1]-[top3 ] ) .
they are all locally equivalent to each other under a coordinate transformation , and so it is tempting to conclude that they are globally equivalent since the @xmath157-coordinate is not periodic .
however at fixed @xmath158 the @xmath159 section describes ( for @xmath106 ) a two - dimensional anti de sitter spacetime .
this spacetime can have non - trivial identifications @xcite and so can describe either a two - dimensional black hole @xcite as in eqs .
( [ top1],[top3 ] ) ( the latter being a special case of a massless extremal black hole ) or else standard ads@xmath160 , as in eq .
( [ top2 ] ) .
another clear difference between this case and the previous cases is that the quartic function @xmath161 in ( [ ftopnut ] ) can have four real roots , which we label @xmath162 in increasing order of size . from this perspective
these metrics can represent two different situations .
first , for @xmath163  ( or @xmath164 ) , these spacetime represent static bubbles .
the coordinate @xmath36 has the range @xmath165 ( or alternatively @xmath166 ) , and the @xmath113-coordinate takes on all possible values .  at large @xmath167 it is straightforward to show that these metrics are asymptotically ads .  both invariants @xmath168 and @xmath169 are finite everywhere and the spacetime has no singularities .  since @xmath170 is not periodic , the bubble does not close back on itself  in this sense it is not really a bubble , but it should be a valid background for string theory .  in this case
each of the metrics contains an ergoregion given by the expression @xmath171  the structure of the ergoregion depends upon the value of @xmath172 and @xmath173 .
consider first the metric ( [ top3 ] ) , ( @xmath174 ) for which the ergoregion is independent of @xmath170 .
if @xmath175 , then there will be an ergoregion everywhere except for a very narrow strip of infinite length near the bubble where @xmath176 , with @xmath177 .
the strip slightly broadens as @xmath79 increases , as can be seen from figure ( [ fig31 ] ) , where the function @xmath178versus @xmath36 is plotted .
as @xmath179 from above @xmath180 gets larger , diverging at @xmath181 , at which the ergoregion vanishes .
once @xmath182 there is no ergoregion .
for the metric ( [ top2 ] ) , with @xmath183 , there will be no ergoregion for @xmath184 . for @xmath185 @xmath186
, the ergoregion is confined to the right half of the @xmath187 plane bounded by a curve with a vertical asymptote at @xmath188 , and a horizontal asymptote at @xmath189 .
this curve has a mirror image about @xmath190 and the ergoregion for @xmath191 has a mirror image for negative @xmath170 .
a typical ergoregion is plotted in figure ( [ fig32 ] ) for special values of @xmath192 in this figure the vertical axis is @xmath170 and horizontal axis is @xmath193     for the metric ( [ top1 ] ) , with @xmath194 , the ergoregion is confined to a strip outside the event horizon at @xmath195 bounded by a curve and its mirror image in @xmath170 that asymptote for large @xmath36 at @xmath196 for @xmath182 .
this curve smoothly intersects the location of the bubble at the event horizon .
a typical ergoregion is plotted in figure ( [ fig33 ] ) for special values of @xmath197    as @xmath179 from below , the strip becomes increasingly wide , and at @xmath198 the boundary curve asymptotes to @xmath199 . for @xmath185 @xmath186 , the boundary of the ergoregion has a vertical asymptote at @xmath188 , and covers the remainder of the @xmath187 plane for larger values of @xmath36 .
if @xmath200 has 4 real roots ,  a second possibility emerges : the coordinate @xmath36 can have the range @xmath201 , and can be transformed to a periodic variable .
removal of conical singularities at @xmath202 is satisfied by requiring @xmath203 , which is @xmath204 or alternatively , if @xmath205 @xmath206   \label{toproots2}\ ] ] since for the function @xmath200 we must have @xmath207 .
( [ toproots2 ] ) has no solutions since its right hand side is always negative .
however if @xmath208 then conical singularities can be removed .  if @xmath209 then @xmath62 , whereas if @xmath210 we can attempt to obtain a new spacetime by setting @xmath211 , @xmath212@xmath213 and then taking the limit as @xmath214 .
this is possible only if @xmath106 , since a rescaling of the coordinates @xmath215 by factors of @xmath216 yields a singularity in the metric .
however if @xmath106 one of the roots becomes @xmath217 at which a singularity is located .     if @xmath62 then @xmath218 in order that @xmath200 be real .
the metric can be written as *  * @xmath219 where @xmath220 $ ] , @xmath16 has the period @xmath221 , and @xmath222 strictly speaking this is not a bubble spacetime : instead the @xmath223 section is the surface of a deformed 2-sphere .  the ergoregion is now located in a band about the equator at @xmath224 , whose specific boundary curves are complicated functions of @xmath172 .  if @xmath174 this band is very narrow and exists only for a very small range of @xmath47 , and it remains constant everywhere .
however if @xmath225 the band will change shape as @xmath170 increases .  note that if we set @xmath226 then we must have at least one conical singularity , since the condition  @xmath227 ca nt be satisfied unless @xmath106 , in which case there is again a singularity . the black hole ( @xmath194 ) becomes a compact region , surrounded by a deformed ` ergosphere ' for small values of @xmath170 .
finally , we can continue to the metric @xmath228 which locally can be transformed to one of the metrics ( [ top1]-[top3 ] ) via a change of coordinates ; it can be regarded as the extension of the @xmath229 case if @xmath230 is not periodic .
however if the @xmath230-coordinate in ( [ topnutperiodic ] ) is identified then these metrics are not globally equivalent , yielding new bubble solutions .
if @xmath231  ( or @xmath164 ) , the metric ( [ topnutperiodic ] ) describe a bubble spacetime in which the bubble expands from zero size to the finite size @xmath232 and then contracts back to zero size again .
all spacetime events are causally connected to one another , and near the initial expansion ( or final contraction ) the scale factor is linear in time and the spacetime behaves like a milne universe .     if @xmath201 then the metric ( [ topnutperiodic ] ) describes a compact time - dependent 3-geometry of a twisted circle over a 2-sphere . from the previous analysis ,
conical singularities are avoided only if either @xmath62 or @xmath106 . in each case
the spacetime oscillates , with the @xmath55 represented by the @xmath230 coordinate repeatedly expanding and contracting .
there are again no event horizons .
these spacetimes are not really bubbles because the @xmath9  section is a compact space .
however they are candidate time - dependent backgrounds . since all of the topological backgrounds ( [ top1])-([top3 ] ) and ( [ topnutperiodic ] ) are  four dimensional they are stable , similar to the other nutty bubbles ( [ tnbub1])-([tnbub3 ] ) .
we close this section by noticing that in the case of  ( [ topnutperiodic ] ) , null geodesics relevant for s - waves behaviour satisfy : @xmath233 the integration involves the square root of a quartic function and so finding a closed - form relation between coordinates @xmath36  and @xmath157  for a null geodesic is quite problematic .
the relation between these two coordinates is essential to obtain a relation between phase differences of  two geodesics on @xmath234 and @xmath235 .
this relation can be used as an evidence for particle creation .
we shall not consider this complicated problem here .
however in the case of kerr - nut bubbles , the integration can be done and a simple relation is obtained for null geodesics , eq .
( [ trnull ] ) , as we shall later consider .
another class of nut - charged spacetimes has the metric @xmath236 where @xmath237 has only two real roots , and the @xmath238 section describes a flat two - dimensional space , possibly with identifications .
there are two distinct double analytic continuations : @xmath239 where @xmath240 .
the @xmath241 section is now a two - dimensional milne / rindler spacetime , which can also be written as @xmath242 which is a nut - charged generalization of the ads soliton @xcite ; we shall continue to refer to it as a bubble . without loss of generality the bubble ( [ nutsoliton ] ) can be located at the positive root of @xmath243 since @xmath79 can take on all real values .
the cases @xmath240 are locally equivalent under a coordinate transformation , but will be globally distinct if @xmath230 is periodically identified .
when @xmath244 the bubble undergoes a milne - type evolution . at any fixed @xmath245
the geometry of the bubble is that of a twisted torus , with the coordinate @xmath230 parametrizing one circle , and with the other twisted circle parameterized by the shifted coordinate @xmath246 .
*  *    turning next to the causal structure of the bubble ( [ nutrind ] ) , we consider two observers at different coordinates @xmath247  in the background of bubble ( [ nutrind ] ) , a null curve satisfies @xmath248  and so @xmath249 as @xmath250 , we have @xmath251 , and so observers at two points with different @xmath230  can communicate to each other for arbitrarily large @xmath252 .
the penrose diagram of the bubble spacetime ( [ nutrind ] ) with @xmath244 is given in figure ( [ fig34 ] ) .
any radial null ray leaving the surface of the bubble reaches infinity in a finite amount of coordinate time @xmath253 , and there is always an infinite amount of proper distance between the bubble and spatial infinity.the bubble is located on the surface of the shaded region at @xmath22 , where @xmath15  is the root of @xmath254    when @xmath194 the geometry is that of a static bubble with an ergoregion .
writing the metric ( [ nutrind ] ) for this case as @xmath255 we see that the lapse function is @xmath256 and so vanishes at the rindler horizon @xmath190 .
the ergoregion is described by @xmath257 and at large @xmath36 includes the entire @xmath187 plane except for a strip bounded by a pair of curves opposite the @xmath36-axis that asymptote to @xmath258 .
as @xmath36 decreases , this strip widens , becoming infinitely wide at the bubble edge .
the ergoregion is similar to that given in fig .
( [ fig32 ] ) , except that the vertical asymptotes are at the bubble edge @xmath259 .
using the same coordinates as ( [ nutrinderg ] ) , we have @xmath260 and so there is no restriction as to the maximum change of coordinate @xmath170 for points on the null curve as @xmath261 .
so observers at any two points with different values of @xmath170  can communicate with each other .
in this section we consider the possibility of extending our nut - charged bubble solutions to higher dimensions .
the six dimensional taub - nut - ads is a fibration over @xmath262 @xcite : @xmath263 if @xmath264 is the largest root of @xmath265 , then singularities in ( [ tneuc ] ) are avoided by choosing the mass parameter and the function @xmath265 by    @xmath266    and    @xmath267    the coordinate @xmath268 must be periodic with the period @xmath269    in performing a single analytic continuation from the euclidean metric ( [ tneuc ] ) we are forced to introduce two timelike directions unless we analytically continue @xmath16 .
the reason is that there is a single nut charge @xmath270 and all analytic continuations in one of the @xmath7 sections forces @xmath271 .
this in turn forces an additional analytic continuation in the other @xmath7 section , yielding two timelike directions .     for example performing the following analytic continuation @xmath272 on the metric ( [ tneuc ] ) yields    @xmath273    where    @xmath274   \label{vtnb1}\ ] ]    to find the structure of the bubble , we set @xmath275 where @xmath276 then the induced metric on the bubble is    @xmath277   \label{nut1}\ ] ]    note that the angle @xmath278 will be restricted to within a certain range in order for @xmath279 to be a timelike killing vector .
another continuation is    @xmath280    and this gives    @xmath281    with the induced bubble metric :    @xmath282   \label{nut3}\ ] ]    similarly , the following continuation :    @xmath283    gives    @xmath284    with the induced bubble metric :    @xmath285 \label{nut4}\ ] ]    which again has two timelike directions .
similar considerations apply for all higher - dimensional taub - nut - ads solutions have @xmath286 fibrations over a product of @xmath7 s .
for example the 8-dimensional solutions has a @xmath286 fibration over @xmath287 @xcite : @xmath288 where @xmath289 in this case , the mass parameter @xmath79 takes the value @xmath290 to have a nut solution .
the coordinate @xmath268 has a period @xmath291 , to avoid a conical singularity at @xmath292 similarly , the ten dimensional taub - nut - ads solution is a result of a @xmath286 fibration over @xmath293 @xmath294 @xmath295 where for a nut solution , the mass parameter and function @xmath4 are given by : @xmath296 and @xmath297    the coordinate @xmath268 must has a period @xmath298 , to avoid a conical singularity at @xmath292    the form for a general fibration over @xmath299 is the same .
all analytic continuations in the @xmath7 sections force the introduction of more than one timelike direction because the nut parameter must be analytically continued to avoid a complex metric .
another possibility is a over @xmath300 . in six dimensions
this is @xcite : @xmath301 where @xmath147 is given by @xmath302 and @xmath303 is the metric over the @xmath300 , which is given by @xmath304    the function @xmath265 in the metric is the same as function ( [ vtneuc ] ) , appearing in the taub - nut - ads fibered over @xmath305 in the limit @xmath306 the metric ( [ tn6cp2 ] ) approaches the metric considered in @xcite .
the situation here is somewhat different in that we have not found any analytic continuation of the coordinates of the @xmath307 section that yields a single timelike direction .
further higher - dimensional generalizations are products of @xmath7 and @xmath308 and the same considerations apply , namely more than one timelike direction appears after analytic continuation unless @xmath16 is the coordinate continued .
we consider next rotating nut solutions in four dimensions .
such solutions are given by the class of kerr - nut metrics , which have the ( euclidean ) form @xcite @xmath309^{2}+\sin ^{2}\theta \lbrack ad\chi -(r^{2}-n^{2}-a^{2})d\phi ] ^{2}}{[r^{2}-(n+a\cos \theta ) ^{2 } ] } \\
+ [ r^{2}-(n+a\cos \theta ) ^{2}](\frac{dr^{2}}{v(r)}+d\theta ^{2 } ) \end{array } \label{kn}\ ] ] where @xmath310 applying the following analytic continuations    @xmath311    on the metric ( [ kn ] ) yields : @xmath312^{2}+\cosh ^{2}t[ad\chi
-(r^{2}+n^{2}-a^{2})d\phi ] ^{2}}{[r^{2}+(n - a\sinh t)^{2 } ] } \\   + [ r^{2}+(n - a\sinh t)^{2}](\frac{dr^{2}}{\widetilde{v}(r)}-dt^{2 } ) \end{array } \label{kn1}\ ] ] where the function @xmath313 @xmath314is given by : @xmath315 the largest root of @xmath313 is @xmath316 and so the induced metric on the bubble is @xmath317dt^{2}+\frac{\cosh ^{2}td% \tilde{\chi}^{2}}{[\widetilde{r}_{+}^{2}+(n - a\sinh t)^{2 } ] }   \label{kn1i}\ ] ] where @xmath318 .
for all values of @xmath319 the coefficient of time is negative definite .
so , unlike the case of kerr - ads @xcite , the induced bubble metric obtained through the analytic continuation of the kerr - nut metric in four dimension has a well defined behaviour . in the special case of @xmath320 ( [ kn1i ] )
reduces to the kerr bubble obtained in @xcite . for @xmath321
the bubble evolves like ds@xmath160 .
this situation is similar to the case of schwarzschild and kerr bubbles @xcite .
however the presence of a nut charge causes the time evolution to become asymmetric , with the bubble approaching a minimal radius and then later a maximal radius near @xmath322 ( or the time - reversal of this , depending on the signs of the parameters ) . for very large
@xmath319 the bubble metric becomes @xmath323 which is independent of the nut charge . in this limit
, the circle parametrized by @xmath324 remains at a fixed radius , so we do not have any time evolution for the bubble .
so when the bubble evolves like ds@xmath325 the effect of nut charge is to shift the value of the coordinate @xmath157 to satisfy the above mentioned inequality . on the other hand , for large values of @xmath319
the nut charge does not have any physical effect on the induced bubble metric ( [ kn1i ] ) , and the spacetime evolves like a milne universe at late times .
the  penrose diagram of the kerr - nut bubble ( [ kn1i ] ) is given in figure ( [ fig4 ] ) .
for the four dimensional kerr - nut bubble ( [ kn1 ] ) , the null geodesics relevant for the s - waves satisfy : @xmath326 integrating the above equation , we find the relation @xmath327 where @xmath147  is a constant labeling the different null geodesics and @xmath328refers to outgoing and incoming modes .
an s - wave at @xmath234  is given by a purely positive frequency @xmath329 , where @xmath330  here @xmath331  and @xmath332  are the minkowski coordinates in the time and radial directions at infinity . for this s - wave , the phase difference between two geodesics at  @xmath333  and @xmath334  on @xmath335
i.e. @xmath336can be related to the phase difference traced back to @xmath337  and @xmath338on @xmath339  where @xmath340  the bubble metric in the @xmath341section for very large @xmath36
( such that @xmath342 ) and very large @xmath319  ( @xmath343 ) , is given by    @xmath344    where in terms of the original coordinates @xmath345 using relation ( [ trnull ] ) , we find that the future coordinate @xmath346  of the null geodesic is @xmath347  similarly , the past coordinate @xmath348  is @xmath349 as with the @xmath106 case @xcite , we can find @xmath350  between two geodesics .
we solve the equation ( [ vv ] )  for @xmath351 in terms of @xmath352  ( taking geodesic @xmath353  as the reference geodesic ) and then substitute the result in eq.([uu ] ) .
this gives @xmath354 indicating that a positive frequency s - wave at future can not continue back to a purely positive frequency wave at past and so this means that we have a particle creation in bubble spacetime ( [ kn1 ] ) . from the above
obtained result , one can calculate the bogoliubov coefficients ( which determines the extent of particle creation ) , the total number of particles and other qualitative results about the particle creation .
another analytic continuation of the metric ( [ kn ] ) is obtained through @xmath355 which yields : @xmath356^{2}+\sinh ^{2}t[ad\chi -(r^{2}+n^{2}+a^{2})d\phi ] ^{2}}{[r^{2}+(n+a\cosh t)^{2 } ] } \\   + [ r^{2}+(n+a\cosh t)^{2}](\frac{dr^{2}}{v^{\ast } ( r)}-dt^{2 } )
\end{array } \label{kn2}\ ] ] where the function @xmath357 is given by : @xmath358 in this case @xmath359 @xmath360 , the induced metric on the bubble is given by : @xmath361dt^{2}+\frac{\sinh ^{2}td% \widehat{\chi } ^{2}}{[r_{+}^{\ast 2}+(n+a\cosh t)^{2 } ] }   \label{kn2i}\ ] ] where @xmath362 so again we have a bubble solution . for @xmath363
the bubble evolves like ds@xmath160 , now with symmetric time evolution .
however in this case there are two possible time - evolution scenarios .
if @xmath47 and @xmath364 are both of the same sign , the bubble has only one minimal vanishing radius .
if @xmath47 and @xmath365 have opposite sign then the bubble will expand to a maximum , shrink to zero , expand again to the same maximum before shrinking once again .
for large @xmath319 the bubble metric becomes @xmath366 again , independent of the nut charge .
hence in this limit , the circle parametrized by @xmath367 remains at a fixed radius , so we do not have any time evolution for the bubble . as with the last case , when the bubble evolves like ds@xmath325 the effect of the nut charge is to shift
the value of coordinate @xmath157 to satisfy the above mentioned inequality .
on the other hand , for large values of coordinate @xmath319 the nut charge does not have any physical effect on the induced bubble metric ( [ kn2i ] ) , whose late - time evolution is also milne - like .
a typical penrose diagram of the kerr - nut bubble ( [ kn2i ] ) , when @xmath47  and @xmath365  have the same sign , is given by figure ( [ fig5 ] ) .
also , by considering an s - wave in the bubble background ( [ kn2 ] ) , we can show that particle creation can occur .
we can describe the causal structure of kerr - nut bubbles ( [ kn1 ] ) and ( [ kn2 ] ) as follows .
if we consider two observers at different coordinates @xmath368  they are causally connected even for times @xmath369  to see this , we consider a null curve which satisfies @xmath370^{2}+\sinh ^{2}t[a\stackrel{\cdot } { \chi } -(r^{2}+n^{2}+a^{2})% \stackrel{\cdot } { \phi } ] ^{2}}{[r^{2}+(n+a\cosh t)^{2 } ] } & &   \nonumber \\ + [ r^{2}+(n+a\cosh t)^{2}](\frac{\stackrel{\cdot } { r}^{2}}{v_{0}(r)}-% \stackrel{\cdot } { t}^{2 } ) & = & 0\end{aligned}\ ] ] where @xmath371 or @xmath357 and the overdot represents derivative with respect to an affine parameter labeling the curve . since we assumed @xmath372  we conclude that @xmath373 hence , for points on the null curve , when time coordinate @xmath374 , we do not have any restriction on the maximum change in the coordinate @xmath36 .
so , two observers at two points with different values of @xmath36  can communicate with each other .
the two killing vectors @xmath375 and @xmath376 yield two constants of the motion @xmath377 we can use the equation ( [ com ] ) to eliminate @xmath378 and @xmath379 from the null geodesics .
in this case , for a fixed time @xmath319 the null geodesic equation reduces to that of motion of a particle in an effective potential :    @xmath380    where the effective potential is : @xmath381 so geodesics exist in the region for which @xmath382 for example , on a constant time slice @xmath383 if we set @xmath384 the effective potential is plotted in figure ( [ fig6 ] ) .
it is infinite at @xmath385 for @xmath386 the potential is positive definite and so a geodesic can not penetrate into this region .
the four dimensional euclidean kerr - nut - ads spacetime is given by @xcite :
@xmath387^{2}+h(\theta ) \sin ^{2}\theta \lbrack ad\chi -(r^{2}-n^{2}-a^{2})d\phi ] ^{2}}{\sigma ^{4}[r^{2}-(n+a\cos \theta ) ^{2 } ] } \\   + [ r^{2}-(n+a\cos \theta ) ^{2}](\frac{dr^{2}}{v(r)}+\frac{d\theta ^{2}}{% h(\theta ) } ) \end{array } \label{kerrnut}\ ] ] where the functions @xmath265 and @xmath388 are given by the following expressions : @xmath389 and @xmath390 there is a bolt solution @xmath391 given by @xmath392 , with @xmath393 ; noanalogue of the pure nut solution exists for this case @xcite . to avoid conical singularities in different sections of the spacetime
, we must impose some constraints .
for example , avoiding conical singularities in the @xmath394 section implies @xmath395 the period of the coordinate @xmath16 is @xmath396 where the left and right hand sides of the equality come from avoiding conical singularities in the @xmath397 section and ensuring regularity of the metric along the @xmath398-axis . in the previous equality
, @xmath399 and @xmath400 is the real root of @xmath401 applying the following analytic continuations    @xmath402    on the metric ( [ kerrnut ] ) yields : @xmath403^{2}+% \widetilde{h}(t)\cosh ^{2}t[ad\chi -(r^{2}+n^{2}-a^{2})d\phi ] ^{2}}{\sigma ^{4}[r^{2}+(n - a\sinh t)^{2 } ] } \\   + [ r^{2}+(n - a\sinh t)^{2}](\frac{dr^{2}}{\widetilde{v}(r)}-\frac{dt^{2}}{% \widetilde{h}(t ) } ) \end{array } \label{kerrnut1}\ ] ] where the functions @xmath313 and @xmath404 are given by : @xmath405 and @xmath406 at @xmath407 ( where @xmath408 ) , and for the coordinate @xmath157 in the range of @xmath409 @xmath410 the function @xmath411 is positive definite , and the induced metric on the bubble is given by : @xmath317\frac{dt^{2}}{\widetilde{% h}(t)}+\frac{\widetilde{h}(t)\cosh ^{2}t[ad\chi -(\widetilde{r}% _ { + } ^{2}+n^{2}-a^{2})d\phi ] ^{2}}{\sigma ^{4}[\widetilde{r}% _ { + } ^{2}+(n - a\sinh t)^{2 } ] }   \label{kerrnut1i}\ ] ] where the critical coordinate @xmath412 is : @xmath413 on the other hand , for the range @xmath414 the function @xmath404 is negative definite and so the induced metric ( [ kerrnut1i ] ) does not describe a well defined bubble .
in fact , the situation is similar to the induced bubble metrics obtained through analytic continuation of the kerr - ads metric in four dimensions @xcite . in that case
the induced metric also changes signature at some critical time .
the result of @xcite is a special case of our result when @xmath106 .
another special case of the metric ( [ kerrnut ] ) is when @xmath415 in this case , the functions @xmath416 and from relation ( [ thetacriticalkerrnut ] ) , we have @xmath417 so we have a well behaved induced bubble metric @xmath418   \label{nuti}\ ] ] where we derived this form of the induced metric in section one .    the other analytic continuation of the metric ( [ kerrnut ] ) is obtained through @xmath419 which yields @xmath356^{2}+h^{\ast } ( t)\sinh ^{2}t[ad\chi -(r^{2}+n^{2}+a^{2})d\phi ] ^{2}}{\sigma ^{\ast 4}[r^{2}+(n+a\cosh t)^{2 } ] } \\   + [ r^{2}+(n+a\cosh t)^{2}](\frac{dr^{2}}{v^{\ast } ( r)}-\frac{dt^{2}}{h^{\ast } ( t ) } ) \end{array } \label{kerrnut2}\ ] ] where the functions @xmath357 and @xmath420 are given by @xmath421 and @xmath422 and @xmath423 in this case at point @xmath359 ( for some values of parameters @xmath424 , we have @xmath425 ) , and for the coordinate @xmath157 in the range of @xmath409 @xmath426 the function @xmath420 is positive definite , and the induced metric on the bubble is given by : @xmath361\frac{dt^{2}}{h^{\ast } ( t)}+% \frac{h^{\ast } ( t)\sinh ^{2}t[ad\chi -(r_{+}^{\ast 2}+n^{2}+a^{2})d\phi ] ^{2}% } { \sigma ^{\ast 4}[r_{+}^{\ast 2}+(n+a\cosh t)^{2 } ] }   \label{kerrnut2i}\ ] ] where the critical coordinate @xmath427 is now    @xmath428    on the other hand , for the range @xmath429 the function @xmath430 is negative definite and so the induced metric ( [ kerrnut2i ] ) does not describe a well defined bubble .
our consideration of bubble spacetimes with nut charge has yielded several new candidate four - dimensional time - dependent backgrounds for string theory whose properties merit further exploration .
we have found that the presence of a nut charge can considerably modify the evolution of the spacetime , particularly in the topological and rotating cases .
the nut charge induces the presence of an ergoregion in the latter case , although for the spacetimes in which this occurs the bubble is not compact .
the bubble ( [ nutrind ] ) is particularly intriguing , since its evolution is milne - like throughout . in higher dimensions
the presence of a nut charge appears to forbid the existence of any bubble spacetime , since the double analytic continuation yields more than one timelike direction .
we close with a few comments about de sitter spacetime .
all of the above results in the ads case can , by the simple continuation @xmath431 be changed to taub - nut - ds bubbles .
for example the tn - de sitter space can be obtained from ( [ tnlo ] ) , and is given by :      where @xmath433 these spacetimes are expected to be better behaved , at least in the rotating case , since the functions @xmath434 and @xmath420 in eqs .
( [ kerrnut1 ] ) and ( [ hkerrnut2 ] ) respectively need never vanish if the parameters are appropriately chosen .
we leave a study of this case for future work .        99 a. strominger , jhep * 0110 * ( 2001 ) 034 , _ hep - th/0106113 _ ; e. witten , _ hep - th/0106109 _ ; m. li , jhep * 0204 * ( 2002 ) 005 , _ hep - th/0106148 _ ; d. klemm , _ nucl.phys . _ * b625 * ( 2002 ) 295 , _ hep - th/0106247 _ ; v. balasubramanian , p. horava and d. minic , jhep * 0105 * ( 2001 ) 043 , _ hep - th/0103171 _ ; a. chamblin and n. d. lambert , _ phys .
* b508 * ( 2001 ) 369 , _ hep - th/0107031 _ ; y. gao , _ hep - th/0107067 _ ; j. bros , h. epstein and u. moschella , _ phys .
* d65 * ( 2002 ) 084012 , _ hep - th/0107091 _ ; s. nojiri and s. d. odintsov , jhep * 0112 * ( 2001 ) 033 , _ hep - th/0107134 _ ; e. halyo ,   _ hep - th/0107169 _ ; i. sachs and s. n. solodukhin , _ phys .
rev . _ * d64 * ( 2001 ) 124023 , _
hep - th/0107173 _ ; a. j. tolley and n. turok , _ hep - th/0108119 _ ; s. nojiri and s. d. odintsov , _ phys . lett . _ * b519 * ( 2001 ) 145 , _ hep - th/0106191 _ ; c. m. hull , jhep * 0111 * ( 2001 ) 012 , _ hep - th/0109213 _ ; s. cacciatori and d. klemm , _ class .
* 19 * ( 2002 ) 579 , _ hep - th/0110031 _ ; b. mcinnes , _ nucl . phys_. * b627 * ( 2002 ) 311 , _ hep - th/0110062 _ ; s. nojiri and s. d. odintsov , _ phys . lett . _ * b523 * ( 2001 ) 165 , _ hep - th/0110064 _ ; a. strominger , jhep * 0111 * ( 2001 ) 049 , _ hep - th/0110087 _ ; v. balasubramanian , j. de boer and d. minic , _ phys . rev . _ * d65 * ( 2002 ) 123508 , _ hep - th/0110108 _ ; y. s. myung , _ mod .
* a16 * ( 2001 ) 2353 , _ hep - th/0110123 _ ; b. g. carneiro da cunha , _ phys . rev .
_ * d65 * ( 2002 ) 104025 , _ hep - th/0110169 _ ; r. cai , y. s. myung , y. zhang , _ phys .
_ * d65 * ( 2002 ) 084019 , _ hep - th/0110234 _ ; u. h. danielsson , jhep * 0203 * ( 2002 ) 020 , _ hep - th/0110265 _ ; s. ogushi ,   _ mod .
lett . _ * a17 * ( 2002 ) 51 , _ hep - th/0111008 _ ; a. c. petkou and g. siopsis , jhep *  0202 * ( 2002 ) 045 , _ hep - th/0111085 _ ; r. cai , _ phys .
_ * b525 * ( 2002 ) 331 , _ hep - th/0111093 _ ; t. shiromizu , d. ida and t. torii , jhep * 11 * ( 2001 ) 010 , _ hep - th/0109057 _ ; a. m. ghezelbash and r. b. mann , jhep * 0201 * ( 2002 ) 005 , _ hep - th/0111217 _ ; m. h. dehghani ,  _ phys .
_ * d65 * ( 2002 ) 104003 , _ hep - th/0112002_.