coordinated motions and pattern formation have been studied for a wide range of biological organisms , from bacteria and amoebae to fish , from birds and wildebeest to humans ( ben - jacob et al . , 2000 ;
deneubourg et al .
, 2001 ; parrish et al . , 2002 ; partridge , 1982 ; wilson , 1975 ) .
animal groups often behave as if they have a single mind , displaying remarkable self - organized behavior . at one extreme ,
the individuals seem to need little information transfer ( e.g. fish schools ) , while at the other end the information exchange occurs in highly integrated ways through long - term associations among the individuals ( e.g. honeybee hives and human communities ) .
controlling such an organized behavior in groups of artificial objects , including autonomous underwater vehicles ( leonard et al . , 2006 ) and groups of autonomous agents ( jadbabaie et al . , 2003 ) , has received extensive attention in contemporary control theory .
challenges in both natural and engineering settings involve understanding which patterns emerge from the interaction among individual agents .
select laboratory experiments have shed some light on the schooling mechanism ( hunter , 1966 ; van olst and hunter , 1970 ; partridge and pitcher , 1979 , 1980 ; pitcher et al . , 1976 ) .
it still remains unclear , however , how the individual - level behavior and group - level ( `` macroscopic '' , or coarse - grained ) patterns are related .
more precise experiments using three - dimensional tracking of every individual in a population should lead to better understanding of this linkage .
an ultimate experimental study with precise control of every relevant detail may not be possible , yet appropriate mathematical models would provide a venue to establish behavioral cause , as one can consider different hypothetical individual - level interaction rules selectively ( see e.g. , flierl et al . , 1999 ) .
several different individual - based models have been proposed , which reproduce certain types of collective behavior in animal groups ( e.g. , see aoki , 1982 ; reynolds , 1987 ; deneubourg and goss , 1989 ) .
self - organization emerges also in a wide spectrum of physical and chemical systems , some of which ( e.g. , crystals and ferromagnetic materials ) exhibit apparent similarities with emergent patterns observed in animal groups .
vicsek et al .
( 1995 ) have introduced a discrete - time model of self - driven particles , or self - propelled particles ( spp ) , based on near - neighbor rules that are similar with those in the ferromagnetic xy model ( kosterlitz and thouless , 1973 ) .
the authors analyzed statistical properties of the model , including phase transition and scaling ( vicsek et al . , 1995 ) .
a long - range interaction has been incorporated into the spp model ( mikhailov and zanette , 1999 ) , and continuum , `` hydrodynamic '' versions of this model have been introduced ( toner and tu , 1995 , 1998 ; topaz et al . , 2006 )
. recently , couzin et al .
( 2002 , 2005 ) have introduced a model to provide insights into the mechanism of decision making in biological systems , which reproduces many important observations made in the field , and provides new insights into these phenomena . a review for various models can be found in parrish et al .
( 2002 ) and czirk and vicsek ( 2001 ) .
the models of couzin et al .
( 2005 ) , and most other such models , incorporate various detailed mechanistic steps .
these shed light on the role of leadership and imitation , and produce a number of surprising results , such as the influence that a few `` informed '' individuals can have on large collectives .
what is needed now are efforts to simplify those models , and to show especially what properties of the microscopic simulators are essential to explain that behavior .
for some models , closure schemes are available ( flierl et al . , 1999 ) ; but more generally , though we may suspect that closures exist , we can not derive explicit expressions for them . in such circumstances
, we need methods such as those used in this paper ; we perform the coarse - grained dynamical analysis by circumventing the derivation of governing equations , using an equation - free computational approach ( theodoropoulos et al . , 2000 ;
kevrekidis et al . , 2003 ) .
a particular goal is to understand how much of the specific spatial detail is fundamental to the behavior . but
turning to the kuramoto - type approximation , where the interaction is assumed to be global , we deliberately ignore local effects . to the extent that the models fail to explain observed types of behavior
, we will need to turn next to more detailed models .
most of previously proposed models concern populations of _ homogeneous _ ( or indistinguishable ) individuals .
furthermore , the dynamical analysis in the literature is often limited to a small subset of the entire parameter space , and a systematic classification of possible global dynamics remains elusive . in the current paper
, we study the _ coarse - grained _ alignment dynamics of individual - based animal group models .
the measurement of the mean angular deviation of fish schools ( e.g. clupeids and scombroids ; see atz , 1953 ; hunter 1966 ) showed that it varies continuously from no alignment to practically perfect alignment .
we account for this continuous change by heterogeneity ( `` quenched noise '' ; characterized by parameters of random variables drawn from a prescribed distribution function ) and the coupling strength .
our approach is flexible in that the heterogeneity can be introduced in various places in the model , and the way we analyze different heterogeneity cases does not require any significant modification .
the rest of the paper is organized as follows : models for homogeneous and heterogeneous animal groups are described in secs .
[ model1 ] and [ model_extension ] , and our approach , equation - free polynomial chaos , is explained in secs .
[ method1 ] and [ method2 ] .
coarse - grained dynamical analysis and its comparison with fine - scale dynamics , for a system of two informed individuals and a large number of heterogeneous uninformed individuals , are presented in sec .  [ results ] .
the case of two groups of heterogeneous informed individuals is presented in sec .
[ results2 ] .
we conclude with a brief discussion in sec .
[ conclusions ] .
we briefly discuss a `` minimal '' model proposed by nabet et al .
( 2006 ) , which we extend in our study .
it concerns the alignment dynamics of a homogeneous population of _ indistinguishable _ @xmath0 individuals with two subgroups of informed individuals ( `` leaders '' ) with populations @xmath1 and @xmath2 respectively and @xmath3 uninformed individuals ( `` followers '' ) , where @xmath4 : @xmath5 here @xmath6 characterizes the average direction of the individuals in each of the two informed subgroups for @xmath7 and the average direction of the uninformed individuals for @xmath8 .
@xmath9 is the corresponding informed , preferred direction ( @xmath10 can be set to zero without loss of generality ) and @xmath11 is the coupling strength .
this minimal model corresponds to the reduced system of the following system of @xmath0 individuals ( nabet et al . , 2006 ) :
@xmath12 where the angle @xmath13 characterizes the direction in which the @xmath14th individual is heading ( we will refer to it as `` orientation '' ) . the average direction @xmath6 is defined as the angle of the average of the phasors ( when each individual s dynamical state is considered as a phasor of unit radius and a phase angle ) of the individuals in the @xmath15th subgroup ; @xmath16 is the magnitude of the average of the phasors .
formally , this is written as @xmath17 the large population model in eq .
( [ m3 ] ) has a separation of time scales .
individuals within each subgroup synchronize quickly , i.e. , @xmath16 quickly converges to 1 .
the slow dynamics are described by the reduced system ( eq .
( [ reduced ] ) ) where the variables @xmath6 characterize the _ lumped behavior _ of each of the three subgroups .
it is assumed that the alignment ( orientational ) dynamics are independent of the translational counterpart ( sepulchre et al . ,
2005 ) ; hence , the dynamical state of an individual can be characterized by its orientation .
the functional form for mutual interaction is borrowed from the well - known kuramoto model ( kuramoto , 1984 ) , a prototypical model for coupled nonlinear oscillators .
this simplified global interaction model is consistent with an observation that the strongest correlations are observed between the ( speed and ) direction of the individual and the average ( speed and ) direction of the entire school ( partridge , 1982 ) : in the mean - field form of the kuramoto model , the interaction term can be rewritten as @xmath18 where @xmath19 in this alternate expression , the dependence on the difference between the individual direction and the average direction stands out explicitly . in the absence of coupling ( @xmath20 ) , each leader eventually heads for its preferred direction .
nontrivial dynamical behavior for the minimal model ( eq .  ( [ reduced ] ) ) are studied in nabet et al .
( 2006 ) ; bifurcations are analyzed for the global phase space in the case @xmath21 and @xmath22 .
the aforementioned models concern populations of _ homogeneous _ subgroups , where the individuals in each subgroup quickly synchronize , nearly perfectly ( @xmath23 ) , during the initial transients ( nabet et al . , 2006 ) . in the more general case ,
the mean angular deviation of fish schools is finite ( atz , 1953 ; hunter , 1966 ) , which is not captured in this `` minimal '' model .
we extend the model to account for the distribution of directions within schools , assuming it arises from the _ heterogeneity _ among the group members .
we introduce the heterogeneity in the following two ways : + ( i ) _ two leaders and many heterogeneous followers  _
first we consider the cases when the population consist of two leaders ( which possibly represent lumped behavior of groups of homogeneous leaders ) and @xmath24 followers :      where the heterogeneity is accounted for through the tendency to deviate from the average direction , characterized by @xmath26 , an i.i.d .
random variable drawn from a prescribed distribution function @xmath27 ( of standard deviation @xmath28 with mean value zero ) . for notational convenience ,
we drop a subscript of a variable to represent a random variable of a proper length ( _ cf .
_ @xmath26 and @xmath29 ) .
as @xmath10 can be set to zero without loss of generality , @xmath30 and @xmath31 are control parameters . in the current study
, we consider @xmath27 to be gaussian , but our analysis is not limited to this particular choice .
+ ( ii ) _ two groups of heterogeneous leaders  _ secondly , we consider two groups of heterogeneous leaders without any followers , focusing only on the dynamics among leaders .
the heterogeneity is accounted for by introducing randomness in the angles preferred by the leaders .
the orientations of the leaders in each group are denoted by @xmath32 s and @xmath33 s ( of sizes @xmath1 and @xmath2 ) respectively :      where the preferred angles @xmath35 and @xmath36 are randomly drawn from prescribed distributions @xmath37 and @xmath38 ( i.e. , i.i.d .
random variables of standard deviations @xmath39 and @xmath40 ) , respectively .
we set @xmath41 , and will vary @xmath42 and @xmath43 @xmath44)$ ] as control parameters ( and investigate some cases of different values of @xmath39 and @xmath40 in sec .
[ results2b ] ) .
direct integration of a system of two leaders ( open circles ; dashed lines indicate preferred angles ) and 300 followers ( dots ) , initialized from uniformly distributed orientations with randomly assigned heterogeneity variable ( i.e. , no initial correlations between @xmath45 and @xmath29 ) , is shown for an initial transient [ ( a ) to ( c ) ] , and for much longer time scales
[ ( d ) to ( f ) ] .
insets illustrate time evolution of the followers orientations on the @xmath46 plane , where strong correlations develop during a short time @xmath47 .
after that , the leaders and followers , the latter effectively as a `` unit '' , slowly drift to the stable steady state .
it takes of the order of @xmath48 for the system to asymptotically converges to this final state .
( @xmath49 ) .
, title="fig : " ]   direct integration of a system of two leaders ( open circles ; dashed lines indicate preferred angles ) and 300 followers ( dots ) , initialized from uniformly distributed orientations with randomly assigned heterogeneity variable ( i.e. , no initial correlations between @xmath45 and @xmath29 ) , is shown for an initial transient [ ( a ) to ( c ) ] , and for much longer time scales [ ( d ) to ( f ) ] .
insets illustrate time evolution of the followers orientations on the @xmath46 plane , where strong correlations develop during a short time @xmath47 .
after that , the leaders and followers , the latter effectively as a `` unit '' , slowly drift to the stable steady state .
it takes of the order of @xmath48 for the system to asymptotically converges to this final state .
( @xmath49 ) .
, title="fig : " ]   direct integration of a system of two leaders ( open circles ; dashed lines indicate preferred angles ) and 300 followers ( dots ) , initialized from uniformly distributed orientations with randomly assigned heterogeneity variable ( i.e. , no initial correlations between @xmath45 and @xmath29 ) , is shown for an initial transient [ ( a ) to ( c ) ] , and for much longer time scales [ ( d ) to ( f ) ] .
insets illustrate time evolution of the followers orientations on the @xmath46 plane , where strong correlations develop during a short time @xmath47 .
after that , the leaders and followers , the latter effectively as a `` unit '' , slowly drift to the stable steady state .
it takes of the order of @xmath48 for the system to asymptotically converges to this final state .
( @xmath49 ) .
, title="fig : " ]     direct integration of a system of two leaders ( open circles ; dashed lines indicate preferred angles ) and 300 followers ( dots ) , initialized from uniformly distributed orientations with randomly assigned heterogeneity variable ( i.e. , no initial correlations between @xmath45 and @xmath29 ) , is shown for an initial transient [ ( a ) to ( c ) ] , and for much longer time scales [ ( d ) to ( f ) ] .
insets illustrate time evolution of the followers orientations on the @xmath46 plane , where strong correlations develop during a short time @xmath47 .
after that , the leaders and followers , the latter effectively as a `` unit '' , slowly drift to the stable steady state .
it takes of the order of @xmath48 for the system to asymptotically converges to this final state .
( @xmath49 ) .
, title="fig : " ]   direct integration of a system of two leaders ( open circles ; dashed lines indicate preferred angles ) and 300 followers ( dots ) , initialized from uniformly distributed orientations with randomly assigned heterogeneity variable ( i.e. , no initial correlations between @xmath45 and @xmath29 ) , is shown for an initial transient [ ( a ) to ( c ) ] , and for much longer time scales [ ( d ) to ( f ) ] .
insets illustrate time evolution of the followers orientations on the @xmath46 plane , where strong correlations develop during a short time @xmath47 .
after that , the leaders and followers , the latter effectively as a `` unit '' , slowly drift to the stable steady state .
it takes of the order of @xmath48 for the system to asymptotically converges to this final state .
( @xmath49 ) .
, title="fig : " ]   direct integration of a system of two leaders ( open circles ; dashed lines indicate preferred angles ) and 300 followers ( dots ) , initialized from uniformly distributed orientations with randomly assigned heterogeneity variable ( i.e. , no initial correlations between @xmath45 and @xmath29 ) , is shown for an initial transient [ ( a ) to ( c ) ] , and for much longer time scales [ ( d ) to ( f ) ] .
insets illustrate time evolution of the followers orientations on the @xmath46 plane , where strong correlations develop during a short time @xmath47 .
after that , the leaders and followers , the latter effectively as a `` unit '' , slowly drift to the stable steady state .
it takes of the order of @xmath48 for the system to asymptotically converges to this final state .
( @xmath49 ) .
, title="fig : " ]      the kuramoto model , a paradigm for all - to - all , phase - coupled oscillator models , has been extensively studied and used to shed light on many synchronization phenomena ( kuramoto , 1984 ; acebrn et al . , 2005 , and references therein ) .
this model has the property that , in the full synchronization regime ( of large enough @xmath31 values ) , phase angles become quickly correlated with ( or  sorted " according to ) the natural frequencies during the initial short transients ( moon et al . ,
similar correlations ( between the angles and heterogeneity random variables ) are expected to arise in the current model ( which is indeed the case , as will be shown later in fig .
[ correlation ] ) , since the coupling term is qualitatively similar .
as in moon et al .
( 2006 ) , we choose expansion coefficients in wiener polynomial chaos as coarse - grained `` observables '' , to explore low - dimensional , coarse - grained dynamics .
wiener(-hermite ) polynomial chaos was introduced by wiener ( 1938 ) , who represented a random process in terms of functional expansions of wiener process ( historically , this method has been termed as polynomial `` chaos '' , because of its initial usage on homogeneous chaos , such as turbulence and brownian motion , rather than the nature of the method ) .
ghanem and spanos ( 1991 ) later extended this idea to treat random processes as functional expansions of random variables , or elements in the hilbert space of random functions , in which a spectral representation in terms of polynomial chaos is identified .
the projections ( or coefficients ) on the polynomial base then can be determined through a galerkin approach .
this method was subsequently applied in uncertainty quantification of various problems ( e.g. , see ghanem , 1999 ; jardak et al . , 2002 ) , and has been extended to general situations using the askey scheme ( xiu et al . , 2002 ; now known as generalized polynomial chaos ) .    in this method ,
dependent random variables ( @xmath45 of the followers for the case ( i ) , and @xmath50 and @xmath51 for the case ( ii ) ) are expanded in polynomials of independent random variables ( @xmath29 , or @xmath52 and @xmath53 ) using appropriately chosen basis functions .
details for the two cases are as follows : + ( i ) _ two leaders and many heterogeneous followers  _ for convenience , we introduce the unit gaussian random variable @xmath54 . using this newly defined variable , we expand @xmath55 ( i.e. @xmath56 ) in hermite polynomials of @xmath57 [ @xmath58 = 1 , @xmath59 , @xmath60 , @xmath61 : @xmath62 where @xmath63 is the highest order retained in the truncated series , @xmath64 is the @xmath65th hermite polynomial , and the @xmath66 s are the expansion coefficients which will be referred to simply as `` chaos coefficients '' in this paper . wiener polynomial chaos , utilizing hermite polynomials as basis functions , is the appropriate choice for gaussian random variables that we consider in the present study .
the probability density function of the gaussian random variables appears as the weighting function of hermite polynomials , and the hermite polynomial expansion is suggested to converge exponentially for gaussian processes ( lucor et al . , 2001 ) . for other random variables ,
use of different basis functions ( for instance , legendre polynomials for uniform random variables ) has been suggested for fast convergence , which is the basis of the development of the generalized polynomial chaos ( xiu et al . , 2002 ) .
we choose the first few nonvanishing chaos coefficients @xmath66 s , as well as the orientations of the leaders ( @xmath67 and @xmath68 ) , to be the coarse - grained `` observables '' . due to symmetry ,
all the even order @xmath66 s vanish , except for the zeroth order @xmath69 that corresponds to the average direction of the followers .
geometrically , @xmath70 and @xmath71 respectively represent a measure for the linear order spread of the angles ( the `` slope '' between @xmath45 and @xmath29 ) and the cubic order measure . in the continuum limit ( @xmath72 )
, the chaos coefficients can be exactly determined using the orthogonality relations for hermite polynomials . however , in the _ finite _ cases of single realization we consider , @xmath73 , those relations hold only approximately , and the coefficients are evaluated using least squares fitting , following moon et al .
+ ( ii ) _ two groups of heterogeneous leaders  _ in the second case , we expand @xmath50 and @xmath51 in terms of @xmath52 and @xmath53 , respectively : @xmath74 where the chaos coefficients and are the coarse `` observables '' of our choice , @xmath64 s are hermite polynomials ( for gaussian @xmath75 and @xmath76 ) , and @xmath77 and @xmath78 are unit gaussian random variables . +      a prerequisite to _ coarse - grained _ dynamical analysis ( which is the main goal of the current study ) is , in a traditional sense , an explicit derivation of coarse - grained governing equations . in principle
, such equations for chaos coefficients , in the continuum limit ( @xmath72 ) , might be obtained through a stochastic galerkin method ( ghanem and spanos , 1991 ) .    in the present study , we do not even attempt to derive such equations . we circumvent their derivation by using an equation - free multiscale computational approach ( theodoropoulos et al . , 2000 ; kevrekidis et al . , 2003 , 2004 ) .
this approach enables us to explore the coarse - grained dynamics without the assumption of the continuum limit .
the premise of this approach is that coarse - grained governing equations conceptually _ exist _ , but are not explicitly available in closed form .
this approach is based on the recognition that short bursts of appropriately initialized microscopic ( fine - scale ) simulations during a time horizon @xmath79 and the projection of the results onto coarse - grained variables , say , result in time - steppers ( mappings ) for those variables @xmath80 ( which is effectively the same as the discretization of unavailable equations ) : @xmath81 one then processes the results of the short simulations to estimate various coarse - grained quantities ( such as time derivatives , action of jacobians , residuals ) to perform relevant coarse - grained level numerical computations , as if those quantities were obtained from coarse - grained governing equations . for instance , one can _ integrate _
unavailable governing equations in time ( coarse projective integration ; see below ) , or compute the steady states of the above coarse time - stepper , by utilizing fixed point algorithms ( such as newton - raphson or newton - gmres ) .
equation - free computations consist of the following steps :    * identify coarse - grained variables ( `` coarse observables '' ) that sufficiently describe both the micro- and macroscopic dynamics ; in our study , they are @xmath66 s ( and @xmath82 s ) . for convenience , we denote the microscopic ( macroscopic ) descriptions by @xmath83 ( @xmath84 ) .
* choose an appropriate _ lifting _ operator @xmath85 , which maps @xmath84 to one ( or more ) consistent description(s ) @xmath83 ( for the purposes of variance reduction and ensemble - averaging ) .
in the current study , this can be achieved by using the relations in eqs .
( [ expansion ] ) and ( [ expansion2 ] ) ; once random variables are drawn , these relations are used to obtain corresponding @xmath83 . * starting from lifted initial condition(s ) @xmath86 , run the microscopic simulator to obtain @xmath87 at a later time ( @xmath88 ) .
* use an appropriate _ restriction _
operator @xmath89 ( least squares fitting , in the current study ) which maps the microscopic state(s ) to the macroscopic description @xmath90 , which effectively results in time series of coarse observables , or their coarse time - stepper @xmath80 ; @xmath91 . * apply desired numerical techniques using the coarse - grained variables obtained from the step 4 . and
repeat some of the above steps as needed .
an extensive discussion can be found in kevrekidis et al .
( 2003 , 2004 ) .
direct integration of the `` fine - scale '' model of eq .
( [ microeq1 ] ) in the strong coupling regime ( @xmath92 ) , started from randomly assigned orientations and the heterogeneity variable ( the latter is a gaussian random variable ) , illustrates that a strong correlation between @xmath45 and @xmath29 develops during a short , initial transient time ; the orientations of the followers quickly become a monotonically increasing function of their heterogeneity variable ( fig .
[ correlation ] ) , after which they slowly drift as a `` unit '' until they settle down in the final steady state . during the latter slow drift
, the system can be described as two leaders and a _ single _
`` clump '' of followers , whose coarse - grained states can be successfully described by a small number of chaos coefficients .
a similar time scale separation exists in the model of homogeneous populations . in this case
, followers quickly collapse asymptotically to the _ same _ direction ( nabet et al . , 2006 ) .
( color online ) accelerated computation of stable steady states via coarse projective integration using five coarse - grained variables , shown here for two different time scales ( @xmath49 ) .
initially all the values are assigned to be 0 .
both @xmath70 and @xmath71 reach their steady state values relatively quickly ( see ( a ) ) , while the others are slowly varying ( see ( b ) ; they are still varying at @xmath93 ) .
dots represent the time intervals during which short direct integration is performed ( and restricted ) , in the course of the projective integration using forward euler method .
solid lines represent the trajectories of direct full integration during the entire time
. higher efficiency can be achieved by optimally choosing the time horizon for the direct integration , the projection stepsize , and projection method .
, title="fig : " ]   ( color online ) accelerated computation of stable steady states via coarse projective integration using five coarse - grained variables , shown here for two different time scales ( @xmath49 ) .
initially all the values are assigned to be 0 .
both @xmath70 and @xmath71 reach their steady state values relatively quickly ( see ( a ) ) , while the others are slowly varying ( see ( b ) ; they are still varying at @xmath93 ) .
dots represent the time intervals during which short direct integration is performed ( and restricted ) , in the course of the projective integration using forward euler method .
solid lines represent the trajectories of direct full integration during the entire time .
higher efficiency can be achieved by optimally choosing the time horizon for the direct integration , the projection stepsize , and projection method .
, title="fig : " ]      we begin by accelerating the approach to a stable steady state using an equation - free algorithm , the coarse projective integration method ( gear and kevrekidis , 2003 ) .
in contrast to a conventional , direct integration of the full fine - scale model during the _ entire _ time ( until sufficient convergence to stable , final states ) , this method exploits smoothness in the coarse variables ( estimated through a direct integration during a _ short _ time ) , in order to extrapolate and take a large projective time - step ( compared to the original integration time - step size ) .
this saves computational effort .
the procedure consists of ( @xmath94 ) _ lifting _ ( appropriate initialization of the fine - scale simulator , an integrator of eq .
( [ microeq1 ] ) , consistent with prescribed coarse - grained values ) , ( @xmath95 ) _ direct integration _ of the microscopic simulator during a relatively short time interval ( but long enough to accurately estimate local coarse - grained time derivatives ) , ( @xmath96 ) _ restriction _ ( of fine - scale description onto coarse - grained variables ) , and ( @xmath97 ) _ taking a projective step _ ( using a traditional numerical integration scheme such as forward euler ) .
the computational payoff of this method depends on the ratio between a short direct integration time interval , the projective time - step size , and the computational effort required for lifting / restriction procedures ( see e.g. , rico - martinez et al . , 2004 ) .
more importantly , successful computation of steady states through this method naturally attests to the _ validity _ of the chosen coarse - grained observables in describing _ both _ fine - scale and coarse - grained states .
( left panel ) a bifurcation diagram observed on one leader @xmath68 , computed using auto2000 ( @xmath98 ) . solid (
dashed ) lines represent stable ( unstable ) branches . there exist a few other unstable branches that are not shown here . at some critical value of @xmath30 ,
an unstable state in the upper branch undergoes a forward pitchfork bifurcation ; two unstable states coincide .
the lower branch of `` trivial '' solutions does not exhibit any bifurcation .
( right panels ) snapshots of two ( symmetric ) stable states in the bistable regime ( @xmath99 ) , marked by dots in the left panel . ]    [ cols="^,^,^,^,^,^,^",options="header " , ]     projective integration using five coarse - grained variables ( @xmath100 , and the first three non - vanishing @xmath66 s ; @xmath69 , @xmath70 , and @xmath71 ) follows virtually the same trajectories of the full , direct integration ( fig .  [ pi ] ) , even if @xmath29 is _ newly _ drawn at each lifting ; the agreement is even better if the same @xmath29 were used ( hence the dynamics are fully deterministic ) . both _ lifting _ ( simply using eq .
( [ expansion ] ) ) and _ restriction _ ( least squares fitting ) operations require minimal computational efforts .
therefore , the computational efficiency in the present case is nearly exclusively determined by the projective step size , which is a factor of about four in fig .
[ pi ] ; with a more sophisticated projection algorithm , a higher efficiency can be obtained .
we see that both @xmath70 and @xmath71 reach their steady state values quickly ( @xmath101 ) , showing that the correlation between @xmath45 and @xmath29 are fully developed by then .
however , the other chaos coefficient @xmath69 ( representing the average direction ) slowly drifts towards the steady state , and so do @xmath67 and @xmath68 ( note that it is still varying at @xmath93 ) ; the computation of an asymptotic , steady state requires a very long time integration .    direct integration ( including projective integration )
can not compute unstable steady states and are inappropriate for stability computations and parametric bifurcation studies .
both stable and _ unstable _ steady state values can be systematically ( and much more efficiently than the projective integrations ) computed by applying coarse - grained fixed point algorithms to the steady state condition of eq .
( [ coarse_stepper ] ) , i.e. , @xmath102 , in _ much lower _ dimension than that of individual level .
we use the coarse newton - gmres ( kelley , 1995 ) , a matrix - free , method to compute coarse - grained fixed points .
we observe that the algorithm accurately converges within a few steps ( tab .
[ gmres ] ) ; the converged values are accurately consistent with the restricted values of the fixed point solution of the detailed ( i.e. , ( @xmath0 + 2)-dimensional ) problem , within prescribed convergence tolerance . by combining a coarse fixed point algorithm with pseudo - arclength continuation ( keller , 1987 ) ,
we numerically compute coarse - grained bifurcation diagrams in the following sections .
the computational efficiency of the coarse fixed point algorithm varies with the choice of the initial guess for the iteration . with a totally uneducated guess
, it could take even longer than the direct integration ( note that the latter never computes the exact solutions ) ; however , during the continuation computation shown below , a good initial guess is always available from the previous parameter value(s ) , and even several orders of magnitude of computational efficiency can be achieved .
( a ) a bifurcation diagram observed on one ( arbitrarily chosen ) follower , as a function of @xmath30 ( @xmath98 ) , computed using auto2000 ( the same case as in fig .
[ micro ] ) .
superscript ` f ' has been added to emphasize that this is the orientation of _ a follower_. a few other existing unstable branches are not included here .
the upper branch undergoes a pitchfork bifurcation and becomes stable .
( b ) a coarse bifurcation diagram observed on @xmath69 ( average direction ) , obtained by the coarse newton - gmres method with pseudo arc - length continuation . only a blowup around the bifurcation point
coarse - grained dynamics exhibit the same structure as in the fine - scale level .
filled ( open ) circles represent stable ( unstable ) steady states .
, title="fig : " ]   ( a ) a bifurcation diagram observed on one ( arbitrarily chosen ) follower , as a function of @xmath30 ( @xmath98 ) , computed using auto2000 ( the same case as in fig .  [ micro ] ) .
superscript ` f ' has been added to emphasize that this is the orientation of _ a follower_. a few other existing unstable branches are not included here .
the upper branch undergoes a pitchfork bifurcation and becomes stable .
( b ) a coarse bifurcation diagram observed on @xmath69 ( average direction ) , obtained by the coarse newton - gmres method with pseudo arc - length continuation . only a blowup around the bifurcation point
coarse - grained dynamics exhibit the same structure as in the fine - scale level .
filled ( open ) circles represent stable ( unstable ) steady states .
, title="fig : " ]      we first analyze the detailed @xmath103-dimensional fine - scale model in the full synchronization regime , in order to obtain insights on fine - scale dynamics to be compared with our coarse - grained analysis below .
we use auto2000 ( doedel et al . , 2000 ) to compute the fine - scale bifurcation diagrams as functions of @xmath30 at a fixed value of @xmath31 ; only projections for one leader ( @xmath68 ) are shown in fig .
[ micro ] and for one follower in fig .
[ macroth2 ] ( a ) .
all the other followers exhibit essentially the same dynamical behavior as the one shown here ( except for some quantitative differences ) .
the interaction between the individuals causes the steady state directions of the leaders to deviate from the preferred angles @xmath104 and @xmath30 , respectively .
such deviation can occur in two directions , either toward the region bounded by @xmath105 $ ] ( an `` obvious '' steady state where followers are directed in between the directions of the leaders ; see fig .
[ micro ] ( b ) ) or the other way around ( e.g. , fig .  [ micro ] ( a ) ) .
the analysis shows that for small @xmath30 values only the former state is stable , while for large values , both of these states become stable .
the branches for `` obvious '' stable steady states , which correspond to lower straight solid lines , exhibit no bifurcations ( see figs .  [ micro ] and [ macroth2 ] ( a ) ) . on the other branches , forward pitchfork bifurcations at some critical value of @xmath30 give birth to another stable branch ( a state on this stable branch
is shown in fig .
[ micro ] ( a ) ) , as well as two unstable asymmetric solution branches , hence the population becomes bistable .
the critical value of @xmath30 for the onset of the bistability depends on @xmath31 ( precisely speaking , @xmath106 ) ; the critical value is @xmath107 at @xmath108 .
as @xmath31 decreases further , the critical value monotonically increases until fully synchronized steady states lose stability at some critical value of @xmath31 .
we now compute coarse - grained steady state solutions . a coarse - grained bifurcation diagram for @xmath69 ( representing the average direction of the followers )
is compared with the corresponding diagram observed for _ one _ follower , in figs .
[ macroth2 ] ( b ) and ( a ) ; ( b ) is a blowup of the region around the bifurcation .
both sides of the bifurcation point can be described by the _
same _ set of coarse - grained observables , which clearly summarize group level dynamical behavior of the followers before and after the bifurcation .
as @xmath31 decreases in the kuramoto model , oscillators get desynchronized ( kuramoto , 1984 ) , starting with the oscillator with the maximum value of @xmath109 ( the `` extreme '' oscillator ) through a saddle - node ( actually a `` sniper '' ) bifurcation on a limit cycle ( moon et al . , 2006 ) .
we expect the same type of bifurcation to occur in this model . however , when we try to compute the coarse - grained steady states as functions of @xmath31 using the previously mentioned five coarse variables ( via coarse newton - gmres method and pseudo arc - length continuation , neither a bifurcation nor an unstable branch is appropriately identified .
the computation , initialized at large @xmath31 steady states , accurately follows stable branches down to some critical value of @xmath31 ( where the transition occurs ) , and then fails to converge .
our coarse - grained observables are not sufficient to describe the states on the `` other side '' of the bifurcation point , as we will explain below .
( color online ) ( a ) a bifurcation diagram observed on a few followers , including the one with the maximum of @xmath109 ( the `` extreme '' follower ) , as a function of @xmath31 ( @xmath110 ) , computed using auto2000 . a critical value where the extreme individual gets desynchronized corresponds to a saddle - node bifurcation point on a limit cycle ( a `` sniper '' bifurcation ) . except for the extreme follower ,
stable and unstable branches nearly coincide ( see the inset ) .
( b ) in order to capture the fine - scale bifurcation , the angle of the extreme follower has to be discounted from the chaos expansion and considered as an _ extra _ coarse - grained variable ( see text ) .
we distinguish these chaos coefficients ( from the ones used so far ) by adding a prime .
it was computed via the coarse newton - gmres method with continuation .
, title="fig : " ]   ( color online ) ( a ) a bifurcation diagram observed on a few followers , including the one with the maximum of @xmath109 ( the `` extreme '' follower ) , as a function of @xmath31 ( @xmath110 ) , computed using auto2000 . a critical value where the extreme individual gets desynchronized corresponds to a saddle - node bifurcation point on a limit cycle ( a `` sniper '' bifurcation ) . except for the extreme follower ,
stable and unstable branches nearly coincide ( see the inset ) .
( b ) in order to capture the fine - scale bifurcation , the angle of the extreme follower has to be discounted from the chaos expansion and considered as an _ extra _ coarse - grained variable ( see text ) .
we distinguish these chaos coefficients ( from the ones used so far ) by adding a prime .
it was computed via the coarse newton - gmres method with continuation .
, title="fig : " ]    a fine - scale bifurcation diagram ( computed using auto2000 ) obtained by starting from a stable steady state on the lower branch in fig .
[ micro ] is shown in fig .
[ turningpt ] ( a ) . here
the diagrams for two leaders and only a few followers , including the extreme one , are shown .
we find that both stable and unstable branches for each angle nearly coincide for all the individuals ( see inset of fig .  [ turningpt ] ( a ) ) , except for the extreme one .
as the difference between stable and unstable branches ( at the same value of @xmath31 ) is appreciable _ only _ when observed on this extreme oscillator , a smooth mapping between @xmath45 and @xmath29 does not prevail for unstable states , and the previously used chaos coefficients are not appropriate any more .    taking these observations into account , it is easy to remedy the situation as follows : the fact that stable and unstable branches nearly coincide , discounting the extreme follower , suggests that all the individuals _ except for the extreme follower _ can be again described by the same set of chaos coefficients .
thus we treat the orientation of the extreme one separately ( introducing it as an additional coarse - grained variable ) , and discount it from the polynomial chaos expansion .
( from the fact that the extreme follower gets desynchronized at the transition , one can also intuitively see that followers have to be considered as a combination of a clump of synchronized `` bulk '' and a separate , extreme one . )
we compute the solutions with continuation , using this new set of _ six _ coarse variables , which captures the bifurcation and appropriately describes the unstable steady states ( fig .
[ turningpt ] ( b ) ) ; we have analyzed exactly the same realization used in fig .
[ turningpt ] ( a ) for direct comparison .
when bifurcation diagrams are computed for ensembles of many realizations , an uncertainty will arise in the exact quantification of the bifurcation point , due to the fluctuation of finite - dimensional random variables among realizations , while the results are qualitatively the same as those of a single realization ( xiu et al . , 2005 ) .
the coarse bifurcation results shown in fig .
[ turningpt ] ( b ) illustrate that the steady state directions of the leaders and the average direction of followers ( @xmath111 , discounting the extreme one ; a prime is added to distinguish it from the previously used notation ) are virtually the same for a range of @xmath31 .
only higher order chaos coefficients ( only @xmath112 is shown in fig .
[ turningpt ] ) appreciably vary as a function of @xmath31 , which means that individuals spread more widely as @xmath31 decreases , until the extreme one eventually starts to oscillate freely , while the average steady state direction remains the same .
lines : a bifurcation diagram observed on @xmath67 in the minimal model ( the first two odes in eq .
( [ reduced ] ) with @xmath113 ) , for varying @xmath30 at a fixed value of @xmath114 , obtained by auto2000 . for large enough preferred angles ( @xmath115 ) , the system becomes bistable through a forward pitchfork bifurcation .
circles : a coarse bifurcation diagram near the pitchfork bifurcation , observed on @xmath69 , as a function of @xmath116 , computed by the coarse newton - gmres method . filled ( open )
circles represent stable ( unstable ) steady states . ]
here we explore both the fine - scale and coarse - grained dynamics of a model for two groups of _ heterogeneous _ leaders ( with no followers ) shown in eq .
( [ for_many ] ) , and compare the results of the two different scales .
one notable difference from the kuramoto model is that `` oscillators '' in eq .
( [ for_many ] ) do not have finite random variables ( natural frequencies ) , hence there is no onset of the synchronization that occurs at a finite value of @xmath31 ( or , they can be alternatively seen as kuramoto - like oscillators of zero natural frequencies , which result in the onset at @xmath117 , hence they get synchronized for all @xmath31 values ) .
the analysis of the minimal model ( the first two of eq .
( [ reduced ] ) with @xmath118 ) reveals that for large enough @xmath119 the system exhibits bistability for a certain range of @xmath31 ( nabet et al . , 2006 ) , as in the previous case in sec .
[ results ] . here
we will vary @xmath43 as the main parameter for two different values of @xmath31 . for large coupling strengths ( @xmath120 ) ,
the bistability in the minimal model appears through a forward pitchfork bifurcation , when @xmath30 is varied as a parameter ( fig .
[ microk2.4 ] ) .
this minimal model can be seen as a special case of the current model , where both @xmath52 and @xmath53 are assumed to be delta functions and each group consists of identical individuals .
coarse - grained bifurcation diagrams observed on the first two chaos coefficients : ( a ) @xmath69 ; the average direction of the first group of leaders , and ( b ) @xmath70 ; the `` slope '' between @xmath50 and @xmath52 , as functions of @xmath43 .
these are blowups of the region around the forward pitchfork bifurcation point in fig .
[ microk2.4 ] .
, title="fig : " ]   coarse - grained bifurcation diagrams observed on the first two chaos coefficients : ( a ) @xmath69 ; the average direction of the first group of leaders , and ( b ) @xmath70 ; the `` slope '' between @xmath50 and @xmath52 , as functions of @xmath43 .
these are blowups of the region around the forward pitchfork bifurcation point in fig .
[ microk2.4 ] .
, title="fig : " ]    we begin by asking whether our model for heterogeneous groups exhibits similar types of dynamical behavior .
one can also do accelerated computations of steady states using the coarse projective integration , but here we skip such computations and present only the coarse bifurcation analysis results . _
coarse _ bifurcation diagrams obtained through the coarse newton - gmres method ( kelley , 1995 ) and pseudo arc - length continuation ( keller , 1987 ) ( for gaussian distributions of @xmath52 and @xmath53 ; @xmath121 , @xmath122 ) show that the heterogeneous groups indeed exhibit the same qualitative type of coarse dynamical behavior around the pitchfork bifurcation point ( fig .
[ macrok2.4 ] ) . as we consider symmetric unimodal distribution functions , all the even order chaos coefficients ( except for @xmath69 and @xmath123 ) virtually vanish .
the diagram for @xmath69 of the first group @xmath124 ( average direction ) exhibits reasonably good quantitative agreement with the corresponding diagram for the minimal model , within fluctuations of finite - size random variables , shown in figs .
[ microk2.4 ] and [ micromacrok1.8 ] .
it is interesting to note that at the critical point , all the followers are headed for the same direction ( @xmath125 , which corresponds to the `` slope '' between @xmath52 and @xmath50 ) ; see fig .
[ macrok2.4 ] ( b ) .
lines : a bifurcation diagram for the minimal model of two leaders , for varying @xmath30 at a fixed value of @xmath126 , obtained by auto2000 . for large enough preferred angles ( @xmath127 ) ,
the system becomes bistable , but the nature of the bifurcation is different from that of higher @xmath31 value cases ( a saddle - node vs. a pitchfork bifurcation ; see fig .
[ microk2.4 ] ) .
circles : a coarse bifurcation diagram observed on the average direction of the first group of leaders ( @xmath69 ) around the saddle - node bifurcation , as a function of @xmath116 , computed via the coarse newton - gmres method with continuation . filled ( open ) circles represent stable ( unstable ) steady states . ]
( color online ) coarse - grained bifurcation diagrams near a turning point in fig .
[ micromacrok1.8 ] , for @xmath52 distributions of three different widths ( the standard deviation @xmath39 = 0.1 for circles ; 0.2 for squares ; 0.3 for triangles ) , obtained via the coarse newton - gmres method and continuation .
the standard deviation for the second group , @xmath40 , is kept the same at 0.1 ( @xmath128 ) . filled and open symbols represent stable and unstable states , respectively .
( a ) the first chaos coefficients @xmath69 ( average direction of the first group ) are nearly the same for the three cases .
the difference between the cases becomes apparent in higher order coefficients that reflect the degree of spreading ; see @xmath70 in ( b ) .
, title="fig : " ]   ( color online ) coarse - grained bifurcation diagrams near a turning point in fig .
[ micromacrok1.8 ] , for @xmath52 distributions of three different widths ( the standard deviation @xmath39 = 0.1 for circles ; 0.2 for squares ; 0.3 for triangles ) , obtained via the coarse newton - gmres method and continuation .
the standard deviation for the second group , @xmath40 , is kept the same at 0.1 ( @xmath128 ) . filled and open symbols represent stable and unstable states , respectively .
( a ) the first chaos coefficients @xmath69 ( average direction of the first group ) are nearly the same for the three cases .
the difference between the cases becomes apparent in higher order coefficients that reflect the degree of spreading ; see @xmath70 in ( b ) .
, title="fig : " ]    the hermite polynomial expansion converges so quickly that the expansions can be accurate even when truncated at the third order .
due to the reflection symmetry ( about @xmath43/2 ) , coefficients have similar structures as the ones , after proper reflection and translation . only results on
are presented here . as the coupling strength decreases across @xmath129 , the nature of the bifurcation changes ( from a pitchfork ) to a saddle - node bifurcation ( fig .
[ micromacrok1.8 ] ) at @xmath129 , which also occurs in the model for homogeneous populations ; the nature of the transition between these different bifurcations , a higher codimension bifurcation , has been discussed in nabet et al .
( 2006 ) .
so far we have considered statistically similar groups , namely @xmath21 and @xmath130 ; they differed only by average preferred directions .
it is natural to ask how the dynamics change as the parameters concerned with the distributions ( for the preferred directions ) are varied .
it is readily expected that the essential dynamics of two different - size groups can be reflected in the minimal model using two different coupling strengths , which is considered in nabet et al .
here we consider only the cases with varying width of the distributions ( @xmath131 ) , which has no analog in the minimal model .
coarse bifurcation diagrams for three different gaussian distributions for @xmath52 ( @xmath39 is varied while @xmath40 is kept at 0.1 ; see fig .
[ vary_std ] ) show that the average directions ( @xmath69 s ) hardly vary with the width of the distributions ; the primary parameter that affects on the average direction is the _ group size_. for the distributions of different widths , the fixed point computation with continuation fails to converge at different values of @xmath69 s ; points marked by arrows in fig .
[ vary_std ] are the last points the newton - gmres computations converged in each case , when approached from the stable branches . such a failure of convergence can be expected , because the steady states on this unstable branch overlap with another nearby unstable branch ( which is not shown in this figure , but was shown in fig .
[ micromacrok1.8 ] ) ; characterizing the distribution with a few wiener chaos coefficients does not provide an accurate description any more .
the differences between the three cases ( of different distribution widths ) manifest themselves clearly in higher order chaos coefficients .
while the average behavior remains nearly the same ( fig .  [ vary_std ] ( a ) ) , individuals in the group spread more widely ( as reflected in @xmath70 and higher order coefficients ; fig .
[ vary_std ] ( b ) ) , as the width of the random variable ( distribution ) increases .
we have demonstrated a computational venue ( an equation - free polynomial chaos approach ) to study coarse - grained dynamics of individual - based models accounting for the _ heterogeneity _ among the individuals in animal group alignment models .
we considered _ finite _ populations of ( i ) two `` leaders '' ( which have direct knowledge on preferred directions ) and @xmath132 uninformed , heterogeneous `` followers '' , and ( ii ) two groups of heterogeneous `` leaders '' .
we explored the _ coarse - grained _ , group level ( low - dimensional ) dynamics using the polynomial chaos expansion coefficients as coarse - grained observables ; these observables account for rapidly developing correlations between random variables , and sufficiently specify both fine - scale and coarse - grained ( group - level ) dynamical states .
all the analysis in our study was done expressively avoiding the derivation of coarse - grained governing equations , following a _ nonintrusive _ , equation - free computational approach wrapped around the direct system simulator .
it should be noted that we have not assumed that @xmath0 is infinitely large ( so - called the `` continuum limit '' ) .
our approach can be used for systems of any finite , large number of populations , and it can be equally applied to various types of random variables ( following generalized polynomial chaos ) and/or various heterogeneity .
we compared our results with those of minimal models that do not account for heterogeneity among the individuals .
they show good agreement in the lowest order ( i.e. , average directions ) , which clearly highlights the correspondence between the individual- and group - level dynamics ( figs .
[ microk2.4 ] and [ micromacrok1.8 ] ) .
indeed this implies that the results in nabet et al .
( 2006 ) , where no heterogeneity is explicitly accounted for , are more robust than demonstrated in that paper alone .    in order to analyze different coarse - grained bifurcations
, it became necessary to use different sets of coarse - grained variables , even if the model is _ the same _ in the fine - scale level ( fig .
[ turningpt ] ) .
this clearly shows that an appropriate choice of coarse - grained observables ( in terms of which one can obtain useful closures ) is an essential step ; different coarse - grained observables are required , as the same fine - scale model closes differently .    in the present study
, we assumed that the orientational dynamics can be separated from their translational counterpart , and considered the simplest nontrivial cases of all - to - all ( `` all - visible '' ) , sinusoidal coupling .
our future work will involve the incorporation of the translational dynamics and more complicated coupling / network topology , including heterogeneous couplings .
our work presented here is the first step of our effort toward the development of more detailed ( and biologically more plausible ) models and their _ coarse - graining_.
s.j.m . and i.g.k .
were financially supported by doe and nsf grant ef0434319 .
s.a.l . was supported in part by nsf grant ef0434319 and darpa grant hr00110510057 .
b.n . and n.e.l .
were supported in part by onr grants n000140210826 and n000140410534 .
atz , j. w. , 1953 .
orientation in schooling fishes , in proceedings of a conference on orientation in animals .
office of naval research , pp .
103 - 130 .
department of the navy , washington , d. c. , section 2 .
kevrekidis , i. g. , gear , c. w. , hyman , j. m. , kevrekidis , p. g. , runborg , o. , theodoropoulos , c. , 2003 .
equation - free coarse - grained multiscale computation : enabling microscopic simulators to perform system - level tasks .
sciences 1 ( 4 ) , 715 - 762 ; e - print physics/0209043 .
leonard , n. e. , paley , d. , lekien , r. , sepulchre , r. , fratantoni , d. m. , davis , r. , 2006 .
collective motion , sensor networks and ocean sampling .
proceedings of the ieee , special issue on `` the emerging technology of networked control systems '' , to appear .
lucor , d. , xiu , d. , karniadakis , g. , 2001 .
spectral representations of uncertainty in simulations : algorithms and applications . in proceedings of the international conference on spectral and high order methods ( icosahom-01 ) , uppsala , sweden .
nabet , b. , leonard , n. e. , couzin , i. d. , levin , s. a. , 2006 .
leadership in animal group motion : a bifurcation analysis .
proceedings of the 17th international symposium on mathematical theory of networks and systems , kyoto , japan , july 2006 .
sepulchre , r. , paley , d. , leonard , n. e. , 2005 .
collective motion and oscillator synchronization , in cooperative control : a post - workshop volume , 2003 block island workshop on cooperative control , eds . v. kumar , n. leonard and a.s .
morse , springer - verlag .


a nucleus is at the heart of every atom , and loosely speaking , is a collection of protons and neutrons that interact pairwise , with much smaller , but significant , three - body interactions .
we are fortunate to know that the underlying laws governing the strong interactions result from a quantum field theory called quantum chromodynamics ( qcd ) .
it is constructed in terms of quark and gluon fields with interactions determined by a local su(3 ) gauge - symmetry and , along with quantum electrodynamics ( qed ) , underpins all of nuclear physics when the five relevant input parameters , the scale of strong interactions @xmath1 , the three light - quark masses @xmath2 , @xmath3 and @xmath4 , and the electromagnetic coupling @xmath5 , are set to their values in nature .
it is remarkable that the complexity of nuclei emerges from `` simple '' gauge theories with just five input parameters .
perhaps even more remarkable is that nuclei resemble collections of nucleons and not collections of quarks and gluons . by solving qcd , we are expecting to predict , with arbitrary precision , nuclear processes and the properties of multi - baryon systems .
the fine - tunings observed in the structure of nuclei , and in the interactions between nucleons , are peculiar and fascinating aspects of nuclear physics .
for the values of the input parameters that we have in our universe , the nucleon - nucleon ( nn ) interactions are fine - tuned to produce unnaturally large scattering lengths in both s - wave channels ( described by non - trivial fixed - points in the low - energy effective field theory ( eft ) ) , and the energy levels in the @xmath6be - system , @xmath7c and @xmath8 are in `` just - so '' locations to produce enough @xmath7c to support life , and the subsequent emergence and evolution of the human species . at a fundamental level
it is important for us to determine the sensitivity of the abundance of @xmath7c to the light - quark masses and to ascertain the degree of their fine - tuning .
being able to solve qcd for the lightest nuclei , using the numerical technique of lattice qcd ( lqcd ) , would allow for a partial unification of nuclear physics .
it would be possible to `` match '' the traditional nuclear physics techniques - the solution of the quantum many - body problem for neutrons and protons using techniques such as no - core shell model ( ncsm ) , greens function monte carlo ( gfmc ) , and others , to make predictions for the structure and interactions of nuclei for larger systems than can be directly calculated with lqcd . by placing these calculations on a fundamental footing , reliable predictions with quantifiable uncertainties
can then be made for larger systems .
during the 1990 s , the nuclear forces were systematized by the hierarchy emerging from the spontaneously broken chiral symmetries of qcd .
the resulting small expansion parameters are powers of the external momenta and powers of the light - quark masses normalized to the scale of chiral symmetry breaking , as pioneered by weinberg , first in the meson sector and then the multi - nucleon sector  @xcite . in addition to generating nuclear forces that are consistent with qcd
, this construction provides the calculational advantage of parametric estimates of the systematic uncertainty introduced by the truncation of the nuclear interactions at a given order in the expansion .
the actual ordering of contributions remains a subject of debate even today , with weinberg s chiral expansion of the potential having its peculiar difficulties , as does the ksw expansion of scattering amplitudes  @xcite .
calculations are being performed at a sufficiently high order where the size of truncation errors are quite small .
weinberg s ordering of operators based upon a chiral expansion of the n - body potentials between nucleons has been carried out to n@xmath9lo , which includes contributions to the three - body ( starting at n@xmath10lo ) and the leading four - body interactions ( starting at n@xmath9lo ) ( for a recent review see ref .
@xcite ) .    during the last several years
, nuclear structure calculations have been performed with the chiral nuclear forces , leading to both postdictions and predictions for nuclei to a given order in the expansion , and compared with experiment , e.g. see fig .
[ fig : ncsma ] .
li and @xmath6b using chiral nuclear forces  @xcite .
[ image is reproduced with the permission of p. maris.],title="fig:",scaledwidth=50.0%]li and @xmath6b using chiral nuclear forces  @xcite .
[ image is reproduced with the permission of p. maris.],title="fig:",scaledwidth=50.0% ]    the nuclear forces that are presently used in such calculations are constrained by experimental measurements of nn scattering and light nuclei . as the desired precision increases , which requires working to higher orders in the expansion , the number of required experimental constraints increases .
eventually , there are too few experimental constraints to practically reduce the systematic uncertainty below some level in any given calculation . however , lqcd calculations are expected to provide a way to constrain the nuclear forces beyond what is possible with experiment , and hence to further reduce the systematic uncertainties in nuclear structure calculations . beyond providing direct calculations of important quantities , lqcd calculations of the light nuclei and nuclear forces can    1 .   verify experimental constraints and/or reduce the uncertainties in the constraints imposed by experiment , 2 .
constrain components of the nuclear forces that are inaccessible to experiment , for instance the light - quark mass dependences which dictates some of the multi - pion vertices , and multi - neutron forces , 3 .   constrain counterterms at higher orders in the expansion to further reduce the systematic uncertainties
lqcd is a technique in which space - time is discretized into a four - dimensional grid and the qcd path integral over the quark and gluon fields at each point in the grid is performed in euclidean space - time using monte carlo methods .
a lqcd calculation of a given quantity will deviate from its value in nature because of the finite volume of the space - time ( with @xmath11 lattice points ) over which the fields exist , and the finite separation between space - time points ( the lattice spacing , @xmath12 ) .
however , such deviations can be systematically removed by performing calculations in multiple volumes with multiple lattice spacings , and extrapolating using the theoretically known functional dependences on each .
supercomputers are needed for such calculations due to the number of space - time points and the monte carlo evaluation of the path integral over the dynamical fields . in order for a controlled continuum extrapolation ,
the lattice spacing must be small enough to resolve structures induced by the strong dynamics , encapsulated by @xmath13 where @xmath14 is the scale of chiral symmetry breaking .
further , in order to have the hadron masses , and also the scattering observables , exponentially close to their infinite - volume values , the lattice volume must be large enough to contain the lightest strongly interacting particle , encapsulated by @xmath15 where @xmath16 is the mass of the pion and @xmath17 is the extent of the spatial dimension of the cubic lattice volume ( this , of course , can be generalized to non - cubic volumes ) .
effective field theory ( eft ) descriptions of these observables exist for @xmath18 ( the symanzik action and its translation into chiral perturbation theory ( @xmath19pt ) and other frameworks ) and @xmath15 ( the p - regime of @xmath19pt and other frameworks ) .
the low - energy constants in the appropriate eft are fit to the results of the lqcd calculations , which are then used to take the limit @xmath20 and @xmath21 .
computational resources devoted to lqcd calculations are becoming sufficient to be able to perform calculations at the physical values of the light quark masses in large enough volumes and at small enough lattice spacings to be relevant , but the majority of present day calculations are performed with pion masses of @xmath22 . therefore , most calculations require the further extrapolation of @xmath23 , but do not yet include strong isospin breaking or electromagnetism . in principle , the gauge - field configurations that are generated in lqcd calculations can be used to calculate an enormous array of observables , spanning the range from particle to nuclear physics . in practice , this is becoming less common , largely due to the different scales relevant to particle physics and to nuclear physics .
calculations of quantities involving the pion with a mass of @xmath24 are substantially different from those of , say , the triton with a mass of @xmath25 , and with the typical scale of nuclear excitations being @xmath26 .
present day dynamical lqcd calculations of nuclear physics quantities are performed with @xmath27 , lattice spacings of @xmath28 and volumes with spatial extent of @xmath29 .
lqcd calculations are approached in the same way that experimental efforts use detectors to measure one or more quantities - the computer is equivalent to the accelerator and the algorithms , software stack , and parameters of the lqcd calculation(s ) are the equivalent of the detector .
the parameters , such as lattice spacing , quark masses and volume , are selected based upon available computational resources , and simulations of the precision of the calculation(s ) required to impact the physical quantity of interest , i.e. simulations of the lqcd monte carlos are performed .
the size of the computational resources required for cutting edge calculations are such that you only get `` one shot at it '' . a typical work - flow of a lqcd calculation consists of three major components .
the first component is the production of an ensemble of gauge - field configurations which contain statistically independent samplings of the gluon fields resulting from the lqcd action .
the production of gauge - fields requires the largest partitions on the leadership class computational facilities , typically requiring @xmath30k compute cores .
present - day calculations have @xmath31 dynamical light - quark flavors and use the wilson , @xmath32-improved - wilson , staggered ( kogut - susskind ) , domain - wall or overlap discretizations , each of which have their own `` features '' .
it is the evaluation of the light - quark determinant ( the determinant of a sparse matrix with dimensions @xmath33 ) that consumes the largest fraction of the resources .
roughly speaking , @xmath34 hybrid monte carlo ( hmc ) trajectories are required to produce an ensemble of @xmath35 decorrelated gauge fields , but in many instances this is an under estimate . for observables involving quarks ,
a second component of production is the determination of the light - quark propagators on each of the configurations .
the light - quark propagator from a given source point ( an example of which is shown in fig .
[ fig : prop ] ) is determined by an iterative inversion of the quark two - point function , using the conjugate - gradient ( cg ) algorithm or variants thereof such as bicgstab , or the most recently developed multi - grid ( mg ) .
[ fig : prop],scaledwidth=60.0% ]    during the last couple of years , the propagator production codes have been ported to run on gpu machines in parallel .
gpu s can perform propagator calculations faster than standard cpu s by an order of magnitude , and have led to a major reduction in the statistical uncertainties in many calculations .
there have been numerous algorithm developments that have also reduced the resources required for propagator production , such as the implementation of deflation techniques and the use of multi - grid methods .
the third component of a lqcd calculation is the production of correlation functions from the light - quark propagators .
this involves performing all of the wick contractions that contribute to a given quantity .
the number of contractions required for computing a single hadron correlation function is small .
however , to acquire long plateaus in the effective mass plots ( emps ) that persist to short times , lscher - wolff type methods involve the computation of a large number of correlation functions resulting from different interpolating operators , and the number of contractions can become large .
in contrast , the naive number of contractions required for a nucleus quickly becomes astronomically large ( @xmath36 for uranium ) , but symmetries in the contractions , and new algorithms ( e.g. ref .
@xcite ) greatly reduce the number of operations that must be performed . a further consequence of the hierarchy of mass scales is that there is an asymptotic signal - to - noise problem in nuclear correlation functions .
the ratio of the mean value of the correlation function to the variance of the sample from which the mean is evaluated degrades exponentially at large times .
however , this is absent at short and intermediate times and the exponential degradation of the signal - to - noise in the correlation functions can be avoided .
capability computing resources provided by leadership class computing facilities are used to produce ensembles of gauge - field configurations , while capacity computing resources , both those operated by usqcd and elsewhere are used to perform observable - dependent calculations of correlation functions using these configurations .
thus the capability resources enable a multitude of physics calculations to be accomplished with the capacity resources . in the area of cold nuclear physics
there is currently a well - defined set of goals , and a program in place to accomplish these goals , as described in one of the 2013 usqcd whitepapers  @xcite : * hadron structure * , * hadron spectroscopy * , * hadronic interactions , nuclear forces and nuclei * , and * fundamental symmetries*.      before calculations of nuclei can be sensibly undertaken , the mass and structure of the nucleon must be reproduced in lqcd calculations .
the spectrum of the lowest - lying hadrons calculated with lqcd is shown in fig .
[ fig : spectrum ] , from which we observe that indeed lqcd is postdicting all of the light - hadron masses within uncertainties .
[ fig : spectrum],scaledwidth=70.0% ]    beyond its mass , one property of the nucleon that is well known experimentally is the forward - matrix element of the isovector axial current , @xmath37 .
significant effort has been put into calculating @xmath37 with lqcd , a summary of which is shown in fig .
[ fig : ga ] , but the extrapolated lqcd value has consistently been smaller than the experimental value . with calculations beginning to be performed at the physical pion mass ,
the community is focused on understanding and quantifying the systematic uncertainties in these calculations .    .
[ image is reproduced with the permission of h .- w .
[ fig : ga ] , scaledwidth=70.0% ]    a central element of the physics program at jlab is to determine the excited spectra of mesons and baryons , including searching for exotic states that are beyond the naive nonrelativistic quark model of hadrons , but arise naturally in qcd .
a critical component of this program is the lqcd calculations of the spectra
. they will play a central role in interpreting and understanding the experimental measurements .
the spectra of such states is complicated by the presence of open multi - hadron channels and significant formal developments remain to be put in place before rigorous statements about the spectra can be made .
calculations at unphysical pion masses have been performed by the jlab lqcd group , examples of which are shown in fig .
[ fig : isoscalars ] , and remarkable progress has been made in the identification of states in these calculations .     by the jlab lqcd group  @xcite
[ image is reproduced with the permission of r. edwards . ]
[ fig : isoscalars ] , scaledwidth=85.0% ]    the aim is to have lqcd predict the exotic spectra of hadrons before , or at the same time as , the gluex experiment at jlab runs , targeting the 2018 milestone hp15 .
multi - hadron lqcd calculations are significantly more challenging than single - hadron calculations for a number of reasons , and systems involving baryons are even more challenging .
meson - meson systems are the simplest multi - hadron systems , and impressive progress has been made in the recent past , particularly when the lqcd calculations are combined with @xmath19pt .
there is little or no signal - to - noise problem in such calculations and therefore highly accurate lqcd calculations of stretched - isospin states can be performed with modest computational resources . moreover , the efts which describe the low - energy interactions of pions and kaons , including lattice - spacing and finite - volume effects , have been developed to non - trivial orders in the chiral expansion .     scattering  @xcite .
[ image in the left panel is reproduced with the permission of h. leutwyler . ]
[ fig : pipi2 ] , title="fig:",scaledwidth=56.0% ]   scattering  @xcite .
[ image in the left panel is reproduced with the permission of h. leutwyler . ]
[ fig : pipi2 ] , title="fig:",scaledwidth=34.0% ]    the @xmath38 pion - pion ( @xmath39 ) scattering length serves as a benchmark calculation with an accuracy that can only be aspired to in other systems .
the scattering lengths for @xmath40 scattering in the s - wave are uniquely predicted at lo in @xmath19pt  @xcite : @xmath41 while experiments do not directly provide stringent constraints on the scattering lengths , a determination of s - wave @xmath40 scattering lengths using the roy equations has reached a remarkable level of precision  @xcite : @xmath42 the roy equations  @xcite use dispersion theory to relate scattering data at high energies to the scattering amplitude near threshold . at present , lqcd can compute @xmath40 scattering only in the @xmath38 channel with precision as the @xmath43 channel contains disconnected diagrams which require large computational resources .
it is of great interest to compare the precise roy equation predictions with lqcd calculations , and fig .
[ fig : pipi2 ] summarizes theoretical and experimental constraints on the s - wave @xmath40 scattering lengths  @xcite .
this is clearly a strong - interaction process for which theory has somewhat out - paced the challenging experimental measurements .
mixed - action @xmath44 lqcd calculations , employing domain - wall valence quarks on a rooted staggered sea and combined with mixed - action @xmath19pt , have predicted  @xcite @xmath45 at the physical pion mass . the agreement between this result and the roy equation determination is a striking confirmation of the lattice methodology , and a powerful demonstration of the constraining power of chiral symmetry in the meson sector . however , lqcd calculations at one or more smaller lattice spacings , and with different discretizations , are required to verify and further refine this calculation .
the etm collaboration has performed a @xmath46 calculation of the @xmath38 @xmath40 scattering length  @xcite , producing a result extrapolated to the physical pion mass of @xmath47    it is interesting to compare the pion mass dependence of the meson - meson scattering lengths to the current algebra predictions . in fig .
[ fig : caplots ] ( left panel ) one sees that the @xmath38 @xmath40 scattering length is consistent with the current algebra result up to pion masses that are expected to be at the edge of the chiral regime in the two - flavor sector . while in the two - flavor theory one expects fairly good convergence of the chiral expansion and
, moreover , one expects that the effective expansion parameter is small in the channel with maximal isospin , the lqcd calculations clearly imply a degree of cancellation between chiral logs and counterterms .
however , as one sees in fig .
[ fig : caplots ] ( right panel ) , the same phenomenon occurs in @xmath48 where the chiral expansion is governed by the strange quark mass and is therefore expected to be much more slowly converging .
this remarkable conspiracy between chiral logs and counterterms for the meson - meson scattering lengths remains mysterious .
vs. @xmath49 ( left panel ) and @xmath50 vs. @xmath51 ( right panel ) .
the solid ( red ) curves are the current algebra predictions .
[ fig : caplots ] , title="fig:",scaledwidth=45.0% ]   vs. @xmath49 ( left panel ) and @xmath50 vs. @xmath51 ( right panel ) . the solid ( red ) curves are the current algebra predictions . [
fig : caplots ] , title="fig:",scaledwidth=45.0% ]    lqcd calculations of the meson - meson scattering phase - shifts are much less advanced than of the scattering length .
this is because the calculation of the phase shift , @xmath52 , at a given energy , @xmath53 , requires a lqcd calculation of the two - meson correlation function at the energy @xmath53 .
generally speaking , a given calculation can determine the lowest few two - hadron energy eigenvalues for a given momentum of the center - of - mass , and that multiple lattice volumes will allow for additional values of @xmath53 at which to determine @xmath52 .
the first serious calculation of the s - wave ( @xmath54 ) @xmath38 @xmath40 phase - shift was done by the cp - pacs collaboration with @xmath46 at a relatively large pion mass  @xcite , and more recently two groups have performed calculations at lower pion masses  @xcite , the results of which are shown in fig .
[ fig : pipidelta ] .     scattering phase - shift .
the left panel shows the results of the lqcd calculations below the inelastic threshold ( @xmath55 ) at a pion mass of @xmath56  @xcite .
the vertical ( blue ) line denotes the start of the t - channel cut .
the shaded region in the right panel shows the results of the lqcd calculation extrapolated to the physical pion mass using nlo @xmath19pt , while the points and uncertainties corresponds to the existing experimental data .
the vertical ( red ) line corresponds to the inelastic threshold .
[ fig : pipidelta ] , title="fig:",scaledwidth=45.0% ]   scattering phase - shift .
the left panel shows the results of the lqcd calculations below the inelastic threshold ( @xmath55 ) at a pion mass of @xmath56  @xcite .
the vertical ( blue ) line denotes the start of the t - channel cut .
the shaded region in the right panel shows the results of the lqcd calculation extrapolated to the physical pion mass using nlo @xmath19pt , while the points and uncertainties corresponds to the existing experimental data .
the vertical ( red ) line corresponds to the inelastic threshold .
[ fig : pipidelta ] , title="fig:",scaledwidth=45.0% ]    further , in some nice work by the hadron spectrum collaboration ( hsc ) , the first efforts have been made to extract the d - wave ( @xmath57 ) @xmath38 @xmath40 phase shift  @xcite .
one of the more exciting recent results is the mapping out of the @xmath58-resonance at @xmath56 from the @xmath59 energy - levels using luscher s method , as shown in fig .
[ fig : rho ]  @xcite .
-resonance at a pion mass of @xmath56  @xcite [ image is reproduced with the permission of r. edwards . ]
[ fig : rho ] , scaledwidth=70.0% ]      calculations of the nucleon - nucleon scattering lengths have been successfully underway for the last decade  @xcite for a range of pion masses .
recently , lqcd calculations have been performed at @xmath60 that also provide the effective ranges  @xcite , the results of which are shown in fig .
[ fig : nnscattnplqcd ] .     channel extracted from lqcd calculations at the su(3 ) symmetric point , including the fit to the ere at n@xmath10lo .
the right panel shows the ratio of the scattering length to effective range , a quantity that is a measure of the naturalness of the system . [
fig : nnscattnplqcd ] , title="fig:",scaledwidth=50.0% ]    channel extracted from lqcd calculations at the su(3 ) symmetric point , including the fit to the ere at n@xmath10lo .
the right panel shows the ratio of the scattering length to effective range , a quantity that is a measure of the naturalness of the system .
[ fig : nnscattnplqcd ] , title="fig:",scaledwidth=50.0% ]    also shown are fits to the effective range expansion ( ere ) , including the shape parameter .
the scattering length and effective range in the @xmath61 channel determined from the nlo fit to the ere are @xmath62 the shape parameter obtained from the nnlo fit to the ere expansion is : @xmath63 .
an interesting aspect of this result is that the ratio of scattering length to effective range , a measure of the naturalness of the system , is @xmath64 , which is to be compared with @xmath65 at the physical quark masses .
this leads one to speculate that the deuteron might be unnatural over a large range of quark masses and not just close to the physical values , indicating that it is not finely tuned .
this speculation requires precise calculations at lighter quark masses to determine if this is , in fact , the situation .
perhaps some of the most important lqcd calculations of late are those of the ground states of the light nuclei , including the deuteron , @xmath9he , @xmath66he and light hypernuclei .
[ fig : ball ] shows the binding energy of the deuteron , @xmath9he and @xmath66he  @xcite as a function of the pion mass .
he ( middle panel ) and @xmath66he ( right panel ) binding energies from @xmath44 lqcd calculations  @xcite .
, title="fig:",scaledwidth=32.0% ]   he ( middle panel ) and @xmath66he ( right panel ) binding energies from @xmath44 lqcd calculations  @xcite .
, title="fig:",scaledwidth=32.0% ]   he ( middle panel ) and @xmath66he ( right panel ) binding energies from @xmath44 lqcd calculations  @xcite .
, title="fig:",scaledwidth=32.0% ]    not only is it exciting to see nuclei emerge from qcd for a range of the light - quark masses , such calculations are crucial in dissecting and refining the chiral nuclear forces . however , it is clear that calculations at lighter pion masses are required , including at the physical pion mass .    )
calculated by nplqcd  @xcite at a pion mass of @xmath60 . , scaledwidth=70.0% ]    a summary of the energy - levels at the flavor su(3 ) symmetry point found in the s - shell nuclei and hypernuclei  @xcite is shown in fig .  [
fig : nuclearsummary ] .
these energy levels are elements of su(3 ) irreps . which allowed , in some cases , e.g. the h - dibaryon , the hypertriton and @xmath67 , for distinct energy levels with the same spin and parity to be determined .
such calculations will become somewhat more complicated at lighter quark masses when the up and down quarks are not degenerate with the strange quark .
the calculations of nplqcd and those of yamazaki _ et al _ are already shedding light on how the ground - state energies of the light nuclei approach their values at the physical light - quark masses .
they are all bound at the heavier light - quark masses and become less bound as the quarks become lighter . in the case of the dineutron , which is bound at @xmath60 ,
it becomes unbound at some intermediate value of the pion mass , giving rise to a neutron - neutron system with an infinite scattering length .
one of the interesting aspects of the nuclear forces to explore is the tensor interaction , responsible for the mixing between the s - wave and d - wave channels in the deuteron channel .
there are a series of lqcd calculations that can be performed that will permit an extraction of the sd mixing parameter , @xmath68 using lscher s method  @xcite , see ref .
@xcite .
it is important to understand how a few of the past roadblocks to progress in this area have been recently overcome .
one of the roadblocks of the past was / is the `` signal - to - noise problem '' that afflicts states other than the pion .
this problem is seen most simply in the single - nucleon correlation function , generated with a three - quark source and a three - quark sink .
the variance of this correlation function is dictated by a 3-quark 3-anti - quark source and a 3-quark 3-anti - quark sink , which overlaps with both the @xmath69 and @xmath70 intermediate states ( and all others with the appropriate quantum numbers ) . at large times , the variance correlation function is dominated by the @xmath70 intermediate state , while the single nucleon correlation function is dominated by the single nucleon , giving rise to an exponentially degrading signal . however , at intermediate times , the behavior of the `` signal - to - noise '' is determined by the overlap of the variance sinks and sources onto the intermediate hadronic states .
the momentum projection onto single nucleon blocks , that nplqcd is currently using , provides a volume suppression of the @xmath70 intermediate state compared to the @xmath69 state .
thus , there is an intermediate time interval in which the signal - to - noise ratio is not exponentially degrading .
it is in this time interval , dubbed the `` golden window '' , that plateaus for the low - lying energy levels in light nuclei can be identified .
unfortunately , the window shrinks as the number of nucleons is increased , and so further developments will be required to go to much larger nuclei .
a second roadblock that inhibited progress in lqcd calculations of nuclei was the number of wick contractions required to form a correlation function .
a system containing @xmath71 up quarks and @xmath72 down quarks requires @xmath73 wick contractions , which is a rapidly growing number as one moves beyond the nucleon .
it was recognized that recursion relations relating the wick contractions in systems with n mesons can be related to those with @xmath74 mesons  @xcite .
further , somewhat more sophisticated algorithms  @xcite have been developed for the multi - baryon systems that greatly reduce the computing resources required to perform the contractions .
these have led to very efficient calculations of the s - shell nuclei and hypernuclei , moving beyond the s - shell requires extensions of these works , and new ideas are required to calculate heavier nuclei .
one of the points of discussion that came up during this presentation was how to optimally couple the results of lqcd calculations to nuclear structure calculations .
given the expertise in the nuclear structure community , it makes little sense for lqcd theorists to `` go it alone '' and attempt to calculate the entire periodic table .
it makes much more sense for the lqcd theorists to produce sets of quantities that can be handed to the nuclear structure theorists who use them in their machinery to determine the periodic table .
the question is what are the optimal quantities to pass along from lqcd .
it seems that the minimal set of quantities that could be passed along are the energy eigenvalues for a given system .
lqcd calculations of the energy spectrum of an a - nucleon system could be performed in multiple lattice volumes , with multiple lattice spacings and at multiple light - quark masses , and handed to the the nuclear structure theorists who in turn reproduce the energies by tuning the chiral interactions .
these tuned interactions are then used to calculate processes in the continuum .
this methodology was used to calculate the @xmath75 interactions at the physical pion mass using @xmath19pt  @xcite .
the chiral interactions were tuned to reproduce the finite - volume energy levels determined in a series of lqcd calculations , and then used to calculate the scattering phase shift at the physical pion mass .
progress in this direction is starting to be made , as demonstrated in recent calculations by nir barnea and collaborators  @xcite , by using the ground state energies of the deuteron , dineutron and @xmath9he at @xmath60 to reproduce the @xmath66he ground state using the pionless eft .
i have summarized the rapid progress that is being made in developing lqcd into a reliable calculational tool for low - energy nuclear physics .
it holds the promise to directly connect the structure and properties of nuclei with qcd , and to enable a refinement of the chiral nuclear forces that are used as input into nuclear structure calculations . at present , the ground states of the s - shell nuclei and hypernuclei are being calculated at unphysically heavy light - quark masses , but within the next few years , such calculations at @xmath76 will be performed ( if hardware and software resources increase as expected ) . within the next five years
, the spectrum and interactions of the lightest nuclei and hypernuclei will be postdicted or predicted with fully - quantified uncertainties .
it is worth emphasizing that the lqcd effort in the us relies heavily on scidac funding to support the scientists who develop and optimize the software to run on the rapidly evolving computational hardware , e.g. gpu - accelerated compute nodes that comprise titan at ornl , or the bg / qs at anl and llnl .
further , the effort requires ongoing access to both capability computing resources on leadership - class computing facilities , and capacity computing obtained from nersc , xsede , through usqcd and at local compute clusters . ongoing software ( see fig .  [
fig : mgsoftware ] ) and hardware support are critical to progress in this area .
ideally , one would start with a lqcd calculation and predict all of the quantities of interest in low - energy nuclear physics .
presently , we are not in a position to do this , even if significantly more computing resources were provided to the program .
while lscher provided the formalism to relate the two - body s - matrix directly to two - particle energy levels inside a cubic volume with the fields subject to periodic boundary conditions  @xcite , which has since been understood and generalized to the two - nucleon systems , e.g. ref .
@xcite , such formalism is complicated to apply in coupled - channels systems  @xcite .
further , the formalism is not in place for the three- and higher - body sectors , but progress is being made in such systems  @xcite .
_ * happy birthday james * : james vary is one of the first nuclear theorists i met when i arrived in the united states to enter the phd program at caltech in the mid 1980 s .
i recall james taking the time to talk physics with me during his stay .
his detailed knowledge of , and passion for , important problems of the day left a lasting impression on me . despite having been able to chat with , and even collaborate with , james since that time ,
when i learned that this conference was in part to celebrate james s 70th birthday , i was taken aback as it seems like yesterday that he was in his early 40 s ( and i was in my early 20 s ) , and he has retained the same passion and energy for science .
i should also add that james is responsible for me remembering the value of @xmath77 ! happy 70@xmath78 ! ! _
d.  b.  kaplan , m.  j.  savage and m.  b.  wise , phys .
b * 424 * , 390 ( 1998 ) [ nucl - th/9801034 ] .
d.  b.  kaplan , m.  j.  savage and m.  b.  wise , nucl .
b * 534 * , 329 ( 1998 ) [ nucl - th/9802075 ] .
e.  epelbaum and u.  -g .
meissner , ann .
nucl .  part .
* 62 * , 159 ( 2012 ) [ arxiv:1201.2136 [ nucl - th ] ] .
p.  maris , j.  p.  vary and p.  navratil , phys .
c * 87 * , 014327 ( 2013 ) [ arxiv:1205.5686 [ nucl - th ] ] .
w.  detmold and k.  orginos , phys .
d * 87 * , 114512 ( 2013 ) [ arxiv:1207.1452 [ hep - lat ] ] .
a.  s.  kronfeld ,
arxiv:1209.3468 [ physics.hist-ph ] . h.  -w .
lin , pos lattice * 2012 * , 013 ( 2012 ) [ arxiv:1212.6849 [ hep - lat ] ] .
j.  j.  dudek , r.  g.  edwards , b.  joo , m.  j.  peardon , d.  g.  richards and c.  e.  thomas , phys .
d * 83 * , 111502 ( 2011 ) [ arxiv:1102.4299 [ hep - lat ] ] .
s.  weinberg , phys .
* 17 * , 616 ( 1966 ) .
g.  colangelo , j.  gasser and h.  leutwyler , nucl .
b * 603 * , 125 ( 2001 ) [ arxiv : hep - ph/0103088 ] .
h.  leutwyler , pos c * confinement8 * , 068 ( 2008 ) [ arxiv:0812.4165 [ hep - ph ] ] .
s.  m.  roy , phys .
b * 36 * , 353 ( 1971 ) .
s.  r.  beane , t.  c.  luu , k.  orginos , a.  parreno , m.  j.  savage , a.  torok and a.  walker - loud , phys .
d * 77 * , 014505 ( 2008 ) [ arxiv:0706.3026 [ hep - lat ] ] .
x.  feng , k.  jansen and d.  b.  renner , phys .
b * 684 * , 268 ( 2010 ) [ arxiv:0909.3255 [ hep - lat ] ] .
t.  yamazaki _ et al .
_ [ cp - pacs collaboration ] , phys .
d * 70 * , 074513 ( 2004 ) [ arxiv : hep - lat/0402025 ] .
j.  j.  dudek , r.  g.  edwards , m.  j.  peardon , d.  g.  richards and c.  e.  thomas , phys .
d * 83 * , 071504 ( 2011 ) [ arxiv:1011.6352 [ hep - ph ] ] .
s.  r.  beane _ et al . _
[ nplqcd collaboration ] , arxiv:1107.5023 [ hep - lat ] .
j.  j.  dudek , r.  g.  edwards and c.  e.  thomas , phys .
d * 87 * , no . 3 , 034505 ( 2013 ) [ arxiv:1212.0830 [ hep - ph ] ] .
s.  r.  beane , e.  chang , s.  d.  cohen , w.  detmold , p.  junnarkar , h.  w.  lin , t.  c.  luu and k.  orginos _ et al .
_ , phys .
c * 88 * , 024003 ( 2013 ) [ arxiv:1301.5790 [ hep - lat ] , arxiv:1301.5790 [ hep - lat ] ] .
m.  fukugita , y.  kuramashi , h.  mino , m.  okawa and a.  ukawa , phys .
* 73 * , 2176 ( 1994 ) [ arxiv : hep - lat/9407012 ] .
m.  fukugita , y.  kuramashi , m.  okawa , h.  mino and a.  ukawa , phys .
d * 52 * , 3003 ( 1995 ) [ arxiv : hep - lat/9501024 ] .
s.  r.  beane , p.  f.  bedaque , k.  orginos and m.  j.  savage , phys .
lett .  * 97 * , 012001 ( 2006 ) [ arxiv : hep - lat/0602010 ] .
s.  r.  beane _ et al . _
[ nplqcd collaboration ] , phys .
d * 81 * , 054505 ( 2010 ) [ arxiv:0912.4243 [ hep - lat ] ] .
n.  ishii , s.  aoki and t.  hatsuda , phys .
lett .  * 99 * , 022001 ( 2007 ) [ arxiv : nucl - th/0611096 ] .
s.  aoki , t.  hatsuda and n.  ishii , comput .
* 1 * , 015009 ( 2008 ) [ arxiv:0805.2462 [ hep - ph ] ] .
s.  aoki , t.  hatsuda and n.  ishii , prog .
theor .  phys .
* 123 * , 89 ( 2010 ) [ arxiv:0909.5585 [ hep - lat ] ] .
t.  yamazaki , y.  kuramashi , a.  ukawa , arxiv:1105.1418 [ hep - lat ] .
t.  yamazaki , y.  kuramashi , a.  ukawa , phys .
d * 81 * , 111504 ( 2010 ) [ arxiv:0912.1383 [ hep - lat ] ] .
de forcrand and m.  fromm , phys .
lett .  * 104 * ,
112005 ( 2010 ) [ arxiv:0907.1915 [ hep - lat ] ] . s.  r.  beane _ et al . _
[ nplqcd collaboration ] , phys .
d * 85 * , 054511 ( 2012 ) [ arxiv:1109.2889 [ hep - lat ] ] .
t.  inoue _ et al .
[ hal qcd collaboration ] , nucl .
a * 881 * , 28 ( 2012 ) [ arxiv:1112.5926 [ hep - lat ] ] . s.  r.  beane , e.  chang , s.  d.  cohen , w.  detmold , h.  w.  lin , t.  c.  luu , k.  orginos and a.  parreno _ et al .
_ , arxiv:1206.5219 [ hep - lat ]
. t.  yamazaki , k.  -i .
ishikawa , y.  kuramashi and a.  ukawa , phys .
d * 86 * , 074514 ( 2012 ) [ arxiv:1207.4277 [ hep - lat ] ] .
m.  luscher , commun .
phys .   * 105 * , 153 ( 1986 ) .
m.  luscher , nucl .
b * 354 * , 531 ( 1991 ) .
r.  a.  briceno , z.  davoudi and t.  c.  luu , phys .
d 88 , * 034502 * ( 2013 ) [ phys .
d * 88 * , 034502 ( 2013 ) ] [ arxiv:1305.4903 [ hep - lat ] ] .
w.  detmold and m.  j.  savage , phys .
d * 82 * , 014511 ( 2010 ) [ arxiv:1001.2768 [ hep - lat ] ] .
t.  doi and m.  g.  endres , comput .
commun .   * 184 * , 117 ( 2013 ) [ arxiv:1205.0585 [ hep - lat ] ] .
s.  r.  beane , e.  chang , s.  d.  cohen , w.  detmold , h.  -w .
lin , t.  c.  luu , k.  orginos and a.  parreno _ et al .
_ , phys .
lett .   * 109 * , 172001 ( 2012 ) [ arxiv:1204.3606 [ hep - lat ] ] .
r.  a.  briceno , z.  davoudi , t.  luu and m.  j.  savage , arxiv:1309.3556 [ hep - lat ] . n. barnea , _ nuclear physics with heavy pions - eft for lqcd _ , a talk presented at the institute for nuclear theory , july 31 , 2013 .
j.  c.  osborn , r.  babich , j.  brannick , r.  c.  brower , m.  a.  clark , s.  d.  cohen and c.  rebbi , pos lattice * 2010 * , 037 ( 2010 ) [ arxiv:1011.2775 [ hep - lat ] ] . at the computational nuclear physics meeting , washington d.c .
, july 23 - 24 , 2012 , by b. joo .
m.  t.  hansen and s.  r.  sharpe , phys .
d * 86 * , 016007 ( 2012 ) [ arxiv:1204.0826 [ hep - lat ] ] .
r.  a.  briceno and z.  davoudi , phys .
d * 87 * , 094507 ( 2013 ) [ arxiv:1212.3398 [ hep - lat ] ] .
p.  guo , j.  dudek , r.  edwards and a.  p.  szczepaniak , phys .
d * 88 * , 014501 ( 2013 ) [ arxiv:1211.0929 [ hep - lat ] ] .
r.  a.  briceno and z.  davoudi , arxiv:1204.1110 [ hep - lat ] .
k.  polejaeva and a.  rusetsky , eur .
phys .  j.  a * 48 * , 67 ( 2012 ) [ arxiv:1203.1241 [ hep - lat ] ] .