The equations of motion of planar elliptic restricted three body problem are transformed to four decoupled Hill's equations. By using the Floquet theorem analytic solution to the oscillator equations with time dependent periodic coefficients are presented. We show that the new analytic approach is valid for system parameters $0 < e \leq 0.05$ and $0 < \mu \leq 0.01$ where $e$ denotes the eccentricity of primaries while $\mu$ is the mass parameter, respectively. We also clarify the transformation details that provide the applicability of the method.

Wide field surveys for transiting planets are well suited to searching diverse stellar populations, enabling a better understanding of the link between the properties of planets and their parent stars. We report the discovery of HAT-P-69b (TOI 625.01) and HAT-P-70b (TOI 624.01), two new hot Jupiters around A stars from the HATNet survey which have also been observed by the Transiting Exoplanet Survey Satellite (TESS). HAT-P-69b has a mass of 3.58 +0.58/-0.58 MJup and a radius of 1.676 +0.051/-0.033 RJup, residing in a prograde 4.79-day orbit. HAT-P-70b has a radius of 1.87 +0.15/-0.10 RJup and a mass constraint of < 6.78 (3 sigma) MJup, and resides in a retrograde 2.74-day orbit. We use the confirmation of these planets around relatively massive stars as an opportunity to explore the occurrence rate of hot Jupiters as a function of stellar mass. We define a sample of 47,126 main-sequence stars brighter than Tmag=10 that yields 31 giant planet candidates, including 18 confirmed planets, 3 candidates, and 10 false positives. We find a net hot Jupiter occurrence rate of 0.41+/-0.10 % within this sample, consistent with the rate measured by Kepler for FGK stars. When divided into stellar mass bins, we find the occurrence rate to be 0.71+/-0.31% for G stars, 0.43+/-0.15% for F stars, and 0.26+/-0.11% for A stars. Thus, at this point, we cannot discern any statistically significant trend in the occurrence of hot Jupiters with stellar mass.

The particle dynamics of the recently observed ring around dwarf planet Haumea is numerically investigated. The point mass gravitational force, a second degree and order gravity field, and the solar radiation pressure as the main perturbations are considered. The quasi-stationary state of the ring varies for different micron-sized grains and depends also on the spin-orbit resonances between the rotation rate of the main body and the orbital period of the dust particles. The simulations confirm the variable radial width of the ring observed during the transit ingress and egress. Results show that the micron sized grains, initially on circular orbits, become eccentric and form an apse-aligned ring at the observed radial distance near to the 3:1 spin-orbit resonance. It is also demonstrated that this coincidence is only apparent and independent of 3:1 resonance.

We formulate a model that describes the escape dynamics in a leaky chaotic system in which the size of the leak depends on the number of the in-falling particles. The basic motivation of this work is the astrophysical process which describes the planetary accretion. In order to study the dynamics generally, the standard map is investigated in two cases when the dynamics is fully hyperbolic and in the presence of KAM islands. In addition to the numerical calculations, an analytic solution to the temporal behavior of the model is also derived. We show that in the early phase of the leak expansion, as long as there are enough particles in the system, the number of survivors deviates from the well-known exponential decay. Furthermore, the analytic solution returns the classical result in the limiting case when the number of particles does not affect the leak size.

We present the first simultaneous photometric and spectroscopic investigation of a large set of RR Lyrae variables in a globular cluster. The radial-velocity data presented comprise the largest sample of RVs of RR Lyrae stars ever obtained. The target is M3; $BVI_{\mathrm{C}}$ time-series of 111 and $b$ flux data of further 64 RRab stars, and RV data of 79 RR Lyrae stars are published. Blazhko modulation of the light curves of 47 percent of the RRab stars are detected. The mean value of the center-of-mass velocities of RR Lyrae stars is $-146.8$ km s$^{-1}$ with 4.52 km s$^{-1}$ standard deviation, which is in good agreement with the results obtained for the red giants of the cluster. The ${\Phi_{21}}^{\mathrm RV}$ phase difference of the RV curves of RRab stars is found to be uniformly constant both for the M3 and for Galactic field RRab stars; no period or metallicity dependence of the ${\Phi_{21}}^{\mathrm RV}$ is detected. The Baade-Wesselink distances of 26 non-Blazhko variables with the best phase-coverage radial-velocity curves are determined; the corresponding distance of the cluster, $10480\pm210$ pc, agrees with the previous literature information. A quadratic formula for the $A_{\mathrm{puls}}-A_V$ relation of RRab stars is given, which is valid for both OoI and OoII variables. We also show that the $(V-I)_0$ of RRab stars measured at light minimum is period dependent, there is at least 0.1 mag difference between the colours at minimum light of the shortest- and longest-period variables.

We present the discovery of the transiting exoplanets HAT-P-65b and HAT-P-66b, with orbital periods of 2.6055 d and 2.9721 d, masses of $0.527 \pm 0.083$ M$_{J}$ and $0.783 \pm 0.057$ M$_{J}$ and inflated radii of $1.89 \pm 0.13$ R$_{J}$ and $1.59^{+0.16}_{-0.10}$ R$_{J}$, respectively. They orbit moderately bright ($V=13.145 \pm 0.029$, and $V=12.993 \pm 0.052$) stars of mass $1.212 \pm 0.050$ M$_{\odot}$ and $1.255^{+0.107}_{-0.054}$ M$_{\odot}$. The stars are at the main sequence turnoff. While it is well known that the radii of close-in giant planets are correlated with their equilibrium temperatures, whether or not the radii of planets increase in time as their hosts evolve and become more luminous is an open question. Looking at the broader sample of well-characterized close-in transiting giant planets, we find that there is a statistically significant correlation between planetary radii and the fractional ages of their host stars, with a false alarm probability of only 0.0041%. We find that the correlation between the radii of planets and the fractional ages of their hosts is fully explained by the known correlation between planetary radii and their present day equilibrium temperatures, however if the zero-age main sequence equilibrium temperature is used in place of the present day equilibrium temperature then a correlation with age must also be included to explain the planetary radii. This suggests that, after contracting during the pre-main-sequence, close-in giant planets are re-inflated over time due to the increasing level of irradiation received from their host stars. Prior theoretical work indicates that such a dynamic response to irradiation requires a significant fraction of the incident energy to be deposited deep within the planetary interiors.

HBC 722 (V2493 Cyg) is a young eruptive star in outburst since 2010. It is an FU Orionis-type object with an atypically low outburst luminosity. Because it was well characterized in the pre-outburst phase, HBC 722 is one of the few FUors where we can learn about the physical changes and processes associated with the eruption. We monitored the source in the BVRIJHKs bands from the ground, and at 3.6 and 4.5 $\mu$m from space with the Spitzer Space Telescope. We analyzed the light curves and the spectral energy distribution by fitting a series of steady accretion disk models at many epochs. We also analyzed the spectral properties of the source based on new optical and infrared spectra. We also mapped HBC 722 and its surroundings at millimeter wavelengths. From the light curve analysis we concluded that the first peak of the outburst in 2010 September was due to an abrupt increase of the accretion rate in the innermost part of the system. This was followed by a long term process, when the brightening was mainly due to a gradual increase of the accretion rate and the emitting area. Our new observations show that the source is currently in a constant plateau phase. We found that around the peak the continuum was bluer and the H$\alpha$ profile changed significantly between 2012 and 2013. The source was not detected in the millimeter continuum, but we discovered a flattened molecular gas structure with a diameter of 1700 au and mass of 0.3 M$_{\odot}$ centered on HBC 722. While the first brightness peak could be interpreted as a rapid fall of piled-up material from the inner disk onto the star, the later monotonic flux rise suggests the outward expansion of a hot component according to the theory of Bell & Lin (1994). Our study of HBC 722 demonstrated that accretion-related outbursts can occur in young stellar objects even with very low mass disks, in the late Class II phase.

We report the discovery of two new transiting extrasolar planets orbiting moderately bright (V = 10.7 and 12.2 mag) F stars (masses of 1.39 Msun and 1.10 Msun, respectively). The planets have periods of P = 4.7322 d and 4.4087 d, and masses of 0.21 MJ and 0.17 MJ which are almost half-way between those of Neptune and Saturn. With radii of 1.31 RJ and 1.13 RJ, these very low density planets are the two lowest mass planets with radii in excess that of Jupiter. Comparing with other recent planet discoveries, we find that sub-Saturns (0.18MJ < Mp < 0.3MJ) and super-Neptunes (0.05MJ < Mp < 0.18MJ) exhibit a wide range of radii, and their radii exhibit a weaker correlation with irradiation than higher mass planets. The two planets are both suitable for measuring the Rossiter-McLaughlin effect and for atmospheric characterization. Measuring the former effect would allow an interesting test of the theory that star-planet tidal interactions are responsible for the tendency of close-in giant planets around convective envelope stars to be on low obliquity orbits. Both planets fall on the edge of the short period Neptunian desert in the semi-major axis-mass plane.

We report the discovery of a new transiting extrasolar planet, HAT-P-55b. The planet orbits a V = 13.207 +/- 0.039 sun-like star with a mass of 1.013 +/- 0.037 solar masses, a radius of 1.011 +/- 0.036 solar radii and a metallicity of -0.03 +/- 0.08. The planet itself is a typical hot Jupiter with a period of 3.5852467 +/- 0.0000064 days, a mass of 0.582 +/- 0.056 Jupiter masses and a radius of 1.182 +/- 0.055 Jupiter radii. This discovery adds to the increasing sample of transiting planets with measured bulk densities, which is needed to put constraints on models of planetary structure and formation theories.

The overtone and multi-mode RR Lyrae stars in the globular cluster M3 are studied using a 200-d long, $B,V$ and $I_{\mathrm C}$ time-series photometry obtained in 2012. 70\% of the 52 overtone variables observed show some kind of multi-periodicity (additional frequency at ${f_{0.61}}={f_{\mathrm {1O}}}/0.61$ frequency ratio, Blazhko effect, double/multi-mode pulsation, period doubling). A signal at 0.587 frequency ratio to the fundamental-mode frequency is detected in the double-mode star, V13, which may be identified as the second radial overtone mode. If this mode-identification is correct, than V13 is the first RR Lyrae star showing triple-mode pulsation of the first three radial modes. Either the Blazhko effect or the ${f_{0.61}}$ frequency (or both of these phenomena) appear in 7 double-mode stars. The $P_{\mathrm{1O}}/P_{\mathrm{F}}$ period ratio of RRd stars showing the Blazhko effect are anomalous. A displacement of the main frequency component at the fundamental-mode with the value of modulation frequency (or its half) is detected in three Blazhko RRd stars parallel with the appearance of the overtone-mode pulsation. The ${f_{0.61}}$ frequency appears in RRc stars that lie at the blue side of the double-mode region and in RRd stars, raising the suspicion that its occurrence may be connected to double-mode pulsation. The changes of the Blazhko and double-mode properties of the stars are also reviewed using the recent and archive photometric data.

We report the discovery and characterization of four transiting exoplanets by the HATNet survey. The planet HAT-P-50b has a mass of 1.35 M_J and a radius of 1.29 R_J, and orbits a bright (V = 11.8 mag) M = 1.27 M_sun, R = 1.70 R_sun star every P = 3.1220 days. The planet HAT-P-51b has a mass of 0.31 M_J and a radius of 1.29 R_J, and orbits a V = 13.4 mag, M = 0.98 M_sun, R = 1.04 R_sun star with a period of P = 4.2180 days. The planet HAT-P-52b has a mass of 0.82 M_J and a radius of 1.01 R_J, and orbits a V = 14.1 mag, M = 0.89 M_sun, R = 0.89 R_sun star with a period of P = 2.7536 days. The planet HAT-P-53b has a mass of 1.48 M_J and a radius of 1.32 R_J, and orbits a V = 13.7 mag, M = 1.09 M_sun, R = 1.21 R_sun star with a period of P = 1.9616 days. All four planets are consistent with having circular orbits and have masses and radii measured to better than 10% precision. The low stellar jitter and favorable R_P/R_star ratio for HAT-P-51 make it a promising target for measuring the Rossiter-McLaughlin effect for a Saturn-mass planet.

The concept of Space Manifold Dynamics is a new method of space research. We have applied it along with the basic idea of the method of Ott, Grebogi and York (OGY method) to stabilize the motion of a spacecraft around the triangular Lagrange point L5 of the Earth--Moon system. We have determined the escape rate of the trajectories in the general three- and four-body problem and estimated the average lifetime of the particles. Integrating the two models we mapped in detail the phase space around the L5 point of the Earth--Moon system. Using the phase space portrait our next goal was to apply a modified OGY method to keep a spacecraft close to the vicinity of L5. We modified the equation of motions with the addition of a time dependent force to the motion of the spacecraft. In our orbit--keeping procedure there are three free parameters: (i) the magnitude of the thrust, (ii) the start time and (iii) the length of the control. Based on our numerical experiments we were able to determine possible values for these parameters and successfully apply a control phase to a spacecraft to keep it on orbit around L5.

The circular restricted three body problem is investigated in the context of accretion and scattering processes. In our model a large number of identical non-interacting mass-less planetesimals are considered in planar case orbiting a star-planet system. This description allows us to investigate in dynamical systems approach the gravitational scattering and possible captures of the particles by the forming planetary embryo. Although the problem serves a large variety of complex motion, the results can be easily interpreted because of the low dimensionality of the phase space. We show that initial conditions define isolated regions of the disk, where accretion or escape of the planetesimals occur, these have, in fact, a fractal structure. The fractal geometry of these "basins" implies that the dynamics is very complex. Based on the calculated escape rates and escape times, it is also demonstrated that the planetary accretion rate is exponential for short times and follows a power-law for longer integration. A new numerical calculation of the maximum mass that a planet can reach (described by the expression of the isolation mass) is also derived.

We report the discovery of HAT-P-54b, a planet transiting a late K dwarf star in field 0 of the NASA K2 mission. We combine ground-based photometric light curves with radial velocity measurements to determine the physical parameters of the system. HAT-P-54b has a mass of 0.760 $\pm$ 0.032 $M_J$, a radius of 0.944 $\pm$ 0.028 $R_J$, and an orbital period of 3.7998 d. The star has V = 13.505 $\pm$ 0.060, a mass of 0.645 $\pm$ 0.020 $M_{\odot}$, a radius of 0.617 $\pm$ 0.013 $R_{\odot}$, an effective temperature of Teff = 4390 $\pm$ 50K, and a subsolar metallicity of [Fe/H] = -0.127 $\pm$ 0.080. HAT-P-54b has a radius that is smaller than 92% of the known transiting planets with masses greater than that of Saturn, while HAT-P-54 is one of the lowest-mass stars known to host a hot Jupiter. Follow-up high-precision photometric observations by the K2 mission promise to make this a well-studied planetary system.

HD 15082 (WASP-33) is the hottest and fastest rotating star known to harbor a transiting extrasolar planet (WASP-33b). The lack of high precision radial velocity (RV) data stresses the need for precise light curve analysis and gathering further RV data. By using available photometric and RV data, we perform a blend analysis, compute more accurate system parameters, confine the planetary mass and attempt to cast light on the observed transit anomalies. We combine the original HATNet observations and various followup data to jointly analyze the signal content and extract the transit component and use our RV data to aid the global parameter determination. The blend analysis of the combination of multicolor light curves yields the first independent confirmation of the planetary nature of WASP-33b. We clearly identify three frequency components in the 15-21 1/day regime with amplitudes 7-5 mmag. These frequencies correspond to the delta Scuti-type pulsation of the host star. None of these pulsation frequencies or their low-order linear combinations are in close resonance with the orbital frequency. We show that these pulsation components explain some but not all of the observed transit anomalies. The grand-averaged transit light curve shows that there is a ~1.5 mmag brightening shortly after the planet passes the mid-transit phase. Although the duration and amplitude of this brightening varies, it is visible even through the direct inspections of the individual transit events (some 40-60% of the followup light curves show this phenomenon). We suggest that the most likely explanation of this feature is the presence of a well-populated spot belt which is highly inclined to the orbital plane. This geometry is consistent with the inference from the spectroscopic anomalies. Finally, we constrain the planetary mass to M_p=3.27+/-0.73 M_J by using our RV data collected by the TRES spectrograph.

Context. V2492 Cyg is a young eruptive star that went into outburst in 2010. The near-infrared color changes observed since the outburst peak suggest that the source belongs to a newly defined sub-class of young eruptive stars, where time-dependent accretion and variable line-of-sight extinction play a combined role in the flux changes. Aims. In order to learn about the origin of the light variations and to explore the circumstellar and interstellar environment of V2492 Cyg, we monitored the source at ten different wavelengths, between 0.55 \mu m and 2.2 \mu m from the ground and between 3.6 \mu m and 160 \mu m from space. Methods. We analyze the light curves and study the color-color diagrams via comparison with the standard reddening path. We examine the structure of the molecular cloud hosting V2492 Cyg by computing temperature and optical depth maps from the far-infrared data. Results. We find that the shapes of the light curves at different wavelengths are strictly self-similar and that the observed variability is related to a single physical process, most likely variable extinction. We suggest that the central source is episodically occulted by a dense dust cloud in the inner disk, and, based on the invariability of the far-infrared fluxes, we propose that it is a long-lived rather than a transient structure. In some respects, V2492 Cyg can be regarded as a young, embedded analog of UX Orionis-type stars. Conclusions. The example of V2492 Cyg demonstrates that the light variations of young eruptive stars are not exclusively related to changing accretion. The variability provided information on an azimuthally asymmetric structural element in the inner disk. Such an asymmetric density distribution in the terrestrial zone may also have consequences for the initial conditions of planet formation.

We investigate the phase space structure of the relativistic Sitnikov problem in the first post-Newtonian approximation. The phase space portraits show a strong dependence on the gravitational radius which describes the strength of the relativistic pericentre advance. Bifurcations appearing at increasing the gravitational radius are presented. Transient chaotic behavior related to escapes from the primaries are also studied. Finally, the numerically determined chaotic saddle is investigated in the context of hyperbolic and non-hyperbolic dynamics as a function of the gravitational radius.