thin shell components can be found in many large built - up mechanical structures such as cars , ships , and aeroplanes .
the prediction of the mid- and high- frequency vibrational properties of these structures becomes computationally prohibitive for standard element - based methods , such as the finite element method @xcite .
the main reasons for this limitation are : firstly , very fine meshes are required for an adequate representation of the highly oscillatory wave solutions and the computational complexity grows with frequency raised to the power of the dimension of the space being modelled . secondly , small uncertainties originating from the manufacturing process lead to a much larger variability in vibro - acoustic responses in the high frequency range @xcite , meaning that the response of any individual manufactured structure is of less interest to computer aided engineering practitioners than the average responses .    methods such as statistical energy analysis ( sea ) @xcite and ray tracing @xcite are more commonly applied for modelling wave problems at high frequencies .
sea has traditionally proved more popular for structural vibration problems with low damping @xcite , whereas ray tracing has found its niche in applications where relatively few reflections need to be tracked in computer graphics @xcite , room acoustics @xcite and seismology @xcite .
ray tracing has also been applied to elastic wave transmission problems on shells and plates @xcite . in this context
the ring frequency , that is the frequency above which longitudinal waves in a curved shell behave as they would in a flat plate , provides a useful point of reference . beyond the ring frequency
the ray dynamics in a curved shell are relatively simple following geodesic paths , but below the ring frequency one finds that asymptotic ( ray ) theories show richer features .
the dispersion relations become highly anisotropic with likewise anisotropic propagation @xcite .
an sea treatment would fail to capture the non - trivial way in which the curved shell geometry influences the wave and ray propagation below the ring frequency @xcite , and hence ray methods can provide useful insight @xcite .
wave scattering from discontinuous line joints below the ring frequency has been considered in @xcite .
however , manufacturers of large built - up mechanical structures are increasingly developing larger and lighter sub - components , whereby large thin shell structures are replacing more traditional plate - beam and multi - plate assemblies .
the manufacturing process for such thin shells often entails casting molten metal ( for example aluminium ) , which gives rise to curved components with smooth transitions between flat and curved regions .
this raises the question of how the ray and wave scattering , and hence the vibrational properties of the structure , are modified in these smooth designs . in this work
we study the case of a singly curved shell , chosen here as an assembly of two plates joined smoothly with a quarter of a cylinder .
this simplified assembly represents a typical curved region within one of the larger sub - structures described above .
we shall go beyond plane waves and ray tracing calculations by solving the full wave scattering problem numerically .
the numerical solution to the full wave scattering problem will be discussed in comparison with the corresponding ray tracing calculations . in both cases
we find effective laws for the scattering properties , which may be inserted into ray or wave propagation modelling techniques such as dynamical energy analysis ( dea ) @xcite or the wave and finite element method @xcite .
in particular , combining these local scattering models within a larger model of a built - up structure will lead to a hybrid method for structures including curved thin shell components in the mid - frequency regime . here
, the natural definition of the mid - frequency regime is given by the range of frequencies that are high enough for a pure fem analysis to be impractical , but low enough so that a simple geodesic description of the trajectory evolution is invalid . in the high - frequency case ,
dea @xcite can be applied to model the vibrational energy transport of a built - up structure along geodesic paths using the mesh data from a fem analysis @xcite .
in fact , dea presents a link between ray tracing and sea by casting the wave or ray problem into the language of evolution operators .
we note that an equivalent operator formalism has also been long known in computer graphics @xcite , although the theory underlying dea arose from the more general setting of evolution operators for transporting flows in dynamical systems @xcite .
the organisation of the article is as follows : we introduce the necessary shell theory and derive a wave scattering model for a singly curved shell in sect .
we then present two approaches for solving the wave scattering model ; a short wavelength asymptotic ray tracing model based on this shell theory is detailed in sect .
[ sec : raytrac ] , and a finite difference discretisation of the full wave model is described in sect .
[ sec : scatscheme ] .
we then discuss and compare numerical results for both the wave and ray scattering models , and the resulting reflection and transmission laws in sect .
[ sec : numer ] .
the thin shell theory of donnell is one of the simplest and widely adopted models @xcite . in this theory ,
moments and transverse forces are expressed by the displacement @xmath0 of the middle surface as known from the theory of laterally loaded plates . as with other theories of continuum mechanics ,
shell theory is formulated in tensor form @xcite .
some properties of tensors are summarised in [ sec : tensors ] .
we assume an isotropic shell of thickness @xmath1 , young s modulus @xmath2 , density @xmath3 and with poisson ratio @xmath4 .
the displacement vector of a point originally on the mid - surface of the shell is decomposed into tangential and normal components thus @xmath5^\text{t}$ ] .
the following tensor equation for the normal displacement @xmath0 may be derived @xcite @xmath6 where @xmath7 are the bending and extensional stiffness , respectively .
all greek alphabet indices take values from the set @xmath8 .
also , the membrane strain is @xmath9 with @xmath10 the second fundamental form and @xmath11 the covariant derivative .
these are discussed further in the next section , where they are simplified for a singly curved shell .
the tangential displacements @xmath12 in the directions @xmath13 , respectively , satisfy @xcite @xmath14 the term @xmath15 represents the inverse of the first fundamental form as explained in [ sec : tensors ] .
note that the above formulation is stated for a shell _ in vacuo _ and , although beyond the scope of the present study , coupling to an acoustic fluid may be included via an additional term in equation ( [ eq : pdesys1 ] ) to express the pressure difference on either side of the shell ; see for example ref .
in addition , structural damping may be incorporated in the usual way by replacing the young s modulus @xmath2 in the equations above with @xmath16 , where @xmath17 is the damping loss factor .
however , we proceed with the lossless case @xmath18 in order to isolate the effect of curvature on the reflection and transmission properties of the shell , whilst also allowing us to check that our results conserve energy in the correct manner .          ) and exterior ( @xmath19 ) regions for the scattering problem .
the interfaces between @xmath20 and @xmath19 lie with the flat regions where @xmath21 . for @xmath22 ,
the curvature increases smoothly to @xmath23 at @xmath24 .
the incident , reflected and transmitted waves in @xmath19 are also indicated . ]    in this section we simplify the shell theory presented above for the case of a cylindrical ridge as shown in fig .
[ fig : ridgegeom ] .
the only simplifying assumption is that the shell is not curved with respect to the direction @xmath25 .
we define the principal curvature @xmath26 ( in the direction @xmath27 ) as @xmath28 where @xmath29 is a constant corresponding to the maximum curvature in the cylindrical region .
a smoothly varying curvature with respect to @xmath27 is obtained via the interpolation function @xmath30 given by @xmath31 where @xmath32 are the centres of the transitions between the two flat regions where @xmath33 and the cylindrical region where @xmath34 , see fig .
[ fig : curv ] . also , @xmath35 is the width of the transition region .
note that the case of a non - smooth joint between a cylindrical ridge and a pair of connecting plates is obtained in the limit @xmath36 .
a cylindrical ridge with a prescribed bending angle @xmath37 can be designed by integrating the angular increment @xmath38 to find the total length @xmath39 of the curved section of the geometry such that @xmath40 for the desired angle .
the length @xmath39 is then found by numerically solving the implicit equation ( [ eq : angleinc ] ) . in this work
we take a quarter cylinder with @xmath41 as shown in fig .
[ fig : ridgegeom ] .    in the above described geometrical setting ,
the shell theory presented in the last section is simplified considerably if we choose an orthonormal basis @xmath42 along the axes @xmath43 .
we then measure all tensor fields relative to this basis , that is , we work in physical components .
the metric becomes @xmath44 and so the first fundamental form and its dual both reduce to the identity matrix @xmath45 therefore the raising and lowering of indices becomes trivial .
the second fundamental form @xmath10 in the basis @xmath42 is simply @xmath46 and zero for all other combinations of @xmath47 and @xmath48 .
the covariant derivative ( see [ sec : tensors ] ) , likewise reduces considerably since all derivatives vanish for the metric ( [ eq : metricid ] ) .
hence , all christoffel symbols vanish and the covariant derivative is identical to the directional derivative , i.e. @xmath49 . applying the above simplifications in the shell equation for the normal displacement ( [ eq : pdesys1 ] ) and assuming that material constants and the thickness are constant yields the following equation @xmath50 where we write @xmath51 and @xmath52 for simplicity of notation .
we adopt this notation for the remainder of the paper since we are no longer considering differential equations written in tensor form . repeating the simplification process for the in - plane wave equation system ( [ eq : pdesys2 ] ) gives @xmath53 \label{eq : pdesys2b}\frac{\rho h}{c }   \frac{\partial^2 u^2}{\partial t^2 } = & \:\frac{\partial^2 u^2}{\partial y^2}+\frac{(1+\nu)}{2}\frac{\partial^2u^1}{\partial x \partial y}+\frac{(1-\nu)}{2}\frac{\partial^2u^2}{\partial x^2}+\nu\kappa\frac{\partial w}{\partial y}.\end{aligned}\ ] ] hence , a simplified set of partial differential equations ( pdes ) describe the wave motion in the cylindrical ridge shell geometry under consideration here .
in the next section we discuss how under certain modelling assumptions , this system can be reduced further to a set of ordinary differential equations ( odes ) .
in this work we are primarily interested in how the scattering properties of a thin shell cylindrical ridge depend on the frequency and the direction of a prescribed incident plane wave originating in the flat region ( where @xmath54 ) . as such it makes sense to consider time - harmonic waves with angular frequency @xmath55 .
further , we can extract scattering properties that are independent of the position along the ridge by exploiting the translational symmetry and assuming that the ridge is of infinite extent in the @xmath56direction . under these assumptions
we may write @xmath57 where @xmath58 is the component of the wavenumber in the @xmath56direction and @xmath59^\text{t}$ ] are the coefficients of ( [ eq : odeansatz ] ) in the in - plane directions @xmath60 and @xmath61 , and the normal direction , respectively .
substitution of the ansatz ( [ eq : odeansatz ] ) into the pde system ( [ eq : pdesys1a ] ) - ( [ eq : pdesys2b ] ) yields the following fourth - order ode system in the variable @xmath60 : @xmath62\label{eq : odesys2 } c_s^2\:\frac{{\mathrm{d}}^2 v}{{\mathrm{d}}x^2}+\left(\omega^2-c_p^2k_y^2\right)v+\mathrm{i } k_y\left((c_p^2-c_s^2)\frac{{\mathrm{d}}u}{{\mathrm{d}}x}+\nu c_p^2\kappa \hat{w}\right)&=0,\\[2pt]\label{eq : odesys3 } \frac{b}{\rho h}\left(\frac{{\mathrm{d}}^4 \hat{w}}{{\mathrm{d}}x^4}-2k_y^2\frac{{\mathrm{d}}^2 \hat{w}}{{\mathrm{d}}x^2}+k_y^4\hat{w}\right)+(c_p^2\kappa^2-\omega^2)\hat{w}+c_p^2\kappa\left(\frac{{\mathrm{d}}u}{{\mathrm{d}}x}+i\nu k_y v\right)&=0.\end{aligned}\ ] ] note that the constants in the above system have been simplified by writing them in terms of the pressure and shear wave velocities , @xmath63 respectively .
in particular , we have made use of the following easily verified relations @xmath64 later in the paper we will write the above ode system in the shorthand form @xmath65 for brevity of exposition .
we connect the curved region to flat plates on each side using interfaces .
imposing a set of conditions at these interfaces will enable us to formulate the set of odes ( [ eq : odesysb ] ) as a boundary value problem , and then by considering the flux of the incoming and outgoing wave fields at these interfaces we are able to formulate a scattering problem .
we assume that the interfaces reside in the asymptotically flat regions and as such the interfaces themselves do not give rise to reflection / transmission phenomena , only the interior region between the interfaces governs the scattering properties .
each interface is assumed to satisfy the continuity conditions given in ref .
@xcite , that is , continuity of displacement , rotation , traction , moment , and normal shear stress . since the material properties are constant throughout the entire geometry , then the interface conditions may be written simply as @xmath66    the superscripts specify the value of the quantity as we approach the interface from either the interior region @xmath67 containing the ridge , or the exterior flat region @xmath68 beyond the interfaces on either side of the ridge . in the sequel
we will refer to the interior region as @xmath69 , and the union of the exterior regions as @xmath19 ( see fig .
[ fig : curv ] ) .
the waves in @xmath19 are precisely those of classical plate theory and therefore the wave modes in @xmath19 that are scattered by @xmath20 will be one ( or more ) of bending , pressure , shear or evanescent bending type . in relation to the vector @xmath70
we have that @xmath71 describes the sum of the bending wave contributions and that the in - plane wave types will each be given by a linear combination of @xmath72 and @xmath73 . in this work
we consider only incident bending modes @xmath74 originating in the exterior @xmath19 and being scattered by the ridge in @xmath69 . however , the extension to other incident wave types is straightforward . for the purpose of investigating directional properties
we consider plane wave scattering .
we introduce the notation that the interface to the left of @xmath69 is located at @xmath75 and the interface to the right is correspondingly at @xmath76 .
then an incident wave of unit amplitude travelling in from the left of @xmath20 can be written in the form @xmath77 where @xmath78 is the wave - number associated with the incident bending mode in the @xmath60-direction . in order to write the resulting scattered waves in a concise manner we introduce @xmath79 and @xmath80 .
then the scattered waves in @xmath19 may be written @xmath81 where @xmath82 designates the scattered wave type as either bending ( @xmath83 ) , evanescent bending ( @xmath84 ) , pressure ( @xmath85 ) or shear ( @xmath86 ) .
the symbol @xmath87 designates whether the scattered wave is reflected ( @xmath88 ) or transmitted ( @xmath89 ) .
that is , whether the scattered wave emerges in @xmath19 on the same side of the ridge as the incident wave is sent in , or on the other side .
the value of @xmath48 also prescribes the sign in the @xmath90 as negative for @xmath91 and positive for @xmath92 .
the coefficient @xmath93 denotes the corresponding wave amplitude , which are related to @xmath70 via the following relations @xmath94 v^{\beta}(x)=&a^p_\beta\exp(\pm \mathrm{i } k^p_x x_\beta)\sin\theta\pm a^s_\beta\exp(\pm \mathrm{i } k^s_x x_\beta ) \cos\theta,\\[2pt]\label{eq : wrt } \hat{w}^{\beta}(x)=&a^b_\beta\exp(\pm \mathrm{i } k^b_x x_\beta)+a^e_\beta\exp(-k^b_x |x_\beta|).\end{aligned}\ ] ] here @xmath95 is the angle between the scattered wave directions and the @xmath96axis ( positive or negatively oriented depending on the direction of propagation ) .
hence , the total wave field in @xmath19 to the left of @xmath20 is given by @xmath97^\text{t},\ ] ] and the total wave field to the right of @xmath20 may be written @xmath98^\text{t}.\ ] ]    the scattering problem is then formed by connecting the interior problem for @xmath99 in @xmath20 with the plane wave ansatz in the exterior regions for @xmath100 and @xmath101 described above , via the interface conditions ( [ eq : if1 ] ) to ( [ eq : if4 ] ) . here @xmath100 and @xmath101 take the role of @xmath102 on the left and right interfaces , respectively .
when this scattering problem is solved , we can extract the scattering solutions @xmath100 and @xmath101 and the scattering coefficients @xmath93 , the latter of which are of primary interest for this study .
one obtains a ray tracing model from the pde model ( [ eq : pdesys1a ] ) - ( [ eq : pdesys2b ] ) by moving to short wavelength asymptotics using the ansatz @xmath103 where @xmath104 is a phase function and @xmath105 is a small parameter . the choice of the parameters @xmath106 and @xmath107 determines the wave type , bending or in - plane .
then we define the frequency @xmath108 and the wavenumber vector @xmath109 .
note that the pre - factors of @xmath110 vanish after applying the asymptotic scaling and setting @xmath111 .
assuming that @xmath112 is large in comparison to the curvature , pierce derived a general dispersion relation @xcite that was later presented in a simpler form by norris and rebinsky @xcite .
in the next section we apply this dispersion relation to generate the hamiltonian dynamics for our ray tracing model on a singly curved shell .      the ring frequency for the cylinder
is defined relative to the pressure mode and corresponds to the frequency above which a longitudinal wave can traverse around the cylinder .
this extensional motion gives rise to breathing modes in cylinders , which also result in a radial motion . in the vicinity of the ring frequency
, approximate forms of the dispersion relations are obtained for the in - plane and bending modes using different scalings . these dispersion relations each reduce to those for a flat plate when the frequency sufficiently exceeds the ring frequency . hence in this study
, we refer to the mid - frequency case as the the frequency range in the vicinity of the ring frequency and the high frequency regime as the frequencies that are sufficiently large for the dispersion relations to reduce to those of a flat plate . defining the longitudinal wavenumber to be @xmath113 , with @xmath114 the longitudinal plate wave speed as before , and denoting the corresponding wavenumber at the ring frequency as @xmath115 , we find that @xmath116 for the singly curved shell configuration shown in fig .  [
fig : ridgegeom ] .     at zeros of the hamiltonian @xmath117 .
( a ) anisotropic dispersion relation for a cylinder . ( b ) isotropic dispersion relation for a flat plate . ]    below the ring frequency , the full dispersion relationship is usually considered . for the configuration in fig.[fig : ridgegeom ] with zero curvature in the @xmath56 direction we obtain the following expression from ref .
@xcite in physical coordinates : @xmath118 two extreme cases of the dispersion curve are depicted in fig.[fig : disp ] , which shows the respective dispersion relations for a cylinder and for a flat plate using the parameters given in [ sec : parameters ] . the outer part of the cylinder dispersion curve has also been observed experimentally @xcite .
the hamiltonian @xmath117 ( [ ham ] ) gives rise to the following ode system , which governs the ray dynamics on a singly curved shell @xmath119 in a full time - domain simulation one would also have two additional equations @xmath120 in this work however , the frequency is constant due to the fact our hamiltonian @xmath117 is time - independent .
furthermore , since we also take the material parameters to be constants throughout the shell , equation ( [ eq : hex1 ] ) simply represents a re - parametrisation of time to the fictitious time used here .
equations ( [ eq : hex1 ] ) and ( [ eq : hex2 ] ) are therefore not needed in the time - harmonic description .
finally , due to the translational invariance in the @xmath61-direction , @xmath58 is constant and it suffices to study the above hamiltonian system in the @xmath121 phase - plane only .
in order to solve the full wave scattering problem we discretise the ode system in @xmath20 using finite difference methods and couple this to the scattering of incoming and outgoing waves in @xmath19 using the interface coupling conditions as described below .
the differential operator @xmath122 in the ode system ( [ eq : odesysb ] ) includes both bending and in - plane waves , along with the coupling between them .
each of the equations in the system is discretised using second order accurate centered finite difference ( fd ) formulae on a set of equi - spaced grid points @xmath123 with @xmath124 .
we denote the centered finite difference matrices for approximating derivatives of order @xmath125 by @xmath126 and define @xmath127\hspace{5 mm } \text { with } \:\ : i , j=0,\dots , n \,\ ] ] to be a diagonal matrix containing the curvature values .    to create a vectorial ode discretisation we use the kronecker product @xmath128 , see for example ref .
@xcite , along with the following @xmath129 projector matrices : @xmath130 \:\text { for } \ : a , b=1,2,3 , \\ { \mathbf{\pi}}_i & = { \mathbf{e}}_{i i } , \\ { \mathbf{h}}_{2 3 } & = \left [ \begin{array}{ccc }   0 & 0 & 0 \\ 0 & 0 & \mathrm{i}\\ 0 & - \mathrm{i } & 0 \end{array } \right ] \,.\end{aligned}\ ] ] using the above notation we can write the various operators of our ode in a compact way .
the in - plane equations ( [ eq : odesys1 ] ) and ( [ eq : odesys2 ] ) will have a diagonal ( self - interaction ) and frequency independent part @xmath131 where @xmath132 is the @xmath129 identity matrix .
the coupling between the @xmath72 and @xmath73 components is symmetric and is given by @xmath133 the bending wave equation ( [ eq : odesys3 ] ) has a diagonal ( self - interaction ) and frequency independent part @xmath134 and the coupling between the bending and in - plane components is given by @xmath135 the shell operator @xmath122 then has the following discrete representation : @xmath136    the fd operators will extend beyond the grid over @xmath20 as the stencil reaches into the scattering region .
this leads to a coupling between the degrees of freedom in the finite difference approximation and the unknown scattering amplitudes @xmath93 .
conversely , the incident wave in @xmath19 will produce a forcing term which drives the finite difference calculation in @xmath20 .
this coupling of the fd solution in @xmath20 and the plane wave description in @xmath19 is discussed further below .
the scattering solution and its derivatives in @xmath19 are matched to the fd solution in @xmath20 at the interfaces between @xmath20 and @xmath19 according to the conditions ( [ eq : if1 ] ) to ( [ eq : if4 ] ) .
the incident wave gives rise to inhomogeneous terms for the excitation of the fd model at the left interface .
the derivatives appearing in the coupling conditions are implemented in the fd - solution in @xmath20 using one - sided finite difference operators with second order accuracy .
forward difference formulae are employed at the left interface and backward difference formulae are employed at the right interface so that the stencils extend only into @xmath20 .
the extremal points of these stencils are taken at the scattering boundaries , where the value of the scattering solution is applied as a boundary condition .
the interfaces are thus taken to be at positions @xmath137 away from the extremal inner grid points , where @xmath137 is the discretisation step in the fd approximation .    as an illustrative example we consider the condition on @xmath138 ( [ eq : if4 ] ) at the interface receiving the incident wave ( [ eq : winc ] ) .
the reflected bending modes include both a propagative contribution @xmath139 and an evanescent contribution @xmath140 as shown in equation ( [ eq : wrt ] ) .
the forward difference formula for the third derivative takes the form @xmath141 matching the scattering and finite difference solutions using the interface condition ( [ eq : if4 ] ) leads to @xmath142 splitting the discretised boundary operator @xmath143 into a part that acts only on the left - most entry at the interface and one that acts on the remaining nodes in @xmath20 , we obtain @xmath144 at @xmath145 , the left hand side of ( [ eq : ifdisc ] ) evaluates to @xmath146 and using the splitting described above , the right hand side may be written as @xmath147 writing @xmath148 for @xmath149 and using the definition of @xmath143 yields a linear equation @xmath150 in terms of the unknowns @xmath151 , @xmath152 and @xmath153 , @xmath149 .
continuing for the remaining boundary conditions at the left interface @xmath154 yields similar relations , all following the same pattern as in ( [ eq : vcond ] ) .
the boundary conditions at the right interface give similar results for the transmitted terms , but without the incident wave terms . finally , we combine the discretised interface coupling conditions with the interior finite difference equations .
this leads to the matrix problem : @xmath155 \left[\begin{matrix } \mathbf{r}\\ \hat{\mathbf{u}}_\delta\\ \mathbf{t } \end{matrix } \right]= \left[\begin{matrix } * \\
\mathbf{0}\\ \mathbf{0 } \end{matrix } \right]\ ] ] with scattering coefficients @xmath156^\text{t},\ ] ] and @xmath157^\text{t}.\ ] ] here , the finite difference solution in the interior of @xmath20 is represented by @xmath158^\text{t},\:\ : i=1,\dots , n-1.\ ] ] note that the stencil for the second derivative fd operator within @xmath159 takes boundary data from the plane wave terms at the interface points @xmath160 and @xmath161 , whereas the fourth derivative operator extends beyond this and also takes boundary data from within the scattering region at @xmath162 and @xmath163 .
in this section we discuss the numerical solution of the full wave scattering problem derived in sect .
[ sec2 ] and compare the results to those obtained using an ode time - stepper for the hamiltonian system presented in sect .
[ sec : raytrac ] .
the parameters chosen for all calculations are summarised in [ sec : parameters ] . for the cylindrical region
, these parameters correspond to a circularly cylindrical steel shell as considered in refs .
@xcite . note that whilst the width of the cylindrical ridge studied here is on the millimeter scale , it only represents a small region of the larger structures and components that serve to motivate this study .
an incident bending mode in the left part of @xmath19 is used to excite the system and is sent in to the interface with @xmath20 at various angles , corresponding to a variation in the trace wave number @xmath164 .
the latter is chosen according to the values given by the dispersion curve . in all cases we set @xmath165 rad / s , which is large enough so that waves will always transmit straight through the cylindrical region ( as though it were a flat plate ) when
@xmath166 , that is when waves approach the cylindrical region directed parallel to the @xmath60-axis .
we note that in the high frequency regime this behaviour will be preserved for almost all incident waves when both @xmath167 and @xmath168 are positive .
however , the choice of @xmath55 here corresponds to the more interesting mid - frequency case where reflections are also possible for a range of sufficiently large @xmath169 .
phase - plane showing two fixed points at @xmath24 and the separatrix ( solid bold curve ) which connects them as it traverses the cylindrical ridge .
( b ) a close up view of the upper fixed point and separatrix region .
the upper fixed point divides the incoming rays from the left hand side ( dashed lines ) into reflecting or transmitting trajectories .
the bold dashed line depicts the extremal transmitting trajectory closest to the set of reflecting rays . in both subplots ,
the solid line trajectories are those which are trapped inside the ridge . ]    in this section we consider the reflection/ transmission behaviour of rays corresponding to an incident bending mode ( see fig .
[ fig : phasespace4 ] ) .
a range of incoming wavenumbers ( directions ) are used corresponding to the strip in the upper left corner of the figure .
we find that the rays transmit for sufficiently large @xmath78 , whereas for smaller positive values of @xmath78 the rays reflect .
note that the symmetry of the problem means that rays also transmit for @xmath78 negative when the incoming ray is from the right hand side and @xmath170 is large enough .
the threshold value of @xmath171 for the change from reflectance to transmittance corresponds to a hyperbolic fixed point and its location can be found conveniently by considering the sizes of the dispersion curves as discussed in @xcite . the upper fixed point shown in both parts of fig .
[ fig : phasespace4 ] gives a threshold value of @xmath172 .
we will also investigate the existence and location of such a threshold incoming wave direction in the finite difference solution of the full wave problem described in the next section .
depending on the angle of incidence , fig .
[ fig : phasespace4 ] shows that the point of reflection varies from the centre of the ridge at @xmath24 to slightly towards the flat region to the left .
the smooth curvature model here therefore differs from the discontinuous curvature models considered in refs .
@xcite , for which the reflection takes place off - centre at a fixed location only .
that is , for the models presented in refs .
@xcite the reflection would take place at the location of the jump in the curvature between the flat region and the cylindrical region .
this would be off - centre in the example here , since the centre point @xmath24 corresponds to the centre of the cylindrical region .    in the next section we consider an equivalent scattering problem to the one above , but instead using the finite difference model described in the previous section to numerically solve the full wave problem .
the full wave model will include the phase information omitted in the pure ray approximation applied in this section , but at the cost of a computational expense which scales with frequency .
such a study is feasible up to reasonably high frequencies due to the one - dimensional setup of the problem derived in sect .
however , the high frequency purely transmissive behaviour for the problem here is relatively straightforward to predict , and we find that both methods may be used for the more interesting mid - frequency case studied here .
we restrict the study to bending excitations as in the previous section . for smooth joints
this means that only the bending mode is active , with negligible conversion to in - plane modes .
that is , the scattering probabilities become @xmath173 ) with very small @xmath35 in the interpolation function ( [ eq : interp ] ) , mode conversions appeared in the numerics at almost normal incidence .
we defer the study of this case to future work and note the possibility of using discontinuous joints instead for this case .
note that we use the term probabilities to describe the coefficients defined in ( [ eq : prob ] ) since for an incoming bending wave of unit amplitude ( [ eq : winc ] ) , conservation of energy gives that @xmath174 . in the full wave picture considered here these coefficients ( [ eq : prob ] )
actually give the proportions of reflected and transmitted wave energy .
however , for comparison with the ray tracing results in the last section , they give the probabilities of reflection and transmission for a particular incoming trajectory .     plotted on a linear scale .
the plot shows anti - resonances close to @xmath175 , @xmath176 and @xmath177 . there is a switch in the dominant behaviour from transmission to reflection close to @xmath178.,scaledwidth=90.0% ]     plotted on a logarithmic scale .
the logarithmic scale emphasises the anti - resonance reflection peaks and the switch from transmission to reflection that were evident from fig .
[ fig : coeff].,scaledwidth=90.0% ]    the dependence of the reflection and transmission probabilities on the value of the trace wave number @xmath168 is shown in fig .
[ fig : coeff ] .
the transmission and reflection probabilities have been computed using the finite difference discretisation described in sect .
[ sec : scatscheme ] for both @xmath179 and @xmath180 grid points in @xmath20 , corresponding to @xmath181 m and @xmath182 m , respectively .
the ratio of wavelength to @xmath137 is @xmath183 for the case @xmath179 and hence the rule of thumb requiring 6 points per wavelength for reliable results suggests that in both cases our results should be reasonably well converged .
this convergence is also evident from fig .
[ fig : coeff ] , since it is difficult to detect differences between the plots for @xmath179 and @xmath180 .    fig .
[ fig : coeff ] shows that the scattering coefficients of the full wave problem include discrete anti - resonance points ( points of perfect reflection ) that were not present in the ray tracing calculation .
in addition , the threshold behaviour from transmission to reflection has been smoothed in the full wave calculations , whereas for the ray tracing solution there is a sudden jump from transmission to reflection corresponding to the location of the hyperbolic fixed point .
hence , for the calculations in this section there is a region of @xmath168 values where both reflection and transmission take place at the same time .
this region lies roughly between @xmath184 and @xmath185 , but one observes an obvious switch from dominant transmission behaviour to dominant reflection behaviour at a threshold value around @xmath186 , which is close to the threshold estimate ( @xmath187 ) from the corresponding ray tracing calculation .
it appears therefore that the asymptotic ray model slightly over predicts the transmission / reflection threshold compared with the full wave solution , but still produces a reliable estimate within around @xmath188 of the total @xmath168 range .
a more complete picture appears when the scattering coefficients are plotted on a logarithmic scale , see fig .
[ fig : logcoeff ] .
we can then clearly identify the discrete reflection points in the region of perfect transmission . the dips or peaks
can be arbitrarily close to zero , or unity , depending on the resolution of the plot .
the discrepancy between the finite difference solution for the reflection coefficient with @xmath179 and @xmath180 is also more evident in this figure since the logarithmic scaling amplifies the differences at very small values of @xmath189 below @xmath190 .
however , the positions of the anti - resonances and the switch from transmission to reflection remain in excellent agreement for both @xmath179 and @xmath180 .
note also that the values of @xmath191 computed with @xmath179 and @xmath180 are in good agreement with one another , even for very small values .
$ ] for @xmath192 ( plots ( a ) to ( c ) ) , @xmath177 ( plots ( d ) to ( f ) ) , @xmath193 ( plots ( g ) to ( i ) ) and @xmath194 ( plots ( j ) to ( l ) ) .
each plot shows a displacement in metres against its corresponding @xmath60-coordinate , also in metres .
plots ( c ) , ( f ) , ( i ) and ( l ) show the bending mode @xmath71 which exhibits transmission in plot ( c ) , anti - resonance in plot ( f ) , mixed reflection and transmission for the threshold region in plot ( i ) and reflection in plot ( l ) .
the remaining plots show the in - plane wave solution , which remains confined inside the cylindrical region in all cases .
solid and dashed lines show the real and imaginary parts , respectively . ]
figure [ fig : wfmiddle ] shows the behaviour of the wave function solutions @xmath59^\text{t}$ ] for various values of the trace wavenumber @xmath168 .
the finite difference computations for these plots were all carried out using @xmath180 grid points .
note that there is no mode conversion and so the in - plane contributions @xmath72 and @xmath73 only have support on the curved region and vanish in @xmath19 . in the case of reflection ( @xmath194 ) and transmission ( @xmath192 )
the expected behaviour is observed .
that is , for the reflective case the bending wave function @xmath195 localises to the left as shown in plot ( l ) of fig .
[ fig : wfmiddle ] . in the ray
tracing model one can pinpoint an exact reflection point , but for the wave problem one observes a decay in amplitude as the wave enters the cylindrical region rather than a hard wall reflection .    in the case of transmission shown in plot ( c ) of fig .
[ fig : wfmiddle ] , the bending wave continues through the ridge without significant deformation from the shape of the initial plane wave excitation .
this is equivalent to the corresponding ray tracing result .
however , if @xmath168 is chosen coincide with one of the anti - resonance peaks in the transmission region that were not present in the ray tracing result , then the bending wave field is localised in the central cylindrical region ; for example , at @xmath177 as shown in plot ( f ) of fig .
[ fig : wfmiddle ] . the bending wave function @xmath195 also localises in the cylindrical region when @xmath168 is chosen in the threshold region as shown in plot ( i ) of fig .
[ fig : wfmiddle ] , which shows a plot of @xmath195 at @xmath193 .
we have investigated bending wave scattering across a smoothed plate and quarter cylinder configuration in the interesting mid - frequency case , close to the ring frequency for the cylinder .
results have been obtained in a simplified and effectively one - dimensional setting using both a high - frequency ray - tracing approximation and a finite difference discretisation of the equivalent full - wave problem .
previous studies on plate and cylinder connections have concentrated on joints having discontinuous curvature where reflections can occur , however the smoothed joints considered here also give rise to reflections within the cylindrical part of the structure .
furthermore , the ray tracing calculations suggest the existence of a threshold incident wave direction which separates waves or rays that exhibit reflective or transmissive behaviour .
this threshold direction is also evident from the smoothed transition between reflection and transmission observed in the full wave calculations .
hence , relatively simple scattering laws can be employed to model the propagation of structure - borne noise in shells , and ultimately in built - up structures containing thin shell components .
the full wave solution shows that in addition to the switch from transmission to reflection as the incident wave direction becomes increasingly oblique , there are also anti - resonances giving rise to perfect reflection at a discrete set of directions where transmission would typically be expected . in these cases
the wave functions appear as trapped modes that localise in the cylindrical part of the configuration .
there are several avenues for further research , including extensions to multiply - curved and fluid - loaded shells .
in addition , the scattering laws could be incorporated within computed aided engineering simulations of built - up structures via wave methods , such as the wave and finite element method or dynamical energy analysis .
support from the eu ( fp7-people-2013-iapp grant no .
612237 ( mhivec ) ) is gratefully acknowledged .
we also wish to thank dr gregor tanner for stimulating discussions and dr jonathan crofts for carefully reading the manuscript
.    99 o.c .
zienkiewicz , r.l .
taylor , j.z .
zhu , _ the finite element method : its basis and fundamentals _ , 7th edition , butterworth - heinemann , oxford , 2003 .
lyon , r.g .
dejong , _ theory and application of statistical energy analysis _
, 2nd edition , butterworth - heinemann , boston ma , 1995 .
glasser , _ an introduction to ray tracing _ , academic press , new york , 1989 .
a. le bot , _ foundation of statistical energy analysis in vibroacoustics _ , oxford university press , oxford , 2015 .
h. kuttruff , _ room acoustics _ , 4th edition , spon , london , 2000 .
v. cerveny , _ seismic ray theory _ , cambridge university press , cambridge uk , 2001 .
norris , d.a .
rebinsky , membrane and flexural waves on thin shells .
_ journal of vibration and acoustics _ 116 ( 1994 ) 457467 .
doi:10.1115/1.2930449 e.g. williams , b.h.houston , j.a .
bucaro , experimental investigation of the wave propagation on a point - driven , submerged capped cylinder using k - space analysis .
_ journal of the acoustical society of america _
87(2 ) ( 1990 ) 513522 .
doi:10.1121/1.398922 r.s .
langley , wave motion and energy flow in cylindrical shells .
_ journal of sound and vibration _
169(1 ) ( 1994 ) 2942 .
doi:10.1006/jsvi.1994.1004 a.n .
norris , reflection and transmission of structural waves at an interface between doubly curved shells .
_ acta acustica united with acustica _
84(6 ) ( 1998 ) 10661076
. g. tanner , dynamical energy analysis - determining wave energy distributions in vibro - acoustical structures in the high - frequency regime .
_ journal of sound and vibration _ 320 ( 2009 ) 10231038 .
doi:10.1016/j.jsv.2008.08.032 d.j .
chappell , d. lchel n. sndergaard , g. tanner , dynamical energy analysis on mesh grids : a new tool for describing the vibro - acoustic response of engineering structures . _ wave motion _
51(4 ) ( 2014 ) 589597 .
doi:10.1016/j.wavemoti.2014.01.004 e. manconi , b.r .
mace , wave characterization of cylindrical and curved panels using a finite element method . _ journal of the acoustical society of america _
125(1 ) ( 2009 ) 154163 .
doi:10.1121/1.3021418 j.m .
renno , b.r .
mace , calculation of reflection and transmission coefficients of joints using a hybrid finite element / wave and finite element approach .
_ journal of sound and vibration _
332(9 ) ( 2013 ) 21492164 .
doi:10.1016/j.jsv.2012.04.029 j.t .
kayija , the rendering equation .
_ acm siggraph computer graphics _
20(4 ) ( 1986 ) 143150 .
doi:10.1145/15886.15902 p. gaspard .
_ chaos , scattering and statistical mechanics _ , cambridge university press , cambridge uk , 1998 .
p. cvitanovi , r. artuso , r. mainieri , g. tanner , g. vattay _ chaos : classical and quantum _ , chaosbook.org , niels bohr institute , copenhagen , denmark , 2012 .
pierce , waves on fluid - loaded inhomogeneous elastic shells of arbitrary shape .
_ journal of vibration and acoustics _
115(4 ) ( 1993 ) 384390 .
doi:10.1115/1.2930361 w.  flgge , _ tensor analysis and continuum mechanics _ , springer , berlin , 1972 .
van loan , the ubiquitous kronecker product .
_ journal of computational and applied mathematics _ 123 ( 2000 ) 85100 .
doi:10.1016/s0377 - 0427(00)00393 - 9 d.a .
rebinsky , a.n .
norris , dispersion of flexural waves on shells .
_ journal of vibration and acoustics _ 118(3 ) ( 1996 ) 526529 .
doi:10.1115/1.2888218 n. sndergaard , d.j .
chappell , g. tanner , tracking vibrational energy on curved shell structures in the mid - high frequency limit - a ray - tracing approach . _ proceedings of noise and vibration - emerging technologies ( novem ) 2015 _ , dubrovnik , croatia , 2015 .
m. do carmo , _ differential geometry of curves and surfaces _ ,
prentice hall , new jersey , 1976 .
t.  frankel , _ the geometry of physics _ , cambridge university press , new york , 2004 .
several physical theories are concisely written in tensor form @xcite .
the tensor formalism has also proved to be useful in continuum mechanics .
the components of tensors are with respect to a given choice of coordinates @xmath196 .
a change of coordinates from @xmath196 to @xmath197 leads to expressions for components in the new coordinates in relation to those of the old .
for example , given a tensor of type @xmath198 ( that is , with one superscript and two subscript indices ) the transformation takes the form @xmath199 where jacobians and inverse jacobians have been used for the change of coordinates .
we have adopted the important summation convention for repeated indices ; this indicates that a summation is to take place over the repeated index ( unless otherwise stated ) .    the metric tensor @xmath200 , known as the first fundamental form , measures distances via @xmath201 the inverse of the metric tensor @xmath15 obeys @xmath202 and is used to alter the type of a tensor .
for example , we can introduce a new tensor of type ( 1,1 ) @xmath203 from a tensor of type ( 0,2 ) by raising the indices of the latter tensor .
likewise @xmath200 lowers indices .
the directional derivative in arbitrary coordinates is generalized to the covariant derivative @xmath204 .
this derivative may be introduced from an embedding using projection of gradients @xcite or intrinsically with quantities only related to the curved space itself @xcite .
for example , in a coordinate basis , the covariant derivative of a tensor of type ( 1,1 ) becomes @xmath205 with @xmath206 the christoffel symbols .
the computations throughout this work were done using the following parameter choices @xcite :


it is well known that any link can be represented by the closure of a braid .
the minimum number of strands needed in a braid whose closure represents a given link is called the _ braid index _ of the link .
defined as the extreme value of a quantity over an infinite family of links that are topologically equivalent like other link invariants ( such as the minimum crossing number ) , the braid index of a link is hard to compute  @xcite in general . in the case of the minimum crossing number , there was a long standing conjecture which states that for a reduced alternating link diagram , the number of crossings in the diagram is equal to the minimum crossing number of the link .
it is known that the span of the jones polynomial @xcite of a link gives a lower bound on the crossing number of the link . in the case of a reduced alternating link diagram
, it was shown that the span of the jones polynomial of a link equals the number of crossings in the diagram , which leads to the proof of the conjecture @xcite . in the case of braid index
, there is a similar inequality relating the braid index of a link to the @xmath1-span of its homfly polynomial ( which is a polynomial of two variables @xmath2 and @xmath1 to be defined in the next section ) .
s.  yamada proved that any link diagram of a given link @xmath3 with @xmath4 _ seifert circles _ can be realized as the closure of a braid on @xmath4 strands , which implies that the braid index of an oriented link @xmath3 equals the minimum number of seifert circles of all link diagrams of @xmath3  @xcite . in  @xcite
, h.  morton showed that the number of seifert circles of a link @xmath3 , hence the braid index of @xmath3 ( in light of yamada s result ) , is bounded from below by @xmath1-span@xmath5 ( which is called the morton - frank - williams inequality , or mfw inequality for short ) .
in analogy to the crossing number conjecture for a reduced alternating link diagram , k.  murasugi conjectured that the number of seifert circles , hence the braid index , in such a diagram equals @xmath1-span@xmath5 ( the murasugi conjecture ) @xcite .
although this conjecture turned out to be false in general  @xcite ( for example the knot @xmath6 has 4 seifert circles , but the @xmath1-span of its homfly polynomial is 4 so @xmath1-span@xmath7 and the braid index of @xmath6 is also 3 ) , researchers had shown that the mfw inequality is sharp for many classes of links ( including some non - alternating ones ) hence the @xmath1-span of the homfly polynomial for these links can be used to determine their braid index .
examples include the closed positive braids with a full twist ( in particular the torus links )  @xcite , the 2-bridge links and fibered alternating links  @xcite , and a new class of links discussed in a more recent paper by s.  y.  lee and m.  seo  @xcite . for more readings on related topics , interested readers can refer to j.s .
birman and w.w .
menasco  @xcite , p.r .
cromwell  @xcite , e.a .
elrifai  @xcite , h.  morton , h.  b.  short  @xcite , t.  nakamura  @xcite and a.  stoimenow  @xcite .    motivated by this question , in this paper the authors seek a special and explicit formulation of the homfly polynomial for certain classes of links where the explicit forms of the homfly polynomial would allow us to analyze and derive the @xmath1-spans of the homfly polynomials of these links .
our main result expresses the homfly polynomial of a link presented in a closed braid diagram form in two explicit formulas .
we show that one of our formulas of the homfly polynomial is equivalent to the expansion derived by f.  jaeger  @xcite . however , our approach is combinatorial in nature and the proof of the formula is shorter than the proof in  @xcite .
we show that the morton - frank - williams inequality is an immediate consequence of these two homfly polynomial formulas . as a more significant application of our result , we use it to show that if a link has an @xmath0-string closed braid diagram that is also reduced and alternating , then the braid index of the link is exactly @xmath0 .
the proof of this is direct and short .
our approach is very different from the proof of murasugi s more general results  @xcite on oriented alternating fibered links , which is inductive in nature .
this paper is structured as follows . in section  [ s2 ]
, we introduce basic definitions and terminology about link diagrams , braids and the homfly polynomial .
we introduce two special classes of resolving trees for closed braids , in which every vertex is a closed braid and the leaf vertices are closed braids that represent trivial links .
we call these resolving trees descending trees and ascending trees respectively . in section  [ s3 ] , we state and
prove our main result , namely two formulas expressing the homfly polynomial of a closed braid as a total contribution of all leaf vertices in our trees and show that the morton - frank - williams inequality  @xcite is a direct consequence of these two formulas . in section  [ s4 ] , we show that the @xmath1-span of the homfly polynomial of a reduced alternating braid on @xmath0 strands is exactly @xmath8 and an application of the the morton - frank - williams inequality shows that the braid index of a reduced alternating braid equals the number of strands in the braid . as another application of our main result from section  [ s3 ] , we also give a short proof that the leading coefficient of the alexander polynomial of such a closed braid equals @xmath9 .
finally , in section  [ acp ] , we show that one of our formulas is equivalent to the expansion of the homfly polynomial derived by f.  jaeger  @xcite based on the concept of admissible circuit partitions .
we assume that the reader has the basic knowledge about the definition of a link and refer a reader without such knowledge to a textbook such as @xcite .
for the convenience of the reader , however , we will review one important result that is needed in section [ s3 ] .
figure [ r - moves ] defines three moves one can make on a link diagram without changing its topology , and these are called reidemeister moves of type i , ii and iii . in 1926 , k.  reidemeister  @xcite proved that two link diagrams represent the same link if and only if one diagram can be changed to the other through a finite sequence of reidemeister moves .        for a given oriented link diagram @xmath3 , we assign each crossing @xmath10 or @xmath11 according to its sign as defined in figure  [ fig : cross ] .
the writhe of @xmath3 , written as @xmath12 , is the sum of these @xmath10 s and @xmath11 s over all crossings of @xmath3 .
however , if we only sum the @xmath10 s and @xmath11 s over crossings between two different components @xmath13 and @xmath14 of the link , then this number is always even and half of it is called the _ linking number _ between @xmath13 and @xmath14 . using the reidemeister moves , it is easy to show that the linking number is a link invariant . in particular ,
if @xmath13 and @xmath14 are separated by a topological plane ( in which case @xmath13 and @xmath14 are said to be _
splittable _ ) , then the linking number between @xmath13 and @xmath14 is zero .
we make the following remark for future reference .
its proof is trivial and is left to the reader .
[ r_remark ] _ reidemeister moves of type ii and iii do not change the writhe of a link diagram hence link diagrams related by a finite sequence of reidemeister moves ii and iii have the same writhe . in particular , if a link diagram @xmath15 is obtained through @xmath16 by a finite sequence of moves that involve only deformation of segments of the link within planes that are parallel to the projection plane of the link diagram , then @xmath17 since such moves do not introduce reidemeister moves of type i. also , if @xmath3 is a link with splittable components @xmath13 , @xmath14 ,  , @xmath18 , then @xmath19 .
_    a crossing in a link diagram is called _ nugatory _ if there is a simple closed curve such that the link diagram intersects the curve only at the crossing and the link diagram intersects both components of the complement of the curve .
a link diagram is said to be _ reduced _ if no crossing in the diagram is nugatory .
a link diagram is _ alternating _ if as we travel through the link diagram by any given orientation , the strands go through the crossings alternately between overpasses and underpasses .
for example , the closure of the braid in figure [ fig : bd ] is an alternating link diagram .
as we mentioned in the introduction , reduced alternating link diagrams are special since the number of crossings in a reduced alternating link diagram is the minimum crossing number of the link @xcite .
consider @xmath20 as the standard euclidean @xmath21-plane . a _ braid diagram _
( or just a _ braid _ ) on @xmath0 strands is a set @xmath22 $ ] consisting of @xmath0 curves called @xmath23 of @xmath24 such that the following four conditions are met .
first , each strand is a monotonic curve in the @xmath25 direction .
second , every point of @xmath26 is a starting point of a unique strand and every point of @xmath27 is an end point of a unique strand .
third , every point of @xmath28 belongs to at most two strands . a point that belongs to two strands
is called a _
crossing_. at each crossing , one strand is identified as an _ overpass _ and the other is as an _ underpass _  @xcite .
fourth , there is at most one crossing in @xmath29 for each @xmath30 $ ] .
note that the second condition gives the braid diagram a downward orientation , so the closure of a braid diagram is an oriented link diagram .
treated as topological objects , one can speak of topological equivalence of braid diagrams , and such equivalence relations provide the foundations for one to treat the braids as elements in the algebraic objects called the _
braid groups_. not to deviate from our main task , we will only point out that a braid group @xmath31 is a group with @xmath32 generators @xmath33 satisfying certain relations .
an element of @xmath31 is a word of these generators and each letter in the word corresponds to a crossing in the braid .
an example of a braid on 5 strands and its counterpart in @xmath34 is given in figure [ fig : braid ] .
please refer to  @xcite for more details on braid groups . in this paper
we are only interested in closed braids as topological objects , not the braids in the algebraic sense .
for this reason , _ we will not distinguish a braid and its closure _ , that is , the word braid ( or a symbol of it ) can either represent the braid itself or its closure .
the reader should rely on the context to determine its meaning ( in many cases it really does not matter ) .
we define the _ label _ of a strand by the @xmath35 coordinate of its starting point and we denote the corresponding mapping by @xmath36 , that is , if a strand @xmath37 starts at @xmath38 , then @xmath39 .
on the other hand , the mapping @xmath40 that takes the label of a strand to the @xmath35 coordinate of its ending point defines a permutation of the labels ( namely the integers from @xmath41 to @xmath0 ) .
denote this permutation by @xmath42 and write it as a product of disjoint cycles , we have @xmath43 where @xmath44 is the label of the first strand in the @xmath45-th cycle , @xmath46 is the label of the second strand in the @xmath45-th cycle , ... , and so on , @xmath47 is the number of cycles in the permutation .
furthermore , we can re - arrange the orders of the cycles and the numbers in each cycle so that @xmath44 is the smallest integer in each cycle for each @xmath45 , and @xmath48 if @xmath49 .
we call this special form of @xmath42 the _ standard form_. from now on , @xmath42 will always be expressed in its standard form . note that the standard form of @xmath42 defines a total order among the strand labels , namely @xmath50 we call this order the _ return order _ of the strands in the braid diagram @xmath24 .
we call each @xmath44 in @xmath42 the _ pivot label _ within its corresponding cycle and @xmath51 the _ pivot point _ of the cycle when @xmath42 is expressed in its standard form .
note that each cycle in @xmath42 corresponds to a connected component in the closed braid diagram and we can travel through @xmath24 by traveling through each such component .
we say that we travel through @xmath24 _ naturally _ if we travel along the strands of @xmath24 in its return order , starting from the pivot point at each component and follow the orientation of the braid @xmath24 .     on 5 strands that has 8 crossings .
]    _ consider the braid diagram @xmath24 shown in figure  [ fig : braid ] whose corresponding word in @xmath34 is @xmath52 .
the dotted curves show how the braid is closed to form a link diagram .
the circled points on top are pivot points of the three components .
we have @xmath53 ( in its standard form ) , so its return order is @xmath54 .
one can see that if we close the braid using disjoint curves starting and ending at the points @xmath55 , @xmath56 for each @xmath57 ( shown in figure [ fig : braid ] with dotted curves ) , strands @xmath41 and @xmath58 are in the same component , strand @xmath59 is a component by itself and strands @xmath60 and @xmath61 are in the same component .
the pivot labels of these components are @xmath41 , @xmath59 and @xmath60 respectively .
_    let @xmath62 be a crossing in a braid diagram @xmath24 , @xmath63 be the overpassing strand at @xmath62 and @xmath64 be the underpassing strand at @xmath62 .
we say that @xmath62 is _ descending _ if @xmath65 and _ ascending _ if @xmath66 .
if all crossings in @xmath24 are descending ( ascending ) , then we say that @xmath24 is a descending ( ascending ) braid diagram .
for example , all crossings except the circled one in figure [ fig : braid ] are descending crossings . switching the circled crossing to a descending crossing
will then results in a descending braid .
[ descending_remark ] _ descending ( ascending ) braid diagrams have the easy and well known property that the closure of any such braid diagram is a trivial link , that is , a link topologically equivalent ( ambient isotopic ) to the disjoint union of several circles contained in the same plane .
take the braid given in figure  [ fig : braid ] as an example by flipping the circled crossing to make the braid descending . using the order @xmath54 , we can actually move strands 1 , 4 , 2 , 3 and 5 into planes @xmath67 , 4 , 3 , 2 and 1 respectively so that their projections would still give the descending braid , and that this move will not change the topology of the corresponding link ( since there are no crossing changes in the process ) .
it is then easy to see that each component corresponding to a cycle @xmath68 , @xmath69 or @xmath70 is an unknot and these knots are splittable since they are separated by planes . _
let @xmath71 , @xmath72 , and @xmath73 be oriented link diagrams that coincide except at a small neighborhood of a crossing where the diagrams are presented as in figure  [ fig : cross ] : the crossing in @xmath71 ( @xmath72 ) is positive ( negative ) and is assigned @xmath10 ( @xmath11 ) in the calculation of the writhe of the link diagram .
we say the crossing presented in @xmath71 is of a _ positive _ sign and the crossing presented in @xmath72 is of a _ negative _ sign .
the following result appears in  @xcite .
[ ho ] there is a unique function that maps each oriented link diagram @xmath3 to a two - variable laurent polynomial with integer coefficients @xmath74 such that    1 .   if @xmath16 and @xmath15 are ambient isotopic , then @xmath75 .
2 .   @xmath76 .
3 .   if l is an unknot , then @xmath77 .
the laurent polynomial @xmath74 is called the _ homfly polynomial _ of the oriented link @xmath3 .
the second condition in the proposition is called the _ skein relation _ of the homfly polynomial . with conditions ( 2 ) and ( 3 ) above , one can easily show that if @xmath3 is a trivial link with @xmath0 connected components , then @xmath78 ( by applying these two conditions repeatedly to a simple closed curve with @xmath32 twists in its projection ) . for our purposes
, we will actually be using the following two equivalent forms of the skein relation : @xmath79    a rooted and edge - weighted binary tree @xmath80 is called a _ resolving tree _ of an oriented link diagram @xmath3 ( for the homfly polynomial ) if the following conditions hold .
first , every vertex of @xmath80 corresponds to an oriented link diagram .
second , the root vertex of @xmath80 corresponds to the original link diagram @xmath3 .
third , each leaf vertex of @xmath80 corresponds to a trivial link .
fourth , if we direct @xmath80 using the directions of the paths from the root vertex to the leaf vertices , then under this direction any internal vertex has exactly two children vertices and the corresponding link diagrams of these three vertices are identical except at one crossing and they are related by one of the two possible relations at that crossing as shown in figure  [ fig : rtree ] , where the edges are weighted and the directions of the edges coincide with the direction of @xmath80 .
[ tree_formula_remark ] _ if @xmath3 admits a resolving tree @xmath80 , then one can easily show that @xmath74 is a summation in which each leaf vertex of @xmath80 contributes exactly one term in the following way .
let @xmath81 be the trivial link corresponding to a leaf vertex in @xmath80 and let @xmath82 be the unique path from the root ( @xmath3 ) to the leaf vertex ( @xmath81 ) . then the contribution of the leaf vertex is simply @xmath83 multiplied by the weights of the edges in @xmath82 , where @xmath84 is the number of components in @xmath81 .
it is known ( and not hard to prove ) that resolving trees exist for any given oriented link diagram @xmath3 , and that they are not unique in general .
if @xmath85 is the mirror image of @xmath3 , a resolving tree for @xmath85 can be obtained from a resolving tree of @xmath3 by taking mirror images of all link diagrams in it and replacing @xmath1 by @xmath86 and @xmath2 by @xmath87 in the edge weights .
it follows the well known fact that @xmath88 . _    ) and
( [ skein2 ] ) . ]
let @xmath89 be a braid on @xmath0 strands and let us define an algorithm ( let us call it algorithm d ) that operates on @xmath89 as follows .
if @xmath89 is descending , then the algorithm does not do anything and simply returns @xmath89 .
if @xmath89 is not descending , then it contains at least one ascending crossing . in this case , let us travel through @xmath89 naturally , until we run into the first ascending crossing @xmath62 , which can be a positive or a negative crossing as shown in either the left hand side or the right hand side of figure  [ fig : rtree ] .
the algorithm then returns the two corresponding braids split from the original one as shown in figure  [ fig : cpd ] , which we will name as @xmath90 and @xmath91 ( @xmath92 for  flipping " the crossing and @xmath37 for  smoothing " the crossing ) .
notice that @xmath90 has one less ascending crossing than @xmath89 does since @xmath90 and @xmath89 are identical ( including the return order of the strands ) except at the crossing @xmath62 .
but one can not say the same for @xmath91 since it does not share the same strands ( and the return order of the strands ) with @xmath89 .
however , @xmath89 and @xmath91 share the same strands and crossings up to crossing @xmath62 ( which are all descending ) while travel through them naturally , and @xmath91 has one less crossing than @xmath89 .
so if algorithm d is repeatedly applied to a braid and the resulting braids of this operation , then this process will end after a finite number of repeats ( in fact this number is bounded above by the number of crossings in @xmath24 ) .
it follows that we can construct a special resolving tree @xmath80 for @xmath93 ( as a link diagram ) as follows .
we apply algorithm d to @xmath24 first if it is not descending , and then apply the algorithm again to the two resulting braids ( if they are not descending ) , and so on , until all the leaf vertices are descending braids .
the closures of the braids involved in the process are the vertices of @xmath80 .
in particular , the closures of the resulting descending braids form the leaf vertices of @xmath80 since they are all trivial links . by assigning appropriate weights to the edges of @xmath80
, one can easily verify that @xmath80 is indeed a resolving tree of @xmath93 . by the way it is constructed
, @xmath80 is unique . in a similar fashion
, we can also construct another ( unique ) resolving tree of @xmath93 by replacing  descending " with  ascending " in the above ( the algorithm corresponding to algorithm d would be called algorithm a ) . to distinguish the two
, we will call the first the  descending tree " and the later the  ascending tree " of @xmath93 , and denote them by @xmath94 and @xmath95 respectively .
an example of a descending tree is shown in the right hand side of figure  [ fig : cpd ] .    ;
right : the descending tree of ( the closure of ) the braid on the top left ( only the braids are shown ) . ]    [ r_cor ] _ for a given braid @xmath81 obtained from @xmath24 by flipping and smoothing some of its crossings , there is a simple way to check whether @xmath96 by the way @xmath97 is constructed .
we travel through @xmath81 naturally and visit each crossing of @xmath24 exactly twice . for each crossing of @xmath24
that we encounter in this process for the first time ( including the smoothed ones , these are marked by small circles in figure [ fig : cpd ] ) we perform the following test .
if this crossing is smoothed in @xmath81 , we check whether we would be approaching it from its underpass if the crossing ( in its original form in @xmath24 ) were not smoothed . on the other hand ,
if the crossing is also a crossing in @xmath81 ( which may or may not have been flipped from its original form in @xmath24 ) , we check whether we are approaching it from its overpass . if all crossings pass this check , then @xmath98 , otherwise @xmath99 . _
in this section , we derive the main result of this paper , namely the two formulas of the homfly polynomial of a braid @xmath93 based on @xmath94 and @xmath95 respectively , given in theorem [ p3 ] .
let @xmath89 be a vertex in @xmath94 ( or @xmath95 ) , @xmath100 be the writhe of @xmath89 , @xmath101 be the number of components in @xmath102 .
note that @xmath89 is obtained from @xmath24 by applying algorithm d ( or algorithm a ) repeatedly , and in this process some crossings of @xmath24 may have been smoothed .
let @xmath103 be the number of smoothed crossings and @xmath104 be the number of smoothed crossings that are negative .
it is trivial to note that @xmath103 is simply the difference between the number of crossings in @xmath24 and the number of crossings in @xmath89 .
[ l3 ] if @xmath81 is a descending braid on @xmath0 strands , then @xmath105 . on the other hand ,
if @xmath106 is an ascending braid on @xmath0 strands , then @xmath107 .
since @xmath81 is descending , by remark  [ descending_remark ] , @xmath108 can be realized by such space curves that its components are separated by planes that are parallel to the @xmath21-plane , hence the writhe contribution of crossings whose strands belong to different components is zero .
let @xmath109 , @xmath13 , @xmath14 ,
, @xmath18 be the cycles of @xmath110 and @xmath111 be the length of the cycle @xmath112 .
we have @xmath113 and @xmath114 .
we claim that for each component @xmath115 of @xmath108 , @xmath116 . for simplicity
let @xmath117 and let the labels of the strands in @xmath115 be @xmath118 ,
, @xmath119 , ordered by their return order in @xmath81 . without loss of generality , assume that strand @xmath120 is in the plane @xmath121 , strand @xmath122 is in the plane @xmath123 , ... and strand @xmath119 is in the plane @xmath124 .
furthermore , the curve connecting the ending point of @xmath125 to the starting point of strand @xmath126 is bounded between the two planes @xmath127 and @xmath128 for @xmath129 .
we first observe that each strand and part of the curves connecting to its ends can be deformed to straight line segments ( within the plane that it is in ) , in a form as shown in the left hand side of figure [ fig : layers ] , where each strand resides in a plane parallel to @xmath130 ( and the equation of the plane is marked in the figure ) .
the part of the connecting curve attached to the strand residing in the same plane is marked by solid lines and the dotted curves have their end points on different planes parallel to the @xmath21-plane whose @xmath2 coordinates differ by exactly one , with the exception of the dotted curve on the far left ( which is between the planes @xmath124 and @xmath121 ) .     with its strands and
the connecting curves straightened .
right : a special re - arrangement of the strands of @xmath115 through a finite sequence of reidemeister moves of type ii and iii . ]    by remark [ r_remark ] , this deformation does not change the writhe .
for the same reason , these straight line segments can freely slide within the plane they reside in since there are no other curves in that plane .
see figure [ slide ] for one such move .
in particular , nothing can preventing us from sliding them into the position where the strands @xmath131 are parallel lines arranged in this order from left to right , as shown in the right hand side of figure  [ fig : layers ] .
since all moves are made within the planes where these curves reside in , the writhe does not change by remark [ r_remark ] .
we see that @xmath132 from the right hand side of figure  [ fig : layers ] since there are exactly @xmath133 crossings in the projection and all of them are negative .
thus @xmath134 , _
i.e. _ , @xmath105
. an ascending braid diagram @xmath106 is the mirror image of a descending braid diagram @xmath81 .
it is known that @xmath135 and it follows that @xmath107 .
is the only strand in the top plane and can be slid freely within the plane without causing self intersection of the link ; right : the same strand after a sliding move . ]
[ p3 ] let @xmath24 be a braid on @xmath0 strands , @xmath136 and @xmath137 be the set of leaf vertices of @xmath94 and @xmath95 respectively , then @xmath138 @xmath139 where @xmath140 is the number of crossings in @xmath24 that are smoothed in obtaining @xmath81 and @xmath141 is the number of negative crossings among these smoothed crossings .
let us consider the descending tree first . by remark  [ tree_formula_remark ] ,
the contribution of @xmath142 to @xmath143 is @xmath83 multiplied by the weights of the edges on the unique path of @xmath94 from @xmath93 to @xmath108 .
as shown in figure  [ fig : rtree ] , the degree of @xmath1 in the weight of an edge is exactly the change of writhe from the starting vertex of the edge ( remember that it is directed from the root to the leaf ) to the ending vertex of the edge , whereas a @xmath2 term in the weight of the edge indicates that the ending vertex is obtained from the starting vertex by a crossing smoothing and a negative sign in the weight indicates that the smoothed crossing is a negative crossing .
it follows that the total contribution of @xmath81 is @xmath144\end{aligned}\ ] ] by lemma  [ l3 ] .
this proves ( [ e1 ] ) .
equation ( [ e2 ] ) can be proved in a similar fashion and is left to the reader .
_ let @xmath3 be a link with braid index @xmath0 and @xmath93 a braid representation of @xmath3 where @xmath24 is a braid of @xmath0 strands .
let @xmath145 and @xmath146 be the maximum and minimum degrees of @xmath1 in @xmath74 .
the @xmath1-span of @xmath74 is defined as @xmath147 . since @xmath148 for any @xmath142 , formulas ( [ e1 ] ) and ( [ e2 ] ) imply that @xmath149 and @xmath150 .
it follows that @xmath1-span@xmath151 , that is , the morton - frank - williams inequality is a direct consequence of theorem [ p3 ] . _
a link is _ splittable _ if there exist components of the link that lie on different sides of a topological plane in @xmath152 and a link is _ non - splittable _ if it is not splittable  @xcite .
a braid diagram is _ reduced _ if its closure is a reduced link diagram .
a braid diagram is _ alternating _ if its closure is an alternating link diagram .
an oriented link is a _ reduced alternating braid _ on @xmath0 strands if it can be represented by a reduced alternating braid diagram on @xmath0 strands .
a reduced alternating braid is an alternating link but the converse is not true in general .
figure [ fig : bd ] is an example a non - splittable , reduced alternating braid diagram . in this section ,
we prove the following theorem with a simple and direct proof based on our results from the last section .
@xcite[theorem4.1 ] let @xmath24 be the closure of a reduced alternating braid on @xmath0 strands , then the braid index of @xmath24 is @xmath0 .    _
the above theorem is a special case of a more general theorem on a class of oriented alternating fibered links due to murasugi @xcite where the proof relies on using a sequence of lemmas shown by induction .
the alternating links in this class are the @xmath153-products of @xmath154 torus links ( which include the alternating closed braids ) @xcite and are also known to be fibered @xcite .
the fact that the reduced alternating closed braids are fibered can be established directly by using another result due to murasugi , which states that an alternating link is fibered if the leading coefficient of its alexander polynomial is @xmath9 @xcite , and by proving that the leading coefficient of the alexander polynomial of a reduced alternating closed braid is indeed @xmath9 .
this latter fact seems to be known but we failed to find a direct proof of it in the literature so we will provide a short one at the end of this section , which also serves as another application of theorem [ p3 ] and our method .
we also note that fact that the reduced alternating closed braids are fibered is a direct consequence of a result due to stallings in which he proved that all homogeneous closed braids ( which include the alternating ones ) are fibered @xcite . _    before we proceed to the proof of the theorem , we note that it suffices for us to prove this result for reduced alternating braids that are non - splittable .
since if not , say @xmath155 with the @xmath156 s being the non - splittable components of @xmath24 , we can simply apply our result to each @xmath156 .
since the braid index of a link equals the sum of the braid indices of its non - splittable components , and @xmath1-span@xmath5 of @xmath157 is the sum of the ( @xmath1-span@xmath5 ) s of the @xmath158 s as one can easily check , the general result then follows .
let @xmath24 be a braid diagram on @xmath159 strands .
consider @xmath24 as a word of generators of the braid group @xmath31 .
note that each crossing in the braid diagram corresponds to a letter @xmath160 or @xmath161 in the word @xmath24 . in a standard drawing of @xmath24 as a braid ,
a crossing corresponding to @xmath160 or @xmath161 is drawn in the space between the vertical line @xmath162 and @xmath163 in the @xmath21-plane .
we call this space a _ gap_. a gap is _ odd _
if @xmath45 is odd , and _ even _ if @xmath45 is even .
thus , a crossing corresponds to @xmath160 or @xmath161 is in an odd gap ( even gap ) if @xmath45 is odd ( even ) .
if @xmath24 is non - splittable , then each gap in the braid diagram contains at least one crossing .
if @xmath24 is also reduced , then each gap in the diagram contains at least two crossings .
suppose @xmath24 is alternating , the overpass ( underpass ) of a crossing @xmath62 must be the underpass ( overpass ) at next crossing @xmath164 .
so if @xmath164 is in the same gap , @xmath164 must have the same sign as @xmath62 and if @xmath164 is in an adjacent gap , @xmath164 must have the opposite sign . therefore ,
if @xmath24 is alternating , all odd gaps contain crossings of the same sign and all even gaps contain crossings of the opposite sign . if the first gap contains positive crossings , we say the alternating braid diagram @xmath24 is _ positive - leading _ and if the first gap contains negative crossings , we say the alternating braid diagram @xmath24 is _ negative - leading_.    [ l4 ]
let @xmath24 be a non - splittable reduced alternating braid diagram on @xmath0 strands , @xmath145 be the highest degree of @xmath1 in @xmath157 and @xmath146 the lowest degree of @xmath1 in @xmath157 . then @xmath165 and @xmath166 .
it follows that @xmath1-span@xmath167 .
a simple application of the morton - frank - williams inequality and lemma  [ l4 ] , together with the note in the second paragraph of this section , immediately shows that if @xmath24 is a reduced alternating braid on @xmath0 strands , then the braid index of @xmath24 is exactly @xmath0 .
we now proceed to prove lemma [ l4 ] .
we will consider the positive - leading alternating braid diagrams first
. bear in mind that for a positive - leading alternating braid diagram , all the odd gaps contain only positive crossings and all the even gaps contain only negative crossings .
consider @xmath157 as a laurent polynomial of @xmath1 with coefficients in the ring of laurent polynomials in the variable @xmath2 . by
( [ e1 ] ) , the highest possible degree of @xmath1 is @xmath168 and the only leaf vertices @xmath98 that can make contributions to the term of @xmath157 with this @xmath1 degree must satisfy the condition @xmath169 .
let us consider the braid @xmath170 obtained from @xmath24 by first smoothing all crossings in the odd gaps , then smoothing all crossings in even gaps except the first one and the last one and finally flip the sign of the last crossings in the even gaps .
this is obvious by the checking method in remark [ r_cor ] .
part of the @xmath170 contribution to @xmath172 is a term of the form @xmath173 .
proof of claim 2 . by ( [ e1 ] ) , the contribution of @xmath170 to @xmath172 is @xmath174 we have @xmath175 and the result follows .
claim 3 . for any @xmath98 , @xmath176 , the contribution of @xmath81 to @xmath172 either has its maximum degree in the variable @xmath1 less than @xmath168 , or has its degree in @xmath2 less than @xmath177 .
proof of claim 3 .
if @xmath178 then there is nothing to prove . assuming that @xmath169 , then the contribution of @xmath81 to @xmath157 is @xmath179 so the degree of @xmath2 is @xmath180 .
we need to show that @xmath181 , that is , @xmath182 .
since @xmath169 , the return order of @xmath81 is @xmath183 .
it follows that each gap contains either no crossings or at least two crossings of @xmath81 since each strand , say it has label @xmath45 with starting point @xmath55 , has to go through a gap an even number of times in order for it to end at the point @xmath56 .
we claim that @xmath81 contains the first crossing of @xmath24 in each even gap in its original sign .
if this is not the case , let @xmath62 be the first such crossing that has been changed ( so it is either smoothed or flipped in @xmath81 ) .
assume that @xmath62 is in the @xmath45-th gap ( with @xmath45 being even ) .
notice that any strand with label @xmath184 less than @xmath45 can only cross the @xmath45-th gap into gaps to the right of the @xmath45-th gap and must return to the point @xmath185 at the end because the return order of @xmath81 is @xmath183 and the strands can only travel downward .
it follows that the strand entering @xmath62 from its right side must have label greater than @xmath45 .
now consider the strand of @xmath81 that first enters @xmath62 as we travel through @xmath81 naturally .
this strand has to come from the left side of @xmath62 by the above observation .
but then it fails the test given in remark [ r_cor ] since @xmath62 is descending in @xmath24 and can not be changed in @xmath81 .
this gives us the needed contradiction .
it now follows that @xmath81 has at least two crossings in each even gap .
if @xmath81 also contains some crossings in the odd gaps , or contains more than two crossings in some even gaps , then we already have @xmath182 and there is nothing left to prove .
so the only case left is the case when @xmath81 contains no crossings in the odd gaps and exactly two crossings in each even gap .
claim 4 . if @xmath98 and it contains no crossings in the odd gaps and exactly two crossings in each even gap , then @xmath186 , that is , @xmath170 is the only element in @xmath187 with this property .
proof of claim 4 . by the proof of claim 3 , if @xmath188 , then there exists an even gap such that @xmath81 contains the first crossing of @xmath24 in this gap in its original sign , and exactly one other crossing @xmath62 of @xmath24 in this gap which is not the last crossing in this gap . by remark [ r_cor ]
again , the sign of @xmath62 in @xmath24 has to be changed to make it descending in @xmath81 .
but as we travel through @xmath81 naturally past @xmath62 and encounter the first crossing of @xmath24 in the same gap below @xmath62 ( keep in mind that this crossing exists because @xmath62 is not the last crossing of @xmath24 in this gap , and the other crossings of @xmath24 in the adjacent odd gap have all been smoothed in @xmath81 ) .
this crossing has been smoothed in @xmath81 but we are now approaching it from its overpass , so it fails the check in remark [ r_cor ] and @xmath189 , contradicting to @xmath190 .
the consequence of claims 1 to 4 is that if we write @xmath157 as a laurent polynomial of @xmath1 whose coefficients are laurent polynomials of @xmath2 , then it contains a nontrivial term of the form @xmath191 and all other terms have degrees less than @xmath168 , that is , @xmath165 . to obtain @xmath146
, we will use @xmath192 and ( [ e2 ] ) , where @xmath193 is obtained from @xmath24 by keeping the first crossing and flipping the last crossing in each odd gap , and smoothing all other crossings .
the details are left to the reader .
finally , if @xmath24 is a non - splittable reduced negative - leading alternating braid diagram , then its mirror image @xmath194 is a non - splittable reduced positive - leading alternating braid diagram and we have @xmath195 .
let @xmath196 and @xmath197 be the highest and lowest degrees of @xmath1 in @xmath198 .
then by the first part of the lemma , we have @xmath199 and @xmath200 . therefore by remark  [ tree_formula_remark ] we have @xmath201 and @xmath202 .
[ remark4.4 ] _ another way to handle the case when @xmath24 is a non - splittable reduced negative - leading alternating braid diagram is to take its connected sum with a positive hopf link , thus creating a new first gap with a single pair of positive crossings .
the statement now follows the positive leading case and from the fact that the braid index is additive under taking the connected sum .
_    finally , as another application of theorem [ p3 ] , we provide a short proof to the following theorem .
this result seems to be known @xcite , although we failed to find a specific proof in the literature .
[ t4.4 ] let @xmath24 be the closure of a reduced alternating braid on @xmath0 strands , then the leading coefficient in the alexander polynomial of @xmath24 is @xmath9 .
let @xmath203 be the alexander polynomial of @xmath24 , then @xmath204 @xcite .
assume that @xmath24 is positive leading , substituting @xmath205 , @xmath206 in ( [ e1 ] ) leads to @xmath207 construct @xmath208 by the following procedure : ( i ) for the first gap of @xmath24 , smooth all crossings ( which are positive ) except the last one , which we will flip and cross it into the second gap .
once there , we have no choice but to keep the first ( negative ) crossing we encounter which will lead us into the third gap . ( ii ) now we will smooth all crossings we encounter ( negative and positive ones in the second and the third gaps ) until we encounter the last positive crossing in the third gap .
notice that in doing so we may have reached the bottom of the braid at @xmath209 and returned to the top of the braid at @xmath210 , so this last positive crossing we encounter may not be the last positive crossing in the third gap . at this positive
crossing we flip it and cross it into the next gap and repeat ( ii ) .
this process is repeated until we reach the last gap and we can and will smooth all crossings in the last gap except one , no matter it is even or odd . at this point , we have exactly one positive crossing in each odd gap but may have negative crossings in the even gaps that we have not visited .
however in the last gap we only have one crossing left and we will be traveling back through this crossing .
as one can easily verify , each time we travel back to an even gap it is through the only crossing left in the odd gap to its right , so if there are any negative crossings that we have not visited before , these crossings will be between the two positive crossings in the two odd gaps to the left and right of the said even gap . by the descending rule , we can and will smooth all these crossings .
thus we have shown that @xmath208 , @xmath211 and @xmath212 where @xmath213 is the total number of crossings in @xmath24 .
furthermore , it is also rather easy to see that @xmath214 implies @xmath215 .
finally we leave it to our reader to verify that @xmath216 is the only element in @xmath136 with the property @xmath217 .
the result of the theorem then follows .
if @xmath24 is negative leading , we will simply apply the above argument to ( [ e2 ] ) instead .
in this section we show that the homfly polynomial formation ( [ e1 ] ) given in theorem  [ p3 ] is equivalent to the expansion derived by f.  jaeger for braids in @xcite .
thus our approach used to prove proposition  [ p3 ] provides an alternative ( and in fact shorter ) proof of his expansion . unlike our approach ( which is combinatorial in nature ) , jaeger stated his formula using a concept called the _ admissible circuit partitions _ of a braid diagram .
we establish this equivalence by showing that for any braid @xmath24 , there is a bijection between the leaf vertices in the descending resolving tree @xmath94 and the admissible circuit partitions of @xmath24 , such that each leaf vertex and its corresponding admissible partition under this bijection make the same contributions to their respective homfly polynomial formulas .
we begin with a brief review of the essential concepts introduced by f.  jaeger with some small modification of the terminology to make the terms more consistent with our current paper .
interested readers please refer to @xcite for the details and his original terminologies .
given a braid diagram @xmath24 , a circuit partition @xmath218 of @xmath24 is a braid diagram obtained from @xmath24 by smoothing every crossing in a subset @xmath219 of its crossings , while leaving the other crossings unchanged .
the smoothed crossings are denoted by small dotted circles in figure 5 .
we can identify a circuit partition @xmath218 by the ordered pair @xmath220 where @xmath24 is the original braid diagram and @xmath219 the set of crossings of @xmath24 to be smoothed .
two circuit partitions @xmath220 and @xmath221 are defined to be equal if and only if @xmath222 and @xmath223 .
let @xmath62 be a crossing in @xmath24 .
if @xmath62 is smoothed , the upper left part of the strand entering @xmath62 is connected to the lower left part of the strand exiting @xmath62 and the resulting curve is called a _ left tangence _ at @xmath62 .
similarly one can define the _
right tangence _ at @xmath62 .
while traveling through @xmath218 naturally , we will meet each crossing in @xmath24 exactly twice ( including the crossings that have been smoothed ) .
a smoothed crossing @xmath62 in @xmath218 is said to be _
admissible _ if it has the following property .
if @xmath62 is a positive crossing , then the first passage at @xmath62 is a left tangence , if @xmath62 is negative in @xmath24 , then the first passage at @xmath62 is a right tangence .
a circuit partition is _ admissible _ if all the smoothed crossings in @xmath218 are admissible .
in particular , @xmath24 itself is an admissible circuit partition , in which no crossing is smoothed .
we denote the set of admissible circuit partitions of @xmath24 by @xmath224 .
then in @xcite jaeger showed that for any braid diagram @xmath24 on @xmath0 strands , @xmath225 where @xmath226 , @xmath227 is the number of crossings in @xmath219 , @xmath228 is the number of negative crossings in @xmath219 and @xmath229 is the number of components in the closure of @xmath218 .
[ thm : da ] for any braid diagram @xmath24 , there exists a bijection @xmath230 such that for each @xmath142 , @xmath81 and @xmath231 are both obtained from @xmath24 by smoothing crossings from the same set .
define a mapping @xmath232 from @xmath136 to the set of all circuit partitions of @xmath24 as follows . for any @xmath233 ,
let @xmath234 be the circuit partition @xmath220 where @xmath219 is the set of crossings in @xmath24 that are smoothed in the process of obtaining @xmath81 .
we claim that @xmath232 is a bijection between @xmath136 and @xmath224 .
we proceed by induction on @xmath4 , the number of vertices in the descending tree @xmath94 .
if @xmath235 , then the tree consists of only the root vertex , that is , @xmath24 is descending .
so @xmath236 .
assume that there is an admissible circuit partition @xmath237 , and let @xmath62 be the first crossing encountered as we travel through @xmath218 naturally .
up to the crossing @xmath62 , traveling through @xmath218 naturally is the same as traveling through @xmath24 naturally since nothing has been changed up to that point . since @xmath62 is descending , the first strand entering @xmath62 is the top strand at @xmath62 , and smoothing @xmath62 results in a right tangence if @xmath62 is positive and in a left tangence if @xmath62 is negative .
this contradicts the definition of an admissible circuit partition .
thus the only admissible circuit partition of @xmath24 is itself .
so @xmath232 is a bijection .
assume that the statement is true for all @xmath24 whose descending tree has at most @xmath238 vertices and let @xmath24 be such that @xmath94 has at most @xmath239 vertices .
since @xmath240 , @xmath24 contains at least one ascending crossing .
let @xmath62 be the first ascending crossing of @xmath24 encountered when we travel through @xmath24 naturally .
all crossings preceding @xmath62 in the return order being descending , they are not switched or smoothed in any vertex ( braid ) of @xmath94 .
let @xmath241 and @xmath242 be the children of the root vertex of @xmath94 , which are obtained by flipping and smoothing @xmath62 respectively .
the descendants of @xmath241 and @xmath242 respectively form the rooted trees @xmath243 and @xmath244 respectively .
the set of leaves @xmath136 is the union of the sets @xmath245 and @xmath246 .
note that this union is disjoint , since @xmath62 is not smoothed in elements of @xmath245 but is smoothed in elements of @xmath246 . using an argument similar to the one used in the case @xmath235 , we see that for every admissible circuit partition @xmath247 , @xmath219 does not contain any crossing preceding @xmath62 when we travel along the strands of @xmath24 naturally thus each @xmath248 can be identified with @xmath249 if @xmath250 , and with @xmath251 if @xmath252 .
let @xmath253 be the inverse of this identifying map ( which is a bijection , of course ) .
since @xmath241 and @xmath242 are the children of the root vertex of @xmath94 , @xmath254 and @xmath255 both have at least one less vertex than @xmath94 does . by the induction hypothesis , @xmath256 and @xmath257 are bijections that map @xmath258 to @xmath259 and @xmath260 to @xmath261 respectively .
it follows that the mapping @xmath262 defined by @xmath263 if @xmath264 and @xmath265 if @xmath266 is a bijection .
since @xmath267 , @xmath232 is a bijection between @xmath136 and @xmath268 .
this concludes our proof .
_ for the descending tree represented in figure  [ fig : cpd ] , the admissible circuit partitions corresponding to the leaf vertices are depicted on the left hand side of the picture . for each leaf
the corresponding admissible circuit partition is at the same level .
note that for each corresponding pair the same crossings of the braid at the root are smoothed . _    the bijection introduced in proposition  [ thm : da ] has the property that @xmath142 and @xmath269 are identical except the signs at some crossings .
hence we have @xmath270 , @xmath271 and @xmath272 .
this leads to the following result .
f.  jaeger s expansion of the homfly polynomial ( [ j_equation ] )  ( * ? ? ?
* proposition 3 ) is an immediate consequence of ( [ e1 ] ) .    as a final note to this section
, we point out that if we change the definition of admissible circuit partition by defining a smoothed crossing @xmath62 in @xmath218 to be admissible if the first passage at @xmath62 is a left tangence when @xmath62 is negative , and the first passage at @xmath62 is a right tangence if @xmath62 is positive , then we can use ( [ e2 ] ) to show that f.  jaeger s expansion of the homfly polynomial ( [ j_equation ] ) becomes @xmath273 where @xmath274 is the set of admissible circuit partitions under this new definition .
we end this paper by noting the potential application of theorem [ p3 ] to other classes of links that are presented in a closed braid form .
it is also an interesting yet challenging question to explore whether there exist other classes of links that allow formulations similar to the ones in theorem [ p3 ]
. these will be the future research directions of the authors .
the research of the third author was partially supported by a grant from the simons foundation ( # 245153 to gbor hetyei ) .
the authors thank alex stoimenow for his insightful comments .
they are also indebted to an anonymous referee for the careful reading of this manuscript and for suggesting several substantial improvements .
the alternative argument given in remark [ remark4.4 ] is also due to the referee .