We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 from June 2000 to June 2019 and derived physical and dynamical properties of the binary systems. We found that the data for 1999 KW4 are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of -0.65 +/- 0.16 deg/yr2 (all quoted uncertainties are 3sigma). This means that the semimajor axis of the mutual orbit of the components of this binary system increases in time with a mean rate of 1.2 +/- 0.3 cm/yr. The data for 2001 SL9 are also inconsistent with a constant orbital period and we obtained two solutions for the quadratic drift of the mean anomaly: 2.8 +/- 0.2 and 5.2 +/- 0.2 deg/yr2, implying that the semimajor axis of the mutual orbit of the components decreases in time with a mean rate of -2.8 +/- 0.2 or -5.1 +/- 0.2 cm/yr for the two solutions, respectively. The expanding orbit of 1999 KW4 may be explained by mutual tides interplaying with binary YORP (BYORP) effect (McMahon and Scheeres, 2010). However, a modeling of the BYORP drift using radar-derived shapes of the binary components predicted a much higher value of the orbital drift than the observed one. It suggests that either the radar-derived shape model of the secondary is inadequate for computing the BYORP effect, or the present theory of BYORP overestimates it. It is possible that the BYORP coefficient has instead an opposite sign than predicted; in that case, the system may be moving into an equilibrium between the BYORP and the tides. In the case of 2001 SL9, the BYORP effect is the only known physical mechanism that can cause the inward drift of its mutual orbit.

We present the first-ever measurements of the surface colors and rotational properties of a recently discovered distant object in an unusual orbit 2016 ND21. In 2017, this object was close to its perihelion, which allowed us to observe it with a 2.0-m telescope at the Peak Terskol observatory. Broadband photometric observations of 2016 ND21 were carried out in October and December of 2017 in standard BVR filters of the Johnson-Cousins photometric system. During our observations, we did not detect any cometary activity. We found the rotational period of 17.53{\pm}0.02 hr, while another a little longer value of 17.65{\pm}0.02 hr is also possible. Assuming an equatorial aspect of observations, a peak-to-peak amplitude of A = 0.31{\pm}0.05 mag (or even higher since only one maximum and one minimum were well-measured) corresponds to an elongated body with an axis ratio a/b ~1.3. The lightcurve behavior indicates a complex, possibly non-convex, shape of this object. The visible absolute magnitude is HV =12.4{\pm}0.1 mag, which was estimated by using the linear phase slope 0.04 mag/deg as the most probable value from our observations. Assuming a surface albedo in the range of 0.04-0.10, the size of 2016 ND21 should be about 15-23 km. From our multi-color observations, we determined surface colors V-R = 0.69{\pm}0.04 mag, B-R = 1.79{\pm}0.08 mag, and B-V = 1.10{\pm}0.08 mag. The measured colors indicate an extremely red surface of this object. A very red surface is unusual for comets, which is in agreement with the fact, that no cometary activity was detected for 2016 ND21. The B-R color is higher than the typical B-R colors of the red D-type asteroids, but it is consistent with colors of the red Centaurs and TNOs classified as RR type in TNOs classification. This result gives a first evidence of a possible outer belt origin of this small body.

We studied the rotational properties of the dwarf planet Makemake. The photometric observations were carried out at different telescopes between 2006 and 2017. Most of the measurements were acquired in BVRI broad-band filters of a standard Johnson-Cousins photometric system. We found that Makemake rotates more slowly than was previously reported. A possible lightcurve asymmetry suggests a double-peaked period of P = 22.8266$\pm$0.0001~h. A small peak-to-peak lightcurve amplitude in R-filter A = 0.032$\pm$0.005 mag implies an almost spherical shape or near pole-on orientation. We also measured BVRI colours and the R-filter phase-angle slope and revised the absolute magnitudes. The absolute magnitude of Makemake has remained unchanged since its discovery in 2005. No direct evidence of a newly discovered satellite was found in our photometric data; however, we discuss the possible existence of another larger satellite.

We studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process.

Asteroid (3200) Phaethon is a Near-Earth Apollo asteroid with an unusual orbit that brings it closer to the Sun than any other known asteroid. Its last close approach to the Earth was in mid-December 2017 and the next one will be on October 2026. Previous rotationally time-resolved spectroscopy of Phaethon showed that its spectral slope is slightly bluish, in agreement with its B/F taxonomic classification, but at some rotational phases, it changes to slightly reddish. Motivated by this result we performed time-resolved imaging polarimetry of Phaethon during its recent close approach to the Earth. Phaethon has a spin period of 3.604 hours and we found a variation of the linear polarisation with rotation. This seems to be a rare case in which such variation is unambiguously found, also a consequence of its fairly large amplitude. Combining this new information with the brightness and colour variation, as well as previously reported results from Arecibo radar observations, we conclude that there is no variation of the mineralogy across the surface of Phaeton. However, the observed change in the linear polarisation may be related to differences in the thickness of the surface regolith in different areas or local topographic features.

A multi-colour phase-polarization curve of asteroid (3200)~Phaethon has been obtained during the December 2017 apparition by merging measurements taken at the observing station of Calern (France) and at the Rhozen observatory (Bulgaria). All the observations were obtained in the positive polarization branch, the phase angle ranging from 36$^\circ$ to 116$^\circ$. The measured values of linear polarization are among the highest ever observed for a Solar system body. The covered interval of phase angle was not sufficiently extended to derive a firm determination of the $P_{\rm max}$ parameter, but this appears to occur at a phase angle around 130$^\circ$ and reaches more than 45\% of linear polarization. Phaethon is the parent body of the Geminid meteor shower, and the real physical nature of this object (asteroid or comet) has been a long-debated subject. Our polarimetric measurements seem to support the asteroid hypothesis with a phase-polarization curve similar to the asteroid (2)~Pallas, but further observations at smaller phase angles are needed to draw definitive conclusions.

The rotation states of small asteroids are affected by a net torque arising from an anisotropic sunlight reflection and thermal radiation from the asteroids' surfaces. On long timescales, this so-called YORP effect can change asteroid spin directions and their rotation periods. We analyzed lightcurves of four selected near-Earth asteroids with the aim of detecting secular changes in their rotation rates that are caused by YORP. We use the lightcurve inversion method to model the observed lightcurves and include the change in the rotation rate $\mathrm{d} \omega / \mathrm{d} t$ as a free parameter of optimization. We collected more than 70 new lightcurves. For asteroids Toro and Cacus, we used thermal infrared data from the WISE spacecraft and estimated their size and thermal inertia. We also used the currently available optical and radar astrometry of Toro, Ra-Shalom, and Cacus to infer the Yarkovsky effect. We detected a YORP acceleration of $\mathrm{d}\omega / \mathrm{d} t = (1.9 \pm 0.3) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ for asteroid Cacus. For Toro, we have a tentative ($2\sigma$) detection of YORP from a significant improvement of the lightcurve fit for a nonzero value of $\mathrm{d}\omega / \mathrm{d} t = 3.0 \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$. For asteroid Eger, we confirmed the previously published YORP detection with more data and updated the YORP value to $(1.1 \pm 0.5) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. We also updated the shape model of asteroid Ra-Shalom and put an upper limit for the change of the rotation rate to $|\mathrm{d}\omega / \mathrm{d} t| \lesssim 1.5 \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. Ra-Shalom has a greater than $3\sigma$ Yarkovsky detection with a theoretical value consistent with observations assuming its size and/or density is slightly larger than the nominally expected values.

We propose an analytic model for the normal YORP and diurnal Yarkovsky effects experienced by a convex asteroid. Both the YORP torque and the Yarkovsky force are expressed as integrals of a universal function over the surface of an asteroid. Although in general this function can only be calculated numerically from the solution of the heat conductivity equation, approximate solutions can be obtained in quadratures for important limiting cases. We consider 3 such simplified models: Rubincam's approximation (zero heat conductivity), low thermal inertia limit (including the next order correction and thus valid for small heat conductivity), and high thermal inertia limit (valid for large heat conductivity). All three simplified models are compared with the exact solution.

We present optical observations of SN 2013dx, related to the Fermi burst GRB 130702A occurred at a redshift z = 0.145. It is the second-best sampled GRB-SN after SN~1998bw: the observational light curves contain more than 280 data points in uBgrRiz filters until 88 day after the burst, and the data were collected from our observational collaboration (Maidanak Observatory, Abastumani Observatory, Crimean Astrophysical Observatory, Mondy Observatory, National Observatory of Turkey, Observatorio del Roque de los Muchachos) and from the literature. We model numerically the multicolour light curves using the one-dimensional radiation hydrodynamical code STELLA, previously widely implemented for the modelling of typical non-GRB SNe. The best-fitted model has the following parameters: pre-supernova star mass M = 25 M_Sun, mass of a compact remnant M_CR = 6 M_Sun, total energy of the outburst E_oburst = 3.5 x 10^(52) erg, pre-supernova star radius R = 100 R_Sun, M_56Ni = 0.2 M_Sun which is totally mixed through the ejecta; M_O = 16.6 M_Sun, M_Si = 1.2 M_Sun, and M_Fe = 1.2 M_Sun, and the radiative efficiency of the SN is 0.1 per cent.

Context. Tangential YORP is a thermophysical effect that can alter the rotation rate of asteroids and is distinct from the "normal" YORP effect, but to date has only been studied for asteroids with zero obliquity. Aims. The tangential YORP force produced by spherical boulders on the surface of an asteroid with an arbitrary obliquity is studied. Methods. A finite element method is used to simulate heat conductivity inside a boulder, to find the recoil force experienced by it. Then an ellipsoidal asteroid uniformly covered by such boulders is considered and the torque is numerically integrated over its surface. Results. Tangential YORP is found to operate on non-zero obliquities and decreases by a factor of 2 for increasing obliquity.

Asteroid modeling efforts in the last decade resulted in a comprehensive dataset of almost 400 convex shape models and their rotation states. This amount already provided a deep insight into physical properties of main-belt asteroids or large collisional families. We aim to increase the number of asteroid shape models and rotation states. Such results are an important input for various further studies such as analysis of asteroid physical properties in different populations, including smaller collisional families, thermophysical modeling, and scaling shape models by disk-resolved images, or stellar occultation data. This provides, in combination with known masses, bulk density estimates, but constrains also theoretical collisional and evolutional models of the Solar System. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as an input for the convex inversion method, and derive 3D shape models of asteroids, together with their rotation periods and orientations of rotation axes. The key ingredient is the support of more that one hundred observers who submit their optical data to publicly available databases. We present updated shape models for 36 asteroids, for which mass estimates are currently available in the literature or their masses will be most likely determined from their gravitational influence on smaller bodies, which orbital deflection will be observed by the ESA Gaia astrometric mission. This was achieved by using additional optical data from recent apparitions for the shape optimization. Moreover, we also present new shape model determinations for 250 asteroids, including 13 Hungarias and 3 near-Earth asteroids.

Astrometric observations performed by the Gaia Follow-Up Network for Solar System Objects (Gaia-FUN-SSO) play a key role in ensuring that moving objects first detected by ESA's Gaia mission remain recoverable after their discovery. An observation campaign on the potentially hazardous asteroid (99 942) Apophis was conducted during the asteroid's latest period of visibility, from 12/21/2012 to 5/2/2013, to test the coordination and evaluate the overall performance of the Gaia-FUN-SSO . The 2732 high quality astrometric observations acquired during the Gaia-FUN-SSO campaign were reduced with the Platform for Reduction of Astronomical Images Automatically (PRAIA), using the USNO CCD Astrograph Catalogue 4 (UCAC4) as a reference. The astrometric reduction process and the precision of the newly obtained measurements are discussed. We compare the residuals of astrometric observations that we obtained using this reduction process to data sets that were individually reduced by observers and accepted by the Minor Planet Center. We obtained 2103 previously unpublished astrometric positions and provide these to the scientific community. Using these data we show that our reduction of this astrometric campaign with a reliable stellar catalog substantially improves the quality of the astrometric results. We present evidence that the new data will help to reduce the orbit uncertainty of Apophis during its close approach in 2029. We show that uncertainties due to geolocations of observing stations, as well as rounding of astrometric data can introduce an unnecessary degradation in the quality of the resulting astrometric positions. Finally, we discuss the impact of our campaign reduction on the recovery process of newly discovered asteroids.

Using our photometric observations taken between 1996 and 2013 and other published data, we derived properties of the binary near-Earth asteroid (175706) 1996 FG3 including new measurements constraining evolution of the mutual orbit with potential consequences for the entire binary asteroid population. We also refined previously determined values of parameters of both components, making 1996 FG3 one of the most well understood binary asteroid systems. We determined the orbital vector with a substantially greater accuracy than before and we also placed constraints on a stability of the orbit. Specifically, the ecliptic longitude and latitude of the orbital pole are 266{\deg} and -83{\deg}, respectively, with the mean radius of the uncertainty area of 4{\deg}, and the orbital period is 16.1508 +/- 0.0002 h (all quoted uncertainties correspond to 3sigma). We looked for a quadratic drift of the mean anomaly of the satellite and obtained a value of 0.04 +/- 0.20 deg/yr^2, i.e., consistent with zero. The drift is substantially lower than predicted by the pure binary YORP (BYORP) theory of McMahon and Scheeres (McMahon, J., Scheeres, D. [2010]. Icarus 209, 494-509) and it is consistent with the theory of an equilibrium between BYORP and tidal torques for synchronous binary asteroids as proposed by Jacobson and Scheeres (Jacobson, S.A., Scheeres, D. [2011]. ApJ Letters, 736, L19). Based on the assumption of equilibrium, we derived a ratio of the quality factor and tidal Love number of Q/k = 2.4 x 10^5 uncertain by a factor of five. We also derived a product of the rigidity and quality factor of mu Q = 1.3 x 10^7 Pa using the theory that assumes an elastic response of the asteroid material to the tidal forces. This very low value indicates that the primary of 1996 FG3 is a 'rubble pile', and it also calls for a re-thinking of the tidal energy dissipation in close asteroid binary systems.

This article discusses how re-emission of absorbed solar light by centimeter- to decimeter-sized structures on the surface of an asteroid can create a component of the recoil force parallel to the surface. Under certain conditions the west side of stones appear to be on average slightly warmer than their east sides, thus experiencing a stronger recoil force and increasing the rotation rate of the asteroid. We study this effect, called the Tangential YORP effect, in a toy model, replacing stones with walls and simulating heat conductivity in them. We discuss general trends of the effect, estimate its magnitude, and find it to be comparable to the normal YORP effect determined by gross-scale asymmetry of the asteroid. The existence of this effect would modify the predictions of the YORP acceleration of asteroids.

Tangential YORP, or TYORP, has recently been demonstrated to be an important factor in the evolution of an asteroid's rotation state. It is complementary to normal YORP, or NYORP, which used to be considered previously. While NYORP is produced by non-symmetry in the large-scale geometry of an asteroid, TYORP is due to heat conductivity in stones on the surface of the asteroid. Yet to date TYORP has been studied only in a simplified 1-dimensional model, substituting stones by high long walls. This article for the first time considers TYORP in a realistic 3-dimensional model, also including shadowing and self-illumination effects via ray tracing. TYORP is simulated for spherical stones lying on regolith. The model includes only 5 free parameters, and the dependence of the TYORP on each of them is studied. The TYORP torque appears to be smaller than previous estimates from 1-dimensional model, but still comparable to the NYORP torques. These results can be used to estimate TYORP of different asteroids, and also as a basis for more sophisticated models of TYORP.

The spin state of small asteroids can change on a long timescale by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric lightcurves. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the aim to detect possible deviations from the constant rotation caused by the YORP effect. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as a linear change in the rotation rate in time d\omega /dt. Values of d\omega/ dt derived from observations were compared with the values predicted by theory. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger d\omega / dt = (1.4 +/- 0.6) x 10^{-8} rad/d (3\sigma error), which corresponds to a decrease in the rotation period by 4.2 ms/yr. The photometry of Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could only constrain maximum values of |d\omega / dt| < 8 x 10^{-9} rad/d for Cerberus, and |d\omega / dt| < 3 x 10^{-8} rad/d for Ra-Shalom.

Main belt asteroids (6070) Rheinland and (54827) 2001NQ8 belong to a small population of couples of bodies which reside on very similar heliocentric orbits. Vokrouhlicky & Nesvorny (2008, AJ 136, 280) promoted a term "asteroid pairs", pointing out their common origin within the past tens to hundreds of ky. Previous attempts to reconstruct the initial configuration of Rheinland and 2001NQ8 at the time of their separation have led to the prediction that Rheinland's rotation should be retrograde. Here we report extensive photometric observations of this asteroid and use the lightcurve inversion technique to directly determine its rotation state and shape. We confirm the retrograde sense of rotation of Rheinland, with obliquity value constrained to be >= 140 deg. The ecliptic longitude of the pole position is not well constrained as yet. The asymmetric behavior of Rheinland's lightcurve reflects a sharp, near-planar edge in our convex shape representation of this asteroid. Our calibrated observations in the red filter also allow us to determine $H_R = 13.68\pm 0.05$ and $G = 0.31\pm 0.05$ values of the H-G system. With the characteristic color index $V-R = 0.49\pm 0.05$ for the S-type asteroids, we thus obtain $H = 14.17\pm 0.07$ for the absolute magnitude of (6070) Rheinland. This a significantly larger value than previously obtained from analysis of the astrometric survey observations. We next use the obliquity constraint for Rheinland to eliminate some degree of uncertainty in the past propagation of its orbit. This is because the sign of the past secular change of its semimajor axis due to the Yarkovsky effect is now constrained. Determination of the rotation state of the secondary component, asteroid (54827) 2001NQ8, is the key element in further constraining the age of the pair and its formation process.

Asteroid pairs sharing similar heliocentric orbits were found recently. Backward integrations of their orbits indicated that they separated gently with low relative velocities, but did not provide additional insight into their formation mechanism. A previously hypothesized rotational fission process4 may explain their formation - critical predictions are that the mass ratios are less than about 0.2 and, as the mass ratio approaches this upper limit, the spin period of the larger body becomes long. Here we report photometric observations of a sample of asteroid pairs revealing that primaries of pairs with mass ratios much less than 0.2 rotate rapidly, near their critical fission frequency. As the mass ratio approaches 0.2, the primary period grows long. This occurs as the total energy of the system approaches zero requiring the asteroid pair to extract an increasing fraction of energy from the primary's spin in order to escape. We do not find asteroid pairs with mass ratios larger than 0.2. Rotationally fissioned systems beyond this limit have insufficient energy to disrupt. We conclude that asteroid pairs are formed by the rotational fission of a parent asteroid into a proto-binary system which subsequently disrupts under its own internal system dynamics soon after formation.

A wide observational campaign was carried out in 2004-2009 aimed to complete the ground-based investigation of Lutetia prior to the Rosetta fly-by in July 2010. We have obtained BVRI photometric and V-band polarimetric measurements over a wide range of phase angles, and visible and infrared spectra in the 0.4-2.4 micron range. We analyzed them together with previously published data to retrieve information on Lutetia's surface properties. Values of lightcurve amplitudes, absolute magnitude, opposition effect, phase coefficient and BVRI colors of Lutetia surface seen at near pole-on aspect have been determined. We defined more precisely parameters of polarization phase curve and showed their distinct deviation from any other moderate-albedo asteroid. An indication of possible variations both in polarization and spectral data across the asteroid surface was found. To explain features found by different techniques we propose that (i) Lutetia has a non-convex shape, probably due to the presence of a large crater, and heterogeneous surface properties probably related to surface morphology; (ii) at least part of the surface is covered by a fine-grained regolith with particle size less than 20 microns; (iii) the closest meteorite analogues of Lutetia's surface composition are particular types of carbonaceous chondrites or Lutetia has specific surface composition not representative among studied meteorites.

We present observations of the afterglow of GRB 080319B at optical, mm and radio frequencies from a few hours to 67 days after the burst. Present observations along with other published multi-wavelength data have been used to study the light-curves and spectral energy distributions of the burst afterglow. The nature of this brightest cosmic explosion has been explored based on the observed properties and it's comparison with the afterglow models. Our results show that the observed features of the afterglow fits equally good with the Inter Stellar Matter and the Stellar Wind density profiles of the circum-burst medium. In case of both density profiles, location of the maximum synchrotron frequency $\nu_m$ is below optical and the value of cooling break frequency $\nu_c$ is below $X-$rays, $\sim 10^{4}$s after the burst. Also, the derived value of the Lorentz factor at the time of naked eye brightness is $\sim 300$ with the corresponding blast wave size of $\sim 10^{18}$ cm. The numerical fit to the multi-wavelength afterglow data constraints the values of physical parameters and the emission mechanism of the burst.

The analysis of the density and brightness of big planets satellites, main asteroid belt objects, Kuiper belt objects and centaurs has been carried out as well as the analysis of suspected unseen satellites of the stars. According to the date on the first of January 2006 the catalogue of planetary objects has been compiled.