-
Octofitter: Fast, Flexible, and Accurate Orbit Modelling to Detect Exoplanets
Authors:
William Thompson,
Jensen Lawrence,
Dori Blakely,
Christian Marois,
Jason Wang,
Mosé Giordano,
Timothy Brandt,
Doug Johnstone,
Jean-Baptiste Ruffio,
S. Mark Ammons,
Katie A. Crotts,
Clarissa R. Do Ó,
Eileen C. Gonzales,
Malena Rice
Abstract:
As next-generation imaging instruments and interferometers search for planets closer to their stars, they must contend with increasing orbital motion and longer integration times. These compounding effects make it difficult to detect faint planets but also present an opportunity. Increased orbital motion makes it possible to move the search for planets into the orbital domain, where direct images…
▽ More
As next-generation imaging instruments and interferometers search for planets closer to their stars, they must contend with increasing orbital motion and longer integration times. These compounding effects make it difficult to detect faint planets but also present an opportunity. Increased orbital motion makes it possible to move the search for planets into the orbital domain, where direct images can be freely combined with the radial velocity and proper motion anomaly, even without a confirmed detection in any single epoch. In this paper, we present a fast and differentiable multimethod orbit-modeling and planet detection code called Octofitter. This code is designed to be highly modular and allows users to easily adjust priors, change parameterizations, and specify arbitrary function relations between the parameters of one or more planets. Octofitter further supplies tools for examining model outputs including prior and posterior predictive checks and simulation-based calibration. We demonstrate the capabilities of Octofitter on real and simulated data from different instruments and methods, including HD 91312, simulated JWST/NIRISS aperture masking interferometry observations, radial velocity curves, and grids of images from the Gemini Planet Imager. We show that Octofitter can reliably recover faint planets in long sequences of images with arbitrary orbital motion. This publicly available tool will enable the broad application of multiepoch and multimethod exoplanet detection, which could improve how future targeted ground- and space-based surveys are performed. Finally, its rapid convergence makes it a useful addition to the existing ecosystem of tools for modeling the orbits of directly imaged planets.
△ Less
Submitted 2 February, 2024;
originally announced February 2024.
-
Effects of the chromospheric Lyα line profile shape on the determination of the solar wind HI outflow velocity using the Doppler dimming technique
Authors:
G. E. Capuano,
S. Dolei,
D. Spadaro,
S. L. Guglielmino,
P. Romano,
R. Ventura,
V. Andretta,
A. Bemporad,
C. Sasso,
R. Susino,
V. Da Deppo,
F. Frassetto,
S. M. Giordano,
F. Landini,
G. Nicolini,
M. Pancrazzi,
M. Romoli,
L. Zangrilli
Abstract:
The determination of solar wind outflow velocity is fundamental in order to probe the mechanisms of wind acceleration in the corona. We aim to study, via the Doppler dimming technique, the effects that the chromospheric Lyα line profile shape causes on the determination of the outflow speed of coronal HI atoms. The Doppler dimming technique takes into account the decrease of coronal Lyα radiation…
▽ More
The determination of solar wind outflow velocity is fundamental in order to probe the mechanisms of wind acceleration in the corona. We aim to study, via the Doppler dimming technique, the effects that the chromospheric Lyα line profile shape causes on the determination of the outflow speed of coronal HI atoms. The Doppler dimming technique takes into account the decrease of coronal Lyα radiation in regions where HI atoms flow out in the solar wind. Starting from UV observations (UVCS/SOHO) of the coronal Lyα line and simultaneous measurements of pB (LASCO/SOHO and Mk3/MLSO), we studied the effect of the pumping chromospheric Lyα line profile through measurements from SOHO/SUMER, UVSP/SMM and LPSP/OSO-8, taken from representative on-disk regions and as a function of time during the solar activity cycle. In particular, we considered the effect of four chromospheric line parameters: line width, depth of the central reversal, asymmetry and distance of the peaks. We find that the range of variability of these parameters is of about 50% for the width, 69% for the depth of the central reversal, 35% for the asymmetry, and 50% for the distance of the peaks. Then, we find that the variability of the pumping Lyα profile affects the estimates of the coronal HI velocity by about 9-12%. Therefore, this uncertainty is smaller than other physical quantities uncertainties, and a constant in time and unique shape of the Lyα profile over the solar disk can be adopted in order to estimate the solar wind outflow velocity.
△ Less
Submitted 12 August, 2021;
originally announced August 2021.
-
Contribution of polar plumes to fast solar wind
Authors:
L. Zangrilli,
S. M. Giordano
Abstract:
Context. Several physical properties of solar polar plumes have been identified by different published studies, however such studies are rare and sometimes in disagreement. Aims. The purpose of the present work is to analyze a set of SOHO/UVCS data dedicated to the observation of plumes and to obtain a picture of the physical properties of plumes in the intermediate solar corona through a self-con…
▽ More
Context. Several physical properties of solar polar plumes have been identified by different published studies, however such studies are rare and sometimes in disagreement. Aims. The purpose of the present work is to analyze a set of SOHO/UVCS data dedicated to the observation of plumes and to obtain a picture of the physical properties of plumes in the intermediate solar corona through a self-consistent analysis. Methods. We applied the Doppler Dimming technique to data acquired by SOHO/UVCS in April 1996, which was during the very early phases of the mission. From this we derived outflow speeds and electron densities. We used SOHO/LASCO images as context data in order to better identify plume and interplume regions in the UVCS field of view. Results. The results we obtain demonstrate that in three cases out of four plumes expand with outflow speeds comparable to those of interplumes, and in a single case with lower speeds. We estimate that the contribution of plumes to the wind coming from the solar poles is about 20%, and that different plumes provide a different contribution, possibly according to different stages of their evolution. Conclusions. We conclude that plumes are not static structures, and that they contribute significantly to the wind coming from the solar poles.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Towards new solutions for scientific computing: the case of Julia
Authors:
Maurizio Tomasi,
Mosè Giordano
Abstract:
This year marks the consolidation of Julia (https://julialang.org/), a programming language designed for scientific computing, as the first stable version (1.0) has been released, in August 2018. Among its main features, expressiveness and high execution speeds are the most prominent: the performance of Julia code is similar to statically compiled languages, yet Julia provides a nice interactive s…
▽ More
This year marks the consolidation of Julia (https://julialang.org/), a programming language designed for scientific computing, as the first stable version (1.0) has been released, in August 2018. Among its main features, expressiveness and high execution speeds are the most prominent: the performance of Julia code is similar to statically compiled languages, yet Julia provides a nice interactive shell and fully supports Jupyter; moreover, it can transparently call external codes written in C, Fortran, and even Python and R without the need of wrappers. The usage of Julia in the astronomical community is growing, and a GitHub organization named JuliaAstro takes care of coordinating the development of packages. In this paper, we present the features and shortcomings of this language and discuss its application in astronomy and astrophysics.
△ Less
Submitted 4 December, 2018;
originally announced December 2018.
-
On the relevance of chaos for halo stars in the solar neighbourhood II
Authors:
Nicolas P. Maffione,
Facundo A. Gómez,
Pablo M. Cincotta,
Claudia M. Giordano,
Robert J. J. Grand,
Federico Marinacci,
Rüdiger Pakmor,
Christine M. Simpson,
Volker Springel,
Carlos S. Frenk
Abstract:
In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of c…
▽ More
In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of chaotic mixing using multicomponent Galactic potential models and solar neighbourhood-like volumes extracted from fully cosmological hydrodynamic simulations, thus naturally accounting for the gravitational potential associated with baryonic components, such as the bulge and disc. Despite the strong change in the global Galactic potentials with respect to those obtained in dark matter only simulations, our results confirm that a large fraction of halo particles evolving on chaotic orbits exhibit their chaotic behaviour after periods of time significantly larger than a Hubble time. In addition, significant diffusion in phase space is not observed on those particles that do exhibit chaotic behaviour within a Hubble time.
△ Less
Submitted 15 May, 2018; v1 submitted 11 January, 2018;
originally announced January 2018.
-
Astrometric microlensing
Authors:
A. A. Nucita,
F. De Paolis,
G. Ingrosso,
M. Giordano,
L. Manni
Abstract:
Astrometric microlensing will offer in the next future a new channel for investigating the nature of both lenses and sources involved in a gravitational microlensing event. The effect, corresponding to the shift of the position of the multiple image centroid with respect to the source star location, is expected to occurr on scales from micro-arcoseconds to milli-arcoseconds depending on the charac…
▽ More
Astrometric microlensing will offer in the next future a new channel for investigating the nature of both lenses and sources involved in a gravitational microlensing event. The effect, corresponding to the shift of the position of the multiple image centroid with respect to the source star location, is expected to occurr on scales from micro-arcoseconds to milli-arcoseconds depending on the characteristics of the lens-source system. Here, we consider different classes of events (single/binary lens acting on a single/binary source) also accounting for additional effects including the finite source size, the blending and orbital motion. This is particularly important in the era of Gaia observations which is making possible astrometric measurements with unprecedent quality.
△ Less
Submitted 4 May, 2017;
originally announced May 2017.
-
A catalogue sample of low mass galaxies observed in X-rays with central candidate black holes
Authors:
A. A. Nucita,
L. Manni,
F. De Paolis,
M. Giordano,
G. Ingrosso
Abstract:
We present a sample of $X$-ray selected candidate black holes in 51 low mass galaxies with $z\le 0.055$ {and mass up to $10^{10}$ M$_{\odot}$} obtained by cross-correlating the NASA-SLOAN Atlas with the 3XMM catalogue. {We have also searched in the available catalogues for radio counterparts of the black hole candidates and find that 19 of the previously selected sources have also a radio counterp…
▽ More
We present a sample of $X$-ray selected candidate black holes in 51 low mass galaxies with $z\le 0.055$ {and mass up to $10^{10}$ M$_{\odot}$} obtained by cross-correlating the NASA-SLOAN Atlas with the 3XMM catalogue. {We have also searched in the available catalogues for radio counterparts of the black hole candidates and find that 19 of the previously selected sources have also a radio counterpart.} Our results show that about $37\%$ of the galaxies of our sample host { an $X$-ray source} (associated to a radio counterpart) spatially coincident with the galaxy center, in agreement with { other recent works}. For these {\it nuclear} sources, the $X$-ray/radio fundamental plane relation allows one to estimate the mass of the (central) candidate black holes which results to be in the range $10^{4}-2\times10^{8}$ M$_{\odot}$ (with median value of $\simeq 3\times 10^7$ M$_{\odot}$ and eight candidates having mass below $10^{7}$ M$_{\odot}$). This result, while suggesting that $X$-ray emitting black holes in low-mass galaxies may have had a key role in the evolution of such systems, makes even more urgent to explain how such massive objects formed in galaxies. {Of course, dedicated follow-up observations both in the $X$-ray and radio bands, as well as in the optical, are necessary in order to confirm our results
△ Less
Submitted 11 April, 2017;
originally announced April 2017.
-
Timing analysis in microlensing
Authors:
Mosè Giordano,
Achille Nucita,
Francesco De Paolis,
Gabriele Ingrosso
Abstract:
Timing analysis is a powerful tool used to determine periodic features of physical phenomena. Here we review two applications of timing analysis to gravitational microlensing events. The first one, in particular cases, allows the estimation of the orbital period of binary lenses, which in turn enables the breaking of degeneracies. The second one is a method to measure the rotation period of the le…
▽ More
Timing analysis is a powerful tool used to determine periodic features of physical phenomena. Here we review two applications of timing analysis to gravitational microlensing events. The first one, in particular cases, allows the estimation of the orbital period of binary lenses, which in turn enables the breaking of degeneracies. The second one is a method to measure the rotation period of the lensed star by observing signatures due to stellar spots on its surface.
△ Less
Submitted 20 December, 2016;
originally announced December 2016.
-
Parallax and orbital effects in astrometric microlensing with binary sources
Authors:
A. A. Nucita,
F. De Paolis,
G. Ingrosso,
M. Giordano,
L. Manni
Abstract:
In gravitational microlensing, binary systems may act as lenses or sources. Identifying lens binarity is generally easy especially in events characterized by caustic crossing since the resulting light curve exhibits strong deviations from smooth single-lensing light curve. On the contrary, light curves with minor deviations from a Paczyński behaviour do not allow one to identify the source binarit…
▽ More
In gravitational microlensing, binary systems may act as lenses or sources. Identifying lens binarity is generally easy especially in events characterized by caustic crossing since the resulting light curve exhibits strong deviations from smooth single-lensing light curve. On the contrary, light curves with minor deviations from a Paczyński behaviour do not allow one to identify the source binarity. A consequence of the gravitational microlensing is the shift of the position of the multiple image centroid with respect to the source star location - the so called astrometric microlensing signal. When the astrometric signal is considered, the presence of a binary source manifests with a path that largely differs from that expected for single-source events. Here, we investigate the astrometric signatures of binary sources taking into account their orbital motion and the parallax effect due to the Earth motion, which turn out not to be negligible in most cases. We also show that considering the above-mentioned effects is important in the analysis of astrometric data in order to correctly estimate the lens-event parameters.
△ Less
Submitted 7 June, 2016;
originally announced June 2016.
-
An optical transmission spectrum of the transiting hot Jupiter in the metal-poor WASP-98 planetary system
Authors:
L. Mancini,
M. Giordano,
P. Molliere,
J. Southworth,
R. Brahm,
S. Ciceri,
Th. Henning
Abstract:
The WASP-98 planetary system represents a rare case of a hot Jupiter hosted by a metal-poor main-sequence star. We present a follow-up study of this system based on multi-band photometry and high-resolution spectroscopy. Two new transit events of WASP-98b were simultaneously observed in four passbands (g,r,i,z), using the telescope-defocussing technique, yielding eight high-precision light curves…
▽ More
The WASP-98 planetary system represents a rare case of a hot Jupiter hosted by a metal-poor main-sequence star. We present a follow-up study of this system based on multi-band photometry and high-resolution spectroscopy. Two new transit events of WASP-98b were simultaneously observed in four passbands (g,r,i,z), using the telescope-defocussing technique, yielding eight high-precision light curves with point-to-point scatters of less than 1 mmag. We also collected three spectra of the parent star with a high-resolution spectrograph, which we used to remeasure its spectral characteristics, in particular its metallicity. We found this to be very low, Fe/H]=-0.49, but larger than was previously reported, [Fe/H]=-0.60. We used these new photometric and spectroscopic data to refine the orbital and physical properties of this planetary system, finding that the stellar and planetary mass measurements are significantly larger than those in the discovery paper. In addition, the multi-band light curves were used to construct an optical transmission spectrum of WASP-98b and probe the characteristics of its atmosphere at the terminator. We measured a lower radius at z compared with the other three passbands. The maximum variation is between the r and z bands, has a confidence level of roughly 6 sigma and equates to 5.5 pressure scale heights. We compared this spectrum to theoretical models, investigating several possible types of atmospheres, including hazy, cloudy, cloud-free, and clear atmospheres with titanium and vanadium oxide opacities. We could not find a good fit to the observations, except in the extreme case of a clear atmosphere with TiO and VO opacities, in which the condensation of Ti and V was suppressed. As this case is unrealistic, our results suggest the presence of an additional optical-absorbing species in the atmosphere of WASP-98b, of unknown chemical nature.
△ Less
Submitted 1 June, 2016;
originally announced June 2016.
-
The Scales of Gravitational Lensing
Authors:
Francesco De Paolis,
Mosè Giordano,
Gabriele Ingrosso,
Luigi Manni,
Achille Nucita,
Francesco Strafella
Abstract:
After exactly a century since the formulation of the general theory of relativity, the phenomenon of gravitational lensing is still an extremely powerful method for investigating in astrophysics and cosmology. Indeed, it is adopted to study the distribution of the stellar component in the Milky Way, to study dark matter and dark energy on very large scales and even to discover exoplanets. Moreover…
▽ More
After exactly a century since the formulation of the general theory of relativity, the phenomenon of gravitational lensing is still an extremely powerful method for investigating in astrophysics and cosmology. Indeed, it is adopted to study the distribution of the stellar component in the Milky Way, to study dark matter and dark energy on very large scales and even to discover exoplanets. Moreover, thanks to technological developments, it will allow the measure of the physical parameters (mass, angular momentum and electric charge) of supermassive black holes in the center of ours and nearby galaxies.
△ Less
Submitted 22 April, 2016;
originally announced April 2016.
-
First X-ray dectection of the young variable V1180 Cas
Authors:
S. Antoniucci,
A. A. Nucita,
T. Giannini,
D. Lorenzetti,
B. Stelzer,
D. Gerardi,
S. Delle Rose,
A. Di Paola,
M. Giordano,
L. Manni,
F. Strafella
Abstract:
V1180 Cas is a young variable that has shown strong photometric fluctuations (Delta_I~6mag) in the recent past, which have been attributed to events of enhanced accretion. The source has entered a new high-brightness state in Sept.2013, which we have previously analyzed through optical and near-IR spectroscopy. To investigate the current active phase of V1180 Cas, we performed observations with th…
▽ More
V1180 Cas is a young variable that has shown strong photometric fluctuations (Delta_I~6mag) in the recent past, which have been attributed to events of enhanced accretion. The source has entered a new high-brightness state in Sept.2013, which we have previously analyzed through optical and near-IR spectroscopy. To investigate the current active phase of V1180 Cas, we performed observations with the Chandra satellite to study the X-ray emission from the object and its connection to accretion episodes. Chandra observations were performed in early Aug.2014. Complementary JHK photometry and J-band spectra were taken at our Campo Imperatore facility to relate the X-ray and near-IR emission from the target. We observe a peak of X-ray emission at the nominal position of V1180 Cas. This signal corresponds to an X-ray luminosity L_X(0.5-7 kev) in the range 0.8-2.2e30 erg/s. Based on the relatively short duration of the dim states in the light curve and on stellar luminosity considerations, we explored the possibility that the brightness minima of V1180 Cas are driven by extinction variations. From the analysis of the spectral energy distribution of the high state we infer a stellar luminosity of 0.8-0.9 Lsun and find that the derived L_X is comparable to the average X-ray luminosities of T Tauri stars. Moreover, the X-ray luminosity is lower than the X-ray emission levels of 5e30 -1e31 erg/s detected at outbursts in similar low-mass objects. Our analysis suggests that at least part of the photometric fluctuations of V1180 Cas might be extinction effects rather than the result of accretion excess emission. However, as the source displays spectral features indicative of active accretion, we speculate that its photometric variations might be the result of a combination of accretion-induced and extinction-driven effects, as suggested for other young variables, such as V1184 Tau and V2492 Cyg.
△ Less
Submitted 25 September, 2015;
originally announced September 2015.
-
Starspot induced effects in microlensing events with rotating source star
Authors:
Mosè Giordano,
Achille A. Nucita,
Francesco De Paolis,
Gabriele Ingrosso
Abstract:
We consider the effects induced by the presence of hot and cold spots on the source star in the light curves of simulated microlensing events due to either single or binary lenses taking into account the rotation of the source star and the orbital motion of the lens system. Our goal is to study the anomalies induced by these effects on simulated microlensing light curves.
We consider the effects induced by the presence of hot and cold spots on the source star in the light curves of simulated microlensing events due to either single or binary lenses taking into account the rotation of the source star and the orbital motion of the lens system. Our goal is to study the anomalies induced by these effects on simulated microlensing light curves.
△ Less
Submitted 18 August, 2015;
originally announced August 2015.
-
On the relevance of chaos for halo stars in the Solar Neighbourhood
Authors:
Nicolás P. Maffione,
Facundo A. Gómez,
Pablo M. Cincotta,
Claudia M. Giordano,
Andrew P. Cooper,
Brian W. O'Shea
Abstract:
We show that diffusion due to chaotic mixing in the Neighbourhood of the Sun may not be as relevant as previously suggested in erasing phase space signatures of past Galactic accretion events. For this purpose, we analyse Solar Neighbourhood-like volumes extracted from cosmological simulations that naturally account for chaotic orbital behaviour induced by the strongly triaxial and cuspy shape of…
▽ More
We show that diffusion due to chaotic mixing in the Neighbourhood of the Sun may not be as relevant as previously suggested in erasing phase space signatures of past Galactic accretion events. For this purpose, we analyse Solar Neighbourhood-like volumes extracted from cosmological simulations that naturally account for chaotic orbital behaviour induced by the strongly triaxial and cuspy shape of the resulting dark matter haloes, among other factors. In the approximation of an analytical static triaxial model, our results show that a large fraction of stellar halo particles in such local volumes have chaos onset times (i.e., the timescale at which stars commonly associated with chaotic orbits will exhibit their chaotic behaviour) significantly larger than a Hubble time. Furthermore, particles that do present a chaotic behaviour within a Hubble time do not exhibit significant diffusion in phase space.
△ Less
Submitted 3 August, 2015;
originally announced August 2015.
-
Signatures of rotating binaries in micro-lensing experiments
Authors:
A. A. Nucita,
M. Giordano,
F. De Paolis,
G. Ingrosso
Abstract:
Gravitational microlensing offers a powerful method with which to probe a variety of binary-lens systems, as the binarity of the lens introduces deviations from the typical (single-lens) Paczyński behaviour in the event light curves. Generally, a static binary lens is considered to fit the observed light curve and, when the orbital motion is taken into account, an oversimplified model is usually e…
▽ More
Gravitational microlensing offers a powerful method with which to probe a variety of binary-lens systems, as the binarity of the lens introduces deviations from the typical (single-lens) Paczyński behaviour in the event light curves. Generally, a static binary lens is considered to fit the observed light curve and, when the orbital motion is taken into account, an oversimplified model is usually employed. In this paper, we treat the binary-lens motion in a realistic way and focus on simulated events that are fitted well by a Paczyński curve. We show that an accurate timing analysis of the residuals (calculated with respect to the best-fitting Paczyński model) is usually sufficient to infer the orbital period of the binary lens. It goes without saying that the independently estimated period may be used to further constrain the orbital parameters obtained by the best-fitting procedure, which often gives degenerate solutions. We also present a preliminary analysis of the event OGLE-2011-BLG-1127 / MOA-2011-BLG-322, which has been recognized to be the result of a binary lens. The period analysis results in a periodicity of \simeq 12 days, which confirms the oscillation of the observed data around the best-fitting model. The estimated periodicity is probably associated with an intrinsic variability of the source star, and therefore there is an opportunity to use this technique to investigate either the intrinsic variability of the source or the effects induced by the binary-lens orbital motion.
△ Less
Submitted 24 January, 2014;
originally announced January 2014.
-
Comparative study of variational chaos indicators and ODEs' numerical integrators
Authors:
Luciano A. Darriba,
Nicolás P. Maffione,
Pablo M. Cincotta,
Claudia M. Giordano
Abstract:
The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficie…
▽ More
The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed.
This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D Hénon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the
Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate
Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-Stöer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.
△ Less
Submitted 8 May, 2012; v1 submitted 4 May, 2012;
originally announced May 2012.
-
A comparison of different indicators of chaos based on the deviation vectors. Application to symplectic mappings
Authors:
N. P. Maffione,
L. A. Darriba,
P. M. Cincotta,
C. M. Giordano
Abstract:
The aim of this research work is to compare the reliability of several variational indicators of chaos on mappings. The Lyapunov Indicator (LI); the Mean Exponential Growth factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI); the Fast Lyapunov Indicator (FLI); the Dynamical Spectra of stretching numbers (SSN) and the corresponding Spectral Distance (D); and the Relative Lyapunov In…
▽ More
The aim of this research work is to compare the reliability of several variational indicators of chaos on mappings. The Lyapunov Indicator (LI); the Mean Exponential Growth factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI); the Fast Lyapunov Indicator (FLI); the Dynamical Spectra of stretching numbers (SSN) and the corresponding Spectral Distance (D); and the Relative Lyapunov Indicator (RLI), which is based on the evolution of the difference between two close orbits, have been included.
The experiments presented herein allow us to reliably suggest a group of chaos indicators to analyze a general mapping. We show that a package composed of the FLI and the RLI (to analyze the phase portrait globally) and the MEGNO and the SALI (to analyze orbits individually) is good enough to make a description of the systems' dynamics.
△ Less
Submitted 10 August, 2011;
originally announced August 2011.
-
Global Dynamics in Galactic Triaxial Systems I
Authors:
Pablo M. Cincotta,
Claudia M. Giordano,
Josefa Perez,
.
Abstract:
In this paper we present a theoretical analysis of the global dynamics in a triaxial galactic system using a 3D integrable Hamiltonian as a simple representation. We include a thorough discussion on the effect of adding a generic non--integrable perturbation to the global dynamics of the system. We adopt the triaxial Stackel Hamiltonian as the integrable model and compute its resonance structure…
▽ More
In this paper we present a theoretical analysis of the global dynamics in a triaxial galactic system using a 3D integrable Hamiltonian as a simple representation. We include a thorough discussion on the effect of adding a generic non--integrable perturbation to the global dynamics of the system. We adopt the triaxial Stackel Hamiltonian as the integrable model and compute its resonance structure in order to understand its global dynamics when a perturbation is introduced. Also do we take profit of this example in order to provide a theoretical discussion about diffussive processes taking place in phase space.
△ Less
Submitted 21 April, 2006;
originally announced April 2006.