Source author record

Cezary Migaszewski

Cezary Migaszewski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2016arXiv

A resonant chain of four transiting, sub-Neptune planets

Surveys have revealed many multi-planet systems containing super-Earths and Neptunes in orbits of a few days to a few months. There is debate whether in situ assembly or inward migration is the dominant mechanism of the formation of such planetary systems. Simulations suggest that migration creates tightly packed systems with planets whose orbital periods may be expressed as ratios of small integers (resonances), often in a many-planet series (chain). In the hundreds of multi-planet systems of sub-Neptunes, more planet pairs are observed near resonances than would generally be expected, but no individual system has hitherto been identified that must have been formed by migration. Proximity to resonance enables the detection of planets perturbing each other. Here we report transit timing variations of the four planets in the Kepler-223 system, model these variations as resonant-angle librations, and compute the long-term stability of the resonant chain. The architecture of Kepler-223 is too finely tuned to have been formed by scattering, and our numerical simulations demonstrate that its properties are natural outcomes of the migration hypothesis. Similar systems could be destabilized by any of several mechanisms, contributing to the observed orbital-period distribution, where many planets are not in resonances. Planetesimal interactions in particular are thought to be responsible for establishing the current orbits of the four giant planets in the Solar System by disrupting a theoretical initial resonant chain similar to that observed in Kepler-223.

preprint2016arXiv

On the migration of three planets in a protoplanetary disc and the formation of chains of mean motion resonances

We study the migration of three-planet systems in an irradiated 1+1D $α$-disc with photoevaporation. We performed $2700$ simulations with various planets' masses and initial orbits. We found that most of the systems which ended up as compact configurations form chains of mean motion resonances (MMRs) of the first and higher orders. Most of the systems involved in chains of MMRs are periodic configurations. The period ratios of such system, though, are not necessarily close to exact commensurability. If a given system resides in a divergent migration zone in the disc, the period ratios increase and evolve along resonant divergent migration paths at ($P_2/P_1, P_3/P_2)$-diagram, where $P_1, P_2, P_3$ are the orbital periods of the first, second and third planet, respectively. The observed systems, though, do not lie on those paths. We show that an agreement between the synthetic and the observed systems distributions could be achieved if the orbital circularization was slower than it results from models of the planet-disc interactions. Therefore, we conclude that most of those systems unlikely formed as a result of divergent migration out of nominal chains of MMRs.

preprint2015arXiv

On the migration of two planets in a disc and the formation of mean motion resonances

We study the dynamics of a system of two super-Earths embedded in a protoplanetary disc. We build a simple model of an irradiated viscous disc and use analytical prescriptions for the planet-disc interactions which lead to migration. We show that depending on the disc parameters, planets' masses and their positions in the disc, the migration of each planet can be inward or outward and the migration of a two-planet system can be convergent (which may lead to formation of a first order mean motion resonance, MMR) or divergent (a system moves away from MMR). We performed $3500$ simulations of the migration of two-planet systems with various masses and initial orbits. Almost all of them end up as resonant configurations, although the period ratios may be very distant from the nominal values of a given MMR. We found that almost all the systems resulting from the migration are periodic configurations.

preprint2015arXiv

The Laplace resonance in the Kepler-60 system

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.

preprint2014arXiv

Mechanic: the MPI/HDF code framework for dynamical astronomy

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).

preprint2014arXiv

Multiple mean motion resonances in the HR 8799 planetary system

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.

preprint2013arXiv

A linear distribution of orbits in compact planetary systems?

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.

preprint2012arXiv

A dynamical analysis of the Kepler-11 planetary system

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.

preprint2012arXiv

Mechanic: a new numerical MPI framework for the dynamical astronomy

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.

preprint2012arXiv

Testing a hypothesis of the νOctantis planetary system

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.

preprint2012arXiv

The generalized non-conservative model of a 1-planet system - revisited

We study the long-term dynamics of a planetary system composed of a star and a planet. Both bodies are considered as extended, non-spherical, rotating objects. There are no assumptions made on the relative angles between the orbital angular momentum and the spin vectors of the bodies. Thus, we analyze full, spatial model of the planetary system. Both objects are assumed to be deformed due to their own rotations, as well as due to the mutual tidal interactions. The general relativity corrections are considered in terms of the post-Newtonian approximation. Besides the conservative contributions to the perturbing forces, there are also taken into account non-conservative effects, i.e., the dissipation of the mechanical energy. This dissipation is a result of the tidal perturbation on the velocity field in the internal zones with non-zero turbulent viscosity (convective zones). Our main goal is to derive the equations of the orbital motion as well as the equations governing time-evolution of the spin vectors (angular velocities). We derive the Lagrangian equations of the second kind for systems which do not conserve the mechanical energy. Next, the equations of motion are averaged out over all fast angles with respect to time-scales characteristic for conservative perturbations. The final equations of motion are then used to study the dynamics of the non-conservative model over time scales of the order of the age of the star. We analyze the final state of the system as a function of the initial conditions. Equilibria states of the averaged system are finally discussed.

preprint2010arXiv

The non-resonant, relativistic dynamics of circumbinary planets

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.

preprint2009arXiv

First Results from the CHEPS: Exoplanets and the Discovery of an Eccentric Brown Dwarf in the Desert

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.

preprint2009arXiv

Is the HR 8799 extrasolar system destined for planetary scattering?

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.

preprint2009arXiv

Modeling the RV and BVS of active stars

We present a method of modeling the radial velocity (RV) measurements which can be useful in searching for planets hosted by chromospherically active stars. We assume that the observed RV signal is induced by the reflex motion of a star as well as by distortions of spectral line profiles, measured by the Bisector Velocity Span (BVS). The RVs are fitted with a common planetary model including RV correction term depending linearly on the BVS, which accounts for the stellar activity. The coefficient of correlation is an additional free parameter of the RV model. That approach differs from correcting the RVs before or after fitting the "pure" planetary model. We test the method on simulated data derived for single-planet systems. The results are compared with the outcomes of algorithms found in the literature.

preprint2008arXiv

Relativistic Lidov-Kozai resonance in binaries

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.