Source author record

Valeri V. Makarov

Valeri V. Makarov 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

27works
6topics
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

27 published item(s)

preprint2025arXiv

Distributions of wide binary stars in theory and in Gaia data: III. Orbital momenta, masses, and manifestations of MOND

Using the censored catalog of 103,169 resolved Gaia DR3 binary stars with accurate astrometric data for each component, a new observable, object-specific parameter is computed for each pair: the projected orbital momentum. This parameter is the product of four functions of physical characteristics: total mass, semimajor axis, eccentricity, and inclination angle. Using the previously estimated marginal probability densities of eccentricity and semimajor axis, and assuming an isotropic orientation of binary systems, the sample distribution of mass was adjusted using a concordance metric of the observed and synthetic distributions of orbital momenta and an ad hoc functional model. The best-fitting mass density model is found to faithfully reproduce the observed dependence of orbital momenta on apparent separation, although the absolute luminosity distributions indicate a tendency of the widest systems to more frequently include solar-type primaries. The anticipated manifestation of MOND is computed in the investigated parameter space \{separation, momentum\}. This effect is absent in the given data. The median total mass of the widest Gaia binaries is found to be somewhat higher than that of the tighter pairs, which is interpreted as a dynamical age effect.

preprint2022arXiv

Chaos over Order: Mapping 3D Rotation of Triaxial Asteroids and Minor Planets

Celestial bodies approximated with rigid triaxial ellipsoids in a two-body system can rotate chaotically due to the time-varying gravitational torque from the central mass. At small orbital eccentricity values, rotation is short-term orderly and predictable within the commensurate spin-orbit resonances, while at eccentricity approaching unity, chaos completely takes over. Here, we present the full 3D rotational equations of motion around all three principle axes for triaxial minor planets and two independent methods of numerical solution based on Euler rotations and quaternion algebra. The domains of chaotic rotation are numerically investigated over the entire range of eccentricity with a combination of trial integrations of Euler's equations of motion and the GALI($k$) method. We quantify the dependence of the order--chaos boundaries on shape by changing a prolateness parameter, and find that the main 1:1 spin-orbit resonance disappears for specific moderately prolate shapes already at eccentricities as low as 0.3. The island of short-term stability around the main 1:1 resonance shrinks with increasing eccentricity at a fixed low degree of prolateness and completely vanishes at approximately 0.8. This island is also encroached by chaos on longer time scales indicating longer Lyapunov exponents. Trajectories in the close vicinity of the 3:2 spin-orbit resonance become chaotic at smaller eccentricities, but separated enclaves of orderly rotation emerge at eccentricities as high as 0.8. Initial perturbations of rotational velocity in latitude away from the exact equilibrium result in a spectrum of free libration, nutation, and polar wander, which is not well matched by the linearized analysis omitting the inertial terms.

preprint2022arXiv

Quasars with Proper Motions and the Link to Double and Multiple AGNs

Gaia used a large sample of photometrically selected active galactic nuclei (AGNs) and quasars to remove the residual spin of its global proper motion system in order to achieve a maximally inertial reference frame. A small fraction of these reference objects have statistically significant astrometric proper motions in Gaia EDR3. We compile a source sample of $105,593$ high-fidelity AGNs with accurate spectroscopically determined redshifts above 0.5 from the SDSS and normalized proper motions below 4. The rate of genuinely perturbed proper motions is at least 0.17\%. A smaller high completeness sample of 152 quasars with excess proper motions at a confidence level of 0.9995 is examined in detail. Pan-STARRS images and Gaia-resolved pairs reveal that 29\% of the sample are either double sources or gravitationally lensed quasars. An Anderson--Darling test on parameters of a smaller high-reliability sample and their statistical controls reveals 17 significant factors that favor multiplicity and multi-source structure as the main cause of perturbed astrometry. Using a nearest neighbor distance statistical analysis and counts of close companions in Gaia on a much larger initial sample of AGNs, an excess of closely separated sources in Gaia is detected. At least 0.33\% of all optical quasars are genuinely double or multiply imaged. We provide a list of 44 candidate double or multiple AGNs and four previously known gravitational lenses. Many proper motion quasars may be more closely separated, unresolved doubles exhibiting the variability imposed motion (VIM) effect, and a smaller fraction may be chance alignments with foreground stars causing weak gravitational lensing.

preprint2022arXiv

Testing the Cosmological Principle: Astrometric Limits on Systemic Motion of Quasars at Different Cosmological Epochs

A sample of $60,410$ bona fide optical quasars with astrometric proper motions in Gaia EDR3 and spectroscopic redshifts above 0.5 in an oval 8400 square degree area of the sky is constructed. Using orthogonal Zernike functions of polar coordinates, the proper motion fields are fitted in a weighted least-squares adjustment of the entire sample and of six equal bins of sorted redshifts. The overall fit with 37 Zernike functions reveals a statistically significant pattern, which is likely to be of instrumental origin. The main feature of this pattern is a chain of peaks and dips mostly in the R.A. component with an amplitude of 25~$μ$as yr$^{-1}$. This field is subtracted from each of the six analogous fits for quasars grouped by redshifts covering the range 0.5 through 7.03, with median values 0.72, 1.00, 1.25, 1.52, 1.83, 2.34. The resulting residual patterns are noisier, with formal uncertainties up to 8~$μ$as yr$^{-1}$ in the central part of the area. We detect a single high-confidence Zernike term for R.A. proper motion components of quasars with redshifts around 1.52 representing a general gradient of 30 $μ$as yr$^{-1}$ over $150\degr$ on the sky. We do not find any small- or medium-scale systemic variations of the residual proper motion field as functions of redshift above the $2.5\,σ$ significance level.

preprint2022arXiv

The tidal quality of the hot Jupiter WASP-12b

WASP-12b stands out among the planets of its class of hot Jupiters because of the observed fast orbital decay attributed to tidal dissipation. The measured rate of the orbital period is $\stackrel{\bf\centerdot}{\textstyle{P}}_{\rm orb}\,=\,-\,29\pm3\;\mbox{ms/yr}\;$=$\;(9.2\pm1.0)\times10^{-10}\;\mbox{s/s}$. In the literature heretofore, all attempts to explain this high rate were based on the assumption that the orbital evolution is dominated by the tides in the star. Since the modified tidal quality factor in yellow dwarfs is insufficient to warrant such a decay rate, a hypothesis was put forward that the star may actually be a subgiant. Using the latest data from the Gaia mission, we deduce that WASP-12a at $1.36\,{M_{Sun}}$ is an evolving dwarf at an early stage of post-turn-off evolution that has not yet depleted hydrogen in its core. Its unremarkable position in the color-magnitude diagram and the existence of close planets orbiting red giants of similar mass contradict the hypothesis of an abrupt boost of tidal quality due to structural internal changes. On the other hand, the previous research neglected the tidal dissipation in the planet assuming it to be negligible due to the likely synchronisation of its rotation and a presumed high quality factor. We critically reassess this assumption in the light of recent astrometric results for Jupiter and Saturn. Assuming that the structure of WASP-12b is similar to that of our Jupiter and Saturn, we find that the observed orbital decay is well explained by the tides in the planet. The estimated value of the planet's modified quality factor coincides almost precisely with that of our Jupiter.

preprint2022arXiv

USNO Bright Star Catalog, version 1

USNO Bright Star Catalog (UBSC) is a new astrometric catalog of 1423 brightest stars covering the entire sky, which is published online. It is nearly complete to $V=3$ mag except for three stellar systems. A combined astrometric solution of the Hipparcos Intermediate Astrometry Data and two dedicated ground-based campaigns in 2013 -- 2020 is the basis for this catalog. The astrometric parameters for each star include position coordinates, parallax, proper motion components, and their covariances on the Hipparcos mean epoch 1991.25. 64 percent of the catalog are flagged as known or suspected double or binary stars. UBSC lists 68 stars missing in Gaia EDR3 and another 114 stars without Gaia parallaxes or proper motions. The formal precision achieved for proper motions is comparable to that of Gaia.

preprint2021arXiv

Mass ratios of long-period binary stars resolved in precision astrometry catalogs of two epochs

Mass ratios of widely separated, long-period, resolved binary stars can be directly estimated from the available data in major space astrometry catalogs, such as the ESA's Hipparcos and Gaia mission results. The method is based on the universal principle of inertial motion of the system's center of mass in the absence of external forces, and is independent of any assumptions about the physical parameters or stellar models. The application is limited by the precision of input astrometric data, the orbital period and distance to the system, and possible presence of other attractors in the vicinity, such as in triple systems. A generalization of this technique to triples is proposed, as well as approaches to estimation of uncertainties. The known long-period binary HIP 473 AB is discussed as an application example, for which a $m_2/m_1=0.996^{+0.026}_{-0.026}$ is obtained.

preprint2020arXiv

A $5/4$ commensurability of KIC 5773205, the smallest eclipsing red dwarf detected by the Kepler mission

KIC 5773205 is the least luminous eclipsing M dwarf found in the Villanova catalog of eclipsing binaries detected by the {\it Kepler} mission. We processed and analyzed the three available quarters of mission data for this star and discovered a persistent periodic variation of the light curve with a period, which is in exact 4:5 commensurability to the orbital period. Three routes of interpretation are considered: 1) non-radial pulsations excited by the tidal interaction at a specific eigenfrequency; 2) a high-order spin-orbit resonance caused by the tides; 3) an ellipsoidal deformation caused by an outer orbiting companion in a mean motion resonance. All three explanations meet considerable difficulties, but the available facts seem to favor the tidally driven pulsation scenario. The star may represent a new type of heartbeat binary with tidally excited pulsations that are close to the orbital motion in frequency.

preprint2020arXiv

Constraining the origin of the planetary debris surrounding ZTF J0139+5245 through rotational fission of a triaxial asteroid

White dwarfs containing orbiting planetesimals or their debris represent crucial benchmarks by which theoretical investigations of post-main-sequence planetary systems may be calibrated. The photometric transit signatures of likely planetary debris in the ZTF J0139+5245 white dwarf system has an orbital period of about 110 days. An asteroid which breaks up to produce this debris may spin itself to destruction through repeated close encounters with the star without entering its Roche radius and without influence from the white dwarf's luminosity. Here, we place coupled constraints on the orbital pericentre ($q$) and the ratio ($β$) of the middle to longest semiaxes of a triaxial asteroid which disrupts outside of this white dwarf's Roche radius ($r_{\rm Roche}$) soon after attaining its 110-day orbit. We find that disruption within tens of years is likely when $β\lesssim 0.6$ and $q\approx 1.0-2.0r_{\rm Roche}$, and when $β\lesssim 0.2$ out to $q\approx 2.5r_{\rm Roche}$. Analysing the longer-timescale disruption of triaxial asteroids around ZTF J0139+5245 is desirable but may require either an analytical approach relying on ergodic theory or novel numerical techniques.

preprint2020arXiv

Gyr-timescale destruction of high-eccentricity asteroids by spin and why 2006 HY51 has been spared

Asteroids and other small celestial bodies have markedly prolate shapes, and the perturbative triaxial torques which are applied during pericenter passages in highly eccentric orbits trigger and sustain a state of chaotic rotation. Because the prograde spin rate around the principal axis of inertia is not bounded from above, it can accidentally reach the threshold value corresponding to rotational break-up. Previous investigations of this process were limited to integrations of $\sim 10^3$ orbits because of the stiff equation of motion. We present here a fast 1D simulation method to compute the evolution of this spin rate over $\sim 10^9$ orbits. We apply the method to the most eccentric solar system asteroid known, 2006 HY51 (with $e = 0.9684$), and find that for any reasonably expected shape parameters, it can never be accelerated to break-up speed. However, primordial solar system asteroids on more eccentric orbits may have already broken up from this type of rotational fission. The method also represents a promising opportunity to investigate the long-term evolution of extremely eccentric triaxial exo-asteroids ($e > 0.99$), which are thought to be common in white dwarf planetary systems

preprint2019arXiv

Chaotic rotation and evolution of asteroids and small planets in high-eccentricity orbits around white dwarfs

Observed planetary debris in white dwarf atmospheres predominately originate from the destruction of small bodies on highly eccentric ($>0.99$) orbits. Despite their importance, these minor planets have coupled physical and orbital evolution which has remained largely unexplored. Here, we present a novel approach for estimating the influence of fast chaotic rotation on the orbital evolution of high-eccentricity triaxial asteroids, and formally characterize the propagation of their angular rotation velocities and orbital elements as random time processes. By employing the impulse approximation, we demonstrate that the violent gravitational interactions during periastron passages transfer energy between the orbit and asteroid's rotation. If the distribution of spin impulses were symmetric around zero, then the net result would be a secular decrease of the semimajor axis and a further increase of the eccentricity. We find evidence, however, that the chaotic rotation may be self-regulated in such a manner that these effects are reduced or nullified. We discover that asteroids on highly eccentric orbits can break themselves apart --- in a type of YORP-less rotational fission --- without actually entering the Roche radius, with potentially significant consequences for the distribution of debris and energy requirements for gravitational scattering in metal-polluted white dwarf planetary systems. This mechanism provides a steady stream of material impacting a white dwarf without rapidly depleting the number of small bodies in the stellar system.

preprint2016arXiv

Variability induced motion in Kepler data

Variability induced motion (VIM) is an observable effect in simultaneous astrometric and photometric measurements caused by brightness variation in one of the components of a double source or blended image, which manifests itself as a strongly correlated shift of the optical photocenter. We have processed the entire collection of the Kepler long-cadence light curve data looking for correlated signals in astrometry and photometry on the time basis of a quarter year. Limiting the VIM correlation coefficient to 0.3, VIM events are detected for 129,525 Kepler stars at least in one quarter. Of 7305 Kepler objects of interest (KOI), 4440 are detected as VIM at least once. Known variable stars and resolved double stars have elevated rates of VIM detection. Confident VIM occurrences are found for stars with suggested superflare events, indicating possible signal contamination. We present a complete catalog of all quarterly VIM detections. This catalog should be checked for such astrophysically significant events as transits of exoplanets, new eclipsing stars, and superflares of solar-type stars.

preprint2015arXiv

Equilibrium rotation of semiliquid exoplanets and satellites

A wide range of exoplanet and exomoon models are characterized by a finite average rigidity and a viscosity much lower than the typical values for terrestrials. Such semiliquid bodies may or may not have rigid crusts with permanent figures. Unlike planets with solid mantles and Earth-like rheology, semiliquid bodies can be captured into stable pseudosynchronous spin resonance, where the average rate of rotation is higher than the synchronous 1:1 resonance. Two basic conditions are derived for capture of planets with a triaxial figure into pseudosynchronous rotation, one related to the characteristic tidal wave number (the product of the tidal frequency by the Maxwell time), and the other to the orbital eccentricity. If a semiliquid object does not satisfy either of the two conditions, it is captured into the synchronous resonance. For nearly axially symmetric bodies, only the first condition is in place, and the other is much relaxed, so they should predominantly be pseudosynchronous. It is also pointed out that the equilibrium pseudosychronous rotation rate can not reach the widely used asymptotic value from the constant time lag model but is in reality closer to the synchronous spin.

preprint2015arXiv

Tidal dissipation in a homogeneous spherical body. I. Methods

A formula for the tidal dissipation rate in a spherical body is derived from first principles, to correct some mathematical inaccuracies found in the literature. The development is combined with the Darwin-Kaula formalism for tides. Our intermediate results are compared with those by Zschau (1978) and Platzman (1984). When restricted to the special case of an incompressible spherical planet spinning synchronously without libration, our final formula can be compared with the commonly used expression from Peale & Cassen (1978, Eqn. 31). The two turn out to differ. In our expression, the contributions from all Fourier modes are positive-definite, this not being the case of the formula from Ibid. (The presence of negative terms in their formula was noticed by Makarov 2013.) Examples of application of our expression for the tidal damping rate are provided in the work by Makarov and Efroimsky (2014).

preprint2015arXiv

UrHip Proper Motion Catalog

Proper motions are computed and collected in a catalog using the Hipparcos positions (epoch 1991.25) and URAT1 positions (epoch 2012.3 to 2014.6). The goal is to obtain a significant improvement on the proper motion accuracy of single stars in the northern hemisphere, and to identify new astrometric binaries perturbed by orbital motion. For binaries and multiple systems, the longer baseline of Tycho2 (~ 100 yr) makes it more reliable despite its larger formal uncertainties. The resulting proper motions obtained for 67,340 stars have a consequent gain in accuracy by a factor of ~ 3 compared to Hipparcos. Comparison between UrHip and Hipparcos shows that they are reasonably close, but also reveals stars with large discrepant proper motions, a fraction of which are potential binary candidates.

preprint2014arXiv

Radial velocities and binarity of southern SIM grid stars

We present analysis of precision radial velocities (RV) of 1134 mostly red giant stars in the southern sky, selected as candidate astrometric grid objects for the Space Interferometry Mission (SIM). Only a few (typically, 2 or 3) spectroscopic observations per star have been collected, with the main goal of screening binary systems. The estimated rate of spectroscopic binarity in this sample of red giants is 32% at the 0.95 confidence level, and 46% at the 0.75 confidence. The true binarity rate is likely to be higher, because our method is not quite sensitive to very wide binaries and low-mass companions. The estimated lower and upper bounds of stellar RV jitter for the entire sample are 24 and 51 m/s, respectively; the adopted mean value is 37 m/s. A few objects of interest are identified with large variations of radial velocities, implying abnormally high mass ratios.

preprint2014arXiv

Tidal dissipation in a homogeneous spherical body. II. Three examples: Mercury, Io, and Kepler-10 b

In Efroimsky & Makarov (2014), we derived from the first principles a formula for the tidal heating rate in a tidally perturbed homogeneous sphere. We compared it with the formulae used in the literature, and pointed out the differences. Using this result, we now present three case studies - Mercury, Kepler-10b, and a triaxial Io. A very sharp frequency-dependence of k2/Q near spin-orbit resonances yields a similarly sharp dependence of k2/Q on the spin rate. This indicates that physical libration may play a major role in tidal heating of synchronously rotating bodies. The magnitude of libration in the spin rate being defined by the planet's triaxiality, the latter should be a factor determining the dissipation rate. Other parameters equal, a synchronously rotating body with a stronger triaxiality should generate more heat than a similar body of a more symmetrical shape. Further in the paper, we discuss scenarios where initially triaxial objects melt and lose their triaxiality. Thereafter, dissipation in them becomes less intensive; so the bodies freeze. The tidal bulge becomes a new permanent figure, with a new triaxiality lower than the original. In the paper, we also derive simplified, approximate expressions for dissipation rate in a rocky planet of the Maxwell rheology, with a not too small Maxwell time. The three expressions derived pertain to the cases of a synchronous spin, a 3:2 resonance, and a nonresonant rotation; so they can be applied to most close-in super-Earth exoplanets detected thus far. In such bodies, the rate of tidal heating outside of synchronous rotation is weakly dependent on the eccentricity and obliquity, provided both these parameters are small or moderate. According to our calculation, Kepler-10b could hardly survive the great amount of tidal heating without being synchronised, circularised and also reshaped through a complete or partial melt-down.

preprint2013arXiv

Dynamical evolution and spin-orbit resonances of potentially habitable exoplanets. The case of GJ 581d

GJ 581d is a potentially habitable super-Earth in the multiple system of exoplanets orbiting a nearby M dwarf. We investigate this planet's long-term dynamics, with an emphasis on its probable final rotation states acquired via tidal interaction with the host. The published radial velocities for the star are re-analysed with a benchmark planet detection algorithm, to confirm that there is no evidence for the recently proposed two additional planets (f and g). Limiting the scope to the four originally detected planets, we assess the dynamical stability of the system and find bounded chaos in the orbital motion. For the planet d, the characteristic Lyapunov time is 38 yr. Long-term numerical integration reveals that the system of four planets is stable, with the eccentricity of the planet d changing quasi-periodically in a tight range around 0.27, and with its semimajor axis varying only a little. The spin-orbit interaction of GJ 581d with its host star is dominated by the tides exerted by the star on the planet. We model this interaction, assuming a terrestrial composition of the mantle. Besides the customarily included secular parts of the triaxiality-caused and tidal torques, we also include these torques' oscillating components. It turns out that, dependent on the mantle temperature, the planet gets trapped into the 2:1 or an even higher spin-orbit resonance. It is very improbable that the planet could have reached the 1:1 resonance. This enhances the possibility of the planet being suitable for sustained life.

preprint2013arXiv

Dynamical evolution and spin-orbit resonances of potentially habitable exoplanets. The case of GJ 667C

We investigate the dynamical evolution of the potentially habitable super-earth GJ 667Cc in the multiple system of at least two exoplanets orbiting a nearby M dwarf, paying special attention to its spin-orbital state. The published radial velocities for this star are re-analyzed and evidence is found for additional periodic signals, which could be taken for two additional planets on eccentric orbits. Such systems are not dynamically viable and break up quickly in numerical integrations. Limiting the scope to the two originally detected planets, we assess the dynamical stability of the system and find no evidence for bounded chaos in the orbital motion. The orbital eccentricity of the planets b and c is found to change cyclicly in the range 0.06 - 0.28 and 0.05 - 0.25, respectively, with a period of approximately 0.46 yr. Taking the eccentricity variation into account, numerical integrations are performed of the differential equations modeling the spin-orbit interaction of the planet GJ 667Cc with its host star assuming a terrestrial composition of its mantle. Depending on the interior temperature of the planet, it is likely to be entrapped in the 3:2 (probability 0.51) or even higher spin-orbit resonance. It is less likely to reach the 1:1 resonance (probability 0.24). The estimated characteristic spin-down times are quite short for the two planets, i.e., within 1 Myr for planet c and even shorter for planet b. Both planet arrived at their current and, most likely, ultimate spin-orbit states a long time ago. The rate of tidal dissipation of energy is formidably high in the planets of GJ 667, estimated at $10^{23.7}$ and $10^{26.7}$ J yr$^{-1}$ for c and b, respectively. This raises a question of how such relatively massive, close super-Earths could survive overheating and destruction.

preprint2013arXiv

No pseudosynchronous rotation for terrestrial planets and moons

We reexamine the popular belief that a telluric planet or satellite on an eccentric orbit can, outside a spin-orbit resonance, be captured in a quasi-static tidal equilibrium called pseudosynchronous rotation. The existence of such configurations was deduced from oversimplified tidal models assuming either a constant tidal torque or a torque linear in the tidal frequency. A more accurate treatment requires that the torque be decomposed into the Darwin-Kaula series over the tidal modes, and that this decomposition be combined with a realistic choice of rheological properties of the mantle. This development demonstrates that there exist no stable equilibrium states for solid planets and moons, other than spin-orbit resonances.

preprint2013arXiv

Tidal Friction and Tidal Lagging. Applicability Limitations of a Popular Formula for the Tidal Torque

Tidal torques play a key role in rotational dynamics of celestial bodies. They govern these bodies' tidal despinning, and also participate in the subtle process of entrapment of these bodies into spin-orbit resonances. This makes tidal torques directly relevant to the studies of habitability of planets and their moons. Our work begins with an explanation of how friction and lagging should be built into the theory of bodily tides. Although much of this material can be found in various publications, a short but self-consistent summary on the topic has been lacking in the hitherto literature, and we are filling the gap. After these preparations, we address a popular concise formula for the tidal torque, which is often used in the literature, for planets or stars.We explain why the derivation of this expression, offered in the paper by Goldreich (1966; AJ 71, 1 - 7) and in the books by Kaula (1968, eqn. 4.5.29), and Murray & Dermott (1999, eqn. 4.159), implicitly sets the time lag to be frequency independent. Accordingly, the ensuing expression for the torque can be applied only to bodies having a very special (and very hypothetical) rheology which makes the time lag frequency independent, i.e, the same for all Fourier modes in the spectrum of tide. This expression for the torque should not be used for bodies of other rheologies. Specifically, the expression cannot be combined with an extra assertion of the geometric lag (or the phase lag) being constant, because at finite eccentricities the said assumption is incompatible with the constant-time-lag condition.

preprint2013arXiv

Why is the Moon synchronously rotating?

If the Moon's spin evolved from faster prograde rates, it could have been captured into a higher spin-orbit resonance than the current 1:1 resonance. At the current value of orbital eccentricity, the probability of capture into the 3:2 resonance is as high as 0.6, but it strongly depends on the temperature and average viscosity of the Moon's interior. A warmer, less viscous Moon on a higher-eccentricity orbit is even more easily captured into supersynchronous resonances. We discuss two likely scenarios for the present spin-orbit state: a cold Moon on a low-eccentricity orbit and a retrograde initial rotation.

preprint2012arXiv

How terrestrial planets traverse spin-orbit resonances: A camel goes through a needle's eye

The dynamical evolution of terrestrial planets resembling Mercury in the vicinity of spin-orbit resonances is investigated using comprehensive harmonic expansions of the tidal torque taking into account the frequency-dependent quality factors and Love numbers. The torque equations are integrated numerically with a small step in time, includng the oscillating triaxial torque components but neglecting the layered structure of the planet and assuming a zero obliquity. We find that a Mercury-like planet with its current value of orbital eccentricity (0.2056) is always captured in the 3:2 resonance. The probability of capture in the higher 2:1 resonance is approximately 0.23. These results are confirmed by a semi-analytical estimation of capture probabilities as functions of eccentricity for both prograde and retrograde evolution of spin rate. As follows from analysis of equilibrium torques, entrapment in the 3:2 resonance is inevitable at eccentricities between 0.2 and 0.41. Considering the phase space parameters at the times of periastron, the range of spin rates and phase angles, for which an immediate resonance passage is triggered, is very narrow, and yet, a planet like Mercury rarely fails to align itself into this state of unstable equilibrium before it traverses the 2:1 resonance.

preprint2012arXiv

Revealing companions to nearby stars with astrometric acceleration

A subset of 51 Hipparcos astrometric binaries among FG dwarfs within 67pc has been surveyed with the NICI adaptive optics system at Gemini-S, directly resolving for the first time 17 sub-arcsecond companions and 7 wider ones. Using these data together with published speckle interferometry of 57 stars, we compare the statistics of resolved astrometric companions with those of a simulated binary population. The fraction of resolved companions is slightly lower than expected from binary statistics. About 10% of astrometric companions could be "dark" (white dwarfs and close pairs of late M-dwarfs). To our surprise, several binaries are found with companions too wide to explain the acceleration. Re-analysis of selected intermediate astrometric data shows that some acceleration solutions in the original Hipparcos catalog are spurious.

preprint2011arXiv

Stability, chaos and entrapment of stars in very wide pairs

The relative motion of stars and other celestial objects in very wide pairs, separated by distances of the order of 1 pc, is strongly influenced by the tidal gravitational potential of the Galaxy. The Coriolis component of the horizontal tidal force in the rotating reference frame tends to disrupt such marginally bound pairs. However, even extremely wide pairs of bodies can be bound over intervals of time comparable to the Hubble time, under appropriate initial conditions. Here we show that for arbitrary chosen initial coordinates of a pair of stars, there exists a volume of the space of initial velocity components where the orbits remain bound in the planar tidal field for longer than 10 Gyr, even though the initial separation is well outside the Jacobi radius. The boundary of this phase space of stable orbits is fractal, and the motion at the boundary conditions is clearly chaotic. We found that the pairs may remain confined for several Gyr, and then suddenly disintegrate due to a particularly close rendezvous. By reversing such long-term stable orbits, we find that entrapment of unrelated stars into wide pairs is possible, but should be quite rare. Careful analysis of precision astrometry surveys revealed that extremely wide pairs of stars are present in significant numbers in the Galaxy. These results are expected to help discriminating the cases of genuine binarity and chance entrapment, and to make inroads in testing the limits of Newtonian gravitation.

preprint2011arXiv

The Worst Distortions of Astrometric Instruments and Orthonormal Models for Rectangular Fields of View

The non-orthogonality of algebraic polynomials of field coordinates traditionally used to model field-dependent corrections to astrometric measurements, gives rise to subtle adverse effects. In particular, certain field dependent perturbations in the observational data propagate into the adjusted coefficients with considerable magnification. We explain how the worst perturbation, resulting in the largest solution error, can be computed for a given non-orthogonal distortion model. An algebraic distortion model of full rank can be converted into a fully orthonormal model based on the Zernike polynomials for a circular field of view, or a basis of functions constructed from the original model by a variant of the Gram-Schmidt orthogonalization process for a rectangular field of view. The relative significance of orthonormal distortion terms is assessed simply by the numerical values of the corresponding coefficients. Orthonormal distortion models are easily extendable when the distribution of residuals indicate the presence of higher order terms.

preprint2010arXiv

Variability of surface flows on the Sun and the implications for exoplanet detection

The published Mount Wilson Doppler-shift measurements of the solar velocity field taken in 1967--1982 are revisited with a more accurate model, which includes two terms representing the meridional flow and three terms corresponding to the convective limb shift. Integration of the recomputed data over the visible hemisphere reveals significant variability of the net radial velocity at characteristic time scales of 0.1--10 years, with a standard deviation of 1.4 \ms. This result is supported by independent published observations. The implications for exoplanet detection include reduced sensitivity of the Doppler method to Earth-like planets in the habitable zone, and an elevated probability of false detections at periods of a few to several years.