Chaos in cuspy triaxial galaxies with a supermassive black hole: a simple toy model

This paper summarises an investigation of chaos in a toy potential which mimics much of the behaviour observed for the more realistic triaxial generalisations of the Dehnen potentials, which have been used to model cuspy triaxial galaxies both with and without a supermassive black hole. The potential is the sum of an anisotropic harmonic oscillator potential, ${\text{V}}_{\text{0}} = \frac{1}{2}\left( {a^2 x^2 + b^2 y^2 + c^2 z^2 } \right)$ , and aspherical Plummer potential, ${\text{V}}_{\text{P}} = M_{BH} /\sqrt {r^2 + \varepsilon ^2 } $ , with $r^2 = x^2 + y^2 + z^2$ . Attention focuses on three issues related tothe properties of ensembles of chaotic orbits which impact on chaotic mixing and the possibility of constructing self-consistent equilibria:(1) What fraction of the orbits are chaotic? (2) How sensitive are the chaotic orbits, that is, how large are their largest (short time) Lyapunov exponents? (3) To what extent is the motion of chaotic orbits impeded by Arnold webs, that is, how 'sticky' are the chaotic orbits? These questions are explored as functions of the axis ratio a: b: c, black hole mass M _BH, softening length ε, and energy E with the aims of understanding how the manifestations of chaos depend onthe shape of the system and why the black hole generates chaos. The simplicity of the model makes it amenable to a perturbative analysis. That it mimics the behaviour of more complicated potentials suggests that much of this behaviour should be generic.     

The motion around the libration points in the restricted three-body problem with the effect of radiation and oblateness

This paper deals with the existence of libration points and their linear stability when the more massive primary is radiating and the smaller is an oblate spheroid. Our study includes the effects of oblateness of $\bar{J}_{2i}$ (i=1,2) with respect to the smaller primary in the restricted three-body problem. Under combining the perturbed forces that were mentioned before, the collinear points remain unstable and the triangular points are stable for 0<μ<μ _c , and unstable in the range $\mu_{c} \le\mu\le\frac{1}{2}$ , where $\mu_{c} \in(0,\frac{1}{2})$ , it is also observed that for these points the range of stability will decrease. The relations for periodic orbits around five libration points with their semimajor, semiminor axes, eccentricities, the frequencies of orbits and periods are found, furthermore for the orbits around the triangular points the orientation and the coefficients of long and short periodic terms also are found in the range 0<μ<μ _c .     

Secular resonances between bodies on close orbits II: prograde and retrograde orbits for irregular satellites

In extending the analysis of the four secular resonances between close orbits in Li and Christou (Celest Mech Dyn Astron 125:133–160, 2016) (Paper I), we generalise the semianalytical model so that it applies to both prograde and retrograde orbits with a one-to-one map between the resonances in the two regimes. We propose the general form of the critical angle to be a linear combination of apsidal and nodal differences between the two orbits \( b_1 \Delta \varpi + b_2 \Delta \varOmega \), forming a collection of secular resonances in which the ones studied in Paper I are among the strongest. Test of the model in the orbital vicinity of massive satellites with physical and orbital parameters similar to those of the irregular satellites Himalia at Jupiter and Phoebe at Saturn shows that \({>}20\) and \({>}40\%\) of phase space is affected by these resonances, respectively. The survivability of the resonances is confirmed using numerical integration of the full Newtonian equations of motion. We observe that the lowest order resonances with \(b_1+|b_2|\le 3\) persist, while even higher-order resonances, up to \(b_1+|b_2|\ge 7\), survive. Depending on the mass, between 10 and 60% of the integrated test particles are captured in these secular resonances, in agreement with the phase space analysis in the semianalytical model.     

Constraining the Angular Momentum of the Sun with Planetary Orbital Motions and General Relativity

The angular momentum of a star is an important astrophysical quantity related to its internal structure, formation, and evolution. Helioseismology yields $S_{\odot}= 1.92\times10^{41}\ \mathrm{kg\ m^{2}\ s^{-1}}$ for the angular momentum of the Sun. We show how it should be possible to constrain it in a near future by using the gravitomagnetic Lense?CThirring effect predicted by General Relativity for the orbit of a test particle moving around a central rotating body. We also discuss the present-day situation in view of the latest determinations of the supplementary perihelion precession of Mercury. A fit by Fienga et al. (Celestial Mech. Dynamical Astron. 111, 363, 2011) of the dynamical models of several standard forces acting on the planets of the solar system to a long data record yielded milliarcseconds per century. The modeled forces did not include the Lense?CThirring effect itself, which is expected to be as large as from helioseismology-based values of S _??. By assuming the validity of General Relativity, from its theoretical prediction for the gravitomagnetic perihelion precession of Mercury, one can straightforwardly infer $S_{\odot}\leq0.95\times10^{41}\ \mathrm{kg\, m^{2}\, s^{-1}}$ . It disagrees with the currently available values from helioseismology. Possible sources for the present discrepancy are examined. Given the current level of accuracy in the Mercury ephemerides, the gravitomagnetic force of the Sun should be included in their force models. MESSENGER, in orbit around Mercury since March 2011, will collect science data until 2013, while BepiColombo, to be launched in 2015, should reach Mercury in 2022 for a year-long science phase: the analysis of their data will be important in effectively constraining S _?? in about a decade or, perhaps, even less.     

Stability of the triangular Lagrange points beyond Gascheau��s value

We examine the stability of the triangular Lagrange points L ₄ and L ₅ for secondary masses larger than the Gascheau??s value ${\mu_{\rm G}= (1-\sqrt{23/27}/2)= 0.0385208\ldots}$ (also known as the Routh value) in the restricted, planar circular three-body problem. Above that limit the triangular Lagrange points are linearly unstable. Here we show that between??? _G and ${\mu \approx 0.039}$ , the L ₄ and L ₅ points are globally stable in the sense that a particle released at those points at zero velocity (in the corotating frame) remains in the vicinity of those points for an indefinite time. We also show that there exists a family of stable periodic orbits surrounding L ₄ or L ₅ for ${\mu \ge \mu_G}$ . We show that??? _G is actually the first value of a series ${\mu_0 (=\mu_G), \mu_1,\ldots, \mu_i,\ldots}$ corresponding to successive period doublings of the orbits, which exhibit ${1, 2, \ldots, 2^i,\ldots}$ cycles around L ₄ or L ₅. Those orbits follow a Feigenbaum cascade leading to disappearance into chaos at a value ${\mu_\infty = 0.0463004\ldots}$ which generalizes Gascheau??s work.     

Disk precession and quasi-periodic brightness oscillations of V603 Aql in 2001–2002

We present the photometric observations of the old nova V603 Aql with the RTT 150 Russian-Turkish telescope during eleven nights of 2001–2002. We show that the star at this time was in a state with positive superhumps and its photometric period of \(0\mathop .\limits^d 144 - 0\mathop .\limits^d 145\) was longer than the orbital period. We found night-to-night variations in the mean brightness of the system that are consistent with disk precession periods of \(3\mathop .\limits^d 3\) and \(3\mathop .\limits^d 0\) in 2001 and 2002, respectively. Analysis of the results and their comparison with the results of other authors using current theoretical models for disk precession lead us to suggest that the change in the disk precession period was caused by a change in the accretion rate in the system. V603 Aql in a state with negative superhumps was found to be brighter than it is in a state with positive superhumps by \(0\mathop .\limits^m 2 - 0\mathop .\limits^m 3\). We hypothesize that the transition between these states could also be caused by a change in the accretion rate. Quasi-periodic oscillations (QPOs) of the brightness with typical time scales of 9–70 min were detected on each observing night. These time scales were found to change from night to night. The detection of QPOs with a period of about 0.05 of the orbital period and its multiples on certain nights provides evidence for the model of QPO generation through accretion-rate modulation by ionization-front oscillations on the surface of the donor star near the inner Lagrangian point.     

Secular period decreasing of detached chromospherically active binaries

The long-term orbital period changes of a large sample of detached chromospherically active binaries (CABs) were studied. Eleven such systems were found to be undergoing secular period decreases with the rates of ?6.3×10^?9 to ?1.1×10^?6 days per year. The period decreasing rates are found to vary depending on the orbital period. The longer the orbital period is, the more rapidly the period decreases. Following Stepien (Mon. Not. R. Astron. Soc. 274:1019, 1995), the period decreasing rate predicted by angular momentum loss (AML) caused by magnetic wind is computed for each system. A comparison between the observed and calculated period decreasing rates shows that the former values are obviously larger than the latter by 1–3 orders of magnitude. It suggests that the magnetic wind is not likely the determinant mechanism driving the AML in these systems. Finally, the orbital angular momentum (AM) and the rate of AML, $\dot{J}$ , are computed for each system. It shows that the AM have a similar change with the orbital period like dP/dt does, but $\log\dot{J}/J$ presents no strict changing with the kinematical ages.     

The Poynting-Robertson effect in the Newtonian potential with a Yukawa correction

We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.     

Changes in binding energy of interacting galaxies

It is shown that the fractional increase in binding energy of a galaxy in a fast collision with another galaxy of the same size can be well represented by the formula $$\xi _2 = 3({G \mathord{\left/ {\vphantom {G {M_2 \bar R}}} \right. \kern-\nulldelimiterspace} {M_2 \bar R}}) ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {V_p }}} \right. \kern-\nulldelimiterspace} {V_p }})^2 e^{ - p/\bar R} = \xi _1 ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {M_2 }}} \right. \kern-\nulldelimiterspace} {M_2 }})^3 ,$$ whereM ₁,M ₂ are the masses of the perturber and the perturbed galaxy, respectively,V _p is the relative velocity of the perturber at minimum separationp, and \(\bar R\) is the dynamical radius of either galaxy.     

UGC 5119: A galaxy with a stellar polar ring?

We present our photometric BV R_c observations of UGC 5119, a candidate polar-ring galaxy. We have determined its absolute magnitude, \(M_{0,B} = - 20\mathop m\limits_. 3\), and total color indices, \((B - V)_t^0 = + 0\mathop m\limits_. 73 \pm 0\mathop m\limits_. 10\) and \((V - R_c )_t^0 = + 0\mathop m\limits_. 54 \pm 0\mathop m\limits_. 10\). A Fourier analysis of the shape of its isophotes shows that UGC 5119 is most likely an elliptical galaxy with a disk component in the central part and a “boxy” feature on the periphery. At distances larger than 8″, the galaxy exhibits a turn of its major axis and a change in the phase of the fourth harmonic. Assuming the position angle of the major axis to be constant, a stellar ringlike structure is distinguished in the galaxy. The age of the ring stars is the same as that of the stars in the host galaxy. The distinguished ringlike structure cannot be attributed to typical polar rings rich in gas and in young stars.     

Secular dynamics of a lunar orbiter: a global exploration using Prony’s frequency analysis

We study the secular dynamics of lunar orbiters, in the framework of high-degree gravity models. To achieve a global view of the dynamics, we apply a frequency analysis (FA) technique which is based on Prony’s method. This allows for an extensive exploration of the eccentricity ( $e$ )—inclination ( $i$ ) space, based on short-term integrations ( $\sim $ 8 months) over relatively high-resolution grids of initial conditions. Different gravity models are considered: 3rd, 7th and 10th degree in the spherical harmonics expansion, with the main perturbations from the Earth being added. Since the dynamics is mostly regular, each orbit is characterised by a few parameters, whose values are given by the spectral decomposition of the orbital elements time series. The resulting frequency and amplitude maps in ( $e_0,i_0$ ) are used to identify the dominant perturbations and deduce the “minimum complexity” model necessary to capture the essential features of the long-term dynamics. We find that the 7th degree zonal harmonic ( $J_7$ term) is of profound importance at low altitudes as, depending on the initial secular phases, it can lead to collision with the Moon’s surface within a few months. The 3rd-degree non-axisymmetric terms are enough to describe the deviations from the 1 degree-of-freedom zonal problem; their main effect is to modify the equilibrium value of the argument of periselenium, $\omega $ , with respect to the “frozen” solution ( $\omega =\pm 90^{\circ }, \forall \Omega $ , where $\Omega $ is the nodal longitude). Finally, we show that using FA on a fine grid of initial conditions, set around a suitably chosen ‘first guess’, one can compute an accurate approximation of the initial conditions of a periodic orbit.     

A simple procedure to extend the Gauss method of determining orbital parameters from three to N points

A simple procedure is developed to determine orbital elements of an object orbiting in a central force field which contribute more than three independent celestial positions. By manipulation of formal three point Gauss method of orbit determination, an initial set of heliocentric state vectors r _i and $\dot{\mathbf{r}}_{i}$ is calculated. Then using the fact that the object follows the path that keep the constants of motion unchanged, I derive conserved quantities by applying simple linear regression method on state vectors r _i and $\dot{\mathbf{r}}_{i}$ . The best orbital plane is fixed by applying an iterative procedure which minimize the variation in magnitude of angular momentum of the orbit. Same procedure is used to fix shape and orientation of the orbit in the plane by minimizing variation in total energy and Laplace Runge Lenz vector. The method is tested using simulated data for a hypothetical planet rotating around the sun.     

Analysis of the exoplanet containing system Kepler-13

We have applied the close binary system analysis program WinFitter, with its physically detailed fitting function, to an intensive study of the complex multiple system Kepler-13 using photometry data from all 13 short cadence quarters downloaded from the NASA Exoplanet Archive (NEA) (http://exoplanetarchive.ipac.caltech.edu). The data-point error of our normalized, phase-sequenced and binned (380 points per bin: 0.00025 phase interval) flux values, at 14 ppm, allows the model’s specification for the mean reference flux level of the system to a precision better than 1 ppm. Our photometrically derived values for the mass and radius of KOI13.01 are \(6.8\pm0.6~\mbox{M}_{\mathrm{J}}\) and \(1.44\pm0.04~\mbox{R}_{\mathrm{J}}\). The star has a radius of \(1.67\pm0.05~\mbox{R}_{\odot}\). Our modelling sets the mean of the orbital inclination \(i\) at \(94.35\pm0.14^{\circ}\), with the star’s mean precession angle \(\phi_{p}\)—\(49.1\pm5.0^{\circ}\) and obliquity \(\theta_{o}\)\(67.9 \pm 3.0^{\circ}\), though there are known ambiguities about the sense in which such angles are measured.Our findings did not confirm secular variation in the transit modelling parameters greater than their full correlated errors, as argued by previous authors, when each quarter’s data was best-fitted with a determinable parameter set without prejudice. However, if we accept that most of the parameters remain the same for each transit, then we could confirm a small but steady diminution in the cosine of the orbital inclination over the 17 quarter timespan. This is accompanied by a slight increase of the star’s precession angle (less negative), but with no significant change in the obliquity of its spin axis. There are suggestions of a history of strong dynamical interaction with a highly distorted planet rotating in a 3:2 resonance with its revolution, together with a tidal lag of \(\sim30~\mbox{deg}\). The mean precessional period is derived to be about 1000 y, but at the present time the motion of the star’s rotation axis appears to be supporting the gravitational torque, rather than providing the balance against it that would be expected over long periods of time.The planet has a small but detectable backwarming effect on the star, which helps to explain the difference in brightness just after transit and just before occultation eclipses. In assessing these findings it is recognized that sources of uncertainty remain, notably with possible inherent micropulsational effects, variations from other components of the multiple star, stellar activity, differential rotation and the neglect of higher order terms (than \(r_{1}^{5}\)) in the fitting function, where \(r_{1}\) is the ratio of the radius of the star to the mean orbital separation of planet and host star.     

The theory of the Trojan asteroids,part V

E.W. Brown conjectured (1911) that the family of the long-periodic orbits in the Troian case of the restricted problem of three bodies terminates in an asymptotic orbit passing through the Lagrangian point L₃ at t=±∞. In 1977 the author showed that such an orbit deviates from L₃ by the epicyclic term mg (±∞). It is shown here that $$g\left( { \pm \infty } \right) = 0,$$ so that the Brown conjecture regarding L₃ is false. Contrary to what Brown believed, there is an entire family ofhomoclinic orbits, doubly asymptotic to short-periodic orbits around L₃. In the complex z-plane of the Poincaré eccentric variables, the latter orbits are circles of radius mR, with R bounded away from zero. The kinematics of the homoclinic family is investigated here in some detail.     

Cosmic-ray antimatter: A primary origin hypothesis

We examine the possibility that the observed cosmic-ray protons are of primary extragalactic origin. The present \(\bar p\) data are consistent with a primary extragalactic component having \(\bar p\) /p?3.2±0.7 x 10^-4 independent of energy. Following the suggestion that most extragalactic cosmic rays are from active galaxies, we propose that most of the observed \(\bar p\) 's are alos from the same sites. This would imply the possibility of destroying the corresponding \(\bar \alpha \) 'sat the source, thus leading to a flux ratio \(\bar \alpha \) /α< \(\bar p\) /p. We further predict an estimate for \(\bar \alpha \) α～10^-5, within the range of future cosmic-ray detectors. the cosmological implications of this proposal are discussed.     

Sur un modele explicitant les differents types morphologiques des galaxies

In the now classical Lindblad-Lin density-wave theory, the linearization of the collisionless Boltzmann equation is made by assuming the potential functionU expressed in the formU=U ₀ + \(\tilde U\) +... WhereU ₀ is the background axisymmetric potential and \(\tilde U<< U_0 \) . Then the corresponding density distribution is \(\rho = \rho _0 + \tilde \rho (\tilde \rho<< \rho _0 )\) and the linearized equation connecting \(\tilde U\) and the component \(\tilde f\) of the distribution function is given by $$\frac{{\partial \tilde f}}{{\partial t}} + \upsilon \frac{{\partial \tilde f}}{{\partial x}} - \frac{{\partial U_0 }}{{\partial x}} \cdot \frac{{\partial \tilde f}}{{\partial \upsilon }} = \frac{{\partial \tilde U}}{{\partial x}}\frac{{\partial f_0 }}{{\partial \upsilon }}.$$ One looks for spiral self-consistent solutions which also satisfy Poisson's equation $$\nabla ^2 \tilde U = 4\pi G\tilde \rho = 4\pi G\int {\tilde f d\upsilon .} $$ Lin and Shu (1964) have shown that such solutions exist in special cases. In the present work, we adopt anopposite proceeding. Poisson's equation contains two unknown quantities \(\tilde U\) and \(\tilde \rho \) . It could be completelysolved if a second independent equation connecting \(\tilde U\) and \(\tilde \rho \) was known. Such an equation is hopelesslyobtained by direct observational means; the only way is to postulate it in a mathematical form. In a previouswork, Louise (1981) has shown that Poisson's equation accounted for distances of planets in the solar system(following to the Titius-Bode's law revised by Balsano and Hughes (1979)) if the following relation wasassumed $$\rho ^2 = k\frac{{\tilde U}}{{r^2 }} (k = cte).$$ We now postulate again this relation in order to solve Poisson's equation. Then, $$\nabla ^2 \tilde U - \frac{{\alpha ^2 }}{{r^2 }}\tilde U = 0, (\alpha ^2 = 4\pi Gk).$$ The solution is found in a classical way to be of the form $$\tilde U = cte J_v (pr)e^{ - pz} e^{jn\theta } $$ wheren = integer,p =cte andJ _v (pr) = Bessel function with indexv (v ² =n ² + α²). By use of the Hankel function instead ofJ _v (pr) for large values ofr, the spiral structure is found to be given by $$\tilde U = cte e^{ - pz} e^{j[\Phi _v (r) + n\theta ]} , \Phi _v (r) = pr - \pi /2(v + \tfrac{1}{2}).$$ For small values ofr, \(\tilde U\) = 0: the center of a galaxy is not affected by the density wave which is onlyresponsible of the spiral structure. For various values ofp,n andv, other forms of galaxies can be taken into account: Ring, barred and spiral-barred shapes etc. In order to generalize previous calculations, we further postulateρ ₀ =kU ₀/r ², leading to Poisson'sequation which accounts for the disc population $$\nabla ^2 U_0 - \frac{{\alpha ^2 }}{{r^2 }}U_0 = 0.$$ AsU ₀ is assumed axisymmetrical, the obvious solution is of the form $$U_0 = \frac{{cte}}{{r^v }}e^{ - pz} , \rho _0 = \frac{{cte}}{{r^{2 + v} }}e^{ - pz} .$$ Finally, Poisson's equation is completely solvable under the assumptionρ =k(U/r ². The general solution,valid for both disc and spiral arm populations, becomes $$U = cte e^{ - pz} \left\{ {r^{ - v} + } \right.\left. {cte e^{j[\Phi _v (r) + n\theta ]} } \right\},$$ The density distribution along the O _z axis is supported by Burstein's (1979) observations.     

On the period of the periodic orbits of the restricted three body problem

We will show that the period T of a closed orbit of the planar circular restricted three body problem (viewed on rotating coordinates) depends on the region it encloses. Roughly speaking, we show that, \(2 T=k\pi +\int _\Omega g\) where k is an integer, \(\Omega \) is the region enclosed by the periodic orbit and \(g:{\mathbb {R}}^2\rightarrow {\mathbb {R}}\) is a function that only depends on the constant C known as the Jacobian constant; it does not depend on \(\Omega \). This theorem has a Keplerian flavor in the sense that it relates the period with the space “swept” by the orbit. As an application we prove that there is a neighborhood around \(L_4\) such that every periodic solution contained in this neighborhood must move clockwise. The same result holds true for \(L_5\).     

Orbites periodiques de satellites

In a previous paper (Stellmacher, 1981, hereafter mentioned as Paper I), we have given an algorithm for the construction of periodic orbits in a rotating frame, for satellites around an oblate planet. In the present paper, we apply this theory to the Mimas-Tethys case; we obtain the following results:

1. Without resonance, it is possible to find a rotating system in which the solution is a periodic one. The angular velocity of this rotating frame is calculated as function of the masses of the two satellites.
2. Including the resonant terms and assuming an exact commensurability of the implied frequencies, we demonstrate that the condition for periodic solutions in the rotating system as defined in (a) is: the initial position of the satellites at conjunction lies on an axis defined by (Ω₁+Ω₂)/2 or (Ω₁+Ω₂)/2 + π/2;Ω₁ and Ω₂ are the longitudes of the ascending nodes of the satellite's orbits. The solution still is a periodic one, thus all the conjunction occur in either axis.
3. In the Mimas Tethys case there is only approximately commensurability between these frequencies. The two satellites are considered as oscillators whose amplitudes and phases are functions of time. The equation of the libration can be established; we find the usual form, but for each satellite the generating solution is a periodic solution (as defined in Paper I), but not a Keplerian one. It follows a determination of the masses which slightly differs from that given by Kozai (1957), when the same values of the observed quantities are used for calculations.
4. The equation of the libration is: $$\ddot z + n_1^2 h^2 \sin z + n_1 q\dot z\sin z = 0$$

     

Comparative analysis of photometric variability of TT ARI in the years 1994–1995 and 2001, 2004

We present the results of photometric observations of a bright cataclysmic variable TT Ari with an orbital period of 0.13755 days. CCD observations were carried out with the Russian-Turkish RTT 150 telescope in 2001 and 2004 (13 nights). Multi-color photoelectric observations of the system were obtained with the Zeiss 600 telescope of SAO RAS in 1994–1995 (6 nights). In 1994–1995, the photometric period of the system was smaller than the orbital one (0 _. ^d 132 and 0 _. ^d 134), whereas it exceeded the latter (0 _. ^d 150 and 0 _. ^d 148) in 2001, 2004. An additional period exceeding the orbital one (0 _. ^d 144) is detected in 1995 modulations. We interpret it as indicating the elliptic disc precession in the direction of the orbital motion. In 1994, the variability in colors shows periods close to the orbital one (0 _. ^d 136, b-v), as well as to the period indicating the elliptic disk precession (0 _. ^d 146, w-b). We confirm that during the epochs characterized by photometric periods shorter than the orbital one, the quasi-periodic variability of TT Ari at time scales about 20 min is stronger than during epochs with long photometric periods. In general, the variability of the system can be described as a “red” noise with increased amplitudes of modulations at characteristic time scales of 10–40 min.