An intriguing question in the context of dynamics arises: Could a moon possess a moon itself? Such a configuration does not exist in the Solar System, although this may be possible in theory. Kollmeier et al. (2019) determined the critical size of a satellite necessary to host a long-lived sub-satellite, or submoon. However, the orbital constraints for these submoons to exist are still undetermined. Domingos et al. (2006) indicated that moons are stable out to a fraction of the host planet Hill radius $R_{H,p}$, which in turn depends on the eccentricity of its host's orbit. Motivated by this, we simulate a system of exomoons and submoons for $10^5$ planetary orbits, while considering many initial orbital phases to obtain the critical semimajor axis in terms of $R_{H,p}$ or the hosts satellite's Hill radius $R_{H,sat}$, respectively. We find that, assuming circular coplanar orbits, the stability limit for exomoons is 0.40 $R_{H,p}$ and for a submoon is 0.33 $R_{H,sat}$. Additionally, we discuss the observational feasibility of detecting these sub-satellites through photometric, radial velocity, or direct imaging observations using the Neptunes-sized exomoon candidates Kepler 1625b-I (Teachey et al. 2018) and identify how stability can shape the identification of future candidates.

We investigate the Interface Region Imaging Spectrograph (IRIS) observations of the quiet-Sun (QS) to understand the propagation of acoustic waves in transition region (TR) from photosphere. We selected a few IRIS spectral lines, which include the photospheric (Mn~{\sc i} 2801.25~{\AA}), chromospheric (Mg~{\sc ii} k 2796.35~{\AA}) and TR (C~{\sc ii} 1334.53~{\AA}), to investigate the acoustic wave propagation.The wavelet cross-spectrum reveals significant coherence (about 70\% locations) between photosphere and chromosphere. Few minutes oscillations (i.e., period range from 1.6 to 4.0 minutes) successfully propagate into chromosphere from photosphere, which is confirmed by dominance of positive phase lags. However, in higher period regime (i.e., greater than $\approx$ 4.5 minutes), the downward propagation dominates is evident by negative phase lags. The broad spectrum of waves (i.e., 2.5-6.0 minutes) propagates freely upwards from chromosphere to TR. We find that only about 45\% locations (out of 70\%) show correlation between chromosphere and TR. Our results indicate that roots of 3 minutes oscillations observed within chromosphere/TR are located in photosphere. Observations also demonstrate that 5 minute oscillations propagate downward from chromosphere. \textbf{However, some locations within QS also show successful propagation of 5 minute oscillations as revealed by positive phase lags, which might be the result of magnetic field}. In addition, our results clearly show that a significant power, within period ranging from 2.5 to 6.0 minutes, of solar chromosphere is freely transmitted into TR triggering atmospheric oscillations. Theoretical implications of our observational results are discussed.

Stability of planetary orbits around the GJ 832 star system, which contains inner (GJ 832c) and outer (GJ 832b) planets, is investigated numerically and a detailed phase-space analysis is performed. A special emphasis is given to the existence of stable orbits for a planet less than 15M$_\oplus$ which is injected between the inner and outer planets. Thus, numerical simulations are performed for three and four bodies in elliptical orbits (or circular for special cases) by using a large number of initial conditions that cover the selected phase-spaces of the planet's orbital parameters. The results presented in the phase-space maps for GJ 832c indicate the least deviation of eccentricity from its nominal value, which is then used to determine its inclination regime relative to the star-outer planet plane. Also, the injected planet displays stable orbital configurations for at least one billion years. Then, the radial velocity curves based on the signature from the Keplerian motion are generated for the injected planets with masses 1M$_\oplus$ to 15M$_\oplus$ in order to estimate their semimajor axes and mass-limit. The synthetic RV signal suggests that an additional planet of mass $\le$ 15M$_\oplus$ with dynamically stable configuration may be residing between 0.25 - 2.0 AU from the star. We have provided an estimated number of RV observations for the additional planet that is required for further observational verification.

We investigate numerically Alfv\'en waves propagating along an axisymmetric and non-isothermal solar flux tube embedded in the solar atmosphere. The tube magnetic field is current-free and diverges with height, and the waves are excited by a periodic driver along the tube magnetic field lines. The main results are that the two wave variables, the velocity and magnetic field perturbations in the azimuthal direction, behave differently as a result of gradients of physical parameters along the tube. To explain these differences in the wave behavior, the time evolution of the wave variables and the resulting cutoff period for each wave variable are calculated, and used to determine regions in the solar chromosphere where strong wave reflection may occur.

Recently Wi\'sniewska et al. demonstrated observationally how the acoustic cutoff frequency varies with height in the solar atmosphere including the upper photosphere and the lower and middle chromosphere, and showed that the observational results cannot be accounted for by the existing theoretical formulas for the acoustic cutoff. In order to reproduce the observed variation of the cutoff with atmospheric height, numerical simulations of impulsively generated acoustic waves in the solar atmosphere are performed, and the spectral analysis of temporal wave profiles is used to compute numerically changes of the acoustic cutoff with height. Comparison of the numerical results with the observational data shows good agreement, which clearly indicates that the obtained results may be used to determine the structure of the background solar atmosphere.

The idea of single exomoon detection due to the radio emissions caused by its interaction with the host exoplanet is extended to multiple-exomoon systems. The characteristic radio emissions are made possible in part by plasma from the exomoons own ionosphere (Noyola et. al. 2014). In this work, it is demonstrated that neighboring exomoons and the exoplanetary magnetosphere could also provide enough plasma to generate a detectable signal. In particular, the plasma-torus-sharing phenomenon is found to be particularly well suited to facilitate the radio detection of plasma-deficient exomoons. The efficiency of this process is evaluated, and the predicted power and frequency of the resulting radio signals are presented.

The three-body problem, which describes three masses interacting through Newtonian gravity without any restrictions imposed on the initial positions and velocities of these masses, has attracted the attention of many scientists for more than 300 years. In this paper, we present a review of the three-body problem in the context of both historical and modern developments. We describe the general and restricted (circular and elliptic) three-body problems, different analytical and numerical methods of finding solutions, methods for performing stability analysis, search for periodic orbits and resonances, and application of the results to some interesting astronomical and space dynamical settings. We also provide a brief presentation of the general and restricted relativistic three-body problems, and discuss their astronomical applications.

The rise of oxygen in Earth's atmosphere that occurred 2.4 to 2.2 billion years ago is known as the Earth's Great Oxidation, and its impact on the development of life on Earth has been profound. Thereafter, the increase in Earth's oxygen level persisted, though at a more gradual pace. The proposed underlying mathematical models for these processes are based on physical parameters whose values are currently not well-established owing to uncertainties in geological and biological data. In this paper, a previously developed model of Earth's atmosphere is modified by adding different strengths of noise to account for the parameters' uncertainties. The effects of the noise on the time variations of oxygen, carbon and methane for the early Earth are investigated by using fractal and multifractal analysis. We show that the time variations following the Great Oxidation cannot properly be described by a single fractal dimension because they exhibit multifractal characteristics. The obtained results demonstrate that the time series as obtained exhibit multifractality caused by long-range time correlations.

We perform numerical simulations of impulsively generated Alfv\'en waves in an isolated photospheric flux tube, and explore the propagation of these waves along such magnetic structure that extends from the photosphere, where these waves are triggered, to the solar corona, and analyze resulting magnetic shells. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical model and specifying the curved magnetic field lines that constitute the magnetic flux tube which is rooted in the solar photosphere. The evolution of the solar atmosphere is described by 3D, ideal magnetohydrodynamic equations that are numerically solved by the FLASH code. Our numerical simulations reveal, based on the physical properties of the multi-shell magnetic twisters and the amount of energy and momentum associated with them, that these multi-shell magnetic twisters may be responsible for the observed heating of the lower solar corona and for the formation of solar wind. Moreover, it is likely that the existence of these twisters can be verified by high-resolution observations.

Aims: Propagation and energy transfer of torsional Alfv\'en waves in solar magnetic flux tubes of axial symmetry is studied. Methods: An analytical model of a solar magnetic flux tube of axial symmetry is developed by specifying a magnetic flux and deriving general analytical formulae for the equilibrium mass density and a gas pressure. The main advantage of this model is that it can be easily adopted to any axisymmetric magnetic structure. The model is used to simulate numerically the propagation of nonlinear Alfv\'en waves in such 2D flux tubes of axial symmetry embedded in the solar atmosphere. The waves are excited by a localized pulse in the azimuthal component of velocity and launched at the top of the solar photosphere, and they propagate through the solar chromosphere, transition region, and into the solar corona. Results: The results of our numerical simulations reveal a complex scenario of twisted magnetic field lines and flows associated with torsional Alfv\'en waves as well as energy transfer to the magnetoacoustic waves that are triggered by the Alfv\'en waves and are akin to the vertical jet flows. Alfv\'en waves experience about 5 % amplitude reflection at the transition region. Magnetic (velocity) field perturbations experience attenuation (growth) with height is agreement with analytical findings. Kinetic energy of magnetoacoustic waves consists of 25 % of the total energy of Alfv\'en waves. The energy transfer may lead to localized mass transport in the form of vertical jets, as well as to localized heating as slow magnetoacoustic waves are prone to dissipation in the inner corona.

We perform numerical simulations of impulsively generated Alfv\'en waves in an isolated solar arcade, which is gravitationally stratified and magnetically confined. We study numerically the propagation of Alfv\'en waves along such magnetic structure that extends from the lower chromosphere, where the waves are generated, to the solar corona, and analyze influence of the arcade size and width of the initial pulses on the wave propagation and reflection. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical VAL-C model and specifying the curved magnetic field lines that constitute the asymmetric magnetic arcade. The propagation and reflection of Alfv\'en waves in this arcade is described by 2.5D magnetohydrodynamic equations that are numerically solved by the FLASH code. Our numerical simulations reveal that the Alfv\'en wave amplitude decreases as a result of a partial reflection of Alfv\'en waves in the solar transition region, and that the waves which are not reflected leak through the transition region and reach the solar corona. We also find the decrement of the attenuation time of Alfv\'en waves for wider initial pulses. Moreover, our results show that the propagation of Alfv\'en waves in the arcade is affected by spatial dependence of the Alfv\'en speed, which leads to phase-mixing that is stronger for more curved and larger magnetic arcades. We discuss processes that affect the Alfv\'en wave propagation in an asymmetric solar arcade and conclude that besides phase-mixing in the magnetic field configuration, plasma properties of the arcade and size of the initial pulse as well as structure of the solar transition region all play a vital role in the Alfv\'en wave propagation.

In the Jupiter-Io system, the moon's motion produces currents along the field lines that connect it to Jupiter's polar regions. The currents generate, and modulate radio emissions along their paths via the electron-cyclotron maser instability. Based on this process, we suggest that such modulation of planetary radio emissions may reveal the presence of exomoons around giant planets in exoplanetary systems. A model explaining the modulation mechanism in the Jupiter-Io system is extrapolated, and used to define criteria for exomoon detectability. A cautiously optimistic scenario of possible detection of such exomoons around Epsilon Eridani b, and Gliese 876 b is provided.

We present results of 3-D numerical simulations of a fast magnetic twister excited above a foot-point of the potential solar coronal arcade that is embedded in the solar atmosphere with the initial VAL-IIIC temperature profile, which is smoothly extended into the solar corona. With the use of the FLASH code, we solve 3-D ideal magnetohydrodynamic equations by specifying a twist in the azimuthal component of magnetic field in the solar chromosphere. The imposed perturbation generates torsional Alfv\'en waves as well as plasma swirls that reach the other foot-point of the arcade and partially reflect back from the transition region. The two vortex channels are evident in the generated twisted flux-tube with a fragmentation near its apex that results from the initial twist as well as from the morphology of the tube. The numerical results are compared to observational data of plasma motions in a solar prominence. The comparison shows that the numerical results and the data qualitatively agree even though the observed plasma motions occur over comparatively large spatio-temporal scales in the prominence.

We present simulation results of the impulsively generated linear and non-linear Alfv\'en waves in the weakly curved coronal magnetic flux-tubes (coronal funnels) and discuss their implications for the coronal heating and solar wind acceleration. We solve numerically the time-dependent magnetohydrodynamic (MHD) equations to obtain the temporal signatures of the small (linear) and large-amplitude (non-linear) Alfv\'en waves in the model atmosphere of expanding open magnetic field configuration (e.g., coronal funnels) by considering a realistic temperature distribution. We compute the maximum transversal velocity of both linear and non-linear Alfv\'en waves at different heights in the coronal funnel, and study their response in the solar corona during the time of their propagation. We infer that the pulse-driven non-linear Alfv\'en waves may carry sufficient wave energy fluxes to heat the coronal funnels and also to power the solar wind that originates in these funnels. Our study of linear Alfv\'en waves show that they can contribute only to the plasma dynamics and heating of the funnel-like magnetic flux-tubes associated with the polar coronal holes.

The propagation of linear transverse waves along a thin isothermal magnetic flux tube is affected by a global cutoff frequency that separates propagating and non-propagating waves. In this paper, wave propagation along a thin but non-isothermal flux tube is considered and a local cutoff frequency is derived. The effects of different temperature profiles on this local cutoff frequency are studied by considering different power-law temperature distributions as well as the semi-empirical VAL C model of the solar atmosphere. The results show that the conditions for wave propagation strongly depend on the temperature gradients. Moreover, the local cutoff frequency calculated for the VAL C model gives constraints on the range of wave frequencies that are propagating in different parts of the solar atmosphere. These theoretically predicted constraints are compared to observational data and are used to discuss the role played by transverse tube waves in the atmospheric heating and dynamics, and in the excitation of solar atmospheric oscillations.

We study the impulsively generated non-linear Alfv\'en waves in the solar atmosphere, and describe their most likely role in the observed non-thermal broadening of some spectral lines in solar coronal holes. We solve numerically the time-dependent magnetohydrodynamic equations to find temporal signatures of large-amplitude Alfv\'en waves in the model atmosphere of open and expanding magnetic field configuration, with a realistic temperature distribution. We calculate the temporally and spatially averaged, instantaneous transversal velocity of non-linear Alfv\'en waves at different heights of the model atmosphere, and estimate its contribution to the unresolved non-thermal motions caused by the waves. We find that the pulse-driven nonlinear Alfv\'en waves with the amplitude $A_{\rm v}$=50 km s$^{-1}$ are the most likely candidates for the non-thermal broadening of Si VIII $\lambda$1445.75 \AA\ line profiles in the polar coronal hole as reported by Banerjee et al. (1998). We also demonstrate that the Alfv\'en waves driven by comparatively smaller velocity pulse with its amplitude $A_{\rm v}$=25 km s$^{-1}$ may contribute to the spectral line width of the same line at various heights in coronal hole without any significant broadening. The main conclusion of this paper is that non-linear Alfv\'en waves excited impulsively in the lower solar atmosphere are responsible for the observed spectral line broadening in polar coronal holes. This is an important result as it allows us to conclude that such large amplitude and pulse-driven Alfv\'en waves do indeed exist in solar coronal holes. The existence of these waves and their undamped growth may impart the required momentum to accelerate the solar wind.

Our aim is to identify and classify mean-motion resonances (MMRs) for the coplanar circular restricted three-body problem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what kind of resonances are present, and the inspection of the orbital elements to classify the periodicity. We show that the 2:1 resonance occurs the most frequently. Among other resonances, the 3:1 resonance is the second most common, and furthermore both 3:2 and 5:3 resonances occur more often than the 4:1 resonance. Moreover, the resonances in the coplanar CR3BP are classified based on the behaviour of the orbits. We show that orbital stability is ensured for high values of resonance (i.e., high ratios) where only a single resonance is present. The resonances attained are consistent with the previously established resonances for the Solar System, i.e., specifically, in regards to the asteroid belt. Previous work employed digital filtering and Lyapunov characteristic exponents to determine stochasticity of the eccentricity, which is found to be consistent with our usage of Lyapunov exponents as an alternate approach based on varying the mass ratio instead of the eccentricity. Our results are expected to be of principal interest to future studies, including augmentations to observed or proposed resonances, of extra-solar planets in binary stellar systems.

We demonstrate that habitable Earth-mass planets and moons can exist in the Kepler-16 system, known to host a Saturn-mass planet around a stellar binary, by investigating their orbital stability in the standard and extended habitable zone (HZ). We find that Earth-mass planets in satellite-like (S-type)orbits are possible within the standard HZ in direct vicinity of Kepler-16b, thus constituting habitable exomoons. However, Earth-mass planets cannot exist in planetary-like (P-type) orbits around the two stellar components within the standard HZ. Yet, P-type Earth-mass planets can exist superior to the Saturnian planet in the extended HZ pertaining to considerably enhanced back-warming in the planetary atmosphere if facilitated. We briefly discuss the potential detectability of such habitable Earth-mass moons and planets positioned in satellite and planetary orbits, respectively. The range of inferior and superior P-type orbits in the HZ is between 0.657 to 0.71 AU and 0.95 to 1.02 AU, respectively.

We provide a detailed theoretical study aimed at the observational finding about the nu Octantis binary system that indicates the possible existence of a Jupiter-type planet in this system. If a prograde planetary orbit is assumed, it has earlier been argued that the planet, if existing, should be located outside the zone of orbital stability. However, a previous study by Eberle & Cuntz (2010) [ApJ 721, L168] concludes that the planet is most likely stable if assumed to be in a retrograde orbit with respect to the secondary system component. In the present work, we significantly augment this study by taking into account the observationally deduced uncertainty ranges of the orbital parameters for the stellar components and the suggested planet. Furthermore, our study employs additional mathematical methods, which include monitoring the Jacobi constant, the zero velocity function, and the maximum Lyapunov exponent. We again find that the suggested planet is indeed possible if assumed to be in a retrograde orbit, but it is virtually impossible if assumed in a prograde orbit. Its existence is found to be consistent with the deduced system parameters of the binary components and of the suggested planet, including the associated uncertainty bars given by observations.

We establish a criterion for the stability of planetary orbits in stellar binary systems by using Lyapunov exponents and power spectra for the special case of the circular restricted 3-body problem (CR3BP). The centerpiece of our method is the concept of Lyapunov exponents, which are incorporated into the analysis of orbital stability by integrating the Jacobian of the CR3BP and orthogonalizing the tangent vectors via a well-established algorithm originally developed by Wolf et al. The criterion for orbital stability based on the Lyapunov exponents is independently verified by using power spectra. The obtained results are compared to results presented in the two previous papers of this series. It is shown that the maximum Lyapunov exponent can be used as an indicator for chaotic behaviour of planetary orbits, which is consistent with previous applications of this method, particularly studies for the Solar System. The chaotic behaviour corresponds to either orbital stability or instability, and it depends solely on the mass ratio of the binary components and the initial distance ratio of the planet relative to the stellar separation distance. Our theoretical results allow us to link the study of planetary orbital stability to chaos theory noting that there is a large array of literature on the properties and significance of Lyapunov exponents. Although our results are given for the special case of the CR3BP, we expect that it may be possible to augment the proposed Lyapunov exponent criterion to studies of planets in generalized stellar binary systems, which is strongly motivated by existing observational results as well as results expected from ongoing and future planet search missions.

We investigate the possibility of habitable Trojan planets in the HD 23079 star-planet system. This system consists of a solar-type star and a Jupiter-type planet, which orbits the star near the outer edge of the stellar habitable zone in an orbit of low eccentricity. We find that in agreement with previous studies Earth-mass habitable Trojan planets are possible in this system, although the success of staying within the zone of habitability is significantly affected by the orbital parameters of the giant planet and by the initial condition of the theoretical Earth-mass planet. In one of our simulations, the Earth-mass planet is captured by the giant planet and thus becomes a habitable moon.

About half of all known stellar systems with Sun-like stars consist of two or more stars, significantly affecting the orbital stability of any planet in these systems. This observational evidence has prompted a large array of theoretical research, including the derivation of mathematically stringent criteria for the orbital stability of planets in stellar binary systems, valid for the "coplanar circular restricted three-body problem". In the following, we use these criteria to explore the validity of results from previous theoretical studies.

About half of all known stellar systems with Sun-like stars consist of two or more stars, significantly affecting the orbital stability of any planet in these systems. Here we study the onset of instability for an Earth-type planet that is part of a binary system. Our investigation makes use of previous analytical work allowing to describe the permissible region of planetary motion. This allows us to establish a criterion for the orbital stability of planets that may be useful in the context of future observational and theoretical studies.

Propagation of torsional Alfv\'en waves along magnetic flux tubes has been extensively studied for many years but no conclusive results regarding the existence of a cutoff frequency for these waves have been obtained. The main purpose of this paper is to derive new wave equations that describe the propagation of linear torsional Alfv\'en waves along thin and isothermal magnetic flux tubes, and use these wave equations to demonstrate that the torsional wave propagation is not affected by any cutoff frequency. It is also shown that this cutoff-free propagation is independent of different choices of the coordinate systems and wave variables adopted in the previous studies. A brief discussion of implications of this cutoff-free propagation of torsional tube waves on theories of wave heating of the solar and stellar atmospheres is also given.

The vector Burgers equation is extended to include pressure gradients and gravity. It is shown that within the framework of the Cole-Hopf transformation there are no physical solutions to this problem. This result is important because it clearly demonstrates that any extension of Burgers equation to more interesting physical situations is strongly limited.