In this Review, I discuss recent developments on the long-term dynamical evolution of exoplanet systems, focusing on how distinctive dynamical processes may have shaped the orbital architectures of observed populations. I include three applications that highlight part of my own work. First, I examine the high-eccentricity tidal migration of hot Jupiters from a phase of dynamical instability and subsequent secular interactions in two-planet systems. Second, secular chaos as the origin of ultra-short-period planets with extreme period ratios. Third, secular resonance sweeping driven by a dispersing protoplanetary disk as the origin hot Neptunes residing in polar orbits. Finally, I discuss how upcoming observations will allow further constraining the prevalence of these dynamical processes.

In the last five years, the number of periodic variable stars has increased by two million. We used the ZTF DR2 data to find and build a catalog that includes 780,000 periodic variable stars. These periodic variable stars were classified into 11 types, which greatly complemented the variable stars in Galactic disk. Based on the latest ZTF DR16 data, we found 2 million variable candidates. We trained a machine learner to classify variable stars, and the learner had a prediction accuracy of 94%. Using millions of variable stars, we carried out studies to optimize the period–luminosity relations and the Galactic structure and the extinction law. With the future China Space Station Telescope, millions of variable stars in the Local Group will be discovered. They help to study the structure of our Local Group and also to cross-check the distance ladders based on different variable stars.

The lightcurves of strongly lensed AGNs get distorted due to gravitational microlensing, which differently magnifies the emission regions of AGNs depending on their size. This effect has been used to measure the size of the AGN accretion disc, but high photometric accuracy lightcurves reveal coherent variations on short time-scales that are not expected by standard accretion disc models. I show that this signal can be produced by emission from the Broad Line Region (BLR) but also by extended (diffuse) continuum emission. I explain how these features can be used to measure the size of the BLR but also reveal additional sources of emission. The multi-colour lightcurves of lensed AGNs, such as those to be obtained with the Vera Rubin Observatory, may become a powerful new tool to reveal the sub-parsec structure of AGNs, and shed light on elusive AGN emitting regions such as the one producing diffuse continuum emission.

We are investigating the extended outflow from G25.82–W1, which is one of the members of the high-mass protocluster G25.82–0.17. The aim is to study the star-forming environment of G25.82–W1. To identify the outflow, we obtained CO 2-1 data using the Atacama Large Millimeter/submillimeter Array.

We have identified several spatial and spectral outflows, including: 1) an extended N1–S1 CO outflow, driven by a high-mass young stellar object (HM-YSO) named G25.82–W1; 2) an elongated SE–NW outflow powered by G25.82–W2; 3) a compact and curved N2–S2 CO outflow originating from G25.82–E; and 4) a pair of knotty lobes centered on G25.82–W.

Furthermore, the innermost region of the N1–S1 CO outflow, traced by the 22 GHz H2O maser, reveals a complex spatial and velocity structure within a 2” from its launching point.

To accurately calculate the properties of the N1–S1 CO outflow, we have utilized an accurate distance measurement of d=4.5 kpc, derived from the annual parallax of the H2O masers. The outflow rate and force are comparable to those observed in outflows from other HM-YSOs. The physical properties of the N1–S1 CO outflow follow a trend connecting the low and high-mass regimes, supporting the idea that the star-forming mode in G25.82–W1 is likely a scaled-up version of low-mass star formation.

We present a progress report of our project aiming to increase the number of known Cepheids in double-lined binary (SB2) systems from six to 100 or more. This will allow us, among other goals, to accurately measure masses for a large sample of Cepheids. Currently, only six accurate Cepheid masses are available, which hinders our understanding of their physical properties and renders the Cepheid mass–luminosity relation poorly constrained. At the same time, Cepheids are widely used for essential measurements (e.g., extragalactic distances, the Hubble constant). To examine Cepheid period–luminosity relations, we selected as binary candidates Cepheids that are too bright for their periods. To date, we have confirmed 56 SB2 systems, including the detection of significant orbital motions of the components for 32. We identified systems with orbital periods up to five times shorter than the shortest reported period to date, as well as systems with mass ratios significantly different from unity (suggesting past merger events). Both features are essential to understand how multiplicity affects the formation and destruction of Cepheid progenitors and what effect this has on global Cepheid properties. We also present eight new systems composed of two Cepheids (only one such system was known before). Among confirmed SB2 Cepheids, there are also several wide-orbit systems. In the future, these may facilitate independent accurate geometric distance measurements to the Large and Small Magellanic Clouds.

We discuss the impact of Gaia, the cornerstone mission of the European Space Agency (ESA), on the calibration of the period–luminosity and luminosity–metallicity relations of Cepheids and RR Lyrae stars, with specific reference to data published as part of the most recent Gaia releases: Early Data Release 3 (EDR3), on 19 December 2020, and Data Release 3 (DR3) on 13 June 2022. We provide future perspectives for the Gaia mission, including extensions approved by ESA and a tentative schedule of the data releases that will take place in the next few years. We briefly present plans for cross-Coordination Unit processing of Gaia data of Cepheids and RR Lyrae stars for DR4 and conclude by outlining the expected improvement in astrometry at the end of the extended Gaia mission, which will help to further strengthen the calibration of the first rung of the cosmic distance ladder.

We present a method to estimate distances to AGB stars, utilizing the rich infrared data sets available for these infrared-bright targets. The method is based on the assumption that stars with intrinsically similar properties (metallicity, initial mass, etc.) produce similar spectral energy distributions (SEDs) and similar luminosities. We here discuss the results for AGB stars belonging to the BAaDE survey sample whose distances were calibrated using the template SEDs of stars with their VLBI parallaxes. As VLBI parallaxes are only known for a handful of sources, the resulting templates only cover a small subset of the BAaDE sample. Additional methods to derive suitable templates will therefore also be required. The work on expanding the template set is promising, although more fine tuning is still needed.

In the context of a perturbed two body problem, in which the Keplerian motion of the small object (the satellite) is perturbed by the oblateness of the central body (the asteroid) and the attraction of a third body (the Sun), we discuss the long-term evolution of the orbital elements of a satellite orbiting an oblate body, with a particular focus on the behavior of the inclination and the longitude of the ascending node. We derive analytically the position of the Laplace plane as a function of several parameters and use this solution to analyse the long-term evolution of distant circular orbits. The analytical study is complemented by numerical tests, performed in the context of both Cartesian and Hamiltonian frameworks. The results give a description of the orbital dynamical environment of asteroids and reveal the parameters that play a key role in the long-term stability of distant circular orbits.

The MAximum-entropy ReconStruction (MARS) method is a free-form strong-lensing (SL) reconstruction algorithm, which adopts the maximum cross-entropy as a regularization. MARS shows remarkable convergence of multiple images in both source (∼0.”02) and image planes (∼0.”05 – 0.”1) while suppressing spurious fluctuations. Although the reconstruction requires a large number of free parameters exceeding ∼19,000, our implementation through PyTorch can obtain the reconstruction within hours. From our test using the publicly available synthetic clusters, we have verified that the reconstructed radial mass profiles are consistent with the truth within 1 percent. This makes MARS one of the best-performing SL reconstruction methods. We apply MARS to the six Hubble Frontier Fields clusters and present new mass reconstruction results. We also reconstruct a mass model of Abell 2744 using both weak-lensing (WL) and SL data from the JWST observations, with the largest dataset of Abell 2744, including 286 SL multiple images and ∼350 arcmin−2 WL constraints.

Forthcoming data from the Vera Rubin Observatory, Euclid and Roman telescopes are expected to increase the number of strong lenses by two orders of magnitude. With current discovery methods these would be accompanied by an even greater number of false positives. In that context we find that using an ensemble of classifiers would provide a more complete sample of high-purity lenses and present methods to post-process the outputs of such classifiers to give reliable probabilities that a given image contains a lens.

A key ingredient in the earliest evolutionary phase of high-mass (M>8 M⊙) star formation (HMSF) is the presence of a jet/outflow system. To study its role in HMSF, we have carried out high resolution (0.1″) VLA K-band (18-26.5 GHz) observations toward IRAS 19035+0641 A, identified as a high-mass protostellar jet candidate based on previous cm continuum data. Our observations resolve the continuum emission into an elongated structure in the NE-SW direction, confirming that the K-band continuum arises from an ionized jet. Furthermore, we detected several 22.2 GHz H2O maser spots aligned in a direction consistent with the jet axis. Zeeman splitting was detected in the strongest maser spot. In this paper, we present our results and discuss the implications of our findings.

Recent ALMA observations detected protostellar outflows in 70-μm dark infrared dark clouds (IRDCs). These sources are candidates for the initial stages of high-mass star formation. We launched a new survey for free-free emission from outflow shocks using the Yamaguchi Interferometer (YI) at 8 GHz. We aim to catalog “proto-high-mass protostar” candidates that are still in the low to intermediate-mass phase. We selected starless-like clumps without any 70-μm point source from Traficante et al. (2015). We currently detected 82 sources from 167 clumps. 37 of them are fainter than 20 mJy (down to a few mJy). They tend to associate with colder and denser clumps that are suitable for star formation. This fact suggests that, at least, some of them trace star-formation activities. The highest-density clumps are, in fact, associated with several masers and molecular outflows. Furthermore, some of them have already shown a signature of ongoing cluster formation.

We briefly consider the history of maser variability, and of flaring variability specifically. We consider six proposed flare generation mechanisms, and model them computationally with codes that include saturation and 3-D structure (the last mechanism is modelled in 1-D). Fits to observational light curves have been made for some sources, and we suggest that a small number of observational parameters can diagnose the flare mechanism in many cases. The strongest flares arise from mechanisms that can increase the number density of inverted molecules in addition to by geometrical effects, and in events where unsaturated quiescent masers become saturated during the flare.

We suggest an advanced algorithm for semi-analytical calculation of orbital perturbations of Earth artificial satellites caused by the gravity attraction of the “3rd-bodies” (the Moon, the Sun, major planets). A new accurate analytical series for the relevant perturbation function is developed. It is obtained through a careful spectral analysis of the long-term DE406 planetary/lunar ephemerides and valid over 2000 years, 1000-3000. The series is used in the author’s semi-analytical model of satellite motion. The results of the motion prediction of several Earth satellites obtained by means of the semi-analytical model and a numerical integration method are compared.

Systems near mean-motion resonances (MMRs) are subject to large transit-timing variations (TTVs). The amplitude and period of the TTVs strongly depend on the distance to exact MMR and the planetary eccentricities which are shaped during the formation and long-term evolution of the system. For close-in planets, the tides raised by the star provide a source of dissipation, placing the planets further away from the MMR. In this work, we will discuss how the tidal interactions with the central star play an important role in shaping the period ratios and resonant angles in resonant chains. Moreover, we will show how they can impact the TTVs and therefore how the TTVs could serve as a means to put constraints on the tidal history of planetary systems. The study will focus on the four-planet resonant chain of Kepler-80.

Close-in planets undergo strong tidal effects with their host stars that modify their spins and orbits. Adopting a Maxwell rheology, it has been shown that for the 5/2 and 7/2 spin-orbit resonances, the obliquity of these planets can stabilise at a high value. Here, we show that these high obliquity metastable states can also be observed for the same spin-orbit resonances considering the Andrade rheology.

The database of circumstellar OH masers by Engels & Bunzel (2015) was updated to include new 1612, 1665, and 1667 MHz OH maser observations published between 2015 and 2022. A cross-correlation of the database was made with infrared catalogues (AllWISE and 2MASS) and with GAIA DR3. This led frequently to significantly improved coordinates and identified contaminations with non-stellar sources. About 40% of all OH maser-detected stars were not detected by GAIA. These are mostly representatives of the population of highly obscured stars at the end of AGB or at the beginning of post-AGB evolution.

Cepheids have been the cornerstone of the extragalactic distance scale for a century. With high-quality data, these luminous supergiants exhibit a small dispersion in their Leavitt (period–luminosity) relation, particularly at longer wavelengths, and few methods rival the precision possible with Cepheid distances. In these proceedings, we present an overview of major observational programs pertaining to the Cepheid extragalactic distance scale, its progress and remaining challenges. In addition, we present preliminary new results on Cepheids from the James Webb Space Telescope (JWST). The launch of JWST has opened a new chapter in the measurement of extragalactic distances and the Hubble constant. JWST offers a resolution three times that of the Hubble Space Telescope (HST) with nearly 10 times the sensitivity. It has been suggested that the discrepancy in the value of the Hubble constant based on Cepheids compared to that inferred from measurements of the cosmic microwave background requires new and additional physics beyond the standard cosmological model. JWST observations will be critical in reducing remaining systematics in the Cepheid measurements and for confirming if new physics is indeed required. Early JWST data for the galaxy, NGC 7250 show a decrease in scatter in the Cepheid Leavitt law by a factor of two relative to existing HST data and demonstrate that crowding/blending effects are a significant issue in a galaxy as close as 20 Mpc.