Research article Special Issues

Consistency in Physical Activity and Increase in Mental Health in Elderly over a Decade: Are We Achieving Better Population Health?

  • Objective: Over the past century, advances in medicine and public health have resulted in an extraordinary increase in life expectancy. As a result, focus has shifted from infectious to chronic diseases. Though current guidelines for healthy behaviors among the elderly exist, it remains unclear whether this growing segment of the population has shifted their behaviors in response to public health campaigns. The objective of this study was to investigate mental health and physical activity trends that may be leading indicators for healthier living and increased life expectancy. Methods: Using nearly a decade of continuous serial cross-sectional data collected in the nationwide Behavioral Risk Factor Surveillance System, this study investigated trends of health behaviors and mental health in a population of nearly 750,000 who were 65 or older from 2003 through 2011. Weighted univariate and multivariable analyses were utilized including investigation of trend analyses over the decade, producing adjusted annual odds of physical activity and mental health. Results: After controlling for demographic and other factors, higher education and income, lower BMI, and current or previous smoking was associated with higher odds of adverse mental health and lower odds of physical activity engagement. Adjusted odds of adverse mental health climbed over the decade of observation whereas the odds of physical activity remained static. Conclusions: These data, encompassing a very large population over a decade of time, suggest that physical activity is stable though mental health challenges are on the rise in this older population. Public health campaigns may face greater barriers in an elderly population due to lifelong habits, dissemination and educational approaches, or decreasing gains. Further research should be conducted to identify more effective approaches towards increasing physical activity in this important and growing subset of the population and towards transforming behaviors earlier in life.

    Citation: Tyler C. Smith, Besa Smith. Consistency in Physical Activity and Increase in Mental Health in Elderly over a Decade: Are We Achieving Better Population Health?[J]. AIMS Medical Science, 2016, 3(1): 147-161. doi: 10.3934/medsci.2016.1.147

    Related Papers:

    [1] Helen Christodoulidi, Christos Efthymiopoulos . Stages of dynamics in the Fermi-Pasta-Ulam system as probed by the first Toda integral. Mathematics in Engineering, 2019, 1(2): 359-377. doi: 10.3934/mine.2019.2.359
    [2] Lorenzo Pistone, Sergio Chibbaro, Miguel D. Bustamante, Yuri V. Lvov, Miguel Onorato . Universal route to thermalization in weakly-nonlinear one-dimensional chains. Mathematics in Engineering, 2019, 1(4): 672-698. doi: 10.3934/mine.2019.4.672
    [3] Camilla Nobili . The role of boundary conditions in scaling laws for turbulent heat transport. Mathematics in Engineering, 2023, 5(1): 1-41. doi: 10.3934/mine.2023013
    [4] Raimondo Penta, Ariel Ramírez-Torres, José Merodio, Reinaldo Rodríguez-Ramos . Effective governing equations for heterogenous porous media subject to inhomogeneous body forces. Mathematics in Engineering, 2021, 3(4): 1-17. doi: 10.3934/mine.2021033
    [5] Federico Cluni, Vittorio Gusella, Dimitri Mugnai, Edoardo Proietti Lippi, Patrizia Pucci . A mixed operator approach to peridynamics. Mathematics in Engineering, 2023, 5(5): 1-22. doi: 10.3934/mine.2023082
    [6] Jana Kopfová, Petra Nábělková . Thermoelastoplastic oscillator with Prandtl-Ishlinskii operator. Mathematics in Engineering, 2025, 7(3): 264-280. doi: 10.3934/mine.2025012
    [7] Daniele Agostinelli, Roberto Cerbino, Juan C. Del Alamo, Antonio DeSimone, Stephanie Höhn, Cristian Micheletti, Giovanni Noselli, Eran Sharon, Julia Yeomans . MicroMotility: State of the art, recent accomplishments and perspectives on the mathematical modeling of bio-motility at microscopic scales. Mathematics in Engineering, 2020, 2(2): 230-252. doi: 10.3934/mine.2020011
    [8] Paola F. Antonietti, Chiara Facciolà, Marco Verani . Unified analysis of discontinuous Galerkin approximations of flows in fractured porous media on polygonal and polyhedral grids. Mathematics in Engineering, 2020, 2(2): 340-385. doi: 10.3934/mine.2020017
    [9] Antonello Gerbi, Luca Dedè, Alfio Quarteroni . A monolithic algorithm for the simulation of cardiac electromechanics in the human left ventricle. Mathematics in Engineering, 2019, 1(1): 1-37. doi: 10.3934/Mine.2018.1.1
    [10] Luca Azzolin, Luca Dedè, Antonello Gerbi, Alfio Quarteroni . Effect of fibre orientation and bulk modulus on the electromechanical modelling of human ventricles. Mathematics in Engineering, 2020, 2(4): 614-638. doi: 10.3934/mine.2020028
  • Objective: Over the past century, advances in medicine and public health have resulted in an extraordinary increase in life expectancy. As a result, focus has shifted from infectious to chronic diseases. Though current guidelines for healthy behaviors among the elderly exist, it remains unclear whether this growing segment of the population has shifted their behaviors in response to public health campaigns. The objective of this study was to investigate mental health and physical activity trends that may be leading indicators for healthier living and increased life expectancy. Methods: Using nearly a decade of continuous serial cross-sectional data collected in the nationwide Behavioral Risk Factor Surveillance System, this study investigated trends of health behaviors and mental health in a population of nearly 750,000 who were 65 or older from 2003 through 2011. Weighted univariate and multivariable analyses were utilized including investigation of trend analyses over the decade, producing adjusted annual odds of physical activity and mental health. Results: After controlling for demographic and other factors, higher education and income, lower BMI, and current or previous smoking was associated with higher odds of adverse mental health and lower odds of physical activity engagement. Adjusted odds of adverse mental health climbed over the decade of observation whereas the odds of physical activity remained static. Conclusions: These data, encompassing a very large population over a decade of time, suggest that physical activity is stable though mental health challenges are on the rise in this older population. Public health campaigns may face greater barriers in an elderly population due to lifelong habits, dissemination and educational approaches, or decreasing gains. Further research should be conducted to identify more effective approaches towards increasing physical activity in this important and growing subset of the population and towards transforming behaviors earlier in life.


    This paper is dedicated to our colleague and friend Antonio Giorgilli, in occasion of his 70$ ^{\mathrm{th}} $ birthday. The scientific contribution of Antonio, in different branches of Mathematical Physics, is well known and is invaluable. Concerning in particular FPU, beyond the interest of specific results, Antonio had a special propulsive role in Italy, and influenced very much our own work. "Understanding the FPU state...", which is the title of the present paper, started for us from studying [14] and the scaling law (3) there illustrated. Thanks, Antonio!

    This paper is devoted to the Fermi-Pasta-Ulam (FPU) problem [1], namely the problem of energy sharing in weakly nonlinear chains of oscillators. The Hamiltonian has the form

    $ H(p,q)=12Ni=1p2i+Ni=0V(qi+1qi) ,
    $
    (1)

    $ V $ being some nearest-neighbors potential with single minimum in zero,

    $ V(r)=r22+αr33+βr44+ ,β>0 .
    $
    (2)

    In (1) the boundary conditions are still not specified; they are generally either fixed ends, i.e., $ q_0 = q_{N+1} = 0 $ (like in the original FPU paper), or periodic, $ q_{N} = q_0 $. The fixed-ends case can be viewed as a submanifold of the periodic case with $ N' = 2(N+1) $ degrees of freedom, with symmetry $ q_{N'-j} = -q_{j} $.

    The value of $ \alpha $, if different from zero, is irrelevant, since a trivial rescaling reports it to any prefixed value; the effective parameters determining the dynamics are indeed $ |\alpha|\sqrt{\varepsilon } $, $ \varepsilon = E/N $ being the specific energy, and then $ \beta/\alpha^2 $, … The choice of $ \alpha $ actually fixes the energy scale, as well as the scale of $ \beta $ and of possibly further coefficients in (2).

    As is well known, Fermi, Pasta and Ulam aimed to investigate how the system reaches the statistical equilibrium, identified with the equipartition of energy among normal modes, if started very far from equilibrium, the whole energy being initially given to only one or two long-wavelength normal modes. With great surprise they found, within the time scale accessible to their computer, * no equilibrium at all: the system apparently reached a stationary state, different both from the initial state and from equipartition, in which only a few normal modes significantly share energy. Such a peculiar state is often called "the FPU state". Besides the lack of equipartition, the FPU state is characterized by long time recurrences, made especially evident a few years later in [3].

    *The celebrated MANIAC I of the Los Alamos Laboratories, capable to process approximately $ 10^4 $ instructions per second [2]. Compare with modern computers, having hundreds or thousands of CPUs with clock at several GHz.

    The literature on FPU, 65 years after the original FPU paper, is so huge, that it is not even conceivable to summarize it through a reasonable number of main references. Two collections of papers [4,5], produced in the occasion of the 50th anniversary of FPU, may give an idea of the status of the art around 2004–2005, as well as of the widespread attitudes of different authors, or different schools, in front of such non easy subject which touches the very foundations of classical statistical mechanics.

    Two main questions have been extensively studied in the literature, since the very beginning:

    ⅰ. Understanding the unexpected FPU state: its characteristic features, its long time recurrences, how is it produced from initial data, how is it possible that, in spite of the nonlinearity which couples normal modes, energy is shared among only a few of them.

    ⅱ. Understanding whether, possibly on a much larger time scale, the FPU state further evolves, eventually reaching thermalization and normal statistical behavior, as in the original FPU expectation.

    Both questions have been studied rather soon, independently. On the one hand, the recurrences observed in [1,3] suggested the presence of a close underlying integrable dynamics, which was initially searched, starting with [6] in 1965, among nonlinear wave equations, and later on, starting with [8] and [9], in the Toda model. On the other hand, almost simultaneously in 1966, the possibility was outlined [7] that normal statistical behavior is recovered if the specific energy $ \varepsilon = E/N $ exceeds a threshold, which in turn, in [7], was conjectured to vanish for $ N\to\infty $.

    Ref. [7], as well as many later papers investigating the conditions for thermalization in FPU, deal with the so-called $ \beta $ model, having $ \alpha = 0 $. The two models in fact (together with other variants of FPU, like Lennard-Jones) quite often have been considered equivalent and interchangeable, without much criticism. Certainly, they are not for what concerns the FPU state, which is peculiar of the $ \alpha $-models. The long time behavior, with thermalization on a suitably long time scale, is instead rather similar; see for a comparison [10].

    The literature on FPU accumulated, along these lines, without much order, and it is not so easy to follow it. Even the basic distinction between formation of the FPU state and possible later evolution to normal statistical behavior, with different phenomena occurring on different time scales, was not always clear at the beginning, and started being emphasized only in 1982 in [11,12], where a possible connection between the FPU state and metastability was outlined. Not so many papers are based on this distinction; among them, [10,13,14,15,16,17,18].

    In this paper we focus the attention on the former problem, namely understanding as much as possible the FPU state. The aim is to revisit the existing literature, tracing possible paths through it, and then to provide a further understanding of the FPU state by looking more closely at the Toda model.

    Browsing the literature, we extract two main ideas (see Section 2 for references):

    $ – $ The FPU state is characterized by a few simple scaling laws; the most relevant one is that the number $ M $ of normal modes effectively involved in the energy sharing, when one or a few long-wavelength normal modes are initially excited, scales as

    $ M  Nε1/4 .
    $
    (3)

    $ – $ The reference integrable dynamics underlying FPU, is not the linear one, made of free normal modes, rather the essentially nonlinear Toda dynamics.

    Our contribution, in this background, is the following:

    $ – $ We check in Toda the scaling laws known to hold for FPU, and find they are fully satisfied. Scaling laws appear thus to be proper of Toda, and FPU obeys them as long as its dynamics remains close to the Toda dynamics.

    $ – $ We exploit the numerical procedure explained in [9] to compute numerically the Toda actions; the procedure allows to determine, in any given state $ (p, q) $, the number of nonvanishing actions, that is the effective dimensionality of the corresponding Toda torus. We consider states in which only one or a few low frequency modes are excited, and observe that the dimensionality $ {\cal M} $ of the torus scales exactly as in (3),

    $ M  Nε1/4 .
    $
    (4)

    This means: the scaling laws characterizing the FPU state are contained in the nontrivial correspondence between normal modes and Toda actions. This is a more innovative view, in which dynamics (integration of Hamilton equations) does not play a role. By the way, the accuracy of the scaling law (4) is the best we ever met (see in particular Figures 8 and 9 below).

    The paper is organized as follows: Section 2, after a few preliminaries, proposes our personal way to look at the literature on the FPU state, and also includes (Sect. 2.4) the new results concerning the scaling laws in the Toda model. Section 3 is instead devoted to the computation of the Toda actions. An Appendix follows, in which the procedure to compute the Toda actions is shortly recalled.

    Consider, to be definite, the case of fixed ends. Normal modes are then

    $ Q_k = \sqrt{\frac{2}{N+1}}\,\sum\limits_{i = 1}^N q_i\,\sin\frac{\pi ki}{N+1}\ , \qquad P_k = \sqrt{\frac{2}{N+1}}\,\sum\limits_{i = 1}^Np_i\,\sin\frac{\pi ki}{N+1}\ ; $

    their energies $ E_k $ and frequencies $ \omega _k $ are

    $ E_k = \frac12(P_k^2+\omega _k^2Q_k^2)\ , \qquad \omega _k = 2\,\sin\frac{\pi k}{2(N+1)}\ . $

    For given initial data let

    $ \overline E_k(T) = \frac1{\frac13 T}\int_{\frac23 T}^T E_k(P(t),Q(t)) \,\mathrm{d} t\ . $

    Restricting the average to the running window $ \frac23 T\le t \le T $ does not change the asymptotics and makes a little quicker the convergence: time averages from the beginning are somehow lazy.

    A preliminary qualitative view of the formation of the FPU state, and of its later slow evolution to equipartition, is provided by Figure 1 (a remake of Figure 1 in [19]) which refers to $ N = 1023 $, $ \alpha = -1 $, $ \beta = 2 $, specific energy $ \varepsilon = 10^{-4} $, and shows $ \overline E_k(T) $ vs. $ k/n $ at selected times $ T $ (marked in the figure) in geometric progression. The energy was initially equidistributed among the lowest $ 10\% $ of modes (the rectangular profile in the figure). One can see that quite soon, already at $ T\simeq10^3 $, a well defined profile is formed, in which only some low frequency modes effectively take part to the energy sharing, the energies of the remaining ones decaying exponentially with $ k/N $. The energy profile keeps its form nearly unchanged for a rather large time scale, $ T\simeq10^5 $ or $ 10^6 $, definitely much larger than the time needed to form it. Afterwords, on a much larger time scale, the dynamics slowly evolves towards energy equipartition, the high-frequency modes being progressively involved into the energy-sharing game; in the above conditions, equipartition requires $ T\simeq10^{10} $. The initial phases (this will have importance, see later) were chosen randomly; moreover, to improve a little the figure reducing the fluctuations, an average of $ E_k(T) $ on 24 different random choices of the phases was performed.

    Figure 1.  The shape of the averaged energy spectrum of normal modes $ \overline E_k(T) $ plotted vs. $ k/N $, at selected times $ T $ (marked in the figure) in geometric progression. Fixed ends FPU, $ N = 1023 $, $ \alpha = 1 $, $ \beta = 2 $, $ \varepsilon = 10^{-4} $. Energy initially equidistributed among modes with $ 0 < k/N < 0.1 $, see the rectangle marked $ t = 0 $. Each point is the average over 24 random extractions of the initial phases.

    If the initial energy distribution is narrower, then the formation of the FPU state is slower, but the state looks the same, as well at the later evolution; see Figure 2, left panel, where the initial distribution is $ 0 < k/N < 0.03 $. If energy is raised, the width of the FPU state grows, see the right panel of Figure 2 where $ \varepsilon = 4\times10^{-4} $.

    Figure 2.  Left: same as in Figure 1, with narrower initial energy distribution, $ 0 < k/N < 0.03 $. Right: same as in Figure 1, at higher energy $ \varepsilon = 4\times10^{-4} $.

    We shall now try to order the literature on the FPU state, by tracing two possible paths through it.

    The FPU state needs to be characterized quantitatively, looking in particular at its width——how many modes are effectively involved——and at the time scale for its formation, to be thought of as depending on $ N $ and $ \varepsilon $ as well as on the kind of the initial excitation.

    Let us start by numerical evidence. Both in [20] and [14] a convincing evidence is produced (particularly accurate in [14]) that the width of the FPU state, namely $ f = M/N $ where $ M $ is the (suitably defined) number of modes effectively sharing energy, for large $ N $ depends only on the specific energy $ \varepsilon $ and not on $ N $, and scales as

    $ f  ε1/4 ;
    $
    (5)

    correspondingly, $ M $ obeys (3). One should remark that [20] deals with a periodic pure $ \alpha $ model, and considers an initial excitation flat on a fraction $ f_0 $ of modes already proportional to $ \varepsilon ^{1/4} $, while [14] deals with a fixed ends model with $ \beta = 4\alpha^2 $, and only the first mode is initially excited. The notion of "number of modes effectively shearing energy", that is the definition of $ M $, is also different, see the original papers; it is remarkable that the scaling law (5) is the same——the phenomenon looks, so to speak, robust. Some numerical evidence of (5) was already present in [21]. Concerning the time of formation of the FPU state, in [20] we find the estimate

    $ t  ε3/4 ,
    $
    (6)

    holding for large $ N $ and initial data as before, that is for $ f_0 $ obeying (5).§

    §No explicit estimate of $ t $ as function of $ \varepsilon $ is provided in [14], but an indication is present, if we correctly interpret, $ t\sim\varepsilon ^{-1/2} $. We conjecture the authors have seen the correct law $ t\sim\varepsilon ^{-3/8} $ valid for fixed $ M_0 $, see below, which is close to $ \varepsilon ^{-1/2} $ and easily confused with it in a non systematic study.

    As is remarkable, the scaling laws (5) and (6) have been also justified theoretically, looking at the way normal modes interact. More precisely, (5) is deduced in [21] by adapting to the $ \alpha $ FPU model the Chirikov method of overlapping of resonances, originally formulated for the $ \beta $ model [7]; both (5) and (6) are instead deduced in [13], on which [20] is based, by adapting to FPU some ideas from the theory of weak turbulence. In both cases, the reference uncoupled model is the harmonic chain.

    The scaling laws entering the formation of the FPU state have been later studied in [22] in a more systematic way. The conventional definition of $ M $——actually the same as in [20]——is as follows: let

    $ \eta = -\sum\limits_{k = 1}^{N}F_k\log F_k\ , \qquad \hbox{where} \quad F_k = \frac{ \overline E_k}{\sum_j \overline E_j}\ , $

    denote the "entropy" of the spectrum [23]; then

    $ M = e^\eta\ . $

    Let $ \overline E_k = E/M $ for $ k\le M $, $ \overline E_k = 0 $ for $ k > M $. Then it is $ \eta = \log M $.

    Consider initial data with the first lowest $ M_0 $ modes initially excited; think $ f $ to be a function of $ t $, $ N $, $ \varepsilon $ and $ f_0 = M_0/N $. The results in [22] can be summarized in two general scaling laws for $ f(t, N, \varepsilon , f_0) $, both holding for large $ N $, depending on whether the phases of the initially excited modes are chosen randomly or are coherent.

    (a) For random initial phases, $ f $ is independent of $ N $ and scales as

    $ f(t,N,ε,f0)=ε1/4G(ε3/8f3/20t) ,
    $
    (7)

    where $ {\cal G} $ is a kind of universal function of sigmoidal shape, reaching a plateau when its argument grows to infinity. According to (7), the width of the FPU state obeys (5), no matter which is $ f_0 $. The time scale for the formation of the state scales instead as

    $ t  ε3/8f3/20 .
    $
    (8)

    The main interest, having in mind thermodynamics, is for fixed $ f_0 $, that is $ M_0 $ proportional to $ N $ (modes excited up to a given frequency). The scaling law (8) however also holds for $ f_0\sim\varepsilon ^{1/4} $, as in [20], giving in this case (6), as well as down to small fixed $ M_0 $, as in [14], giving in such case times growing with $ N $,

    $ t  ε3/8N3/2 ;
    $
    (9)

    for fixed small $ M_0 $ the assumption that the initial phases are random is meaningless, and is in fact irrelevant to deduce (9).

    (b) If instead many modes are initially excited with coherent initial phases, for example equal phases, or phases alternating in some easy pattern, then the total energy $ E $, rather than the specific energy $ \varepsilon $, appears to be the the relevant variable, and (7) is replaced by

    $ f(t,N,ε,f0)=E1/4f1/40G(E3/8f15/80t) .
    $
    (10)

    For large $ E $, in particular for fixed $ \varepsilon $ and large $ N $, the intermediate FPU state practically confuses with equipartition among all modes: unless of course $ f_0\sim 1/N $, in which case, with interesting continuity, (9) is recovered.

    It is worthwhile to remark that the choice of equal initial phases, although not well motivated, is rather common in the literature: it occurs indeed whenever, for example, either the $ q_i $'s or the $ p_i $'s initially vanish, and correspondingly energy is purely kinetic or purely potential (the easiest initial data to be handled numerically). With such initial data, unless $ f_0\sim N^{-1} $, the FPU state for large $ N $ disappears. The lack of attention to initial phases, in our opinion, is among the reasons why it not so easy, sometimes, to compare papers.

    In [22] the general scaling laws (7) and (10) are also motivated theoretically, but in a way substantially different from [21] and [13], namely using, as a normal form for FPU, the KdV nonlinear wave equation. Considering FPU as a perturbation of a nonlinear integrable system is indeed an old idea, that we shall shortly revisit in the next subsection.

    As already remarked in the Introduction, to explain the presence of long time recurrences in FPU it was soon suggested that the FPU dynamics is close to the dynamics of a nonlinear but nevertheless integrable system, and this was first searched among nonlinear wave equations. The first attempt in this direction, going back to 1965, is [6], a fundamental paper which is at the basis of modern theory of nonlinear integrable wave equations, in which FPU appears as a main motivation to study again, after years, the KdV equation,

    $ u_t = \tfrac{\alpha}{\sqrt{2}} uu_x+\tfrac{1}{24}u_{xxx} \ . $

    The progress in the field is then rapid: a couple of years later the method of inverse scattering [24] and Lax pairs [25] are introduced, and the presence of infinitely many constants of motion is established [26]. Finally, in 1971, KdV is shown to be an infinite dimensional completely integrable Hamiltonian system [27].

    In fact, the nonlinear wave equation which is more immediately related to FPU, if one searches for a continuum limit in which the first nonlinear and the first dispersive terms beyond the wave equation are kept, is the Boussinesq equation; a possible form is

    $ u_{tt} = u_{xx}+2\alpha u_x u_{xx}+\frac{1}{12}u_{xxxx}\ . $

    From Boussinesq it is possible to deduce, in a convenient nontrivial limit, KdV, but the Boussinesq equation itself was soon proved, in 1973, to be integrable [28]. In [28] the connection between FPU and Boussinesq is quite emphasized, the lack of integrability of FPU being explained as due to the discretization.

    In the meanwhile, statistical physicists become interested in the Toda model [29,30,31,32]. As is well known, this is a Hamiltonian system with the same form as (1), $ V $ being the Toda exponential potential

    $ V_T(r) = \frac{1}{\lambda^2}(e^{\lambda r}-1-\lambda r)\ . $

    In 1974 the Toda model was proved to be completely integrable, remarkably in three independent papers [8,33,34]. Reference [8] is particularly important for FPU, because the connection with FPU is there emphasized. Indeed, for $ \lambda = 2\alpha $ it is

    $ V_T(r) = \frac12 r^2+\frac13 \alpha r^3+\frac14\beta_T r^4 +\frac15\gamma_Tr^5+\cdots\ , $

    with

    $ \beta_T = \tfrac23\alpha^2\ ,\qquad \gamma_T = \tfrac13 \alpha^3\ , \qquad \ldots \ , $

    so the model is tangent to FPU up to the third order included and provides an integrable approximation to FPU better than the harmonic chain. In [8] the slow stochastization of FPU is no more attributed to the discretization with respect to an integrable wave equation, but to the small difference, with dominating term $ \frac14(\beta-\beta_T)r^4 $ in the potential, between FPU and Toda. The reference to Toda as the best integrable approximation for FPU, is a considerable change of paradigm in the FPU problem.

    The connection between FPU and Toda was proposed again, and emphasized, in 1982, in [9]. In this paper, on the one hand a striking evidence is provided that the dynamics of FPU, in the FPU state before the slow drift towards equipartition, is hardly distinguishable from the Toda dynamics; on the other hand, the integrability of Toda is studied very constructively, and an algorithm is proposed to compute numerically the Toda actions. We shall come back to this point in Section 3.

    The obvious conjecture, at this point, is that the scaling laws (7) and (10), characterizing the FPU state and implying in particular (5) and (6), are not proper of FPU but come from the underlying Toda integrable dynamics. This is what we shall see in this section, where we revisit and remake for Toda some numerical work originally done in [22] for FPU.

    The scaling laws (7) and (10), reducing the variables from four to two, are obtained in [22] as combinations of three independent homogeneity relations characterizing the function $ f(t, N, \varepsilon , f_0) $ for large $ N $. Two of them do not depend on the choice of the phases of the excited modes and are

    $ f(λ3/4t,λ1/4N,λε,λ1/4f0)=λ1/4f(t,N,ε,f0)
    $
    (11)
    $ f(μ3/2t,μN,ε,μ1f0)=f(t,N,ε,f0) .
    $
    (12)

    We cannot enter here the way (11) and (12), which might appear obscure, are (heuristically) deduced from KdV, demanding for this to [22].

    Figure 3 provides a test of the scaling law (11). The left panel shows $ f $ vs. $ t $ for eight different choices of $ N $, $ \varepsilon $ and $ f_0 $ obeying (11), namely $ N = 2896-9742 $ with step $ 2^{-1/4} $, $ \varepsilon = 10^{-3}-7.8\times10^{-6} $ with step $ 1/2 $, and $ f_0 = M_0/N $ with constant $ M_0 = 102 $. The right panel shows the same curves rescaled according to (11), namely rescaled horizontally by a factor $ (\varepsilon /\varepsilon _*)^{-3/4} $ and vertically by a factor $ (\varepsilon /\varepsilon _*)^{1/4} $, with $ \varepsilon _* = 2.5\times10^{-4} $ (the value of $ \varepsilon $ of the third curve from above, which does not move). One can see the rescaled curves superimpose almost exactly. The phases were chosen randomly, and each curve is in fact the average over 24 different random choices of the phases.

    Figure 3.  A test of the scaling law (11) in the Toda model. Parameters: $ N = 2896-9742 $, step $ 2^{-1/4} $; $ \varepsilon = 10^{-3}-7.8\times10^{-6} $, step $ 1/2 $; $ M_0 = 102 $. Random phases (24 choices). The plotted quantity is $ f $ vs. $ t $ as it is (left) and rescaled (right).

    Figure 4 provides similarly a test of the scaling law (12). Here $ \varepsilon = 2.4\times10^{-4} $ is kept fixed, while $ N $ ranges between $ 1023 $ and $ 16383 $; $ f_0 = M_0/N $ with fixed $ M_0 = 26 $. The left panel shows the curves as they are, while in the right panel $ t $ is scaled by a factor $ (N/N_*)^{-3/2} $, with $ N_* = 4095 $ (corresponding to the central curve, which stays unchanged). As before, each curve is the average of 24 different random choices of the initial phases.

    Figure 4.  A test of the scaling law (12) in the Toda model. Parameters: $ N = 1023-16383 $, $ \varepsilon = 2.5\times10^{-4} $, $ M_0 = 26 $. Random phases (24 choices). The plotted quantity is $ f $ vs. $ t $ as it is (left) and rescaled (right).

    Putting together (11) and (12) one gets

    $ f(t,N,ε,f0)=ε1/4F(ε3/8N3/2t,Nf0) ,
    $
    (13)

    with a convenient $ {\cal F} $ depending only on two variables. Such a relation can be further tested by considering combinations of $ N $, $ \varepsilon $ and $ f_0 $ different from the previous ones, for example one can keep $ N $ constant and move $ \varepsilon $; the corresponding scaling factors are in this case $ \varepsilon ^{-3/8} $ on the horizontal axis, $ \varepsilon ^{1/4} $ on the vertical axis. Figure 5 shows the result for $ N = 1023 $, seven values of $ \varepsilon $ in the range $ 10^{-3}-1.56\times 10^{-5} $, $ f_0 = M_0/N $ with small $ M_0 = 3 $. As before the left panel shows $ f $ vs. $ t $ as it is, while the right panel reports the rescaled curves.

    Figure 5.  A test of the scaling law (13) in the Toda model. Parameters: $ N = 1023 $, $ \varepsilon = 10^{-3}-1.56\times10^{-5} $, $ M_0 = 3 $. The plotted quantity is $ f $ vs. $ t $ as it is (left) and rescaled (right).

    The third relation which completes (11) and (12) depends instead in an essential way on the choice of the phases. For random phases, it simply tells that, for large $ N $ and given $ f_0 $, $ f $ is independent of $ N $. This means $ {\cal F}(x, y) = {\cal G}(xy^{3/2}) $, and (7) is deduced. A test of the independence of $ f $ on $ N $ for large $ N $ at fixed $ f_0 $ is provided in Figure 6, where $ \varepsilon $ is either $ 10^{-3} $ or $ 2.5\times10^{-4} $, while $ N $ assumes five different values in the range $ 1023-16383 $; $ f_0 = 0.025 $ is fixed. The left panel displays $ f $ vs. $ t $; quite clearly, $ f $ depends on $ \varepsilon $ but not on $ N $. The further rescaling of the vertical axis by a factor $ \varepsilon ^{1/4} $, shown in the right panel, makes all curves superimpose.

    Figure 6.  A test of the independence of $ f $ on $ N $, at fixed $ f_0 $, for random initial phases. Parameters: $ \varepsilon $ either $ 10^{-3} $ or $ 2.5\times10^{-4} $, five different values of $ N $ in the range $ 1023-16383 $, $ f_0 = 0.025 $. Left panel: $ f $ vs. $ t $ as it is; right panel: vertical rescaling proportional to $ \varepsilon ^{-1/4} $.

    For coherent phases instead the additional scaling law is

    $ f(t,ρN,ρ1ε,f0)=f(t,N,ε,f0) ;
    $
    (14)

    that is, $ f $ is independent of $ N $ at fixed $ f_0 $ only if the total energy $ E = N\varepsilon $ is kept fixed. From (14) and (9), (10) is deduced. A test of (14) is exhibited in Figure 7, left panel, where $ E = 0.265 $ and $ f_0 = 0.1 $ are kept fixed, while $ N $ assumes six values in the range $ 511-16383 $. The right panel shows that if $ \varepsilon $, instead of $ E $, is fixed, then the independence of $ N $ is lost and for large $ N $ the FPU state does not exist any more. The phases are here constant, namely $ \phi_k = \pi/2 $ for all excited modes (purely potential initial energy), but a qualitative similar picture is obtained by choosing for example $ \phi_k = 0 $, or $ \phi_k = k\pi/2 $, or similar regular patterns.

    Figure 7.  Left: $ f(t, N, \varepsilon , f_0) $ vs.$ t $ for seven values of $ N $ in the range $ 255-16383 $, fixed $ E = 0.256 $ and $ \varepsilon = E/N $, fixed $ f_0 = 0.1 $. Right: fixed $ \varepsilon = 2.4\times10^{-4} $ instead of fixed $ E $.

    In this section we shall investigate the basic scaling law (5) following a substantially different road, namely looking at the correspondence between linear normal modes and action variables in the Toda model. The aim is understanding the scaling law (5) governing the FPU state in Toda by looking only at the initial state——how many actions are excited, which is the effective dimensionality of the torus supporting the dynamics——with no need to integrate the equations of motion.

    Let us remark that computing the Toda actions for large $ N $ is a rather nontrivial work, facing delicate questions of numerical precision. We present here only a few preliminary results, without entering in details, with the only aim to show this is a promising way (see also [35]). We shall come back soon to the subject, in a wider devoted paper.

    The appropriate definition of Toda actions, as well as the numerical procedure to compute them, are provided in [9]; for the convenience of the unfamiliar reader, however, the essentials are reported in the Appendix. The natural model to look at, in this connection, is a periodic chain with $ N $ particles, that is, eliminating the barycenter, with $ N-1 $ degrees of freedom. Let us recall that for the corresponding harmonic periodic chain, normal modes are not stationary waves but traveling waves; a convenient choice is

    $ q_i(t) = \sum\limits_{k = 1}^{N-1}A_k \, q^{(k)}_i(t)\ , \qquad i = 1,\ldots,N\ , $

    where

    $ q^{(k)}_i(t) = \sqrt{\frac{2}{N}}\,\frac{1}{\omega _k}\, \sin\Big(\frac{2\pi ik}{N}+\omega _k t+\theta_k\Big) \ , \qquad \omega _k = 2\sin\frac{\pi k}{N} $

    (similarly for $ p_i(t) $). The harmonic energies and harmonic actions are respectively

    $ E_k = A_k^2 \ , \qquad I_k = E_k/\omega _k\ . $

    Let $ {\cal I}_1, \ldots, {\cal I}_{N-1} $ denote the Toda actions, as defined in [9]. They turn out to be perturbations of the harmonic actions $ I_1, \ldots, I_n $, and in particular, if a single harmonic action $ I_k $ is different from zero, then in the limit of vanishing energy $ E $, it is

    $ IkIk1 ,IkIk0forkk .
    $
    (15)

    The limit however, see [36,37,38], is highly non uniform in $ N $, and for small $ k $, to get convergence, $ EN^3 $, equivalently $ \varepsilon N^4 $, needs to be small. The perturbative approach instead fails, and the correspondence between normal modes and actions gets out of theoretical control, for large $ N $ at fixed $ \varepsilon $, that is in the thermodynamic limit.

    Having in mind the FPU state and the scaling law (5) it satisfies, we naturally expect that for large $ N $, if one or a small number $ M_0 $ of waves of low $ k $ are present, then the number of Toda actions that are substantially different from zero scales as $ N\varepsilon ^{1/4} $. This is precisely what we found in our numerical computations.

    To define and compute the effective number $ {\cal M} $ of non vanishing actions in any given state of the system, we found convenient to resort to the analog of spectral entropy, namely

    $ {\cal H} = -{ \sum\limits_k} {\cal F}_k\log {\cal F}_k\ , \qquad \hbox{where} \quad {\cal F}_k = \frac{ {\cal E}_k}{\sum\limits_j {\cal E}_j}\ ,\quad {\cal E}_k = \omega _k {\cal I}_k\ ; $

    then we let, in complete analogy, $ {\cal M} = e^{ {\cal H}} $.

    Figure 8 reports $ {\cal M} $ vs. $ \varepsilon $ in the wide range $ 10^{-15}\le\varepsilon \le1 $, for several values of $ N $ between $ 32 $ and $ 8192 $; the state is a single traveling wave with $ k = 1 $. The solid lines drawn in the figure are obtained by a least square fit, and all slopes, denoted $ a $ and reported in the figure aside $ N $, turn out to be very close to $ 1/4 $. So, the behavior is as expected: in agreement with (15), all curves converge to $ {\cal M} = 1 $ for $ \varepsilon \to0 $, while for larger $ \varepsilon $, with crossover depending on $ N $, $ {\cal M} $ grows proportionally to $ \varepsilon ^{1/4} $.

    Figure 9, left panel, is the same as Figure 8, but the fraction $ f = {\cal M}/(N-1) $ of nonvanishing actions, instead of $ {\cal M} $, is there reported vs. $ \varepsilon $. Quite clearly all curves flatten, by growing $ \varepsilon $, on the same line. The solid line drawn in the figure interpolates the data for $ N = 8192 $ and has slope $ a = 0.25 $. An overall check that the scaling law is

    $ M  (N1)ε1/4
    $
    (16)

    is provided by the right panel of Figure 9, where $ {\cal M} $ is reported vs. $ (N-1)\varepsilon ^{1/4} $. Here too, the line interpolates data for $ N = 8192 $; its computed slope is $ 0.998 $. In conclusion, in the state of single traveling wave with $ k = 1 $, the expected scaling law is satisfied with quite great accuracy.

    Such an accuracy was indeed a nice surprise, and was not really expected, the numerical procedure to compute the actions being new for us and somehow problematic (see the Appendix, Section C). It is hard to say how reliable are the values of the slopes reported in Figure 8. The only easy analysis is the uncertainty of the least-squares fit giving $ a $. Figure 10 reports the computed slopes vs. $ N $, the error bars corresponding, as is typical, to $ \pm3 $ std. deviations; the numerical results are clearly compatible with $ a = 0.25 $. Such an analysis of the accuracy does not take into consideration other possible sources of error: At the moment we are not able to draw error bars around the symbols of Figure 8, and the choice of the intervals of $ \varepsilon $ where performing the best-fit analysis, after the crossover with the horizontal line $ N = 1 $, is also somehow arbitrary. Anyhow, the overall precision looks remarkable.

    Figure 8.  The effective number $ {\cal M} $ of nonvanishing actions vs. $ \varepsilon $, for several values of $ N $ in the range $ 32-8192 $. Single traveling wave, $ k = 1 $. Solid lines are best fit interpolations, slopes are reported in the figure aside $ N $.
    Figure 9.  Left: as in Figure 8, $ f = {\cal M}/(N-1) $ in place of $ {\cal M} $. Right: as in Figure 8, $ (N-1)\varepsilon ^{1/4} $ in the abscissa.
    Figure 10.  The slopes of the lines appearing in Figure 8 vs. $ N $, with error bars given by the best-fit analysis ($ \pm 3 $ std. deviations).

    Nothing essentially changes if, in place of a single traveling wave with $ k = 1 $, one considers, for example, three waves with $ k = 1, 2, 3 $, having the same $ {\cal E}_k $. Figure 11 shows, in the left panel, the analog of Figure 8, that is $ {\cal M} $ vs. $ \varepsilon $, and in the right panel the same as Figure 9, right panel, that is $ {\cal M} $ vs. $ N\varepsilon ^{1/4} $. Now the limit for $ N\varepsilon ^{1/4}\to 0 $ is of course $ {\cal M} = 3 $. The figure shows a curious phenomenon, that we did not try to investigate: before growing, here too as $ \varepsilon ^{1/4} $, $ {\cal M} $ decreases from $ 3 $ to about $ 2 $, and this for any $ N $. The accuracy of the scaling law is not as good as for a single wave, but remains satisfactory.

    Figure 11.  Left: same as Figure 8, for three traveling waves with $ k = 1, \, 2, \, 3 $. Right: same as Figure 9-left.

    Considering long waves with small $ k $ is essential to produce a FPU state obeying the scaling law (16). If instead we consider a single short wave excitation with, for example, $ k = 0.4N $ (recall the maximal frequency is at $ k = N/2 $), then the limit $ {\cal M}\to 1 $ for $ \varepsilon \to0 $ looks uniform in $ N $, and up to rather large $ \varepsilon \simeq0.1 $, see Figure 12, only one action is appreciably different from zero. Note that the vertical scale is quite enlarged with respect to Figure 8: passing from $ {\cal M} = 1 $ to $ 2 $ requires $ \varepsilon \simeq1 $.

    Figure 12.  Same as Figure 8, for a single traveling wave with $ k = 0.4N $.

    We feel indebted to Martino Cecchetto (Padova) for helping us in understanding the computation of the Toda actions.

    The authors declare no conflict of interest.

    We recall here the procedure to compute the actions proposed in [9]. Let us stress that the procedure is very clear in principle, but not so easy to implement in practice for large $ N $. To understand, ref. [35] was very helpful to us.

    A. The theoretical frame. Consider a periodic Toda chain with $ N $ particles, assume the barycenter is at rest. Denote

    $ a_i = e^{\lambda(q_{i+1}-q_i)/2}\ , \qquad b_i = \lambda p_i\ , $

    so that

    $ \prod_i a_i = 1\ ,\qquad \sum_i b_i = 0 \ , $

    and let

    $ L^\pm = (b1a1±aNa1b2a2a2b3a3bN1aN1±aNaN1bN)
    \qquad ; $

    $ L^+ $ is the Lax matrix associated to the system, used to prove integrability in [8,34], $ L^- $ is an auxiliary matrix. Following [9], consider the characteristic polynomials

    $ P^\pm(x) = \det(x\,I-L^\pm) $

    of $ L^\pm $, which turn out to satisfy, for any $ x $,

    $ P^-(x)-P^+(x) = 4\ . $

    Though the matrices change during the evolution, the characteristic polynomials do not (this is indeed the essence of Lax method to prove integrability). In particular, the eigenvalues $ \lambda^\pm_1, \ldots, \lambda^\pm_N $ of $ L^\pm $ stay constant, and for either choice of the sign they provide an independent set of $ N $ constants of motion. The eigenvalues will be supposed to be ordered decreasingly.

    To define the actions it is convenient to introduce the discriminant $ \Delta $, defined as

    $ \Delta(x) = P^+(x)+2 = P^-(x)-2\ , $

    which then satisfies

    $ \Delta(\lambda^+_i) = 2 , \qquad \Delta(\lambda^-_i) = -2\ . $

    We recall here some basic properties of $ \Delta $, for details and discussion see [9]. To be definite we assume $ N $ is even.

    (a) At mechanical equilibrium ($ a_i = 1 $, $ b_i = 0 $), the shape of $ \Delta $ is as in Figure 13, left panel, which refers to $ N = 16 $; $ \Delta $ oscillates exactly between $ \pm 2 $ and all eigenvalues $ \lambda^\pm_i $, but the extremal ones $ \lambda^+_1 $ and $ \lambda^+_N $, are double. All eigenvalues stay between|| $ \pm2 $.

    Figure 13.  Left: the discriminant for $ N = 16 $, at rest. Right: the discriminant (blue) and the arcs $ \pm(2+\rho) $ (red), for $ \varepsilon = 10^{-2} $; several traveling waves are present.

    ||Between $ \pm1 $ in [9]; for a certain convenience, our notation differs from [9] by a few sparse factors $ 2 $.

    $ (b) $ For nonvanishing energy, one or more pairs of degenerate eigenvalues $ \lambda^\pm_k, \lambda^\pm_{k+1} $ open, producing a "gap", and correspondingly there are maxima of $ \Delta $ exceeding $ 2 $ or minima exceeding $ -2 $. The shape of $ \Delta $ looks as in the right panel of Figure 13, blue curve. Denote by $ J_1, \ldots, J_{N-1} $, right to left, the intervals between nearby eigenvalues:

    $ {J_k} = {\mkern 1mu} \left\{ {\begin{array}{*{20}{l}} \begin{array}{l} {\mkern 1mu} (\lambda _{k + 1}^ - ,\lambda _k^ - )\\ (\lambda _{k + 1}^ + ,\lambda _k^ + ) \end{array}
    &koddkeven
    &{\;.} \end{array}} \right. $

    To each nonvanishing interval $ J_k $ corresponds a positive action $ {\cal I}_k $, defined by

    $ {\cal I}_k = \frac{1}{\pi}\,\int_{J_k}\rho(x)\, \,\mathrm{d} x\ ,\qquad \rho(x) = {\rm{acosh}}\frac{|\Delta(x)|}{2}\ . $

    Following [9], one can trace inside each open gap an arc $ \rho'(x) = \pm (2+\rho(x)) $, see the red curves of Figure 13. The actions are precisely, up to a trivial factor, the areas between $ \rho'(x) $ and the lines $ \pm2 $. Notice that if $ |\Delta| $ is large, then $ \rho\simeq\log|\Delta| $.

    B. An easy example. To exemplify the computation in an easy case, consider a state with a single sinusoidal wave with not small $ k $. Figure 14 refers to $ N = 32 $, $ k = 12\simeq0.4N $, $ \varepsilon $ in the range $ 10^{-4}-10^{-1} $. Quite evidently, the action $ {\cal I}_{12} $ is dominant, although it is clear that by increasing $ \varepsilon $ other actions start to be appreciable. The peak of $ \Delta $, already for $ \varepsilon = 0.01 $, is large, approximately $ 130 $, and grows to approximately $ 4\times 10^6 $ for $ \varepsilon = 0.1 $.

    Figure 14.  The discriminant (blue) and the arcs $ \pm(2+\rho) $ (red) for $ N = 32 $, $ \varepsilon = 10^{-4} $, $ 10^{-3} $, $ 10^{-2} $, $ 10^{-1} $; single traveling wave, $ k = 12 $.

    By increasing $ N $ the peaks of $ \Delta $ grow explosively, and this is the main problem in the numerical computation of the actions. Table 1 shows the peaks for some values of $ N $ up to $ 1024 $, $ \varepsilon $ between $ 10^{-4} $ and $ 10^{-1} $. High peaks should not surprise: in a state with a single traveling wave it is $ I_k = E_k/\omega _k = \varepsilon N/\omega _k $; should it be $ {\cal I}_k\simeq I_k $, then (assuming the gap, i.e., the base of the lobe formed by $ \rho'(x) $, is large) $ \rho $ is of order $ \varepsilon N $ and $ \Delta $ is its exponential.**

    Table 1.  The maximum of $ |\Delta| $ for different values of $ N $ and $ \varepsilon $; single left traveling wave, $ k\simeq0.4N $.
    $ \varepsilon =10^{-4} $ $ 10^{-3} $ $ 10^{-2} $ $ 10^{-1} $
    $ N=32 $ $ 2.5 $ $ 4.9 $ $ 1.3\times10^2 $ $ 3.5\times10^6 $
    $ 64 $ $ 3.0 $ $ 21 $ $ 1.4\times10^4 $ $ 7.3\times10^{12} $
    $ 128 $ $ 6.9 $ $ 4.1\times10^2 $ $ 1.8\times10^8 $ $ 3.4\times 10^{25} $
    $ 256 $ $ 45 $ $ 1.5\times10^5 $ $ 3.2\times10^{16} $ $ 1.2\times10^{51} $
    $ 512 $ $ 2.0\times10^3 $ $ 2.9\times10^{10} $ $ 1.0\times10^{33} $ $ 1.4\times10^{102} $
    $ 1024 $ $ 4.1\times10^6 $ $ 8.3\times10^{20} $ $ 9.3\times10^{65} $ $ 1.2\times10^{204} $

     | Show Table
    DownLoad: CSV

    **In [9], when peaks as large as $ 10^{15} $ are met, the comment is: The value of the discriminant at its maximum $ (10^{15}) $ is entirely typical, and perhaps surprisingly large for someone who has done inverse spectral theory with paper and pencil only.

    C. Technical questions. Problems of theoretical precision, in presence of such large numbers, are serious. The results reported in Section 3 should be consider as highly tentative, subject to criticism and justified mainly by their overall coherence.

    All computations have been made in the so-called quadruple precision, that is representing numbers with approximately $ 33 $ decimal digits. This is fairly enough for small $ N $, but questions arise for, say, $ N = 1024 $ or larger, unless $ \varepsilon $ is very small.

    To compute the actions one needs to compute with the best possible precision the eigenvalues $ \lambda^\pm_i $. Since this is done only once per each state, computer time can be invested. The procedure we found convenient is:

    $ – $ Using standard routines, namely routines (in double precision, $ 16 $ decimal digits) designed for sparse symmetric matrices, to compute the eigenvalues of $ L^\pm $. This provides a first approximation of $ \lambda^\pm_i $, with errors of the order, say, $ 10^{-14} $.

    $ – $ Improving the precision by searching the eigenvalues as zeros of $ P^{\pm}(x) $, with the Newton method. Easy recursive expressions are available for the determinant of three-diagonal periodic matrices, that is for $ P^{\pm}(x) $, as well as for the derivative. A little problem arises when the zeros of $ P^{\pm}(x) $ are tangent or nearly tangent, but this is solved by approximating $ P^{\pm}(x) $ with a parabola, and moreover, the corresponding gaps $ J_k $ and actions $ {\cal I}_k $ in this case are small, so one should not worry too much. The serious numerical problem instead is the usual loss of precision, when a zero results from a huge sum of large numbers. Here the only way out seems to us to be the use of multiprecision arithmetics, as we plan to do in a nearby future.

    Once the $ \lambda^{\pm}_i $ are obtained, the characteristic polynomials can be computed in the standard way,

    $ P^\pm(x) = { \prod_i}\,(x-\lambda^\pm_i)\ ; $

    $ \Delta(x) $, $ \rho(x) $ and the actions are immediately obtained. When $ |P^\pm(x)| $ starts to be large, one can directly compute $ \rho(x)\simeq\log|\Delta(x)|\simeq\log|P^\pm(x)| $,

    $ \rho(x)\simeq { \sum_i}\,\log|x-\lambda^\pm_i|\ . $

    Everything was somehow unusual in our experience, and a little painful.

    [1] Armstrong GL, Conn LA, Pinner RW (1999) Trends in infectious disease mortality in the United States during the 20th century. JAMA 281: 61-66.
    [2] Hamlin C, Sheard S (1998) Revolutions in public health: 1848, and 1998? BMJ 317: 587-591.
    [3] Leading Causes of Death, 1900-1998.
    [4] Breslow L (1999) From disease prevention to health promotion. JAMA 281: 1030-1033. doi: 10.1001/jama.281.11.1030
    [5] Institute of Medicine. Committee on Quality Care in America. (2001) CROSSING THE QUALITY CHASM: A New Health System for the 21st Century. Washington, DC: Institute of Medicine.
    [6] Winslow CE (1920) The Untilled Fields of Public Health. Science 51: 23-33.
    [7] Centers for Disease C, Prevention. (2011) Ten great public health achievements--United States, 2001-2010. MMWR. Morbidity and mortality weekly report 60: 619-623.
    [8] Centers for Disease C, Prevention. (1999) Ten great public health achievements--United States, 1900-1999. MMWR. Morbidity and mortality weekly report 48: 241-243.
    [9] Arias E (2015) United States Life Tables, 2011. National vital statistics reports : from the Centers for Disease Control and Prevention, National Center for Health Statistics, National Vital Statistics System. 64: 1-63.
    [10] Olshansky SJ, Passaro DJ, Hershow RC, et al. (2005) A potential decline in life expectancy in the United States in the 21st century. The New England Journal of Medicine 352: 1138-1145.
    [11] Oeppen J, Vaupel JW (2002) Demography. Broken limits to life expectancy. Science 296: 1029-1031.
    [12] Medicine IIo. A Nationwide Framework for Surveillance of Cardiovascular and Chronic Lung Disease. Washington, DC2011.
    [13] Nelson DE, Holtzman D, Bolen J, et al. (2001) Reliability and validity of measures from the Behavioral Risk Factor Surveillance System (BRFSS). Sozial- und Praventivmedizin 46: S3-42.
    [14] Stein AD, Lederman RI, Shea S (1993) The Behavioral Risk Factor Surveillance System questionnaire: its reliability in a statewide sample. Am J Public Health 83: 1768-1772.
    [15] Li C, Balluz LS, Ford ES, et al. (2012) A comparison of prevalence estimates for selected health indicators and chronic diseases or conditions from the Behavioral Risk Factor Surveillance System, the National Health Interview Survey, and the National Health and Nutrition Examination Survey, 2007-2008. Prev Med 54: 381-387.
    [16] Arday DR, Tomar SL, Nelson DE, et al. (1987) State smoking prevalence estimates: a comparison of the Behavioral Risk Factor Surveillance System and current population surveys. Am J Public Health 87: 1665-1669.
    [17] Hu SS, Balluz L, Battaglia MP, et al. (2011) Improving public health surveillance using a
    dual-frame survey of landline and cell phone numbers. Am J Epidemiol 173: 703-711.
    [18] A Nationwide Framework for Surveillance of Cardiovascular and Chronic Lung Diseases. Washington (DC) 2011.
    [19] Carlson SA, Densmore D, Fulton JE, et al. (2009) Differences in physical activity prevalence and trends from 3 U.S. surveillance systems: NHIS, NHANES, and BRFSS. J Physi Act Health 6: S18-27.
    [20] Aging. USDoHaHSAo. Aging statistics. http://www.aoa.acl.gov/aging_statistics/index.aspx. Accessed January 13, 2016.
    [21] Nelson ME, Rejeski WJ, Blair SN, et al. (2007) Physical activity and public health in older adults: recommendation from the American College of Sports Medicine and the American Heart Association. Circulation 116: 1094-1105. doi: 10.1161/CIRCULATIONAHA.107.185650
    [22] 2008 Physical Activity Guidelines for Americans. US Department of Health and Human Services; 2008.
    [23] Larson EB, Wang L, Bowen JD, et al. (2006) Exercise is associated with reduced risk for incident dementia among persons 65 years of age and older. An Int Med 144: 73-81.
    [24] DiPietro L (2001) Physical activity in aging: changes in patterns and their relationship to health and function. The journals of gerontology. Series A, Biological sciences and medical sciences. 56: 13-22.
    [25] Jones DW, Peterson ED, Bonow RO, et al. (2008) Translating research into practice for healthcare providers: the American Heart Association's strategy for building healthier lives, free of cardiovascular diseases and stroke. Circulation 118: 687-696. doi: 10.1161/CIRCULATIONAHA.108.189934
    [26] Clays E, Lidegaard M, De Bacquer D, et al. (2014) The combined relationship of occupational and leisure-time physical activity with all-cause mortality among men, accounting for physical fitness. Am J Epidemiol 179: 559-566.
    [27] Kampert JB, Blair SN, Barlow CE, et al. (1996) Physical activity, physical fitness, and all-cause and cancer mortality: a prospective study of men and women. Ann Epidemiol 6: 452-457.
    [28] Blair SN, Kohl HW, Paffenbarger RS, et al. (1989) Physical fitness and all-cause mortality. A prospective study of healthy men and women. JAMA 262: 2395-2401.
    [29] Smith TC, Wingard DL, Smith B, et al. (2007) Walking decreased risk of cardiovascular disease mortality in older adults with diabetes. J Clin Epidemiol 60: 309-317.
    [30] Prohaska T, Belansky E, Belza B, et al. (2006) Physical activity, public health, and aging: critical issues and research priorities. The journals of gerontology. Series B, Psychological sciences and social sciences 61: S267-273.
    [31] Slingerland AS, van Lenthe FJ, Jukema JW, et al. (2007) Aging, retirement, and changes in physical activity: prospective cohort findings from the GLOBE study. Am J Epidemiol 165 :1356-1363.
    [32] Littman A, Jacobson IG, Boyko EJ, et al. (2015) Changes in Meeting Vigorous Physical Activity Guidelines After Discharge From the Military. J Physi Act Health 12: 666-674.
    [33] Reilly T, Waterhouse J, Atkinson G. (1997) Aging, rhythms of physical performance, and adjustment to changes in the sleep-activity cycle. Occu Envir Med 54: 812-816.
    [34] Reuter I (2012) Aging, physical activity, and disease prevention. J Aging Res 2012: 373294.
    [35] Martin LG, Schoeni RF, Andreski PM (2010) Trends in health of older adults in the United States: past, present, future. Demography 47: S17-40.
    [36] Seeman TE, Merkin SS, Crimmins EM, et al. (2010) Disability trends among older Americans: National Health And Nutrition Examination Surveys, 1988-1994 and 1999-2004. Am J Public Health 100: 100-107. doi: 10.2105/AJPH.2008.157388
    [37] Lim K, Taylor L (2005) Factors associated with physical activity among older people—a population-based study. Prev Med 40: 33-40.
    [38] Florindo AA, Guimaraes VV, Cesar CL, et al. (2009) Epidemiology of leisure, transportation, occupational, and household physical activity: prevalence and associated factors. J Physi Act Health 6: 625-632.
    [39] Todt K, Skargren E, Jakobsson P, et al. (2015) Factors associated with low physical activity in patients with chronic obstructive pulmonary disease: a cross-sectional study. Scandinavian J Caring Scie 29: 697-707.
    [40] Healthy People 2020 Washington, DC: U.S. Department of Health and Human Services, Office of Disease Prevention and Health Promotion.
    [41] Piane GM, Smith TC (2014) Building an evidence base for the co-occurrence of chronic disease and psychiatric distress and impairment. Prevent Chronic Dis11: E188.
    [42] Kessler RC, Chiu WT, Demler O, et al. (2005) Prevalence, severity, and comorbidity of
    12-month DSM-IV disorders in the National Comorbidity Survey Replication. Arch General Psychiatry 62: 617-627.
    [43] Kessler RC, Berglund P, Demler O, et al. (2005) Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the National Comorbidity Survey Replication. Arch General Psychiatry 62: 593-602.
    [44] Karel MJ, Gatz M, Smyer MA (2012) Aging and mental health in the decade ahead: what psychologists need to know. The Am Psycholo 67: 184-198.
    [45] Institute of Medicine (IOM). The Mental Health and Substance Use Workforce for Older Adults: In Whose Hands? 2012.
    [46] U.S. Department of Health and Human Services. Healthy People 2010: Understanding and Improving Health. 2nd ed. Vol 2005. Washington, DC: U.S. Government Printing Office; November 2000.
    [47] Koh HK, Piotrowski JJ, Kumanyika S, et al. (2011) Healthy people: a 2020 vision for the social determinants approach. Health education & behavior: the official publication of the Society for Public Health Education 38: 551-557.
    [48] Koh HK (2010) A 2020 vision for healthy people. Eng J Med 362: 1653-1656.
    [49] Checkoway H, Pearce N, Kriebel D (2007) Selecting appropriate study designs to address specific research questions in occupational epidemiology. Occu Envir Med 64: 633-638.
    [50] Pearce N (2012) Classification of epidemiological study designs. Int J Epide 41: 393-397.
    [51] Rothman K, Greenland S (1998) Modern Epidemiology. Second ed. Philadelphia, PA: Lippincott-Raven.
    [52] Bowling A (2005) Mode of questionnaire administration can have serious effects on data quality. J Pub Health 27: 281-291. doi: 10.1093/pubmed/fdi031
    [53] Krebs NF, Himes JH, Jacobson D, et al. (2007) Assessment of child and adolescent overweight and obesity. Pediatrics 120: S193-228.
  • This article has been cited by:

    1. Chiara Caracciolo, Ugo Locatelli, Elliptic tori in FPU non-linear chains with a small number of nodes, 2021, 97, 10075704, 105759, 10.1016/j.cnsns.2021.105759
    2. Matteo Gallone, Antonio Ponno, Bob Rink, Korteweg–de Vries and Fermi–Pasta–Ulam–Tsingou: asymptotic integrability of quasi unidirectional waves, 2021, 54, 1751-8113, 305701, 10.1088/1751-8121/ac0a2e
    3. Wilberd van der Kallen, Ferdinand Verhulst, Explorations for alternating FPU-chains with large mass, 2021, 137, 00207462, 103809, 10.1016/j.ijnonlinmec.2021.103809
    4. Anna Vainchtein, Solitary waves in FPU-type lattices, 2022, 434, 01672789, 133252, 10.1016/j.physd.2022.133252
    5. Giancarlo Benettin, Antonio Ponno, 2023, Chapter 3, 978-981-19-6461-9, 21, 10.1007/978-981-19-6462-6_3
    6. M. Onorato, Y.V. Lvov, G. Dematteis, S. Chibbaro, Wave Turbulence and thermalization in one-dimensional chains, 2023, 1040, 03701573, 1, 10.1016/j.physrep.2023.09.006
    7. Giancarlo Benettin, Giuseppe Orsatti, Antonio Ponno, On the Role of the Integrable Toda Model in One-Dimensional Molecular Dynamics, 2023, 190, 1572-9613, 10.1007/s10955-023-03147-x
  • Reader Comments
  • © 2016 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6400) PDF downloads(1498) Cited by(1)

Figures and Tables

Figures(1)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog