Numerical Study of Colliding Winds in Massive Stars

We run a numerical experiment ejecting stellar winds in a very massive binary system measuring the properties of the resulting colliding wind structure and accreted mass onto the companion under different conditions. Colliding massive binaries interact and create a colliding wind structure with a shape that depends on the momentum ratio, orbital motion, distance between the stars, and other factors. We run simulations of a static LBV-WR binary and in each simulation abruptly varying the mass loss rate of the LBV from the fiducial value. The modified wind front propagates and interacts with the previous colliding wind structure, and modifies its shape. We calculate the emitted X-ray from the interaction and investigate the proprieties of the new shape. We derive the mass accretion rate onto the secondary, and find that it depends on the momentum ratio of the winds. We then add orbital velocity that reduces the mass accretion rate, a similar behaviour as the analytical estimates based on modified Bondi–Hoyle–Lyttleton. Creating a large set of simulations like those presented here can allow constraining parameters for specific colliding wind binaries and derive their stellar parameters and orbital solution.

Milankovic Pseudo-cycles Recorded in Sediments and Ice Cores Extracted by Singular Spectrum Analysis

Milankovic cycles describe the changes in the Earth's orbit and rotation axis and their impact on its climate over thousands of years. Singular Spectrum Analysis (SSA) is a signal processing method that is best known for its ability to find and extract pseudo-cycles in complex signals. In this short paper, we propose to apply it to three time series that have been proposed as geological reference time scales, in order to retrieve, compare and identify their Milankovic periodicities: (1) LR04, a stack of Plio-Pleistocene benthic microfossil records (Lisiecki and Raymo, 2005), (2) the CO2 and CH4 records from the Vostok ice core (Petit et al, 1999) and (3) the long-term orbital solution La04 for the insolation of Laskar et al (2004). The Vostok CO2 and CH4 series share the first 7 SSA components, three main ones (98, 104, 39 kyr), and four smaller ones (18, 22, 65, 180 kyr). CO2 displays a component at 28 kyr and a doublet at 61 and 62 kyr. CH4 displays a doublet near 50 kyr. 18/22 ky is a precession doublet, 62 kyr an insolation component, and 95/105 kyr an insolation/eccentricity doublet. The 49/50 kyr doublet in CH4 is not found in the orbital model. The SSA results for the La04 orbital solution are in excellent agreement with the values obtained by Laskar et al (2004). Four SSA components of obliquity are almost identical (rounded figures are 41, 54, 29 and 39 kyr). As far as eccentricity is concerned, the first five components are 404, 95, 124, 99, and 132 kyr. The next components are not found in our list of components for eccentricity, but they are in the SSA of insolation, at 2338, 970, 488 and 684 kyr. With more than 20 components, the LR04 stack is the richest series. In order of decreasing amplitude, one encounters 41, 95 and 75 kyr components. Next are smaller 39.5 and 53.6 kyr components, and a 22.4 kyr component. One recognizes one of the two main precession components, the doublet of obliquity components, a line at 47.4 kyr that is not found in any of the other spectra, and a doublet at 53.6 and 55.7 kyr, corresponding to the line at 54 kyr found in all four orbital quantities. Next comes a line at 63.6 kyr that may correspond to a line in insolation, CH4 and CO2. Then come components from eccentricity variations at 75.2, 94.5, 107.2, 132.1, 198.6 and 400.9 kyr. The remaining components of LR04 show up in La04. The “elusive ~200 kyr eccentricity cycle” of Hilgen et al (2020) is actually present in all three series, in the La04 orbital model as a 195 ± 6 kyr component of eccentricity and in LR04 as a 198.6 ± 5.6 kyr component. Finding not only the main expected Milankovic periodicities but also many “secondary” components with much smaller amplitudes gives confidence in our iterative SSA method (iSSA), on the quality of the La04 model and on the remarkable LR04 sedimentary stack, with more than 15 “ Milankovic periods”.

BRITE photometry and STELLA spectroscopy of bright stars in Auriga: Rotation, pulsation, orbits, and eclipses

Context. Knowing rotational and pulsational periods across the Hertzsprung-Russell diagram is of top priority for understanding stellar activity as a function of time. Aims. We aim to determine periods for bright stars in the Auriga field that are otherwise not easily accessible for ground-based photometry. Methods. Continuous photometry with up to three BRITE satellites was obtained for 12 targets and subjected to a period search. Contemporaneous high-resolution optical spectroscopy with STELLA was used to obtain radial velocities through cross correlation with template spectra as well as to determine astrophysical parameters through a comparison with model spectra. Results. The Capella red light curve was found to be constant over 176 days with a root mean square of 1 mmag, but the blue light curve showed a period of 10.1 ± 0.6 d, which we interpret to be the rotation period of the G0 component. From STELLA we obtained an improved orbital solution based on 9600 spectra from the previous 12.9 yr. We derive masses precise to ≈0.3% but 1% smaller than previously published. The BRITE light curve of the F0 supergiant ε Aur suggests 152 d as its main pulsation period, while the STELLA radial velocities reveal a clear 68 d period. An ingress of an eclipse of the ζ Aur binary system was covered with BRITE and a precise timing for its eclipse onset derived. A possible 70 d period fits the proposed tidal-induced, nonradial pulsations of this ellipsoidal K4 supergiant. η Aur is identified as a slowly pulsating B (SPB) star with a main period of 1.29 d and is among the brightest SPB stars discovered so far. The rotation period of the magnetic Ap star θ Aur is detected from photometry and spectroscopy with a period of 3.6189 d and 3.6177 d, respectively, likely the same within the errors. The radial velocities of this star show a striking non-sinusoidal shape with a large amplitude of 7 km s−1. Photometric rotation periods are also confirmed for the magnetic Ap star IQ Aur of 2.463 d and for the solar-type star κ1 Cet of 9.065 d, and also for the B7 HgMn giant β Tau of 2.74 d. Revised orbital solutions are derived for the eclipsing SB2 binary β Aur, which replaces the initial orbit dating from 1948 for the 27-year eclipsing SB1 ε Aur, and for the RS CVn binary V711 Tau, for which a spot-corrected orbital solution was achieved. The two stars ν Aur and ι Aur are found to be long-term, low-amplitude RV and brightness variables, but provisional orbital elements based on a period of 20 yr and an eccentricity of 0.7 could only be extracted for ν Aur. The variations of ι Aur are due to oscillations with a period of ≈4 yr.

The orbit of Triton with new precise observations and the INPOP19a ephemeris

Aims. The Gaia catalogue brings new opportunities and challenges to high-precision astronomy and astrometry. The precision of data reduction is therefore improved by a large number of reference stars with high-precision positions and proper motions. Numerous precise positions for Triton are obtained from the latest observations using the Gaia catalogue. Furthermore, the new INPOP19a planetary ephemeris, which also fits the observations from the Gaia Data Release 2, has recently become available. In this paper, a new orbit of Triton is calculated using the latest precise charge-coupled device (CCD) observations and the INPOP19a ephemeris. Methods. Triton’s orbital solution is calculated using a numerical integrator, while the orientation of Neptune’s pole in particular is obtained by integrating the simplified Euler’s equations of motion. We determine the orbit of Triton over 170 yr based on 11 040 Earth-based observations made between 1847 and 2016 and on Voyager 2 data. The positions of the Sun and planets are provided by the INPOP19a ephemeris. We compare our results to those from other previous works to check the influences on Triton’s orbit from different planetary ephemerides. Results. A new orbit of Triton is provided here. The root-mean-square of the residuals for the Earth-based CCD absolute observations are 0.102″ in right ascension and 0.142″ in declination. Although most different planetary ephemerides have large differences in Neptune’s position, the orbits of Triton using different planetary ephemerides are still close, under similar dynamical models. The Voyager 2 data add a constraint on Triton’s orbit here.

Dwarfs in the Milky Way halo outer rim: first infall or backsplash satellites?

ABSTRACT Leo T is a gas-rich dwarf located at $414\, {\rm kpc}$ (1.4Rvir) distance from the Milky Way (MW) and it is currently assumed to be on its first approach. Here, we present an analysis of orbits calculated backwards in time for the dwarf with our new code delorean, exploring a range of systematic uncertainties, e.g. MW virial mass and accretion, M31 potential, and cosmic expansion. We discover that orbits with tangential velocities in the Galactic standard-of-rest frame lower than $| \vec{u}_{\rm t}^{\rm GSR}| \le 63^{+47}_{-39}\, {\rm km}\, {\rm s}^{\rm -1}$ result in backsplash solutions, i.e. orbits that entered and left the MW dark matter halo in the past, and that velocities above $| \vec{u}_{\rm t}^{\rm GSR}| \ge 21^{+33}_{-21}\, {\rm km}\, {\rm s}^{\rm -1}$ result in wide-orbit backsplash solutions with a minimum pericentre range of $D_{\rm min} \ge 38^{+26}_{-16}\, {\rm kpc}$, which would allow this satellite to survive gas stripping and tidal disruption. Moreover, new proper motion estimates overlap with our orbital solution regions. We applied our method to other distant MW satellites, finding a range of gas stripped backsplash solutions for the gasless Cetus and Eridanus II, providing a possible explanation for their lack of cold gas, while only first infall solutions are found for the H i-rich Phoenix I. We also find that the cosmic expansion can delay their first pericentre passage when compared to the non-expanding scenario. This study explores the provenance of these distant dwarfs and provides constraints on the environmental and internal processes that shaped their evolution and current properties.

Resonance in the K2-19 system is at odds with its high reported eccentricities

ABSTRACT K2-19 hosts a planetary system composed of two outer planets, b and c, with size of 7.0 ± 0.2 R⊕ and 4.1 ± 0.2 R⊕, and an inner planet, d, with a radius of 1.11 ± 0.05 R⊕. A recent analysis of Transit-Timing Variations (TTVs) suggested b and c are close to but not in 3:2 mean motion resonance (MMR) because the classical resonant angles circulate. Such an architecture challenges our understanding of planet formation. Indeed, planet migration through the protoplanetary disc should lead to a capture into the MMR. Here, we show that the planets are in fact, locked into the 3:2 resonance despite circulation of the conventional resonant angles and aligned periapses. However, we show that such an orbital configuration cannot be maintained for more than a few hundred million years due to the tidal dissipation experienced by planet d. The tidal dissipation remains efficient because of a secular forcing of the innermost planet eccentricity by planets b and c. While the observations strongly rule out an orbital solution where the three planets are on close to circular orbits, it remains possible that a fourth planet is affecting the TTVs such that the four planet system is consistent with the tidal constraints.

Timing of the accreting millisecond pulsar IGR J17591–2342: evidence of spin-down during accretion

ABSTRACT We report on the phase-coherent timing analysis of the accreting millisecond X-ray pulsar IGR J17591–2342, using Neutron Star Interior Composition Explorer (NICER) data taken during the outburst of the source between 2018 August 15 and 2018 October 17. We obtain an updated orbital solution of the binary system. We investigate the evolution of the neutron star spin frequency during the outburst, reporting a refined estimate of the spin frequency and the first estimate of the spin frequency derivative ($\dot{\nu }\sim -7\times 10^{-14}$ Hz s−1), confirmed independently from the modelling of the fundamental frequency and its first harmonic. We further investigate the evolution of the X-ray pulse phases adopting a physical model that accounts for the accretion material torque as well as the magnetic threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. From this analysis we estimate the neutron star magnetic field Beq = 2.8(3) × 108 G. Finally, we investigate the pulse profile dependence on energy finding that the observed behaviour of the pulse fractional amplitude and lags as a function of energy is compatible with the down-scattering of hard X-ray photons in the disc or the neutron star surface.

The M31/M33 tidal interaction: a hydrodynamic simulation of the extended gas distribution

ABSTRACT We revisit the orbital history of the Triangulum galaxy (M33) around the Andromeda galaxy (M31) in view of the recent Gaia Data Release 2 proper motion measurements for both Local Group galaxies. Earlier studies consider highly idealized dynamical friction, but neglect the effects of dynamical mass loss. We show the latter process to be important using mutually consistent orbit integration and N-body simulations. Following this approach, we find an orbital solution that brings these galaxies to within ∼50 kpc of each other in the past, ∼6.5 Gyr ago. We explore the implications of their interaction using an N-body/hydrodynamical simulation with a focus on the origin of two prominent features: (1) M31’s Giant Stellar Stream; and (2) the M31–M33 H i filament. We find that the tidal interaction does not produce a structure reminiscent of the stellar stream that survives up to the present day. In contrast, the M31–M33 H i filament is likely a fossil structure dating back to the time of the ancient encounter between these galaxies. Similarly, the observed outer disc warp in M33 may well be a relic of this past event. Our model suggests the presence of a tidally induced gas envelope around these galaxies, and the existence of a diffuse gas stream, the ‘Triangulum stream’, stretching for tens of kpc from M33 away from M31. We anticipate upcoming observations with the recently commissioned, Five-hundred-metre Aperture Spherical radio Telescope that will target the putative stream in its first years of operation.