The Media, Accountability and Civic Engagement in Africa

human development, democracy

Variational Integrators for the Gravitational N-Body Problem

This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the algorithm to handle individual time steps; this introduces fifth-order errors in angular momentum conservation and symplecticity. We show that using adaptive block power of two timesteps does not increase the error in symplecticity. In contrast to other high-order, symplectic, individual timestep, momentum-preserving algorithms, the algorithm takes only forward timesteps. We compare a code integrating an N-body system using the algorithm with a direct-summation force calculation to standard stellar cluster simulation codes. We find that our algorithm has about 1.5 orders of magnitude better symplecticity and momentum conservation errors than standard algorithms for equivalent numbers of force evaluations and equivalent energy conservation errors.Comment: 31 pages, 8 figures. v2: Revised individual-timestepping description, expanded comparison with other methods, corrected error in predictor equation. ApJ, in pres

The WISDOM Study: breaking the deadlock in the breast cancer screening debate.

There are few medical issues that have generated as much controversy as screening for breast cancer. In science, controversy often stimulates innovation; however, the intensely divisive debate over mammographic screening has had the opposite effect and has stifled progress. The same two questions-whether it is better to screen annually or bi-annually, and whether women are best served by beginning screening at 40 or some later age-have been debated for 20 years, based on data generated three to four decades ago. The controversy has continued largely because our current approach to screening assumes all women have the same risk for the same type of breast cancer. In fact, we now know that cancers vary tremendously in terms of timing of onset, rate of growth, and probability of metastasis. In an era of personalized medicine, we have the opportunity to investigate tailored screening based on a woman's specific risk for a specific tumor type, generating new data that can inform best practices rather than to continue the rancorous debate. It is time to move from debate to wisdom by asking new questions and generating new knowledge. The WISDOM Study (Women Informed to Screen Depending On Measures of risk) is a pragmatic, adaptive, randomized clinical trial comparing a comprehensive risk-based, or personalized approach to traditional annual breast cancer screening. The multicenter trial will enroll 100,000 women, powered for a primary endpoint of non-inferiority with respect to the number of late stage cancers detected. The trial will determine whether screening based on personalized risk is as safe, less morbid, preferred by women, will facilitate prevention for those most likely to benefit, and adapt as we learn who is at risk for what kind of cancer. Funded by the Patient Centered Outcomes Research Institute, WISDOM is the product of a multi-year stakeholder engagement process that has brought together consumers, advocates, primary care physicians, specialists, policy makers, technology companies and payers to help break the deadlock in this debate and advance towards a new, dynamic approach to breast cancer screening

Variational Integrators for Almost-Integrable Systems

We construct several variational integrators--integrators based on a discrete variational principle--for systems with Lagrangians of the form L = L_A + epsilon L_B, with epsilon << 1, where L_A describes an integrable system. These integrators exploit that epsilon << 1 to increase their accuracy by constructing discrete Lagrangians based on the assumption that the integrator trajectory is close to that of the integrable system. Several of the integrators we present are equivalent to well-known symplectic integrators for the equivalent perturbed Hamiltonian systems, but their construction and error analysis is significantly simpler in the variational framework. One novel method we present, involving a weighted time-averaging of the perturbing terms, removes all errors from the integration at O(epsilon). This last method is implicit, and involves evaluating a potentially expensive time-integral, but for some systems and some error tolerances it can significantly outperform traditional simulation methods.Comment: 14 pages, 4 figures. Version 2: added informative example; as accepted by Celestial Mechanics and Dynamical Astronom

The role of chaotic resonances in the solar system

Our understanding of the Solar System has been revolutionized over the past decade by the finding that the orbits of the planets are inherently chaotic. In extreme cases, chaotic motions can change the relative positions of the planets around stars, and even eject a planet from a system. Moreover, the spin axis of a planet-Earth's spin axis regulates our seasons-may evolve chaotically, with adverse effects on the climates of otherwise biologically interesting planets. Some of the recently discovered extrasolar planetary systems contain multiple planets, and it is likely that some of these are chaotic as well.Comment: 28 pages, 9 figure

Forward Symplectic Integrators and the Long Time Phase Error in Periodic Motions

We show that when time-reversible symplectic algorithms are used to solve periodic motions, the energy error after one period is generally two orders higher than that of the algorithm. By use of correctable algorithms, we show that the phase error can also be eliminated two orders higher than that of the integrator. The use of fourth order forward time step integrators can result in sixth order accuracy for the phase error and eighth accuracy in the periodic energy. We study the 1-D harmonic oscillator and the 2-D Kepler problem in great details, and compare the effectiveness of some recent fourth order algorithms.Comment: Submitted to Phys. Rev. E, 29 Page

Sums and differences of four k-th powers

We prove an upper bound for the number of representations of a positive integer $N$ as the sum of four $k$-th powers of integers of size at most $B$, using a new version of the Determinant method developed by Heath-Brown, along with recent results by Salberger on the density of integral points on affine surfaces. More generally we consider representations by any integral diagonal form. The upper bound has the form $O_{N}(B^{c/\sqrt{k}})$, whereas earlier versions of the Determinant method would produce an exponent for $B$ of order $k^{-1/3}$ in this case. Furthermore, we prove that the number of representations of a positive integer $N$ as a sum of four $k$-th powers of non-negative integers is at most $O_{\epsilon}(N^{1/k+2/k^{3/2}+\epsilon})$ for $k \geq 3$, improving upon bounds by Wisdom.Comment: 18 pages. Mistake corrected in the statement of Theorem 1.2. To appear in Monatsh. Mat

Pseudo-High-Order Symplectic Integrators

Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly with order in timestep, rendering higher-order methods impractical. However, symplectic integrators are often applied to systems in which perturbations between bodies are a small factor of the force due to a dominant central mass. In this case, it is possible to create optimized symplectic algorithms that require fewer substeps per timestep. This is achieved by only considering error terms of order epsilon, and neglecting those of order epsilon^2, epsilon^3 etc. Here we devise symplectic algorithms with 4 and 6 substeps per step which effectively behave as 4th and 6th-order integrators when epsilon is small. These algorithms are more efficient than the usual 2nd and 4th-order methods when applied to planetary systems.Comment: 14 pages, 5 figures. Accepted for publication in the Astronomical Journa

A precise modeling of Phoebe's rotation

Although the rotation of some Saturn's satellites in spin-orbit has already been studied by several authors, this is not the case of the rotation of Phoebe, which has the particularity of being non resonant. The purpose of the paper is to determine for the first time and with precision its precession-nutation motion. We adopt an Hamiltonian formalism of the motion of rotation of rigid celestial body set up by Kinoshita (1977) based on Andoyer variables and canonical equations. First we calculate Phoebe's obliquity at J2000,0 from available astronomical data as well as the gravitational perturbation due to Saturn on Phoebe rotational motion. Then we carry out a numerical integration and we compare our results for the precession rate and the nutation coefficients with pure analytical model. Our results for Phoebe obliquity (23{\deg}95) and Phoebe precession rate (5580".65/cy) are very close to the respective values for the Earth. Moreover the amplitudes of the nutations (26" peak to peak for the nutaton in longitude and 8" for the nutation in obliquity) are of the same order as the respective amplitudes for the Earth. We give complete tables of nutation, obtained from a FFT analysis starting from the numerical signals. We show that a pure analytical model of the nutation is not accurate due to the fact that Phoebe orbital elements e, M and Ls are far from having a simple linear behaviour. The precession and nutation of Phoebe have been calculated for the first time in this paper. We should keep on the study in the future by studying the additional gravitational effects of the Sun, of the large satellites as Titan, as well as Saturn dynamical ellipticity.Comment: 11 pages,15 figures, accepted for publication in A&

Psychosocial Aspects of Physical Activity and Fitness In Special-Population, Minority Middle School Children

Special-population research predicting physical activity (PA) and fitness with minority middle school children from at-risk environments is rare. Hence, the purpose of our investigation was to evaluate the ability of important social cognitive and environment-based measures to predict PA and fitness with children with developmental delay, cognitive, and emotional impairments. Children (N = 89, ages 11-15) completed questionnaires assessing social cognitive and environment-based constructs, self report PA, and completed fitness testing. Correlational results supported some hypotheses. The descriptive and correlational results also indicated commonalities with similar research on non special-population minority middle school children from at-risk environments