Search

In this work, we investigate the near 2:1 mean motion resonance (MMR) capture scenario and tidal evolution for the hot-Neptune system HD 106315 by theoretical analysis and numerical simulations. To make a deeper understanding of migration for this system, we extensively explore the initial semi-major axis (SMA), eccentricity, as well as the attenuation coefficient of the eccentricity of HD 106315c that affect the orbital configuration of the system. The numerical simulations show that when the initial orbital elements of two planets are a b ∼ 0.4 au, a c ∼ 0.8 au, e b and e c lower than 0.03, HD 106315b and HD 106315c can migrate to currently observed location, and form near 2:1 mean motion resonance in about 65000 years. This scenario takes into account the gravitational effect of the central star and the viscosity of disk. In addition, the tidal dissipation that originates from the host star, can influence the MMR configuration for planets. From the point of theoretical analysis, the results show that for the tidal dissipation factor Q = 100 , the ratio of orbital period of two planets will change over the evolution. This may be indicative of the resonant departure for two planets in the system. However, according to numerical simulations, the tidal dissipation factor Q may be larger than 10 3 for the two planets, therefore we may safely conclude that tidal effect cannot play a major role in driving two Neptunian planets of the HD 106315 system move out of 2:1 mean motion resonance within the residual lifetime of the star.

The highly-substructured outskirts of the Magellanic Clouds provide ideal locations for studying the complex interaction history between both Clouds and the Milky Way (MW). In this paper, we investigate the origin of a >20$^\circ$ long arm-like feature in the northern outskirts of the Large Magellanic Cloud (LMC) using data from the Magellanic Edges Survey (MagES) and Gaia EDR3. We find that the arm has a similar geometry and metallicity to the nearby outer LMC disk, indicating that it is comprised of perturbed disk material. Whilst the azimuthal velocity and velocity dispersions along the arm are consistent with those in the outer LMC, the in-plane radial velocity and out-of-plane vertical velocity are significantly perturbed from equilibrium disk kinematics. We compare these observations to a new suite of dynamical models of the Magellanic/MW system, which describe the LMC as a collection of tracer particles within a rigid potential, and the SMC as a rigid Hernquist potential. Our models indicate the tidal force of the MW during the LMC's infall is likely responsible for the observed increasing out-of-plane velocity along the arm. Our models also suggest close LMC/SMC interactions within the past Gyr, particularly the SMC's pericentric passage ~150 Myr ago and a possible SMC crossing of the LMC disk plane ~400 Myr ago, likely do not perturb stars that today comprise the arm. Historical interactions with the SMC prior to ~1 Gyr ago may be required to explain some of the observed kinematic properties of the arm, in particular its strongly negative in-plane radial velocity., Accepted by MNRAS

A galaxy moving through a background of dark matter particles induces an overdensity of these particles or a wake behind it. The back reaction of this wake on the galaxy is a force field that can be decomposed into an effective deceleration (called dynamical friction) and a tidal field. In this paper, we determine the tidal forces, thus generated on the galaxy, and the resulting observables, which are shown to be warps, lopsidedness, and/or kinematic-photometric position angle misalignments. We estimate the magnitude of the tidal-like effects needed to reproduce the observed warp and lopsidedness on the isolated galaxy IC 2487. Within a realistic range of dark matter distribution properties, the observed, warped, and lopsided kinematical properties of IC 2487 are possible to reproduce (the background medium of dark matter particles has a velocity dispersion of $\lesssim 80\, {\rm km\, s^{-1}}$ and the density $10^4{\!-\!}10^5\, {\rm M_\odot \, kpc^{-3}}$, more likely at the lower end). We conclude that the proposed mechanism can generate warps, lopsidedness, and misalignments observed in isolated galaxies or galaxies in loose groups. The method can be used also to constrain dark matter spatial and velocity distribution properties.

We use the astrometric and photometric data from Gaia Data Release 2 and line-of-sight velocities from various other surveys to study the 3d structure and kinematics of the Sagittarius dwarf galaxy. The combination of photometric and astrometric data makes it possible to obtain a very clean separation of Sgr member stars from the Milky Way foreground; our final catalogue contains ~2.6e5 candidate members with magnitudes G, matches published version

Large dips in the brightness for a number of stars have been observed, for which the tentative explanation is occultation of the star by a transiting circumplanetary disk or ring system. In order for the circumplanetary disk/rings to block the host star's light, the disk must be tilted out of the planet's orbital plane, which poses stability problems due to the radial extent of the disk required to explain the brightness dip durations. This work uses N-body integrations to study the structure and stability of circumplanetary disk/ring systems tilted out of the planet's orbital plane by the spinning planet's mass quadrupole. Simulating the disk as a collection of test particles with orbits initialized near the Laplace surface (equilibrium between tidal force from host star and force from planet's mass quadrupole), we find that many extended, inclined circumplanetary disks remain stable over the duration of the integrations (~3-16 Myr). Two dynamical resonances/instabilities excite the particle eccentricities and inclinations: the Lidov-Kozai effect which occurs in the disk's outer regions, and ivection resonance which occurs in the disk's inner regions. Our work places constraints on the maximum radial extent of inclined circumplanetary disk/ring systems, and shows that gaps present in circumplanetary disks do not necessarily imply the presence of exomoons., Comment: 14 pages, 11 figures. Accepted to MNRAS (2020 July 20)

A star coming too close to a supermassive black hole gets disrupted by the tidal force of the compact object in a tidal disruption event, or TDE. Following this encounter, the debris evolves into an elongated stream, half of which coming back to pericentre. Relativistic apsidal precession then leads to a self-crossing shock that initiates the formation of an accretion disc. We perform the first simulation of this process considering a parabolic encounter with a supermassive black hole, which has so far eluded investigations for computational reasons. This numerical issue is alleviated by using as initial conditions the outflow launched by the self-crossing shock according the local simulation of Lu & Bonnerot (2020). We find that the gas leaving the intersection point experiences numerous secondary shocks that result in the rapid formation of a thick and marginally bound disc. The mass distribution features two overdensities identified as spiral shocks that drive slow gas inflow along the mid-plane. Inward motion primarily takes place along the funnels of the newly formed torus, from which a fraction of the matter can get accreted. Further out, the gas moves outward forming an extended envelope completely surrounding the accretion flow. Secondary shocks heat the debris at a rate of a few times $10^{44} \, \rm erg\, s^{-1}$ with a large fraction likely participating to the bolometric luminosity. These results pave the way towards a complete understanding of the early radiation from TDEs that progressively becomes accessible from observations.

The variation of the orbital periods found in high-mass X-ray binaries occurring over time can be explained by two models. According to the first model, orbital decay occurs due to dynamic tidal forces. The second model shows that the change occurs due to the inertia moments of a massive star during its evolution. These models are applied to LMC X-4 and Cen X-3.

In this paper, we discuss the possible disruption of a prolate ellipsoidal small body by the Earth’s tidal forces during a close flyby. An exact expression for the stress tensor under the action of gravity, rotation, and tidal forces is obtained using the analytical solution of the elastic problem. Constraints for the critical distance are derived taking into account the size, shape, density, mechanical properties (material strength and Poisson ratio), rotation, and two variants of the body orientation. Asteroid 433 Eros is considered as a case study. The results are important for methods for assessing the asteroid hazard of near-Earth objects (NEO) and the evolution of terrestrial planets.

Its extraordinary level of geologic activity, its potential for habitability, and the prospects of returning samples from its plume of erupting water ice make Saturn's small (∼500 km diameter) moon Enceladus a high priority target for future exploration and a key to our developing understanding of icy ocean worlds. The structure of its outer ice shell is particularly important as it relates to the global heat budget, the global-scale response to tidal forces, and the nature of the ongoing eruptions. It is also diagnostic of how and where heat is dissipated internally. Here, using the most recent shape model and a new approach to modeling isostasy, we obtain a shell structure that simultaneously accommodates the shape, gravity, and libration observations and suggests that tidal dissipation near the base of the ice shell is likely an important mode of internal heating. The implied conductive heat loss is greater than the heat loss associated with the eruptions but is nevertheless compatible with the condition of steady state.

We take a closer and new look at the effects of tidal forces on the free fall of a quantum particle inside a spherically symmetric gravitational field. We derive the corresponding Schr\"odinger equation for the particle by starting from the fully relativistic Klein-Gordon equation in order (i) to briefly discuss the issue of the equivalence principle and (ii) to be able to compare the relativistic terms in the equation to the tidal-force terms. To the second order of the nonrelativistic approximation, the resulting Schr\"odinger equation is that of a simple harmonic oscillator in the horizontal direction and that of an inverted harmonic oscillator in the vertical direction. Two methods are used for solving the equation in the vertical direction. The first method is based on a fixed boundary condition, and yields a discrete-energy spectrum with a wavefunction that is asymptotic to that of a particle in a linear gravitational field. The second method is based on time-varying boundary conditions and yields a quantized-energy spectrum that is decaying in time. Moving on to a freely-falling reference frame, we derive the corresponding time-dependent energy spectrum. The effects of tidal forces yield an expectation value for the Hamiltonian and a relative change in time of a wavepacket's width that are mass-independent. The equivalence principle, which we understand here as the empirical equivalence between gravitation and inertia, is discussed based on these various results. For completeness, we briefly discuss the consequences expected to be obtained for a Bose-Einstein condensate or a superfluid in free fall using the nonlinear Gross-Pitaevskii equation., Comment: 33 pages. References added and typos corrected. Accepted for publication in IJMPD

Apsidal motion is a gradual shift in the position of periastron. The impact of dynamic tides on apsidal motion has long been debated, because the contribution could not be quantified due to the lack of high quality observations. KIC 4544587 with tidally excited oscillations has been observed by \textit{Kepler} high-precision photometric data based on long time baseline and short-cadence schema. In this paper, we compute the rate of apsidal motion that arises from the dynamic tides as $19.05\pm 1.70$ mrad yr$^{-1}$ via tracking the orbital phase shifts of tidally excited oscillations. We also calculate the procession rate of the orbit due to the Newtonian and general relativistic contribution as $21.49 \pm 2.8$ and $2.4 \pm 0.06$ mrad yr$^{-1}$, respectively. The sum of these three factors is in excellent agreement with the total observational rate of apsidal motion $42.97 \pm 0.18$ mrad yr$^{-1}$ measured by eclipse timing variations. The tidal effect accounts for about 44\% of the overall observed apsidal motion and is comparable to that of the Newtonian term. Dynamic tides have a significant contribution to the apsidal motion. The analysis method mentioned in this paper presents an alternative approach to measuring the contribution of the dynamic tides quantitatively., 10 pages, 5 figures, accepted by ApJ

We present high spatial resolution (12pc) Atacama Large Millimeter/sub-millimeter Array CO(J=3-2) observations of the nearby lenticular galaxy NGC4429. We identify 217 giant molecular clouds within the 450pc radius molecular gas disc. The clouds generally have smaller sizes and masses but higher surface densities and observed linewidths than those of Milky Way disc clouds. An unusually steep size - line width relation and large cloud internal velocity gradients (0.05 - 0.91 km s^-1 pc^-1) and observed Virial parameters (alpha_obs,vir = 4.0) are found, that appear due to internal rotation driven by the background galactic gravitational potential. Removing this rotation, an internal Virial equilibrium appears to be established between the self-gravitational (Usg) and turbulent kinetic (Eturb) energies of each cloud, i.e. alpha_sg,vir=Usg/Eturb = 1.3. However, to properly account for both self and external gravity (shear and tidal forces), we formulate a modified Virial theorem and define an effective Virial parameter alpha_eff,vir = alpha_sg,vir + Usg/Eext (and associated effective velocity dispersion). The NGC4429 clouds then appear to be in a critical state in which the self-gravitational energy and the contribution of external gravity to the cloud's energy budget (Eext) are approximately equal, i.e. Eext/Usg~1. As such, alpha_eff,vir = 2.2 and most clouds are not virialised but remain marginally gravitationally bound. We show this is consistent with the clouds having sizes similar to their tidal radii and being generally radially elongated. External gravity is thus as important as self-gravity to regulate the clouds of NGC4429., Comment: Accepted by MNRAS

Despite numerous attempts, no exomoon has firmly been confirmed to date. New missions like CHEOPS aim to characterize previously detected exoplanets, and potentially to discover exomoons. In order to optimize search strategies, we need to determine those planets which are the most likely to host moons. We investigate the tidal evolution of hypothetical moon orbits in systems consisting of a star, one planet and one test moon. We study a few specific cases with ten billion years integration time where the evolution of moon orbits follows one of these three scenarios: (1) "locking", in which the moon has a stable orbit on a long time scale ($\gtrsim$ 10$^9$ years); (2) "escape scenario" where the moon leaves the planet's gravitational domain; and (3) "disruption scenario", in which the moon migrates inwards until it reaches the Roche lobe and becomes disrupted by strong tidal forces. Applying the model to real cases from an exoplanet catalogue, we study the long-term stability of moon orbits around known exoplanets. We calculate the survival rate which is the fraction of the investigated cases when the moon survived around the planet for the full integration time (which is the age of the star, or if not known, then the age of the Sun).The most important factor determining the long term survival of an exomoon is the orbital period of the planet. For the majority of the close-in planets (, 65 pages: 18 pages of text with figures, 47 pages of table in appendix. Accepted for publication in PASP

In nuclei of galaxies strong tidal forces can destroy stars passing within a critical distance from the central super-massive black hole (SMBH). Observational signatures of tidal disruption events (TDEs) depend on the environment around the SBMH horizon and the level of its accretion activity. Evidence for optical and UV spectral features has been reported in TDE flares; nevertheless, to test the effects of General Relativity in the immediate vicinity of SMBH, the relativistically broadened and skewed X-ray line would tell us significantly more useful information. This will require to proceed beyond inactive nuclei. To this end we consider a system where the material from a disrupted star forms a gaseous ring that circularises near the tidal radius around SMBH, it gradually spreads in radius by viscous processes, and resides embedded within a hot corona. In our test calculation the remnant trail is assumed to be fully circularized and embedded in a hot environment and illuminated by X-rays from a surrounding corona or a jet base of (mild) AGN activity. We show the expected effects on the observed profile and the centroid energy. In future the evolving spectral features can enhance the diagnostic capability and provide a novel way to reveal the parameters of TDEs in such sources; namely, the distance of the remnant gas from the SMBH, the radial extent of the gaseous trail, and the spin of the SMBH could be measured., Comment: 8 pages, 4 figures, to appear in Astronomische Nachrichten, proceedings of the conference IBWS 2019, Karlovy Vary

The recent hypothesis of a giant supercluster binary-like structure formed by the Corona Borealis and its close companion Abell 2142 (supercluster) belongs to a little known area of investigation as the dynamics of gravitationally interacting galaxy supercluster pairs. From the observational point of view this structure approximates the configuration of a binary-like system in linear orbit interconnected by a huge filamentary structure which, if confirmed, it would be the first case to date observed at z ≥ 0.07. Given the importance to disentangle this issue, a follow-up analysis has been performed on the region constrained by the common envelop of the two superclusters in order to search for new hints to confirm their mutual gravitational interaction. Observational signatures of that interaction have been found mapping the inner peculiar motions where the observed negative peculiar velocities measured within the A2142 (supercluster) region suggest a general matter flow toward the Corona Borealis supercluster. Besides, analyzing the effects on both superclusters due to the mutual impact of the external tidal forces, turns out that their inner dynamics remain unperturbed up to the turnaround radii. Outside, where the binding forces are overlapped by the tidal ones, the outskirts of both superclusters should be unstable and subject to fragmentation. Such a scenario indicates that both superclusters interact with comparable and reciprocal tidal perturbations leaving the whole system in a substantial dynamical equilibrium. The origin of such a dynamical dichotomy would be explained either by a much more massive Corona Borealis supercluster than that estimated in the present work or by a selection effect biasing the small sample of peculiar velocities due to the remoteness of the system worsened by the large uncertainty on their measurements.

In past studies, the reorientation or true polar wander (TPW) of visco-elastic bodies has been studied with approximated solutions. Two types of methods are commonly adopted: those based on the quasi-fluid approximation e.g. Ricard et al., 1993 and those based on the fluid limit approximation which only consider the final orientation (e.g. Matsuyama and Nimmo, 2007). Recently, Hu et al., (2017b) established a method which provides a dynamic solution for calculating the reorientation of tidally deformed bodies. However, they did not provide the links between the complete solution and the fluid limit solution. This paper provides a semi-analytical method for calculating the reorientation of tidally deformed bodies and shows the relation between the complete and the approximated solutions. Furthermore, we provide a criterion, the fluid limit process number F , to test for a given model and estimated TPW speed if the quasi-fluid approximation or fluid limit solution is valid. This number is a quantitative description of how close the body stays in hydrostatic equilibrium during a reorientation process. We use this number to obtain the largest allowed TPW speed of Mars as a function of viscosity, for which the approximated solutions may be used.

The dynamical environment on and about the Martian moon Phobos is explored. This planetary moon provides a unique dynamical environment in the solar system, being subject to extreme tidal forces and having a characteristically non-spherical shape. Further, it is not in a fully circular orbit, meaning that it has librations that arise from its eccentricity, contributing to a periodic forcing environment. Thus, to plan and implement missions in the vicinity of and on Phobos will require these considerations be taken into account. In this paper the latest published models of the Phobos shape and dynamics are used to characterize its dynamical environment in close proximity orbit about the body, for motion across its surface and for controlled hovering motion in its vicinity. It is found that surface motion is subject to a number of “speed limits” that can cause a moving vehicle to leave the surface and to possibly escape the moon and enter orbit about Mars. In terms of orbital stability, the existence of libration orbit families are characterized down to the surface using an exact potential, and the known stable QSO orbits are shown to be associated with families of stable quasi-periodic orbits.

X-ray irradiation in X-ray binaries is thought to control the behaviour at the outer disc, which is observable mainly at optical wavelengths. It is generally parametrized phenomenologically, but it can also be predicted from theoretical models of irradiated discs and their coronae/winds. We test these models using five multiwavelength Hubble Space Telescopeand quasi simultaneous Rossi X-Ray Timing Explorer (RXTE) data sets from the black hole binary XTE J1859 + 226. These show how the reprocessed optical emission changes during outburst as the source fades from the very high/intermediate state at ∼0.4LEdd down through the high/soft state towards the transition to the hard state at ∼0.02LEdd. The models are able to reproduce the small change in reprocessing efficiency as the source flux decreases by a factor of two, and the spectrum softens from the very high/intermediate state to the bright high/soft state. However, the low luminosity high/soft state as well as the transition spectrum show more complex behaviour that is not well described by current models. We suggest the disc geometry has changed drastically during the outburst, probably due to tidal forces, and that the disc is no longer in steady state at the late stage of the outburst.

Close-in planets evolve under extreme conditions, raising questions about their origins and current nature. Two predominant mechanisms are orbital migration, which brings them close to their star, and atmospheric escape under the resulting increased irradiation. Yet, their relative roles remain unclear because we lack models that couple the two mechanisms with high precision on secular timescales. To address this need, we developed the JADE code, which simulates the secular atmospheric and dynamical evolution of a planet around its star, and can include the perturbation induced by a distant third body. On the dynamical side, the 3D evolution of the orbit is modeled under stellar and planetary tidal forces, a relativistic correction, and the action of the distant perturber. On the atmospheric side, the vertical structure of the atmosphere is integrated over time based on its thermodynamical properties, inner heating, and the evolving stellar irradiation, which results, in particular, in photo-evaporation. The JADE code is benchmarked on GJ436 b, prototype of evaporating giants on eccentric, misaligned orbits at the edge of the hot Neptunes desert. We confirm that its orbital architecture is well explained by Kozai migration and unveil a strong interplay between its atmospheric and orbital evolution. During the resonance phase, the atmosphere pulsates in tune with the Kozai cycles, which leads to stronger tides and an earlier migration. This triggers a strong evaporation several Gyr after the planet formed, refining the paradigm that mass loss is dominant in the early age of close-in planets. This suggests that the edge of the desert could be formed of warm Neptunes whose evaporation was delayed by migration. It strengthens the importance of coupling atmospheric and dynamical evolution over secular timescales, which the JADE code will allow simulating for a wide range of systems., 20 pages, 2 figures, accepted in A&A

The standard model of stellar evolution in close binary systems assumes that during mass transfer episodes, the system is in a synchronized and circularized state. Remarkably, the redback system PSR J1723–2837 has an orbital period derivative P˙orb too large to be explained by this model. Motivated by this fact, we investigate the action of tidal forces in between two consecutive mass transfer episodes for a system under irradiation feedback, which is a plausible progenitor for PSR J1723–2837. We base our analysis on Hut’s treatment of equilibrium tidal evolution, generalized by considering the donor as a two layers object that may not rotate as a rigid body. We also analyse three different relations for the viscosity with the tidal forcing frequency. We found that the large value measured for P˙orb can be reached by systems where the donor star rotates slower (by few per cent) than the orbit just after mass transfer episodes. Van Staden & Antoniadis have observed this object and reported a lack of synchronism, opposite to that required by the Hut’s theory to account for the observed P˙orb⁠. Motivated by this discrepancy, we analyse photometric data obtained by the spacecraft Kepler second mission K2, with the purpose of identifying the periods present in PSR J1723–2837. We notice several periods close to those of the orbit and the rotation. The obtained periods pattern reveals the presence of a more complex phenomenology, which would not be well described in the frame of the weak friction model of equilibrium tides., Instituto de Astrofísica de La Plata

We explore the effect of anisotropic wind driving on the properties of accretion onto black holes in close binaries. We specifically focus on line-driven winds, which are common in high-mass X-ray binaries. In close binary systems, the tidal force from the companion star can modify the wind structure in two different ways. One is the reduction of wind terminal velocity due to the weaker effective surface gravity. The other is the reduction in mass flux due to gravity darkening. We incorporate these effects into the so-called CAK theory in a simple way and investigate the wind flow around the accretor on the orbital scale. We find that a focused accretion stream is naturally formed when the Roche lobe filling factor is $\gtrsim0.8$-0.9, analogous to that of wind Roche lobe overflow, but only when the velocity reduction is taken into account. The formation of a stream is necessary to bring in sufficient angular momentum to form an accretion disc around the black hole. Gravity darkening effects reduce the amount of accreted angular momentum, but not enough to prevent the formation of a disc. Based on these results, we expect there to be a discrete step in the observability of high-mass X-ray binaries depending on whether the donor Roche lobe filling factor is below or above $\sim$0.8-0.9., Comment: accepted for publication in the Publications of the Astronomical Society of Australia

Both classical and relativistic weak-field and slow-motion perturbations to planetary orbits can be treated as perturbative corrections to the Keplerian model. In particular, tidal forces and General Relativity (GR) induce small precession rates of the apsidal line. Accurate measurements of these effects in transiting exoplanets could be used to test GR and to gain information about the planetary interiors. Unfortunately, models for transiting planets have a high degree of degeneracy in the orbital parameters that, combined to the uncertainties of photometric transit observations, results in large errors on the determinations of the argument of periastron and precludes a direct evaluation of the apsidal line precession. Moreover, tidal and GR precession time-scales are many order of magnitudes larger than orbital periods, so that on the observational time-spans required to cumulate a precession signal enough strong to be detected, even small systematic errors in transit ephemerides add up to cancel out the tiny variations due to precession. Here we present a more feasible solution to detect tidal and GR precession rates through the observation of variations of the time interval ($\Delta \tau$) between primary and secondary transits of hot Jupiters and propose the most promising target for such detection, WASP-14 b. For this planet we expect a cumulated $\Delta \tau$ $\approx$ -250 s, due to tidal and relativistic precession, since its first photometric observations., Comment: This paper has been accepted for publication by MNRAS

We present a series of chemodynamical simulations of Magellanic-like systems consisting of two interacting, equal-mass dwarf galaxies orbiting a massive host galaxy, including feedback and star formation, tides, and ram pressure. We study the star formation and chemical enrichment history of the dwarfs, and the production of a Magellanic Stream analogue. The dwarfs interact with each other through tidal forces, distorting their morphologies and triggering star formation. A stream is naturally produced as outflows, induced by feedback and interactions, are stretched by tidal forces. Counter to some recent simulations, we find that the leading arm persists even in the presence of ram pressure from the host galaxy. Interactions between the dwarfs and the host galaxies produce multiple kinematic components in the stream, as observed. A combination of ongoing star-formation and entrained low-metallicity gas causes the stream to have a complex chemical structure, with an average metallicity that is generally lower than that of the dwarfs., Accepted to ApJ

Stellar components in binaries are subject to tidal forces which influence asteroseismic properties. Tidally pertubed pulsations have been reported for different objects but none of these are in their pre-main sequence phase of evolution. This makes RS Cha, consisting of two $\delta$ Scuti stars and with pulsational characteristics influenced by tidal effects , the first such object observed. We aim to investigate the pulsational properties of the eclipsing binary RS Cha in terms of the theory of tidally perturbed pulsations. Based on photometric time series obtained from the TESS satellite, we performed binary modelling using PHOEBE to interpret the binary light curve and to allow the investigation of the pulsations of both components in RS Cha. We modelled the detrended light curve with the superposition of linear modes. The frequencies were then interpreted as self excited modes perturbed by tidal forces. We find evidence for tidally perturbed modes, which enables the identification of pulsation modes. RS Cha mainly exhibits dipole modes, while one prominent $l=2$ or $l=3$ mode is also inferred. The latter verifies previous results from spectroscopic time series. This work shows that RS Cha is an ideal candidate to test the theory of tidally perturbed pulsations within the framework of asteroseismic modelling. The identification of multiple pulsation modes using this theory is unprecedented and will be a keystone in the future of pre-main sequence asteroseismology. However, amplitude modulation caused by the changing light ratio during the orbital phase in an eclipsing binary also plays a significant role, which can complicate mode identification., Comment: 11 pages plus 7 for the appendix, 8 pages plus 8 for the appendix, accepted for publication in Astronomy & Astrophysics

Aims. Planets with large moon(s) or those in the habitable zone of low-mass stars may experience much stronger tidal force and tide-induced ocean mixing than that on Earth. Thus, the vertical diffusivity (or, more precisely, diapycnal diffusivity) on such planets, which represents the strength of vertical mixing in the ocean, would be greater than that on Earth. In this study, we explore the effects of extremely high diffusivity on the ocean circulation and surface climate of Earth-like planets in one asynchronous rotation orbit. Methods. The response of planetary climate to 10 and 100 times greater vertical diffusivity than that found on Earth is investigated using a fully coupled atmosphere–ocean general circulation model. In order to perform a clear comparison with the climate of modern Earth, Earth’s orbit, land–sea configuration, and present levels of greenhouse gases are included in the simulations. Results. We find that a larger vertical diffusivity intensifies the meridional overturning circulation (MOC) in the ocean, which transports more heat to polar regions and melts sea ice there. Feedback associated with sea ice, clouds, and water vapor act to further amplify surface warming. When the vertical diffusivity is 10 (100) times the present-day value, the magnitude of MOC increases by ≈3 (18) times, and the global-mean surface temperature increases by ≈4 °C (10 °C). This study quantifies the climatic effect of an extremely strong vertical diffusivity and confirms an indirect link between planetary orbit, tidal mixing, ocean circulation, and surface climate. Our results suggest a moderate effect of varying vertical ocean mixing on planetary climate.

Phobos is one of the two moons of Mars, the only other terrestrial planet to have companions. In spite of certain common orbital characteristics between Earth’s Moon and Phobos, such as almost circular, low-inclination, and synchronous orbits, the origin of the Martian moons remains to be understood. This unsettled state of affairs reflects a scarcity of data pertinent to its interior structure and bulk composition, both of which hold the key to solving the question of the provenance of the satellites. Drawing on the importance of interior structure as a means of providing further constraints on this problem, we construct a series of models of the interior of Phobos in accordance with current observations and determine their tidal deformation. The models include a homogeneous body, ice-rock mixtures, models with positive and negative density gradients, and layered models. To compute the deformation of Phobos precisely, we solve the three-dimensional elastostatic problem to obtain the full displacement field. For this, we rely on a higher-order spectral-element method and to account for a correct representation of shape and resulting displacement field, we accurately mesh the figure of Phobos by employing the digital terrain model of Willner et al. (2014). To enable us to further distinguish between models, we also rely on currently available geophysical data (e.g., magnitude of libration in longitude, mean density, gravity field, and moments of inertia). Our results show that the homogeneous, ice-rock mixture, and negative density gradient models are largely degenerate and therefore difficult to separate, but that models with a density increase with depth can be differentiated from the former model families. From a purely deformational point of view, we find that for positive gradient and layered models, the largest deformation occurs on the rim of the Stickney crater, rather than immediately at the sub-Mars point on Phobos, as in the case of the other models, where tidal forces are largest. These observations will be of considerable interest and should provide important anchoring points for the future exploration of Phobos as envisaged with the Martian Moons Exploration mission.

During close approaches to planets or stars, the morphological and dynamical properties of rubble-pile small bodies can be modified, and some may catastrophically break up. This phenomenon is of particular interest for understanding the evolution and population of small bodies, and for making predictions for future encounters. Previous numerical explorations typically used methods that do not adequately represent the nature of rubble piles. The encounter outcomes and influence factors are still poorly constrained. Based on recent advances in modeling rubble piles, we aim to provide a better understanding of the tidal encounter processes through soft-sphere discrete element modeling (SSDEM) and to establish a database of encounter outcomes and the dependencies on encounter conditions and rubble-pile properties. We performed thousands of numerical simulations using the SSDEM implemented in the N-body code pkdgrav to study the dynamical evolution of rubble piles during close encounters with the Earth. The effects of encounter conditions, material strength, arrangement, and resolution of constituent particles are explored. Three typical tidal encounter outcomes are classified, i.e., deformation, mass shedding, and disruption, ranging from mild modifications to severe damages of the progenitor. The encounter speed and distance required for causing disruption events are much smaller than those predicted by previous studies, indicating a smaller creation rate of tidally disrupted small body populations. Extremely elongated fragments with axis ratios ~1:6 can be formed by moderate encounters. Our analyses of the spin-shape evolution of the largest remnants reveal reshaping mechanisms of rubble piles in response to tidal forces, consistent with stable rubble-pile configurations derived by continuum theory. A case study for SL9 suggests a low bulk density (0.2-0.3 g/cc) for its progenitor., 15 pages, 15 figures, Accepted for publication in A&A, author version

Context. The strong tidal dissipation in the couple Jupiter-Io is spread to all the moons involved in the Laplace resonance (Io, Europa, and Ganymede), leading to a migration of their orbits. Aims. We aim to characterize the future behavior of the Galilean satellites over the Solar System lifetime and to quantify the stability of the Laplace resonance. Since tidal dissipation makes possible the exit from the current resonances or capture into new ones, we investigate the capture of Callisto into resonance. Methods. We perform hundreds of propagations using an improved version of a recent semi-analytical model. As Ganymede moves outwards, it approaches the 2:1 resonance with Callisto, inducing a temporary chaotic motion in the system. For this reason, we draw a statistical picture of the outcome of the resonant encounter. Results. The system can settle into two distinct outcomes: A) a chain of three 2:1 two-body resonances (Io-Europa, Europa-Ganymede and Ganymede-Callisto), or B) a resonant chain involving the 2:1 two-body resonance Io-Europa plus at least one pure 4:2:1 three-body resonance, most frequently between Europa, Ganymede and Callisto. In case A (56\% of the simulations), the Laplace resonance is always preserved and the eccentricities remain confined to small values below 0.01. In case B (44\% of the simulations), the Laplace resonance is generally disrupted and the eccentricities of Ganymede and Callisto can increase up to about 0.1, making this configuration unstable and driving the system into new resonances. Conclusion. From our results, the capture of Callisto into resonance appears to be extremely likely (100\% of our simulations). Assuming the most recent estimate of the dissipation between Io and Jupiter, the resonant encounter happens at about 1.5 Gyrs from now. Therefore, the stability of the Laplace resonance is guaranteed at least up to about 1.5 Gyrs., Comment: 22 pages, 13 figures, 6 sections. Accepted for publication in Astronomy & Astrophysics

White dwarf stars that enter the tidal radius of black holes with masses $\lesssim$ $10^5$ M$_{\odot}$ are doomed to be ripped apart by tidal forces. Black holes in this mass range between stellar black holes and supermassive black holes have not been conclusively identified so the detection of a tidal disruption of a white dwarf would provide clear evidence for the existence of intermediate-mass black holes. In this review, we present a theoretical and observational overview of the transient events that result from the tidal disruptions of white dwarfs by intermediate-mass black holes. This includes discussion of the latest simulations and predicted properties, the results of observational searches, as well as a summary of the potential for gravitational wave emission to be detected with upcoming missions., Accepted for publication in Springer Space Science Reviews. Chapter in review, 'The Tidal Disruption of Stars by Massive Black Holes', editors: P. G. Jonker, S. Phinney, E. M. Rossi, S. v. Velzen, I. Arcavi & M. Falanga

The long-wavelength coherent overdensity and tidal force, which are not direct observables for a finite-volume survey, affect time evolution of cosmic structure formation and therefore clustering observables through the mode coupling. In this paper we develop an "anisotropic" separate universe (SU) simulation technique to simulate large-scale structure formation taking into account the effect of large-scale tidal force into the anisotropic expansion of local background. We modify the TreePM N-body simulation code to implement the anisotropic SU simulations, and then study the "response" function of matter power spectrum that describes how the matter power spectrum responds to the large-scale tidal effect as a function of wavenumber and redshift for a given global cosmology. We test and validate the SU simulation results from the comparison with the perturbation theory predictions and the results from high-resolution PM simulations. We find that the response function displays characteristic scale dependences over the range of scales down to nonlinear scales, up to k ~ 6 h/Mpc., 15 pages, 10 figures, 1 table, accepted for publication in MNRAS, added a figure and a subsection, changed to double column

When a star passes close to a supermassive black hole (BH), the BH's tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code H-AMR. We carry out the first grid-based simulation of a deep-penetration TDE ($\beta$=7) with realistic system parameters: a black-hole-to-star mass ratio of $10^6$, a parabolic stellar trajectory, and a nonzero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disk forms due to the dissipation of orbital energy with $\sim$20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream-disk interactions near the pericenter. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 hours and a flaring in the accretion rate. We also find that the disk is eccentric with mean eccentricity e$\approx$0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericenter. Although these partial intersections eject gas out of the orbital plane, an accretion disk still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disk formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57., Comment: Accepted to MNRAS on November 23rd, 2021, 23 pages, 25 figures, uses mnras.cls. Comments welcome. Movies available at https://www.youtube.com/playlist?list=PL7YbfRC6zxzAqIPXWbJgXpr4Wq_nwHiDu

When a star comes too close to a supermassive black hole, it gets torn apart by strong tidal forces in a tidal disruption event, or TDE. Half of the elongated stream of debris comes back to the stellar pericenter where relativistic apsidal precession induces a self-crossing shock. As a result, the gas gets launched into an outflow that can experience additional interactions, leading to the formation of an accretion disc. We carry out the first radiation-hydrodynamics simulations of this process, making use of the same injection procedure to treat the self-crossing shock as in our previous adiabatic study (Bonnerot & Lu 2020). Two sets of realistic parameters of the problem are considered that correspond to different strengths of this initial interaction. In both cases, we find that the injected matter has its trajectories promptly circularized by secondary shocks taking place near the black hole. However, the generated internal energy efficiently diffuses away in the form of radiation, which results in a thin vertical profile of the formed disc. The diffusing photons promptly irradiate the surrounding debris until they emerge with a bolometric luminosity of $L\approx 10^{44} \, \rm erg\, s^{-1}$. Towards the self-crossing shock, diffusion is however slowed that results in a shallower luminosity increase, with a potentially significant component in the optical band. Matter launched to large distances continuously gains energy through radiation pressure, which can cause a significant fraction to become unbound. This work provides direct insight into the origin of the early emission from TDEs, which is accessed by a rapidly increasing number of observations., 23 pages, 19 figures, accepted by MNRAS. Movies of the simulation are available at http://www.tapir.caltech.edu/~bonnerot/first-light.html

A tidal disruption event (TDE) occurs when a star passes too close to a supermassive black hole and gets torn apart by its gravitational tidal field. After the disruption, the stellar debris form an expanding gaseous stream. The morphology and evolution of this stream is particularly interesting as it ultimately determines the observational properties of the event itself. In this work we perform 3D hydrodynamical simulations of the TDE of a star modelled as a polytropic sphere of index {\gamma} = 5/3, and study the gravitational stability of the resulting gas stream. We provide an analytical solution for the evolution of the stream in the bound, unbound and marginally bound case, that allows us to describe the stream properties and analyse the time-scales of the physical processes involved, applying a formalism developed in star formation context. Our results are that, when fragmentation occurs, it is fueled by the failure of pressure in supporting the gas against its self-gravity. We also show that a stability criterion that includes also the stream gas pressure proves to be far more accurate than one that only considers the black hole tidal forces, giving analytical predictions of the time evolution of the various forces associated to the stream. Our results point out that fragmentation occurs on timescales longer compared with the observational windows of these events and is thus not expected to give rise to significant observational features., Comment: 12 pages,15 figures, accepted for publication in MNRAS

An origin of the Optical/UV radiation from tidal disruption events (TDEs) has recently been discussed for different scenarios, but observational support is generally missing. In this Letter, we test applicability of the `Wind-Driven model' (Uno & Maeda 2020) to a sample of UV/Optical TDEs. With the model, we aim to derive the physical properties of the Optical/UV TDEs, such as mass-loss rates and characteristic radii. The model assumes optically thick continuous outflows like stellar winds, and one key question is how the wind-launched radius is connected to physical processes in TDEs. As one possibility, through a comparison between the escape velocities estimated from their black-hole masses and the wind velocities estimated from observed line widths, we propose that the outflow is launched from the self-interaction radius ($R_{\rm SI}$) where the stellar debris stretched by the tidal force intersects; we show that the escape velocities at $R_{\rm SI}$ are roughly consistent with the wind velocities. By applying the model to a sample of Optical/UV TDE candidates, we find that explosive mass ejections ($\gtrsim 10 ~M_{\odot}{\rm yr^{-1}}$) from $R_{\rm SI}$ ($\sim 10^{14}{\rm ~cm}$) can explain the observed properties of TDEs around peak luminosity. We also apply the same framework to a peculiar transient, AT2018cow. The model suggests that AT2018cow is likely a TDE induced by an intermediate-mass black hole ($M_{\rm BH} \sim 10^{4}~M_{\odot}$)., 8 pages, 1 table, 4 figures, accepted for publication in the ApJL

The frequent detection of binary mergers of ∼30 M ⊙ black holes (BHs) by the Laser Interferometer Gravitational-Wave Observatory (LIGO) rekindled researchers’ interest in primordial BHs (PBHs) being dark matter (DM). In this work, we investigated PBHs distributed as DM with a monochromatic mass of 30 M ⊙ and examined the encounter-capture scenario of binary formation, where the densest central region of DM halo dominates. Thus, we paid special attention to the tidal effect by the supermassive black hole (SMBH) present. In doing so, we discovered a necessary tool called loss zone that complements the usage of loss cone. We found that the tidal effect is not prominent in affecting binary formation, which also turned out to be insufficient in explaining the totality of LIGO’s event rate estimation, especially due to a microlensing event constraining the DM fraction in PBH at the mass of interest from near unity to an order smaller. Meanwhile, an early-universe binary formation scenario proves so prevailing that the LIGO signal in turn constrains the PBH fraction below one percent. Thus, people should put more faith in alternative PBH windows and other DM candidates.

An extended gravitational object embedded in a parent system comprised of gas and collisionless particles may undergo both dynamical friction (DF) and mass loss by tidal forces. If the object is compact enough, it can increase its mass through accretion of material from the surrounding medium. We extend the classical linear analysis of DF on a constant-mass body in a gaseous medium to the case where its mass changes with time. We show that the structure of the wake may differ significantly from the constant-mass case. For instance, the front-back symmetry of density about subsonic constant-mass perturbers is broken down for variable-mass perturbers. The density wake keeps a memory of the past mass history of the perturber. For dissolving perturbers, the density field is more dense than expected using the instantaneous mass of the perturber in the classical formula. As a consequence, the instantaneous-mass approximation underestimates the drag force for mass-losing perturbers and overestimates it for mass-gaining perturbers. We present cases in which the percentage error in the drag force using the instantaneous-mass approximation is greater than 50%., Comment: 16 pages, 15 figures, accepted for publication in MNRAS

In hierarchical triple systems, the inner binary is perturbed by a distant companion. For large mutual inclinations, the Lidov-Kozai mechanism secularly excites large eccentricity and inclination oscillations of the inner binary. The maximal eccentricity attained, $e_{\rm max}$ is simply derived and widely used. However, for mildly hierarchical systems (i.e. the companion is relatively close and massive), non-secular perturbations affect the evolution. Here we account for fast non-secular variations and find new analytic formula for $e_{{\rm max}}$, in terms of the system's hierarchy level, correcting previous work and reproducing the orbital flip criteria. We find that $e_{{\rm max}}$ is generally enhanced, allowing closer encounters between the inner binary components, thus significantly changing their interaction and its final outcome. We then extend our approach to include additional relativistic and tidal forces. Using our results, we show that the merger time of gravitational-wave (GW) sources orbiting massive black-holes in galactic nuclei is enhanced compared with previous analysis accounting only for the secular regime. Consequently, this affects the distribution and rates of such GW sources in the relevant mild-hierarchy regime. We test and confirm our predictions with direct N-body and 2.5-level Post-Newtonian codes. Finally, we calculate the formation and disruption rates of hot-Jupiters (HJ) in planetary systems using a statistical approach, which incorporates our novel results for $e_{{\rm max}}$. We find that more HJ migrate from further out, but they are also tidally disrupted more frequently. Remarkably, the overall formation rate of HJs remains similar to that found in previous studies. Nevertheless, the different rates could manifest in different underlying distribution of observed warm-Jupiters., Comment: Accepted to MNRAS

We present a sample of 446 galaxy pairs constructed using the cosmological simulation IllustrisTNG-100 at z = 0, with M$_{FoF, dm}$ = 10$^{11}$-10$^{13.5}$ M$_{\odot}$. We produce ideal mock SDSS g-band images of all pairs to test the reliability of visual classification schema employed to produce samples of interacting galaxies. We visually classify each image as interacting or not based on the presence of a close neighbour, the presence of stellar debris fields, disturbed discs, and/or tidal features. By inspecting the trajectories of the pairs, we determine that these indicators correctly identify interacting galaxies $\sim45\%$ of the time. We subsequently split the sample into the visually identified interacting pairs (VIP; 38 pairs) and those which are interacting but are not visually identified (nonVIP; 47 pairs). We find that VIP have undergone a close passage nearly twice as recently as the nonVIP, and typically have higher stellar masses. Further, the VIP sit in dark matter haloes that are approximately 2.5 times as massive, in environments nearly 2 times as dense, and are almost a factor of 10 more affected by the tidal forces of their surroundings than the nonVIP. These factors conspire to increase the observability of tidal features and disturbed morphologies, making the VIP more likely to be identified. Thus, merger rate calculations which rely on stellar morphologies are likely to be significantly biased toward massive galaxy pairs which have recently undergone a close passage., 20 pages, 11 figures, thumbnail catalog of interacting pairs sample

The disruption of substructure in galaxy clusters likely plays an important role in shaping the cluster population as a significant fraction of cluster galaxies today have spent time in a previous host system, and thus may have been pre-processed. Once inside the cluster, group galaxies face the combined environmental effects from group and cluster - so called 'post-processing'. We investigate these concepts, by tracking the evolution of satellites and their hosts after entering the cluster and find that tidal forces during their first pericentric passage are very efficient at breaking up groups, preferentially removing satellites at larger distances from their hosts. 92.2% of satellites whose host has passed pericentre will leave their host by $z=0$, typically no later than half a Gyr after pericentric passage. We find satellites leave with high velocities, and quickly separate to large distances from their hosts, making their identification within the cluster population challenging. Those few satellites ($\sim$7.8%) that remain bound to their hosts after a pericentric passage are typically found close to their host centres. This implies that substructure seen in clusters today is very likely on first infall into the cluster, and yet to pass pericentre. This is even more likely if the substructure is extended, with satellites beyond R$_{200}$ of their host. We find the cluster dominates the tidal mass loss and destruction of satellites, and is responsible for rapidly halting the accretion of new satellites onto hosts once they reach 0.6-0.8 R$_{200}$ radii from the cluster., Comment: 14 pages, 14 figures, accepted for publication in MNRAS on 04.10.19

We introduce a new hybrid method to perform high-resolution tidal disruption simulations, at arbitrary orbits. An sph code is used to simulate tidal disruptions only in the immediate spatial domain of the star, namely, where the tidal forces dominate over gravity, and then during the fragmentation phase in which the emerging tidal stream may collapse under its own gravity to form fragments. Following each hydrodynamical simulation, an analytical treatment is then applied to instantaneously transfer each fragment back to the tidal sphere for its subsequent disruption, in an iterative process. We validate the hybrid model by comparing it to both an analytical impulse approximation model of single tidal disruptions, as well as full-scale sph simulations spanning the entire disc formation. The hybrid simulations are essentially indistinguishable from the full-scale sph simulations, while computationally outperforming their counterparts by orders of magnitude. Thereby our new hybrid approach uniquely enables us to follow the long-term formation and continuous tidal disruption of the planet/planetesimal debris, without the resolution and orbital configuration limitation of previous studies. In addition, we describe a variety of future directions and applications for our hybrid model, which is in principle applicable to any star, not merely white dwarfs.

The absorption of stellar radiation observed by the HD209458b in resonant lines of OI and CII has not yet been satisfactorily modeled. In our previous 2D simulations we have shown that the hydrogen-dominated upper atmosphere of HD209458b, heated by XUV radiation, expands supersonically beyond the Roche lobe and drags the heavier species along with it. Assuming solar abundances, OI and CII particles accelerated by tidal forces to velocities up to 50 km/s should produce the absorption due to Doppler resonance mechanism at the level of 6-10%, consistent with the observations. Since the 2D geometry does not take into account the Coriolis force in the planet reference frame, the question remained to which extent the spiraling of the escaping planetary material and its actually achieved velocity may influence the conclusions made on the basis of 2D modeling. In the present paper we apply for the first time in the study of HD209458b a global 3D hydrodynamic multi-fluid model that self-consistently describes the formation and expansion of the escaping planetary wind, affected by the tidal and Coriolis forces, as well as by the surrounding stellar wind. The modeling results confirm our previous findings that the velocity and density of the planetary flow are sufficiently high to produce the absorption in HI, OI, and CII resonant lines at the level close to the in-transit observed values. The novel finding is that the matching of the absorption measured in MgII and SiIII lines requires at least 10 times lower abundances of these elements than the Solar system values.

A terrestrial planet is molten during formation and may remain so if subject to intense insolation or tidal forces. Observations continue to favour the detection and characterisation of hot planets, potentially with large outgassed atmospheres. We aim to determine the radius of hot Earth-like planets with large outgassed atmospheres and explore differences between molten and solid silicate planets and their influence on the mass-radius relationship and transmission and emission spectra. An interior-atmosphere model, combined with static structure calculations, tracks the evolving radius of a rocky mantle that is outgassing CO$_2$ and H$_2$O. Synthetic emission and transmission spectra are generated for CO$_2$ and H$_2$O dominated atmospheres. Atmospheres dominated by CO$_2$ suppress the outgassing of H$_2$O to a greater extent than previously realised, as previous studies have applied an erroneous relationship between volatile mass and partial pressure. We therefore predict more H$_2$O can be retained by the interior during the later stages of magma ocean crystallisation. Furthermore, formation of a lid at the surface can tie outgassing of H$_2$O to the efficiency of heat transport through the lid, rather than the atmosphere's radiative timescale. Contraction of the mantle as it solidifies gives $\sim5\%$ radius decrease, which can partly be offset by addition of a relatively light species to the atmosphere. We conclude that a molten silicate mantle can increase the radius of a terrestrial planet by around $5\%$ compared to its solid counterpart, or equivalently account for a $13\%$ decrease in bulk density. An outgassing atmosphere can perturb the total radius according to its speciation. Atmospheres of terrestrial planets around M-stars that are dominated by CO$_2$ or H$_2$O can be distinguished by observing facilities with extended wavelength coverage (e.g., JWST)., Comment: 19 pages, published in A&A, abstract shortened

The planets with a radius l 4 R⊕ observed by the Kepler mission exhibit a unique feature, and propose a challenge for current planetary formation models. The tidal effect between a planet and its host star plays an essential role in reconfiguring the final orbits of the short-period planets. In this work, based on various initial Rayleigh distributions of the orbital elements, the final semi-major axis distributions of the planets with a radius l 4 R⊕ after suffering tidal evolutions are investigated. Our simulations have qualitatively revealed some statistical properties: the semi-major axis and its peak value all increase with the increase of the initial semi-major axis and eccentricity. For the case that the initial mean semi-major axis is less than 0.1 au and the mean eccentricity is larger than 0.25, the results of numerical simulation are approximately consistent with the observation. In addition, the effects of other parameters, such as the tidal dissipation coefficient, stellar mass and planetary mass, etc., on the final semi-major axis distribution after tidal evolution are all relatively small. Based on the simulation results, we have tried to find some clues for the formation mechanism of low-mass planets. We speculate that these low-mass planets probably form in the far place of protoplanetary disk with a moderate eccentricity via the type I migration, and it is also possible to form in situ.

Observational studies have revealed that galaxy pairs tend to have lower gas-phase metallicity than isolated galaxies. This metallicity deficiency can be caused by inflows of low-metallicity gas due to the tidal forces and gravitational torques associated with galaxy mergers, diluting the metal content of the central region. In this work we demonstrate that such metallicity dilution occurs in state-of-the-art cosmological simulations of galaxy formation. We find that the dilution is typically 0.1 dex for major mergers, and is noticeable at projected separations smaller than $40$ kpc. For minor mergers the metallicity dilution is still present, even though the amplitude is significantly smaller. Consistent with previous analysis of observed galaxies we find that mergers are outliers from the \emph{fundamental metallicity relation}, with deviations being larger than expected for a Gaussian distribution of residuals. Our large sample of mergers within full cosmological simulations also makes it possible to estimate how the star formation rate enhancement and gas consumption timescale behave as a function of the merger mass ratio. We confirm that strong starbursts are likely to occur in major mergers, but they can also arise in minor mergers if more than two galaxies are participating in the interaction, a scenario that has largely been ignored in previous work based on idealised isolated merger simulations., Submitted to MNRAS

The Small Magellanic Cloud (SMC) has recently been found to harbour more than two hundred per cent increase of its known cluster population. We provide here with solid evidence that such an unprecedented number of clusters could be largely overestimated. On the one hand, the fully-automatic procedure used to identify such an enormous cluster candidate sample did not recover ~ 50 per cent, in average, of the known relatively bright clusters located in the SMC main body. On the other hand, the number of new cluster candidates per time unit as a function of time results noticeably different to the intrinsic SMC cluster frequency (CF), which should not be the case if these new detections were genuine physical systems. We additionally found that the SMC CF varies spatially, in such a way that it resembles an outside-in process coupled with the effects of a relatively recent interaction with the Large Magellanic Cloud. By assuming that clusters and field stars share the same formation history, we showed for the first time that the cluster dissolution rate also depends on the position in the galaxy. The cluster dissolution results higher as the concentration of galaxy mass increases or external tidal forces are present., Comment: 8 pages, 5 figures. Accepted for publication in Monthly Noticies of the Royal Astronomical Society

We run a suite of hydrodynamics simulations of tidal disruption events (TDEs) of a white dwarf (WD) by a black hole (BH) with a wide range of WD/BH masses and orbital parameters. We implement nuclear reactions to study nucleosynthesis and its dynamical effect through release of nuclear energy. The released nuclear energy effectively increases the fraction of unbound ejecta. This effect is weaker for a heavy WD with 1.2 $\mathrm{M}_{\odot}$, because the specific orbital energy distribution of the debris is predominantly determined by the tidal force, rather than by the explosive reactions. The elemental yield of a TDE depends critically on the initial composition of a WD, while the BH mass and the orbital parameters also affect the total amount of synthesized elements. Tanikawa et al. (2017) find that simulations of WD-BH TDEs with low resolution suffer from spurious heating and inaccurate nuclear reaction results. In order to examine the validity of our calculations, we compare the amounts of the synthesized elements with the upper limits of them derived in a way where we can avoid uncertainties due to low resolution. The results are largely consistent, and thus support our findings. We find particular TDEs where early self-intersection of a WD occurs during the first pericentre passage, promoting formation of an accretion disc. We expect that relativistic jets and/or winds would form in these cases because accretion rates would be super-Eddington. The WD-BH TDEs result in a variety of events depending on the WD/BH mass and pericentre radius of the orbit., 13 pages, 13 figures, 2 tables, accepted to MNRAS

Formation of supermassive stars (SMSs) with mass ~10^4 Msun is a promising pathway to seed the formation of supermassive black holes in the early universe. The so-called direct-collapse (DC) model postulates that such an SMS forms in a hot gas cloud irradiated by a nearby star-forming galaxy. We study the DC SMS formation in a fully cosmological context using three-dimensional radiation hydrodynamics simulations. We initialize our simulations using the outputs of the cosmological simulation of Chon et al. (2016), where two DC gas clouds are identified. The long-term evolution over a hundred thousand years is followed from the formation of embryo protostars through their growth to SMSs. We show that the strength of the tidal force by a nearby galaxy determines the multiplicity of the formed stars and affects the protostellar growth. In one case, where a collapsing cloud is significantly stretched by strong tidal force, multiple star-disk systems are formed via filament fragmentation. Small-scale fragmentation occurs in each circumstellar disk, and more than 10 stars with masses of a few times 10^3 Msun are finally formed. Interestingly, about a half of them are found as massive binary stars. In the other case, the gas cloud collapses nearly spherically under a relatively weak tidal field, and a single star-disk system is formed. Only a few SMSs with masses ~ 10^4 Msun are found already after evolution of a hundred thousand years, and the SMSs are expected to grow further by gas accretion and to leave massive blackholes at the end of their lives., 19 pages, 20 figures, Submitted to MNRAS

A repeating fast radio burst (FRB), FRB 20180916B (hereafter FRB 180916), was reported to have a 16.35-day period. This period might be related to a precession period. In this paper, we investigate two precession models to explain the periodic activity of FRB 180916. In both models, the radio emission of FRB 180916 is produced by a precessing jet. For the first disk-driven jet precession model, an extremely low viscous parameter (i.e., the dimensionless viscosity parameter $\alpha \lesssim 10^{-8}$) is required to explain the precession of FRB 180916, which implies its implausibility. For the second tidal force-driven jet precession model, we consider a compact binary consists of a neutron star/black hole and a white dwarf; the white dwarf fills its Roche lobe and mass transfer occurs. Due to the misalignment between the disk and orbital plane, the tidal force of the white dwarf can drive jet precession. We show that the relevant precession periods are several days to hundreds of days, depending on the specific accretion rates and component masses. The duration of FRB 180916 generation in the binary with extremely high accretion rate will be several thousand years., Comment: 12 pages, 4 figures, accepted for publication in ApJ