-
TOI-2458 b: A mini-Neptune consistent with in situ hot Jupiter formation
Authors:
Ján Šubjak,
Davide Gandolfi,
Elisa Goffo,
David Rapetti,
Grzegorz Nowak,
Toshiyuki Mizuki,
Fei Dai,
Luisa M. Serrano,
Thomas G. Wilson,
Dawid Jankowski,
Krzysztof Goździewski,
Jon M. Jenkins,
Joseph D. Twicken,
Joshua N. Winn,
Allyson Bieryla,
William D. Cochran,
Karen A. Collins,
Hans J. Deeg,
Rafael A. García,
Eike W. Guenther,
Artie P. Hatzes,
Petr Kabáth,
Judith Korth,
David W. Latham,
John H. Livingston
, et al. (9 additional authors not shown)
Abstract:
We report on the discovery and spectroscopic confirmation of TOI-2458 b, a transiting mini-Neptune around an F-type star leaving the main-sequence with a mass of $M_\star=1.05 \pm 0.03$ M$_{\odot}$, a radius of $R_\star=1.31 \pm 0.03$ R$_{\odot}$, an effective temperature of $T_{\rm eff}=6005\pm50$ K, and a metallicity of $-0.10\pm0.05$ dex. By combining TESS photometry with high-resolution spectr…
▽ More
We report on the discovery and spectroscopic confirmation of TOI-2458 b, a transiting mini-Neptune around an F-type star leaving the main-sequence with a mass of $M_\star=1.05 \pm 0.03$ M$_{\odot}$, a radius of $R_\star=1.31 \pm 0.03$ R$_{\odot}$, an effective temperature of $T_{\rm eff}=6005\pm50$ K, and a metallicity of $-0.10\pm0.05$ dex. By combining TESS photometry with high-resolution spectra acquired with the HARPS spectrograph, we found that the transiting planet has an orbital period of $\sim$3.74 days, a mass of $M_p=13.31\pm0.99$ M$_{\oplus}$ and a radius of $R_p=2.83\pm0.20$ R$_{\oplus}$. The host star TOI-2458 shows a short activity cycle of $\sim$54 days revealed in the HARPS S-index time series. We took the opportunity to investigate other F stars showing activity cycle periods comparable to that of TOI-2458 and found that they have shorter rotation periods than would be expected based on the gyrochronology predictions. In addition, we determined TOI-2458's stellar inclination angle to be $i_\star\,=\,10.6_{-10.6}^{+13.3}$ degrees. We discuss that both phenomena (fast stellar rotation and planet orbit inclination) could be explained by in situ formation of a hot Jupiter interior to TOI-2458 b. It is plausible that this hot Jupiter was recently engulfed by the star. Analysis of HARPS spectra has identified the presence of another planet with a period of $P\,=\,16.55\pm0.06$ days and a minimum mass of $M_p \sin i=10.22\pm1.90$ M$_{\oplus}$.
△ Less
Submitted 26 September, 2024;
originally announced September 2024.
-
Tracking Advanced Planetary Systems (TAPAS) with HARPS-N VIII. A wide-orbit planetary companion in the hot-Jupiter system HD 118203
Authors:
G. Maciejewski,
A. Niedzielski,
K. Gozdziewski,
A. Wolszczan,
E. Villaver,
M. Fernandez,
M. Adamow,
J. Sierzputowska
Abstract:
The star HD 118203, classified as a K0 subgiant, was known to harbour a transiting hot Jupiter planet on a 6.1-day eccentric orbit. Previous studies also revealed a linear trend in the radial velocity (RV) domain, indicative of a companion on a wide orbit. Such a hierarchical orbital architecture could be helpful in studies of the origins of hot Jupiters. We acquired precise RV measurements over 1…
▽ More
The star HD 118203, classified as a K0 subgiant, was known to harbour a transiting hot Jupiter planet on a 6.1-day eccentric orbit. Previous studies also revealed a linear trend in the radial velocity (RV) domain, indicative of a companion on a wide orbit. Such a hierarchical orbital architecture could be helpful in studies of the origins of hot Jupiters. We acquired precise RV measurements over 17 years using the 9.2 m Hobby-Eberly Telescope and the 3.6 m Telescopio Nazionale Galileo. Combining these observations with space-born photometric time series from the Transiting Exoplanet Survey Satellite, we constructed a two-planetary model for the system. Astrometric observations from Hipparcos and Gaia were used to constrain the orbital inclination of the wide-orbit companion and its mass. Numerical simulations were used to investigate the dynamics of the system. The photometric data were searched for additional transit-like flux drops. We found that the additional companion is an 11-Jupiter mass planet orbiting HD 118203 on a 14-year moderately eccentric orbit, constituting a hierarchical planetary system with the hot Jupiter. Both planets were found to be dynamically decoupled mainly due to the general relativistic apsidal precession of the inner planet, marginalising secular interactions. The orbits of both planets might have a relatively low mutual inclination unless the longitudes of the ascending node differ substantially. This configuration favours the coplanar high-eccentricity migration as a path to the present-day orbital configuration. No other transiting planets with radii down to 2 Earth radii and orbital periods less than 100 days were found in the system.
△ Less
Submitted 16 July, 2024;
originally announced July 2024.
-
The orbital architecture and stability of the $μ$ Arae planetary system
Authors:
Krzysztof Goździewski
Abstract:
We re-analyze the global orbital architecture and dynamical stability of the $μ$ Arae planetary system. We have updated the best-fit elements and minimal masses of the planets based on literature radial velocity (RV) measurements, now spanning 15 years. This is twice the RVs interval used for the first characterization of the system in 2006. It consists of a Saturn- and two Jupiter-mass planets in…
▽ More
We re-analyze the global orbital architecture and dynamical stability of the $μ$ Arae planetary system. We have updated the best-fit elements and minimal masses of the planets based on literature radial velocity (RV) measurements, now spanning 15 years. This is twice the RVs interval used for the first characterization of the system in 2006. It consists of a Saturn- and two Jupiter-mass planets in low-eccentric orbits resembling the Earth-Mars-Jupiter configuration in the Solar system, as well as the close-in warm Neptune with a mass of ~14 Earth masses. Here, we constrain this early solution with the outermost period to be accurate to one month. The best-fit Newtonian model is characterized by moderate eccentricities of the most massive planets below 0.1 with small uncertainties ~0.02. It is close but meaningfully separated from the 2e:1b mean motion resonance of the Saturn-Jupiter-like pair, but may be close to weak three-body MMRs. The system appears rigorously stable over a wide region of parameter space covering uncertainties of several $σ$. The system stability is robust to a five-fold increase in the minimal masses, consistent with a wide range of inclinations, from 20 to 90 deg. This means that all planetary masses are safely below the brown dwarf mass limit. We found a weak statistical indication of the likely system inclination I~20-30 deg. Given the well constrained orbital solution, we also investigate the structure of hypothetical debris disks, which are analogs of the Main Belt and Kuiper Belt, and may naturally occur in this system.
△ Less
Submitted 22 September, 2022; v1 submitted 9 September, 2022;
originally announced September 2022.
-
Orbital and dynamical analysis of the system around HR 8799. New astrometric epochs from VLT/SPHERE and LBT/LUCI
Authors:
A. Zurlo,
K. Gozdziewski,
C. Lazzoni D. Mesa,
P. Nogueira,
S. Desidera,
R. Gratton,
F. Marzari,
E. Pinna,
G. Chauvin,
P. Delorme,
J. H. Girard,
J. Hagelberg,
Th. Henning,
M. Janson,
E. Rickman,
P. Kervella,
H. Avenhaus,
T. Bhowmik,
B. Biller,
A. Boccaletti,
M. Bonaglia,
M. Bonavita,
M. Bonnefoy,
F. Cantalloube,
A. Cheetham
, et al. (22 additional authors not shown)
Abstract:
HR\,8799 is a young planetary system composed of 4 planets and a double debris belt. Being the first multi-planetary system discovered with the direct imaging technique, it has been observed extensively since 1998. This wide baseline of astrometric measurements, counting over 50 observations in 20 years, permits a detailed orbital and dynamical analysis of the system. To explore the orbital parame…
▽ More
HR\,8799 is a young planetary system composed of 4 planets and a double debris belt. Being the first multi-planetary system discovered with the direct imaging technique, it has been observed extensively since 1998. This wide baseline of astrometric measurements, counting over 50 observations in 20 years, permits a detailed orbital and dynamical analysis of the system. To explore the orbital parameters of the planets, their dynamical history, and the planet-to-disk interaction, we made follow-up observations of the system during the VLT/SPHERE GTO program. We obtained 21 observations, most of them in favorable conditions. In addition, we observed HR\,8799 with the instrument LBT/LUCI. All the observations were reduced with state-of-the-art algorithms implemented to apply the spectral and angular differential imaging method. We re-reduced the SPHERE data obtained during the commissioning of the instrument and in 3 open-time programs to have homogeneous astrometry. The precise position of the 4 planets with respect to the host star was calculated by exploiting the fake negative companions method. To improve the orbital fitting, we also took into account all of the astrometric data available in the literature. From the photometric measurements obtained in different wavelengths, we estimated the planets' masses following the evolutionary models. We obtained updated parameters for the orbits with the assumption of coplanarity, relatively small eccentricities, and periods very close to the 2:1 resonance. We also refined the dynamical mass of each planet and the parallax of the system (24.49 $\pm$ 0.07 mas). We also conducted detailed $N$-body simulations indicating possible positions of a~putative fifth innermost planet with a mass below the present detection limits of $\simeq 3$~\MJup.
△ Less
Submitted 21 July, 2022;
originally announced July 2022.
-
Planet-planet scattering in presence of a companion star
Authors:
Francesco Marzari,
Makiko Nagasawa,
Krzyszof Goździewski
Abstract:
Planet Planet scattering is a leading dynamical mechanism invoked to explain the present orbital distribution of exoplanets. Many stars belong to binary systems, therefore it is important to understand how this mechanism works in presence of a companion star. We focus on systems of three planets orbiting the primary star and estimate the timescale for instability finding that it scales with the ke…
▽ More
Planet Planet scattering is a leading dynamical mechanism invoked to explain the present orbital distribution of exoplanets. Many stars belong to binary systems, therefore it is important to understand how this mechanism works in presence of a companion star. We focus on systems of three planets orbiting the primary star and estimate the timescale for instability finding that it scales with the keplerian period for systems that have the same ratio between inner planet and binary semimajor axes. An empirical formula is also derived from simulations to estimate how the the binary eccentricity affects the extent of the stability region. The presence of the secondary star affects the Planet Planet scattering outcomes causing a broadening of the final distribution in semimajor axis of the inner planet as some of the orbital energy of the planets is absorbed by the companion star. Repeated approaches to the secondary star causes also a significant reduction in the frequency of surviving two planet systems in particular for larger values of the inner planet semimajor axis. The formation of Kozai states with the companion star increases the number of planets which may be tidally circularized. To predict the possible final distribution of planets in binaries we have performed a large number of simulations where the initial semimajor axis of the inner planets is chosen randomly. For small values of the binary semimajor axis, the higher frequency of collision alter the final planet orbital distributions which, however, beyond 50 au appear to be scalable to wider binary separations.
△ Less
Submitted 9 December, 2021;
originally announced December 2021.
-
An exact, generalised Laplace resonance in the HR 8799 planetary system
Authors:
Krzysztof Gozdziewski,
Cezary Migaszewski
Abstract:
A system of four super-Jupiter planets around HR 8799 is the first multi-planet configuration discovered via the direct imaging technique. Despite over decade of research, the system's architecture remains not fully resolved. The main difficulty comes from still narrow observing window of ~20 years that covers small arcs of orbits with periods from roughly 50 to 500 years. Soon after the discovery…
▽ More
A system of four super-Jupiter planets around HR 8799 is the first multi-planet configuration discovered via the direct imaging technique. Despite over decade of research, the system's architecture remains not fully resolved. The main difficulty comes from still narrow observing window of ~20 years that covers small arcs of orbits with periods from roughly 50 to 500 years. Soon after the discovery it became clear that unconstrained best-fitting astrometric configurations self-disrupt rapidly, due to strong mutual gravitational interactions between ~10-Jupiter-mass companions. Recently, we showed that the HR 8799 system may be long term stable when locked in a generalized Laplace 8:4:2:1 mean motion resonance (MMR) chain, and we constrained its orbits through the planetary migration. Here we qualitatively improve this approach by considering the MMR in terms of an exactly periodic configuration. This assumption enables us to construct for the first time the self-consistent $N$-body model of the long-term stable orbital architecture, using only available astrometric positions of the planets relative to the star. We independently determine planetary masses, which are consistent with thermodynamic evolution, and the parallax overlapping to $1σ$ with the most recent GAIA DR2 value. We also determine the global structure of the inner and outer debris discs in the [8, 600] au range, consistent with the updated orbital solution.
△ Less
Submitted 15 September, 2020;
originally announced September 2020.
-
Two Super-Earths in the 3:2 MMR around KOI-1599
Authors:
F. Panichi,
C. Migaszewski,
K. Goździewski
Abstract:
We validate the planetary origin of the KOI-1599 transit time variations (TTVs) with statistical and dynamical tests. We re-analysed KEPLER Q1-Q17 light-curves of the star, and we independently derived the TTVs. They appear as strongly anti-correlated, suggestive of two mutually interacting planets. We found similar radii of the candidates, $1.9 \pm 0.2\, \mbox{R}_{\oplus}$ for the inner KOI-1599.…
▽ More
We validate the planetary origin of the KOI-1599 transit time variations (TTVs) with statistical and dynamical tests. We re-analysed KEPLER Q1-Q17 light-curves of the star, and we independently derived the TTVs. They appear as strongly anti-correlated, suggestive of two mutually interacting planets. We found similar radii of the candidates, $1.9 \pm 0.2\, \mbox{R}_{\oplus}$ for the inner KOI-1599.02, and $1.9 \pm 0.3 \, \mbox{R}_{\oplus}$ for the outer KOI-1599.01. The standard MCMC TTV analysis constrains the planet masses safely below the dynamical instability limit of $\simeq 3\, \mbox{M}_{\mbox{Jup}} $. The best-fitting MCMC model yields $(9.0\pm 0.3)\, \mbox{R}_{\oplus}$, and $(4.6\pm0.3),\ \mbox{R}_{\oplus}$, for the inner and the outer planet, respectively. The planets are trapped in 3:2 mean motion resonance (MMR) with anti-aligned apsides ($Δ\varpi = 180^{\circ}$) at low-eccentric ($e\simeq 0.01)$ orbits. However, we found that the TTV mass determination depends on eccentricity priors with the dispersion in the (0.01,0.05) range. They permit a second family of TTV models with smaller masses of $\simeq 7\,\mbox{R}_{\oplus}$, and $\simeq 3.6\,\mbox{R}_{\oplus}$, respectively, exhibiting two modes of $Δ\varpi = 0^{\circ},180^{\circ}$ librations. The 3:2 MMR is dynamically robust and persists for both modes. In order to resolve the mass duality, we re-analysed the TTV data with a quasi-analytic model of resonant TTV signals. This model favours the smaller masses. We also reproduced this model with simulating the migration capture of the system into the 3:2 MMR.
△ Less
Submitted 24 March, 2019; v1 submitted 5 January, 2019;
originally announced January 2019.
-
The orbital architecture and debris disks of the HR 8799 planetary system
Authors:
Krzysztof Gozdziewski,
Cezary Migaszewski
Abstract:
The HR8799 planetary system with four ~10 mJup planets in wide orbits up to 70 au, and periods up to 500 yr has been detected with the direct imaging. Its intriguing orbital architecture is not fully resolved due to time-limited astrometry covering ~20 years. Earlier, we constructed a heuristic model of the system based on rapid, convergent migration of the planets. We developed a better structure…
▽ More
The HR8799 planetary system with four ~10 mJup planets in wide orbits up to 70 au, and periods up to 500 yr has been detected with the direct imaging. Its intriguing orbital architecture is not fully resolved due to time-limited astrometry covering ~20 years. Earlier, we constructed a heuristic model of the system based on rapid, convergent migration of the planets. We developed a better structured and CPU-efficient variant of this model. We re-analyzed the self-consistent, homogeneous astrometric dataset in Konopacky (2016). The best-fitting configuration agrees with our earlier findings. The planets are likely involved in dynamically robust Laplace resonance chain. Hypothetical planets with masses below the current detection limit of 0.1-3 Jupiter masses, within the observed inner, or beyond the outer orbit, respectively, do not influence the long term stability of the system. We predict positions of such non-detected objects. The long-term stable orbital model of the observed planets helps to simulate the dynamical structure of debris disks in the system. A CPU-efficient fast indicator technique makes it possible to reveal their complex, resonant shape in 10^6 particles scale. We examine the inner edge of the outer disk detected between 90-145 au. We also reconstruct the outer disk assuming that it has been influenced by convergent migration of the planets. A complex shape of the disk strongly depends on various dynamical factors, like orbits and masses of non-detected planets. It may be highly non-circular and its models are yet non-unique, regarding both observational constraints, as well as its origin.
△ Less
Submitted 16 July, 2018;
originally announced July 2018.
-
A periodic configuration of the Kepler-25 planetary system?
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski
Abstract:
We study proximity of the Kepler-25 planetary system to a periodic configuration, which is known to be the final state of a system that undergoes smooth migration resulting from the planet-disc interaction. We show that the system is close to the periodic configuration of 2:1 mean motion resonance (MMR) what indicates that its past migration was neither disturbed significantly by turbulence in the…
▽ More
We study proximity of the Kepler-25 planetary system to a periodic configuration, which is known to be the final state of a system that undergoes smooth migration resulting from the planet-disc interaction. We show that the system is close to the periodic configuration of 2:1 mean motion resonance (MMR) what indicates that its past migration was neither disturbed significantly by turbulence in the disc nor the orbits were perturbed by planetesimals that left after the disc dispersal. We show that, because of the TTV model degeneracy, a periodic configuration is difficult to be found when the standard modelling of the transit timing variations (TTVs) is used. The TTV signal of a periodic configuration (with anti-aligned apsidal lines) may be misinterpreted as an aligned non-resonant system. We demonstrate that the standard MCMC modelling of the Kepler-25 TTVs is very sensitive to an a~priori information on the eccentricities (prior probability distributions). Wide priors (of the order of the ones typically used in the literature) result in favouring the aligned non-resonant configurations with small planets' masses and moderate eccentricities, while for the narrower priors the most likely are the anti-aligned resonant systems with larger masses and very low eccentricities.
△ Less
Submitted 20 July, 2018; v1 submitted 27 March, 2018;
originally announced March 2018.
-
Another look at AM Herculis $-$ radio-astrometric campaign with the e-EVN at 6 cm
Authors:
M. P. Gawroński,
K. Goździewski,
K. Katarzyński,
G. Rycyk
Abstract:
We conducted radio-interferometric observations of the well known binary cataclysmic system AM Herculis. This particular system is formed by the magnetic white dwarf (primary) and the red dwarf (secondary), and is the prototype of so-called polars. Our observations were conducted with the European VLBI Network (EVN) in the e-EVN mode at 5 GHz. We obtained six astrometric measurements spanning one…
▽ More
We conducted radio-interferometric observations of the well known binary cataclysmic system AM Herculis. This particular system is formed by the magnetic white dwarf (primary) and the red dwarf (secondary), and is the prototype of so-called polars. Our observations were conducted with the European VLBI Network (EVN) in the e-EVN mode at 5 GHz. We obtained six astrometric measurements spanning one year, which make it possible to update the annual parallax for this system with the best precision to date ($π=11.29\pm0.08$ mas), equivalent to the distance of $88.6\pm0.6$ pc. The system was observed mostly in the quiescent phase (visual magnitude $m_v\sim15.3$), when the radio emission was at the level of about $300~μ$Jy. Our analysis suggests that the radio flux of AM Herculis is modulated with the orbital motion. Such specific properties of the radiation can be explained using the emission mechanism similar to the scenario proposed for V471 Tauri and, in general, for RS CVn type stars. In this scenario the radio emission arises near the surface of the red dwarf, where the global magnetic field strength may reach a few kG. We argue that the quiescent radio emission distinguishes AM Herculis together with AR Ursae Majoris (the second known persistent radio polar) from other polars, as the systems with a magnetized secondary star.
△ Less
Submitted 17 January, 2018;
originally announced January 2018.
-
The architecture and formation of the Kepler-30 planetary system
Authors:
Federico Panichi,
Krzysztof Goździewski,
Cezary Migaszewski,
Ewa Szuszkiewicz
Abstract:
We study the orbital architecture, physical characteristics of planets, formation and long-term evolution of the Kepler-30 planetary system, detected and announced in 2012 by the KEPLER team. We show that the Kepler-30 system belongs to a particular class of very compact and quasi-resonant, yet long-term stable planetary systems. We re-analyse the light curves of the host star spanning Q1-Q17 quar…
▽ More
We study the orbital architecture, physical characteristics of planets, formation and long-term evolution of the Kepler-30 planetary system, detected and announced in 2012 by the KEPLER team. We show that the Kepler-30 system belongs to a particular class of very compact and quasi-resonant, yet long-term stable planetary systems. We re-analyse the light curves of the host star spanning Q1-Q17 quarters of the KEPLER mission. A huge variability of the Transit Timing Variations (TTV) exceeding 2 days is induced by a massive Jovian planet located between two Neptune-like companions. The innermost pair is near to the 2:1 mean motion resonance (MMR), and the outermost pair is close to higher order MMRs, such as 17:7 and 7:3. Our re-analysis of photometric data allows us to constrain, better than before, the orbital elements, planets' radii and masses, which are $9.2 \pm 0.1$, $536 \pm 5$, and $23.7 \pm 1.3$ Earth masses for Kepler-30b, Kepler-30c and Kepler-30d, respectively. The masses of the inner planets are determined within $\sim 1\%$ uncertainty. We infer the internal structures of the Kepler-30 planets and their bulk densities in a wide range from $ (0.19 \pm 0.01)$ g$\cdot$cm$^{-3}$ for Kepler-30d, $(0.96\pm0.15)$ g$\cdot$cm$^{-3}$ for Kepler-30b, to $(1.71 \pm 0.13)$ g$\cdot$cm$^{-3}$ for the Jovian planet Kepler-30c. We attempt to explain the origin of this unique planetary system and a deviation of the orbits from exact MMRs through the planetary migration scenario. We anticipate that the Jupiter-like planet plays an important role in determining the present dynamical state of this system.
△ Less
Submitted 27 February, 2018; v1 submitted 16 July, 2017;
originally announced July 2017.
-
The Reversibility Error Method (REM): a new, dynamical fast indicator for planetary dynamics
Authors:
Federico Panichi,
Krzyszof Goździewski,
Giorgio Turchetti
Abstract:
We describe the Reversibility Error Method (REM) and its applications to planetary dynamics. REM is based on the time-reversibility analysis of the phase-space trajectories of conservative Hamiltonian systems. The round-off errors break the time reversibility and the displacement from the initial condition, occurring when we integrate it forward and backward for the same time interval, is related…
▽ More
We describe the Reversibility Error Method (REM) and its applications to planetary dynamics. REM is based on the time-reversibility analysis of the phase-space trajectories of conservative Hamiltonian systems. The round-off errors break the time reversibility and the displacement from the initial condition, occurring when we integrate it forward and backward for the same time interval, is related to the dynamical character of the trajectory. If the motion is chaotic, in the sense of non-zero maximal Characteristic Lyapunov Exponent (mLCE), then REM increases exponentially with time, as exp λt, while when the motion is regular (quasi-periodic) then REM increases as a power law in time, as t^α, where α and λ are real coefficients. We compare the REM with a variant of mLCE, the Mean Exponential Growth factor of Nearby Orbits (MEGNO). The test set includes the restricted three body problem and five resonant planetary systems: HD 37124, Kepler-60, Kepler-36, Kepler-29 and Kepler-26. We found a very good agreement between the outcomes of these algorithms. Moreover, the numerical implementation of REM is astonishing simple, and is based on solid theoretical background. The REM requires only a symplectic and time-reversible (symmetric) integrator of the equations of motion. This method is also CPU efficient. It may be particularly useful for the dynamical analysis of multiple planetary systems in the KEPLER sample, characterized by low-eccentricity orbits and relatively weak mutual interactions. As an interesting side-result, we found a possible stable chaos occurrence in the Kepler-29 planetary system.
△ Less
Submitted 30 March, 2017;
originally announced March 2017.
-
Is there a circumbinary planet around NSVS 14256825?
Authors:
Ilham Nasiroglu,
Krzysztof Gozdziewski,
Aga Slowikowska,
Krzysztof Krzeszowski,
Michal Zejmo,
Staszek Zola,
Huseyin Er,
Waldemar Ogloza,
Marek Drozdz,
Dorota Koziel-Wierzbowska,
Bartlomiej Debski,
Nazil Karaman
Abstract:
The cyclic behaviour of (O-C) residuals of eclipse timings in the sdB+M eclipsing binary NSVS 14256825 was previously attributed to one or two Jovian-type circumbinary planets. We report 83 new eclipse timings that not only fill in the gaps in those already published but also extend the time span of the (O-C) diagram by three years. Based on the archival and our new data spanning over more than 17…
▽ More
The cyclic behaviour of (O-C) residuals of eclipse timings in the sdB+M eclipsing binary NSVS 14256825 was previously attributed to one or two Jovian-type circumbinary planets. We report 83 new eclipse timings that not only fill in the gaps in those already published but also extend the time span of the (O-C) diagram by three years. Based on the archival and our new data spanning over more than 17 years we re-examined the up to date system (O-C). The data revealed systematic, quasi-sinusoidal variation deviating from an older linear ephemeris by about 100 s. It also exhibits a maximum in the (O-C) near JD 2,456,400 that was previously unknown. We consider two most credible explanations of the (O-C) variability: the light propagation time due to the presence of an invisible companion in a distant circumbinary orbit, and magnetic cycles reshaping one of the binary components, known as the Applegate or Lanza-Rodono effect. We found that the latter mechanism is unlikely due to the insufficient energy budget of the M-dwarf secondary. In the framework of the third-body hypothesis, we obtained meaningful constraints on the Keplerian parameters of a putative companion and its mass. Our best-fitting model indicates that the observed quasi-periodic (O-C) variability can be explained by the presence of a brown dwarf with the minimal mass of 15 Jupiter masses rather than a planet, orbiting the binary in a moderately elliptical orbit (~ 0.175) with the period of ~ 10 years. Our analysis rules out two planets model proposed earlier.
△ Less
Submitted 18 January, 2017;
originally announced January 2017.
-
Physical properties and astrometry of radio-emitting brown dwarf TVLM513-46546 revisited
Authors:
Marcin P. Gawronski,
Krzysztof Gozdziewski,
Krzysztof Katarzynski
Abstract:
We present multi-epoch astrometric observations of the M9 ultra-cool dwarf TVLM513-46546 that is placed at the brown dwarf boundary. The new observations have been performed with the European VLBI Network (EVN) at 6cm band. The target has been detected at 7~epochs spanning three years, with measured quiescent emission flux in the range 180-300 muJy. We identified four short-duration flaring events…
▽ More
We present multi-epoch astrometric observations of the M9 ultra-cool dwarf TVLM513-46546 that is placed at the brown dwarf boundary. The new observations have been performed with the European VLBI Network (EVN) at 6cm band. The target has been detected at 7~epochs spanning three years, with measured quiescent emission flux in the range 180-300 muJy. We identified four short-duration flaring events (0.5-2 mJy) with very high circular polarization (~75%-100%). Properties of the observed radio flares support the physical model of the source that is characterized by the electron cyclotron maser instability responsible for outbursts of radio emission. Combined with Very Long Baseline Array (VLBA) earlier data, our detections make it possible to refine the absolute parallax pi=93.27(+0.18,-0.17) mas. Our measurements rule out TVLM513-46546 companions more massive than Jupiter in orbits with periods longer than ~1 yr.
△ Less
Submitted 16 December, 2016;
originally announced December 2016.
-
The origin and 9:7 MMR dynamics of the Kepler-29 system
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski,
Federico Panichi
Abstract:
We analyse the Transit Timing Variation (TTV) measurements of a~system of two super-Earths detected as Kepler-29, in order to constrain the planets' masses and orbital parameters. A dynamical analysis of the best-fitting configurations constrains the masses to be $\sim 6$ and $\sim 5$ Earth masses for the inner and the outer planets, respectively. The analysis also reveals that the system is likel…
▽ More
We analyse the Transit Timing Variation (TTV) measurements of a~system of two super-Earths detected as Kepler-29, in order to constrain the planets' masses and orbital parameters. A dynamical analysis of the best-fitting configurations constrains the masses to be $\sim 6$ and $\sim 5$ Earth masses for the inner and the outer planets, respectively. The analysis also reveals that the system is likely locked in the 9:7~mean motion resonance. However, a variety of orbital architectures regarding eccentricities and the relative orientation of orbits is permitted by the observations as well as by stability constraints. We attempt to find configurations preferred by the planet formation scenarios as an additional, physical constraint. We show that configurations with low eccentricities and anti-aligned apsidal lines of the orbits are a natural and most likely outcome of the convergent migration. However, we show that librations of the critical angles are not necessary for the Kepler-29 system to be dynamically resonant, and such configurations may be formed on the way of migration as well. We argue, on the other hand, that aligned configurations with $e \gtrsim 0.03$ may be not consistent with the migration scenario.
△ Less
Submitted 3 November, 2016; v1 submitted 28 September, 2016;
originally announced September 2016.
-
Search for exoplanets and brown dwarfs with VLBI
Authors:
K. Katarzynski,
M. Gawronski,
K. Gozdziewski
Abstract:
The main aim of this work is to estimate possible radio GHz emission of extrasolar planets and brown dwarfs and to check if such radiation can be detected by Very Large Baseline Interferometers (VLBI). In the estimation we assume that the emission may originate in processes similar to those observed in the Jupiter system. The frequency of the radio emission that is produced in this system depends…
▽ More
The main aim of this work is to estimate possible radio GHz emission of extrasolar planets and brown dwarfs and to check if such radiation can be detected by Very Large Baseline Interferometers (VLBI). In the estimation we assume that the emission may originate in processes similar to those observed in the Jupiter system. The frequency of the radio emission that is produced in this system depends mostly on the magnetic field strength. Jupiter's magnetic field ($\sim 9$ G on average) allows for radiation from kHz frequencies up to 40 MHz. This is well below the frequency range of VLBI. However, it was demonstrated that the magnetic field strength in massive and young object may be up to two orders of magnitude higher than for Jupiter, which is especially relevant for planets around short-lived A type stars. This should extend the range of the emission up to GHz frequencies. We calculated expected flux densities of radio emission for a variety of hypothetical young planetary systems. We analysed two different emission scenarios, and found that the radiation induced by moons (process similar to Jupiter-Io interactions) appears to be less efficient than the emission generated by a stellar wind on a planetary magnetosphere. We also estimated hypothetical emission of planets and brown dwarfs located around relatively young and massive main sequence A-type stars. Our results show that the emission produced by stellar winds could be detected by currently operating VLBI networks.
△ Less
Submitted 24 August, 2016;
originally announced August 2016.
-
The Laplace resonance in the Kepler-60 system
Authors:
Krzysztof Gozdziewski,
Cezary Migaszewski,
Federico Panichi,
Ewa Szuszkiewicz
Abstract:
We investigate the dynamical stability of the Kepler-60 planetary system with three super-Earths. We first determine their orbital elements and masses by Transit Timing Variation (TTV) data spanning quarters Q1-Q16 of the KEPLER mission. The system is dynamically active but the TTV data constrain masses to ~4 Earth masses and orbits in safely wide stable zones. The observations prefer two types of…
▽ More
We investigate the dynamical stability of the Kepler-60 planetary system with three super-Earths. We first determine their orbital elements and masses by Transit Timing Variation (TTV) data spanning quarters Q1-Q16 of the KEPLER mission. The system is dynamically active but the TTV data constrain masses to ~4 Earth masses and orbits in safely wide stable zones. The observations prefer two types of solutions. The true three-body Laplace MMR exhibits the critical angle librating around 45 degrees and aligned apsides of the inner and outer pair of planets. In the Laplace MMR formed through a chain of two-planet 5:4 and 4:3 MMRs, all critical angles librate with small amplitudes of ~30 degrees and apsidal lines in planet's pairs are anti-aligned. The system is simultaneously locked in a three-body MMR with librations amplitude of ~10 degrees. The true Laplace MMR can evolve towards a chain of two-body MMRs in the presence of planetary migration. Therefore the three-body MMR formed in this way seems to be more likely state of the system. However, the true three-body MMR cannot be disregarded a priori and it remains a puzzling configuration that may challenge the planet formation theory.
△ Less
Submitted 9 October, 2015;
originally announced October 2015.
-
The HU Aqr planetary system hypothesis revisited
Authors:
K. Gozdziewski,
A. Slowikowska,
D. Dimitrov,
K. Krzeszowski,
M. Zejmo,
G. Kanbach,
V. Burwitz,
A. Rau,
P. Irawati,
A. Richichi,
M. Gawronski,
G. Nowak,
I. Nasiroglu,
D. Kubicki
Abstract:
We study the mid-egress eclipse timing data gathered for the cataclysmic binary HU Aquarii during the years 1993-2014. The (O-C) residuals were previously attributed to a single ~7 Jupiter mass companion in ~5 au orbit or to a stable 2-planet system with an unconstrained outermost orbit. We present 22 new observations gathered between June, 2011 and July, 2014 with four instruments around the worl…
▽ More
We study the mid-egress eclipse timing data gathered for the cataclysmic binary HU Aquarii during the years 1993-2014. The (O-C) residuals were previously attributed to a single ~7 Jupiter mass companion in ~5 au orbit or to a stable 2-planet system with an unconstrained outermost orbit. We present 22 new observations gathered between June, 2011 and July, 2014 with four instruments around the world. They reveal a systematic deviation of ~60 - 120 seconds from the older ephemeris. We re-analyse the whole set of the timing data available. Our results provide an erratum to the previous HU Aqr planetary models, indicating that the hypothesis for a third and fourth body in this system is uncertain. The dynamical stability criterion and a particular geometry of orbits rule out coplanar 2-planet configurations. A putative HU Aqr planetary system may be more complex, e.g., highly non-coplanar. Indeed, we found examples of 3-planet configurations with the middle planet in a retrograde orbit, which are stable for at least 1Gyr, and consistent with the observations. The (O-C) may be also driven by oscillations of the gravitational quadrupole moment of the secondary, as predicted by the Lanza et al. modification of the Applegate mechanism. Further systematic, long-term monitoring of HU Aqr is required to interpret the (O-C) residuals.
△ Less
Submitted 18 December, 2014;
originally announced December 2014.
-
Predicting a third planet in the Kepler-47 circumbinary system
Authors:
Tobias C. Hinse,
Nader Haghighipour,
Veselin B. Kostov,
Krzysztof Goździewski
Abstract:
We have studied the possibility that a third circumbinary planet in the Kepler-47 planetary system be the source of the single unexplained transiting event reported during the discovery of these planets. We applied the MEGNO technique to identify regions in the phase space where a third planet can maintain quasi-periodic orbits, and assessed the long-term stability of the three-planet system by in…
▽ More
We have studied the possibility that a third circumbinary planet in the Kepler-47 planetary system be the source of the single unexplained transiting event reported during the discovery of these planets. We applied the MEGNO technique to identify regions in the phase space where a third planet can maintain quasi-periodic orbits, and assessed the long-term stability of the three-planet system by integrating the entire 5 bodies (binary + planets) for 10 Myr. We identified several stable regions between the two known planets as well as a region beyond the orbit of Kepler-47c where the orbit of the third planet could be stable. To constrain the orbit of this planet, we used the measured duration of the unexplained transit event ($\sim 4.15$ hours) and compared that with the transit duration of the third planet in an ensemble of stable orbits. To remove the degeneracy among the orbits with similar transit durations, we considered the planet to be in a circular orbit and calculated its period analytically. The latter places an upper limit of 424 days on the orbital period of the third planet. Our analysis suggests that if the unexplained transit event detected during the discovery of the Kepler-47 circumbinary system is due to a planetary object, this planet will be in a low eccentricity orbit with a semimajor axis smaller than 1.24 AU. Further constraining of the mass and orbital elements of this planet requires a re-analysis of the entire currently available data including those obtained post-announcement of the discovery of this system. We present details of our methodology and discuss the implication of the results.
△ Less
Submitted 17 November, 2014; v1 submitted 4 September, 2014;
originally announced September 2014.
-
Mechanic: the MPI/HDF code framework for dynamical astronomy
Authors:
Mariusz Słonina,
Krzysztof Goździewski,
Cezary Migaszewski
Abstract:
We introduce the Mechanic, a new open-source code framework. It is designed to reduce the development effort of scientific applications by providing unified API (Application Programming Interface) for configuration, data storage and task management. The communication layer is based on the well-established Message Passing Interface (MPI) standard, which is widely used on variety of parallel compute…
▽ More
We introduce the Mechanic, a new open-source code framework. It is designed to reduce the development effort of scientific applications by providing unified API (Application Programming Interface) for configuration, data storage and task management. The communication layer is based on the well-established Message Passing Interface (MPI) standard, which is widely used on variety of parallel computers and CPU-clusters. The data storage is performed within the Hierarchical Data Format (HDF5). The design of the code follows em core-module approach which allows to reduce the user's codebase and makes it portable for single- and multi-CPU environments. The framework may be used in a local user's environment, without administrative access to the cluster, under the PBS or Slurm job schedulers. It may become a helper tool for a wide range of astronomical applications, particularly focused on processing large data sets, such as dynamical studies of long-term orbital evolution of planetary systems with Monte Carlo methods, dynamical maps or evolutionary algorithms. It has been already applied in numerical experiments conducted for Kepler-11 (Migaszewski et al., 2012), and nuOctantis planetary systems (Goździewski et al., 2013). In this paper we describe the basics of the framework, including code listings for the implementation of a sample user's module. The code is illustrated on a model Hamiltonian introduced by (Froeschle et al., 2000) presenting the Arnold diffusion. The Arnold Web is shown with the help of the MEGNO (Mean Exponential Growth of Nearby Orbits) fast indicator (Goździewski et al., 2008a) applied onto symplectic SABA integrators family (Laskar and Robutel, 2001).
△ Less
Submitted 24 January, 2014;
originally announced January 2014.
-
Revisiting the proposed circumbinary multi-planet system NSVS14256825
Authors:
Tobias C. Hinse,
Jae Woo Lee,
Krzysztof Gozdziewski,
Jonathan Horner,
Robert A. Wittenmyer
Abstract:
In this work we carry out an analysis of the observed times of primary and secondary eclipses of the post-common envelope binary NSVS14256825. Recently, \cite{Almeida2013} proposed that two circumbinary companions orbit this short-period eclipsing binary, in order to explain observed variations in the timing of mutual eclipses between the two binary components. Using a standard weighted least-squa…
▽ More
In this work we carry out an analysis of the observed times of primary and secondary eclipses of the post-common envelope binary NSVS14256825. Recently, \cite{Almeida2013} proposed that two circumbinary companions orbit this short-period eclipsing binary, in order to explain observed variations in the timing of mutual eclipses between the two binary components. Using a standard weighted least-squares minimisation technique, we have extensively explored the topology of $χ^2$ parameter space of a single planet model. We find the data set to be insufficient to reliably constrain a one-companion model. Various models, each with similar statistical significance, result in substantially different orbital architectures for the additional companion. No evidence is seen for a second companion of planetary nature. We suspect insufficient coverage (baseline) of timing data causing the best-fit parameters to be unconstrained.
△ Less
Submitted 13 November, 2013;
originally announced November 2013.
-
Project RISARD
Authors:
Marcin P. Gawronski,
Krzysztof Gozdziewski,
Krzysztof Katarzynski
Abstract:
Red Dwarfs (RDs) are the most common, low-mass stars (<=0.5 M_sun) in the Solar neighbourhood, and probably in the Universe as well. Most (likely all) young RDs are magnetically active, and therefore it is impossible to measure their radial velocities (RVs) with the accuracy required by contemporary planetary searches (5-50 m/s). Some of RDs are known as a source of variable radio emission at cent…
▽ More
Red Dwarfs (RDs) are the most common, low-mass stars (<=0.5 M_sun) in the Solar neighbourhood, and probably in the Universe as well. Most (likely all) young RDs are magnetically active, and therefore it is impossible to measure their radial velocities (RVs) with the accuracy required by contemporary planetary searches (5-50 m/s). Some of RDs are known as a source of variable radio emission at centimetre wavelengths. This radiation is relatively weak (~0.2 to 1 mJy). However, recent development of the EVN and e-VLBI systems in terms of sensitivity gives a new possibility to investigate such stars. Here, we present first results from our ongoing RISARD project (Radio Interferometric Survey of Active Red Dwarfs). The main goal of this project is to detect an exoplanet by direct, precise measurements of a RD position and possible changes to this position caused by the planet.
△ Less
Submitted 18 September, 2013;
originally announced September 2013.
-
Multiple mean motion resonances in the HR 8799 planetary system
Authors:
Krzysztof Gozdziewski,
Cezary Migaszewski
Abstract:
HR 8799 is a nearby star hosting at least four ~10 Jovian mass planets in wide orbits up to ~70au, detected through the direct, high-contrast infrared imaging. Large companions and debris disks reported interior to ~10au, and exterior to ~100au indicate massive protoplanetary disc in the past. The dynamical state of the HR 8799 system is not yet fully resolved, due to limited astrometric data cove…
▽ More
HR 8799 is a nearby star hosting at least four ~10 Jovian mass planets in wide orbits up to ~70au, detected through the direct, high-contrast infrared imaging. Large companions and debris disks reported interior to ~10au, and exterior to ~100au indicate massive protoplanetary disc in the past. The dynamical state of the HR 8799 system is not yet fully resolved, due to limited astrometric data covering tiny orbital arcs. We construct a new, orbital model of the HR 8799 system, assuming rapid migration of the planets after their formation in wider orbits. We found that the HR 8799 planets are likely involved in double Laplace resonance, 1e:2d:4c:8b MMR. Quasi-circular planetary orbits are coplanar with the stellar equator and inclined by ~25 degrees to the sky plane. This best-fit orbital configuration matches astrometry, debris disk models, and mass estimates from cooling models. The multiple MMR is stable for the age of the star ~160Myr, for at least 1 Gyr unless significant perturbations to the N-body dynamics are present. We predict four configurations with the fifth hypothetical innermost planet HR 8799f in ~9.7au, or ~7.5au orbit, extending the MMR chain to triple Laplace resonance 1f:2e:4d:8c:16b MMR or to the 1f:3e:6d:12c:24b MMR, respectively. Our findings may establish strong boundary conditions for the system formation and its early history.
△ Less
Submitted 6 March, 2014; v1 submitted 29 August, 2013;
originally announced August 2013.
-
A linear distribution of orbits in compact planetary systems?
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski,
Mariusz Slonina
Abstract:
We report a linear ordering of orbits in a sample of multiple extrasolar planetary systems with super-Earth planets. We selected 20 cases, mostly discovered by the Kepler mission, hosting at least four planets within \sim 0.5 au. The semi-major axis a_n of an n-th planet in each system of this sample obeys a(n) = a_1 + (n-1) Δa, where a_1 is the semi-major axis of the innermost orbit and Δa is a s…
▽ More
We report a linear ordering of orbits in a sample of multiple extrasolar planetary systems with super-Earth planets. We selected 20 cases, mostly discovered by the Kepler mission, hosting at least four planets within \sim 0.5 au. The semi-major axis a_n of an n-th planet in each system of this sample obeys a(n) = a_1 + (n-1) Δa, where a_1 is the semi-major axis of the innermost orbit and Δa is a spacing between subsequent planets, which are specific for a particular system. For instance, the Kepler-33 system hosting five super-Earth planets exhibits the relative deviations between the observed and linearly predicted semi-major axes of only a few percent. At least half of systems in the sample fulfill the linear law with a similar accuracy. We explain the linear distribution of semi-major axes as a natural implication of multiple chains of mean motion resonances between subsequent planets, which emerge due to planet--disk interactions and convergent migration at early stages of their evolution.
△ Less
Submitted 6 August, 2013; v1 submitted 14 June, 2013;
originally announced June 2013.
-
Multi-site campaign for transit timing variations of WASP-12 b: possible detection of a long-period signal of planetary origin
Authors:
G. Maciejewski,
D. Dimitrov,
M. Seeliger,
St. Raetz,
L. Bukowiecki,
M. Kitze,
R. Errmann,
G. Nowak,
A. Niedzielski,
V. Popov,
C. Marka,
K. Gozdziewski,
R. Neuhaeuser,
J. Ohlert,
T. C. Hinse,
J. W. Lee,
C. -U. Lee,
J. -N. Yoon,
A. Berndt,
H. Gilbert,
Ch. Ginski,
M. M. Hohle,
M. Mugrauer,
T. Roell,
T. O. B. Schmidt
, et al. (16 additional authors not shown)
Abstract:
The transiting planet WASP-12 b was identified as a potential target for transit timing studies because a departure from a linear ephemeris was reported in the literature. Such deviations could be caused by an additional planet in the system. We attempt to confirm the existence of claimed variations in transit timing and interpret its origin. We organised a multi-site campaign to observe transits…
▽ More
The transiting planet WASP-12 b was identified as a potential target for transit timing studies because a departure from a linear ephemeris was reported in the literature. Such deviations could be caused by an additional planet in the system. We attempt to confirm the existence of claimed variations in transit timing and interpret its origin. We organised a multi-site campaign to observe transits by WASP-12 b in three observing seasons, using 0.5-2.6-metre telescopes. We obtained 61 transit light curves, many of them with sub-millimagnitude precision. The simultaneous analysis of the best-quality datasets allowed us to obtain refined system parameters, which agree with values reported in previous studies. The residuals versus a linear ephemeris reveal a possible periodic signal that may be approximated by a sinusoid with an amplitude of 0.00068+/-0.00013 d and period of 500+/-20 orbital periods of WASP-12 b. The joint analysis of timing data and published radial velocity measurements results in a two-planet model which better explains observations than single-planet scenarios. We hypothesize that WASP-12 b might be not the only planet in the system and there might be the additional 0.1 M_Jup body on a 3.6-d eccentric orbit. A dynamical analysis indicates that the proposed two-planet system is stable over long timescales.
△ Less
Submitted 25 January, 2013;
originally announced January 2013.
-
Investigating AM Her Cataclysmic Variables with OPTIMA
Authors:
Aga Słowikowska,
Krzysztof Goździewski,
Ilham Nasiroglu,
Gottfried Kanbach,
Arne Rau,
Krzysztof Krzeszowski
Abstract:
We focus on short--period eclipsing binaries that belong to a class of Cataclysmic Variables (CVs). They are known as polars and intermediate polars, closely resembling their prototype AM Herculis. These binaries consist of a red dwarf and a strongly magnetic white dwarf, having orbital periods of only a few hours. Monitoring eclipses of these typically faint sources demands high-time resolution p…
▽ More
We focus on short--period eclipsing binaries that belong to a class of Cataclysmic Variables (CVs). They are known as polars and intermediate polars, closely resembling their prototype AM Herculis. These binaries consist of a red dwarf and a strongly magnetic white dwarf, having orbital periods of only a few hours. Monitoring eclipses of these typically faint sources demands high-time resolution photometry. We describe the very recent results obtained for two CVs, HU Aqr and DQ Her, which were observed with the Optical Pulsar Timing Analyzer (OPTIMA). The new observations of HU Aqr confirm that the O--C (Observed minus Calculated) diagrams exhibit variations known for this binary which can be explained by a single, massive Jupiter--like planet, possibly accompanied by a very distant companion.
△ Less
Submitted 2 October, 2012;
originally announced October 2012.
-
On the HU Aquarii planetary system hypothesis
Authors:
Krzysztof Gozdziewski,
Ilham Nasiroglu,
Aga Slowikowska,
Klaus Beuermann,
Gottfried Kanbach,
Bartosz Gauza,
Andrzej J. Maciejewski,
Robert Schwarz,
Axel D. Schwope,
Tobias C. Hinse,
Nader Haghighipour,
Vadim Burwitz,
Mariusz Slonina,
Arne Rau
Abstract:
In this work, we investigate the eclipse timing of the polar binary HU Aquarii that has been observed for almost two decades. Recently, Qian et al. attributed large (O-C) deviations between the eclipse ephemeris and observations to a compact system of two massive jovian companions. We improve the Keplerian, kinematic model of the Light Travel Time (LTT) effect and re-analyse the whole currently av…
▽ More
In this work, we investigate the eclipse timing of the polar binary HU Aquarii that has been observed for almost two decades. Recently, Qian et al. attributed large (O-C) deviations between the eclipse ephemeris and observations to a compact system of two massive jovian companions. We improve the Keplerian, kinematic model of the Light Travel Time (LTT) effect and re-analyse the whole currently available data set. We add almost 60 new, yet unpublished, mostly precision light curves obtained using the time high-resolution photo-polarimeter OPTIMA, as well as photometric observations performed at the MONET/N, PIRATE and TCS telescopes. We determine new mid--egress times with a mean uncertainty at the level of 1 second or better. We claim that because the observations that currently exist in the literature are non-homogeneous with respect to spectral windows (ultraviolet, X-ray, visual, polarimetric mode) and the reported mid--egress measurements errors, they may introduce systematics that affect orbital fits. Indeed, we find that the published data, when taken literally, cannot be explained by any unique solution. Many qualitatively different and best-fit 2-planet configurations, including self-consistent, Newtonian N-body solutions may be able to explain the data. However, using high resolution, precision OPTIMA light curves, we find that the (O-C) deviations are best explained by the presence of a single circumbinary companion orbiting at a distance of ~4.5 AU with a small eccentricity and having ~7 Jupiter-masses. This object could be the next circumbinary planet detected from the ground, similar to the announced companions around close binaries HW Vir, NN Ser, UZ For, DP Leo or SZ Her, and planets of this type around Kepler-16, Kepler-34 and Kepler-35.
△ Less
Submitted 18 May, 2012;
originally announced May 2012.
-
Testing a hypothesis of the νOctantis planetary system
Authors:
Krzysztof Gozdziewski,
Mariusz Slonina,
Cezary Migaszewski,
Anna Rozenkiewicz
Abstract:
We investigate the orbital stability of a putative Jovian planet in a compact binary νOctantis reported by Ramm et al. We re-analyzed published radial velocity data in terms of self-consistent Newtonian model and we found stable best-fit solutions that obey observational constraints. They correspond to retrograde orbits, in accord with an earlier hypothesis of Eberle & Cuntz, with apsidal lines an…
▽ More
We investigate the orbital stability of a putative Jovian planet in a compact binary νOctantis reported by Ramm et al. We re-analyzed published radial velocity data in terms of self-consistent Newtonian model and we found stable best-fit solutions that obey observational constraints. They correspond to retrograde orbits, in accord with an earlier hypothesis of Eberle & Cuntz, with apsidal lines anti-aligned with the apses of the binary. The best-fit solutions are confined to tiny stable regions of the phase space. These regions have a structure of the Arnold web formed by overlapping low-order mean motion resonances and their sub-resonances. The presence of a real planet is still questionable, because its formation would be hindered by strong dynamical perturbations. Our numerical study makes use of a new computational Message Passing Interface (MPI) framework MECHANIC developed to run massive numerical experiments on CPU clusters.
△ Less
Submitted 19 December, 2012; v1 submitted 7 May, 2012;
originally announced May 2012.
-
A dynamical analysis of the Kepler-11 planetary system
Authors:
Cezary Migaszewski,
Mariusz Slonina,
Krzysztof Gozdziewski
Abstract:
The Kepler-11 star hosts at least six transiting super-Earth planets detected through the precise photometric observations of the Kepler mission (Lissauer et al.). In this paper, we re-analyze the available Kepler data, using the direct N-body approach rather than an indirect TTV method in the discovery paper. The orbital modeling in the realm of the direct approach relies on the whole data set, n…
▽ More
The Kepler-11 star hosts at least six transiting super-Earth planets detected through the precise photometric observations of the Kepler mission (Lissauer et al.). In this paper, we re-analyze the available Kepler data, using the direct N-body approach rather than an indirect TTV method in the discovery paper. The orbital modeling in the realm of the direct approach relies on the whole data set, not only on the mid-transits times. Most of the results in the original paper are confirmed and extended. We constrained the mass of the outermost planet g to less than 30 Earth masses. The mutual inclinations between orbits b and c as well as between orbits d and e are determined with a good precision, in the range of [1,5] degrees. Having several solutions to four qualitative orbital models of the Kepler-11 system, we analyze its global dynamics with the help of dynamical maps. They reveal a sophisticated structure of the phase space, with narrow regions of regular motion. The dynamics are governed by a dense net of three- and four-body mean motion resonances, forming the Arnold web. Overlapping of these resonances is a main source of instability. We found that the Kepler-11 system may be long-term stable only in particular multiple resonant configurations with small relative inclinations. The mass-radius data derived for all companions reveal a clear anti-correlation between the mean density of the planets with their distance from the star. This may reflect the formation and early evolution history of the system.
△ Less
Submitted 13 September, 2012; v1 submitted 3 May, 2012;
originally announced May 2012.
-
International observational campaigns of the last two eclipses in EE Cephei: 2003 and 2008/9
Authors:
C. Gałan,
M. Mikołajewski,
T. Tomov,
D. Graczyk,
G. Apostolovska,
I. Barzova,
I. Bellas-Velidis,
B. Bilkina,
R. M. Blake,
C. T. Bolton,
A. Bondar,
L. Brát,
T. Brożek,
B. Budzisz,
M. Cikała,
B. Csák,
A. Dapergolas,
D. Dimitrov,
P. Dobierski,
M. Drahus,
M. Dróżdż,
S. Dvorak,
L. Elder,
S. Frcakowiak,
G. Galazutdinov
, et al. (65 additional authors not shown)
Abstract:
Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in…
▽ More
Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete set of observational data collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous with a "bump" that is very clear in all the UBVRI light curves. Variations in the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and strong absorption is superimposed on the H_alpha emission. These observations confirm that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the "bump" observed during the last two eclipses in terms of the different times of obscuration of the hot polar regions of the Be star by the disk; and (iii) our suggested period of the disk precession (~11-12 P_orb) and predicted depth of about 2 mag the forthcoming eclipse in 2014.
△ Less
Submitted 26 September, 2012; v1 submitted 30 April, 2012;
originally announced May 2012.
-
Mechanic: a new numerical MPI framework for the dynamical astronomy
Authors:
Mariusz Slonina,
Krzysztof Gozdziewski,
Cezary Migaszewski
Abstract:
We develop the Mechanic package, which is a new numerical framework for dynamical astronomy. The aim of our software is to help in massive numerical simulations by efficient task management and unified data storage. The code is built on top of the Message Passing Interface (MPI) and Hierarchical Data Format (HDF5) standards and uses the Task Farm approach to manage numerical tasks. It relies on th…
▽ More
We develop the Mechanic package, which is a new numerical framework for dynamical astronomy. The aim of our software is to help in massive numerical simulations by efficient task management and unified data storage. The code is built on top of the Message Passing Interface (MPI) and Hierarchical Data Format (HDF5) standards and uses the Task Farm approach to manage numerical tasks. It relies on the core-module approach. The numerical problem implemented in the user-supplied module is separated from the host code (core). The core is designed to handle basic setup, data storage and communication between nodes in a computing pool. It has been tested on large CPU-clusters, as well as desktop computers. The Mechanic may be used in computing dynamical maps, data optimization or numerical integration. The code and sample modules are freely available at http://git.astri.umk.pl/projects/mechanic.
△ Less
Submitted 29 February, 2012;
originally announced February 2012.
-
The Proposed Quadruple System SZ Herculis: Revised LITE Model and Orbital Stability Study
Authors:
T. C. Hinse,
K. Gozdziewski,
J. W. Lee,
N. Haghighipour,
C. -U. Lee
Abstract:
In a recent study, Lee et al. presented new photometric follow-up timing observations of the semi-detached binary system SZ Herculis and proposed the existence of two hierarchical cirumbinary companions. Based on the light-travel time effect, the two low-mass M-dwarf companions are found to orbit the binary pair on moderate to high eccentric orbits. The derived periods of these two companions are…
▽ More
In a recent study, Lee et al. presented new photometric follow-up timing observations of the semi-detached binary system SZ Herculis and proposed the existence of two hierarchical cirumbinary companions. Based on the light-travel time effect, the two low-mass M-dwarf companions are found to orbit the binary pair on moderate to high eccentric orbits. The derived periods of these two companions are close to a 2:1 mean-motion orbital resonance. We have studied the stability of the system using the osculating orbital elements as presented by Lee et al. Results indicate an orbit-crossing architecture exhibiting short-term dynamical instabilities leading to the escape of one of the proposed companions. We have examined the system's underlying model parameter-space by following a Monte Carlo approach and found an improved fit to the timing data. A study of the stability of our best-fitting orbits also indicates that the proposed system is generally unstable. If the observed anomalous timing variations of the binary period is due to additional circumbinary companions, then the resulting system should exhibit a long-term stable orbital configuration much different from the orbits suggested by Lee et al. We, therefore, suggest that based on Newtonian-dynamical considerations, the proposed quadruple system cannot exist. To uncover the true nature of the observed period variations of this system, we recommend future photometric follow-up observations that could further constrain eclipse-timing variations and/or refine light-travel time models.
△ Less
Submitted 21 February, 2012;
originally announced February 2012.
-
New Light-Travel Time Models and Orbital Stability Study of the Proposed Planetary System HU Aquarii
Authors:
Tobias C. Hinse,
Jae Woo Lee,
Krzysztof Gozdziewski,
Nader Haghighipour,
Chung-Uk Lee,
Eamon M. Scullion
Abstract:
In this work we propose a new orbital architecture for the two proposed circumbinary planets around the polar eclipsing binary HU Aquarii. We base the new two-planet, light-travel time model on the result of a Monte Carlo simulation driving a least-squares Levenberg-Marquardt minimisation algorithm on the observed eclipse egress times. Our best-fitting model with $χ_{r}^2=1.43$ resulted in high fi…
▽ More
In this work we propose a new orbital architecture for the two proposed circumbinary planets around the polar eclipsing binary HU Aquarii. We base the new two-planet, light-travel time model on the result of a Monte Carlo simulation driving a least-squares Levenberg-Marquardt minimisation algorithm on the observed eclipse egress times. Our best-fitting model with $χ_{r}^2=1.43$ resulted in high final eccentricities for the two companions leading to an unstable orbital configuration. From a large ensemble of initial guesses we examined the distribution of final eccentricities and semi-major axes for different $χ_{r}^2$ parameter intervals and encountered qualitatively a second population of best-fitting parameters. The main characteristic of this population is described by low-eccentric orbits favouring long-term orbital stability of the system. We present our best-fitting model candidate for the proposed two-planet system and demonstrate orbital stability over one million years using numerical integrations.
△ Less
Submitted 30 November, 2011;
originally announced December 2011.
-
Modeling the Radial Velocities of HD 240210 with the Genetic Algorithms
Authors:
A. Rozenkiewicz,
K. Gozdziewski,
C. Migaszewski
Abstract:
More than 450 extrasolar planets are known to date. To detect these intriguing objects, many photometric and radial velocity (RV) surveys are in progress. We developed the Keplerian FITting (KFIT) code, to model published and available RV data. KFIT is based on a hybrid, quasi-global optimization technique relying on the Genetic Algorithms and simplex algorithm. Here, we re-analyse the RV data of…
▽ More
More than 450 extrasolar planets are known to date. To detect these intriguing objects, many photometric and radial velocity (RV) surveys are in progress. We developed the Keplerian FITting (KFIT) code, to model published and available RV data. KFIT is based on a hybrid, quasi-global optimization technique relying on the Genetic Algorithms and simplex algorithm. Here, we re-analyse the RV data of evolved K3III star HD 240210. We confirm three equally good solutions which might be interpreted as signals of 2-planet systems. Remarkably, one of these best-fits describes the long-term stable two-planet system, involved in the 2:1~mean motion resonance (MMR). It may be the first instance of this strong MMR in a multi-planet system hosted by evolved star, as the 2:1 MMR configurations are already found around a few solar dwarfs.
△ Less
Submitted 30 June, 2010;
originally announced July 2010.
-
The non-resonant, relativistic dynamics of circumbinary planets
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski
Abstract:
We investigate the non-resonant, 3-D (spatial) model of the hierarchical system composed of point-mass stellar (or sub-stellar) binary and a low-mass companion (a circumbinary planet or a brown dwarf). We take into account the leading relativistic corrections to the Newtonian gravity. The secular model of the system relies on the expansion of the perturbing Hamiltonian in terms of the ratio of sem…
▽ More
We investigate the non-resonant, 3-D (spatial) model of the hierarchical system composed of point-mass stellar (or sub-stellar) binary and a low-mass companion (a circumbinary planet or a brown dwarf). We take into account the leading relativistic corrections to the Newtonian gravity. The secular model of the system relies on the expansion of the perturbing Hamiltonian in terms of the ratio of semi-major axes $α$, averaged over the mean anomalies. We found that the low-mass object in a distant orbit may excite large eccentricity of the inner binary when the mutual inclination of the orbits is larger than about of 60 deg. This is related to strong instability caused by a phenomenon which acts similarly to the Lidov-Kozai resonance (LKR). The secular system may be strongly chaotic and its dynamics unpredictable over the long-term time scale. Our study shows that in the Jupiter-- or brown dwarf-- mass regime of the low-massive companion, the restricted model does not properly describe the long-term dynamics in the vicinity of the LKR. The relativistic correction is significant for the parametric structure of a few families of stationary solutions in this problem, in particular, for the direct orbits configurations (with the mutual inclination less than 90 degrees). We found that the dynamics of hierarchical systems with small $α\sim 0.01$ may be qualitatively different in the realm of the Newtonian (classic) and relativistic models. This holds true even for relatively large masses of the secondaries.
△ Less
Submitted 19 October, 2010; v1 submitted 30 June, 2010;
originally announced June 2010.
-
The Calan-Hertfordshire Extrasolar Planet Search
Authors:
J. S. Jenkins,
H. R. A. Jones,
K. Gozdziewski,
C. Migaszewski,
J. R. Barnes,
M. I. Jones,
P. Rojo,
D. J. Pinfield,
A. C. Day-Jones,
S. Hoyer
Abstract:
The detailed study of the exoplanetary systems HD189733 and HD209458 has given rise to a wealth of exciting information on the physics of exoplanetary atmospheres. To further our understanding of the make-up and processes within these atmospheres we require a larger sample of bright transiting planets. We have began a project to detect more bright transiting planets in the southern hemisphere by…
▽ More
The detailed study of the exoplanetary systems HD189733 and HD209458 has given rise to a wealth of exciting information on the physics of exoplanetary atmospheres. To further our understanding of the make-up and processes within these atmospheres we require a larger sample of bright transiting planets. We have began a project to detect more bright transiting planets in the southern hemisphere by utilising precision radial-velocity measurements. We have observed a constrained sample of bright, inactive and metal-rich stars using the HARPS instrument and here we present the current status of this project, along with our first discoveries which include a brown dwarf/extreme-Jovian exoplanet found in the brown dwarf desert region around the star HD191760 and improved orbits for three other exoplanetary systems HD48265, HD143361 and HD154672. Finally, we briefly discuss the future of this project and the current prospects we have for discovering more bright transiting planets.
△ Less
Submitted 29 December, 2009;
originally announced December 2009.
-
Relative equilibria in the unrestricted problem of a sphere and symmetric rigid body
Authors:
Mikhail Vereshchagin,
Andrzej J. Maciejewski,
Krzysztof Gozdziewski
Abstract:
We consider the unrestricted problem of two mutually attracting rigid bodies, an uniform sphere (or a point mass) and an axially symmetric body. We present a global, geometric approach for finding all relative equilibria (stationary solutions) in this model, which was already studied by Kinoshita (1970). We extend and generalize his results, showing that the equilibria solutions may be found by…
▽ More
We consider the unrestricted problem of two mutually attracting rigid bodies, an uniform sphere (or a point mass) and an axially symmetric body. We present a global, geometric approach for finding all relative equilibria (stationary solutions) in this model, which was already studied by Kinoshita (1970). We extend and generalize his results, showing that the equilibria solutions may be found by solving at most two non-linear, algebraic equations, assuming that the potential function of the symmetric rigid body is known explicitly. We demonstrate that there are three classes of the relative equilibria, which we call "cylindrical", "inclined co-planar", and "conic" precessions, respectively. Moreover, we also show that in the case of conic precession, although the relative orbit is circular, the point-mass and the mass center of the body move in different parallel planes. This solution has been yet not known in the literature.
△ Less
Submitted 11 December, 2009;
originally announced December 2009.
-
First Results from the CHEPS: Exoplanets and the Discovery of an Eccentric Brown Dwarf in the Desert
Authors:
James S. Jenkins,
Hugh R. A. Jones,
Krzysztof Gozdziewski,
Cezary Migaszewski,
John R. Barnes,
Matias I. Jones,
Patricio Rojo,
David J. Pinfield,
Avril C. Day-Jones,
Sergio Hoyer
Abstract:
We report the discovery of a brown dwarf on an eccentric orbit and with a semimajor axis that places it in the brown dwarf desert region around the star HD191760. The star has a spectral type of G3IV/V and a metallicity ([Fe/H]) of 0.29 dex. HD191760 adds to the small number of metal-rich stars with brown dwarf companions. The brown dwarf (HD191760b) is found to have an orbital period of 505.57+…
▽ More
We report the discovery of a brown dwarf on an eccentric orbit and with a semimajor axis that places it in the brown dwarf desert region around the star HD191760. The star has a spectral type of G3IV/V and a metallicity ([Fe/H]) of 0.29 dex. HD191760 adds to the small number of metal-rich stars with brown dwarf companions. The brown dwarf (HD191760b) is found to have an orbital period of 505.57+/-0.40 days and semimajor axis of 1.35+/-0.01 AU, placing it firmly in the brown dwarf desert. The eccentricity of HD191760b is found to be 0.63+/-0.01, meaning it reaches as close as 0.5 AU from the host star. Dynamical simulations indicate that no inner planets could reside at separations beyond ~0.17 AU due to the disastrous gravity imposed by HD191760b. In addition to these first results we also refine the orbits found for the exoplanets around the stars HD48265, HD143361 and HD154672. All 1-planet solutions are in agreement with those previously published by the Magellan Planet Search.
△ Less
Submitted 18 May, 2009;
originally announced May 2009.
-
Is the HR 8799 extrasolar system destined for planetary scattering?
Authors:
Krzysztof Gozdziewski,
Cezary Migaszewski
Abstract:
The recent discovery of a three-planet extrasolar system of HR 8799 by Marois et al. is a breakthrough in the field of the direct imaging. This great achievement raises questions on the formation and dynamical stability of the HR 8799 system, because Keplerian fits to astrometric data are strongly unstable during ~0.2Myr. We search for stable, self-consistent N-body orbits with the so called GAM…
▽ More
The recent discovery of a three-planet extrasolar system of HR 8799 by Marois et al. is a breakthrough in the field of the direct imaging. This great achievement raises questions on the formation and dynamical stability of the HR 8799 system, because Keplerian fits to astrometric data are strongly unstable during ~0.2Myr. We search for stable, self-consistent N-body orbits with the so called GAMP method that incorporates stability constraints into the optimization algorithm. Our searches reveal only small regions of stable motions in the phase space of three-planet, coplanar configurations. Most likely, if the planetary masses are in 10-Jupiter-mass range, they may be stable only if the planets are involved in two- or three-body mean motion resonances (MMRs). We found that 80% systems found by GAMP that survived 30Myr backwards integrations, eventually become unstable after 100Myr. It could mean that the HR 8799 system undergo a phase of planet-planet scattering. We test a hypothesis that the less certain detection of the innermost object is due to a blending effect. In such a case, two-planet best-fit systems are mostly stable, on quasi-circular orbits and close to the 5:2 MMR, resembling the Jupiter-Saturn pair.
△ Less
Submitted 27 April, 2009;
originally announced April 2009.
-
Relativistic Lidov-Kozai resonance in binaries
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski
Abstract:
We consider the secular dynamics of a binary and a planet in terms of non-restricted, hierarchical three-body problem, including the general relativity corrections to the Newtonian gravity. We determine regions in the parameter space where the relativistic corrections may be important for the long-term dynamics. We try to constrain the inclinations of putative Jovian planets in recently announce…
▽ More
We consider the secular dynamics of a binary and a planet in terms of non-restricted, hierarchical three-body problem, including the general relativity corrections to the Newtonian gravity. We determine regions in the parameter space where the relativistic corrections may be important for the long-term dynamics. We try to constrain the inclinations of putative Jovian planets in recently announced binary systems of HD 4113 and HD 156846.
△ Less
Submitted 31 December, 2008;
originally announced January 2009.
-
Equilibria in the secular, non-coplanar two-planet problem
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski
Abstract:
We investigate the secular dynamics of a planetary system composed of the parent star and two massive planets in mutually inclined orbits. The dynamics are investigated in wide ranges of semi-major axes ratios (0.1-0.667), and planetary masses ratios (0.25-2) as well as in the whole permitted ranges of the energy and total angular momentum. The secular model is constructed by semi-analytic avera…
▽ More
We investigate the secular dynamics of a planetary system composed of the parent star and two massive planets in mutually inclined orbits. The dynamics are investigated in wide ranges of semi-major axes ratios (0.1-0.667), and planetary masses ratios (0.25-2) as well as in the whole permitted ranges of the energy and total angular momentum. The secular model is constructed by semi-analytic averaging of the three-body system. We focus on equilibria of the secular Hamiltonian (periodic solutions of the full system), and we analyze their stability. We attempt to classify families of these solutions in terms of the angular momentum integral. We identified new equilibria, yet unknown in the literature. Our results are general and may be applied to a wide class of three-body systems, including configurations with a star and brown dwarfs and sub-stellar objects. We also describe some technical aspects of the semi-numerical averaging. The HD 12661 planetary system is investigated as an example configuration.
△ Less
Submitted 2 February, 2009; v1 submitted 15 December, 2008;
originally announced December 2008.
-
On the habitability of the OGLE-2006-BLG-109L planetary system
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski,
Tobias C. Hinse
Abstract:
We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analog of the Solar system. Our work is inspired by work of Malhotra and Minton (2008). Using the linear Laplace-Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to unce…
▽ More
We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analog of the Solar system. Our work is inspired by work of Malhotra and Minton (2008). Using the linear Laplace-Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the sub-system of Jupiters may be found in dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semi-major axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.
△ Less
Submitted 2 February, 2009; v1 submitted 10 November, 2008;
originally announced November 2008.
-
Secular dynamics of coplanar, non-resonant planetary system under the general relativity and quadrupole moment perturbations
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski
Abstract:
We construct a secular theory of a coplanar system of N-planets not involved in strong mean motion resonances, and which are far from collision zones. Besides the point-to-point Newtonian mutual interactions, we consider the general relativity corrections to the gravitational potential of the star and the innermost planet, and also a modification of this potential by the quadrupole moment and ti…
▽ More
We construct a secular theory of a coplanar system of N-planets not involved in strong mean motion resonances, and which are far from collision zones. Besides the point-to-point Newtonian mutual interactions, we consider the general relativity corrections to the gravitational potential of the star and the innermost planet, and also a modification of this potential by the quadrupole moment and tidal distortion of the star. We focus on hierarchical planetary systems. A survey regarding model parameters (the masses, semi-major axes, spin rate of the star) reveals a rich and non-trivial dynamics of the secular system. Our study is focused on its equilibria. Such solutions predicted by the classic secular theory, which correspond to aligned (mode I) or anti-aligned (mode II) apsides, may be strongly affected by the gravitational corrections. The so called true secular resonance, which is a new feature of the classic two-planet problem discovered by Michtchenko & Malhotra (2004), may appear in other, different regions of the phase space of the generalized model. We found bifurcations of mode II which emerge new, yet unknown in the literature, secularly unstable equilibria and a complex structure of the phase space. These equilibria may imply secularly unstable orbital configurations even for nitially moderate eccentricities. The point mass gravity corrections can affect the long term-stability in the secular time scale, which may directly depend on the age of the host star through its spin rate. We also analyze the secular dynamics of the upsilon Andromede system in the realm of the generalized model. Also in this case of the three-planet system, new secular equilibria may appear.
△ Less
Submitted 30 September, 2008;
originally announced September 2008.
-
The Successful Prediction of the Extrasolar Planet HD 74156 d
Authors:
Rory Barnes,
Krzysztof Gozdziewski,
Sean N. Raymond
Abstract:
Most of the first-discovered extrasolar multi-planet systems were found to lie close to dynamically unstable configurations. However a few observed multi-planet systems (e.g. HD 74156) did not show this trait. Those systems could share this property if they contain an additional planet in between those that are known. Previous investigations identified the properties of hypothetical planets that…
▽ More
Most of the first-discovered extrasolar multi-planet systems were found to lie close to dynamically unstable configurations. However a few observed multi-planet systems (e.g. HD 74156) did not show this trait. Those systems could share this property if they contain an additional planet in between those that are known. Previous investigations identified the properties of hypothetical planets that would place these systems near instability. The hypothetical planet in HD 74156 was expected to have a mass about equal to that of Saturn, a semi-major axis between 0.9 and 1.4 AU, and an eccentricity less than 0.2. HD 74156 d, a planet with a mass of 1.3 Saturn masses at 1.04 AU with an eccentricity of 0.25, was recently announced. We have reanalyzed all published data on this system in order to place tighter constraints on the properties of the new planet. We find two possible orbits for this planet, one close to that already identified and another (with a slightly better fit to the data) at ~0.89 AU. We also review the current status of other planet predictions, discuss the observed single planet systems, and suggest other systems which may contain planets in between those that are already known. The confirmation of the existence of HD 74156 d suggests that planet formation is an efficient process, and planetary systems should typically contain many planets.
△ Less
Submitted 28 April, 2008;
originally announced April 2008.
-
A secular theory of coplanar, non-resonant planetary system
Authors:
Cezary Migaszewski,
Krzysztof Gozdziewski
Abstract:
We present the secular theory of coplanar $N$-planet system, in the absence of mean motion resonances between the planets. This theory relies on the averaging of a perturbation to the two-body problem over the mean longitudes. We expand the perturbing Hamiltonian in Taylor series with respect to the ratios of semi-major axes which are considered as small parameters, without direct restrictions o…
▽ More
We present the secular theory of coplanar $N$-planet system, in the absence of mean motion resonances between the planets. This theory relies on the averaging of a perturbation to the two-body problem over the mean longitudes. We expand the perturbing Hamiltonian in Taylor series with respect to the ratios of semi-major axes which are considered as small parameters, without direct restrictions on the eccentricities. Next, we average out the resulting series term by term. This is possible thanks to a particular but in fact quite elementary choice of the integration variables. It makes it possible to avoid Fourier expansions of the perturbing Hamiltonian. We derive high order expansions of the averaged secular Hamiltonian (here, up to the order of 24) with respect to the semi-major axes ratio. The resulting secular theory is a generalization of the octupole theory. The analytical results are compared with the results of numerical (i.e., practically exact) averaging. We estimate the convergence radius of the derived expansions, and we propose a further improvement of the algorithm. As a particular application of the method, we consider the secular dynamics of three-planet coplanar system. We focus on stationary solutions in the HD 37124 planetary system.
△ Less
Submitted 9 May, 2008; v1 submitted 22 March, 2008;
originally announced March 2008.
-
Stability constraints in modeling of multi-planet extrasolar systems
Authors:
K. Gozdziewski,
C. Migaszewski,
A. Musielinski
Abstract:
We present an analysis of high precision radial velocity (RV) observations of stars hosting multi-planet systems with Jovian companions. We use dynamical stability constraints and quasi-global methods of optimization. As an illustration, we present new results derived for the RV data of the Sun-like dwarfs HD 155358 and $τ^1$ Gruis.
We present an analysis of high precision radial velocity (RV) observations of stars hosting multi-planet systems with Jovian companions. We use dynamical stability constraints and quasi-global methods of optimization. As an illustration, we present new results derived for the RV data of the Sun-like dwarfs HD 155358 and $τ^1$ Gruis.
△ Less
Submitted 2 February, 2008;
originally announced February 2008.
-
The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects
Authors:
Krzysztof Gozdziewski,
Slawomir Breiter,
Wojciech Borczyk
Abstract:
We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincare variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator MEGNO (a measure of the maximal Lyapunov exponent) that helps to distinguish chaotic and regular configur…
▽ More
We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincare variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator MEGNO (a measure of the maximal Lyapunov exponent) that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best fit solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best fit solutions are found in dynamically active region of the phase space. The long term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myrs of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behavior.
△ Less
Submitted 29 October, 2007;
originally announced October 2007.
-
A dynamical analysis of the 14 Her planetary system
Authors:
K. Gozdziewski,
C. Migaszewski,
M. Konacki
Abstract:
Precision radial velocity (RV) measurements of the Sun-like dwarf
14 Herculis published by Naef et. al (2004), Butler et. al (2006) and
Wittenmyer et al (2007) reveal a Jovian planet in a 1760 day orbit and a trend indicating the second distant object. On the grounds of dynamical considerations, we test a hypothesis that the trend can be explained by the presence of an additional giant plane…
▽ More
Precision radial velocity (RV) measurements of the Sun-like dwarf
14 Herculis published by Naef et. al (2004), Butler et. al (2006) and
Wittenmyer et al (2007) reveal a Jovian planet in a 1760 day orbit and a trend indicating the second distant object. On the grounds of dynamical considerations, we test a hypothesis that the trend can be explained by the presence of an additional giant planet. We derive dynamical limits to th orbital parameters of the putative outer Jovian companion in an orbit within ~13 AU. In this case, the mutual interactions between the Jovian planets are important for the long-term stability of the system. The best self-consistent and stable Newtonian fit to an edge-on configuration of Jovian planets has the outer planet in 9 AU orbit with a moderate eccentricity ~0.2 and confined to a zone spanned by the low-order mean motion resonances 5:1 and 6:1. This solution lies in a shallow minimum of \Chi and persists over a wide range of the system inclination. Other stable configurations within 1σconfidence interval of the best fit are possible for the semi-major axis of the outer planet in the range of (6,13) AU and the eccentricity in the range of (0,0.3). The orbital inclination cannot yet be determined but when it decreases, both planetary masses approach ~10 Jupiter masses and for ~30 deg the hierarchy of the masses is reversed.
△ Less
Submitted 5 January, 2008; v1 submitted 13 May, 2007;
originally announced May 2007.
-
On the extrasolar multi-planet system around HD160691
Authors:
Krzysztof Gozdziewski,
Andrzej J. Maciejewski,
Cezary Migaszewski
Abstract:
We re-analyze the precision radial velocity (RV) observations of HD160691 (mu Ara) by the Anglo-Australian Planet Search Team. The star is supposed to host two Jovian companions (HD160691b, HD160691c) in long-period orbits about 630 days and about 2500 days, respectively) and a hot-Neptune (HD160691d) in about 9 days orbit. We perform a global search for the best fits in the orbital parameters s…
▽ More
We re-analyze the precision radial velocity (RV) observations of HD160691 (mu Ara) by the Anglo-Australian Planet Search Team. The star is supposed to host two Jovian companions (HD160691b, HD160691c) in long-period orbits about 630 days and about 2500 days, respectively) and a hot-Neptune (HD160691d) in about 9 days orbit. We perform a global search for the best fits in the orbital parameters space with a hybrid code employing the genetic algorithm and simplex method. The stability of Keplerian fits is verified with the N-body model of the RV signal that takes into account the dynamical constraints (so called GAMP method). Our analysis reveals a signature of the fourth, yet unconfirmed, Jupiter-like planet HD160691e in about 307 days orbit. Overall, the global architecture of four-planet configuration recalls the Solar system. All companions of mu Ara move in quasi-circular orbits. The orbits of two inner Jovian planets are close to the 2:1 mean motion resonance. The alternative three-planet system involves two Jovian planets in eccentric orbits (e about 0.3), close to the 4:1 MMR, but it yields a significantly worse fit to the data. We also verify a hypothesis of the 1:1 MMR in the subsystem of two inner Jovian planets in the four-planet model.
△ Less
Submitted 12 October, 2006; v1 submitted 14 August, 2006;
originally announced August 2006.
-
Planetary nebulae with emission-line central stars
Authors:
K. Gesicki,
A. A. Zijlstra,
A. Acker,
S. K. Gorny,
K. Gozdziewski,
J. R. Walsh
Abstract:
The kinematic structure of a sample of planetary nebulae, consisting of 23 [WR] central stars, 21 weak emission line stars (wels) and 57 non-emission line central stars, is studied. The [WR] stars are shown to be surrounded by turbulent nebulae, a characteristic shared by some wels but almost completely absent from the non-emission line stars. The fraction of objects showing turbulence for non-e…
▽ More
The kinematic structure of a sample of planetary nebulae, consisting of 23 [WR] central stars, 21 weak emission line stars (wels) and 57 non-emission line central stars, is studied. The [WR] stars are shown to be surrounded by turbulent nebulae, a characteristic shared by some wels but almost completely absent from the non-emission line stars. The fraction of objects showing turbulence for non-emission-line stars, wels and [WR] stars is 7%, 24% and 91%, respectively. The [WR] stars show a distinct IRAS 12-micron excess, indicative of small dust grains, which is not found for wels. The [WR]-star nebulae are on average more centrally condensed than those of other stars. On the age-temperature diagram, the wels are located on tracks of both high and low stellar mass, while [WR] stars trace a narrow range of intermediate masses. Emission-line stars are not found on the cooling track. One group of wels may form a sequence wels--[WO] stars with increasing temperature. For the other groups both the wels and the [WR] stars appear to represent several, independent evolutionary tracks. We find a discontinuity in the [WR] stellar temperature distribution and suggest different evolutionary sequences above and below the temperature gap. One group of cool [WR] stars has no counterpart among any other group of PNe and may represent binary evolution. A prime factor distinguishing wels and [WR] stars appears to be stellar luminosity. We find no evidence for an increase of nebular expansion velocity with time.
△ Less
Submitted 20 January, 2006; v1 submitted 19 January, 2006;
originally announced January 2006.