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

Survival of Terrestrial Planets in the Presence of Giant Planet Migration

The presence of ``Hot Jupiters'', Jovian mass planets with very short orbital periods orbiting nearby main sequence stars, has been proposed to be primarily due to the orbital migration of planets formed in orbits initially much further from the parent star. The migration of giant planets would have profound effects on the evolution of inner terrestrial planets in these systems, and previous analyses have assumed that no terrestrial planets survive after migration has occurred. We present numerical simulations showing that a significant fraction of terrestrial planets could survive the migration process, eventually returning to circular orbits relatively close to their original positions. A fraction of the final orbits are in the Habitable Zone, suggesting that planetary systems with close-in giant planets are viable targets for searches for Earth-like habitable planets around other stars.Comment: 5 pages, 3 figures, emulateapj. ApJL in press, referee comments changes and edited for lengt

On data skewness, stragglers, and MapReduce progress indicators

We tackle the problem of predicting the performance of MapReduce applications, designing accurate progress indicators that keep programmers informed on the percentage of completed computation time during the execution of a job. Through extensive experiments, we show that state-of-the-art progress indicators (including the one provided by Hadoop) can be seriously harmed by data skewness, load unbalancing, and straggling tasks. This is mainly due to their implicit assumption that the running time depends linearly on the input size. We thus design a novel profile-guided progress indicator, called NearestFit, that operates without the linear hypothesis assumption and exploits a careful combination of nearest neighbor regression and statistical curve fitting techniques. Our theoretical progress model requires fine-grained profile data, that can be very difficult to manage in practice. To overcome this issue, we resort to computing accurate approximations for some of the quantities used in our model through space- and time-efficient data streaming algorithms. We implemented NearestFit on top of Hadoop 2.6.0. An extensive empirical assessment over the Amazon EC2 platform on a variety of real-world benchmarks shows that NearestFit is practical w.r.t. space and time overheads and that its accuracy is generally very good, even in scenarios where competitors incur non-negligible errors and wide prediction fluctuations. Overall, NearestFit significantly improves the current state-of-art on progress analysis for MapReduce

Adoption and Foster Care by Gay and Lesbian Parents in the United States

Discussion and debate about adoption and foster care by gay, lesbian, and bisexual (GLB) parents occurs frequently among child welfare policymakers, social service agencies, and social workers. They all need better information about GLB adoptive and foster parents and their children as they make individual and policy-level decisions about placement of children with GLB parents. This report provides new information on GLB adoption and foster care from the U.S. Census 2000, the National Survey of Family Growth (2002), and the Adoption and Foster Care Analysis and Reporting System (2004)

Wind-Tunnel Investigation of an Advanced General Aviation Canard Configuration

Wind-tunnel tests of a model of an advanced canard configuration designed for general aviation were conducted in the Langley 30- by 60-Foot Tunnel. The objective of the tests was to determine the aerodynamic stability and control characteristics of the configuration for a large range of angles of attack and sideslip at several power conditions. Analysis of the aerodynamic data indicates significant effects of power and of center-of-gravity location. For forward center-of-gravity locations, the configuration had extremely stall-resistant stability and control characteristics. For aft center-of-gravity locations and high-power conditions, the combined effects of increased pitch control and reduced longitudinal stability overpowered the stall resistance provided by the canard, which led to a high-angle-of-attack, deep-stall trim condition. Other aspects of the aerodynamic characteristics studied include the following: flow-visualization study, effect of negative angles of attack, lateral-directional characteristics, and comparison of the stall characteristics with another canard configuration

The impact of risk management practice upon the implementation of recovery-oriented care in community mental health services: a qualitative investigation

Background: Recovery-oriented care has become guiding principle for mental health policies and practice in the UK and elsewhere. However, a pre-existing culture of risk management practice may impact upon the provision of recovery-oriented mental health services. Aims: To explore how risk management practice impacts upon the implementation of recovery-oriented care within community mental health services. Method: Semi-structured interviews using vignettes were conducted with 8 mental health worker and service user dyads. Grounded theory techniques were used to develop explanatory themes. Results: Four themes arose: 1) recovery and positive risk taking; 2) competing frameworks of practice; 3) a hybrid of risk and recovery; 4) real-life recovery in the context of risk. Discussion: In abstract responses to the vignettes, mental health workers described how they would use a positive risk taking approach in support of recovery. In practice, this was restricted by a risk-averse culture embedded within services. Mental health workers set conditions with which service users complied to gain some responsibility for recovery. Conclusion: A lack of strategic guidance at policy level and lack of support and guidance at practice level may result in resistance to implementing ROC in the context of RMP. Recommendations are made for policy, training and future research

Testing for seasonal unit roots by frequency domain regression

This paper develops univariate seasonal unit root tests based on spectral regression estimators. An advantage of the frequency domain approach is that it enables serial correlation to be treated non-parametrically. We demonstrate that our proposed statistics have pivotal limiting distributions under both the null and near seasonally integrated alternatives when we allow for weak dependence in the driving shocks. This is in contrast to the popular seasonal unit root tests of, among others, Hylleberg et al. (1990) which treat serial correlation parametrically via lag augmentation of the test regression. Moreover, our analysis allows for (possibly infinite order) moving average behaviour in the shocks, while extant large sample results pertaining to the Hylleberg et al. (1990) type tests are based on the assumption of a finite autoregression. The size and power properties of our proposed frequency domain regression-based tests are explored and compared for the case of quarterly data with those of the tests of Hylleberg et al. (1990) in simulation experiments.Seasonal unit root tests; moving average; frequency domain regression; spectral density estimator; Brownian motion

The (In)Stability of Planetary Systems

We present results of numerical simulations which examine the dynamical stability of known planetary systems, a star with two or more planets. First we vary the initial conditions of each system based on observational data. We then determine regions of phase space which produce stable planetary configurations. For each system we perform 1000 ~1 million year integrations. We examine upsilon And, HD83443, GJ876, HD82943, 47UMa, HD168443, and the solar system (SS). We find that the resonant systems, 2 planets in a first order mean motion resonance, (HD82943 and GJ876) have very narrow zones of stability. The interacting systems, not in first order resonance, but able to perturb each other (upsilon And, 47UMa, and SS) have broad regions of stability. The separated systems, 2 planets beyond 10:1 resonance, (we only examine HD83443 and HD168443) are fully stable. Furthermore we find that the best fits to the interacting and resonant systems place them very close to unstable regions. The boundary in phase space between stability and instability depends strongly on the eccentricities, and (if applicable) the proximity of the system to perfect resonance. In addition to million year integrations, we also examined stability on ~100 million year timescales. For each system we ran ~10 long term simulations, and find that the Keplerian fits to these systems all contain configurations which may be regular on this timescale.Comment: 37 pages, 49 figures, 13 tables, submitted to Ap

Crystallization of the oligopeptide-binding protein AppA from Bacillus subtilis

AppA is the membrane-anchored extracellular receptor component of an ABC transporter responsible for the uptake of oligopeptides into Bacillus subtilis. AppA has been overexpressed as a cleavable maltose-binding protein fusion in Escherichia coli. Following removal of the fusion portion, AppA has been crystallized from morpholino-ethanesulfonic acid-buffered solutions at pH 6.5 containing polyethylene glycol and zinc acetate. A complete X-ray diffraction data set extending to 2.3 Angstrom spacing has been collected