Supporting the Perpetuation and Reproducibility of Numerical Method Publications

AbstractThe development of new numerical methods is of great importance in computational science. Due to their many appealing properties, Finite Element (FEMs), Finite Volume (FVMs) and Finite Difference (DFMs) methods are of particular interest, with a very large number of journal articles devoted to them. Unfortunately, these methods can be time consuming to implement when one would like to explore certain numerical properties or make comparisons with other methods. In this paper, we propose a new computational environment in which the specification and implementation of such methods introduced in journal articles can be perpetuated and the experimentation with them can be made reproducible and easily compared with related work

On the dynamics of Extrasolar Planetary Systems under dissipation. Migration of planets

We study the dynamics of planetary systems with two planets moving in the same plane, when frictional forces act on the two planets, in addition to the gravitational forces. The model of the general three-body problem is used. Different laws of friction are considered. The topology of the phase space is essential in understanding the evolution of the system. The topology is determined by the families of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, or close to them, that the planets migrate when dissipative forces act. At the critical points where the stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and the migration process changes route and follows the new stable family up to large eccentricities or to a chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 resonances are studied. The migration to larger or smaller eccentricities depends on the particular law of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For particular laws of friction and for special values of the parameters of the frictional forces, it is possible to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed.Comment: Accepted in Celestial Mechanics and Dynamical Astronom

The Relativistic Factor in the Orbital Dynamics of Point Masses

There is a growing population of relativistically relevant minor bodies in the Solar System and a growing population of massive extrasolar planets with orbits very close to the central star where relativistic effects should have some signature. Our purpose is to review how general relativity affects the orbital dynamics of the planetary systems and to define a suitable relativistic correction for Solar System orbital studies when only point masses are considered. Using relativistic formulae for the N body problem suited for a planetary system given in the literature we present a series of numerical orbital integrations designed to test the relevance of the effects due to the general theory of relativity in the case of our Solar System. Comparison between different algorithms for accounting for the relativistic corrections are performed. Relativistic effects generated by the Sun or by the central star are the most relevant ones and produce evident modifications in the secular dynamics of the inner Solar System. The Kozai mechanism, for example, is modified due to the relativistic effects on the argument of the perihelion. Relativistic effects generated by planets instead are of very low relevance but detectable in numerical simulations

The Science of Sungrazers, Sunskirters, and Other Near-Sun Comets

This review addresses our current understanding of comets that venture close to the Sun, and are hence exposed to much more extreme conditions than comets that are typically studied from Earth. The extreme solar heating and plasma environments that these objects encounter change many aspects of their behaviour, thus yielding valuable information on both the comets themselves that complements other data we have on primitive solar system bodies, as well as on the near-solar environment which they traverse. We propose clear definitions for these comets: We use the term near-Sun comets to encompass all objects that pass sunward of the perihelion distance of planet Mercury (0.307 AU). Sunskirters are defined as objects that pass within 33 solar radii of the Sun’s centre, equal to half of Mercury’s perihelion distance, and the commonly-used phrase sungrazers to be objects that reach perihelion within 3.45 solar radii, i.e. the fluid Roche limit. Finally, comets with orbits that intersect the solar photosphere are termed sundivers. We summarize past studies of these objects, as well as the instruments and facilities used to study them, including space-based platforms that have led to a recent revolution in the quantity and quality of relevant observations. Relevant comet populations are described, including the Kreutz, Marsden, Kracht, and Meyer groups, near-Sun asteroids, and a brief discussion of their origins. The importance of light curves and the clues they provide on cometary composition are emphasized, together with what information has been gleaned about nucleus parameters, including the sizes and masses of objects and their families, and their tensile strengths. The physical processes occurring at these objects are considered in some detail, including the disruption of nuclei, sublimation, and ionisation, and we consider the mass, momentum, and energy loss of comets in the corona and those that venture to lower altitudes. The different components of comae and tails are described, including dust, neutral and ionised gases, their chemical reactions, and their contributions to the near-Sun environment. Comet-solar wind interactions are discussed, including the use of comets as probes of solar wind and coronal conditions in their vicinities. We address the relevance of work on comets near the Sun to similar objects orbiting other stars, and conclude with a discussion of future directions for the field and the planned ground- and space-based facilities that will allow us to address those science topics

Angular momentum exchange during secular migration of two-planet systems

We investigate the secular dynamics of two-planet coplanar systems evolving under mutual gravitational interactions and dissipative forces. We consider two mechanisms responsible for the planetary migration: star-planet (or planet-satellite) tidal interactions and interactions of a planet with a gaseous disc. We show that each migration mechanism is characterized by a specific law of orbital angular momentum exchange. Calculating stationary solutions of the conservative secular problem and taking into account the orbital angular momentum leakage, we trace the evolutionary routes followed by the planet pairs during the migration process. This procedure allows us to recover the dynamical history of two-planet systems and constrain parameters of the involved physical processes.Comment: 20 pages, 9 figures, accepted for publication in Celestial Mechanics and Dynamical Astronomy (special issue on Exoplanets

Search for direct production of charginos and neutralinos in events with three leptons and missing transverse momentum in √s = 7 TeV pp collisions with the ATLAS detector

A search for the direct production of charginos and neutralinos in final states with three electrons or muons and missing transverse momentum is presented. The analysis is based on 4.7 fb−1 of proton–proton collision data delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in three signal regions that are either depleted or enriched in Z-boson decays. Upper limits at 95% confidence level are set in R-parity conserving phenomenological minimal supersymmetric models and in simplified models, significantly extending previous results

Jet size dependence of single jet suppression in lead-lead collisions at sqrt(s(NN)) = 2.76 TeV with the ATLAS detector at the LHC

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sensitivity to the physics of jet quenching. In a sample of lead-lead collisions at sqrt(s) = 2.76 TeV corresponding to an integrated luminosity of approximately 7 inverse microbarns, ATLAS has measured jets with a calorimeter over the pseudorapidity interval |eta| < 2.1 and over the transverse momentum range 38 < pT < 210 GeV. Jets were reconstructed using the anti-kt algorithm with values for the distance parameter that determines the nominal jet radius of R = 0.2, 0.3, 0.4 and 0.5. The centrality dependence of the jet yield is characterized by the jet "central-to-peripheral ratio," Rcp. Jet production is found to be suppressed by approximately a factor of two in the 10% most central collisions relative to peripheral collisions. Rcp varies smoothly with centrality as characterized by the number of participating nucleons. The observed suppression is only weakly dependent on jet radius and transverse momentum. These results provide the first direct measurement of inclusive jet suppression in heavy ion collisions and complement previous measurements of dijet transverse energy imbalance at the LHC.Comment: 15 pages plus author list (30 pages total), 8 figures, 2 tables, submitted to Physics Letters B. All figures including auxiliary figures are available at http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/HION-2011-02

Measurement of the polarisation of W bosons produced with large transverse momentum in pp collisions at sqrt(s) = 7 TeV with the ATLAS experiment

This paper describes an analysis of the angular distribution of W->enu and W->munu decays, using data from pp collisions at sqrt(s) = 7 TeV recorded with the ATLAS detector at the LHC in 2010, corresponding to an integrated luminosity of about 35 pb^-1. Using the decay lepton transverse momentum and the missing transverse energy, the W decay angular distribution projected onto the transverse plane is obtained and analysed in terms of helicity fractions f0, fL and fR over two ranges of W transverse momentum (ptw): 35 < ptw < 50 GeV and ptw > 50 GeV. Good agreement is found with theoretical predictions. For ptw > 50 GeV, the values of f0 and fL-fR, averaged over charge and lepton flavour, are measured to be : f0 = 0.127 +/- 0.030 +/- 0.108 and fL-fR = 0.252 +/- 0.017 +/- 0.030, where the first uncertainties are statistical, and the second include all systematic effects.Comment: 19 pages plus author list (34 pages total), 9 figures, 11 tables, revised author list, matches European Journal of Physics C versio