since the advent of infrared arrays in the early 1990s , there have been many low resolution spectroscopic atlases published covering the @xmath3 , @xmath4 , or @xmath5 bands .
previous efforts that included the full o  m spectral classes and i  v luminosity classes made use of moderate resolution fourier transform spectroscopy techniques ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
an additional window of atmospheric transparency between these traditional near - infrared bands and the optical window , dubbed the @xmath6 band , was highlighted by @xcite . with existing observational capabilities on instruments at the irtf , keck , las campanas , and ukirt telescopes , and with upcoming wide - field survey facilities such as pan - starrs and vista
, the @xmath6 band is becoming an increasingly appreciated and scientifically important near - infrared bandpass . several low resolution ( @xmath7 ) spectral atlases for m ,
l , and t dwarfs have been published including @xmath6 band data , e.g. @xcite , @xcite , and @xcite . however , a high - resolution spectroscopic atlas covering earlier spectral types and the full range of luminosity classes has not yet appeared .    to complement existing moderate resolution catalogs in the @xmath3 , @xmath4 , and @xmath5 bands , we present here a high - resolution infrared spectral atlas for the @xmath6 band .
our data cover the @xmath8 m region sampled by nirspec @xcite at the w. m. keck observatory .
we obtained spectra of 20 mk - classified stars ranging in spectral type from f through m and in luminosity class from i
v for use as spectroscopic standards in future work .
the data were processed as briefly described in section [ sec : data ] .
we present a representative grid of spectroscopic standards , a table of identified absorption features , and the resulting atlas of spectral lines in section [ sec : atlas ] . in section [ sec : analysis ] , we present relevant atomic data and illustrate as well as discuss line sensitivity to stellar effective temperature and surface gravity .
near - infrared echelle spectra of 20 stars with well - established mk classifications were acquired between 2001 and 2009 with nirspec @xcite on keck ii .
the targets , listed in table [ tab : stds ] , were drawn largely from the lists of standards utilized in the @xcite and @xcite studies of stellar photospheres in the @xmath4 and @xmath3 bands , while the remainder were selected by simbad searches on spectral type and luminosity class , brightness , and accessibility .
we used the n1 filter and have wavelength coverage from @xmath9 to @xmath10 @xmath11 m at r = 25,000 ( @xmath12v = 12 km s@xmath13 ) .
the detector is a 1024 @xmath14 1024 pixel insb chip , and the projected spectrograph slit we used was @xmath15 in angular size .
the @xmath6 band is covered in 12 full and 2 partial orders , without gaps .
the wavelength coverage of each order is given in table [ tab : coverage ] .
spectra were generally taken in nodded pairs using an abba pattern , with the astronomical target offset by @xmath16 in each ab pair .
exposure times ranged from 1 to720 s , integrated in the case of each spectral standard to acquire signal - to - noise ratios @xmath17 in the reduced one - dimensional spectra .
data reduction , including wavelength calibration and spatial rectification , extraction of one - dimensional spectra from the images , and removal of telluric emission and absorption features , is discussed in @xcite and @xcite .
we used a modified version of an iraf script  reduceit " provided by marianne takamiya to reduce the 2002 and 2005 data , and we used the idl package redspec by s. kim , l. prato , and i. mclean to reduce data acquired in 2006 and 2009 .
from our 20 spectroscopic standards we identified a grid representing each of the four spectral classes ( f , g , k , m ) and three luminosity classes ( i - i - i , iii , iv
v ; hereafter supergiants , giants , dwarfs ) .
the star used to typify each class is indicated in table [ tab : stds ] and was chosen based on having a good signal - to - noise ratio and relatively low rotational velocity ( @xmath18 ) .
figures [ fig : grid_80][fig : grid_69 ] show the grid of representative spectroscopic standards displayed from shorter to longer wavelengths for orders 80 through 69 .
each spectrum is normalized to the local continuum and shifted to line rest velocity .
orders 81 and 68 at the edges of our wavelength coverage are excluded due to the low signal - to - noise ratio and the relatively high telluric contamination towards the edges of the @xmath6 band window . as can be seen ,
the 1 @xmath11 m spectral region is rich in spectral features .
the absorption lines originate from both atomic and molecular species .      in order to generate a complete list of the atomic and molecular lines seen in the grid of spectroscopic standards
, we began with a candidate list of transitions from the infrared atlas of the arcturus spectrum @xcite .
arcturus is a k2 iii star and so its associated line list should be relevant for the majority of the cool stars in our sample , with some supplementing needed as described below  especially for the hotter stars . of the 786 lines in our final candidate line list , 209 were added beyond those present in the arcturus atlas list . because the arcturus atlas does not contain the complete details of each atomic transition ( only line frequency and species ) , the candidate line list was cross - referenced with the atomic line database @xcite in order to obtain excitation potentials and log @xmath19 values , as well as the atomic line list @xcite , to obtain the transition terms .
the kurucz line list was developed from theoretical computations and it is far more extensive in line transitions ( over 16,000 within our wavelength range ) than necessary given the moderate spectral resolution and signal - to - noise ratio typical of our ( and most ) empirical astrophysical data .
thus great care was taken to properly match the atomic data to the lines identified in the arcturus atlas .
each line in the arcturus atlas was matched by species , ionization , and wavelength within @xmath20 to lines in the kurucz atomic line database ( in vacuum wavelengths ) . in some cases this resulted in multiple potential matches between the lists among lines of the same species but having different ionization states .
taking ionization into account generally resolved the situation but ambiguities often remained for lines such as the numerous cr i , fe i , and ti i transitions . in such cases we chose to associate  at least initially  the arcturus atlas line with the kurucz line possessing the largest log @xmath19 among those with a realistic ( @xmath2110ev ) lower excitation potential .
after the excitation potential , log @xmath19 , and transition term details of the atomic lines were added from the various sources to the candidate line list , the spectra were visually inspected .
we observed that some lines identified in the arcturus atlas did not correspond to any observed spectral features .
lines that were not present at @xmath22 of the continuum level were flagged appropriately ( as having strength 0 on the scale of 03 described below ) , as were any multiply matched to those in the kurucz list but decided not to be the correct identification .
these lines are ignored in our further analysis but remain in the candidate line list for reference .
visual inspection of the spectra relative to the line list also identified strong absorption features seen in the spectra but not present in the assembled atomic line list . for later type stars ( and especially prominent in the supergiants ) such features were mostly identifiable as molecular lines based again on comparison to the arcturus atlas .
we do not list these individually herein , but we found that nearly all unidentified strong lines in the later type stars appear to be molecular . for the earlier type f stars ,
additional atomic lines were identified using the kurucz line list , now allowing for higher excitation potentials . in the case that there were multiple possible lines within the spectral resolution of the data ( @xmath23 ) the line with lowest ground state energy level was chosen .
if this line was a species that had not been previously identified in the spectra , then other line candidates near the observed line center were considered .    after the initial matching using the arcturus atlas ( hinkle ) and atomic line database ( kurucz ) , all lines that were part of the same multiplet
were identified using the on - line atomic line list ( van hoof ) .
any lines in multiplets of those already identified based on the above exercises and within the spectral range of the data were added to the line list if not already present , and are flagged as so in the final catalog .
as above , we considered only the expected range of log @xmath19 and excitation potential values when initially assigning new candidate lines based on multiplets .
we then also considered the relative line strengths among the multiplets .
specifically , for optically thin ( weak ) lines in local thermodynamic equilibrium , the expected line intensity ratio @xmath24 is @xmath25 or @xmath26 where the symbols have their usual meanings . for stronger lines
the line broadening factors should be taken into account but we have ignored this aspect for ease of analysis . in the case of multiplets ( lines arising from the same lower level ) the above expression simplifies because of the very small energy differences , to @xmath27 which is what we use as a check on the line identifications .
if during the multiplet analysis the case arose that an already - identified line could also be part of a multiplet from a different line species , we chose the multiplet identification over the alternative species unless both could be justified via multiplets or previous identification , in which case an additional flag for line blending was added .
we note that the line wavelengths as well as their other data as given among the arcturus atlas ( hinkle ) , the cfa atomic line database ( kurucz ) and the on - line atomic line list ( van hoof ) were not always identical within quoted errors .
further , there are lines appearing in some of the above sources that are not available in the others .
we believe sensible choices for line identifications have been made , and we have indicated remaining uncertainties in the tabulation .
we also note that in some cases we have been unable to discern the transition term from the existing sources of atomic data .
the 1 @xmath11 m region is a relatively unscrutinized wavelength regime in stellar astrophysics and we are offering some of the first data that can help refine the available atomic line lists and parameters .    finally , all lines in the catalog were assigned a numerical grade based on their relative strength in the spectroscopic standard grid .
the grades range from zero to three , where the zero corresponds to lines that are not seen in any of the spectroscopic standards and are left in the catalog only for completeness , while the three is assigned to the deepest 210 lines per order ( often deeper than 70% of the continuum level across multiple spectral types and luminosity classes ) . throughout this paper
, we consider a one to be a weak line ( shallower than 90% of the continuum ) , while a two ( deeper than 90% of continuum ) or a three ( up to 50% of continuum ) is a strong line .
lines that were blended ( at greater than the 50% level ) , or lines that were only seen in the f standards were also noted in the catalog .
we note that for late - type and especially low - gravity stars , some expected atomic lines can be easily lost amidst strong molecular absorption features , generally cn lines which can exceed 80% depth in some orders .
the number of lines in each strength classification , and the fractions that are molecular , blended , or only seen in the f standards are indicated in table [ tab : strengths ] .    the strongest ( deepest ) atomic and ionic lines ,
those given grades two and three above , plus their related atomic data are presented in table [ tab : lines ] .
the log @xmath19 values among this strong line set range from @xmath28 and up for fe i and tii lines , to -2.5 and up for cr i lines , to @xmath29 and up for si i lines , etc .
we show in figure  [ fig : lowerep ] the lower excitation potential as a function of wavelength and of log @xmath19 for these strong lines .
the general trends in figure  [ fig : lowerep ] are as expected from atomic physics considerations . in figure
[ fig : ephisto ] , we show a histogram of the number of atomic lines in 1ev bins of the lower excitation levels , for both the entire line catalog and for just the strong lines ( those presented in table [ tab : lines ] and figure  [ fig : lowerep ] ) .
this figure is discussed in more detail below .
the table and figures include only the 0.951.11@xmath11 m spectral region , where our final spectra are most uniform in signal - to - noise ratio across the standards grid .
the strong atomic lines identified above as well as the molecular lines are illustrated in figures [ fig : fmgiant_pg1][fig : fmdwarf_pg2 ] using the extremes of our spectral type grid ( f and m , giants and dwarfs ) for comparison .
orders 79 and 71 appear to cover particularly useful spectral ranges for discriminating between hotter and cooler stars ( see section [ sec : diagnostic ] ) based on the trade off in strong features as a function of spectral type .    summarizing , notably present in the 1@xmath11 m spectral region are numerous low and moderate excitation lines of neutral ca i , cr i , fe i , k i , mg i , na i , ni i , si i , and ti i along with singly ionized sr ii . lines of higher excitation potential like h i ( and he i , though not within the effective temperature range sampled here ) , c i , n i , p i , and s i and singly ionized lines of al ii , ca ii , fe ii , and mg ii are also pervasive in the spectroscopic standards . based on visual examination of the spectral grid ,
many of these lines are temperature sensitive and a few , such as sr ii ( @xmath30 ) , are gravity sensitive .
such trends among the atomic lines are discussed more quantitatively below .
the spectra also feature many molecular lines , including ( for all g  k  m types ) a very prominent cn 00 transition molecular series with band heads near 1.093 @xmath11 m ( @xmath31 ) and 1.097 @xmath11 m ( @xmath32 ) according to @xcite ; these occur in the red part of order 70 and the series covers much of order 69 . at the later spectral types ( m stars ) ,
the feh molecule has a prominent 00 band head ( the  wing  ford " band which is apparent longward of @xmath33 ) and @xmath34-branch ( at @xmath35 ) .    at even later spectral types ,
for which we do not have data , the @xmath6 band contains tio bands which appear at 0.942 and 0.950 @xmath11 m , feh appearing at 1.025 and 1.064 @xmath11 m @xcite , and vo appearing at 0.971 and 1.064 @xmath11 m through @xmath36 @xcite .
there is also h@xmath37o toward the short- and long - wavelength edges of the @xmath6 band atmospheric window .
late m and l dwarfs also feature very prominent atomic k i , na i , and ca i features in the @xmath6 band that are only weakly ( grade one on our system ) or not present ( grade zero ) yet in the latest types among our spectra .
we examine in figure  [ fig : ephisto ] the distribution of energy levels for the suite of atomic spectral features identified in the @xmath6 band ( complete candidate line list as well as the strong lines of table  [ tab : lines ] ) .
the number of lines as a function of lower excitation potential is peaked at approximately 5 ev .
the distribution declines more rapidly on the high excitation potential side of the peak than the lower , largely due to the presence of many strong moderate excitation si i and fe i transitions throughout the f , g , k , m spectral range .
heavier metals and higher ionization states generally contribute to the high excitation potential tail , while a large number of ti i and some cr i lines contribute to the lower excitation potential tail .
notable exceptions to this general trend include several c i lines at @xmath38ev , and a few h i lines in the upper excitation potential range .
as expected , many of these high excitation potential lines are only seen ( with strength values of 3 or 2 on our system ) among the earlier type f stars within our spectroscopic grid .    in figure
[ fig : aveprofs ] , we have plotted for the spectroscopic standards grid the average line profile for four different bins in excitation potential .
the atomic lines were divided into bins such that roughly the same number of @xmath6 band lines are contained in each bin , but with an attempt to minimize the splitting of lines of the same species into multiple bins .
this resulted in divisions at @xmath39 , @xmath40 , and @xmath41 . only strong
, unblended lines were included . however , hydrogen lines as well as any lines within @xmath42 of it were excluded since these particularly strong / broad lines dominate in the averaged profiles over nearby line features .
there are 35 lines in the 0 - 2.4 ev bin , 41 lines in the 2.44.9 ev bin , 29 lines in the 4.96.8 ev bin , and 26 lines in the @xmath43 ev bin .
as expected , the hotter stars tend to have stronger lines with higher excitation potentials while cooler stars are dominated by low excitation potential lines .
supergiants also tend to have stronger lines overall relative to the dwarfs .
specifically , the f stars are dominated by the 6.8evand higher lines , the g stars are dominated by the 4.96.8ev lines , and the m stars are dominated by the 02.4ev lines .
one would expect that the k stars would lie in sequence and be dominated by the 2.44.9ev lines .
instead , we find that the line strengths are roughly constant from 0 to 6.8ev and comparable to the 4.96.8ev bin in the g standards .
although the k standards do fit within the overall excitation - potential - temperature trend , they do not match the observed pattern in detail in the expected way .      among the strong lines there are some obvious temperature dependencies apparent in figures  [ fig : fmgiant_pg1]-[fig : fmdwarf_pg2 ] .
for example , in order 79 , there are several strong c i lines that appear in the early spectral types and are not seen in the later types . in the same order
there are also tii lines that complement this trend , appearing very strongly in the late spectral types , but only weakly , if at all , in the early spectral types .
a similar trend can be seen in order 73 for si i relative to ti i , and in order 71 for both c i and si i relative to ti i.    we investigated through equivalent width analysis the behaviors with both effective temperature and surface gravity of the various strong lines in the 1 @xmath11 m region .
several strong lines stand out over the others in table  [ tab : lines ] as potentially useful for analytic or diagnostic purposes , the strongest of which have line depths greater than 50% of the continuum level .
figure  [ fig : equivwidths ] plots the equivalent widths of several such lines versus the standard star effective temperature .
in constructing this plot we were particularly sensitive to including those lines with temperature dependencies discernable for dwarf stars .
as discussed above , lines of note include c i and s i multiplets , which show a strong temperature dependence , and several tii multiplets , which show a strong inverse temperature dependence .
this behavior was also highlighted by for the a5f
z5go multiplet of ti i in the 1.04 - 1.07@xmath11 m region , though we find many of these long wavelength lines of the multiplet blended with other strong lines .
ca i , cr i , and fe i have generally weaker ( inverse ) temperature dependence .
several neutral h i lines and higher excitation potential neutral metals such as c i , p i , and s i and the singly ionized metals ca ii , fe ii , and mg ii appear only in the f standards and are therefore also good temperature indicators .
the si i lines peak in strength in the g  k spectral range and become weaker toward both f and m types , with some gravity sensitivity .
in addition to their temperature dependence , the ti i lines are also gravity sensitive at log @xmath44 3.7 as are the c i and s i lines at log @xmath45 .
srii shows gravity dependence at all temperatures but little trend with temperature .
we have identified a few lines having line depth ratio that are useful indicators of spectral and luminosity class .
these ratios versus standard star effective temperature are illustrated in figure  [ fig : diagnostic ] .
the ti i ( @xmath46 ) and sr ii ( @xmath47 ) lines , are the temperature and gravity sensitive lines , respectively , with the strongest trends apparent in figure  [ fig : equivwidths ] . in this case ,
division of the temperature - sensitive line by the gravity - sensitive line spreads the spectroscopic standards such that dwarfs and then giants have increased line strength relative to supergiants at similar effective temperatures .
the ti i ( @xmath46 ) and c i ( @xmath48 ) lines are both temperature sensitive , but the temperature dependence of c i is opposite that of ti i. this creates an amplified trend in line strength versus temperature relative to that of the single lines ( as in figure  [ fig : equivwidths ] ) .
we also investigated lines pairs that were closer together in wavelength and therefore suitable for temperature and gravity designation from single - order moderate to high dispersion spectra .
thus the cai ( @xmath49 ) @xmath50 sr ii ( @xmath51 ) ratio is also illustrated as sensitive to gravity and the ti i ( @xmath52 ) @xmath50 ci ( @xmath53 ) ratio as sensitive to temperature .
either of these line pair sets could be used in comparing the relative spectral typing of a set of previously unclassified stellar spectra , depending on the spectral dispersion    although not present in our effective temperature range , there is a high excitation ( 21ev ) he i line at @xmath54 in addition to the high excitation potential ( 12 ev ) h i lines already noted here .
these he i and h i lines are useful spectral type diagnostics for o , b , and a stars .
both lines are also particularly sensitive probes at all spectral types of energetic circumstellar processes such as chromospheric activity and wind processes @xcite .
we have presented a spectral atlas for the @xmath6 band region of the near infrared portion of the electromagnetic spectrum .
our data set consists of a grid of mk - classified stars for use as spectroscopic standards spanning the f through m spectral classes and i  v luminosity classes observed at r@xmath5525,000 .
the associated catalog of identified lines contains numerous spectral absorption features , both atomic and molecular in origin , that span a wide range of excitation potentials .
many of the atomic features show strong line depth and equivalent width trends across the spectral atlas . this atlas and the identified atomic spectral lines are potentially useful for two - dimensional spectral classification , especially line depth ratios involving ti i and c i or s i , which are temperature sensitive in opposite directions , and sr ii , which is gravity sensitive .
ascii versions of all spectra are available to download with the online version of the journal ( an example of the form and content of these spectra can be found in table [ tab : asciiex ] ) .
we thank the anonymous referee for his / her careful inspection of the manuscript .
we acknowledge with gratitude the builders of nirspec and the hard work and assistance from the staff at w.m .
keck observatory .
the observatory is operated as a scientific partnership among the california institute of technology , the university of california and the national aeronautics and space administration and was made possible by the generous financial support of the w.m .
keck foundation .
the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of mauna kea has always had within the indigenous hawaiian community .
we are most fortunate to have the opportunity to conduct observations from this mountain .
some of the spectra presented herein were obtained by dan stark and richard ellis or colette salyk and geoff blake ; we are grateful for their contributions .
we consulted the atomic line list ( http://www.pa.uky.edu/@xmath56peter/atomic/ ) and the cfa kurucz line list ( http://cfa-www.harvard.edu/amp/ampdata/kurucz23/sekur.html ) while conducting this research .
we are grateful to andrew baker for his tolerance during the final analysis and writing stages of this paper .
llrclrc hd 21770 & hr 1069 & 2006 nov 30 & f4 iii & 7170 & @xcite & x + hd 193370 & hr 7770 & 2006 nov 30 & f5 ib & 6640 & @xcite & + hd 38232 & bd+29 983 & 2005 dec 13 & f5 ii & 6640 & @xcite & x + hd 55052 & hr 2706 & 2002 nov 26 & f5 iii & 6470 & @xcite & + hd 210027 & hr 8430 & 2006 nov 30 & f5 v & 6530 & @xcite & x + hd21019 & hr 1024 & 2006 nov 29 & g2 v & 5830 & @xcite & x + hd 87 & hr 4 & 2009 jul 11 & g5 iii & 5072 & @xcite & x + hd 58367 & hr 2828 & 2005 dec 13 & g6 iib & 4720 & @xcite & x + hd 75935 & bd+27 1682 & 2006 nov 30 & g8 v & 5430 & @xcite & + hd 166 & hr 8 & 2006 nov 30 & k0 v & 5240 & @xcite & + hd 39400 & hr 2037 & 2005 dec 13 & k2 iib & 4330 & @xcite & x + bd + 59 2667 & & 2006 nov 30 & k2 iv & 5010 & @xcite & + hd 61772 & hr 2959 & 2005 dec 13 & k3 iii &
4320 & @xcite & x + hd 219134 & hr 8832 & 2006 nov 30 & k3 v & 4780 & @xcite & x + hd 201091 & 61 cyga & 2006 dec 01 & k5 v & 4340 & @xcite & + hd 201092 & 61 cygb & 2006 dec 01 & k7 v & 4040 & @xcite & + hd 232979 & bd+52 857 & 2006 nov 30 & m0.5 v & 5240 & @xcite & + hd 39225 & hr 2028 & 2005 dec 13 & m1 + iii & 3550 & @xcite & x + hd 37536 & hr 1939 & 2005 dec 13 & m2 iab & 3450 & @xcite & x + hd 1326 & gl 15a & 2006 nov 30 & m2 v & 3530 & @xcite & x +    lrcl 81 & 0.94118 - 0.94813 + 80 & 0.94646 - 0.95985 + 79 & 0.95835 - 0.97192 + 78 & 0.97054 - 0.98543 + 77 & 0.98304 - 0.99696 + 76 & 0.99590 - 1.00998 + 75 & 1.00906 - 1.02332 + 74 & 1.02255 - 1.03709 + 73 & 1.03628 - 1.05122 + 72 & 1.05081 - 1.06578 + 71 & 1.06553 - 1.08065 + 70 & 1.08065 - 1.09604 + 69 & 1.09584 - 1.11184 + 68 & 1.11241 - 1.12224 +      lcrrrllc 0.9504665 & mn i & -1.200 & 5.69633 & 7.00088 & ? ?
& 2b & uncertain but in the arcturus list + 0.950867 & ti i & 0.078 & 3.58290 & 4.88689&y5fo - f5f & 2b & + 0.951329 & ti i & -0.316 & 3.55869 & 4.86196&y5fo - f5f & 2b & + 0.951587 & fe i & -1.491 & 5.03341 & 6.33641 & x5fo - f5f & 2 & + 0.952265 & n ii & -0.960 & 4.16750 & 5.46958&y3do - e3d & 2 & + 0.954859
& h i & -0.540 & 12.08751 & 13.38597 & 3 - 8 & 3f & + 0.954865 & ti i & -1.598 & 0.83605 & 2.13459&a5f - z5fo & 3 & * + 0.957253 & fe i & -0.542 & 4.99161 & 6.28690&x5fo - e5 g & 2 & + 0.957438 & cr i & -1.784 & 2.54444 & 3.83949&a5g - z5fo & 2b & + 0.957692 & cr i & -1.088 & 2.54479 & 3.83949&a5g - z5fo & 3b & + 0.960218 & ti i & -1.440 & 0.82591 & 2.11721&a5f - z5fo & 3 & * + 0.960567 & c i & -0.940 & 7.48090 & 8.77172&3po-3s & 3f & + 0.962342 & c i & -0.516 & 7.48278 & 8.77114&3po-3s & 3f & + 0.962903 & cr i & -3.451 & 2.91531 & 4.20301&a3p - z3po & 2b & + 0.962913 & fe i & -0.610 & 5.03340 & 6.32108&x5fo - e5 g & 2 & + 0.963682 & fe i & -0.607 & 5.06412 & 6.35078&x5fo - e5 g & 2 & + 0.964090 & ti i & -0.638 & 0.84848 & 2.13459&a5f - z5fo & 3b & only weak blend , * + 0.965002 & ti i & -1.462 & 0.81820 & 2.10309&a5f - z5fo & 3 & * + 0.965222 & s i & 0.350 & 8.41201 & 9.69661&3do-3d & 2f & + 0.965576 & fe i & -0.778 & 4.73346 & 6.01759&y3do - e3f & 2 & + 0.966108 & c i & -0.264 & 7.48830 & 8.77172 & 3po-3s & 3f & + 0.967314 & cr i & -1.230 & 2.54400 & 3.82626&a5g - z5fo & 2 & + 0.967494 & s i & -0.090 & 8.40872 & 9.69031&3do-3d & 2f & + 0.967815 & ti i & -0.831 & 0.83605 & 2.11721&a5f - z5fo & 3 & * + 0.968346 & s i & 0.100 & 8.40982 & 9.69028&3do-3d & 2f & + 0.969133 & ca i & -0.422 & 4.73889 & 6.01830&3do-3p & 2b & + 0.969146 & ti i & -1.638 & 0.81300 & 2.09240&a5f - z5fo & 3 & * + 0.969651 & fe i ] & -3.115 & 5.08593 & 6.36467&x5fo - e7f & 2b & + 0.970242 & fe i ] & -2.215 & 5.06412 & 6.34208&x5fo - e7f & 2 & + 0.970459 & ca i & -0.169 & 4.74385 & 6.02151&3do-3p & 2 & + 0.970834 & ti i & -1.037 & 0.82591 & 2.10309&a5f - z5fo & 3 & * + 0.971964 & ti i & -1.590 & 2.30552 & 3.58121&z5do - a5d & 2 & + 0.972167 & ti i & -1.272 & 1.50259 & 2.77801&a1g - z1fo & 2 & + 0.973100 & ti i & -1.234 & 0.81820 & 2.09240&a5f - z5fo & 3 & * + 0.973719 & cr i & -1.377 & 2.54400 & 3.81740&a5g - z5fo & 2 & + 0.974037 & ti i & -1.882 & 2.29697 & 3.56995&z5do - a5d & 2b & + 0.974124 & fe i & 0.037 & 4.99161 & 6.26447&x5fo - e5 g & 2b & + 0.974625 & ti i & -1.334 & 0.81300 & 2.08521&a5f - z5fo & 3 & * + 0.974956 & ti i & -1.258 & 2.31807 & 3.58984&z5do - a5d & 2b & + 0.976610 & fe i & -0.468 & 5.09989 & 6.36951&x5fo - e5 g & 2b & + 0.976658 & fe i ] & -1.123 & 5.03341 & 6.30297&x5fo - e7f & 2b & + 0.977293 & ti i & -1.604 & 0.84848 & 2.11721&a5f - z5fo & 3 & + 0.978602 & ti i & -1.452 & 0.83605 & 2.10309&a5f - z5fo & 3 & * as blend + 0.978623 & ti i & -1.643 & 0.81820 & 2.08521&a5f - z5fo & 3 & * as blend + 0.979030 & ti i & -1.469 & 0.82591 & 2.09240&a5f - z5fo & 3 & * blend + 0.980352 & fe i ] & -2.223 & 5.09989 & 6.36467&x5fo - e7f & 2 & + 0.985746 & ca ii & -0.228 & 7.505644 & 8.763499&2po-2s & 2f & + 0.986444 & fe i & -0.299 & 5.06412 & 6.32109&x5fo - e5 g & 2f & + 0.989175 & fe i & -0.125 & 5.03341 & 6.28690&x5fo - e5 g & 2 & + 0.989334 & ca ii & -0.601 & 8.438551 & 9.69184&2fo-2 g & 3f & + 0.992711 & fe i ] & -3.095 & 3.54676 & 4.79579&a1d - y3do & 2 & + 0.993006 & ti i & -1.580 & 1.87933 & 3.12799&a3g - y3fo & 2 & + 0.993410 & ca ii & 0.072 & 7.515348 & 8.763499 & 2po-2s & 2f & + 0.994412 & ti i & -1.924 & 2.16043 & 3.40732&a3d - x3do & 2 & + 0.995181 & ti i & -1.831 & 2.15370 & 3.39963&a3d - x3do & 2b & + 0.996401 & na i & -0.820 & 3.61721 & 4.86162&2d-2fo & 2b & doublet , @xmath61 + 1.000032 & fe ii & -1.826 & 5.48450 & 6.72439&z4fo - b4 g & 2f & + 1.000070 & ti i & -1.840 & 1.87329 & 3.11313&a3g - y3fo & 2b & + 1.000580 & ti i & -1.135 & 2.16043 & 3.39963&a3d - x3do & 2 & + 1.001460 & ti i & -1.258 & 2.15370 & 3.39182&a3d - x3do & 2 & + 1.003720 & ti i & -2.002 & 1.46024 & 2.69556&b3f - z3go & 3 & + 1.003940 & sr ii & -1.312 & 1.80482 & 3.03988 & 2d-2po & 3 & + 1.005160 & ti i & -2.131 & 1.44325 & 2.67681&b3f - z3go & 3 & + 1.005210
& h i & -0.303 & 12.08830 & 13.32180 & 3 - 7 & 3f & + 1.006050 & ti i & -0.900 & 2.17485 & 3.40732&a3d - x3do & 2 & + 1.006270 & ti i & -2.278 & 1.42985 & 2.66206&b3f - z3go & 2 & + 1.006780 & fe i & -2.098 & 4.83518 & 6.06676&y3do - e3f & 2 & + 1.009498 & mg ii & 1.020 & 11.63047 & 12.85873&2fo-2 g & 2f & + 1.0115253 & n i & 0.595 & 11.75833 & 12.98412 & 4do-4f & 2f & + 1.011680 & fe i & -3.771 & 2.75876 & 3.98438&a3g - z3fo & 2 & + 1.0117413 & n i & 0.755 & 11.76464 & 12.990176&4do-4f & 2f & + 1.012370 & ti i & -1.672 & 2.17485 & 3.39963&a3d - x3do & 2b & + 1.012660 & c i & -0.032 & 8.53767 & 9.76209 & 1p-1po & 3f & + 1.014830 & fe i & -0.298 & 4.79579 & 6.01759&y3do - e3f & 2 & + 1.015790 & fe i & -4.229 & 2.17609 & 3.39674&a5p - z5fo & 2 & + 1.017030 & fe i & -4.061 & 2.19801 & 3.41718&a5p - z5fo & 2 & + 1.017326 & ti i & -3.292 & 1.44325 & 2.66206&b3f - z3go & 2 & + 1.019192 & ti i & -3.264 & 1.46024 & 2.67681&b3f - z3go & 2 & + 1.019600 & n ii & -0.790 & 4.08956 & 5.30564&1po-5/2[5/2 ] & 2 & + 1.019790 & fe i & -3.699 & 2.72773 & 3.94359&a3g - z3fo & 2 & + 1.021910 & fe i & -0.215 & 4.73346 & 5.94680&y3do - e3f & 2 & * + 1.022120 & fe i & -2.760 & 3.07155 & 4.28464&c3p - z3po & 2 & + 1.026800 & fe i & -4.467 & 2.22286 &
3.43042&a5p - z5fo & 2b & + 1.029180 & si i & -1.360 & 4.92042 & 6.12519 & 3po-3s & 2b & + 1.033010 & sr ii & -0.353 & 1.83958 & 3.03988 & 2d-2po & 3 & * + 1.033310 & n ii & -1.228 & 4.10569 & 5.30564 & 3go-3f &
2 & + 1.034370 & fe i & -3.493 & 2.19801 & 3.39674&a5p - z5fo & 3 & * as blend + 1.034670 & ca i & -0.408 & 2.93271 & 4.13109 & 1po-1s & 3 & * as blend , @xmath61 + 1.037411 & sii & -0.447 & 4.92965 & 6.12519 & 3po-3s & 2 & * + 1.038150 &
n ii & -0.965 & 4.08846 & 5.28282 & ? & 2b & + 1.038184 & fe i & -4.065 & 2.22286 & 3.41718&a5p - z5fo & 2b & + 1.039860 & fe i & -3.330 & 2.17609 & 3.36848&a5p - z5fo & 3b & * as blend + 1.039961 & ti i & -1.537 & 0.84848 & 2.04076&a5f - z5go & 3b & * as blend + 1.042590 & fe i & -3.776 & 2.69259 & 3.88187&a3g - z3fo & 2b & + 1.042660 & fe i & -2.916 & 3.07155 & 4.26074&c3p - z3po & 2b & * + 1.045544 & [ fe i ] & ? & 0.99 & 2.17 & a3f - a5p & 2 & in arcturus list + 1.045830 & s i & 0.260 & 6.86061 & 8.04620 & 3so-3p & 3f & + 1.045960 & s i & -0.430 & 6.86061 & 8.04605 & 3so-3p & 2f & + 1.046230 & s i & 0.040 & 6.86061 & 8.04575 & 3so-3p & 3f & + 1.047250 & fe i & -1.229 & 3.88377 & 5.06775&z3do-3p & 2 & * + 1.048910 & cr i & -0.972 & 3.01083 & 4.19294&b5dz5do & 2 & + 1.049900 & ti i & -1.651 & 0.83605 & 2.01705&a5f - z5go & 3 & + 1.050438 & fe ii & -1.997 & 5.54914 & 6.72953&z4fo - b4 g & 2f & + 1.051290 & cr i & -1.558 & 3.01351 & 4.19294&b5dz5do &
2 & + 1.05324 & p i & 0.200 & 6.95482 & 8.13207 & 4p-4do & 2f & + 1.053340 & n ii & -1.189 & 4.10569 & 5.28282 & ? ?
& 2 & + 1.053510 & fe i & -1.532 & 3.92887 & 5.10582&z3do-3p & 2 & + 1.054410 & c i & -1.290 & 8.53767 & 9.71361 & 1p-1po & 2f & + 1.058000 & fe i & -3.268 & 3.30163 & 4.47357&b3h - z3go & 2 & + 1.058450 & p i & 0.480 & 6.98570 & 8.15716&4p-4do & 2f & + 1.058750 & ti i & -1.777 & 0.82591 & 1.99703&a5f - z5go & 3b & * as blend + 1.058800 & si i & -0.020 & 4.95413 & 6.12519 & 3po-3s & 3b & * as blend , @xmath61 + 1.060630 & si i & -0.380 & 4.92998 & 6.09902 & 3po-3p & 3 & * as blend , @xmath61 + 1.061060 & ti i & -2.676 & 0.84848 & 2.01705&a5f - z5go & 3 & * as blend + 1.06146 & fe i & -0.336 & 6.16964 & 7.33777&n7fo-9/2[13/2 ] & 2 & + 1.06196 & fe i & -3.257 & 3.26733 & 4.43491&b3h - z3go & 2 & + 1.063060 & si i ] & 0.000 & 5.86287 & 7.02925 & 1p-3po & 2 & + 1.063883 & al ii & -0.610 & 15.30358 & 16.46905 & 3fo-3d & 2f & + 1.065050 & cr i & -1.613 & 3.01083 & 4.17502&b5dz5do & 2 & + 1.066390 & si i & -0.350 & 4.92042 & 6.08315 & 3po-3p & 3b & * as blend , @xmath61 + 1.066460 & ti i & -1.916 & 0.81820 & 1.98085&a5f - z5go & 3b & * as blend + 1.067047 & cr i & -1.489 & 3.01301 & 4.17502&b5dz5do & 2 & + 1.067510 & cr i & -1.374 & 3.01351 & 4.17502&b5dz5do & 2 & + 1.068000 & ti i & -2.507 & 0.83605 & 1.99703&a5f - z5go & 3 & + 1.068600 & c i & 0.076 & 7.48328 & 8.64360 & 3po-3d & 3f & + 1.068830 & c i & -0.276 & 7.48090 & 8.64098 & 3po-3d & 3f & + 1.069260 &
si i & -0.190 & 5.95413 & 7.11374 & 3d-3fo & 2 & * + 1.069420 & c i & 0.348 & 7.48830 & 8.64774 & 3po-3d & 3f & + 1.069720 &
si i & -0.060 & 5.96433 & 7.12345 & 3d-3fo & 2 & + 1.071030 & c i & -0.401 & 7.48328 & 8.64098 & 3po-3d & 3f & + 1.072920 &
ti i & -2.064 & 0.81300 & 1.98652&a5f - z5go & 3b & + 1.073030 & si i & 0.060 & 5.98441 & 7.13995 & 3d-3fo & 3b & + 1.073250 & c i & -0.401 & 7.48830 & 8.64360 & 3po-3d & 3f & + 1.073590 &
ti i & -2.504 & 0.82591 & 1.98085&a5f - z5go & 3 & + 1.075230 & si i & -0.300 & 4.92998 & 6.08315 & 3po-3p & 3 & + 1.075595 & fe i & -1.902 & 3.95999 & 5.11277 & z3do-3p & 2 & + 1.075690 & c i & -1.590 & 7.48830 & 8.64098 & 3po-3d & 2f & + 1.077770 & ti i & -2.658 & 0.81814 & 1.96852&a5f - z5go & 3 & + 1.07836 & fe i & -3.414 & 3.23694 & 4.38676 & b3h - z3go & 2 & + 1.078600 & fe i & -2.597 & 3.11117 & 4.26074&c3p - z3po & 2 & + 1.078751 & si i & -0.910 & 5.96433 & 7.11374 & 3d-3fo & 2 & + 1.078980 & si i & -0.390 & 4.92998 & 6.07914 & 3po-3p & 3 & + 1.08043 & cr i & -1.715 & 3.01135 & 4.15897 & b5d - z5do & 2 & + 1.081410 & mg i & -0.320 & 5.94632 & 7.09290 & 3d-3fo & 2b & @xmath61 + 1.082130 & fe i & -1.998 & 3.95999 & 5.10582&z3do-3p & 2b & + 1.082470 & cr i & -1.678 & 3.01351 & 4.15897&b5d - z5do & 2 & + 1.083010 & si i & 0.220 & 4.95413 & 6.09902 & 3po-3p & 3 & + 1.083091
& ti i & -3.878 & 0.83605 & 1.98085&a5f - z5go & 2b & + 1.08378 & na i & -0.500 & 3.61721 & 4.76129&2d-2fo & 2 & @xmath61 + 1.084680 & si i & -0.310 & 5.86287 & 7.00599 & 1p-1do & 2 & + 1.086562 & fe ii & -2.121 & 5.58957 & 6.73071&z4fo - b4 g & 2f & + 1.086650 & fe i ] & -0.803 & 4.73346 & 5.87452&y3do - e5f & 2 & + 1.087180 & si i & 0.040 & 6.19142 & 7.33192&3fo-3/2[9/2 ] & 2b & + 1.087250 & si i & 0.320 & 5.08269 & 6.22311 & 1po-1d & 3b & @xmath61 + 1.08847 & fe i ] & -3.713 & 2.84523 & 3.98438 & b3p - z3fo & 2 & + 1.088578 & si i & -0.870 & 5.98441 & 7.12345 & 3d-3fo & 2 & + 1.08883 & si i & 0.040 & 6.18113 & 7.31990&3fo-3/2[7/2 ] & 2 & + 1.089930 & fe i & -2.729 & 3.07155 & 4.20917&c3p - z3po & 2 & + 1.090886 & cr i & -0.647 & 3.43820 & 4.57480&y7po - e7s & 2 & + 1.091723 & mg ii & 0.020 & 8.86425 & 10.00002 & 2d-2po & 2fb & + 1.091790 & sr ii & -0.638 & 1.80482 & 2.94051 & 2d-2po & 3 & + 1.091827 & mg ii & -0.930 & 8.86436 & 10.00002 & 2d-2po & 2fb & + 1.093287 & cr i & -2.330 & 3.01351 & 4.15897 & b5d - z5do & 2 & amidst cn bands + 1.094110 & h i & 0.002 & 12.08830 & 13.22160 & 3 - 6 & 3f & + 1.095477 & mg ii & -0.230 & 8.86436 & 9.99622 & 2d-2po & 2f & + 1.095630 & mg i & -0.852 & 5.93194 & 7.06364 & 3po-3d & 2b & @xmath61 + 1.096023 & cr i & -2.067 & 3.01351 & 4.15897 & b5d - z5do & 2b & + 1.096030 & mg i & -0.690 & 5.93235 & 7.06364 & 3po-3d & 2b & @xmath61 + 1.096845 & mg i & -0.250 & 5.93319 & 7.06363 & 3po-3d & 2b & cn band head confused , @xmath61 + 1.097038
& p i & 0.120 & 8.07889 & 9.20914 & 2d-2do & 2f & + 1.098231 & sii & -0.590 & 4.95413 & 6.08315 & 3po-3p & 3 & + 1.098284 & n ii & -0.929 & 4.15386 & 5.28282&y3do - e3d & 3f & + 1.09851 & si i & -0.050 & 6.19142 & 7.32015&3fo-3/2[7/2 ] & 2 & + 1.09875 &
si i & -0.500 & 6.19142 & 7.31990&3fo-3/2[7/2 ] & 2 & + 1.10163 & fe i ] & -1.334 & 4.79579 & 5.92133&y3do - e5f & 2b & + 1.101870 & cr i & -0.515 & 3.43820 & 4.57480&y7po - e7s & 2 & + 1.102100 & si i & 0.610 & 6.20642 & 7.33148&3fo-3/2[9/2 ] & 3 & + 1.102287 & ki & -0.010 & 2.67014 & 3.79501 & 2d-2fo & 2 & doublet , wavelength off , @xmath61 + 1.104770 & cr i & -2.508 & 3.01135 & 4.13369&b5d - z5do & 2 & +    cccl 1.0963141 & 1.07008 + 1.0963299 & 1.07359 + 1.0963458 & 1.10119 + 1.0963615 & 1.08921 + 1.0963775 & 1.11744 + 1.0963933 & 1.14190 + 1.0964092 & 1.08628 + 1.0964249 & 1.06581 + 1.0964408 & 1.08024 + 1.0964566 & 1.03010 +                               of the continuum level and appearing in table  [ tab : lines ] ( bottom panel ) .
the lines are divided into @xmath62 bins according to the energy of the lower transition level .
most lines are in the @xmath63 range , but there is a tail extending to higher excitation potential , up to @xmath64 . while the strong lines are more uniformly distributed in the peak excitation potential range , they are still representative of the whole of the line list . [
fig : ephisto ] ]     ( vertical dotted lines with element names ) identified relative to the supergiant spectroscopic standards hd 38232 , an f5 ii ( black line ) and hd 37536 , an m2 iab ( green line ) .
weaker atomic lines from our completed catalog are marked at the top ( short black hashes ) as are molecular lines ( short red hashes ) ( a color version of this figure is available in the online journal . ) .
[ fig : fmgiant_pg1],height=672 ]       ( vertical dotted lines with element names ) identified relative to the dwarf spectroscopic standards hd 55052 , an f5 v ( black line ) and hd 1326 , an m2 v ( green line ) .
weaker atomic lines from our completed catalog are marked at the top ( short black hashes ) as are molecular lines ( short red hashes ) ( a color version of this figure is available in the online journal . ) .
[ fig : fmdwarf_pg1],height=672 ]       ( black ; containing 26 lines ) .
the three h lines are excluded from the calculation of the solid black profile but are shown separately for the f stars as the dotted black line ( in the cooler standards the h profile is too heavily blended with neighboring features for the absorption profile to be informative ) .
only lines that are unblended and strong enough to be tabulated in table  [ tab : lines ] are included .
a clear trend can be seen with hotter stars having stronger high - excitation lines compared to cooler stars having stronger low - excitation lines .
additionally , supergiants tend to have stronger lines then the dwarfs in the same spectral class .
[ fig : aveprofs ] ]     as those of strength three with significant temperature and/or gravity sensitivity within our spectral grid .
generally , only one line from among those in a multiplet set is shown ; its behavior is typical of all lines in the multiplet .
the three columns depict lines with higher excitation potential ( left ) , lines with lower excitation potential ( right ) , and lines sensitive to surface gravity ( center ) . within each column
the panels are ordered from top to bottom alphabetically by line species and then by increasing wavelength .
[ fig : equivwidths ] ]    ) @xmath50 sr ii ( @xmath51 ) , the lines with the strongest temperature and gravity trends in our spectra , while the right side shows ca i ( @xmath49 ) @xmath50 sr ii ( @xmath51 ) which are chosen due to their proximity in wavelength .
bottom panel : the ratio of line strengths for two temperature sensitive lines with opposite temperature dependence , ti i ( @xmath65 ) and c i ( @xmath66 ) on the left side and ti i ( @xmath52 ) and c i ( @xmath53 ) on the right side , as a function of stellar effective temperature .
squares represent supergiants , triangles the giants , and crosses the dwarf standards .
although strong dependences are evident , the illustrated lines connect the existing data only and should not be misinterpreted as defining the detailed trends with temperature .
[ fig : diagnostic ] ]


challenging questions appear often at edges of traditional fields . as an example , the new branch of mathematical physics , the `` statistical topology '' emerged recently by absorbing ideas from the statistical physics , theory of integrable systems , and algebraic topology ( see , @xcite for review ) .
the scope of the statistical topology includes , on the one hand , mathematical problems involved in the construction of topological invariants of knots and links based on solvable models and , on the other hand , the physical and statistical problems related to summation over knot ensembles . in this work , we dwell predominantly to problems of the latter kind , demonstrating the emergence of a critical behavior in ensemble of @xmath0 torus knots .
this critical behavior is formulated in terms of knot invariants .
torus knots @xmath3 seem to be among the simplest objects in the knot theory .
it is difficult to overestimate their role in different branches of mathematical physics .
the topology of a torus knot is uniquely determined by the pair @xmath4 , which fixes windings along two torus periods . in the fig.[fig:01 ] few particular examples of torus knots , @xmath5 from the series @xmath0 are depicted .
the closed curves wrap around the torus which is not shown .    :
@xmath5.,width=302 ]    various explicit expressions for knot invariants are known , ranging from classical jones
kauffman polynomials @xcite , to recent superpolynomials of torus knots @xcite .
new approaches to the construction of torus knot invariants for particular knots , based on the application of topological string theories and deformed matrix models , have been formulated relatively recently in @xcite .
much less is known about properties of _ knot ensembles _ , where the particular topology of a knot diagram is considered as a topologically `` quenched '' variable similar to the quenched disorder in statistical physics .
the weighted summation over different torus knot types ( i.e. different pairs @xmath4 ) means the consideration of the grand canonical ensemble ( i.e. of the generating function ) of torus knots .    in this work
we uncover the relation between three physical problems : i ) two - dimensional directed ( i.e. ( 1 + 1)-dimensional ) lattice random walks with fixed area under the curve , ii ) ensemble of @xmath0 torus knots , and iii ) a five - dimensional sqed with the chern
simons term at the level one .
the reason for random walks to appear in ii ) and iii ) can be intuitively explained as follows .
the main tool for the evaluation of the torus knot superpolynomials @xcite and of nekrasov partition function in the susy gauge theory @xcite , is the `` equivariant localization '' approach , which reduces the integral over the particular moduli space to the summation over the young tableau .
the last problem can be reformulated as a weighted sums over directed paths on a square lattice .
the schematic relations between the problems considered in the paper is shown in the flowchart in the fig.[fig:02 ] .
our main goal is the investigation of the critical behavior in these topological ensembles . from the `` random walk side ''
the critical behavior occurs in the space of fugacities and we shall focus at the `` area+length+corners '' statistics of paths .
we derive an explicit expression for the free energy in the random walks problem for the `` area+length+corners '' statistics , which provides the very nontrivial phase structure in the three - dimensional space of fugacities . in the generic case
we find a cascade of phase transitions and in the degenerate cases we reproduce the previously known results .
the relation between the random walks and the torus knot superpolynomials can be traced from the mathematical literature , however the relation with the particular observable in 5d sqed instanton ensemble is new . collecting different points of view , and the results of @xcite , we interpret the generating function of ensemble of @xmath0 knots as the weighted sum of instanton contributions to the particular observable in 5d sqed and analyze the structure of corresponding generating function .
the interpretation of the random walk fugacities in terms of the generating parameters in the torus knot ensemble is quite straightforward , moreover , these fugacities are identified in the 5d sqed as well : the corresponding parameters turn out to be the gauge coupling , the mass of the hypermultiplet , and the parameters of the @xmath1-deformation .    given an explicit expression for the free energy in the random walks problem for the `` area+length+corners '' statistics , we use it to analyze the ensembles of knots and instantons .
we show that at the `` gauge theory side '' the particular third derivative of the instanton partition function with respect to the masses exhibits an unexpected critical behavior at some critical line in the space of parameters .
we provide a physical interpretation of critical behaviors in all three theories considered here . having an exact expression for the @xmath2-narayana numbers ( the generating function of the area- and corner - weighted ( 1 + 1)d brownian excursion ) , we describe explicitly the phase portrait of different ensembles .
we show that there is a 3rd order phase transition in the instanton ensemble , corresponding to the 1st order phase transition in the ensembles of random paths and torus knots .
the paper is organized as follows . in the section 2
we describe the sum over the paths with different statistics and focus at the `` area+length+corners '' ensembles . in the section 3
we discuss the representation of the partition function of ( 1 + 1)d brownian excursions ( dyck paths ) for the superpolynomials of @xmath0 torus knots .
the section 4 is devoted to the identification of the sum over the paths as the instanton contribution to the particular observable in the 5d sqed . in the section 5 we consider the physical interpretation of the critical behavior in the space of fugacities in all models discussed at length of the paper .
the relation with the toda - like integrable system is mentioned in section 6 .
our findings and the open questions are summarized in the conclusion . in the appendix
a , using the combinatorial description , we derive the new explicit expression for the generating function of @xmath6-narayana polynomials ( the generating function of area- and corner - weighted dyck paths ) . in the appendix b we remind for completeness the relation of the critical behavior with the `` hydrodynamic description '' of the edge singularities in gue matrix ensembles .
by definition , the dyck path of length @xmath7 on the square lattice starts at the origin @xmath8 , ends at point @xmath9 and consists of the union of sequential elementary `` @xmath10 '' and `` @xmath11 '' steps , such that the path always stays above the diagonal of the square  see the fig.[fig:03 ] .
the number of all dyck paths of length @xmath7 is given by the catalan number , @xmath12 .
we also denote dyck paths as `` brownian excursions '' ( be ) , having in mind an image of a charged particle on a square lattice in an external transversal magnetic field ( after applied wick rotation ) , where the motion of a particle is subject to two restrictions : it moves only up and right and never intersects the diagonal . calculating the action for such a particle
, we see that @xmath13 is the fugacity of the area , @xmath14 , under the dyck path , and the @xmath15 is the fugacity of the path length , @xmath16 .
the information about statistics of area - weighted dyck paths can be easily extracted from the generating function , which is the sum over all path lengths .
this model can be referred to as the `` chiral hofstadter system '' , considered in @xcite .
-narayana polynomials.,width=226 ]    to proceed , turn the lattice by @xmath17 and write the recursion relation for the partition function @xmath18 on a half - line , @xmath19 , where @xmath20 is the height of the path at the step @xmath21 .
the area , @xmath14 , below the path is counted as a sum of filled plaquettes ( i.e. `` heights '' ) between the path and the @xmath22-axis , as shown in the fig.[fig:04]a .
each plaquette has the weight @xmath2 .
-step dyck path ( the same figure as in the fig.[fig:03 ] ) turned by @xmath23 , the heights are measured along the dotted lines and the @xmath24-corners have fugacity @xmath25.,width=302 ]    the partition function of area - weighted dyck paths , @xmath18 , satisfies the relation \ {
l z_k+1(x;q ) = q^x-1 z_k(x+1;q ) + z_k(x-1;q ) + z_k(0;q)=0 + z_k=0(x;q)=_x,1    .
[ eq:01 ] equivalently , ( [ eq:01 ] ) can be written in a matrix form .
define @xmath26 . then _
n(q ) = t^n(q ) * z*_0 where t(q)= (    ccccc 0 & 1 & 0 & 0 &  + 1 & 0 & q & 0 &  + 0 & 1 & 0 & q^2 & + 0 & 0 & 1 & 0 & + & & & &    ) ; _ 0= (    c 1 + 0 + 0 +    ) [ eq:02 ] we are interested in the value @xmath27 since at the very last step the trajectory returns to the initial point . evaluating powers of the matrix @xmath28 , we can straightforwardly check that the values of @xmath29 are given by carlitz  riordan @xmath2-catalan numbers , namely , z_n(1;q )
= \ {    ll c_n/2(q ) & + 0 &    .
[ eq:03 ] recall that @xmath30 satisfy the recursion c_n(q ) = _ k=0^n-1q^k
c_k(q ) c_n - k-1(q ) [ eq:04 ] which is the @xmath2-extension of the standard recursion for catalan numbers .
the generating function @xmath31 obeys the functional relation f(s , q ) = 1 + s f(s , q ) f(sq , q ) [ eq:05 ] it is known that the solution of ( [ eq:05 ] ) can be written as a continued fraction expansion , f(s , q ) = = [ eq:06 ] where @xmath32 is the @xmath2-airy function , a_q(s)=_k=0^ ; ( t;q)_k=_k=0^k-1 ( 1-t q^k ) [ eq:07 ]    let us describe how the critical behavior emerges at the brownian excursion side . in the works @xcite it has been shown that in the double scaling limit @xmath33 and @xmath34 the function @xmath35 has the following asymptotic form c(z ) ~c_reg+(1-q)^1/3 ( 4z ) ; z= , [ eq : asymp1 ] where @xmath36 is the regular part at @xmath37 and @xmath38 is the airy function .    the function @xmath39 is the generating function for the undeformed catalan numbers : c(s , q=1)= the generating function @xmath39 is defined for @xmath40 , and at the point @xmath41 the first derivative of @xmath39 experiences a singularity which is interpreted as the critical behavior .
the limit @xmath42 , @xmath43 can be read also from the asymptotic expression for @xmath35 , eq.([eq : asymp1 ] ) : c(s , q)|_q1 ^ -~c_reg-2 .
[ eq : limq1 ] where @xmath36 is the regular part of @xmath35 at @xmath44 and @xmath45 and @xmath46 at @xmath41 .
note that the first non - singular term in ( [ eq : limq1 ] ) does not contain @xmath2 , so it is no matter in which order the limit in ( [ eq : asymp1 ] ) is taken .
however to define the double scaling behavior and derive the airy - type asymptotic , the simultaneous scaling in @xmath47 and @xmath2 is required .    the generating function , @xmath48 , for narayana numbers , which count dyck paths with fixed fugacity of corners , @xmath25 , demonstrates the behavior similar to ( [ eq : limq1 ] ) , namely the square - root singularity , f(a , s)= [ eq : nar ] this behavior at @xmath49 coincides with the one of catalans ( [ eq : limq1 ] ) . in appendix
a we have derived the explicit expression for the generating function @xmath50 of @xmath2-narayanas with the airy - type asymptotic ( [ eq : asymp1 ] ) in terms of @xmath2-orthogonal polynomials related to rogers  ramanujan continued fractions , f(q , a , s)= [ eq : crit ] where @xmath51 is the extension of the @xmath2-airy function @xmath32 defined in ( [ eq:07 ] ) .
the function @xmath51 reads a_q(s;s(1-a ) ) = _ k=0^ [ eq : crit2 ] one can immediately see that @xmath52 , where @xmath53 is given by ( [ eq:07 ] ) .    completing this section it is worth reminding that appearance of the singularity of type ( [ eq : asymp1 ] )
is the manifestation of the third - order phase transition . in the seminal paper @xcite
it has been shown that the largest eigenvalue , @xmath54 , of the gaussian @xmath55 random matrix ensemble , converges at @xmath56 to @xmath57 , where the random variable @xmath58 has a limiting @xmath16-independent distribution , @xmath59 , being the so - called tracy - widom distribution for gue ensemble @xcite .
so , the normalized value @xmath60 at large ( but finite ) @xmath16 has an uncertainty ( i.e. the width of the distribution ) of order of @xmath61 , typical for the 3rd order phase transitions . above and below the critical value @xmath62 , the tails of the distribution @xmath63 have different asymptotics , signifying existence of strong ( for @xmath64 ) and weak ( for @xmath65 ) couplings .
in this section we briefly explain the relation between the invariants of knots and the random walks with three fugacities . we restrict ourselves by the torus knots , @xmath66 , and consider the superpolynomial introduced in @xcite .
they depend on three generating parameters , which can be related to the fugacities of 2d directed random walk .
remind that the superpolynomial is the poincar polynomial of the triply graded khovanov homologies , @xmath67 , which in the vector space can be attributed to the knot .
the knot superpolynomials are the generalizations of the homfly knot polynomials and depend on three variables , corresponding to gradings @xcite p_n , m(a , q , t)=_ijka^i q^j t^k h_ijk at @xmath68 the superpolynomial @xmath69 reduces to the standard homfly polynomial .
alternatively , it can be interpreted as the generating function for the multiplicities in the particular sector of bps states in the susy gauge theories @xcite .
we will be interested in the critical behavior of the `` area+length+corner '' type statistics in the ensemble of torus knots and focus on the particular series of @xmath0 uncolored knots parameterized by one integer , @xmath16 .
the main object , as before , is the generating function for superpolynomials , @xmath70 , in the ensemble of @xmath0 knots , where the fugacity , @xmath47 , is conjugated to the index @xmath16 , which defines the winding around the cycle on the solid torus ( see the fig.[fig:01 ] ) , z(s , a , q , t)=_n=0^ p_n , n+1(a , q , t ) s^n    hopefully there is an explicit expression for the superpolynomial of the @xmath0 torus knots obtained in two different ways .
the first one deals with the combinatorics of the young diagrams and can be related to the brownian excurion approach @xcite , while the second approach has been developed in @xcite via the localization on the fixed points of the torus action in the moduli space of @xmath16 points in @xmath71 .
the second approach will be used in the next section to compare the generating function for superpolyniomials of @xmath0 uncolored torus knots with the instanton contribution to the particular observable in 5d sqed with the chern
simons term .
we focus on the be representation and provide a dictionary identifying the fugacities at be side with the ones at the knot side .
this dictionary establishes the translation of the be language to the language of torus knot superpolynomials .
the construction of torus knot invariants involves two more statistics for dyck paths , the statistics of `` up - corners '' , where the path changes direction from `` up '' ( @xmath10 ) to `` right '' ( @xmath11 ) , and the statistics of `` dinv '' , a definition for which can be found in @xcite .
the superpolynomial for @xmath0 knots expressed in terms of dyck paths can be presented by the following partition function : z_n(q , t , a ) = _ _ ndyck  paths q^a t^#dinv a^#corners , [ eq:1 ] where @xmath2 , @xmath25 and @xmath72 are the fugacities of area , corners and dinv correspondingly .
the main object of our study is the generating function , z(q , t , a , s ) = _
n=0^ z_n(q , t , a ) s^n [ eq:1a ] where s is fugacity for the length of the path conjugated to the integer @xmath16 , which weights knot type ( at a knot side ) and dyck path length ( at a be side ) .
as before , the area , dinv , and corner statistics in the ensemble of lattice paths are represented by the corresponding fugacities in the grand ensemble . to use the known critical behavior for the generating function of superpolynomials in the ensemble of torus knots , we have to consider the reduction of paths statistics to the `` length+area+corners '' one , switching off the `` dinv '' fugacity ( i.e. setting @xmath73 in ( [ eq:1 ] ) ) .
due to @xmath74 duality , one can equivalently consider the `` length+area+corners+bounce '' statistics , switching off the `` bounce '' fugacity .
this duality can be seen , for instance , in the lowest row in the expansion in the @xmath25 variable of the superpolynomial : p_n , n+1(a , q , t)= _ k a^k p^k_n , n+1(q , t ) .
it turns out that the corresponding term at @xmath75 in the expansion coincides with the @xmath76-catalan number , namely , @xmath77 , when the duality is well - known .    in the next section
we explain the relation between the generating function for the torus knot superpolynomials and the some observables in the 5d sqed upon the particular identification of parameters . in the section 6
we shall use the generating function for @xmath2-narayanas to analyze the critical behavior in the torus knots ensembles classified by the homfly knot invariants .
consider now the nekrasov - like partition function @xcite in the abelian 5d susy gauge theory with the massless hypermultiplet in fundamental representation and massive multiplet in the antifundamental representation in the @xmath1background .
the coefficient in front of the chern
simons term is fixed , @xmath78 , the coupling constant in 5d theory is dimensionful , and the fifth coordinate is compact .
the nekrasov partition function is trivial in this theory , however we shall be interested in the vacuum matrix element @xmath79 of the particular operator @xmath80 .
its evaluation involves the weighted sum of the integrals over the instanton moduli space , where the parameters of the @xmath1background , @xmath81 , serve as the equivariant parameters of two torus actions for the integration over the moduli space , @xmath82 , of @xmath16 point - like instantons , all located at the origin .
the instanton number is weighted with the counting parameter @xmath83 , where @xmath84 and @xmath85 is the radius of the compact fifth dimension .
the desired operator @xmath80 can be identified as follows .
first , we have to recognize the deformed catalan numbers , @xmath86 , with @xmath79 in the @xmath16instanton sector .
to this aim we use the following important result @xcite ^t(hilb^n(c^2,0 ) , v ^n v ) = c_n(q , t ) , which interprets the @xmath76-deformed catalans as equivariant integrals over the moduli space of @xmath16abelian instantons valued in the @xmath16th power of the @xmath16dimensional tautological bundle , @xmath87 .
looking at the representation of the @xmath76catalans in terms of the paths on the young tableau , the desired operator , up to the normalization , can be written in the @xmath16instanton sector as o _
n ( tr e^)_n , where @xmath88 is the hypermultiplet , @xmath89 substitutes the integral over the @xmath90 in the @xmath1background , and @xmath91 is the `` long scalar '' in the @xmath1background @xcite . hence , the composite operator under consideration , is the product of local and nonlocal 4-observables .
it is more convenient to use the equivalent , however a bit more symmetric formulation of the desired observable in 5d sqed @xcite .
consider the 5d qed with two flavors in the fundamental representation with masses @xmath92 , @xmath93 , and the flavor in the antifundamental representation with a mass @xmath94 .
the following relation between the second derivative of the nekrasov partition function and the generating function for the @xmath0 knot superpolynomials , holds : @xmath95 this is equivalent to the evaluation of the correlator of local and nonlocal operators defined above due to the relation .
|_m = ( tr e^ ) , we keep the parameter @xmath25 in the superpolynomial arbitrary , hence the complete list of the identifications reads as follows a=- e^-m_a , t = e^-_1 , q = e^-_2 thus , as above , we have the generating function depending on four fugacities .
the relation between the derivative of the nekrasov partition function and torus knots superpolynomials is important per se , however in this paper we are focusing at the critical behavior in the topological ensembles .
therefore , what we need is i ) the identification of the parameters of the gauge theory as the fugacities of the random walk , and ii ) expression of the generating function for the @xmath2-narayana numbers . to make use the known expression for @xmath2-narayana s
, we have to switch off one equivariant parameter , @xmath96 .
the critical behavior is now the @xmath97phase space , and in this formulation it is clear that we are dealing with the 3rd order phase transition . however , the new point is that the 3rd derivative of the nekrasov partition function is taken with respect to three distinct variables @xmath98 ( or , equivalently , to @xmath99 ) . the physical interpretation of the critical behavior is given in the next section .
in this section we discuss the physical interpretation of the critical behavior found above .
as we have already argued , the critical behavior is exact due to the explicit expression of the generating function for @xmath2-narayana s , which depends on three fugacities ( chemical potentials ) controlling length of the path , area under the path and corners .
the interpretations of the fugacities and of the critical behavior are different at brownian excursion side , knot side , and 5d sqed side , so we consider them separately .
however , these seemingly different critical behaviors reflect one and the same generic pattern of the phase transition .      in terms of the ( 1 + 1)d
random walks the critical behavior in the @xmath100 space has the following interpretation . in the limit @xmath101
the fugacity of the area is switched off and we deal with the generating function @xmath48 ( see ( [ eq : nar ] ) ) defined in the @xmath102plane . at large finite lengths , @xmath16 , the average number of corners , @xmath103 diverges with the length , @xmath16 , as @xmath104 .
the corresponding contour plot of the generating function @xmath105 is shown in the fig.[fig:05 ] in the @xmath42panel .
the critical lines divide the @xmath102plane in three domains . in the first domain ( `` phase 1 '' ) since @xmath47 is small , one has short and ( in average ) sufficiently wrinkled paths , with varying number of corners controlled by the fugacity @xmath25 . in the second ( intermediate ) domain
the trajectories are long since they exceed the critical value for the fugacity , again with varying number of corners , however since the generating function @xmath106 diverges , currently we have not any physical interpretation of this phase . in the third domain ( `` phase 2 '' ) the trajectories are long and essentially wrinkled , since for any @xmath47
the value of @xmath25 is bounded from below : @xmath107 .
for @xmath108 the behavior of the system becomes much more rich and the ensemble of area- and corner - weighted dyck paths exhibits the cascade of airy - type ( 3rd order ) phase transitions as it is seen from the contour plots for the function @xmath109 drawn at few fixed values of @xmath2 ( @xmath110 )  see the corresponding panels in the fig.[fig:05 ] .     above the @xmath102plane for different fixed values of the area fugacity , @xmath2 . for @xmath42
the phase 1 corresponds to short paths , and the phase 2  to long and rather wrinkled paths . for
@xmath108 one sees in the @xmath102plane the cascade of transitions with the airy - type asymptotics.,width=302 ]    the asymptotic behavior of the function @xmath50 as a function of @xmath47 for fixed values @xmath2 and @xmath25 near the singularities ( transition points ) is better seen in the figure fig.[fig:06 ] .
one clearly distinguishes the finite cascade of phase transitions at different @xmath111 and @xmath25 .     for few fixed values of @xmath2 and @xmath25.,width=302 ]    for completeness , we provide in the fig.[fig:07 ] two contour plots of the function @xmath112 at two fixed values of corner fugacity , @xmath113 and @xmath114 , which correspond to weakly and highly wrinkled paths .     above the @xmath115 plane for fixed values @xmath116 .
the cascade of transitions in clearly seen.,width=302 ]    the nature of existence of the cascade of transitions lies in the lattice origin of the problem and deals with the incommensurability of `` optimal '' trajectories with respect to the lattice spacing for particular values of @xmath25 and @xmath2 and @xmath47 . actually , it is very easy to imagine that there are conflicts between wrinkled ( with large number of corners ) and `` inflated '' ( with extended area ) trajectories at some specific length . on the lattice this competition leads to a finite number of phase transitions , manifested in the divergence of the generating function @xmath50 .      in the ensemble of the @xmath0
uncolored torus knots , the situation when we can use the exact answer for the @xmath2narayana s corresponds to the unrefined case with 3d space of parameters @xmath100 . hence , we consider generating function for homfly polynomials , @xmath117 , for @xmath118 torus knots .    what could be the interpretation of the phase transitions in the ensembles of the torus knots ?
naively , different behavior of the generating function at different ridges of the transition point , means different ability to recognize a particular knot in the ensemble .
speaking more formally , we could try to distinguish phases via the derivatives of the generating function with respect to the particular generating parameters , that is the average values of the corresponding `` topological numbers '' .
the derivative with respect to the variables yields the average number of winding , or in other terms , the type of the torus knot .
however the interpretation of other topological numbers is less evident . on the other hand ,
these quantum numbers are quite transparent if we remember the that knot invariants are related to the enumeration of bps states . along this approach , other generating parameters count the number of the corresponding branes .
let us remind , that in a pure `` knot framework '' , one can obtain the @xmath119 invariants of the knots upon the substitution a = q^n hence , the @xmath120 limit corresponds to the situation when all @xmath119 invariants are `` degenerate '' and can not be distinguished . at the critical line in the @xmath121 plane in the ensemble of torus knots ,
the average winding @xmath122 diverges and the average value @xmath123 diverges as well . on the critical line they are related as n hence , in the ensemble of torus knots , the stable limit ( @xmath56 ) is reached at the critical line , where we also forget about the `` nonabelian character '' of knot invariants .
note that the counting problem can be formulated in terms of the seifert surface , which for the @xmath0 torus knots reads as x^n = y^n+1 where @xmath124 .
for instance , the alexander polynomials can be formulated in terms of the monodromy in the cohomologies related to the seifert surface .
also , the alexander polynomials can be interpreted in terms of the stratification of the algebra of functions on the seifert surface .
hence , the critical behavior can be formulated in terms of the seifert surface as well . whether one can relate the type of the seifert surface with the singularity type near the transition point of the generating function ,
is still an open question .
let us comment on the interpretation of the cascades of the phase transitions at `` torus knot side '' .
assuming that there are infinite number of phases in the generic situation , it is naturally to conjecture the domination of the different values or windings in different phases .
turn finally to the third part of our construction . in the ensemble of instantons
, the account of three fugacities @xmath100 corresponds to the account of the coupling constant , mass and parameter of the @xmath1deformation in the unrefined limit .
the critical cascades in this case can be derived again from the analytic structure of the @xmath2narayana generating function , @xmath109 .
the limit @xmath44 corresponds to switching off the @xmath1deformation .
hence , the @xmath102phase space corresponds to the space of the mass , @xmath94 , of the hypermultiplet in the antifundamental representation @xcite , and the instanton charge counting parameter , @xmath125 .
the critical line in the physical @xmath102variables is defined by the equation ( 1+q - qe^m_a)^2=4q [ eq : q ] there are two critical lines , @xmath126 and @xmath127 , defined by ( [ eq : q ] ) at the particular relation between the gauge coupling and mass of the antifundamental representation .
configurations with different number of instantons dominate in the different regions of the @xmath102space . at the critical line we have the following scaling behavior    when the parameter of the @xmath1deformation is switched back again , the cascade of the singular surfaces emerges , dividing the three - dimensional parameter space in the regions .
since the parameter @xmath2 corresponds in the 5d sqed to the fugacity for the angular momentum in some plane in @xmath128 , the critical surfaces separate the regions with the different responses of our observable to the rotation . when @xmath75 , the antifundamental becomes massless and we have the critical cascades in @xmath129plane ; their scaling is dictated by the analyitic structure of the @xmath2-airy functions .    in the simplest case when only the @xmath47parameter matters , the critical behavior of the @xmath2catalans corresponds to the following values of parameters in the 5d gauge theory _
1 = 0,_20,s = -8 ^ 2 g^-2 = 42    the general issues concerning the 3rd order phase transitions will discussed in the separate publication @xcite .
let us point out that the scaling function @xmath130 , which appeared many times throughout the text ( see , for instance , ( [ eq : asymp1 ] ) , or ( [ ricatti2 ] ) ) , plays also a very important role , connecting be to the integrable systems .
it is known ( see @xcite and the references therein ) that @xmath131 , defined in ( [ ricatti2 ] ) , is itself a generating function : w(z)|_z ~_k=0^ _ k z^-(3k-1)/2 , [ eq : gen2 ] where the coefficients @xmath132 have well - defined physical sense : representing @xmath132 in the form @xmath133 where @xmath134 is the gamma - function .
one can show @xcite that @xmath135 is the @xmath16th moment of the area under the brownian excursion on the unit interval . defining
@xmath136 one sees that @xmath137 satisfy the recursion k_n = k_n-1+_k=1^n-1k_k k_n - k ; k_0=- .
[ eq : rec3 ] on the other hand , the function @xmath138 appears it the theory of integrable systems as the solution of the rational painlev ii equation ( @xmath139 ) at @xmath140 : u(x)=2u^3(x)+ 4x u(x)+ 4(+ ) .
[ eq : painleve ]    the connection of area - weighted generating function @xmath138 with the rational solutions of painlev ii is not restricted exclusively by ( [ eq : painleve ] ) , and can be pushed for any @xmath141 , where @xmath142 .
take into account that the rational solutions of ( [ eq : painleve ] ) can be written ( see @xcite ) as @xmath143 , where @xmath144 is the @xmath145-function of the toda system ( see @xcite ) , written as a hankel determinant _
n= (    cccc a_0 & a_1 & & a_n-1 + a_1 & a_2 & & a_n + & & & + a_n-1 & a_n & & a_2n-2    ) [ eq : det1 ] and the entries @xmath146 satisfy the recursion @xcite a_n = 2(n-2 ) a_n-3 + _
k=0^n-2 a_k a_n - k-2 , [ eq : rec2 ] with @xmath147 . the associated generating function , g(x , s)=_j=0^a_j(x ) ( -2s)^-j , [ eq : gen1 ] obeys the riccati equation @xcite ( compare to ( [ ricatti ] ) ) -2 g + g^2 + _ s g - ( 4s^2+s^-1)g + 4 x s^2=0 , [ ricatti3 ] whose solution is g(x , s ) = 2s^2 + ( s^2-x ) .
[ eq : airy2 ]    the equation ( [ eq : rec2 ] ) at large @xmath16 resembles ( though being different in details ) the recursion ( [ eq : rec3 ] ) for the function @xmath137 .
the connection between ( [ eq : rec3 ] ) and ( [ eq : rec2 ] ) can be set by comparing ( [ eq : gen2 ] ) and ( [ eq : gen1 ] ) .
finally , we get a_j=(-2)^j_k=0^ _
j=0^ (    c + m    ) ( -x)^m _ 3k+2m+1,j , where @xmath148 is the kronecker @xmath149function .
thus , we explicitly see how the linear combinations of moments of area - weighted brownian excursions , @xmath132 , are connected to the coefficients @xmath150 in the expansion of the toda @xmath145function .
let us also point out the striking similarity of the recursion equation ( [ eq : rec3 ] ) for different momenta of area - weighted dyck paths with the summation over genus , @xmath151 , the partition function of the @xmath152 chern  simons  matter theory , also known as the abjm theory @xcite .
we plan to discuss this question in details in the forthcoming paper @xcite .
in this paper , the consideration of the directed random walks with the `` area+length+corners '' statistics allows us to find a new critical behavior in the ensemble of the @xmath0 torus knots and in the ensemble of instantons in 5d sqed .
this critical behavior takes place in the space of fugacities and the corresponding parameters in the gauge theory .
the result is exact due to the explicit answer of the generating function for paths with these statistics .    in the space of fugacities
there are patterns of critical lines , and cascades of the phase transitions .
the physics of the phase transitions at this critical line for the ensembles of torus knots and instantons is very rich . in both cases at the critical line
the contribution to the quantity under consideration with the infinitely large topological number starts to dominate .
these quantities are windings in the knot case , or topological charges in the 5d sqed . in the gauge theory
we have the critical lines in the ( gauge coupling , mass ) parameter plane .
the phase transition in this plane is physically meaningful and certainly deserved detailed investigation .
even more profound behavior is seen when all three fugacities are nontrivial .
the critical cascades are quite unexpected and physically they take place in the unrefined limit of the @xmath1deformed gauge theory . at the knot side ,
the critical behavior occurs in the generating function of the @xmath0 knots characterized by the homfly polynomials .
the following general comment is in order .
one has to distinguish two `` critical phenomena '' : phase transition and the wall  crossing phenomena .
the latter is typical for the supersymmetric gauge theories when the spectrum of bps states jumps at some curves and surfaces in the moduli space . upon the susy breaking , the curves of marginal stability become the curves of phase transitions . hence literally the wall - crossing should not be considered as a phase transition , since only the spectrum of the stable bps states gets rearranged .
however , there is a pattern of real phase transitions , even in susy theory , with the argyres
douglas superconformal point in sqcd as the notable example . in this case
the particular condensates serve as the order parameter .
do we deal with the phase transition or with the wall - crossing phenomena in our case ?
formally , the 3rd derivative of the nekrasov partition function in 5d sqed diverges , hence we could speak about the 3rd order phase transition in the instanton ensemble . moreover , the uv properties of the condensate play the key role , and critical lines involve the hypermultiplet masses and the gauge coupling . hence , we could speak on some similarity with the argyres
douglas scenario . on the other hand , at the torus knot side , we could remind the interpretation of the knot superpolynomial as the index of the bps states in some theory .
hence the jump in the index is better interpreted in terms of wall - crossing . as far as we know , there is no evident notion of the wall - crossing in terms of random walks , however , we cold think about it in terms of the stokes lines in quantum mechanics .
summarizing , we think that the critical behavior we have described , admits the interpretation both as the wall - crossing , and as the phase transition , depending on the viewpoint .    the critical behavior we have considered in 5d sqed is a nontrivial example of the 3rd order phase transition .
the familiar example of the 3rd order phase transition in the gauge theories is the phase transition of douglas  kazakov type @xcite in 2d pure yang - mills theory on the sphere characterized by the singularity ( [ eq : asymp1 ] ) . somewhat similar phase transition of douglas  kazakov type has been observed in @xcite . in all these situations
the phase transitions are driven by instantons .
on the other hand , our study , as well as @xcite , yield the clear relation between the torus knots and instantons .
this suggests more general viewpoint that there is a knotting process behind the 3rd order phase transitions in the general case and the knot invaraints play the role similar to the central charges of the virasoro algebra in the 2nd order phase transitions .
we shall provide the general discussion of the 3rd order phase transition in the gauge theories , knot ensembles and statistical mechanics elsewhere @xcite .
there are many challenging issues , related to our work , which deserve further study .
for example , it would be interesting to extend the analysis to the whole ensemble of @xmath66 torus knots and the whole set of fugacities .
also , it would be important to recognize the counterparts of the `` bounces '' and `` corners '' in the particle path integral in the continuum ( i.e. off lattice ) and obtain the interpretation of the external magnetic field in the be approach as a kind of the berry curvature from a 4d viewpoint
. it would be also very interesting to extend our consideration and to relate the spectrum of the full ( i.e. non - chiral ) hofstadter model at be side with knot invariants and 4d instantons .
we plan to address all these questions in @xcite .
it seems also very important to recognize the critical behavior observed in this paper in the ensembles of branes @xcite , hopfions @xcite and `` ensemble '' of 3d theories classified by the torus knots @xcite .
another interesting question concerns the physical interpretation of the critical behavior in terms of the extended objects involved in the torus knot description along the discussion in @xcite .
this is the explicit realization of the @xmath153-dual magnetic approach to the evaluation of the knot invariants suggested in @xcite .
we are grateful to e. gorsky , a. milekhin and n. nekrasov for the useful discussions .
the work of a.g . and k.b .
was supported in part by grants rfbr-12 - 02 - 00284 and pics-12 - 02 - 91052 .
the work of k.b . was also supported by the dynasty fellowship program and princeton centennial fellowship .
thanks scgp at stony brook university where the part of this work has been done during the simons workshop on mathematics and physics 2014 for the hospitality and support .
let us associate each entrance of a dyck path at a height @xmath20 with the application of a `` generator '' @xmath154 .
define the semigroup @xmath155 with the infinite set of generators @xmath156 , obeying the commutation relations _
j _ k = _ k _ j |k - j|2 [ a:08 ] any dyck path can be uniquely encoded by a `` word '' written in terms of generators @xmath154 ( @xmath157 ) .
for example , the dyck path in the fig.[fig:03 ] corresponds to a word = _ 0_3_2_1_2_1_3_3 _ 2 _ 1 _ 1_0 _ 3_2 _ 1_1_1_0 [ a:09 ] ( see also the fig.[fig:04 ] ) .
define a `` normal order '' representation of words in the semigroup @xmath155 , which consists in pushing the generators with smaller indices _ as left as possible _ if such a reordering does not violate the commutation relations @xcite .
let us note that all dyck paths written in terms of generators @xmath158 ( @xmath157 ) are automatically normally ordered .
thus , we can compute the partition function of all such dyck paths via the transfer matrix approach , where the transfer matrix , @xmath159 , defines which particular generator , @xmath160 , can stay next to the previous one , @xmath154 :    * if @xmath22 then @xmath161 ; * if @xmath162 then @xmath163 ; * if @xmath164 then @xmath165    graphical representation of the rules ( i)-(ii ) is shown in the fig.[fig:08 ] .     if the previous index is @xmath22 ( i ) , @xmath166 ( ii ) and @xmath167 ( iii).,width=226 ]    in order to count heights , identify a weight @xmath168 with the generator @xmath154
. then the weighted transfer matrix , @xmath169 , consistent with the rules ( i)-(iii ) has the following form r(q)= (    ccccc 1 & q & q^2 & q^3 &  + 1 & q & q^2 & q^3 &  + 0 & q & q^2 & q^3 & + 0 & 0 & q^2 & q^3 & + & & & &    ) [ a:10 ] define @xmath170 .
the number of dyck paths , starting at the height @xmath22 and returning to the height @xmath22 is given by the element @xmath171 of the vector @xmath172 , where w_n(0;q ) = * v * r^n(q ) * v*^ ; = ( 1,0,0 , ... ) [ a:11 ] one can straightforwardly check that w_n(0;q ) = c_n(q)=z_2n(1;q ) [ a:12 ] where @xmath30 is the @xmath173 carlitz  riordan @xmath2-catalan number and @xmath27 is given by ( [ eq:03 ] ) .
thus , the magnetic random walks can be equivalently described either by the transfer matrix @xmath174 ( eq.([eq:02 ] ) ) , or by the transfer matrix @xmath169 ( eq.([a:10 ] ) ) .
the last description can be straightforwardly generalized to account for corners .
the approach developed in the section [ sect:2:2 ] can be easily extended on counting of simultaneously area - weighted and @xmath24-corner - weighted @xmath7-step dyck paths shown in the fig.[fig:04 ] .
specifically , we are interested in calculation of the partition function , @xmath176 , defined in ( [ eq:1 ] ) , where the summation runs over the ensemble of @xmath7-step dyck paths enclosing the area @xmath14 , controlled by the fugacity @xmath2 , and having @xmath93 @xmath24-corners , controlled by the fugacity @xmath25 ( bounces are not counted , so we set @xmath73 in ( [ eq:1 ] ) ) .
note that the @xmath24-corner appears in position @xmath21 if and only if the following condition is fulfilled : @xmath177 . in the opposite situation , i.e. for @xmath178
, the @xmath24-corner is not created . in terms of the transfer matrix approach , which counts the normally ordered words in the semigroup @xmath155
, the @xmath24-corner creation condition means that for all steps which do not decrease the height we should apply the weight @xmath25 ( the fugacity of the corner ) .
thus , we get the following @xmath175-deformed transfer matrix @xmath179 , where we have multiplied all entries of matrix @xmath169 , except the lower sub - diagonal , by a factor @xmath25 : r(q , r)= (    ccccc a & aq & aq^2 & aq^3 &  + 1 & aq & aq^2 & aq^3 &  + 0 & q & aq^2 & aq^3 & + 0 & 0 & q^2 & aq^3 & + & & & &    ) [ a:14 ] ( compare to ( [ a:10 ] ) ) . evaluation of the area- and corner - weighted partition function , @xmath180 is similar to ( [ a:11 ] ) , namely : _ n(q , a)z_2n(q , a)=*v * r^n(q , a ) * v*^;=(1,0,0 , ... ) [ a:15 ]    in particular , for @xmath181 we get the narayana polynomial , which is the generating function for the narayana numbers @xmath182 , _
n(q=1,a ) = _
k=1^k n_n , ka^k =  _ 2f_1(1-n ,- n,2,a)a ; where @xmath183 is the hypergeometric function and _ n , k= (    c n + k    ) (    c n + k-1    ) [ a:16 ] for example , for @xmath184 both eq.([a:15 ] ) ( in which we take @xmath42 ) and eq.([a:16 ] ) give z_2n(q=1,a ) = a + 21 a^2 + 105 a^3 + 175 a^4 + 105 a^5 + 21 a^6 + a^7 meaning that in the ensemble of @xmath185step dyck paths has 21 different configurations with 2 @xmath24-corners , 105 configurations with 3 @xmath24-corners etc .    for any @xmath175
the recursion relation for @xmath186 reads _
n(q , a ) = a n_n-1(q , a)+_k=1^n-1q^k n_k(q , a)n_n - k-1(q , a ) [ a:17 ] ( compare to ( [ eq:04 ] ) ) .
the generating function @xmath187 obeys the functional relation f(q , a , s)=1+(a-1)s f(q , a , s)+s f(q , a , s)f(sq , a , s ) [ a:18 ] ( compare to ( [ eq:05 ] ) ) .
these polynomials coincide with the standard @xmath2-narayana polynomials @xcite . for @xmath42
we reproduce the generating function for narayana numbers @xcite : f(a , s)= [ a:18a ]      consider the recursion relation _
j+1(q , a ) = r(q , a ) * w*_j(q , a ) [ a:19 ] for the vector @xmath188 , where @xmath189 and the transfer - matrix @xmath179 is defined in ( [ a:14 ] ) .
one can easily rewrite ( [ a:19 ] ) for components of the vector @xmath190 : @xmath191 for the generating function , @xmath192 we obtain from ( [ a:20 ] ) s^-1 |f(x ) - s^-1 = q^x-2 |f(x-1 ) + a_y=1^l q^y-1 |f(y ) [ a:21 ] which can be rewritten in the local form    relation looks pretty much similar to the recursion relation for orthogonal polynomials associated with rogers
ramanujan continued fraction
see @xcite :    l u_m+1(z;b , c ) = z(1 + b q^m ) u_m(z;b , c ) + - c q^m-1 u_m-1(z;b , c ) ; ( m > 0 ) + u_0(z;b , c)=1 ; u_1(z;b , c)=z(1+b )    [ a:24 ] the function @xmath193 has the following explicit expression u_m(z;b , c ) = _ k=0^[m/2 ] q^k(k-1 ) [ a:25 ] where @xmath194 is the pochhammer symbol . by setting z=1 ; b= ; c= ; m = x [ a:26 ] we get the following expression for the generating function @xmath195 f(x;q , a , s)=u_x(1 ; , ) [ a:27 ] the function @xmath195 can be represented as a continued fraction of rogers
ramanujan @xcite : @xmath196 where @xmath197 is the extension of the @xmath2-airy function @xmath32 defined in ( [ eq:07 ] ) .
the function @xmath51 reads a_q(s;s(1-a ) ) = _ k=0^ [ a:29 ] at @xmath49 eq.([a:29 ] ) coincides with @xmath32 in eq.([eq:07 ] ) .
in this appendix we shall mention that the asymptotics ( [ eq : asymp1 ] ) describes the scaling of top line in a bunch of directed vicious walks . proceeding as in @xcite ,
take the ensemble of @xmath198 vicious walkers , define the averaged position of the top line and consider its fluctuations near the averaged position .
in such a description all vicious walkers lying below the top line play a role of a `` mean field '' , which pushes the top line to some `` atypical '' equilibrium position , around which it fluctuates .
it is naturally to suppose that the fluctuations of the top line in a mean - field approximation have the same scaling as the fluctuations of the `` inflated '' brownian excursion with fixed area under the path .
one actually can show that , following the line of reasoning of the work @xcite .
the solution of the inviscid burgers equation @xmath199 is @xmath200 and gives the wigner semicircle law centered at the point @xmath201 .
one can smear the function @xmath202 near the boundary value , @xmath203 , adding by hands the gaussian fluctuations , i.e. passing to the burgers equation with a weak diffusivity ( @xmath204 ) , @xmath205 seeking for weakly fluctuating solutions of viscous burgers equation near the top line , ( @xmath206 ) , in the form @xcite @xmath207 and substituting the ansatz for @xmath208 into the viscous burgers equation , one gets the equation for @xmath209 , which for @xmath210 , @xmath211 and appropriate boundary conditions is transformed in the limit @xmath212 into the dimensionless ricatti equation @xcite -yt^-3/2 + w^2 + _ y w=0 , [ ricatti ] having the solution ( for @xmath206 )
w(z ) = 2(2 ^ -1/3z ) z = yn^-1/2 .
[ ricatti2 ] in ( [ ricatti2 ] ) one can recognize the singular part of the grand partition function of the `` area+length''weighted brownian excursion . to make this connection precise , define the partition function @xmath213 of @xmath16-step directed 2d random walk in the upper half - plane of the square lattice ( i.e. the brownian excursion ) with the fixed area , @xmath14 .
thus , one can straightforwardly identify @xmath214 with @xmath215 under the following redefinitions : \ {                        l. carlitz and j. riordan , two element lattice permutation numbers and their @xmath2-generalization , duke j. math . * 31 * , 371 - 388 ( 1964 ) ; j. frlinger and j. hofbauer , @xmath2-catalan numbers , j. comb .
a * 40 * , 248 - 264 ( 1985 )            s. mashkevich and s. ouvry , area distribution of two - dimensional random walks on a square lattice j. stat . phys . * 137 * , 71 - 78 ( 2009 ) ; + s. matveenko and s. ouvry , the area distribution of two - dimensional random walks and non - hermitian hofstadter quantum mechanics , j. phys . a : math
. theor . * 47 * 185001 ( 2014 )                                                        e. witten , fivebranes and knots , arxiv:1101.3216 [ hep - th ] + d. gaiotto and e. witten , knot invariants from four - dimensional gauge theory , adv .
* 16 * , 935 ( 2012 ) arxiv:1106.4789 [ hep - th ]      n. nekrasov and a. okounkov , seiberg - witten theory and random partitions , hep - th/0306238 + l. baulieu , a. losev , and n. nekrasov , chern - simons and twisted supersymmetry in various dimensions , nucl .
b * 522 * , 82 ( 1998 ) [ hep - th/9707174 ]              p.l .
ferrari , m. praehofer , and h. spohn , stochastic growth in one dimension and gaussian multi - matrix models , in proceedings of the _ 14th international congress on mathematical physics _ ( icmp 2003 ) , world scientific ( ed . j .-
zambrini ) ( 2006 ) , 404 - 411 , arxiv : math - ph/0310053    j .- p .
blaizot , m.a .
nowak , large @xmath220 confinement , universal shocks and random matrices , lectures given at the 49 cracow school of theoretical physics , may 31 - june 10 ( 2009 ) zakopane , poland , arxiv:0911.3683