Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: homework is a key element in every undergraduate - level physics course for practicing problem solving . students commonly have difficulties combining physics concepts with mathematical calculations and selecting among known mathematical tools . therefore , it is important to explore how students employ mathematics when approaching physics homework problems .    in this paper , we focus on the mathematical tools required in upper - division physics courses , such as evaluating integrals , using differential equations , or choosing appropriate coordinate systems . although these techniques are covered in mathematics classes , where students can successfully solve problems , many of them still struggle to apply these mathematical concepts to problems in physics . this is especially true at the upper - division when problems are more complicated and mathematics more tightly entwined with problem solving . in our study , we use the acer ( _ activation  construction  execution  reflection _ ) framework combined with network analytic methods to investigate how upper - division students written work shows their mathematical thinking on typical homework problems @xcite . the traditional undergraduate classical mechanics class is a textbook based 4 credit hour course with a solid foundation in the basics of theoretical physics . it is usually taken as the fourth physics course for physics majors and minors in the spring of their second year . typically , the students are concurrently enrolled in a differential equations course . the mechanics course under study is a typical exmaple of this course in the us . the textbook used in this course is _ classical mechanics _ by j. r. taylor @xcite . the lectures normally cover the first 10 chapters , including topics such as newton s laws of motion , momentum , angular momentum , energy , oscillations , lagrange s equations , two - body central - force problems and rotational motion of rigid bodies . we collected our data in a traditionally - taught classical mechanics class at a large land grant university . all homework problems were selected from the textbook and were chosen to encourage students to solve more practice problems . in return , students could earn extra credit towards their final grade . students were encouraged to work on problems in groups , but , for pedagogical reasons , they were to write their solutions independently .    in spring 2015 , twelve students were enrolled in the class . they had thirteen weekly homework assignments ( each including 10 - 15 problems ) . a complete solution to each homework problem was expected to be about half a page long and include diagrams , mathematics and verbal statements . we scanned students submitted homework for analysis before passing them to the grader . all students in the course consented to participate in our research study . several attempts to model physics problem solving break it into discrete steps to be followed in a linear order @xcite . however , unless students are explicitly taught to follow these algorithms , their approaches to problem solving rarely fit these linear models . in contrast , studies using knowledge - in - pieces frameworks to understand students non - linear problem solving `` in the wild '' ( e.g. , ref . @xcite ) rely on the richness of interview or video data to make inferences about student reasoning . the acer framework bridges the gap between prescriptive problem solving and knowledge - in - pieces @xcite . it was developed to analyze the use of mathematics in upper - division physics courses . the acer framework is organized around four components : _ activation _ ( determining the proper mathematical tool ) ; _ construction _ ( making a mathematical model ) ; _ execution _ ( performing mathematical steps ) ; and _ reflection _ ( checking the final solution ) . within each of these four components , students may perform multiple specific actions that are then detailed using sub - codes . for example , within construction , a student might choose a coordinate system or set the limits of integration . in solving a particular problem , a student may freely move among these components in any order , visiting each of them several times with different sub - codes . naturally , the details of the sub - codes depend on the specific physics topic and problems at hand . the acer framework was developed using think - aloud interview data on conceptually - rich research - based problems . these problems are different in character from the more traditional problems in our study . our data are wholly written accounts with sparse reasoning evidence . within these constraints , we refine our research question : how , if at all , can the modified acer framework be used to analyze written evidence of typical problem solving in a traditional course ?    prior work using acer framework has focused mostly on electricity and magnetism . we extended the standard sub - codes set to include codes appropriate for mechanics , as well as to cover students errors ( e.g. , mistakes done when taking derivatives , using an incorrect strategy , etc . ) . we coded students problem solutions , iteratively seeking new content codes and combining proposed codes until our code book covered @xmath0 of student work . as the code book approached stability eight people participated in inter - rater reliability testing to assure the code book was a valid representation of student problem solving . when the code book was stable , multiple graders looked at different problems to confirm the reliability of coding . once two graders coded individual problems and achieved @xmath1 agreement on student work , the code book was established . in this paper , we report on 3 problems selected from homework assignments across all students . we compare both across students and across problems to find patterns of mathematical tools use within each problem . to compare patterns in students solutions , we used techniques borrowed from network analysis @xcite . a network graph is a collection of points , called nodes , and a collection of lines connecting these points , called edges . when the direction of a connection matters , we have a _ directed graph_. in such case an edge from node a to b ( @xmath2 ) means something different than an edge from b to a ( @xmath3 ) . if the direction of a connection does not contain any information about the relation between nodes ( i.e , @xmath4 has the same meaning as @xmath5 ) , a graph is called _ undirected_. a graph can be characterized using various measures , such as node and edge count , diameter or density . node and edge count correspond to the number of nodes and edges , respectively . diameter is the longest graph distance between any two nodes in the network and it represents the linear size of a network . density , on the other hand , is the ratio of the number of existing edges to the number of all possible edges , and can be thought of as a measure of network effectiveness .    to determine the most important nodes in a network , one can perform centrality analysis . depending on the nature of a network , as well as the category of the `` importance of a node '' one is interested in , there are multiple centrality measures that can be used . the most commonly used centralities are degree , betweenness , closeness and eigenvector . the degree indicates the importance of a node by the number of nodes connected to it , where the larger the degree , the more important the node is . betweenness counts the number of shortest paths that pass through a given node . the closeness centrality is the inverse of the sum of distances to all other nodes . finally , the eignevector , is the sum of a node s connections to other nodes weighted by their degrees @xcite .    within this framework , in our case nodes represent acer sub - codes and the directed edges represent the order of steps students took to solve a given problem . in particular , we look at the diameter and density of each graph and compare it with these values for an expert solution . we also analyze the degree and eigenvector centralities as the most appropriate for studying our network . when looking at the network created from all solutions to a given problem , we weight the edges of the graph by how many times students connected a given two nodes to emphasize which connections between nodes are important . edges with higher weights represent crucial connections within a network . * problem 1 * a particle of mass @xmath6 is moving on a frictionless horizontal table and is attached to a massless string , whose other end passes through a hole in the table , where i m holding it . initially the particle is moving in a circle of radius @xmath7 with angular velocity @xmath8 , but i now pull the string down through the hole until a length @xmath9 remains between the hole and the particle . ( a. ) what is the particle s angular velocity now ? ( b. ) assuming that i pull the string so slowly that we can approximate the particle s path by circle of slowly shrinking radius , calculate the work i did pulling the string . ( c. ) compare your answer to part ( b. ) with the particle s gain in kinetic energy . students in this class were not taught an explicit procedure to solve problems . rather , they were free to solve them any way that they could which gave us the opportunity to gain insight into their reasoning when left to their own devices . a solution for the sample problem is presented in . this particular problem was not presented with a diagram . we expected that in order to visualize the physical system , the students would start with drawing one . in our acer framework , that would be represented by codes c1 ( pick a coordinate system ) and c2 ( visualize the problem ) . this problem has three parts and for all of them we expected students to begin with a general form of an equation , which would be coded as a2 ( use a general form of an equation ) . in part ( b. ) , as they have to do an integration , we expected codes e1 ( do the integration ) , e10 ( evaluate the integration ) and c5 ( set the limits of the integration ) to occur . both ( b. ) and ( c. ) require substitutions from the answer obtained in part ( a. ) , therefore we expected also codes r4 ( does this answer fit in the next part ) and e3 ( make substitutions ) to be used . as students perform some mathematical operations here , it s obvious to expect codes related to doing algebraic steps as well as doing substitutions , which are e3 and e5 . when we consider the problem as a whole , we expected nodes a2 , c1 , c2 , c5 , e1 , e3 , e5 ( doing algebra ) , e10 , r4 to be dominating codes .    a solution to problem 1 with a response coded according to the modified acer framework . ] when we used the acer framework alone , we saw codes a2 , c1 , c3 [ use an equation ( specific to the particular problem ) ] , c5 , e1 , e3 , e5 , e10 , r4 as the mostly commonly used , dominant codes . surprisingly , despite the lack of a figure in the problem statement , not all the students started by visualizing the problem . also , they used both the general form of the equation and the problem related non - general form . when we applied network analysis ( na ) along with centrality measures to the acer codes , we found that some of the nodes in the network , identified as important by the acer framework , turned out to be not significant when looked at using the na approach . on the other hand , network analysis identified nodes as important which were not significant in the acer framework . network graphs representing two students solutions to the problem presented in ( a. , b. ) , all students responses to this problems merged into a single network ( c. ) and the expert s response to it ( d. ) . ]    shows the problem with a single student s solution coded according to the acer framework . in , the network graph for that student ( fig . 2a ) , one other student ( fig . 2b ) as well as for all students responses combined into single graph ( fig . 2c ) and an expert solution ( fig . 2d ) to this particular problem are presented . one can see that , while for most part students used procedures which resulted in the same set of acer codes , their solving methods differ ( see fig . 2a and fig . moreover , the `` all responses '' graph ( fig . 2c ) is much more dense than each individual solution network , suggesting that , even though it was a well structured homework assignment , students took many different approaches to solving it . presents a summary of the basic descriptives for all networks presented in . when we look at sizes as measured by the diameter , we see that students networks are smaller than the expert s network . this was true in general , as the average diameter for the students was ( @xmath10 @xmath11 @xmath12 ) . the average density of students solutions was ( @xmath13 @xmath11 @xmath14 ) with values for individual student s ranging from @xmath15 to @xmath16 . when we look at centrality measures for this problem , the node e3 has the highest degree , indicating that it most frequently connects other steps in the problem solution . this is evidence that the mathematical operation of substitution is a crucial element in the solving process for this particular problem . moreover , almost all the homework problems had multiple parts within the problem and often subsequent parts required a substitution found in a previous section . when we look at the network created by combining solutions to this problem from all students , we found that the top three nodes for all students are the same as the top nodes for both sample students ( though their order is different ) . they differ , however , from the expert s solution . in particular , the component c3 lost its significance in the `` all students '' network . @lcccc & student a & student b & all students & expert solution + diameter & 6 & 5 & 7 & 7 + density & 0.155 & 0.167 & 0.632 & 0.191 + & c3 and e3 & r4 and e3 & e3 & e3 + & e5 and r4 & c3 & r4 & c3 + & a2 & e5 & e5 & e5 +      comparison between network graphs clearly show that students use two different approaches to solve this problem : 1 ) starting with visualizing the problem followed by either general or non general form of equations ; and 2 ) using general or non general form equations to start the problem , omitting the visualization part . while this particular problem does not explicitly ask for graphical representation of the physical situation , drawing appropriate diagrams is a useful problem solving strategy @xcite . students also tend to skip some mathematical steps , even though they are crucial in the solution . it seems as if some students were trying to solve the problem , but did not have a good grasp of the concepts underlying it . at the same time the expert s solution presents the problem solving process that , for illustrative purposes , includes steps that students might find unnecessary , thus becoming somewhat longer than it needs to be . thus , it may not be accurate to compare network graph for expert s solution with individual student s network graphs . comparison between all students and the expert s solution , on the other hand , will give us an idea about what procedures students used but the expert did not and vice versa . when we compare fig . 2c and fig . 2d , we see some additional nodes in the network for all students together that do not appear in the expert s solution . other than one node in reflection on the final solution most of these nodes belongs to either activation or the construction of the solution . because the approach and the construction of a solution is somewhat person dependent , we could expect those differences and also the way they reflect on their answer could be different . we found that there is a dependence on the original problem statement reflected on the network graphs . in particular , we identified three different types of prompts that suggest to students how they should start the problem :    1 . the prompt directs students towards the physical system / a diagram as the problem statement asks for an explanation of the physical system . an example statement looks like `` start with a diagram / free body diagram '' . 2 .   problem statement directs students to think about what bits / equations should they should bring together to get equivalent expression along with the physical system . to do so students requires physics knowledge along with math steps . an example statement looks like `` show that / prove that ... with the help of a physical system '' . 3 .   problem statement is straightforward and asks them to perform certain mathematical operations . an example statement looks like `` take the derivative / solve the above position equation for acceleration '' .    in order to determine whether students take the same approach when working on problems that are structured similarly , we identified three problems that belong to the same category ( the second kind mentioned above ) and we studied networks created by combining solutions to these three problems . since we compare problems structured similarly , we expect the crucial components of students solutions to be similar across all students . allowing for some slight variations , the core of the solutions should remain unchanged . the problems we chose for our analysis include the one discussed in a previous section and thus we expected nodes a2 , c1 , c2 , c3 , c5 , e1 , e3 , e5 , e10 , r4 to be the dominating codes . @lcccc & student c & student d & student e & all students + diameter & 13 & 11 & 10 & 11 + density & 0.133 & 0.277 & 0.211 & 0.309 + & e5 & e3 & e3 & e5 + & a2 & e5 & c5 & e3 + & c3 & c5 and a2 & e5 & c5 and c3 +    table  [ tabl2 ] shows characteristics of networks created by combining three structurally similar problems for three individual students , as well as of a network representing solutions from all students in the class combined . it s clear that nodes a2 , c3 , c5 , e3 and e5 are most central for these networks . as we expected based on the structure similarly , students used general ( a2 ) or non - general ( c3 ) form equation to start the problem . nodes e3 and e5 again are the key players in the network , but since the problems include integrals , also the node c5 ( setting the integration limits ) becomes a prominent one . by comparing networks across different problems , we can see whether students are consistent in their approaches to solving problems with similar structure . it is clear that in all cases , node e5 was the most central node . also , the four most prominent nodes ( a2 , c3 , e3 , e5 ) are the same for all students and for the merged network , although their order differs slightly . this alone could suggest that students take similar approaches to solve all three problems and that these approaches are quite similar among students . network graphs created by combining students solutions to three problems from three individual students ( a.  c. ) and from all students responses merged into a single network ( d. ) . ]    however , when we look at the network representing student s solutions in we can see that the network presented in fig . 3a is significantly different from the other two networks ( fig . 3b and 3c ) . also the network characteristics in reveal that the network for student c has greater diameter and lower density ( less connected graph with more nodes ) , indicating that this student got sidetracked in his solution . this could not be seen by the acer framework alone , as the significant codes for this students were the same as for students who took more directed approach . here we compare two problems that falls into type one and type three as classified at the beginning of section [ sec : net_dep ] . problem 2 belongs to the third type and students are required to perform trivial math procedures . problem 3 , one the other hand , falls into the first type and required to start with a diagram . * problem 2 * the shortest path between two points on a curved surface , such as the surface of a sphere , is called a geodesic . to find a geodesic , one has first to set up an integral that gives the length of a path on the surface in question . this will always be similar to the integral @xmath17 but may be more complicated ( depending on the nature of the surface ) and may involve different coordinates than @xmath18 and _ y_. to illustrate this , use spherical polar coordinates ( _ r _ , _ @xmath19 _ , _ @xmath20 _ ) to show that the length of a path joining two points on a sphere of radius _ r _ is @xmath21    * problem 3 * consider the pendulum of figure  [ figure4 ] , suspended inside a railroad car that is being forced to accelerate with a constant acceleration @xmath22 . ( a. ) write down the lagrangian for the system and the equation of motion for the angle @xmath20 . use a trick similar to the one used in @xmath23t@xmath24 to write the combination of @xmath25 and @xmath26 as a multiple of @xmath27 ( b. ) find the equilibrium angle @xmath20 at which the pendulum can remain fixed ( relative to the car ) as the car accelerates . use the equation of motion to show that this equilibrium is stable . what is the frequency of small oscillation about this equilibrium position ? ( commenter : have some kind of header break here ? )    diagram for problem 3 . ] the size of a network is a measure of the number of nodes within a network . the network graph discussed in section [ sec : net_dep_sin_stud ] ( problem 1 ) has 16 nodes and 40 edges with 8 strongly connected nodes . the network graph for problem 2 has 14 nodes and 32 edges along with five strongly connected nodes and problem 3 has a network graph with 11 nodes and 31 edges with two strongly connected nodes . the network graph of problem 1 has a relatively large network size compared to other two problems . this may be because problem 1 has three sub sections and other two problems have fewer sub sections . ( a. ) : network graph for all responses to problem 2 . ( b. ) : network graph for all responses to problem 3 . ] network analysis shows that the most central nodes for the network graph ( c  all student responses ) for problem 1 are e3 , r4 , e5 and a2 with a network diameter of 5 . according to the problem statement , first students have to understand the physical system and then they could use relevant equations to solve the problem . the network graph clearly shows that the students start with c- codes ( construction ) , where the students do a visual representation and then use the a- codes ( activation ) , where students use appropriate equations . because the problem has three sub sections , students might recall the answer from the previous part and substitute it in next part thus , we get r4 and e3 as most central nodes . the network graph ( a  all students responses ) for problem 2 have e5 , e3 , c5 and e2 as most central nodes along with network diameter 5 . the phrasing of problem 2 also requires visual representation and the students start with c- codes ( construction ) with a diagram . because we do not have subsections , r4 is not a most central node . the problem 3 ( b  all students responses in problem 3 ) network graph has most central nodes e3 , e5 , e2 and a2 with a network diameter of 4 . clearly some different nodes appear in this network graph compared to previous two problems . e2 becomes a central node because students were asked to find the equation of motion and students take derivatives . also to show that this equilibrium is stable students use boundary conditions by the c- codes c10 ( state boundary conditions ) and c11 ( apply boundary conditions ) . comparison between problem statements and the network graphs clearly show that there is a dependence of problem statement on the shape of the network . the analysis suggests that problems that consist of more sub sections lead to relatively complex network . the network graphs by all students and by each individual student can be directly used to estimate the students understanding of the subject material . more importantly the individual network graphs could be used to help or guide the students who might have stuck at a certain step with in the problem solving process . also , network graphs from both individual and all students together can be used to understand the difficulty level of the problems . if the students are taking more steps to reach the final answer than expected , there might be a issue with the way the problem statement was written . thus , the information about students performance revealed by the network graphs could be used to enhance the quality of the homework or quiz problems . the acer framework to code solutions lays out what steps students take in their solutions . the network analysis gives the most important nodes within the network . using both acer framework and network analysis , we can see what are the most important steps within the solution . network analysis indicates that in most of the analyzed cases , the most important nodes are e3 , e5 , a2 , r4 , and therefore substitution and algebra in general are crucial tools for solving these problems . most of the homework problems had a few parts within the problem and in these cases students may have had to go back and check the previous solution and then substitute it in the next part . while many problems have similar central nodes , the exact nodes and their order are different across students and across problems . the shape of the network graphs varies with different problem statements . we can conclude , even though these homework problems seem to be well designed , different students have different approaches to start the problem and different patterns in their construction section , in most of the cases , though , the math execution steps are the same . further research is needed to see exactly what properties of problem statements prompt the students to use specific mathematical tools . nonetheless , network analysis of students problem solutions reveals patterns in their tool use on typical problems . though our data cover 23 problems overall , they are limited to classical mechanics . to expand our study more broadly across the upper division , we are currently collecting data on an electricity and magnetism ( em ) course , and augmenting students written work with video recordings of them solving homework problems . 10 r.e .  pepper , s.v . chasteen , s.j.pollock , and k.k . perkins , observations on student difficulties with mathematics in upper - division electricity and magnetism , _ phys . * 8 , * 010111 ( 2012 ) . wilcox , m.d . caballero , r.e . pepper , and s.j . pollock , upper - division student understanding of coulomb s law : difficulties with continuous charge distributions , _ proceedings of the physics education research conference 2012 _ p. 1283 , ( 2012 ) . wilcox , m.d . caballero , d.a . rehn , and s.j . pollock , analytic framework for students use of mathematics in upper - division physics , _ phys . _ * 9 * , 020119 ( 2013 ) . kolaczyk and g.  csrdi , statistical analysis of network data with r ( springer , 2014 ) . taylor , classical mechanics ( university science books , 2005 ) . p.  heller , r.  keith , and s.  anderson , teaching problem solving through cooperative grouping . part 1 : group vs. individual problem solving , am . 60 , 627 ( 1992 ) . e.  redish , problem solving and the use of math in physics courses , arxiv : physics/0608268 ( 2005 ) . sayre and m.c . wittmann , plasticity of intermediate mechanics students coordinate system choice , _ phys . res . _ * 4 * , 020105 ( 2008 ) . wallace and s.v . chasteen , upper - division students difculties with ampre s law , _ phys . res . _ * 6 * , 020115 ( 2010 ) . j.  larkin and h.  simon , why a diagram is ( sometimes ) worth ten thousand words , _ cognitive science _ , * 11 * , 65 - 69 ( 1987 ) . j. zhang and d. norman , representations in distributed cognitive tasks , _ cognitive science _ , * 18 * , 87 - 122 ( 1994 ) . thanks to deepa chari for her important contribution to this paper . also we thank participating teachers and students . this work funded by nsf award 1430967 . jpz is supported by nsf phy 1344247 .
Abstract:

Answer: the mathematization project investigates students use of mathematical tools across the undergraduate physics curriculum . as a part of this project , we look at intermediate mechanics students written homework solutions to understand how they use those tools in approaching traditional mechanics problems . we use a modified version of the acer framework to analyze students solutions and to identify patterns of mathematical skills used on traditional problems . we apply techniques borrowed from network analysis and the resources framework to build a `` fingerprint '' of students mathematical tool use . in this paper , we present preliminary findings on patterns that we identified in students problem solving .    _ keywords _ : mathematization , undergraduate homework , problem solving , network analysis


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: one - dimensional quantum systems have always attracted a lot of theoretical interest : their exact solutions give , in fact , useful insights on the role of the interactions and other non - perturbative features @xcite . in this respect , the growing ability to realize and manipulate one - dimensional bose gases @xcite has provided an highly controllable experimental counterpart to these theoretical achievements . quasi - one - dimensional bose gases are obtained by using a cigar - shaped external trapping potential , elongated in a direction , with the other degrees of freedom frozen due to the presence of a tight transverse confinement . in the experiments , several variants of the interacting bose gas in one dimension can be implemented : an optical lattice can be added to detect the mott - superfluid transition in one dimension @xcite , the effective one - dimensional interaction can be tuned @xcite to observe a tonks - girardeau gas of ultracold atoms @xcite , or the effect of the temperature can be studied @xcite . important tools that permit to further control the properties of low - dimensional bose systems are the tuning of an external magnetic field near a feshbach resonance @xcite , and , in perspective , the implementation of the recently proposed schemes to engineer effective three - body interactions @xcite .    the technique of feshbach resonances permits to change the sign of the scattering length : by switching from repulsion to attraction , i.e. from positive to negative scattering length , the homogeneous @xmath2-d gross - pitaevskii equation ( gpe ) admits a solution corresponding to a localized wave - function , the so - called _ bright soliton _ bright matter - wave solitons were created both in bose - einstein condensates of @xmath3li @xcite and @xmath4rb atoms @xcite . various localized states has been also produced in quasi - one - dimensional geometries ( for reviews see the book @xcite ) .    with attractive two - body interactions , a crucial role is played both by the dimension of the system and the trapping potential . in three dimensions , for instance , homogeneous attractive bosons are unstable against the collapse , but the presence of an external harmonic trap can stabilize them : the critical value of the interaction coupling that gives rise to the collapse can be obtained from the gpe @xcite , and the critical particle number is given by @xmath5 where @xmath6 is the total number of particles , @xmath7 is the scattering length and @xmath8 is the harmonic oscillator length @xcite . in the one - dimensional case , the bright soliton solution is the ground state of the homogeneous gpe with negative scattering length . furthermore , the gpe ground state energy is in agreement , in the thermodynamic limit , with the ground state energy obtained by bethe ansatz for the attractive one - dimensional bose gas @xcite ( see more in section ii ) .    in this paper , motivated by the recent papers @xcite in which different schemes have been proposed to realize effective tunable three - body interactions , we consider an attractive three - body contact potential and , more generally , a @xmath0-body contact interaction . we consider the limit of large number of particles , @xmath9 , with the constraint @xmath10 ( @xmath11 being the strength of the @xmath0-body interactions ) so that the energy per particle is finite . since no bethe solution is available in the general case of @xmath0-body interaction , we employ an hartree approximation to study the problem in the limit mentioned above . this means that the ground state energy is estimated by using the bright soliton solution of a generalized mean - field gpe equation . as we will show , the @xmath1 is a special case : for this value , in fact , a localized soliton wavefunction exists only for a critical value of the interaction strength and has an infinite degeneracy . the stabilization of this bound state can be cured by putting the system in an external harmonic trap . the variational approach , that we will also employ , reveals the tendency of the higher body interactions to become more unstable in higher dimensions . it is worth stressing that the case we are considering does not consist of a @xmath0-body interaction added to the @xmath12-body interaction of the bose gas : we are interested , in fact , to the effect of the @xmath0-body in its own , since the coefficient of the two - body interaction can be tuned to be zero @xcite . the plan of the paper is the following : in section ii we introduce the hamiltonian corresponding to @xmath0-body contact attractive interactions and we write the ( mean - field ) generalized gpe . the familiar case @xmath13 is briefly recalled . in section iii the bright soliton solution for the homogeneous limit is obtained by using a mechanical analogy with a fictitious particle moving in a potential , and its properties are investigated . the comparison with the numerical results confirms that for @xmath14 this is the ground state of the generalized gpe , as expected . the ground state energy by varying @xmath0 is also determined . in section iv we consider the effect of an harmonic trap : using a variational ansatz for the ground state we determine the critical value of the interaction needed to stabilize the bound state . in section v there are our conclusions . the general quantum hamiltonian for an homogeneous one - dimensional bose gas with @xmath0-body interactions @xmath15 is @xmath16 where @xmath17 is the bosonic field operator . the lieb - liniger hamiltonian for the interacting one - dimensional bose gas @xcite has the kinetic term plus a density - density term involving pairs of particles interacting via a contact two - body potential ; this corresponds to @xmath13 and @xmath18 : @xmath19 positive ( negative ) corresponds to repulsion ( attraction ) between the bosons . the low - energy properties of the lieb - liniger model can be studied by the luttinger liquid effective description @xcite obtained by bosonization @xcite ( a general discussion of the correlation functions is presented in @xcite ) .    for @xmath0-body attractive contact interactions we set @xmath20 ( @xmath21 ) : the hamiltonian ( [ ham - gen ] ) reads then @xmath22 in the heisenberg representation , the equation of motion for the field operator is given by @xmath23= -\frac{\hbar^2}{2 m } \frac{\partial^2}{\partial x^2 } \hat{\psi } -   c \left ( \hat{\psi}^{\dag } \right)^{n-1 } \left ( \hat{\psi } \right)^{n-1 } \hat{\psi}\,\,\ , . \label{dyn}\ ] ] for @xmath13 , the corresponding lieb - liniger model is integrable and the ground state energy @xmath24 can be determined by bethe ansatz @xcite : the final result is given by @xmath25 where @xmath6 is the total number of particles . for large @xmath6 , from ( [ gs - bethe ] ) it follows that one has to keep the product @xmath26 in order to have a finite ground state energy per particle . using the integrability of the @xmath13 model , the correlation functions of the attractive one - dimensional bose gas at zero temperature were recently calculated in @xcite . for the three - body problem ( @xmath1 ) , no bethe ansatz solution is available , except for a more complicate double-@xmath27 function bose gas which can be mapped in a one - dimensional anyon gas @xcite . hence , to estimate the ground state energy @xmath24 we propose here to employ a mean - field ( hartree ) approach : in this approach , the ground state energy is given in terms of the ground state energy of a generalized gpe . the same procedure will be employed for other values of @xmath0 . before we start the discussion of the general @xmath0-body case , let s briefly remind how this task can be successfully done for @xmath13 @xcite . first of all , in the mean - field approximation the ground - state wavefunction is written as @xmath28 where the function @xmath29 is the ground state of the time - independent homogeneous gpe , i.e. the nonlinear schrdinger equation ( nlse ) , given by @xmath30 where @xmath31 is the chemical potential and the normalization is given by @xmath32 . the energy is expressed as @xmath33 \psi_0(x)\,\,\ , . \label{energy - gpe}\ ] ] the static bright soliton solution of ( [ gpe ] ) is given by @xmath34 with @xmath35 and @xmath36 . substituting this expression in ( [ energy - gpe ] ) one gets @xmath37 i.e. , the exact result ( [ gs - bethe ] ) apart terms @xmath38 . a comment is in order : in the homogeneous one - dimensional interacting case there is , strictly speaking , no condensate . however the condition @xmath39 implies that , for large @xmath6 , the coupling constant should scale to zero , @xmath40 : hence , we are in a weak coupling regime where the mean - field gpe is expected to give reasonable results . in a similar way , for @xmath41 ( repulsive interaction ) the comparison between the exact and the gpe ground state energy shows that the latter gives the correct behaviour for @xmath40 while the bogoliubov approximation gives the exact first - order corrections for small @xmath42 @xcite . based on the analysis above , for general @xmath0 and in the limit @xmath40 we expect that a reasonable description of both the ground state properties and the low - energy dynamics is given by the mean - field generalized homogeneous gpe @xmath43 where the nonlinearity degree @xmath44 of the generalized nonlinear schrdinger ( gnlse ) equation is related to @xmath0 by @xmath45 the mean - field ground state is given by the time - independent gnlse equation @xmath46 where , as before , @xmath31 is chemical potential and @xmath47 is normalized to the total number of particles @xmath6 , i.e. @xmath48 .    with @xmath49 integer , @xmath44 is an even integer ; however , in eqn . ( [ mu ] ) @xmath44 can take any real positive value and , in the following , we will consider this general case . in this respect let s comment that the axial dynamics of a bose - einstein condensate induced by an external potential with cylindrical symmetry in the transverse directions can be studied by introducing an effective one - dimensional gpe equation with @xmath50 @xcite and that for bose - einstein condensates in one - dimensional optical lattices the effective equation has a value of @xmath44 that depends on the details of the trapping potentials and it is , in general , a non - integer value @xcite . in the following we will study the attractive @xmath0-body problem in the thermodynamic limit , defined by @xmath51 , with the product @xmath52 kept fixed . this will ensure the energy per particle of the gnlse bright soliton to be finite . in dimensionless units , rescaling the wave function @xmath53 , eqn . ( [ mu ] ) reads @xmath54 where @xmath55 and @xmath56 are the dimensionless versions of @xmath57 and @xmath31 respectively . for the ground state of this equation we look for a real solution , with the normalization condition @xmath58 obviously , once a static solution @xmath29 of eqn . ( [ mu - ad ] ) has been found , the corresponding soliton wave solution with velocity @xmath59 is given by @xmath60 the solution of eqn . ( [ mu - ad ] ) can be found by using a mechanical analogy . in fact , interpreting @xmath61 as the time variable and @xmath29 as the coordinate of a fictitious particle , eqn . ( [ mu - ad ] ) formally corresponds to the newton s equation of motion of this particle ( of mass @xmath62 ) , subjected to the force @xmath63 this force can be derived by the potential @xmath64 as any motion of a particle in a potential , this is accompanied by the integral of motion that corresponds to its mechanical energy @xmath65     for negative values of @xmath56 . @xmath66 are the inversion points of the motion . ]    following this mechanical analogy , it is easy to see that a non - trivial motion can take place only if @xmath67 , where the typical shape of the potential is similar to the one drawn in fig . [ potentialfigure ] . notice that a solution is always given by the equilibrium configuration of the potential ( [ potential ] ) , i.e. @xmath68 this solution can be normalized only on a finite volume @xmath69 , with the dependence of the chemical potential on the volume determined by the normalization condition ( [ important ] ) , i.e. @xmath70 : given this dependence of the chemical potential , the constant solution is simply @xmath71 to determine the ground state , we have to compare the gpe energy of this constant solution with the one of a localized wavefunction . for this solution , both @xmath29 and @xmath72 should vanish when @xmath73 . this condition fixes the constant value of @xmath74 to be zero ( notice that this value is @xmath75 the gpe energy ) . in this case , the fictitious particle takes off from the origin at @xmath76 , moving to the right ( or , equivalently to the left , since the original equation is invariant under @xmath77 ) , until it reaches the inversion point @xmath78 at the time @xmath79 . once the particle arrives in @xmath78 , it inverts its motion and comes back to the origin with a vanishing velocity . it is clear from this analogy that @xmath78 will be the maximum of the bright soliton solution . what we said , however , is not the end of the story . in fact , the kind of motion we have just described occurs for any potential with the shape shown in fig . [ potentialfigure ] . but we are looking for that particular motion that satisfies the additional constraint ( [ important ] ) and this condition can be fulfilled only for a particular shape of the potential , i.e. for a particular combination of the parameters @xmath56 and @xmath55 : it is as if the solution is looking for its proper potential .    in the following , it is convenient to introduce the quantities @xmath80 using the integral of motion @xmath74 , for generic values of @xmath56 , @xmath55 and @xmath44 ( with @xmath67 , @xmath81 and @xmath82 ) the solution is given by a quadrature @xmath83 where @xmath84 is the inversion point reached by the particle at the `` time '' @xmath85 . using the exact expression of the integral of the left hand side @xmath86 \,\,\ , ,   \label{exactintegral}\ ] ] one gets @xmath87 it remains now to impose the normalization eqn . ( [ important ] ) to the solution ( [ q - result ] ) : this fixes the shape of the potential , i.e. the relation between @xmath56 and @xmath55 @xmath88 when @xmath89 , we can use this equation to express @xmath56 as a function of @xmath55 and , in particular , to write the normalization @xmath78 as @xmath90 in fig . [ gss ] we plot the soliton solution ( [ q - result ] ) for different values of @xmath44 at @xmath91 .    ) for @xmath92 ( dimensionless units are used with @xmath91 ) . ]    however , when @xmath93 , eqn . ( [ relation - mu - g ] ) leaves @xmath56 undetermined : this means that the corresponding wavefunction @xmath94 is solution of the nonlinear schrdinger equation ( [ mu - ad ] ) for _ every _ @xmath56 . in this case , however , only a particular value of @xmath55 , given by @xmath95 guarantees its correct normalization ( [ important ] ) . expressed in more physical terms , the attractive @xmath96-body interaction has the peculiarity that one can arbitrarily vary the chemical potential provided that the coupling constant be fine - tuned to the critical value @xmath97 : increasing or decreasing ( in modulus ) the chemical potential simply results , in this case , in shrinking or enlarging the shape of the soliton . this is shown in fig . [ mu - pot ] where the wavefunction ( [ psialpha4 ] ) is plotted for two different values of @xmath56 : in the inset we plot the corresponding potential ( [ potential ] ) , showing a larger ( smaller ) inversion point corresponding to the smaller ( larger ) width . the fact that @xmath56 is undetermined and the soliton can change arbitrarily its shape does not imply that the gpe energy is undetermined : in fact , the explicit computation of the next subsection shows that , in this case , the energy does not depend on the value of @xmath56 . then , for @xmath1 and @xmath98 , an infinite degeneracy parametrized by the chemical potential @xmath99 occurs .    in fig . [ mu - en](a)-(b ) we plot the chemical potential @xmath56 for two different values of @xmath55 , one smaller than @xmath97 and the other larger : it is seen that for @xmath100 ( @xmath101 ) , then @xmath102 ( @xmath103 ) for @xmath104 while @xmath105 ( @xmath102 ) for @xmath106 . the singular nature of the @xmath96-body interaction can then be recovered by studying the limit @xmath107 of the formulas ( [ q - result ] ) , ( [ relation - mu - g ] ) , ( [ enne ] ) given above : for @xmath108 , if @xmath109 the normalization @xmath78 goes to 1 , while if @xmath110 , @xmath111 and @xmath112 ( i.e. we have a non - localized solution ) whereas if @xmath113 , both @xmath78 and @xmath56 diverge , i.e. the wavefunction collapses to the origin . it is worth to mention that a singular behavior of the nonlinear schrdinger equation relative to the value @xmath114 has been also observed in the dynamical blowing up of the moving wave - packets of this equation : the interested reader is referred to the mathematical literature for a detailed discussion of this issue @xcite . in the present application , this instability means that the local @xmath96-body local attractive interactions can not sustain a bound state unless there is a fine tuning of the interaction . in the next section we will show how an external trap can help to stabilize the bound state for a generic value of the coupling .    ) for @xmath1 and @xmath98 plotted for @xmath115 ( solid line ) and @xmath116 ( dot - dashed line ) . inset : corresponding potential ( [ potential ] ) for @xmath115 ( solid line ) and @xmath116 ( dot - dashed line ) . ]    ) for @xmath117 ( solid line ) and @xmath118 ( dashed line ) for @xmath119 - ( b)-(d ) chemical potential and energy for @xmath120 ( solid line ) and @xmath121 ( dashed line ) for @xmath122 . ]    to understand better the behaviour of the solution @xmath29 as a function of @xmath44 , let s define the width @xmath123 as @xmath124 . one gets @xmath125 where @xmath126 . one finds @xmath127 and @xmath128 . for large @xmath44 one has @xmath129 , while , of course , @xmath130 for @xmath131 ( no localized soliton without interaction ) . for @xmath100 , from eqn . ( [ sigma ] ) one sees that for @xmath104 , @xmath132 , while for @xmath106 , @xmath133 . in fig . [ disp ] we plot @xmath134 for @xmath135 from eqn . ( [ sigma ] ) and , for completeness , also the widths for some value of @xmath44 obtained from the numerical gnlse . a divergence is observed for @xmath104 , corresponding to the @xmath96-body attraction : the bright soliton becomes larger and larger getting close to @xmath93 , while for @xmath44 slightly larger than @xmath136 the soliton becomes extremely narrow . this means that there is a collapse of the solution ( [ q - result ] ) going to @xmath93 from large values of @xmath44 . at variance , for @xmath101 , then for @xmath104 , @xmath133 , while for @xmath106 , @xmath132 . it should be stressed that , for @xmath122 , although ( [ q - result ] ) is a solution of the gnlse ( [ mu - ad ] ) , it is no longer its ground state : the divergence of @xmath123 for @xmath104 is signaling the disappearance of the bound state due to the @xmath96-body interaction .     versus the nonlinearity degree @xmath44 . solid line : eqn . ( [ sigma ] ) ; open circles : @xmath137 from the numerical determination of the ground state of the gnlse ( [ mu - ad ] ) . the value @xmath135 and dimensionless units are used . ]      using the bright soliton solution ( [ q - result ] ) we can now estimate the energy per particle . going back to the physical dimensions of all quantities and normalizing now @xmath47 to @xmath6 , for @xmath89 the chemical potential is given by @xmath138 where @xmath139 the energy per particle is then obtained from the gnlse energy functional @xmath140 \psi_0(x)\,\,\ , . \label{energy - ggpe}\ ] ] using eqn . ( [ q - result ] ) we obtain @xmath141 where @xmath142\,\,\ , . \label{en - fun}\ ] ] for @xmath13 , it is @xmath143 and the previous energy ( [ gs - gpe ] ) is recovered . from eqn . ( [ en - dim ] ) it follows that in order to maintain finite the energy per particle for large @xmath6 one has to keep @xmath57 fixed . by a numerical determination of the ground state of the gnlse , we have verified that ( [ q - result ] ) indeed coincides with the ground state for @xmath144 _ both for @xmath100 and @xmath101_. in fig . [ ener ] we compare for @xmath135 the ground state energy per particle from eqn . ( [ en - dim ] ) with the ground state energy obtained for some values of @xmath44 obtained from the numerical gnlse . for @xmath101 and @xmath119 , a similar agreement is obtained .    for @xmath93 , as discussed in the previous section , the chemical potential is undetermined . however , a direct substitution of ( [ psialpha4 ] ) in ( [ energy - ggpe ] ) reveals that @xmath145 for @xmath98 . since ( [ psialpha4 ] ) is a solution of the gnlse ( [ mu - ad ] ) for arbitrary @xmath146 , and then with arbitrary width , we conclude that an infinite degeneracy - parametrized by a negative chemical potential - occurs .    using the energy ( [ en - dim ] ) we can also estimate the energy of the constant solution in the finite interval @xmath147 $ ] : it is @xmath148 , where @xmath149 is the density . to compare this energy with the previous result ( [ en - dim ] ) for the bright soliton solution we have to choose how to perform the thermodynamic limit : if we choose to keep fixed the quantity @xmath150 , with large but finite value of @xmath6 , sending @xmath69 to infinite the energy of the constant solution vanishes . this means that for @xmath122 the constant solution is the ground state of the system ; @xmath93 and @xmath98 is the case in which both the solutions have zero energy .    . solid line ( dashed line ) : eqn . ( [ energy - ggpe ] ) for @xmath119 ( @xmath122 ) ; open circles : energy of the numerical ground state of the gnlse ( [ mu - ad ] ) . dimensionless units ( with @xmath151 ) and @xmath135 , as well as different scales of the energy for @xmath119 ( left part ) and @xmath122 ( right ) , are used . ] in three dimensions , the presence of an external harmonic trap ( with frequency @xmath152 ) can help to stabilize the attractive two - body interaction : the critical value @xmath153 of the particle number that induces the collapse , obtained from the gpe , is given by @xmath154 @xcite . this critical value can be estimated by a variational method @xcite : a gaussian trial wavefunction , with the width variational parameter , is introduced and the corresponding energy computed @xcite . without external trap ( @xmath155 ) , the energy does not have a minimum . however , with @xmath156 , two situations are possible : for a number of particle @xmath157 a metastable minimum appears , while for @xmath158 there is no minimum .    in this section we consider the corresponding problem in one dimension with a @xmath96-body interaction and we show that there is a critical value @xmath159 of the interaction , such that for @xmath160 , the bound state is stable . to this aim , we use the variational wavefunction @xmath161 normalized to @xmath6 . the energy to be minimized is obtained by inserting the variational wavefunction ( [ var ] ) in the generalized gpe functional @xmath162\psi(x ) \,\,\ , . \label{e - functional}\ ] ]    to better illustrate the peculiarities of the one - dimensional case , it is useful to perform the analysis also in higher dimensions . the energy in @xmath163 is given by @xmath164 where @xmath165    let us consider initially the homogeneous case : for @xmath166 the energy ( [ effe ] ) has a minimum only when @xmath167 . for @xmath168 , the critical value corresponds to @xmath93 : this is in agreement with the result of the previous section , which is now obtained by a variational approach . the critical condition @xmath169 is plotted in fig.[crit - fig ] : this figure shows that the higher @xmath0-body interactions tend to be more unstable in higher dimensions .    ): @xmath163 corresponds respectively to @xmath170 , i.e. @xmath171 . ]    when the harmonic trap is present ( @xmath156 ) , there is still a minimum when the values @xmath44 is pushed to the range of value @xmath172 . when @xmath173 , we can identify the critical value @xmath159 as indicated in @xcite for the @xmath174 and @xmath175 case : since @xmath176 for @xmath177 and @xmath178 for @xmath179 , the critical value is obtained by the conditions @xmath180 . in this way we arrive to the result @xmath181 the instability curve ( [ g_critical ] ) depends on @xmath182 : with dimensionless units ( and @xmath151 ) , in one dimension as @xmath183 the critical value @xmath97 goes to zero for @xmath184 and to infinity otherwise ; in two dimensions the behaviour is similar except that @xmath185 when @xmath186 ; while in three dimensions critical value goes to infinity for @xmath187 and to zero otherwise . a plot of the @xmath188 in @xmath168 for @xmath122 is presented in fig . [ plot - c - crit ] . let s now examine the critical point @xmath189 , which in @xmath168 corresponds to the @xmath96-body interactions . minimizing the energy ( [ energy ] ) with respect to @xmath190 , one finds that there is the critical value @xmath191 for @xmath192 there is a minimum and the system is stable , while for @xmath159 the energy does not have ever a minimum : hence , irrespectively of how large @xmath152 may be , the system always collapses . for @xmath155 ( no trap ) and @xmath193 ( @xmath194 ) , the minimum value of the energy is then obtained for @xmath195 ( @xmath196 ) . notice that the variational approach gives for @xmath168 and @xmath114 , the critical value ( in dimensionless units ) @xmath197 , in good agreement with the analytical value @xmath198 .     vs. @xmath44 for @xmath122 , with @xmath168 . dimensionless units ( with @xmath199 ) are used , with @xmath200 ( dashed line ) and @xmath201 ( solid line ) . ] in this paper we have analyzed the one - dimensional bose gases with @xmath0-body local attractive interactions : by using a mean - field approach , we found that @xmath1 ( i.e. , @xmath93 ) is a critical point , and that the localized solution is possible only for a critical value @xmath202 of the interaction strength . for this critical value , an infinite degeneracy occurs : this degeneracy is parametrized by the chemical potential ( i.e. , eigenfunctions with the same negative @xmath31 has the same energy ) . for @xmath119 , the bright soliton coincides with the ground state of the gross - pitaevskii equation : when @xmath203 ( @xmath194 ) , then @xmath104 gives a width going to diverge ( vanish ) . we have also studied how an harmonic trap can make stable this bound state , pointing out that for @xmath1 there exist a critical value of the @xmath96-body interaction strength , and below such critical value the localized state occurs . above this critical value , the collapse is not prevented even for very large trap frequency . a brief discussion of the role played by the dimension for @xmath0-body local attractive interactions has been also presented , showing that higher body interactions are more unstable in higher dimensions . we also mention that two - body nonlocal attractive interactions has been studied , showing different ranges of stability with respect to the local ones @xcite : we could then expect that for @xmath96- and @xmath0-body interactions this effect could become even more relevant . several proposals have recently addressed the issue of inducing and controlling three - body terms . in @xcite it has been proposed to use cold polar molecules driven by microwave fields to obtain strong three - body interactions , controllable in a separate way from the two - body interactions , which in turn can be switched off @xcite . three - body interactions can be effectively induced in mixtures of bosonic particles and molecules : in @xcite the ground state of rotating bose gases close to a feshbach resonance has been studied , showing that for suitable parameters they are fractional quantum hall states , whose excitations obey non - abelian exchange statistics . in @xcite it was shown that a system of atoms and molecules in a one - dimensional lattice can be effectively modeled by a three - body local ( i.e. , contact ) interaction , characterized by a strength @xmath204 and in the limit @xmath205 ( without a two - body interaction ) the ground state properties were investigated by a pfaffian - like ansatz . the strength @xmath204 of the three - body interaction can be made also negative by using feshbach resonances . one of the main reasons of interest of these proposals relies on the fact that exotic quantum phases , such as topological phases , appear to be ground states of hamiltonian with three or more body interaction terms , an example being the fractional quantum hall states described by the pfaffian wavefunctions @xcite . the excitations of pfaffian states are non - abelian anyons , on which schemes of fault - tolerant topological quantum computation are based @xcite . we think that , in perspective , the possibility to induce and tune effective @xmath96-body interactions could become an important tool to control the nonlinear dynamical properties of localized wave - packets and to induce new exotic strongly correlated phases in ultracold atoms .    _ acknowledgments_. we would like to thank b. dubrovin and p. calabrese for stimulating discussions . this work is partially supported by the esf grant instans and by the miur projects `` quantum field theory and statistical mechanics in low dimensions '' and `` quantum noise in mesoscopic systems '' .
Abstract:

Answer: we study the ground state properties of a one - dimensional bose gas with @xmath0-body attractive contact interactions . by using the explicit form of the bright soliton solution of a generalized nonlinear schrdinger equation , we compute the chemical potential and the ground state energy . for @xmath1 , a localized soliton wave - function exists only for a critical value of the interaction strength : in this case the ground state has an infinite degeneracy that can be parameterized by the chemical potential . the stabilization of the bright soliton solution by an external harmonic trap is also discussed , and a comparison with the effect of @xmath0-body attractive contact interactions in higher dimensions is presented .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the recent observation suggesting that the fine - structure constant was smaller in the past  @xcite than it is today has led to renewed interest in the idea of using time ( and space ) variation of fundamental parameters as a probe of high - energy physics . based on the principle that all that is not forbidden is mandatory , time variation in the electromagnetic force is perhaps not surprising . moreover , one would expect that the other forces of nature vary with time as well . this idea is quantifiable by assuming a suitable grand unified theory ( gut ) at a scale @xmath0 or a particular brane - world scenario  @xcite . the assumption that the short - distance couplings of the theory are related at all times leads to relations between the variation of the strength of the electromagnetic interaction , @xmath1 , and the variation of other standard model parameters , such as the light quark masses , @xmath2 . as a change in @xmath2 and @xmath1 will naively lead to a change in the positions of nuclear energy levels , special interest has been paid to the light - element abundances predicted by big bang nucleosynthesis ( bbn ) and also to the abundance of isotopes produced by the oklo `` natural reactor '' in the hope that these abundances can be used to constrain high - energy physics  @xcite .    in this work we will critically analyze the two - nucleon sector using an effective field theory ( eft ) that respects the approximate @xmath3 chiral symmetry of qcd and has consistent power - counting  @xcite . this theory allows for a systematic study , consistent with qcd , of the bound - state in the @xmath4 coupled - channels , the deuteron , and the scattering amplitude in the @xmath5 channel . this is the first study of the @xmath2-dependence of these quantities that is _ complete _ at next - to - leading order ( nlo ) in a consistent eft . unfortunately , due to the current lack of understanding of the low - energy behavior of qcd , only weak bounds , if any , can be placed on the time - dependence of fundamental couplings from the two - nucleon sector . we will explicitly demonstrate this by finding plausible values for couplings in the eft for which the deuteron binding energy varies relatively little , and for which the di - neutron system in the @xmath5 channel is unbound , over a relatively wide range of quark masses . the deuteron is an interesting object to study for a variety of reasons . from a fundamental point of view it is quite intriguing as it is bound by only @xmath6 , which is much smaller than the typical scale of strong interactions , @xmath7 . clearly , in order to arrive at such a small binding energy , fine - tunings are involved   and the @xmath8 coupled - channels are very close to an unstable ir fixed point of the renormalization group . at this fixed point , the scattering lengths are infinite  @xcite . ] and consequently naive dimensional analysis ( nda ) as applied to the single - nucleon sector or the mesonic sector , which we will call naive naive dimensional analysis ( n@xmath9da ) , on this system is doomed to fail . unfortunately , all previous attempts to extract bounds on the variations of fundamental couplings from the deuteron , and in general the two - nucleon sector , have implemented n@xmath9da . from a more phenomenological standpoint , the smallness of @xmath10 plays a key role in the synthesis of light elements in bbn . the impressive agreement between the predictions of bbn and observation suggests that new physics that would have significantly modified the deuteron at the time of bbn is absent . the @xmath5 channel is quite similar to the @xmath4 coupled - channels in one important way , its scattering length is unnaturally large , @xmath11 . while there is no bound state in this channel for the physical values of the light - quark masses , a small increase in the strength of the nuclear force would bind two nucleons in this channel . the existence of a bound state in this channel , e.g. a di - neutron , @xmath12 , in the nucleosynthesis epoch would be quite profound and could substantially modify the predictions of bbn . during the last decade there has been a significant effort to construct an eft to describe nuclear physics . while it is straightforward to write down all possible terms in the effective lagrange density for two or more nucleons , arriving at the correct power - counting proved to be a difficult task . weinberg s ( w ) original proposal  @xcite for an eft describing multi - nucleon systems was to determine the nucleon - nucleon ( nn ) potentials using the organizational principles of the well - established eft s describing the meson - sector and single - nucleon sector ( chiral perturbation theory , @xmath13pt ) , and then to insert these potentials into the schrdinger equation to solve for nn wavefunctions . observables are computed as matrix elements of operators between these wavefunctions . w power - counting has been extensively and successfully developed during the past decade to study processes in the few - nucleon systems . this method is intrinsically numerical and is similar in spirit to traditional nuclear - physics potential theory . unfortunately , there are formal inconsistencies in w power - counting  @xcite , in particular , divergences that arise at leading order ( lo ) in the chiral expansion can not be absorbed by the lo operators . problems persists at all orders in the chiral expansion , and the correspondence between divergences and counterterms appears to be lost , leading to uncontrolled errors in the predictions for observables . this formal issue was partially resolved by kaplan , savage and wise ( ksw ) who introduced a power - counting in which pions are treated perturbatively  @xcite . the nn phase - shifts and mixing angle in the @xmath5 and @xmath4 coupled - channels have been computed to next - to - next - to - leading order ( n@xmath9lo ) in the ksw expansion by fleming , mehen and stewart ( fms )  @xcite from which it can be concluded that the ksw expansion converges slowly in the @xmath5 channel and does not converge in the @xmath4 coupled - channels . therefore , neither w or ksw power - counting provide a complete description of nuclear interactions ( for recent reviews see ref  @xcite ) . the problems with w and ksw power - counting appear to have been resolved in the work by bedaque , van kolck and the authors  @xcite , which from this point on we will refer to as bbsvk . it was realized in fms that the contributions to the amplitude that lead to non - convergence in the @xmath4 coupled - channels persist in the chiral limit ( it is the chiral limit of iterated one - pion - exchange ( ope ) that is troublesome ) . therefore , in bbsvk power - counting the scattering amplitude is an expansion about the chiral limit . this recovers ksw power - counting in the @xmath5 channel , where fms found it to be slowly converging . however , in the @xmath4 coupled - channels , the chiral limit has contributions from both local four - nucleon operators and from the chiral limit of ope . it is these two contributions that must be resummed using the schrdinger equation to provide the lo scattering amplitude in the @xmath4 coupled - channels   channel , so far no one has managed to dimensionally regulate the @xmath4 coupled - channels due to the presence of the @xmath14 component of the nucleon - nucleon potential . therefore , in this work we will dimensionally regulate the @xmath5 channel and use position - space square - well regularization in the @xmath4 coupled - channels . it has been shown that a singular @xmath15 potential ( @xmath16 ) , can be regulated and renormalized by a single momentum - independent square - well  @xcite . ] .    the lo lagrange density describing the single - nucleon sector and the pseudo - goldstone bosons in two - flavor qcd is @xmath17 \ + \ \lambda { \rm tr}\left[\ m_q\sigma^\dagger + m_q\sigma\ \right ] \nonumber\\ & & \ + \   n^\dagger \left(\ i\partial_0 + { \nabla^2/ 2 m_n^{(0 ) } } \ \right ) n \ + \ g_a^{(0 ) } n^\dagger \sigma\cdot { \bf a } n   \ \ \ , \label{eq : onebody}\end{aligned}\ ] ] where @xmath18 is the pion decay constant in the chiral limit , @xmath19 is the nucleon mass in the chiral limit , @xmath20 is the axial coupling constant in the chiral limit , @xmath21 is the exponential of the pion field , and @xmath22 is the axial - vector field . the properties of the nucleons and mesons have been studied in the first few orders of the chiral expansion ( a good review can be found in ref . if one is interested in the @xmath2-dependence of the nuclear force , as we are , one needs to have the chiral expansion for the nucleon mass , for the axial coupling and for the pion decay constant up to nlo . each of these observables has been studied extensively , the results of which can be found in refs . @xcite , @xmath23 \nonumber\\ m_n & = & m_n^{(0 ) } - 4 m_\pi^2 c_1 + { \cal o}\left(m_\pi^3\right ) \nonumber\\ g_a & = & g_a^{(0 ) } \left [ 1 -   { 2 ( g_a^{(0)})^2 + 1 \over 8\pi^2 ( f_\pi^{(0)})^2 }   m_\pi^2\log\left({m_\pi^2\over{\lambda^2}}\right)+{\cal o}(m_\pi^2)\right ] \ \ \ , \label{eq : snparams}\end{aligned}\ ] ] where @xmath24  @xcite , @xmath25  @xcite are @xmath2-independent constants , and @xmath26 has been used to determine the constants in the chiral limit . we use @xmath27 , @xmath28 and @xmath29 . we have retained only the leading chiral - logarithmic contribution to @xmath30 , and have chosen a renormalization scale of @xmath31 . this point requires discussion . extraction of the counterterm relevant to the @xmath2-dependence of @xmath32 at one - loop order presently yields @xmath33 ( see ref . @xcite and references therein ) . as this anomalously small value seemingly indicates a breakdown of the chiral expansion , we assume that this is a problem with the extraction and use a natural value of @xmath34 for this analysis . this , of course , introduces an uncertainty at nlo in the eft calculation . the interactions between two nucleons arise from pion exchange , resulting from the lagrange density in eq . ( [ eq : onebody ] ) , and also from local four - nucleon interactions , which for s - wave interactions result from a lagrange density of the form , @xmath35 -{1\over 2 } d_{t2 } \left ( n^\dagger \sigma^i n\right ) \left ( n^\dagger \sigma^i n\right){\rm tr}\left [ { \cal m}_{q+ } \right ] \ + \ ... \ \ \ , \label{eq : twobody}\end{aligned}\ ] ] where @xmath36 , and @xmath37 . the ellipses denote operators involving two derivatives  @xcite , which are the same order in the power - counting as a single insertion of @xmath2 , and also higher - dimension operators . the @xmath2-independent coefficients @xmath38 , @xmath39 and so forth , are to be determined from experimental data . an important point to notice is that one can not separate the contributions from the @xmath38 and the @xmath39 using nn scattering data alone , as both operators are momentum independent . however , the @xmath39 operators contribute to interactions between two nucleons and two or more pions , while the @xmath38 operators only contribute to interactions between two nucleons . one might imagine that grave violence would be done to the predictions of bbn if a di - neutron or di - proton were to be stable during the nucleosynthesis epoch . thus it is interesting to know just how much the couplings in the lagrange density describing the interactions between two nucleons in the @xmath5 channel can change before a di - neutron state becomes stable  , @xmath40 and @xmath41 bound states would have to be considered . ] . however , from the point of view of setting bounds on the variation of fundamental couplings , what is important to know is if there is a set of `` reasonable '' couplings for which the scattering length remains unnaturally large and the di - neutron is unbound over a large range of pion masses . in this situation the pion mass can vary substantially , yet still not significantly modify the physics of two nucleons in the @xmath5 channel . indeed , such a parameter set exists . = 2.5 in    0.15 in    .2 in    = 4.5 in    0.15 in    .2 in    the scattering amplitude for two nucleons in the @xmath5 channel has been computed out to n@xmath9lo  @xcite , but for our purposes it will be sufficient to work with the amplitude at nlo  @xcite . the scattering amplitude thus has an expansion of the form @xmath42 where @xmath43 is of order @xmath44 in the dual @xmath2 and momentum expansions . the lo amplitude arises from the diagrams in fig . [ fig:1s0lo ] and is given by    @xmath45 } \ \ \ , \label{eq : amplo}\end{aligned}\ ] ]    while the nlo amplitude arises from the diagrams in fig . [ fig:1s0nlo ] and is given by the sum of @xmath46 ^ 2 } \nonumber\\   { \cal a}_0^{(b ) } & = &   \left({g_a^2\over 2f_\pi^2}\right ) \left(-1 + { m_\pi^2\over 4p^2 } \ln \left ( 1 + { 4p^2\over m_\pi^2}\right)\right ) \nonumber\\ { \cal a}_{0}^{(c ) } & = & { g_a^2\over f_\pi^2 } \left ( { m_\pi m{\cal a}_{-1}\over 4\pi } \right ) \bigg ( - { ( \mu + ip)\over m_\pi } + { m_\pi\over 2p } \left[\tan^{-1 } \left({2p\over m_\pi}\right ) + { i\over 2 } \ln \left(1 + { 4p^2\over m_\pi^2 } \right)\right]\bigg ) \nonumber\\ { \cal a}_0^{(d ) } & = & { g_a^2\over 2f_\pi^2 }   \left({m_\pi m{\cal a}_{-1}\over 4\pi}\right)^2   \bigg(1 -\left({\mu + ip\over m_\pi}\right)^2 + i\tan^{-1 } \left({2p\over m_\pi}\right ) - { 1\over 2 } \ln \left({m_\pi^2 + 4p^2\over\mu^2}\right ) \bigg ) \nonumber\\ { \cal a}_0^{(e ) } & = &   { -d_2 m_\pi^2\over\left[1 + { c_0 m\over 4\pi } ( \mu + ip)\right]^2 } \ \ \ , \label{ampnlo}\end{aligned}\ ] ] where we have chosen to work in the isospin limit , @xmath47 , and we have turned off the electromagnetic interaction . these amplitudes are manifestly renormalization - scale independent order - by - order in the eft expansion . the coefficients @xmath48 , @xmath49 and @xmath50 are the combinations of couplings from the lagrange density in eq . ( [ eq : twobody ] ) appropriate for the @xmath5 channel . @xmath48 corresponds to a four - nucleon operator that is independent of @xmath2 and momentum , @xmath49 corresponds to a four - nucleon operator that has a single insertion of @xmath2 and no derivatives , and @xmath50 corresponds to a four - nucleon operator that is independent of @xmath2 and has two derivatives . the quantity @xmath51 is the renormalization scale , and we have used dimensional regularization and the power - divergence subtraction procedure ( pds )  @xcite to renormalize the theory ; @xmath52 is the magnitude of the three momentum of each nucleon in the center - of - mass frame . in addition , we have used the lo relation between the quark masses and the square of the pion mass so that the amplitude at nlo is written entirely in terms of @xmath53 . it is important to note that the operator with coefficient @xmath49 is required at nlo . it is this operator that absorbs the scale dependence of @xmath54 .    from the nlo amplitude it is easy to construct @xmath55 , which has a well - behaved power - series expansion for @xmath56 , and thus we can determine a linear combination of @xmath48 and @xmath49 in terms of the scattering length @xmath57 , and @xmath50 in terms of the effective range , @xmath58 , at the physical value of the pion mass . once these parameters are fixed , the scattering length , effective range and phase - shift can be determined as a function of the pion mass . the scattering length is given by    @xmath59   -\frac{m_n m_\pi^2}{4\pi}(\gamma -\mu ) ^2 \ d_2   \ \ , \label{eq : scattlength}\end{aligned}\ ] ]    where @xmath60 . given that it is only a combination of @xmath48 and @xmath49 that can be fixed from nn scattering , and at present @xmath49 can not be isolated with pionic processes , unique values for @xmath49 or @xmath48 that are determined experimentally do not exist . thus , one must explore a `` reasonable '' range of values for @xmath49 and @xmath48 in order to gain some understanding of how the @xmath5 channel is modified away from the physical values of the quark masses . it is quite straightforward to identify a set of couplings ,    @xmath61    for which the scattering length in the @xmath5 channel remains unnaturally large and negative over a wide range of pion masses , as shown in fig . [ fig : singscatt ] , and hence a di - neutron is not stable in this region   rather than @xmath62 , which would be more natural from the standpoint of qcd . ] . we have extracted the @xmath38 and @xmath49 using eq . ( [ eq : scattlength ] ) and an analogous equation for the effective range , with @xmath30 , @xmath63 and @xmath64 set to their chiral limit values in eq . ( [ eq : snparams ] ) . in order to determine the impact of a varying pion mass on nuclei more complex than the deuteron , it is useful to know the behavior of the phase - shift as a function of the pion mass over a relatively wide range of momentum . in fig . [ fig : singphase ] we show the @xmath5 phase - shift for @xmath65 , @xmath66 and @xmath67 , for the couplings in eq . ( [ eq : singpars ] ) . one notices that the nlo phase - shift does not reproduce the results of the nijmegen partial - wave analysis  @xcite above @xmath68 , even at the physical value of the pion mass   channel description with eft , see refs . @xcite . ] . this is because we have fixed the parameters , @xmath48 , @xmath49 and @xmath50 to the scattering length and effective range . however , this means that the dual chiral and momentum expansion of the eft at nlo has been matched to the momentum expansion of @xmath55 . consequently , both the scattering length and the effective range have their own chiral expansions that we are truncating at nlo . thus we do not expect to reproduce the phase - shift exactly , but should be perturbatively close , as can be seen to be the case in fig . [ fig : singphase ] . rather than compare the phase - shifts with the nijmegen partial - wave analysis  @xcite , it is perhaps more informative to compare the phase - shifts for @xmath65 and @xmath67 with those at @xmath69 . thus we have identified a set of nlo couplings for which the @xmath5 channel is quite insensitive to moderate changes in the pion mass . it is interesting to ask what one might expect if nature has chosen values for the couplings other than those in eq . ( [ eq : singpars ] ) . choosing some arbitrary values for the couplings , but ones that still respect nda in the two - nucleon sector  @xcite ,    @xmath70    we show the scattering length as a function of the pion mass in fig .  [ fig : singrandom ] . one can see that there are couplings for which the scattering length becomes positive , indicating the presence of a di - neutron that is stable with respect to the strong interaction , for relatively small variations in the pion mass . perhaps the unjustified choice of nn interaction arising from n@xmath9da used in previous works on this subject ( e.g. ref  @xcite ) is in some way comparable to such an ad - hoc choice of couplings in the eft . it is clear that sufficiently little is known about the @xmath2-dependence of the @xmath5 channel that placing bounds on the time - variation of @xmath2 from nuclear processes sensitive to this channel is not possible at this point in time . as we have already discussed , the deuteron plays a central role in bbn and the production of light elements . if the deuteron was deeply bound or unbound in the nucleosynthesis epoch then the abundances of the light elements would look radically different from the predictions of bbn and from those observed in nature . therefore , it is possible that limits can be placed upon the time - variation of fundamental couplings if the dependence of @xmath10 , and relevant nuclear cross - sections on these couplings is known . first realistic attempts to understand @xmath71 as a function of @xmath2 were made in refs  @xcite . one motivation for bbsvk was to understand just what was possible for the @xmath2-dependence of @xmath71 , and in particular to determine if the deuteron was bound or unbound in the chiral limit . this issue could not be resolved due to the lack of information about the @xmath49 operators that contribute to the @xmath4 coupled - channels . a second motivation was to understand what was required to make any sort of contact with quenched lattice qcd calculations  @xcite of the scattering lengths in the @xmath72 channel . = 3.63 in    0.15 in    .2 in    in the @xmath4 coupled - channels , ope generates both central and tensor potentials , @xmath73 where @xmath74 . the constant @xmath75 is uncertain  @xcite but as our results are quite insensitive to this quantity , we set it to zero in what follows . at lo in bbsvk counting we should only keep the chiral limit of these potentials and the @xmath48 operator   channel , in the @xmath4 coupled - channels there is a local four - nucleon interaction from ope that we have included in our definition of @xmath48 . ] that contributes to the @xmath4 coupled - channels , as shown in fig . [ fig:3s1lo ] . deviations from the chiral limit are inserted in perturbation theory . we choose to keep the full potentials in our nlo calculation , meaning that we also include part of the n@xmath9lo calculation , which we find does not modify the divergence structure . in addition , at nlo in bbsvk counting there is a contribution from the chiral limit of two - pion exchange ( tpe ) and from an insertion of @xmath50 and @xmath49 , as shown in fig . [ fig:3s1nlo ] . the tpe potential in coordinate space has been computed in ref . @xcite , and in the chiral limit is given by @xmath76    0.15 in    .2 in    in order to regulate the potentials at the origin , we use a spatial square - well of radius @xmath77  @xcite , where the potential outside the square well is @xmath78 where @xmath79 the potential inside the square well is @xmath80 where @xmath81 , @xmath82 and @xmath83 are constant potentials corresponding to the renormalized local operators with coefficients @xmath48 , @xmath49 and @xmath50 in the @xmath4 coupled - channels , respectively , and again we have used the lo relation between the pion mass and @xmath2 . it is important to recall that there is implicit @xmath2-dependence in this potential due to @xmath30 , @xmath63 and @xmath64 . for the chiral limit of the ope potential , the nlo @xmath2-dependence of these constants must be retained . however , for the deviations of the ope potential from the chiral limit , and also for the chiral limit of the tpe potential it is sufficient ( and in fact necessary to be able to renormalize the theory ) to keep only the chiral limit of these constants . defining @xmath84 to be @xmath85 where @xmath86 is the s - wave wavefunction and @xmath87 is the d - wave wavefunction , the regulated schrdinger equation is @xmath88 \psi ( r ) & = & 0 . \end{aligned}\ ] ] as discussed in bbsvk , we can make an identification between the coefficients of local operators , @xmath38 and @xmath39 , and the constant potentials of the square - wells that enter into eq . ( [ eq : potin ] ) , @xmath89 and @xmath90 , @xmath91    at the physical value of the pion mass , @xmath50 and a linear combination of @xmath48 and @xmath49 are fixed to reproduce the deuteron binding energy and the effective range in the @xmath72-channel . given that it is only a combination of @xmath48 and @xmath49 that we determine , we have freedom to vary @xmath49 and compensate this by a change in @xmath48 , as in the @xmath5 channel . we find couplings @xmath48 , @xmath49 and @xmath50 , renormalized at the scale set by the radius of the square - well , @xmath92 , @xmath93 that are consistent with nda and all the low - energy phase - shift data in the @xmath4 coupled - channels , and for which the deuteron is loosely bound over a large range of pion masses , as shown in fig . [ fig : deutbind ]   dependence in eq . ( [ eq : opepots ] ) is considered  @xcite . ] . the phase - shifts and mixing parameter for the @xmath4 coupled channels are shown in fig .  [ fig : tripphase ] . as expected , the @xmath72 phase - shift falls more steeply as the deuteron becomes more loosely bound . further , as the pion becomes lighter , the higher partial waves and mixing parameter are expected to increase due to the longer range of the potential and their relative insensitivity to short - distance physics .    in bbsvk , the renormalization scale dependence of the theory was investigated  @xcite . the dependence on the width of the square - well , @xmath77 , was found to be small , consistent with the eft expectations . we find the same to be true at nlo .    for purposes of contrast , we present the deuteron binding energy , phase - shifts and mixing parameter for an arbitrary set of couplings that respect nda , @xmath94 for these choices , the variation of @xmath71 with respect to @xmath2 is quite rapid , as shown in fig .  [ fig : deutbindrapid ] , and for @xmath95 , the deuteron is unbound . the phase - shifts and mixing parameter resulting from the couplings in eq . ( [ eq : triprand ] ) are shown in fig . [ fig : tripphasefast ] ,    and as expected , they have a relatively strong dependence upon @xmath2 . it is important to point out that while the deuteron binding energy curves shown in fig . [ fig : deutbind ] and fig . [ fig : deutbindrapid ] yield unbound deuterons in the chiral limit , there exist parameter sets consistent with nda for which the deuteron remains bound in the chiral limit . we have explored the light - quark mass dependence of low - energy nucleon - nucleon interactions . the motivation for this work was to determine if , in fact , bounds could be set on the time - variation of fundamental couplings from nuclear processes , such as those occurring during big bang nucleosynthesis . we have demonstrated the existence of sets of strong interaction couplings in the low - energy effective field theory describing the nucleon - nucleon interaction that are consistent with all available data and with naive dimensional analysis for which the di - neutron remains unbound , and the deuteron remains loosely bound over a wide range of light quark masses . we do not mean to imply that these are the sets of couplings that nature has chosen , but rather that this scenario is not excluded at present . thus , we conclude that bounds that have previously been set on the time - variation of fundamental couplings from processes in the two - nucleon sector are not rigorous and should be discarded .    our calculation does suffer from some limitations . first , we have not included electromagnetism , and thus could not address the @xmath41-system in the @xmath5-channel . however , we do not believe that its inclusion will modify the results we have presented in any significant way . the same we believe to be true of isospin breaking . second , and perhaps the most important limitation , is that we have not included the strange quark in our discussions , and have worked with @xmath3 chiral symmetry . variations in the strange ( and charm , bottom and top ) quark mass will manifest themselves as changes in the values of the coefficients in the lagrange density , @xmath38 and @xmath39 . at this point in time the tools do not exist to analyze this scenario . for instance , one might consider developing an eft with three light - quark flavors . however , since nuclei are such finely - tuned systems , and @xmath96-breaking is not so small , it is likely that a perturbative treatment of nuclei using the approximate @xmath97 chiral symmetry of qcd would converge very slowly , if at all .    from a conventional nuclear physics point of view , the @xmath2-dependence of the nuclear potentials not only requires knowledge of the @xmath2-dependence of one - pion exchange , but also of the @xmath2-dependence of the @xmath98 , @xmath99 , @xmath100 , ... masses and their couplings to nucleons . this would appear to be an intractable problem . it is possible that some of the traditional treatments of the nucleon - nucleon interaction  @xcite , and treatments of light nuclei could mimic these dependences by variations in their ad - hoc short - distance components of the nucleon - nucleon and higher - body potentials . the recently developed effective field theory tools organize this problem in a very simple way , and thus at nlo in the eft expansion , only the @xmath39 are needed . in experimental processes examined to date that , in principle , allow one to separate the @xmath39 from the @xmath38 , it is found that there are other amplitudes contributing to the process that dominate over the @xmath39 contributions , rendering such a separation exceedingly difficult . in the absence of an experimental determination of the @xmath39 , it would appear that the only viable means for determination of this @xmath2-dependence is through lattice qcd simulations  . we find this to be very strong motivation to pursue a lattice qcd program focused on the two - nucleon sector . we thank george fuller for useful discussions and iain stewart for helpful comments on the manuscript . we are especially grateful to evgeni epelbaum for pointing out an error in an earlier version of the manuscript , and for useful discussions . this research was supported in part by the doe grant de - fg03 - 97er41014 . beane , p.f . bedaque , w.c . haxton , d.r . phillips and m.j . savage , essay for the festschrift in honor of boris ioffe , in the encyclopedia of analytic qcd , at the frontier of particle physics , vol .  1 , 133 - 269 , edited by m.  shifman ( world scientific ) ; nucl - th/0008064 .              u .- g . meiner , essay for the festschrift in honor of boris ioffe , in the encyclopedia of analytic qcd , at the frontier of particle physics , vol .  1 , 417 - 505 , edited by m.  shifman ( world scientific ) ; hep - ph/0007092 .
Abstract:

Answer: the dependence of the nuclear force on standard model parameters plays an important role in bounding time and space variations of fundamental couplings over cosmological time scales . we discuss the quark - mass dependence of deuteron and di - neutron binding in a systematic chiral expansion . the leading quark - mass dependence of the nuclear force arises from one - pion exchange and from local quark - mass dependent four - nucleon operators with coefficients that are presently unknown . by varying these coefficients while leaving nuclear observables at the physical values of the quark masses invariant , we find scenarios where two - nucleon physics depends both weakly and strongly on the quark masses . while the determination of these coefficients is an exciting future opportunity for lattice qcd , we conclude that , at present , bounds on time and space variations of fundamental parameters from the two - nucleon sector are much weaker than previously claimed . this brings into question the reliability of coupling - constant bounds derived from more complex nuclei and nuclear processes . * variation of fundamental couplings + and nuclear forces *    * silas r.  beane and martin j.  savage *