Aims: Our goal was to increase number of known doubly eclipsing systems such that the resulting dataset would allow to study them via statistical means, as well as prove that they constitute gravitationally bound 2+2 quadruple system. Methods: We analysed photometric data for eclipsing binaries provided by the OGLE survey in the LMC fields. We found a large number of new doubly eclipsing systems (our discoveries are 3x more numerous than previous studies). With a typical orbital period of days for the binaries, we sought eclipse time variations (ETVs) on the timescale of years. In the cases where we were able to detect the ETV period, the difference between the inner and outer periods in the quadruple system is large enough. This allows us to interpret ETVs primarily as the light-time effect, thus providing an interesting constraint on masses of the binaries. Results: In addition to significantly enlarging the database of known doubly eclipsing systems, we performed a thorough analysis of 72 cases. ETVs for 28 of them (39% of the studied cases) showed evidence of relative motion. We note OGLE BLG-ECL-145467 as the most interesting case; it is bright (I=12.6 mag), consists of two detached binaries with periods of about 3.3 d and 4.9 d (making it a candidate for a 3:2 resonant system), mutual period about 1538 d. Distribution of the orbital period ratio P_A/P_B of binaries in 2+2 quadruples shows statistically significant excess at 1 and 1.5. The former is likely a natural statistical preference in weakly interacting systems with periods within the same range. The latter is thought to be evidence of a capture in the 3:2 mean motion resonance of the two binaries. This sets important constraints on evolutionary channels in these systems. The total number of doubly eclipsing systems increased to 146, more than 90% of which are at low declinations on the southern sky.

The recent close approach of the NEA (3200) Phaethon offered a rare opportunity to obtain high-quality observational data. We used the newly obtained optical light curves to improve the spin and shape model of Phaethon and to determine its surface physical properties derived by thermophysical modeling. We also used the available astrometric observations of Phaethon, including those obtained by the Arecibo radar and the Gaia spacecraft, to constrain the secular drift of the orbital semimajor axis. This constraint allowed us to estimate the bulk density by assuming that the drift is dominated by the Yarkovsky effect. We used the convex inversion model to derive the 3D shape model of Phaethon, and a detailed numerical approach for an accurate analysis of the Yarkovsky effect. We obtained a unique solution for Phaethon's pole orientation at $(318,-47)^{\circ}$ ecliptic longitude and latitude (uncertainty of $5^{\circ}$), and confirm the previously reported thermophysical properties ($D=5.1\pm0.2$ km, $\Gamma=600\pm200$ SI). Phaethon has a top-like shape with possible north-south asymmetry. The characteristic size of the regolith grains is 1-2 cm. The orbit analysis reveals a secular drift of the semimajor axis of $-(6.9\pm1.9)\times 10^{-4}$ au Myr$^{-1}$. With the derived volume-equivalent size of 5.1~km, the bulk density $\rho$ is $1.67\pm0.47$ g cm$^{-3}$. If the size is slightly larger $\sim5.7$ km, as suggested by radar data, $\rho$ would decrease to $1.48\pm0.42$ g cm$^{-3}$. We further investigated the suggestion that Phaethon may be in a cluster with asteroids (155140) 2005 UD and (225416) 1999 YC that was formed by rotational fission of a critically spinning parent body. Phaethon's $\rho$ is consistent with typical values for large ($>100$ km) C-complex asteroids and supports its association with asteroid (2) Pallas. These findings render a cometary origin unlikely for Phaethon.

We present an analysis of the apsidal motion and light curve parameters of 54 never-before-studied galactic Algol-type binaries. This is the first analysis of such a large sample of eccentric eclipsing binaries in our Galaxy, and has enabled us to identify several systems that are worthy of further study. Bringing together data from various databases and surveys, supplemented with new observations, we have been able to trace the long-term evolution of the eccentric orbit over durations extending back up to several decades. Our present study explores a rather different sample of stars to those presented in the previously published catalogue of eccentric eclipsing binaries, sampling to fainter magnitudes, covering later spectral types, sensitive to different orbital periods with more than 50% of our systems having periods longer than six days. The typical apsidal motion in the sample is rather slow (mostly of order of centuries long), although in some cases this is less than 50 yr. All of the systems, except one, have eccentricities less than 0.5, with an average value of 0.23. Several of the stars also show evidence for additional period variability. In particular we can identify three systems in the sample, HD 44093, V611 Pup, and HD 313631, which likely represent relativistic apsidal rotators.

Since its discovery in 2008, the Andromeda galaxy nova M31N 2008-12a has been observed in eruption every single year. This unprecedented frequency indicates an extreme object, with a massive white dwarf and a high accretion rate, which is the most promising candidate for the single-degenerate progenitor of a type-Ia supernova known to date. The previous three eruptions of M31N 2008-12a have displayed remarkably homogeneous multi-wavelength properties: (i) From a faint peak, the optical light curve declined rapidly by two magnitudes in less than two days; (ii) Early spectra showed initial high velocities that slowed down significantly within days and displayed clear He/N lines throughout; (iii) The supersoft X-ray source (SSS) phase of the nova began extremely early, six days after eruption, and only lasted for about two weeks. In contrast, the peculiar 2016 eruption was clearly different. Here we report (i) the considerable delay in the 2016 eruption date, (ii) the significantly shorter SSS phase, and (iii) the brighter optical peak magnitude (with a hitherto unobserved cusp shape). Early theoretical models suggest that these three different effects can be consistently understood as caused by a lower quiescence mass-accretion rate. The corresponding higher ignition mass caused a brighter peak in the free-free emission model. The less-massive accretion disk experienced greater disruption, consequently delaying re-establishment of effective accretion. Without the early refueling, the SSS phase was shortened. Observing the next few eruptions will determine whether the properties of the 2016 outburst make it a genuine outlier in the evolution of M31N 2008-12a.

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.

New CCD photometric observations of fourteen short-period low-mass eclipsing binaries (LMB) in the photometric filters I, R and V were used for the light curve analysis. There still exists a discrepancy between radii as observed and those derived from the theoretical modelling for LMB in general. Mass calibration of all observed LMB was done using only the photometric indices. The light curve modelling of these selected systems were performed, yielding the new derived masses and radii for both components. We compared these systems with the compilation of other known double-lined LMB systems with uncertainties of masses and radii less then 5 \%, which includes 66 components of binaries where both spectroscopy and photometry were combined together. All of our systems are circular short-period binaries, and for some of them the photospheric spots were also used. A purely photometric study of the light curves without spectroscopy seems unable to achieve high enough precision and accuracy in the masses and radii to provide for a meaningful test of the M-R relation for low-mass stars.

The Andromeda Galaxy recurrent nova M31N 2008-12a had been observed in eruption ten times, including yearly eruptions from 2008-2014. With a measured recurrence period of $P_\mathrm{rec}=351\pm13$ days (we believe the true value to be half of this) and a white dwarf very close to the Chandrasekhar limit, M31N 2008-12a has become the leading pre-explosion supernova type Ia progenitor candidate. Following multi-wavelength follow-up observations of the 2013 and 2014 eruptions, we initiated a campaign to ensure early detection of the predicted 2015 eruption, which triggered ambitious ground and space-based follow-up programs. In this paper we present the 2015 detection; visible to near-infrared photometry and visible spectroscopy; and ultraviolet and X-ray observations from the Swift observatory. The LCOGT 2m (Hawaii) discovered the 2015 eruption, estimated to have commenced at Aug. $28.28\pm0.12$ UT. The 2013-2015 eruptions are remarkably similar at all wavelengths. New early spectroscopic observations reveal short-lived emission from material with velocities $\sim13000$ km s$^{-1}$, possibly collimated outflows. Photometric and spectroscopic observations of the eruption provide strong evidence supporting a red giant donor. An apparently stochastic variability during the early super-soft X-ray phase was comparable in amplitude and duration to past eruptions, but the 2013 and 2015 eruptions show evidence of a brief flux dip during this phase. The multi-eruption Swift/XRT spectra show tentative evidence of high-ionization emission lines above a high-temperature continuum. Following Henze et al. (2015a), the updated recurrence period based on all known eruptions is $P_\mathrm{rec}=174\pm10$ d, and we expect the next eruption of M31N 2008-12a to occur around mid-Sep. 2016.

Analyzing the available photometry from the Kepler satellite and other databases, we performed detailed light curve modeling of 10 eclipsing binary systems that were found to exhibit a periodic modulation of their orbital periods. All of the selected systems are detached Algol type, with orbital periods from 0.9 to 2.9 days. In total, 9448 times of minimum for these binaries were analyzed in an attempt to identify the period variations caused by the third bodies in these systems. The well-known method of the light-travel time effect was used for the analysis. The orbital periods of the outer bodies were found to be between 1 and 14 years. This hypothesis makes such systems interesting for future prospective detections of these components, despite their low predicted masses. Considering the dynamical interaction between the orbits, the system KIC 3440230 seems to be the most interesting, in which one would expect the detection of some effects (i.e., changing the inclination) even after a few years or decades of observations.

The minimum timings of eclipsing binaries V418 Aql, SU Boo, RV CVn, CR Cas, GV Cyg, V432 Per, and BD+42 2782 were collected and analyzed. Their long-term behavior was studied via period analysis, revealing a periodic term in eclipse times. We derived 576 new times of minimum. Hence, to describe the periodic variation, a thirdbody hypothesis was proposed and the resulting orbital periods are as follows: 70, 7.4, 53, 37, 27, 53, and 18 yr, respectively. For the system V432 Per an additional 9.5 yr variation was also found. The predicted minimum masses of these distant bodies were calculated and their detectability discussed. The light curves of SU Boo and RV CVn were analyzed using the PHOEBE program, resulting in physical parameters of the components. New variable stars in the field of V418 Aql were discovered.

The photometric and spectroscopic observational campaign organized for the 2008/9 eclipse of EE Cep revealed features, which indicate that the eclipsing disk in the EE Cep system has a multi-ring structure. We suggest that the gaps in the disk can be related to the possible planet formation.

We systematically surveyed period variations of superhumps in SU UMa-type dwarf novae based on newly obtained data and past publications. In many systems, the evolution of superhump period are found to be composed of three distinct stages: early evolutionary stage with a longer superhump period, middle stage with systematically varying periods, final stage with a shorter, stable superhump period. During the middle stage, many systems with superhump periods less than 0.08 d show positive period derivatives. Contrary to the earlier claim, we found no clear evidence for variation of period derivatives between superoutburst of the same object. We present an interpretation that the lengthening of the superhump period is a result of outward propagation of the eccentricity wave and is limited by the radius near the tidal truncation. We interpret that late stage superhumps are rejuvenized excitation of 3:1 resonance when the superhumps in the outer disk is effectively quenched. Many of WZ Sge-type dwarf novae showed long-enduring superhumps during the post-superoutburst stage having periods longer than those during the main superoutburst. The period derivatives in WZ Sge-type dwarf novae are found to be strongly correlated with the fractional superhump excess, or consequently, mass ratio. WZ Sge-type dwarf novae with a long-lasting rebrightening or with multiple rebrightenings tend to have smaller period derivatives and are excellent candidate for the systems around or after the period minimum of evolution of cataclysmic variables (abridged).