The intermediate polar FO Aquarii (FO Aqr) experienced its first-reported low-accretion states in 2016, 2017, and 2018, and using newly available photographic plates, we identify pre-discovery low states in 1965, 1966, and 1974. The primary focus of our analysis, however, is an extensive set of time-series photometry obtained between 2002 and 2018, with particularly intensive coverage of the 2016-2018 low states. After computing an updated spin ephemeris for the white dwarf (WD), we show that its spin period began to increase in 2014 after having spent 27 years decreasing; no other intermediate polar has experienced a sign change of its period derivative, but FO Aqr has now done so twice. Our central finding is that the recent low states all occurred shortly after the WD began to spin down, even though no low states were reported in the preceding quarter-century, when it was spinning up. Additionally, the system's mode of accretion is extremely sensitive to the mass-transfer rate, with accretion being almost exclusively disk-fed when FO Aqr is brighter than V~14 and substantially stream-fed when it is not. Even in the low states, a grazing eclipse remains detectable, confirming the presence of a disk-like structure (but not necessarily a Keplerian accretion disk). We relate these various observations to theoretical predictions that during the low state, the system's accretion disk dissipates into a non-Keplerian ring of diamagnetic blobs. Finally, a new XMM-Newton observation from a high state in 2017 reveals an anomalously soft X-ray spectrum and diminished X-ray luminosity compared to pre-2016 observations.

Data from transit light curves, radial velocity and transit timing observations can be used to probe the interiors of exoplanets beyond the mean density, by measuring the Love numbers $h_2$ and $k_2$. The first indirect estimate of $k_2$ for an exoplanet from radial velocity and transit timing variations observations has been performed by taking advantage of the years-spanning baseline. Not a single measurement of $h_2$ has been achieved from transit light curves, mostly because the photometric precision of current observing facilities is still too low. We show that the Imaging Spectrograph instrument on-board the Hubble Space Telescope could measure $h_2$ of the hot Jupiter WASP-121b if only few more observations were gathered. We show that a careful treatment of the noise and stellar limb darkening must be carried out to achieve a measurement of $h_2$. In particular, we find that the impact of the noise modelling on the estimation of $h_2$ is stronger than the impact of the limb darkening modelling. In addition, we emphasize that the wavelet method for correlated noise analysis can mask limb brightening. Finally, using presently available data, we briefly discuss the tentative measurement of $h_2 = 1.39^{+0.71}_{-0.81}$ in terms of interior structure. Additional observations would further constrain the interior of WASP-121b and possibly provide insights on the physics of inflation. The possibility of using the approach presented here with the Hubble Space Telescope provides a bridge before the high-quality data to be returned by the James Webb Space Telescope and PLATO telescope in the coming decade.

The evolution of the interior of stagnant-lid bodies is comparatively easier to model and predict with respect to the Earth's due to the absence of the large uncertainties associated with the physics of plate tectonics, its onset time and efficiency over the planet's history. Yet, the observational record for these bodies is both scarcer and sparser with respect to the Earth's. It is restricted to a limited number of samples and a variety of remote-sensing measurements of billions-of-years-old surfaces whose actual age is difficult to determine precisely. Combining these observations into a coherent picture of the thermal and convective evolution of the planetary interior represents thus a major challenge. In this chapter, we review key processes and (mostly geophysical) observational constraints that can be used to infer the global characteristics of mantle convection and thermal evolution of the interior of Mercury, the Moon and Mars, the three major terrestrial bodies where a stagnant lid has likely been present throughout most of their history.

We explore the application of machine learning based on mixture density neural networks (MDNs) to the interior characterization of low-mass exoplanets up to 25 Earth masses constrained by mass, radius, and fluid Love number $k_2$. We create a dataset of 900$\:$000 synthetic planets, consisting of an iron-rich core, a silicate mantle, a high-pressure ice shell, and a gaseous H/He envelope, to train a MDN using planetary mass and radius as inputs to the network. For this layered structure, we show that the MDN is able to infer the distribution of possible thicknesses of each planetary layer from mass and radius of the planet. This approach obviates the time-consuming task of calculating such distributions with a dedicated set of forward models for each individual planet. While gas-rich planets may be characterized by compositional gradients rather than distinct layers, the method presented here can be easily extended to any interior structure model. The fluid Love number $k_2$ bears constraints on the mass distribution in the planets' interior and will be measured for an increasing number of exoplanets in the future. Adding $k_2$as an input to the MDN significantly decreases the degeneracy of the possible interior structures.

We are witness to a great and increasing interest in internal structure, composition and evolution of exoplanets. However, direct measurements of exoplanetary mass and radius cannot be uniquely interpreted in terms of interior structure, justifying the need for an additional observable. The second degree fluid Love number, $k_2$, is proportional to the concentration of mass towards the body's center, hence providing valuable additional information about the internal structure. When hydrostatic equilibrium is assumed for the planetary interior, $k_2$ is a direct function of the planetary shape. Previous attempts were made to link the observed tidally and rotationally induced planetary oblateness in photometric light curves to $k_2$ using ellipsoidal shape models. Here, we construct an analytical 3D shape model function of the true planetary mean radius, that properly accounts for tidal and rotational deformations. Measuring the true planetary mean radius is critical when one wishes to compare the measured $k_2$ to interior theoretical expectations. We illustrate the feasibility of our method and show, by applying a Differential Evolution Markov Chain to synthetic data of WASP-121b, that a precision $\leq$ 65 ppm/$\sqrt{2\,min}$ is required to reliably retrieve $k_2$ with present understanding of stellar limb darkening, therefore improving recent results based on ellipsoidal shape models. Any improvement on stellar limb darkening would increase such performance.

Context. The mass and radius of a planet directly provide its bulk density, which can be interpreted in terms of its overall composition. Any measure of the radial mass distribution provides a first step in constraining the interior structure. The fluid Love number $k_2$ provides such a measure, and estimates of $k_2$ for extrasolar planets are expected to be available in the coming years thanks to improved observational facilities and the ever-extending temporal baseline of extrasolar planet observations. Aims. We derive a method for calculating the Love numbers $k_n$ of any object given its density profile, which is routinely calculated from interior structure codes. Methods. We used the matrix-propagator technique, a method frequently used in the geophysical community. Results. We detail the calculation and apply it to the case of GJ 436b, a classical example of the degeneracy of mass-radius relationships, to illustrate how measurements of $k_2$ can improve our understanding of the interior structure of extrasolar planets. We implemented the method in a code that is fast, freely available, and easy to combine with preexisting interior structure codes. While the linear approach presented here for the calculation of the Love numbers cannot treat the presence of nonlinear effects that may arise under certain dynamical conditions, it is applicable to close-in gaseous extrasolar planets like hot Jupiters, likely the first targets for which $k_2$ will be measured.

We describe the current state of knowledge about Mercury's interior structure. We review the available observational constraints, including mass, size, density, gravity field, spin state, composition, and tidal response. These data enable the construction of models that represent the distribution of mass inside Mercury. In particular, we infer radial profiles of the pressure, density, and gravity in the core, mantle, and crust. We also examine Mercury's rotational dynamics and the influence of an inner core on the spin state and the determination of the moment of inertia. Finally, we discuss the wide-ranging implications of Mercury's internal structure on its thermal evolution, surface geology, capture in a unique spin-orbit resonance, and magnetic field generation.

We observed RZ LMi, which is renowned for the extremely (~19d) short supercycle and is a member of a small, unusual class of cataclysmic variables called ER UMa-type dwarf novae, in 2013 and 2016. In 2016, the supercycles of this object substantially lengthened in comparison to the previous measurements to 35, 32, 60d for three consecutive superoutbursts. We consider that the object virtually experienced a transition to the novalike state (permanent superhumper). This observed behavior extremely well reproduced the prediction of the thermal-tidal instability model. We detected a precursor in the 2016 superoutburst and detected growing (stage A) superhumps with a mean period of 0.0602(1)d in 2016 and in 2013. Combined with the period of superhumps immediately after the superoutburst, the mass ratio is not as small as in WZ Sge-type dwarf novae, having orbital periods similar to RZ LMi. By using least absolute shrinkage and selection operator (Lasso) two-dimensional power spectra, we detected possible negative superhumps with a period of 0.05710(1)d. We estimated the orbital period of 0.05792d, which suggests a mass ratio of 0.105(5). This relatively large mass ratio is even above ordinary SU UMa-type dwarf novae, and it is also possible that the exceptionally high mass-transfer rate in RZ LMi may be a result of a stripped core evolved secondary which are evolving toward an AM CVn-type object.

Z Cam stars are a small subset of dwarf novae that exhibit standstills in their light curves. Most modern literature and catalogs of cataclysmic variables quote the number of known Z Cams to be on the order of 30 or so systems. After a four-year observing campaign and an exhaustive examination of the data in the AAVSO International Database we have trimmed that number by a third. One of the reasons for the misclassification of some systems is the fact that the definition of what a Z Cam is has changed over the last 85 years to what it is today. This has caused many stars formerly assumed to be Z Cams or rumored to be Z Cams to be eliminated from the final list. In this paper we present the results of our investigation into 65 stars listed at one time or another in the literature as Z Cams or possible Z Cams.

We carried out the photometric observations of the SU UMa-type dwarf nova ER UMa during 2011 and 2012, which showed the existence of persistent negative superhumps even during the superoutburst. We performed two-dimensional period analysis of its light curves by using a method called "least absolute shrinkage and selection operator" (Lasso) and "phase dispersion minimization" (PDM) analysis, and we found that the period of negative superhumps systematically changed between a superoutburst and the next superoutburst. The trend of the period change can beinterpreted as reflecting the change of the disk radius. This change of the disk radius is in good agreement with the predicted change of the disk radius by the thermal-tidal instability (TTI) model. The normal outbursts within a supercycle showed a general trend that the rising rate to maximum becomes slower as the next superoutburst approaches. The change can be interpreted as the consequence of the increased gas-stream flow onto the inner region of the disk as the result of the tilted disk. Some of the superoutbursts were found to be triggered by a precursor normal outburst when the positive superhumps appeared to develop. The positive and negative superhumps co-existed during the superoutburst. The positive superhumps were prominent only during four or five days after the supermaximum, while the signal of the negative superhumps became strong after the middle phase of the superoutburst plateau. A simple combination of the positive and negative superhumps was found to be insufficient in reproducing the complex profile variation. We were able to detect the developing phase of positive superhumps (stage A superhumps) for the first time in ER UMa-type dwarf novae. Using the period of stage A superhumps, we obtained a mass ratio of 0.100(15), which indicates that ER UMa is on the ordinary evolutional track of CVs.

Continuing the project described by Kato et al. (2009a, arXiv:0905.1757), we collected times of superhump maxima for SU UMa-type dwarf novae mainly observed during the 2012-2013 season. We found three objects (V444 Peg, CSS J203937 and MASTER J212624) having strongly positive period derivatives despite the long orbital period (Porb). By using the period of growing stage (stage A) superhumps, we obtained mass ratios for six objects. We characterized nine new WZ Sge-type dwarf novae. We made a pilot survey of the decline rate of slowly fading part of SU UMa-type and WZ Sge-type outbursts. The decline time scale was found to generally follow the expected Porb^(1/4) dependence and WZ Sge-type outbursts also generally follow this trend. There are some objects which show slower decline rates, and we consider these objects good candidates for period bouncers. We also studied unusual behavior in some objects, including BK Lyn which made a transition from an ER UMa-type state to the novalike (standstill) state in 2013 and unusually frequent occurrence of superoutbursts in NY Ser and CR Boo. We applied least absolute shrinkage and selection operator (Lasso) power spectral analysis, which has been proven to be very effective in analyzing the Kepler data, to ground-based photometry of BK Lyn and detected the dramatic disappearance of the signal of negative superhumps in 2013. We suggested that the mass-transfer rates did not vary strongly between the ER UMa-type state and novalike state in BK Lyn, and this transition was less likely caused by a systematic variation of the mass-transfer rate.

We report on photometric observations of WZ Sge-type dwarf novae, MASTER OT J211258.65+242145.4 and MASTER OT J203749.39+552210.3 which underwent outbursts in 2012. Early superhumps were recorded in both systems. During superoutburst plateau, ordinary superhumps with a period of 0.060291(4) d (MASTER J211258) and of 0.061307(9) d (MASTER J203749) in average were observed. MASTER J211258 and MASTER J203749 exhibited eight and more than four post-superoutburst rebrightenings, respectively. In the final part of the superoutburst, an increase in the superhump periods was seen in both systems. We have made a survey of WZ Sge-type dwarf novae with multiple rebrightenings, and confirmed that the superhump periods of WZ Sge-type dwarf novae with multiple rebrightenings were longer than those of WZ Sge-type dwarf novae without a rebrightening. Although WZ Sge-type dwarf novae with multiple rebrightenings have been thought to be the good candidates for period bouncers based on their low mass ratio (q) from inferred from the period of fully grown (stage B) superhumps, our new method using the period of growing superhumps (stage A superhumps), however, implies higher q than those expected from stage B superhumps. These q values appear to be consistent with the duration of the stage A superoutbursts, which likely reflects the growth time of the 3:1 resonance. We present a working hypothesis that the small fractional superhump excesses for stage B superhumps in these systems may be explained as a result that a higher gas pressure effect works in these systems than in ordinary SU UMa-type dwarf novae. This result leads to a new picture that WZ Sge-type dwarf novae with multiple rebrightenings and SU UMa-type dwarf novae without a rebrightening (they are not period bouncers) are located in the same place on the evolutionary track.