throughout history , the flow of ideas and innovation has led to vast cultural , economic , and political changes .
social scientists have studied this phenomenon in detail in several different settings  @xcite and found that ideas and innovations tend to diffuse along social links .
first , an innovator adopts a novel idea or practice , then people in contact with the innovator adopt it , then people in contact with those people , and so on . in this way , information cascades on a social network .
not surprisingly , information cascades are also common in online social networks .
they are created , for example , when an individual forwards an email she receives to her contacts  @xcite , or retweets a news item to her followers on twitter  @xcite .
understanding how information spreads in online social networks may be indicative of its quality  @xcite . a mathematical tool for analysis of cascades can find extensive use anywhere where
cascades are studied : anomaly and spam detection , information classification , viral marketing , epidemiological studies , computer virus spread , political and social unrest and even power transmission failure  @xcite .
availability of large scale data about human behavior in online social networks has enabled computational scientists to investigate what drives information diffusion and suggest mechanisms to facilitate its spread . however , as in any other field of research , there are two distinct ways of tackling this problem : model - centric or empirical .
model - centric approaches make certain assumptions about how individuals participating in a cascade are affected by their neighbors ( independent cascade or threshold model ) . using these models ,
researchers have tried to infer global properties of information cascades in social networks  @xcite , devise efficient methods to infer the underlying network structure  @xcite or maximize cascade size  @xcite , and identify influential spreaders  @xcite .
however , empirical approaches are needed to validate assumptions made by these models .
we need principled mathematical tools to quantitatively characterize the temporal and spatial properties of cascades as they occur in real - world networks . however , to the best of our knowledge
, no previous work has attempted to quantify the dynamics of information cascades on social networks or characterize their microscopic growth . at most
, researchers have visualized the shape of cascades  @xcite or enumerated their commonly observed patterns  @xcite .
such approaches do not scale to even moderately large cascades .    to address this gap
, we propose a _ practical , general , and scalable quantitative framework _ for the analysis of cascades on social networks that is applicable even to large cascades .
we define a _ cascade generating function _ , which captures the details of the dynamics of information diffusion on networks .
we can use this function to ( 1 ) compute the macroscopic properties of the cascade , such as its size , diameter , average path length , etc .
, ( 2 ) reconstruct the shape of the cascade , and ( 3 ) analyze its microscopic dynamic properties . the cascade generating function is a good _ signature _
@xcite of the contagion process occurring on a network .
it could help us identify patterns , trends , and anomalies within the cascades in near real - time .
it could aid spam filtering , since the flow of spam messages within a network will be different from the flow of valid information .
it could be useful for viral marketing , since it can help us discover the signature of trends that become popular as compared to those which do not .
as the size of cascades grows , storing their complete structure may not be feasible . however , the cascade generating function can approximate the structure of the cascade with very high accuracy , in spite of having pseudo - linear space complexity .
hence , the cascade generating function can provide _ efficient compression _ of the information in a cascade .
this paper makes the following contributions . in section  [ sec : framework ] we describe a general mathematical framework for representing and quantitatively analyzing cascades on social networks .
specifically , in section  [ sec : characterization ] , we define the cascade generating function , which describes how information spreads through the network . we show that this function can be used to compute cascades macroscopic properties , such as its size , diameter , number of paths in the cascade , etc . in section  [ sec : comp ] , we present a fast , efficient algorithm to compute this function , having @xmath0 runtime complexity and @xmath1 space complexity in its naive implementation , where @xmath2 is the number of nodes participating in a cascade , @xmath3 is the maximum degree of any node and @xmath4 is the number of independent cascade seeds .
we demonstrate the use of cascade generating function to study dynamics of cascades in section  [ sec : analysis ] .
we illustrate the framework on simple cascades often observed in online social networks . in section  [ sec : digg ] we also apply it to study large information cascades occurring on a real - world social network of digg ( http://digg.com ) .
this site allows people to submit and vote for news stories , and also to create links to other people in order to see what new stories they have recently voted for .
stories propagate on digg s social network through a series of cascades as users influence their fans to vote for the story  @xcite .
we study the distribution of several macroscopic properties of these cascades .
in addition , we study the microscopic dynamics of their temporal evolution .
time plots of the cascade generating function show several characteristic signatures of cascade growth , such as star - like , chain - like and community - like growth .
[ sec : definitions ]    consider a social network , represented by a graph @xmath5 with @xmath6 nodes and @xmath7 directed edges .
if node @xmath8 wants to watch activities of node @xmath9 , she must create an edge to @xmath9 by designating @xmath9 as a @xmath10 .
we call @xmath8 a _ fan _ ( or follower ) of @xmath9
. figure  [ fig : network](a ) shows a directed network in which node @xmath11 is a fan of @xmath12 and @xmath13 .
we call an edge @xmath14 _ active _ , if node @xmath15 is a fan of node @xmath16 and node @xmath16 is activated before node @xmath15 .
information or influence flows _ from _ activated nodes _ to _ their fans . in the figure above
, information flows from nodes @xmath12 and @xmath13 to @xmath11 .
a _ cascade _ is a sequence of activations generated by a _ contagion process _ , in which nodes cause connected nodes to be activated with some probability . in analogy with the spread of an infectious disease on a network , an _ infected _ ( activated ) node _ exposes _ his fans to the infection .
disease cascades through the network as exposed fans become infected , thereby exposing their own fans to the disease , and so on .
the _ seed _ of a cascade is the node that initiates the cascade . in information
cascades the seed is an independent originator of information , who then influences others to adopt , endorse , or transmit that information .
we call a node that participates in a cascade a _ member _ of the cascade
. a contagion process can generate multiple cascades , and a node can participate in more than one cascade , resulting in a commonly observed `` collision of cascades ''  @xcite phenomenon .
figure  [ fig : network](b ) shows cascades on the network shown in fig .
[ fig : network](a ) , in which nodes are labeled in the order they are activated , with links showing the direction of influence .
as shown , the contagion process generates two cascades whose seeds are nodes @xmath12 and @xmath13 , respectively .
node @xmath11 participates in both cascades .
a _ cascade chain _ is a sequence of connected nodes participating in a cascade .
each node in the cascade chain is influenced by all the nodes in the chain activated before it and influences all the successive nodes in the chain .
the length of the longest chain is the _ diameter _ of the cascade  @xcite .
the _ spread _ of the cascade is the maximal branching number of its participants , i.e. , the maximum number of nodes a single member infects .
the diameter of the contagion process in fig .
[ fig : network](b ) is two ( longest chain is @xmath17 , the spread of cascade 1 ( yellow ) is 4 and of cascade 2 ( in red ) is 2 .
we characterize a cascade mathematically by the _ cascade generating function _ , @xmath18 , which describes how activation spreads through the network .
contagion process is parameterized by the transmission rates @xmath19 @xmath20 $ ] , which give the probability that a node @xmath16 activated at time @xmath21 will activate a connected node @xmath15 at a later time @xmath22 .
though , in principle , @xmath19 could be different for different values of @xmath16 and @xmath15 , for simplicity , we assume that they are all the same , i.e. , @xmath23 . note , that since the nodes are labeled in the temporal order of their activation , @xmath18 characterizes the cascade at time @xmath22 .
we use the contagion process shown in fig .
[ fig : network](b ) to illustrate how the cascade generating function is calculated .
the initial value of the cascade function is some constant . in the example , nodes @xmath12 and @xmath13 are seeds ; therefore , the values of the cascade function at the times they are activated are constant .
while these values may be different , for convenience we set them both to one : @xmath24 .
the value of @xmath25 captures the cumulative effect on node @xmath15 of activated nodes that are connected to @xmath15 .
node @xmath26 is connected to @xmath12 and activated by it with probability @xmath27 ; therefore , @xmath28 . at the time node @xmath11
is activated , cascade function is @xmath29 .
nodes continue to activate others in this fashion . at time @xmath30 ,
the cascade function is @xmath31 .
since @xmath32 only depends on @xmath33 , @xmath34 can be rewritten as @xmath35 .    in general terms , if node @xmath16 is a node activated at time @xmath36 , the value of the cascade generating function at later time @xmath22 when node @xmath15 is activated is : @xmath37 where @xmath38 is a set of nodes connected to node @xmath15 that are activated before it .
since links are directed , without loss of generality , we can assume that there are @xmath39 cascades in a contagion process .
let @xmath40 , @xmath41 , @xmath42 , @xmath43 be the weights of their seeds . then , eq .
[ eq : phi ] reduces to @xmath44 the value of @xmath45 is proportional to the cumulative effect or influence of all cascades on node @xmath15 activated at time @xmath22 and can be described using the vector @xmath46 .
@xmath47 captures the cumulative effect of the cascade generated at seed node @xmath48 on the node @xmath15 where @xmath49 . in fig .
[ fig : network](b ) , at time @xmath50 , @xmath51 @xmath52 where @xmath53 and @xmath54 . at time @xmath55 , @xmath56 . here
@xmath57 and @xmath58 .    if the values of the cascade generating function for nodes @xmath16 and @xmath15 are the same , @xmath59 , the nodes @xmath16 and @xmath15 are _ isomorphic _ with respect to the contagion process .
such nodes are _ structurally similar _ with respect to the cascade ; therefore , the value of the cascade function is independent of the order in which they are activated . by structural similarity , we mean that in a network comprising of only the activated nodes and active edges between them , the topological distance of two isomorphic from all the seeds is the same . here , the topological distance of a node from the seed is measured in terms of the total number of attenuated paths over active edges .
isomorphic nodes can be grouped together in a _ tier _ with its own characteristic @xmath60 . in the contagion process in fig .
[ fig : network](b ) , nodes @xmath26 and @xmath61 are isomorphic and form a tier with value @xmath62
.    * cascade properties .
* we can use the cascade generating function to compute the macroscopic properties of cascades , such as their size , diameter , number of paths , and their average length .    if we take @xmath63 , where @xmath48 is the seed of @xmath64 cascade activated at time @xmath65 , then the total number of paths from @xmath48 to node @xmath15 is equal to @xmath66 in eq .
[ eq : phi1 ]
. the total length of paths from the seed @xmath67 to @xmath15 , @xmath68 , can be obtained by differentiating @xmath25 with respect to @xmath27 and evaluating the derivative at @xmath69 i.e @xmath70 @xmath71 to illustrate this , consider again the contagion process shown in fig .
[ fig : network](b ) . for example
, if we pick node @xmath72 , there are two paths from the seed ( node @xmath12 ) to node @xmath72 : @xmath17 and @xmath73 .
the total length of these paths is three .
there are no paths from the second seed ( node 2 ) to node @xmath72 .
we can also get this answer from @xmath74    we can use similar reasoning to compute other cascade properties . the average path length , @xmath75 is given by :
@xmath76 the diameter of the contagion process is the length of the longest path of any cascade generated by this process .
it is given by modifying eq .
[ eq : av ] : @xmath77}{\frac{d \phi_{min}(j,\alpha)}{d\alpha}\frac { \alpha}{\phi_{min}(j,\alpha ) } } , \label{eq : diam}\ ] ] where @xmath78    [ ! htbp ]    @c|@c@|@c@|@c@|@c|@c@|@c@|@c@| & & & & & & & + &   &   &   &   &   &   &   + ( b ) & & & & &   &   &   + & tot.paths=5&tot.paths=5&tot.paths=31&tot.paths=3&tot.paths=3&tot.paths=3&tot.paths=3 + & tot.len=5&tot.len=15 & tot.len=80 & tot.len=4&tot.len=4 & tot.len=4 & tot.len=5 + & av.len=1&av.len=3 & av.len=2.58&av.len=1.33&av.len=1.33&av.len=1.33 & av.len=1.67 + ( c)&diam.=2 & diam.=5&diam.=5&diam.=2 & diam.=2 & diam.=2 & diam.=2 + & \{1},\{2,3,4,5,6}&\{1},\{2},\{3 } , & & \{1},\{2,3},\{4 } & \{1},\{2,3},\{4 } & \{1},\{2,4},\{3 } & \{1},\{2},\{3,4 } + ( d ) & & \{4},\{5},\{6 } & \{4},\{5},\{6 } & & & & +   + ( a ) &   &   &   &   &   &   &   + ( b ) & & & & &   &   &   + & & & & tot.paths=10 & tot.paths=7 & tot.paths=41 & tot.paths=12 + & tot.len=6 & tot.len=7&tot.len=7&tot.len=15&tot.len=8 & tot.len=154&tot.len=25 + & av.len=1.2&av.len=1.4&av.len=1.4&av.len=1.50&av.len= 1.14&av.len=3.76&av.len=2.08 + ( c ) & diam.=2 & diam.=2&diam.=2&diam.=2&diam.=2 & diam.=8 & diam.=4 + & \{1},\{2,3,4,5}&\{1},\{2,4,6}&\{1},\{2}&\{1},\{2},&\{1,2},\{3,7},&\{1,5,11,12 } , @xmath79 & \{1,3,7},\{2},\{4 } , + ( d ) & \{6 } & \{3,5 } & \{3,4 } & \{3,4,6,7},\{5 } & \{5},\{4},\{6 } & \{8,9},@xmath79,\{18 } & \{5},\{6},\{8 } + & & & & & & & +      [ [ cascade - graph ] ] cascade graph + + + + + + + + + + + + +    for the analysis of the contagion process , we create a _ cascade graph _
@xmath80 from the original network @xmath5 as follows .
let @xmath81 be the number of nodes participating in all cascades .
let a cascade begin at time @xmath82 and end at @xmath83 .
we arrange and label the nodes in the temporal order in which they are activated , e.g. , transmit information : @xmath84 , where node @xmath4 activated at time @xmath85 and @xmath86 .
an edge exists from @xmath15 to @xmath16 in @xmath87 ( i.e. @xmath16 is activated by @xmath15 ) if an edge exists from @xmath16 to @xmath15 in @xmath88 ( @xmath16 is a fan of @xmath15 ) and @xmath89 .
the adjacency matrix of @xmath90 of the cascade graph @xmath80 , the _ cascade matrix _ , is : @xmath91 we break ties randomly . if nodes @xmath8 and @xmath9 receive information at the same time @xmath92 , without loss of generality , we assume @xmath93 and @xmath94 . also , we modify the adjacency matrix @xmath90 , making @xmath95 and @xmath96 , irrespective of whether or not an edge exists between @xmath4 and @xmath97 .
this means that neither node can activate the other , since they are activated at the same time .
we note that @xmath12 is always the seed of a cascade .
the cascade matrix can encode a contagion process that generates multiple cascades .    [ [ contagion - and - length - matrix ] ] contagion and length matrix + + + + + + + + + + + + + + + + + + + + + + + + + + +    in addition to the cascade matrix , we introduce the _ dynamic adjacency matrix _ of the cascade graph , @xmath98 . this is a time - dependent matrix , whose non - zero elements include all nodes that have been activated up to time @xmath99 : @xmath100    the dynamic adjacency matrix allows us to compute connectivity between nodes in a cascade , as measured by the number of paths that exist between them . following @xcite ,
let the attenuation parameter @xmath27 be the probability of transmitting a message or influence along any edge from node @xmath16 at time @xmath36 to node @xmath15 at time @xmath22 .
the _ contagion matrix _ over the time period @xmath101 $ ] : @xmath102 the term @xmath103 gives the number of attenuated paths from node @xmath15 to @xmath16 in @xmath80 and @xmath104 is the identity matrix .
the total length of paths from one node to another can be modeled using a formalism similar to contagion matrix .
we define the _ length matrix _ as : @xmath105 here @xmath106 gives the total length of paths from node @xmath15 to node @xmath16 in @xmath80 .
@xmath107 then gives the _ average length _ of paths from node @xmath15 to node @xmath16 .
the first step towards quantifying cascades is seed identification .
the can be achieved by collecting all the maximal elements of @xmath87 , seen as partially ordered set .
equivalently , if all the elements of the @xmath108 row of @xmath90 are zero , then node @xmath16 is a seed of the cascade .
finding all the seeds gives the total number of cascades , @xmath39 , in the contagion process .
let @xmath109 be the cascade function value of each seed .
the value of the cascade function of node @xmath15 , which is activated after @xmath110 ( but before @xmath111 ) cascade seeds , is @xmath112 .
a non - zero value of @xmath113 indicates that @xmath15 is a member of the cascade initiated by @xmath67 . hence , @xmath114 in eq .
[ eq : phi1 ] .
the efficient design of the cascade generating function is such , that knowing only the columns corresponding to the seeds in the contagion matrix , would help us to characterize the entire matrix .
this along with the triangular nature of the cascade matrix enables us to calculate the contagion and length matrices , and hence the corresponding cascade generating function efficiently using dynamic programming ( algorithm [ alg:2 ] ) .
this algorithm has @xmath1 space and @xmath115 runtime complexity even in its naive implementation , where , @xmath3 is the maximum degree of any node and @xmath2 is the number of nodes in the process .
the _ contagion _ and _ length _ matrices together fully determine @xmath25 and @xmath116 , and therefore , capture the microscopic details of the contagion process .
if vector @xmath117 , @xmath118 , @xmath119 @xmath120 , @xmath121 , @xmath42 , @xmath122 , then @xmath123 and @xmath124 are _ isomorphic _ with respect to the contagion process .
the total number and total length of paths in the cascade from seed @xmath67 to node @xmath15 is given by @xmath125 and @xmath126 in eq .
[ eq : phi1 ] and eq .
[ eq : len ] .
hence , the total number of paths , total length , and average path length for the entire contagion process is given by @xmath127 , @xmath128 and @xmath129 .
as can be seen in eq .
[ eq : diam ] , analogous to the the length matrix , we have devised an efficient algorithm to calculate the diameter .
due to lack of space , we do not provide the algorithm here . since its formalism is very similar to that of the length matrix
, computation has comparable runtime and space complexity .
plotting the cascade generating function @xmath130 _ vs _ time ( @xmath15 ) shows how the structure of the cascade evolves over time .
[ cascade_fig ] illustrates the cascade plots computed for a variety of contagion processes , which include several prototypes of cascades frequently observed in recommendation and blog networks  @xcite .
we label nodes in the order in which they are activated and take @xmath131 when @xmath16 is the seed of the cascade , thus giving equal weights to all cascades in the contagion process . without loss of generality ,
in this study , we set the value of @xmath27 to 0.5 .
future work includes estimation of the transmission rate empirically from the network .
we show that cascade plots contain as much information as cascade graphs , but can be used to analyze the structure and evolution of even large cascades , for which visualization is not feasible .
in addition to showing the cascade plot for each cascade ( row ( a ) ) , fig .
[ cascade_fig ] also reports some of the macroscopic properties of the cascade ( row ( c ) ) , such as total number of paths and their length , average path length , and diameter .
note that this is not the exhaustive list of properties that can be calculated using the cascade characterization function .
row ( d ) lists groups of isomorphic nodes in each cascade .
cascades ( 1)(3 ) in fig .
[ cascade_fig ] are three of the commonly observed patterns , such as a star ( fig .  [ cascade_fig](1 ) ) , a chain ( fig .
[ cascade_fig](2 ) ) , and a community ( clique ) ( fig .
[ cascade_fig](3 ) ) . in the star - like contagion process ,
[ cascade_fig](1 ) , nodes activated by @xmath132 have the same value of @xmath25 , and form an isomorphic group . interchanging the order of their activation does not affect the value of @xmath25 or the cascade plot . in the chain - like contagion process
, cascade function decreases as the chain becomes longer .
there are no isomorphic nodes . in the clique - like contagion process
, the value of the cascade function grows in time as more paths are created in the cascade .
there are also no isomorphic nodes in this cascade .    in the contagion process in fig .
[ cascade_fig](4 ) , nodes activated at @xmath133 and @xmath134 are isomorphic , therefore , the evolution of this cascade is indistinguishable from the cascade shown in fig .
[ cascade_fig](5 ) . however , if the shape of the cascade is the same , but nodes are activated in different order , fig .
[ cascade_fig](6 ) , the cascade plot and its structure are different .
this is because in the contagion processes ( 4 ) and ( 5 ) , cascade widens first ( it is star - like ) , before lengthening , while in the contagion process ( 6 ) , cascade lengthens first ( it is chain - like ) , before widening .
similarly , the cascade ( 7 ) first deepens , then widens , opposite of cascade ( 8) , while cascade ( 9 ) alternates between deepening and widening . in none of these cascades ( except ( 3 ) )
are there multiple paths to a node .
once this happens , as in cascades ( 10 ) and ( 11 ) , the value of the cascade function increases .    we can also disentangle multiple cascades co - occurring in a contagion process .
contagion processes ( 12)(14 ) contain multiple cascades , whose cascade functions are shown in different color .
note that in the contagion process ( 12 ) , node @xmath11 is isomorphic to @xmath26 and @xmath61 with respect to the cascade initiated by @xmath12 , and it is isomorphic to @xmath135 with respect to the cascade initiated by @xmath13 .      given the contagion matrix , it is possible to reconstruct the contagion process with a high level of accuracy .
the cascade generating function @xmath25 , compresses information , and has a space complexity of @xmath136 where @xmath39 is the number of seeds in the contagion process .
how well does the compressed representation @xmath25 capture the contagion process ?    using @xmath25 with @xmath137 , a tier - level reconstruction of the contagion process is possible .
this reconstruction does not remove degeneracy of isomorphic nodes .
temporal ordering of the nodes , help us to fine - tune the tier - level reconstruction .
additional information , such us the number of nodes @xmath110 and their @xmath138 and @xmath139 can help us further improve the approximation . in all the examples shown in fig .
[ cascade_fig ] , using just @xmath25 , we are able to obtain the _ exact tier - level reconstruction _ in all cases . to illustrate , consider fig .
[ fig : network ] , taking @xmath140 , we get @xmath141 , where @xmath27 is the value of the cascade function for the cascade initiated by seed @xmath12 , and @xmath142 is the value of the cascade function for the cascade initiated by second seed , node @xmath13 .
likewise , @xmath143 , @xmath144 and @xmath145 .
hence we can reconstruct that nodes 1 and 2 are independent seeds , 3 and 7 are connected to only node 1 and 5 is connected to only node 2 .
node 4 is connected to both 1 and 2 .
node 6 is connected to 1 and to that _ tier _ of nodes containing node 3 and 7 .
however due to the temporal arrangement of the nodes , we know that 6 is activated before 7 , hence it is necessarily connected to node 3 .
thus we are able to obtain the _ exact _ reconstruction of the cascade .    in fig .
[ cascade_fig ] , using just @xmath25 , we are able to obtain the _ exact tier - level reconstruction _ for cases 1 , 4 , 5 , 7 , 8 , 9,10 , 11 and 13 . in most of these cases
we are also able to disambiguate between isomorphic nodes in the same tier . for cases 2 , 3 , 6 , 12 , and 14
, we are also able to obtain _ exact node - level reconstruction _ of the cascade graph .
[ [ space - and - time - complexity ] ] space and time complexity + + + + + + + + + + + + + + + + + + + + + + + + +    clearly , as demonstrated by the discussion above , knowing the values of @xmath25 at different times allows us to deduce the dynamics of a cascade , and reconstruct its structure ( up to the degeneracy that exists for isomorphic nodes ) .
storing the shape of the cascade has @xmath146 space complexity .
however , as demonstrated above , the cascade generating function can reconstruct this shape with high degree of accuracy . having a pseudo - linear space complexity @xmath136 , it provides an efficient compression of this information .
besides , this model is general , because the same model can be used to investigate cascades in information flow , epidemics , computer viruses , and so on .
this method is fast having @xmath147 runtime complexity even in its naive implementation where @xmath3 is the maximum degree of any node .
moreover , the cascade generating function of a node activated at time @xmath99 , depends only on the cascade value of his friends activated before him .
hence @xmath25 can be calculated real - time and is appropriate even for applications which require streaming , online or near real - time analysis of cascades .
we use the framework described above to study information spread on the social news aggregator digg which allows users to _ post _ and _ vote for _ news stories .
digg users can also create social networks by adding others as friends .
digg highlights the stories a user s friends posted or voted for by marking them with a green ribbon and also displaying them on the _ friends interface _ , a special page for watching friends activity .
a fan may then see the story , and if she decides to vote for it , the story then becomes visible to her own fans , who may in turn vote for it , etc . by voting for a story , a user may influence her fans to also vote it  @xcite .
the spread of a story through the social network of digg is a contagion process that generates many cascades .
the submitter is the _ seed _ of a cascade .
however , there can be other means through which the story can reach a user .
for instance , the user could independently find it on one of digg s web pages or through a link from an external site .
if a user find the story through other means than the friend s interface , she becomes an independent seed for another cascade .
not all seeds , however , generate non - trivial cascades . if a voter is unconnected or does not influence at least one of her fans to vote , the story does not spread .
an independent user who generates a non - trivial cascade is its _ active seed_.    we used digg api to collect data about 3,553 stories promoted to the front page in june 2009 .
the data associated with each story contains its title , i d , link , submitter s name , submission time , list of voters and the time of each vote , and the time the story was promoted to the front page .
in addition , we collected the list of voters friends.lerman/downloads/digg2009.html . ]
we define an _ active user _ as a person who votes for at least one story . in our data set
there are 139,410 active users .
next , we get the connections between the active users .
we say user @xmath8 is @xmath148 to user @xmath9 if she is either a _ friend _ or _ fan _ of user @xmath9 .
we store an active user in _ active users network _ if she is connected to one or more active users .
hence , we are able to determine , whether a user who votes for a particular story is a fan of any previous voter .
out of all active users , 69,524 users @xmath148 to one or more active users .
these 69,524 @xmath148 users form the underlying _
friendship network_. of these , 57,908 users form a giant connected component .
@xmath149 of the @xmath150 distinct submitters belong to this friendship network .
we treat each story as an independent contagion process .
we arrange all voters in the temporal order in which they voted for the story and extract the underlying social network of these voters .
let @xmath151 be the number of _ active _ seeds of the contagion process of a story @xmath152 .
we take each active seed to be independent of other seeds .
therefore , we can quantitatively characterize the cascade by an @xmath153 vector @xmath154 for every node @xmath15 participating in the contagion process .
transmission rate @xmath27 can be derived empirically from the network . in this work , without loss of generality , we set the value of @xmath27 to 0.5 .
we use the framework described above to study the macroscopic properties of cascades on digg , such as the distribution of cascade size , diameter , etc .
we also study the dynamics of evolution of cascades associated with some sample stories .
the stories in our data set generated 216,088 distinct information cascades on the digg social network . using the formalism described above , we calculate global properties of these cascades and plot their distribution .
these properties include cascade size , spread , diameter , etc . due to lack of space ,
we have included in this paper just some examples of the many properties that we can calculate using @xmath25 .    to fit continuous distributions to discrete data , we treat a discrete distribution as if it was generated from a continuous probability density function and then rounded to the nearest integer .
we do not use the commonly used methods such as least square minimization , because the data that spans many orders of magnitude and least square minimization can produce substantially inaccurate estimation of parameters of heavy - tailed distributions like the power - law  @xcite .
we use maximum likelihood parameter estimation ( mle ) to estimate the values of parameters for these distributions and ks statistics to test the goodness of fit .
the closer the ks - statistics to 0 , better the fit .
we study the following distributions : lognormal @xmath155 $ ] , weibull @xmath156 , mixed weibull @xmath157 with @xmath158 , and power - law @xmath159 more often power law applies only for values greater than a certain minimum @xmath160 . in such cases
the _ tail _ of the distribution follows the power law .
using the mle estimates of @xmath160 and scaling parameter @xmath27 , we find what percent of the data comprises this tail of the distribution .
we also investigated distribution fitting using the double pareto lognormal distribution @xcite @xmath161   \\    & - & \frac{1}{\alpha+\beta}(0.5 \beta x^{-\alpha } a(\alpha,\mu,\sigma ) \mbox{erfc}[-\frac{lnx - \mu-\alpha \sigma^2}{\sigma \sqrt{2 } } ] \\   & + & \alpha x^{-\beta } a(-\beta,\mu,\sigma)(1 - 0.5\mbox{erfc}[-\frac{lnx - \mu+\beta \sigma^2}{\sigma \sqrt{2 } } ] )   ) \end{aligned}\ ] ] where @xmath162 .
double pareto lognormal ( dpln ) distribution with @xmath163 , @xmath164 , @xmath165 , @xmath166 gives the best fit for the number of cascades ( better than any of the distributions shown in table [ tbl : mi ] ) with likelihood of @xmath167 and @xmath168 statistic of 0.0109 .
[ htbp ]    fig .
[ fig : cascade_metrics ] shows the distribution of several macroscopic properties of the information cascades on digg , along with functions that best describe them .
table [ tbl : mi ] shows the mle estimates of these distributions .
[ tbl : mi ]    we observe that lognormal or stretched exponential gives a good fit with the observed distribution , and that power law mostly ( if at all ) accounts for a small percentage at the tail of the distribution .
this indicates that a small number of core users may not be driving information propagation in online social networks on the whole .
however , as the cascade size increases , some users may have disproportionate influence on information propagation .
lognormal distribution indicates that the distribution might be generated by a multiplicative effect of many i.i.d random variables .
following the fisher - tippet - gnedenko theorem , the stretched exponential distribution is the limit distribution of properly normalized extrema of a sequence of i.i.d random variables .
hence the distribution may have been generated by the extreme value of the a set of i.i.d random variables .
a very good fit of the distribution of number of cascades with the dpln distribution suggests a possible relationship between the distribution of number of cascades and geometric brownian motion .
future work includes , delving deeper into the probable causes of these distributions .
the cascade generating function @xmath25 is an effective tool not only for computing the global properties of cascades , but also for analyzing their microscopic dynamic signatures . in previous works ,
this was done by visualizing individual cascades  @xcite or by creating a generative model of the contagion process  @xcite .
visualization , however , quickly becomes difficult , even for moderately - sized cascades .
generative models are _ ad hoc _ in nature , and while they are designed to produce cascades with similar macroscopic properties as the observed cascades , they are not guaranteed to reproduce their microscopic characteristics .
the cascade generating function , on the other hand , allows us to study microscopic properties of even very large cascades without the need to visualize them .
[ cols="^,^,^,^ " , ]     we illustrate the use of cascade plots to study microscopic dynamics of cascades with four different stories .
story 1 , titled `` infomercial king billy mays dead at 50 '' was submitted by a user who had 760 fans .
this story was among the most popular in our data set , receiving 8,471 votes , of which 1,244 were from fans .
the contagion process of this story generated 853 cascades .
its diameter was 46 , spread 412 , and the average path length 24 .
[ fig : phi](a ) shows evolution of the cascade function @xmath169 of the top three cascades , ranked by their largest @xmath25 value .
the left - hand set of plots shows the early dynamics of the cascade ( @xmath170 ) , while the right - hand set of plots shows cascade dynamics over the entire time period .
the seed of the cascade is shown in red .
the top cascade attains its largest value of @xmath171 .
this cascade started early in the contagion process .
though the seeds of the next two cascades were also activated within the first 100 votes , these cascades did not start growing until later .
values of @xmath172 imply that the voter is a fan of two or more previous voters .
large values of @xmath25 in fig .
[ fig : phi](a ) indicate a _ community effect _ ( _ cf _ fig .
[ cascade_fig](3 ) ) .
this implies that information is spreading within an interconnected fan network .
though initially the three cascades of story 1 are very different , in their later stages , they become increasingly similar .
this is due to mixing caused by `` collision of cascades , '' which happens when the same nodes participate in different cascades .    the popularity of story 2 , titled `` bender s back , '' is comparable to popularity of story 1 .
story 2 received 8,034 votes of which 1,464 were from fans and generated 722 cascades .
its diameter was 26 , spread 401 , and the average path length 12 .
[ fig : phi](b ) shows both the early and late - stage dynamics of the top three cascades generated by this story .
however , the largest value of @xmath25 attained by any cascade was just @xmath173 , four orders of magnitude smaller than for story 1 .
this indicates a much lower connectivity of the underlying fan network . in the three dominant cascades of this story , @xmath25 does not rise above 2.5 during the first 100 votes .
low values of @xmath25 in the initial stages of cascade evolution imply a _ chaining effect _ ( _ cf _ fig .
[ cascade_fig](1 ) ) , or cascade growth by deepening . unlike story 1 , here the seed of the dominant cascade is the @xmath174 voter .
however , as seen from the larger values of @xmath25 in the cascade plots , in the later stages of information spread , community effect also comes into picture .    the third story in fig .
[ fig : phi](c ) is titled `` play doctor on yourself : 16 things to do between checkups . ''
while this story was submitted by a well - connected user ( with 1,701 fans ) it did not become popular , receiving only 390 votes of which 158 were from fans .
this story generated 11 cascades , and its diameter was 48 , spread 5 , and the average path length 25 .
all of the first 100 voters participated in the dominant cascade , one initiated by the submitter himself .
the maximum @xmath25 value reached by this cascade was very high ( @xmath175 ) , even though this cascade was of short duration .
unlike in previous stories , we observed very high values of @xmath25 already within the first 100 votes , which indicates strong community effect , and high connectivity within the fan network .
for the final illustration we consider the story titled `` apod : 2009 july 1 - three galaxies in draco , '' shown in fig .
[ fig : phi](d ) .
the submitter of this story has only 27 fans .
this story is one of the least popular in our data set , receiving only 199 votes , of which 27 were from fans .
this contagion process generated eight cascades , its diameter was 7 , spread 7 , and the average path length 2.6 .
in the early stages , constant values of @xmath25 in the dominant cascade ( top plot in fig .  [
fig : phi](d ) ) indicate a _ branching effect _ ( _ cf _ fig .
[ cascade_fig](1 ) ) .
this implies that cascade is growing in a star - like fashion , rather than deepening .
the decreasing values of @xmath25 of the third cascade ( bottom plot in fig .
[ fig : phi](d ) ) indicate a _ chaining effect _ , implying that this cascade is deepening .
we do not observe the _ community effect _ either in the initial or later stages of this contagion process .
the maximum value of @xmath25 for this story is two and the average path length is 2.6 , indicating that most of the voters are the fans of the submitter or submitter s fans , but are not themselves interconnected .    in summary , cascade plots can tell us much about the microscopic evolution of information cascade .
popular stories that have large participation also generated many cascades and had high spread .
initially they showed chaining and branching effects , as evidenced by @xmath25 values that are decreasing or staying constant in time , respectively .
the community effect , manifested by growing values of @xmath25 , is visible in later stages when a story penetrates and then spreads through a community .
the trends of the dominant cascades grow increasingly similar with time due to the mixing effect of `` colliding cascades .
'' however stories that do not become popular generate very few cascades and have low spread . when submitter is well connected , the community effect is visible in all stages of the contagion process , implying that the story spreads within submitter s community only . however , when submitter is poorly connected , cascades grow by chaining and branching .
most of the earlier work does not clearly distinguish between cascades and the contagion processes generating these cascades .
we believe that ours is the first work studying large scale cascades without link ambiguity .
though large scale studies of information cascades have been carried out earlier @xcite , the cascades in general were small in size ( o(10 ) ) .
we on the other hand , have very large cascades ( extending up to @xmath176 ) . the quantitative framework for analyzing cascades that we present here
is very scalable and can be easily used to provide a efficient compressed representation of large cascades , when storing the complete information of the entire cascade is no longer trivial . in previous studies
@xcite , the authors characterize the cascades / contagion processes using a multilevel approach comprising of global and local signatures .
cascades are considered to be approximately isomorphic if they have the same global signature .
if the global signatures match , more expensive isomorphism tests based on local signatures are carried out . to aid reasoning about cascades ,
the authors focus on _ local cascades _ , which they define as the ` cascade in the ( undirected ) _ neighborhood _ of the node ' , which for every node is the subgraph induced on the nodes reachable from it .
they enumerate the shapes of these local cascades .
as cascades grow in size , the number of possible shapes increases exponentially and such enumeration becomes infeasible .
note that in this work we provide a scalable , efficient and compressed representation of the observed cascades and make no claims about whether the observed cascades are the actual cascades or fragments of them .    in our study of digg
, we have cascades of size up to @xmath177 .
we aim to deduce their qualitative as well as quantitative properties , such as shape and size .
hence a more formal framework for characterizing cascades is required . in this paper
we provide such a formalism , which not only captures the macro and micro level signatures described above , but much more .
for instance , it captures the similarity between cascades , which were initially similar but later become dissimilar ; or , similarity between cascades that are similar in some stage of their growth .
this formalism enables us to distinguish between cascades and obviates the need for enumerating them or drawing their shape .    in @xcite , the underlying network on which information spreads is not observed , but has to be inferred from the observed cascades .
however such inferences @xcite are based on the hypothesis that the contagion process follows an independent cascade model  @xcite .
our work , on the other hand , focuses on providing a quantitative tool to analyze the trends and patterns of actual contagion processes observed on real - life networks . even when the underlying network is predicted using a different inference methods , e.g. ,  @xcite , the trends of the contagion process occurring on the network can be investigated using the cascade generating function .
future work includes using these tools to aid the verification or rejection of the hypothesis used for modeling information spread  @xcite .
it can also prove to be an effective tool to evaluate the robustness of inferred networks  @xcite .    as demonstrated by the third story in our examples
, we observe that if the submitter is well connected , the community effect is visible at all stages .
however , initial popularity only within the tightly knit community ( shown by a high cascade value and few seeds in the initial stages ) does not ensure global popularity ( large number of votes ) .
in contrast , stories submitted by a not so well connected user , which spreads by branching and deepening initially ( with low cascade values ) , but have larger number of initial active seeds become more popular globally ( as shown by the second story in the example ) .
this observation is in agreement to those reported in @xcite that content diffusing primarily through an interconnected community tends to be confined to that community .
these cascades are also complicated by the interplay between social influence and homophily  @xcite .
future work will address these questions more closely .
in previous works  @xcite , the cascade size was found to be described well by the power - law distribution .
however , we observe that power - law only accounts for a small fraction of cascades at the tail of the distribution .
rather , the entire data can be approximated well with a stretched - exponential ( weibull ) , lognormal or double pareto lognormal distributions , similar to those observed in @xcite .
in this paper we adopt an empirical approach to study cascades on networks .
we believe that our work is first to provide a mathematical framework to quantify and analyze cascades , even for applications requiring real - time or online analysis .
the mathematical framework is based on the cascade generating function , which quantitatively characterizes both the micro and the macroscopic properties of the cascade .
the macroscopic properties that can be efficiently calculated using this tool include the diameter and the spread of the cascades .
this function also provides an efficient compression of the information encoded in cascades . in spite of having pseudo - linear space complexity
, it can be used to reconstruct the shape of the cascade with high degree of accuracy .
although large scale studies of cascades have been carried out , the size of cascades in these studies was relatively small . to the best of our knowledge ,
this is the first study of very large cascades with thousands of participants .
we use this function to study information cascades on an online - social news aggregator digg . for macroscopic properties like number of cascades in a contagion process , cascade size , spread , diameter , average length and so on ,
we observe a stretched exponential ( weibull ) or a lognormal distribution fits well with the observed distribution .
double pareto lognormal distribution gives a very good fit for the distribution of number of cascades .
usually power law accounts ( if at all ) for a small percentage of data in the tail of the distribution .
microscopic analysis also revealed interesting insight to cascades and contagion processes , such as the possible effect of the initial number of seeds and of the branching , chaining and community effect on the initial popularity of news .
10 leskovec j. and mcglohon m. and faloutsos c. and glance n. and hurst m. cascading behavior in large blog graphs , , 2007 .
leskovec j. and singh a. and kleinberg j. patterns of influence in a recommendation network , , 380 - 389 , 2005 .
bailey n. the mathematical theory of infectious diseases and its applications , 1975 .
clauset a. and shalizi c.r and newman m.e.j .
power law distributions in empirical data , , 51,4,661 - 703 , 2009 .
guo l. and tan e. and chen s. and zhang x. and zhao y e. analyzing patterns of user content generation in online social networks , , 369 - 378 , 2009 .
bakshy e. and karrer b. and adamic l. a. social influence and diffusion of user - generated content , , 2009 .
anagnostopoulos a. and kumar r. and mahdian m. influence and correlation in social networks , , 2009 .
dodds p. s. and watts d. j. universal behavior in a generalized model of contagion , , 2004 .
watts d. j. a simple model of global cascades in random networks , , 2002 .
choudhury m. d. and sundaram h. and john a. and seligmann d. d. and kelliher a. `` birds of a feather '' : does homophily among users impact information diffusion in social media ?
gruhl d. and guha r. and liben - nowell d. and tomkins a. information diffusion through blogspace , , 2004 .
ghosh r. and lerman k. predicting influentials in online social networks , , 2010 .
rodriguez m. g. and leskovec j. and krause a. inferring networks of diffusion and influence , , 2010 .
kitsak m. and gallos l. k. and havlin s. and liljeros f.and muchnik l. and stanley h. e. and makse h. a. identifying influential spreaders in complex networks , , 2010 .
rogers e. m. diffusion of innovations , 2003 .
wu f. and huberman b. and adamic l. and tyler j. information flow in social groups , , 2004 .
leskovec j. and backstrom l. and kleinberg j. meme - tracking and the dynamics of the news cycle , , 2009 .
lerman k. and ghosh r. information contagion : an empirical study of the spread of news on digg and twitter social networks , , 2010 .
crane r. and sornette d. viral , quality , and junk videos on youtube : separating content from noise in an information - rich environment , , 2008 .
lerman k. and galstyan a. analysis of social voting patterns on digg , , 2008 .
kempe d. and kleinberg j. and tardos e. maximizing the spread of influence through a social network , , 2003 .
domingos p. and richardson m. mining the network value of customers , , 2001 .
young p.h .
the diffusion of innovations in social networks , , 2003 .
reed w.j . and jorgensen m. the double pareto - lognormal distribution - a new parametric model for size distributions , 2003 .
lerman k. and ghosh r. and kang j.h .
centrality metric for dynamic networks , , 2010 .


recent advances in atom de broglie waveguide technology @xcite and its potential applicability to atom interferometry @xcite and integrated atom optics @xcite create a need for accurate theoretical modelling of such systems in the low temperature , tight waveguide regime where transverse excitations are frozen out and the quantum dynamics becomes essentially one - dimensional ( 1d ) ( tonks - gas limit ) .
it has been shown by olshanii @xcite , and also recently by petrov et al .
@xcite , that at sufficiently low temperatures and densities , high transverse frequencies @xmath0 , and large positive scattering length , where thermal and longitudinal zero - point energies are small compared with @xmath1 , a bose - einstein condensate ( bec ) in a thin cigar - shaped trap has dynamics which reduce to those of a 1d gas of hard core , or impenetrable , point bosons .
this is a model for which the exact many - body energy eigensolutions were found in 1960 using an exact mapping from the hilbert space of energy eigenstates of an _ ideal _ gas of spinless fermions to that of many - body eigenstates of hard core , and therefore _ strongly interacting _ , bosons @xcite .
in this limit there are strong short - range pair correlations which are omitted in the gross - pitaevskii ( gp ) approximation , which assumes that all @xmath2 bosons occupy the same orbital ( complete bec , condensed fraction unity ) . in the absence of a trap potential
it is known @xcite that the occupation of the lowest orbital is of order @xmath3 where @xmath2 is the total number of atoms , in contrast to @xmath2 for the ideal bose gas as well as the gp approximation .
nevertheless , this system exhibits some bec - like behavior such as talbot recurrences following an optical lattice pulse @xcite and dark soliton - like behavior in response to a phase - imprinting pulse @xcite .
the case of harmonically trapped , hard core bosons in 1d is more relevant to recent atom waveguide experiments @xcite .
the spatial profile of the single - particle density is expressible in closed form , and has recently been shown @xcite to be well approximated by a modified 1d effective field theory , although we have recently shown in a numerically accurate time - dependent calculation @xcite that spatial interference fringes of separated and recombined condensates in the exact many - body solution are much weaker than those predicted by the corresponding time - dependent mean field theory @xcite .
although the fermi - bose mapping theorem @xcite implies that all physical properties expressible in terms of spatial configurational probabilities are the same for the actual bosonic system and the fictitious
spinless fermion " system used for the mapping , the momentum distribution of the bosonic system , or more generally its occupation distribution over the relevant orbitals for a given geometry , is very different in the bosonic system .
it is known @xcite that for a spatially uniform system of hard core bosons in 1d , the momentum distribution is strongly peaked in the neighborhood of zero momentum , whereas that of the corresponding fermi system is merely a filled fermi sea . in the case of hard core bosons in a 1d harmonic trap , it has been an interesting question whether the system undergoes true bec or merely an attenuated one such as that in the uniform system .
ketterle and van druten @xcite have shown that true bec occurs for a finite number of atoms in a 1d harmonic oscillator ( ho ) in the case of an _ ideal _ gas , but the behavior turns out to be different for the exact many - body solutions in the tonks ( 1d ) limit . the most fundamental definition of bec and the condensate orbital is based on the large distance behavior of the one - particle reduced density matrix @xmath4 . if off - diagonal long - range order ( odlro ) is present and hence the largest eigenvalue of @xmath5 is macroscopic ( proportional to @xmath2 ) then the system is said to exhibit true bec and the corresponding eigenfunction , the condensate orbital , plays the role of an order parameter @xcite .
although the precise definition of odlro requires a thermodynamic limit not strictly applicable to mesoscopic traps , the gp approximation assumes from the start that odlro and macroscopic occupation of a single orbital are good approximations in a trap , so examination of this assumption is important .
we have recently @xcite used the fermi - bose mapping theorem to determine , for this case of 1d impenetrable bosons in a harmonic trap , the exact many - body ground state and its salient features , including the one - particle reduced density matrix and its eigenvalues ( occupation number distribution function ) and eigenfunctions ( natural orbitals ) , as well as the momentum distribution function .
we find that the largest eigenvalue of the one - particle density matrix is proportional to lower than the first power of @xmath2 , as is the momentum distribution function at @xmath6 , so that the tonks gas does not show true bec .
in this paper we give an overview of our recent work on exact many - body quantum solutions for 1d tonks gases .
in particular , secs .
ii and iii describe the ground state and dynamic properties of tonks gases , respectively . for the most part we investigate 1d tonks gases assuming the conditions are satisfied for realizing them . in sec .
iv we examine the use of bessel optical dipole traps for experimental realization of 1d tonks gases .
we first describe the fermi - bose mapping method for obtaining exact solutions of the time - independent many - body schrdinger equation in the tonks limit ( impenetrable point bosons ) and its application to the spatially uniform case ( no trap potential ) .
this description is very brief and the original literature @xcite should be consulted for details .
some salient features of the recent generalization to the harmonically trapped tonks gas @xcite are then described .
our basic model consists of @xmath2 bosonic atoms at zero temperature moving in 1d ( @xmath7 ) , the quantum motion in the other two spatial dimensions having been frozen out by tight transverse confinement via an atomic waveguide .
the schrdinger hamiltonian is then assumed to have the structure : @xmath8 where @xmath9 is the one - dimensional position of the @xmath10 particle ( atom ) and @xmath11 is symmetric ( invariant ) under permutations of the particles .
the two - particle interaction potential is assumed to contain a hard core of 1d diameter @xmath12 .
this is conveniently treated as a constraint on allowed wave functions @xmath13 : @xmath14 rather than as an infinite contribution to @xmath11 , which then consists of all other ( finite ) interactions and external potentials .
let @xmath15 be a fermionic solution of @xmath16 which is antisymmetric under all particle pair exchanges @xmath17 , hence all permutations .
one can consider @xmath18 to be either the wave function of a fictitious system of `` spinless fermions '' , or else that of a system of real fermions whose spins are all aligned , as in magnetically trapped atomic vapor becs .
define a  unit antisymmetric function " @xcite @xmath19 where @xmath20 is the algebraic sign of the coordinate difference @xmath21 , i.e. , it is + 1(-1 ) if @xmath22(@xmath23 ) . for given antisymmetric @xmath18 ,
define a bosonic wave function @xmath24 by @xmath25 which defines the fermi - bose mapping .
@xmath24 satisfies the hard core constraint ( 2 ) if @xmath18 does , is totally symmetric ( bosonic ) under permutations , obeys the same boundary conditions as @xmath18 , and @xmath26 follows from @xmath27 @xcite . in the case of periodic boundary conditions ( no trap potential , spatially uniform system )
one must add the proviso that the boundary conditions are only preserved under the mapping if @xmath2 is odd , but the case of even @xmath2 is easily accomodated by imposing periodic boundary conditions on @xmath18 but _ anti_periodic boundary conditions on @xmath24 . for a trapped system the boundary condition that wave functions vanish at infinity
is preserved for all @xmath2 .
the mapping theorem leads to explicit expressions for all many - body energy eigenstates and eigenvalues of a 1d scalar condensate ( bosons all of the same spin ) under the assumption that the only two - particle interaction is a zero - range hard core repulsion , represented by the @xmath28 limit of the hard - core constraint .
such solutions were obtained in sec .
3 of the original work @xcite for periodic boundary conditions and no external potential .
such a situation could be realized experimentally , for example , by trapping the tonks gas around the peak of an tight toroidal optical dipole trap formed using a red - detuned laguerre - gaussian laser beam @xcite . in a system of 1d bosons with a hard - sphere interaction of diameter @xmath12 ,
it is sufficient at low densities @xcite to consider the case of impenetrable point particles , the zero - range limit @xmath28 .
since wave functions of  spinless fermions " are antisymmetric under coordinate exchanges , their wave functions vanish automatically whenever any @xmath29 , the constraint has no effect , and the corresponding fermionic ground state is the ground state of the _ ideal _ gas of fermions , a slater determinant of the lowest @xmath2 single - particle plane - wave orbitals .
the exact many body ground state was found to be a pair product of bijl - jastrow form @xmath30| , \ ] ] where @xmath31 is the perimeter of the annular trap . in spite of the very long range of the individual pair correlation factors @xmath32|$ ] , the pair distribution function @xmath33 , which is physically the joint probability density that if one particle is found at @xmath34 a second will be found at @xmath9 , was found to be of short range @xmath35 with @xmath36 the linear number density of the system .
clearly , @xmath37 which reflects the hard core nature of the two - particle interaction . by examination of the excited states the system
was also found to support propagation of sound with speed @xmath38 @xcite .
here we briefly review our recently - obtained exact solution for the many - body ground state of a harmonically trapped tonks gas @xcite .
the hamiltonian of n bosons in a 1d harmonic trap is @xmath39   .\ ] ] we again assume that the two - body interaction potential consists only of a hard core of negligible diameter @xmath40 .
it follows from the mapping theorem that the exact n - boson ground state @xmath41 is @xmath42 where @xmath43 is the ground state of a fictitious system of @xmath2 spinless fermions with the same hamiltonian and constraint .
since the hard - core interaction has no effect in the zero - range limit for fermions , whose wave functions already vanish at contact due to antisymmetry , it follows that the fermionic ground state is a slater determinant of the lowest @xmath2 single - particle eigenfunctions @xmath44 of the harmonic oscillator ( ho ) : @xmath45 the ho orbitals are @xmath46 with @xmath47 the hermite polynomials and @xmath48 being the ground state width of the harmonic trap for a single atom . by factoring the gaussians out of the determinant and carrying out elementary row and column operations
, one can cancel all terms in each @xmath49 except the one of highest degree , with the result @xcite @xmath50 substitution into ( 6 ) then yields a simple but exact analytical expression of bijl - jastrow pair product form for the @xmath2-boson ground state : @xmath51 \prod_{1\le j
< k\le n}|x_{k}-x_{j}|\ ] ] with @xmath52 and normalization constant @xmath53^{-1/2 }   .\ ] ] it is interesting to note the strong similarity between this exact 1d @xmath2-boson wave function and the famous laughlin variational wave function of the 2d ground state for the quantized fractional hall effect @xcite , as well as the closely - related wave functions for bosons with weak repulsive delta - function interactions in a harmonic trap in 2d found recently by smith and wilkin @xcite .
both the single particle density and pair distribution function depend only on the absolute square of the many - body wave function , and since @xmath54 they reduce to standard ideal fermi gas expressions .
the single particle density , normalized to @xmath2 , is @xmath55 we shall not exhibit it here since it has recently been calculated by kolomeisky _
_ @xcite ; see also our recent discussion of the time - dependent case @xcite . the pair distribution function , normalized to @xmath56 , is @xmath57 noting that terms with @xmath58 , which vanish by antisymmetry , can be formally added to the summation ( 13 ) , one can write @xmath59 although the hermite polynomials have disappeared from the expression ( 10 ) for the many - body wave function , they reappear upon integrating @xmath60 over @xmath61 coordinates to get the single particle density @xmath62 and over @xmath63 to get the pair distribution function @xmath64 , and the expressions in terms of the ho orbitals @xmath65 are the most convenient for evaluation .
figure [ fig : one ] shows a gray scale plot of the dimensionless pair distribution function @xmath66 versus the normalized coordinates @xmath67 for a ) @xmath68 , b ) @xmath69 , and c ) @xmath70 .
some qualitative features of the pair distribution function are apparent : in the first place it follows either from the original expression ( 12 ) or from eq .
( 14 ) that @xmath71 vanishes at contact @xmath72 , as it must because of impenetrability of the particles , and we see this to be true in fig . [
fig : one ] .
furthermore , the correlation term @xmath73 is a truncated closure sum and approaches the dirac delta function @xmath74 as @xmath75 , as is to be expected since the healing length in a spatially uniform 1d hard core bose gas varies inversely with particle number @xcite . as a result
the width of the null around the diagonal @xmath76 decreases with increasing @xmath2 , and vanishes in the limit .
away from the diagonal along @xmath77 the pair distribution function rises , exhibits modulations for @xmath78 , due to the oscillatory nature of the ho orbitals , before decreasing back to zero at large distances . for @xmath79
much larger than the healing length , @xmath80 reduces to the uncorrelated density product @xmath81 , so the spatial extent of the pair distribution function is that of the density and varies as @xmath82 @xcite .
the reduced single - particle density matrix with normalization @xmath83 is @xmath84 although this multi - dimensional integral can not be evaluated analytically , it can be evaluated numerically by monte carlo integration for not too large values of @xmath2 ( the computing time scales as @xmath85 ) .
figure [ fig : two ] shows a gray scale plot of the dimensionless reduced single - particle density matrix @xmath86 versus the normalized coordinates @xmath87 and @xmath88 for a ) @xmath68 , b ) @xmath69 , and c ) @xmath70 . along the diagonal @xmath89 reproduces the single - particle density @xcite .
the off - diagonal elements relate to odlro , and it is clear that as @xmath2 increases the off - diagonal elements are decreasing , in contrast with the diagonal .
this is an indication that odlro vanishes for a system of hard core bosons in a 1d ho in the thermodynamic limit , although at present only numerical evidence exists , there being no analytical proof generalizing the result of lenard for the untrapped tonks gas @xcite .    in a macroscopic system ,
the presence or absence of bec is determined by the behavior of @xmath4 as @xmath90 .
off - diagonal long - range order is present if the largest eigenvalue of @xmath5 is macroscopic ( proportional to @xmath2 ) , in which case the system exhibits bec and the corresponding eigenfunction , the condensate orbital , plays the role of an order parameter @xcite .
although this criterion is not strictly applicable to mesoscopic systems , if the largest eigenvalue of @xmath5 is much larger than one then it is reasonable to expect that the system will exhibit some bec - like coherence effects .
thus we examine here the spectrum of eigenvalues @xmath91 and associated eigenfunctions @xmath92 ( `` natural orbitals '' ) of @xmath5 .
although natural orbitals are a much - used tool in theoretical chemistry , they have only recently been applied to mesoscopic atomic condensates @xcite .
the relevant eigensystem equation is @xmath93 @xmath91 represents the occupation of the orbital @xmath94 , and one has @xmath95 .
numerical evaluation of the integral by discretization yields a matrix eigensystem equation giving accurate numerical results for the largest eigenvalues and associated eigenvectors .
we remark that for the corresponding problem of @xmath2 free fermions in a 1d ho , in which case @xmath96 is evaluated using the fermion wave function @xmath97 , the natural orbitals are simply the ho orbitals , and @xmath98 , all higher eigenvalues being zero ( filled fermi sea )
. however , the n - boson wave function is the modulus of the fermion wave function , and this leads to significant differences in the spectrum of natural orbitals and eigenvalues for the hard core bose gas . in fig .
[ fig : three](a ) we show a log - log plot of the fractional occupation of the lowest orbital @xmath99 versus the total particle number @xmath2 ( solid line ) , along with a best fit power - law @xmath100 ( dashed line ) .
this is to be contrasted with the case of a spatially uniform system of hard core bosons for which @xmath101 @xcite . in both cases
the fractional occupation decreases with increasing @xmath2 , and thus does not correspond to a true condensate for which @xmath102 .
nevertheless , the occupation of the lowest orbital may still be large , @xmath103 , and is larger than the spatially uniform case @xmath104 , so macroscopic quantum coherence effects reminiscent of bec can still result @xcite .
figure [ fig : three](b ) shows the distribution of occupations @xmath91 versus orbital number @xmath105 ( the orbitals are ordered according to eigenvalue magnitude , the largest eigenvalue being @xmath106 ) for @xmath68 ( circles ) , @xmath69 ( stars ) , and @xmath70 ( squares ) .
this figure shows that as the lowest orbital occupation @xmath107 increases with increasing @xmath2 so does the range of significantly occupied higher - order orbitals with @xmath108 .
this means that the dominance of the lowest orbital decreases with increasing @xmath2 , so singling out @xmath109 as a macroscopic wave function for the whole system becomes more problematic with increasing @xmath2 @xcite .
next we examine the momentum distribution for the trapped 1d tonks gas . for a spatially uniform system ( no trap )
the natural orbitals are plane waves , so the occupation distribution of the natural orbitals is the same as the momentum distribution .
although this is not the case here due to the effect of the harmonic trap potential , the momentum distribution is still physically important . in terms of the boson annihilation and creation operators in position representation (
quantized bose field operators ) the one - particle reduced density matrix is @xmath110 the momentum distribution function @xmath111 , normalized to @xmath112 , is @xmath113 where @xmath114 is the annihilation operator for a boson with momentum @xmath115 . then @xmath116 the spectral representation of the density matrix then leads to @xmath117 where the @xmath118 are fourier transforms of the natural orbitals : @xmath119 .
figure [ fig : four ] shows the numerically calculated dimensionless momentum spectrum @xmath120 versus normalized momentum @xmath121 , with @xmath122 , for a ) @xmath68 , b ) @xmath69 , and c ) @xmath70 .
the key features are that the momentum spectrum maintains the sharp peaked structure reminiscent of the spatially uniform case @xcite for the 1d ho , and that the peak becomes sharper with increasing atom number @xmath2 .
this is to be expected since as the number of atoms increase the many - body repulsion causes the system to become more spatially uniform within the trap interior .
by way of contrast , for a 1d fermi gas the corresponding momentum spectrum is a filled fermi - sea , here rounded out by the trapping potential .
in particular , the momentum distribution for the 1d trapped fermi gas can be expressed as @xmath123 thus , the momentum spectrum provides a means of distinguishing between the 1d fermi and tonks gases .
figure [ fig : five ] shows an example of the tonks ( dashed line ) and fermi gas ( solid line ) momentum spectra for @xmath70 atoms in a harmonic trap , and the difference is clearly seen .
furthermore , the distinction only grows larger with increasing number of atoms , as the filled fermi - sea broadens whereas the tonks spectrum becomes narrower . in a recent paper
we have devised a scheme to measure the momentum spectrum of trapped gases based on raman outcoupling , and we refer the reader to our paper for details @xcite .
the fermi - bose mapping theorem is very easily generalized so as to yield exact solutions of the time - dependent many - body schrdinger equation ( tdmbse ) by noting that since the mapping function @xmath124 of eq .
( 3 ) is independent of time , one can merely replace @xmath125 by @xmath126 , implying that the relationship of eq .
( 4 ) , with time arguments added to the bose and fermi wave functions , is valid for solutions of the tdmbse . in the special case where the only interatomic interaction is that of hard cores of vanishing diameter , but external potentials including
a trap potential as well as time - dependent fields may be present , the many - fermion solutions of the tdmbse are slater determinants of solutions of the _ single particle _
tdse in the given external potential , and each many - boson solution is obtained by multiplying the corresponding determinant by the mapping function @xmath124 .
some salient features of our recently - obtained dynamical solutions for solitons in a ring geometry and for interference effects in a dynamically split , harmonically trapped tonks gas will be described .
the description will again be brief since details are available in the literature @xcite .
dark and gray solitons are a generic feature of the nonlinear schrdinger equation with repulsive interactions , and several calculations of their dynamics based on the mean - field gross - pitaevskii ( gp ) equation have appeared @xcite , as well as experiments demonstrating their existence in atomic becs @xcite . since the underlying _
many - body _ schrdinger equation is linear , this raises the question of how observed solitonic behavior arises . here
this issue will be examined with the aid of exact many - body solutions for the tonks gas .
the model consists of a 1d hard - core bose gas in a toroidal trap , or ring , with cross section so small that motion is essentially circumferencial .
the fermi - bose mapping is employed to generate exact solutions for this problem .
we identify stationary solutions which reflect some properties of dark solitons from the gp theory when the ring is pierced at a point by an intense blue - detuned laser .
we also present dynamical solutions when half of an initially homogeneous ring bec is phase - imprinted via the light - shift potential of an applied laser , leading to gray soliton - like structures whose velocity depends on the imposed phase - shift @xcite .
such structures are apparent for times less than the echo time @xmath127 , with @xmath31 the ring circumference and @xmath128 the speed of sound in the bec . on
longer time scales the dynamics becomes very complex showing talbot recurrences which are beyond the gp theory .
consider @xmath2 bosons in a tight toroidal trap , for example a toroidal optical dipole potential @xcite , and denote their 1d positions measured around the circumference by @xmath9 .
this is equivalent to the exactly - solved model @xcite of @xmath2 impenetrable point bosons in 1d with wave functions satisfying periodic boundary conditions with period @xmath31 equal to the torus circumference , and the fundamental periodicity cell may be chosen as @xmath129 .
however , the rotationally invariant quantum states of this problem do not reveal any dark soliton - like structures .
to proceed we therefore consider the case that a blue - detuned laser field pierces the ring at @xmath130 by virtue of the associated repulsive dipole force : the light sheet then provides a reference position for the null of the dark soliton .
assume that the light sheet is so intense and narrow that it may be replaced by a constraint that the many - body wave function ( hence the orbitals @xmath131 ) must vanish whenever any @xmath132 .
then the appropriate orbitals @xmath133 are free - particle energy eigenstates vanishing at @xmath130 and periodic with period @xmath31 .
the complete orthonormal set of even - parity eigenstates @xmath134 and odd - parity eigenstates @xmath135 are @xmath136   , \nonumber\\ \phi_{n}^{(-)}(x ) & = & \sqrt{2/l}\sin(2n\pi x / l )        , \end{aligned}\ ] ] with @xmath137 running from @xmath138 to @xmath139 .
the odd eigenstates are the same as those of free particles with no @xmath130 constraint , since these already vanish at @xmath130 . however , the even ones are strongly affected by the constraint , their cusp at @xmath130 being a result of the impenetrable light sheet at that point . if one bends a 1d box @xmath140 with impenetrable walls into a ring , identifying the walls at @xmath141 , then those particle - in - a - box eigenfunctions which are even about the box center become identical with the @xmath134 , and their cusp results from the nonzero slope of these functions at the walls . the @xmath2-fermion ground state is obtained by inserting the lowest @xmath2 of these orbitals into a slater determinant ( filled fermi sea ) . since @xmath142 , the one - particle density @xmath62 of the corresponding many - boson ground state given by eq .
( 4 ) is the same as that of the @xmath2-fermion ground state . in the thermodynamic limit @xmath75 , @xmath143 , @xmath144 for fixed @xmath7 , @xmath145 ,
one finds @xcite @xmath146   .\ ] ] @xmath62 vanishes at @xmath130 and approaches the mean density @xmath147 over a healing length @xmath148 with damped spatial oscillations about its limiting value .
suppose next that the light - sheet is turned off at @xmath149 by removing the constraint that the wave function vanish at @xmath130 .
the solution of the tdmbsb for the many - boson system is then still given by the mapping theorem , but the slater determinant representing the corresponding many - fermion state has to be built from solutions of the time - dependent single - particle schrdinger equation satisfying the initial conditions of eq .
( [ lightsheet ] ) .
the odd - parity solutions are stationary in time since they already vanish at the position of the light sheet .
the even - parity solutions are nontrivial , but are expressible as sums over the space - periodic , even - parity solutions of the time - independent free - particle schrdinger equation @xcite : @xmath150 where @xmath151 and @xmath152 .
the time - dependent density @xmath153 is then the sum of absolute squares of all @xmath2 orbitals in the fermi sea .
it is found @xcite that there are two important time scales : one is the poincar recurrence time @xmath154 . noting that @xmath155 is proportional to @xmath156 , one finds that all terms in the sum are time - periodic with period @xmath157 , which is therefore the recurrence time for the density and in fact all properties of our model @xcite .
the other important time is the echo time @xmath158 , the time for sound to make one circuit around the torus .
recalling that the speed of sound in this system is @xmath38 @xcite , one finds @xmath159 .
for @xmath160 after the constraint is removed , the initial density develops sound waves that propagate around the ring , and that we examine below in the context of phase - imprinting . for @xmath161
the evolution is very complex , but complete recurrences occur for times @xmath162 with fractional revivals in between .
suppose next that the toroidal impenetrable point bose gas is in its ground state for times @xmath163 ( no light sheet obstacle in this case ) , and a phase - imprinting laser pulse is applied as a delta - function pulse over half the ring at @xmath149 .
this is described by the hamiltonian @xmath164\ ] ] where @xmath165 , i.e. , unity for @xmath166 and zero elsewhere .
this is the technique used in recent experiments @xcite , here idealized to a delta - function in time and to sharp spatial edges . before the pulse
the most convenient free - particle orbitals in ( 5 ) are plane waves @xmath167 where @xmath168 and @xmath169 with @xmath170 .
let @xmath171 be the solution of the corresponding single - particle tdse reducing to the above @xmath172 just before the pulse .
then the solutions just after the pulse are @xmath173 .
the potential gradients at the pulse edges impart momentum kicks to the particles there which induce both compressional waves propagating at the speed , @xmath128 , of sound and density dips ( gray solitons ) moving at speeds @xmath174 .
the expansion of @xmath171 in terms of the unperturbed plane waves is evaluated as @xmath175 and the time - dependent density is the sum of the absolute squares of the lowest @xmath2 of these .
figure [ fig : six ] shows numerical simulations for @xmath176 , @xmath177 , and @xmath178 ( solid line ) , and @xmath179 ( dashed line ) : due to symmetry we show only half of the ring @xmath180 , the phase shift being imposed at @xmath181 .
considering times short compared to the echo time means that the corresponding results are not very sensitive to the periodic boundary conditions , and also therefore apply to a linear geometry .
the initial density profile is flat with a value @xmath182 . for both phase shifts two
distinct maxima are seen , which travel at close to the speed of sound @xmath128 , and two distinct minima , which are analogous to gray solitons and travel at velocities @xmath183 .
in addition , there are also high wavevector oscillations which radiate at velocities greater than @xmath128 , analogous to precursors in electromagnetic wave propagation in a medium . in the case of a phase shift @xmath178 ,
the density is symmetric about @xmath181 , whereas for a phase - shift other than a multiple of @xmath184 the evolution is not symmetric , see the dashed line where the global minimum moves to the right in reponse to the phase shift . in fig .
[ fig : seven ] we plot the calculated velocity of the global density minimum relative to the speed of sound for a variety of phase shifts @xmath185 .
the basic trend is that larger phase - shift means lower velocity , in qualitative agreement with recent experiments @xcite , but there is a sharp velocity peak at @xmath186 : this peak results from the cross - over between two local minima in the density .
these general features , the generation of gray solitons and density waves , agree with those of the gp theory , but here arise out of the exact many - body calculation .    a detailed comparison between our results and current experiments is not possible as they do not conform to the conditions for a 1d system .
however , some estimates are in order to set the appropriate time scales : if we consider @xmath187rb with a ring of circumference @xmath188 @xmath189 m , and a high transverse trapping frequency @xmath190 hz , then we are limited to atom numbers @xmath191 @xcite , so these are small condensates .
we then obtain @xmath192 s , and @xmath193 ms for @xmath176 .
finally , we remark that since our approach relied on the mapping between the strongly - interacting bose system and a non - interacting  spinless fermi gas " model , this suggests that dark and gray solitons should also manifest themselves in the density for the 1d fermi system .
although real fermions have spin , the interactions used here to generate solitons were spin - independent .      mean - field theory ( mft ) has proven remarkably successful at predicting both the static and dynamic behavior of bose - einstein condensates ( becs ) in weakly - interacting atomic vapors @xcite , including the ground state properties @xcite , the spectrum of collective excitations @xcite , four - wave mixing @xcite , matter - wave solitons @xcite , and interference between becs @xcite .
the basic notion underlying mft is of a macroscopic wave function @xcite , or order parameter , which defines the spatial mode into which a significant fraction of the atoms condense below the critical temperature .
the macroscopic wave function typically obeys a nonlinear schrdinger equation ( nlse ) , the gross - pitaevskii equation , and is most suitable for dilute bose gases .
however , the success of mft is not assured in all cases .
for example , in one - dimension ( 1d ) a spatially homogeneous ideal gas in its many - body ground state exhibits complete bec into the lowest single - particle state , but no bec at any nonzero temperature .
furthermore , previous exact analysis @xcite of a spatially uniform tonks gas by one of us ( mg ) and its extension by lenard @xcite and vaidya and tracey @xcite have shown that in the many - body ground state the occupation of the lowest single - particle state is of order @xmath3 where @xmath2 is the total number of atoms , in contrast to @xmath2 for usual bec , and similar behavior , but with a slightly larger exponent , occurs in the harmonically trapped tonks gas @xcite . we have recently shown that mean - field theory breaks down in the analysis of the dynamics of 1d atom clouds , in that it predicts interference effects that are absent in the exact theory @xcite .
_ @xcite have proposed a nonlinear schrdinger equation ( nlse ) with a quartic nonlinearity to extend the usual mean - field theory for 1d atom clouds . for a harmonic trap the ground - state density profiles from their theory show excellent agreement with the exact many - body results ( see fig . 1 of their paper ) .
the key question , then , is whether this extended nlse can be used in all circumstances . to address this issue we examine the problem of a 1d atomic cloud in the ground state of a harmonic trap that is split by a blue - detuned laser and recombined , both using an exact many - body treatment based on the fermi - bose mapping and the approximate nlse :
the nlse predicts interference whereas the exact analysis does not .
consider a tonks gas which is initially in its @xmath2-body ground state in a harmonic trap @xcite .
a central gaussian repulsive potential simulating a blue - detuned laser field is turned on quasi - adiabatically at time @xmath149 to split the initial state . after some time @xmath194 both the harmonic trap and repulsive potential are turned off and the two split components allowed to recombine : this is an interference experiment of the cool , cut , interfere variety previously discussed @xcite .
the external potential is taken of the specific form @xmath195 and @xmath196 for @xmath197 .
figure [ fig : eight ] shows an illustrative example for @xmath70 with a repulsive potential of height @xmath198 , and width @xmath199 , with @xmath200 the ground state harmonic oscillator width , and @xmath201 .
this figure shows a gray scale plot of the single - particle density @xmath202 of the @xmath2-boson system as functions of @xmath203 ( horizonatal axis ) and position @xmath204 ( vertical axis ) , with white being the highest density .
the potential height was chosen such that @xmath205 , where @xmath189 is the chemical potential of the @xmath2-particle oscillator ground state @xcite , noting that the top of the @xmath2-particle fermi sea is at @xmath206 .
as expected , as the repulsive potential turns on it splits the initial ground state into two separated components . upon release at @xmath207 the two components
expand and subsequently recombine .
what is noteworthy is that although there is some modulation upon recombination there are no strong interference fringes indicative of the interference provisionally expected for bosons : this was a generic finding from our simulations irrespective of the time scale on which the repulsive potential was turned on @xcite .
in contrast , the density @xmath208 for each individual orbital @xmath209 can show strong fringes , but with a different period in each case .
thus , the minimal inteference is a result of washing out of the individual interferences by summing over n orbitals .
thus , the remnant of any interference fringes decreases with increasing @xmath2 and vanishes in the thermodynamic limit .
physically , it makes sense that the interference fringes are all but absent since the fermi - bose mapping shows that in this 1d limit the system of bosons acts effectively like a system of free fermions insofar as effects expressible only in terms of @xmath210 are concerned , so interference is not expected @xcite .
the lack of interference is therefore a signature of the fermi - bose duality that occurs in 1d systems of impenetrable particles @xcite .
we next turn to the mean - field description proposed by kolomeisky _
@xcite for low - dimensional systems . in particular
, they introduce an order parameter @xmath211 , normalized to the number of particles @xmath2 , for such systems . using energy functional arguments
they deduce the following nlse with quartic nonlinearity for a 1d system of impenetrable bosons : @xmath212\phi   .\ ] ] our goal is to compare the predictions of this nlse for the same cool , cut , interfere simulation in figure [ fig : eight ] , with the initial condition @xmath213 corresponding to the exact many - body solution , all other parameters being equal .
figure [ fig : nine ] shows the corresponding gray - scale plot of the density @xmath214 , and two features are apparent : first , during the splitting phase when the repulsive potential is on there is very good overall agreement between the exact many - body theory and the nlse prediction .
second , when the split components are released and recombine they produce pronounced interference fringes , in contrast to the exact theory .
indeed , this interference in the mft is to be expected on the basis of previous theoretical work @xcite , even though a quartic ( rather than quadratic ) nonlinearity is employed here .
thus the mft can not accurately capture the time - dependent dynamics in all situations .
so far we have assumed that conditions are satisfied for the tonks limit . however , previous analysis by olshanii @xcite and petrov et al .
@xcite show that there are stringent requirements on the temperature , transverse confinement , and linear atom number density to approach the tonks limit .
clearly , what is needed first of all it to create a thin cigar shaped atomic trap so that the transverse mode becomes frozen and the atomic motion becomes essentially 1d .
recent experimental developments suggest that tonks gases may be realizable in magnetic atom waveguides @xcite , and bongs et al .
@xcite have proposed a hybrid trap composed of optical dipole trap formed with a first - order lg beam combined with magnetic longitudinal trapping .
here we examine the use of bessel optical dipole traps as an all - optical means to realize 1d tonks gases .
ideal bessel beams are solutions of the free - space wave equation which propagate with unchanging beam profile along the propagation axis which we take as @xmath7 in cylindrical coordinates @xmath215 . the electric field of a monochromatic , linearly polarized ideal bessel beam of order @xmath216 and frequency @xmath217 is @xcite @xmath218 where @xmath219 is a scale electric field value , @xmath220 is the @xmath216th order bessel function
, @xmath221 is the azimuthal mode number which we take as positive for simplicity in notation , and @xmath222 and @xmath223 are the radial and longitudinal wavevectors such that @xmath224 with @xmath225 .
the zeroth - order solution @xmath226 has a central maximum surrounded by concentric rings of roughly equal power while the higher - order solutions @xmath220 have zero on - axis intensity also with concentric rings .
the ultility of bessel beams for optical dipole traps lies in the fact that they can produce very elongated , and hence low aspect ratio quasi-1d traps . in practice
, of course , one can not produce an ideal bessel optical beams as it carries infinite power .
however , using axiconal optics finite power approximations to bessel beams can be produced which can propagate over significant distances .
the basic scheme is illustrated in figure [ fig : ten ] where an incident laguerre - gaussian beam of azimuthal mode index @xmath216 is transformed into a bessel beam of the same order .
the interested reader is referred to our recent paper @xcite for details , but for an incident gaussian beam on the axiconal optics the resulting zero - order bessel beam intensity profile past the axicon is @xmath227 where @xmath228 is the input spot size , @xmath229 the input power , @xmath230 , @xmath137 and @xmath231 are the index and opening angle of the axicon , and @xmath232 which gives the longitudinal extent of the bessel beam .
this intensity profile of the bessel optical beam is readily converted to an effective dipole potential @xmath233 with @xmath234 the laser detuning from the optical transition frequency @xmath235 , @xmath236 the natural linewidth of the optical transition , @xmath237 is the resonant saturation intensity , and @xmath238 . for a red - detuned laser
the potential is negative and the atoms are attracted to the regions of high intensity , whereas for a blue - detuned laser the atoms are repelled into the low field regions . here for a red - detuned laser we have approximated the bessel optical dipole potential by a parabolic potential near the axis under conditions of tight confinement , where @xmath239 is the longitudinal position of the peak of the bessel beam , and @xmath240 a red - detuned ( @xmath241 ) @xmath226 optical dipole potential therefore provides confinement in both the radial and longitudinal directions . here @xmath242 is the radial oscillation frequency with corresponding ground state oscillator width @xmath243 , and @xmath244 is the ratio between the longitudinal and radial trap frequencies @xmath245 @xcite , which also determines the aspect ratio between the radial and longitudinal ground state widths @xmath246 ( in the absence of many - body repulsion ) .
petrov _ et al .
_ @xcite have theoretically studied the diagram of state for a one - dimensional gas of trapped bosons , assuming @xmath247 , and found that a true bec , or at least a quasi - condensate , with concomitant macroscopic occupation of a single state , is only attained for high enough particle numbers @xmath248 with @xmath249 for @xmath250 and temperatures @xmath251 , one obtains a tonks gas of impenetrable bosons for which hard core repulsion between the bosonic atoms prevents them from penetrating through each other in the one - dimensional system , and the system acquires properties reminiscent of a one - dimensional system of fermionic atoms .
a highly elongated bessel beam discussed , say with parameters @xmath252 nm , @xmath253 w , @xmath254 cm , @xmath255 @xmath189 m , would be an ideal candidate for the experimental realization of a tonks gas . the low aspect ratio @xmath256 and tight radial confinement @xmath257 nm result in a high upper boundary @xmath258 for the particle number of the tonks gas .
for the commonly used @xmath187rb isotope with a scattering length @xmath259 nm one finds @xmath260 . although this is still a fairly low value it should be possible to experimentally realize a small @xmath187rb tonks gas .
however , more promising would be the use of the @xmath261rb isotope , where a feshbach resonance can be used to tune the normally negative scattering length to positive values of several hundred nanometers magnitude @xcite . as @xmath258 is proportional to the square of the scattering length even a moderate increase to @xmath262 nm
would make it possible to create a larger tonks gas , with say @xmath263 atoms , which should be easily detectable .
the bessel beam trap offers some advantages compared to alternative suggested approaches using magnetic waveguides @xcite and a hybrid magnetic - optical trap @xcite .
firstly , it involves only a very simple _ all - optical _ system for which the aspect ratio of the trap may be controlled simply by varying the gaussian spot size incident on the axicon .
more specifically , being all - optical , it does not involve material surfaces as in magnetic waveguides , which can cause matter - wave decoherence @xcite .
furthermore , it allows for the possibility of trapping multiple magnetic sublevels and the investigation of multi - component tonks gases , which would not be possible in the hybrid magnetic - optical trap of bongs _ et al .
in this paper we have given an overview of our work on the quantum dynamics of 1d trapped tonks gases . in particular
, we hope to have shown that tonks gases display a rich variety of behaviors that are worthy of experimental investigation .
the fermi - bose mapping approach to solving for the 1d tonks shows that many features are shared between the tonks gas of impenetrable bosons and a 1d gas of non - interacting fermions : they have the same ground state energies and density profiles , and neither shows true bec . however , the tonks and fermi gases differ markedly in their momentum spectra , and we have devised a scheme to measure this difference @xcite . furthermore , we have shown that the 1d tonks gas can show dark soliton - like solutions which are typically associated with solutions of the mean - field gross - pitaevskii equation .
this is so in spite of the fact that mean - field theory greatly overestimates the coherence properties of the tonks gas , as we demonstrated in our simulations of split and recombined tonks gases .
our work presented here shows that tonks gases have some degree of coherence but much less than a true bec , and this may have implications for proposed atomic interferometers employing high density , eg .
atom laser , sources .
in particular , strong transverse confinement is desirable in atomic interferometers to avoid multi - transverse mode effects , and at the same time one would like to keep the atomic density down to avoid many - body shifts , but these are exactly the requirements for the tonks limit !
therefore it is very important to study the bec - tonks transition @xcite and its effect on the performance of atom interferometers .
we hope that our work is a first step in this direction .
this work was supported at the university of arizona by office of naval research grant n00014 - 99 - 1 - 0806 and also the us army research office .
we thank joe triscari for collaborations on the ground state properties of the tonks gas , and the work on bessel optical dipole traps was conducted in collaboration with j. arlt and k. dholakia from st .
andrews university , scotland , and j. soneson at the university of arizona .
et al . _ , .
d. mller _ et al .
n.h . dekker _
et al . _ , .
thywissen , r.m .
westervelt , and m. prentiss , .
hinds , m.g .
boshier , and i.g .
hughes , . , ed . p.r .
berman ( academic press , boston , 1997 ) .
j. schmiedmayer , .
m. olshanii , .
d. s. petrov , g. v. shlyapnikov , and j. t. m. walraven , .
m. girardeau , .
girardeau , .
see particularly secs .
2 , 3 , and 6 .
a. lenard , .
a. g. rojo , g. l. cohen , and p. r. berman , .
m.d . girardeau and e.m .
k. bongs _ et al .
e. b. kolomeisky _
et al . _ .
m.d . girardeau and e.m .
vaidya and c.a .
w. ketterle and n. j. van druten , . o.
penrose and l. onsager , .
girardeau , e.m .
wright , and j.m .
triscari , .
e.m . wright , j. arlt , and k. dholakia , .
aitken , _ determinants and matrices _ ( oliver and boyd , edinburgh and london , 1951 ) , p. 112 . r.b .
laughlin , .
r.a . smith and n.k .
j. l. dubois and h. r. glyde , .
m.d . girardeau and e.m .
wright , cond - mat/0012365 ( 2000 ) .
w.p . reinhardt and c.w .
r. dum _ et al .
scott , r.j .
ballagh , and k. burnett , .
jackson , g.m .
kavoulakis , and c.j .
pethick , . s. burger _
et al . _ , .
j. denshlag _
et al . _ .
a.e . muryshev _
busch and j.r .
for a recent review see f. dalfovo _
et al . _ , . g. baym and c. j. pethick , .
y. kagan , g. v. shlapnikov , and j. t. m. walraven , .
m. edwards _ et al .
s. stringari , .
e. v. goldstein , k. plttner , and p. meystre , _ quantum .
semiclass . opt . _ * 7 * , 743 ( 1995 ) .
_ , _ science _ * 398 * , 218 ( 1999 ) .
s. a. morgan , r. j. ballagh , and k. burnett , .
w.p . reinhardt and c.w .
m. naraschewski _
et al . _ , .
j. javanainen and m. wilkens , ; _ ibid .
_ * 81 * , 1345 ( 1998 ) .
a. j. legget and f. sols , .
e. m. wright _
et al . _ , .
m.d . girardeau and e.m .
k. e. cahill and r. j. glauber .
j. durnin , j. j. miceli , and j. h. eberly , _ phys .
* 58 * , 1499 ( 1987 ) .
, to be published in _ phys .
a _ ( 2001 ) .
v. m. perez - garcia _ et al .
_ , _ phys .
a _ * 57 * , 3837 ( 1998 ) . y. s. kivshar and t. j. alexander , cond - mat/9905048 ( 1999 ) .
s.  l. cornish _
et  al . _ ,
* 85 * , 1795 ( 2000 ) . c. henkel and m. wilkens , _ europhys .
_ * 47 * , 414 ( 1999 ) . c. henkel , s. poetting , and m. wilkens , _ appl .
b _ * 69 * , 379 ( 1999 ) .
v. dunjko , v. lorent , and m. olshanii , cond - mat/0103085 ( 2001 ) .