Search

There is a lack of exoplanets with sizes similar to Neptune orbiting their host stars with periods $\lesssim 3$ days -- hence the name ``sub-Jovian/Neptune desert''. Recently, several exoplanets have been confirmed to reside in the desert transforming it into a ``savanna'' with several ``giraffe'' planets (such as LTT 9779 b and TOI-674 b). The most prominent scenarios put forward for the explanation of the formation of the desert are related to the stellar irradiation destroying the primary atmosphere of certain specific exoplanets. We aim to present three targets (LTT 9779 b, TOI-674 b and WASP-156 b) which, when observed at wide wavelength ranges in infrared (IR), could prove the presence of these processes, and therefore improve the theories of planetary formation/evolution. We simulate and analyse realistic light curves of the selected exoplanets with PLATO/NCAM and the three narrow-band filters of Ariel (VISPhot, FGS1 and FGS2) based on TESS observations of these targets. We improved the precision of the transit parameters of the three considered planets from the TESS data. We find that the combination of the three narrow-band filters of Ariel can yield inner precision of $\lesssim 1.1\%$ for the planetary radii. Data from the three telescopes together will span decades, allowing the monitoring of changes in the planetary atmosphere through radius measurements. The three selected ``giraffe'' planets can be golden targets for Ariel, whereby the loss of planetary mass due to stellar irradiation could be studied with high precision, multi-wavelength (spectro-)photometry., Comment: Accepted for publication in MNRAS; 16 pages, 17 figures

We present the discovery of a new highly compact quadruple star system, TIC 219006972, consisting of two eclipsing binary stars with orbital periods of 8.3 days and 13.7 days, and an outer orbital period of only 168 days. This period is a full factor of 2 shorter than the quadruple with the shortest outer period reported previously, VW LMi, where the two binary components orbit each other every 355 days. The target was observed by TESS in Full-Frame Images in sectors 14-16, 21-23, 41, 48 and 49, and produced two sets of primary and secondary eclipses. These show strongly non-linear eclipse timing variations (ETVs) with an amplitude of $\sim$0.1 days, where the ETVs of the primary and secondary eclipses, and of the two binaries are all largely positively correlated. This highlights the strong dynamical interactions between the two binaries and confirms the compact quadruple configuration of TIC 219006972. The two eclipsing binaries are nearly circular whereas the quadruple system has an outer eccentricity of about 0.25. The entire system is nearly edge-on, with a mutual orbital inclination between the two eclipsing binary star systems of about 1 degree., 16 pages, 13 figures, 4 tables, MNRAS accepted

Open-set recognition models, in addition to generalizing to unseen instances of known categories, have to identify samples of unknown classes during the training phase. The main reason the latter is much more complicated is that there is very little or no information about the properties of these unknown classes. There are methodologies available to handle the unknowns. One possible method is to construct models for them by using generated inputs labeled as unknown. Generative adversarial networks are frequently deployed to generate synthetic samples representing unknown classes to create better models for known classes. In this paper, we introduce a novel approach to improve the accuracy of recognition methods while reducing the time complexity. Instead of generating synthetic input data to train neural networks, feature vectors are generated using the output of a hidden layer. This approach results in a less complex structure for the neural network representation of the classes. A distance-based classifier implemented by a convolutional neural network is used in our implementation. Our solution’s open-set detection performance reaches an AUC value of 0.839 on the CIFAR-10 dataset, while the closed-set accuracy is 91.4%, the highest among the open-set recognition methods. The generator and discriminator networks are much smaller when generating synthetic inner features. There is no need to run these samples through the first part of the classifier with the convolutional layers. Hence, this solution not only gives better performance than generating samples in the input space but also makes it less expensive in terms of computational complexity.

This paper is the second part of a series of studies discussing a novel attitude determination method for nano-satellites. Our approach is based on the utilization of thermal imaging sensors to determine the direction of the Sun and the nadir with respect to the satellite with sub-degree accuracy. The proposed method is planned to be applied during the Cubesats Applied for MEasuring and LOcalising Transients (CAMELOT) mission aimed at detecting and localizing gamma-ray bursts with an efficiency and accuracy comparable to large gamma-ray space observatories. In our previous work we determined the spherical projection function of the MLX90640 infrasensors planned to be used for this purpose. We showed that with the known projection function the direction of the Sun can be located with an overall accuracy of $$\sim 40^\prime$$ ∼ 40 ′ . In this paper we introduce a simulation model aimed at testing the applicability of our attitude determination approach. Its first part simulates the orbit and rotation of a satellite with arbitrary initial conditions while its second part applies our attitude determination algorithm which is based on a multiplicative extended Kalman filter. The simulated satellite is assumed to be equipped with a GPS system, MEMS gyroscopes and the infrasensors. These instruments provide the required data input for the Kalman filter. We demonstrate the applicability of our attitude determination algorithm by simulating the motion of a nano-satellite on Low Earth Orbit. Our results show that the attitude determination may have a 1$$\sigma$$ σ error of $$\sim 30'$$ ∼ 30 ′ even with a large gyroscope drift during the orbital periods when the infrasensors provide both the direction of the Sun and the Earth (the nadir). This accuracy is an improvement on the point source detection accuracy of the infrasensors. However, the attitude determination error can get as high as 25$$^{\circ }$$ ∘ during periods when the Sun is occulted by the Earth. We show that following an occultation period the attitude information is immediately recovered by the Kalman filter once the Sun is observed again.

The Cubesats Applied for MEasuring and LOcalising Transients (CAMELOT) initiative proposes to deploy a fleet of 3U nanosatellites in order to localise GRBs with all-sky coverage. The operation is based on measuring the time delay of the event trigger between satellites that are otherwise uniformly distributed around the Earth in low-Earth orbit (between 500 - 600 km of altitude). In this design, caesium-iodide crystals interact with soft gamma radiation by emitting optical photons. Utilization of this effect, each member of the fleet is equipped with four of such scintillators and the emitted photons are detected by multi-pixel photon counters (MPPCs). Precise timing is crucial for this concept, the timestamping of the events and the synchronisation is provided by GPS. In order to demonstrate the feasibility of the CAMELOT concept, a single-unit CubeSat, named "GRBAlpha" is currently being developed. GRBAlpha is equipped with a single block of scintillator but the other subsystems are all the same as it will be on the CAMELOT units. We describe this single-unit platform system, focusing on the model versions suitable for high-altitude stratospheric balloon flights. This model has a standardized layout (including pin-out configuration, signalling and bus communication) and compatible with significant proportion of CubeSat system vendors. This system of ours is also capable of hosting multiple payloads at the same time, optimizing the utilization of balloon experiments.

Gaia19fct is one of the Gaia-alerted eruptive young stars that has undergone several brightening events. We conducted monitoring observations using multi-filter optical and near-infrared photometry, as well as near-infrared spectroscopy, to understand the physical properties of Gaia19fct and investigate whether it fits into the historically defined two classes. We present the analyses of light curves, color variations, spectral lines, and CO modeling. The light curves show at least five brightening events since 2015, and the multi-filter color evolutions are mostly gray. The gray evolution indicates that bursts are triggered by mechanisms other than extinction. Our near-infrared spectra exhibit both absorption and emission lines and show time-variability throughout our observations. We found lower rotational velocity and lower temperature from the near-infrared atomic absorption lines than from the optical lines, suggesting that Gaia19fct has a Keplerian rotating disk. The CO overtone features show a superposition of absorption and emission components, which is unlike other young stellar objects. We modeled the CO lines, and the result suggests that the emission and absorption components are formed in different regions. We found that although Gaia19fct exhibits characteristics of both types of eruptive young stars, FU Orionis-type objects (FUors) and EX Lupi-type objects (EXors), it shows more similarity with EXors in general., 29 pages, 14 figures, 7 tables, Accepted for publication in ApJ

The K2 mission of the Kepler Space Telescope allowed the observations of light curves of small solar system bodies throughout the whole Solar system. In this paper we present the results of a collection of K2 transneptunian object observations, between Campaigns C03 (November 2014 -- February 2015) to C19 (August -- September, 2018), which includes 66 targets. Due to the faintness of our targets the detectability rate of a light curve period is $\sim$56%, notably lower than in the case of other small body populations, like Hildas or Jovian trojans. We managed to obtain light curve periods with an acceptable confidence for 37 targets; the majority of these cases are new identifications. We were able to give light curve amplitude upper limits for the other 29 targets. Several of the newly detected light curve periods are longer than $\sim$24 h, in many cases close to $\sim$80 h, i.e., these targets are slow rotators. This relative abundance of slowly rotating objects is similar to that observed among Hildas, Jovian trojans and Centaurs in the K2 mission, and also among main belt asteroids measured with the TESS Space Telescope. Transneptunian objects show notably higher light curve amplitudes at large (D $\gtrsim$ 300 km) sizes than that found among large main belt asteroids, in contrast to the general expectation that due to their lower compressive strength they reach hydrostatic equlibrium at smaller sizes than their inner solar system counterparts., Accepted for publication in ApJS

Earlier reviews of trans-Neptunian light curves reported mean rotation periods of P = 7-8 h (Duffard et al., 2009), and it was also found that the binary trans-Neptunian population rotates slower (Thirouin & Sheppard, 2014), and objects in the cold classical population have larger variability and rotate slower than the non-cold classical TNOs (Benecchi et al., 2013,Thirouin & Sheppard, 2019). While ground-based observations have obvious limitations in detecting long-period light curves the K2 mission of the Kepler Space Telescope allowed long (up to ~80 days), uninterrupted observations of many Solar system objects, including main belt asteroids, Hildas, Jovian trojans, and also the irregular satellites of giant planets. Light curves were also published for a few, selected trans-Neptunian objects based on K2 observations (see Kiss et al., 2020, for a summary). A common outcome of the studies of larger samples, across all dynamical classes, was the identification of an increased number of targets with long rotation periods compared to previous ground-based studies. A similar trend is observed among the data of nearly 10 000 main belt asteroids obtained by the TESS Space Telescope (Pál et al., 2020), and asteroids with long rotation periods were identified in other surveys like the Asteroid Terrestrial-impact Last Alert System (ATLAS), the Zwicky Transient Facility (Erasmus et al., 2021) and the All-Sky Automated Survey for Supernovae (Hanus et al., 2021). We have collected the K2 trans-Neptunian object observations between Campaigns C03 (November 2014 -- February 2015) to C19 (August -- September, 2018), which includes 67 targets. Due to the faintness of our targets the detectability rate of a light curve period is ~57 %, notably lower than in the case of other small body populations, like Hildas or Jovian trojans. We managed to obtain light curve periods with an acceptable confidence for 36 targets; the majority of these cases are new identifications. We were able to give light curve amplitude upper limits for the other 31 targets. Several of the newly detected light curve periods are longer than ~24 h, in many cases close to ~100 h, i.e., slow rotators. There is a very significant difference between the rotation rates of the LCDB and K2 TNO samples (Figs. 1 and 2). The mean LCD spin frequency is 2.71 c/d (8.8 h), while it is 0.87 c/d (27.6 h) in the K2 sample which is more similar to the K2 Hilda and Jovian Trojan spin frequency distrbutions. Thirouin & Sheppard (2019) obtained 9.48±1.53 h and 8.45±0.58 h mean rotation periods for the cold classical and the non-cold classical TNOs. Our mean values for the same dynamical groups (but using different targets) are notably longer: 1.21+1.58-0.63 c/d (19.83 h) and 0.83+1.81-0.23 c/d (P=28.91 h), respectively. The K2 mean frequency is higher than that of the K2 Jovian Trojans and Hildas, but we could not detect the very long period targets that were observed in these other K2 samples. Figure 1: frequency distribution of asteroids. The cyan, magenta, green and blue colours represent the TNOs in the LCDB and Jovian trojans, Hildas and TNOs from K2, respectively. Figure 2: Frequency as a function of absolute magnitude. Big circles with error bars mark the median values standard deviations for the different samples. The horizontal dashed lines mark the spin frequencies of fast, slow, and very slow rotators (Pravec et al., 2002). While there are only three objects with D>500 km in our sample, there are a number of objects -- both with and without detected light curve periods -- that fall in the 300≤D≤500 km transitional zone where asphericity -- hence light curve amplitude -- is expected to drop assuming a single rotating body, assuming main belt composition. Main belt asteroids are already almost extinct in this size range, and so are Centaurs -- for these bodies irregular shapes are expected in most cases. While the general trend is that larger objects have smaller light curve amplitudes among TNOs -- a trend followed both by our sample and the LCDB TNOs -- there are a considerable number of TNOs with high asphericity in the 300≤D≤500 km size range. This contradiction could be resolved if TNOs had higher-than-expected compressive strength and become spherical for sizes larger than their main belt counterparts, and remain 'irregular' in the 300≤D≤500 km range. However, their general low density and high porosity point against this scenario. A notable fraction of contact or semi-contact binary systems in which the members themselves are in hydrostatical equlibrium could produce a population of high-amplitude light curves in this size range (Lacerda et al., 2006, 2014). As some authors pointed out, contact binaries may be very frequent, especially in the plutino population (Thirouin & Sheppard, 2018,2019). The long term stability of such systems against their tidal evolution, however, should be investigated to answer the reliability of this assumption. Finally, spherical (rotationally flattened) bodies with large albedo variegations could also explain the observed amplitudes. While the general expectation in most TNO light curve studies was a double peak light curve, in our sample most light curves were found to be single-peak, after comparing the single-peak and double-peak solutions. Fgiure 3: Light curve amplitude versus the estimated size of the targets in our sample. The region between the vertical dashed lines mark the irregular-to-spherical transition size range in the main belt. Blue and red symbols mark the K2 targets and K2 upper limits, small grey symbols correspond to main belt asteroids. Large gray symbols represent the theoretical maximum light curve amplitudes the of large main belt objects if it was solely caused by the elongated shape of a body with homogeneous albedo.

Eris is currently the most massive known dwarf planet in the Solar system, it has one known satellite, Dysnomia (Brown & Schaller 2007; Holler et al. 2020). There have been several attempts trying to identify the rotation period of Eris from visible ground based measurements which resulted in a wide range of possible values (Duffard et al., 2008). Here we present some new light curve data of Eris, taken with ∼1m-class ground based telescopes, with the GROND instrument at the 2.2m MPG/ESO telescope La Silla, and also with the Transiting Exoplanet Survey Satellite (TESS). Observations of Eris with GROND were made in 3 nights in August 2010. Auxiliary ground based photometry data of Eris from five telescopes were obtained in irregular intervals between 2005 and 2020 (see table 1.) TESS observed Eris in Sector 30 with its Camera 1 and CCD 3. A significant portion of the light curve data had to be excluded due to Eris’ encounter with nearby background sources which left a dominant feature in the background-subtracted image. We note that due to the limited length of the TESS light blocks (2.3 d and 6.6 d) it was not possible to detect light curve periods in to order or close to the orbital period (15.78 d). For most of the ground-based measurements we used our night-averaged values for the 1.5m, La Hita and CA2.2m data), except for the GROND measurements, for which we used the all the individual data points. In addition, we also considered previously published data, including ground-based data from Carraro et al., 2006 and Duffard et al., 2008, and the SWIFT satellite data from Roe et al., 2008. We assumed that i) the light curve amplitude is the same in any of the photometric bands used and ii) the light curve can be described by a simple sinusoidal variation. With these assumptions each model light curve can be described by four parameters: a light curve amplitude, period, phase-shift, and an offset from the photometric zero point. We allowed a different offset for each measurement block even if the data were taken by the same instrument and filter combination (1.5m measurements) due to the year-long gaps between the measurement blocks. We chose an amplitude A and period P, and determined the best fitting light curve phase using a Levenberg–Marquardt minimization algorithm. We expect that the best-fitting period-amplitude values provide the lowest C(P,A) values. We searched the period range P ∈ [1d, 15.88 d], where 15.88 d is the orbital period of Dysnomia, and it would correspond to a synchronised rotation. The C(P,A) contour map shows two minima, one at P ≈ 8 h, and A ≈ 0.05, and another at P ≈ 11.5 h, and A ≈ 0.08. To check the robustness of this result, we repeated the process by modifying the photometric data points by adding a random value with normal distribution and repeating the fitting process several times for the whole data set. The same two minima popped up in all cases. TESS data favours a rotation period of 59.420±0.527 h (See Figure 1.), while ground based data is more compatible with a synchronous or nearly synchronous rotation (See Figure 2.). We inspected the colours calculated from the GROND photometry and compared our results to previous values from the literature (Table 2). The optical colours (B-V, V-R, R-I) are in good agreement with previous colours, but the J-H colour from GROND is different from the ones in the literature. Also, we can see a big scatter in those. (See Table 2.) One possible explanation can be, that there is a big CH4 feature on the surface of Eris, which has strong absorption lines in H band, and it can change the NIR colour, but not the optical ones. To confirm this, more long-term optical and NIR measurements are needed. Finally, we use the simple tidal evolution code by Hastings et al. (2016) to calculate the evolution of the satellite orbit (only the semi-major axis in this approximation) and the spin evolution of Eris and Dysnomia. Although the spin period of Eris is not well defined yet, it is an interesting question whether the possible rotation periods of Eris – ranging from a few hours to the orbit-synchronous state – are feasible in terms of tidal evolution, using the current knowledge on the system components.

Using multi-directional measurements of the Transiting Exoplanet Survey Satellite we have been able to perform analysis on small bodies in the Solar System. In order to validate a purely TESS-based shape and spin-axis reconstruction modelling we determined the spin axis of about forty objects, selected from our large light curve sample. We chose those targets that had been measured in at least three TESS sectors and showed an amplitude variation large enough to constraint the direction of the spin axis. In our simple model the asteroid's shape is approximated with a triaxial ellipsoid, rotating around its shortest axis, with a fixed spin axis direction. Seen from different directions, the light curves and amplitudes would also be different, and we minimize the difference between the predicted and measured amplitudes to obtain the best fitting spin axis direction (lambda_p, beta_p) and shape parameters. Our spin axis solutions are compared to those in the DAMIT database which are based on multi-epoch, mainly ground-based observations. This benchmark study will indicate the capabilities in obtaining shape and spin axis solutions from large surveys. Extracting this information for several such asteroids we can obtain shape and spin axis direction distributions which will be used to infer the formation and dynamics of the asteroid belt using formation and collisional evolution models, for a significantly larger number of targets then before, and using high quality light curve data. Figure 1: Light curves of the asteroid (4717) Kaneko observed in six different sectors of TESS Figure 2: χ2 contours obtained from the comparision of observed and predicted light curve amplitudes for (4717) Kaneko. The lowest contours represent our best fitting spin axis solutions in ecliptic coordinates. Circles mark the spin axis solutions provided by DAMIT.

Hi’iaka is the largest satellite of the dwarf planet Haumea, with an estimated area-equivalent diameter of 300 km (Fernández-Valenzuela et al., 2021). It is the best studied satellite in the trans-Neptunian region. Its rotational light-curve was observed with Hubble, for which an approximate rotation period of 9.8 h was obtained (Hastings et al. 2016). The system is very peculiar because it stands out from all other TNO-binary systems. While all other known satellites are thought to be synchronous, Hi’iaka’s rotation period is fast compared to the 49 days that takes to complete an orbit around Haumea. Therefore, the study of Haumea-Hi’iaka system yields important information about the formation processes of the whole Haumea’s system, which includes another moon (Brown et al. 2006), a ring (Ortiz et al. 2017) and a family of objects (Brown et al. 2007). Our group has been observing Haumea since its discovery, compiling a large database of images since around 20 years ago. Using this set of images we have obtained high accuracy astrometric measurements of the photocenter of the Haumea-Hi’iaka system. We have applied a similar procedure as in Ortiz et al. (2017) to disentangle the position of Haumea from the contribution of Hi'iaka, but for a much larger time span as mentioned above. Therefore, we have been able to determine more accurate orbits for Haumea and Hi'iaka. Additionally, we have carried out two specific observational runs of several days in order to obtain the rotational phase of Hi’iaka at that moment of the two stellar occultations that occurred last year (in April 2021). We used the 1.23-m telescope at Calar Alto Observatory, the Artemis telescope at Teide Observatory and the 1.5-m telescope at Sierra Nevada Observatory to acquire images of the unresolved system. The resulting photometry of these images give rise two rotational light-curves of Haumea in 2021 and 2022. We fitted a fourth-order Fourier function, which represents Haumea’s body-shape contribution to the rotational light-curves. From this fit, we took the residuals of the observational data and searched for periodicities within them. We obtained a rotation period in agreement with the estimations in Hastings et al. (2016), but much more accurate. These residuals, when folded to the resulting period, provide Hi’iaka’s rotational light-curve. The amplitude obtained for Hi’iaka’s rotational light-curve is 0.015 mag, which agrees with the expected signal induced in Haumea’s rotational light-curve when accounting for a variable source as that produced by Hi’iaka, i.e., considering the rotational light-curve obtained in Hastings et al. (2016). We have not detected a change in the amplitude of Hi’iaka’s rotational light-curve when comparing our data, taken in 2021 and 2022, with those from Hastings et al. (2016), taken in 2010. This means that the obliquity of Hi’iaka must be close to 90º in its orbit around Haumea. Brown et al. (2006), The Astrophysical Journal, Volume 639, Issue 1, pp. L43-L46. Brown et al. (2007), Nature, Volume 446, Issue 7133, pp. 294-296. Fernández-Valenzuela et al. (2021), AAS Division of Planetary Science meeting #53, id. 503.05. Bulletin of the American Astronomical Society, Vol. 53, No. 7 e-id 2021n7i503p05. Hastings et al. (2016), The Astronomical Journal, Volume 152, Issue 6, article id. 195, 12 pp. Ortiz et al. (2017), Nature, Volume 550, Issue 7675, pp. 219-223.

The Transiting Exoplanet Survey Satellite (TESS) produces a series of large field-of-view images with a size of 96x24 degress using a cadence of 10 minutes. Even when these field-of-views do not cover directly the ecliptic, thousands of asteroids are available for good quality photometry on a monthly basis. In this poster, we display the structure of the multi-stage pipeline designed to extract precise and accurate photometry from asteroids in an effective manner. The corresponding technical aspects include but not limited to the layout of the backend storage system, staging storages, data product indexing, memory utilization, metadata handling, intermediate products for shift & stack and efficient object querying implementation via SQL. Outputs of these pipeline stages are employed many presentations of this conference as perts of the early second data release.

Classical T Tauri stars are low-mass objects, which are still accreting material from the surrounding circumstellar disk. The accretion process is essential in the formation of Sun-like stars and in setting the properties of the disk at the time when planet formation occurs. We constructed a complex dataset in order to examine the accretion process of VW Cha, a classical T Tauri multiple system with the aim of studying the physical origin of the photometric and spectroscopic variability of the system. The TESS Space Telescope observed VW Cha between 2019 April 22 and June 19, and we complemented these data with contemporaneous ground-based $I_CJHK$ band photometric measurements. In addition, we obtained high-resolution optical spectra with the VLT/ESPRESSO and the 2.2\,m/FEROS instruments. Analyzing these data, we found that the TESS light curve shows photometric variations on timescales from minutes to weeks with a peak-to-peak amplitude of $\sim$0.8 mag. The near-infrared light curves follow the shape of the optical measurements, however, the peak-to-peak amplitudes are slightly increasing towards the longer wavelengths. We took spectra in both fainter and brighter photometric states of the system, allowing us to examine the origin of a photometric brightening event. Our results show that this brightening event can be explained by increased accretion. In addition, our spectroscopic data also suggest that the primary component of VW Cha is a spectroscopic binary, as it was proposed in earlier works., Accepted for publication in ApJ. 18 pages, 10 figures

We describe results of a planetary defense exercise conducted during the close approach to Earth by the near-Earth asteroid (99942) Apophis during 2020 December–2021 March. The planetary defense community has been conducting observational campaigns since 2017 to test the operational readiness of the global planetary defense capabilities. These community-led global exercises were carried out with the support of NASA’s Planetary Defense Coordination Office and the International Asteroid Warning Network. The Apophis campaign is the third in our series of planetary defense exercises. The goal of this campaign was to recover, track, and characterize Apophis as a potential impactor to exercise the planetary defense system including observations, hypothetical risk assessment and risk prediction, and hazard communication. Based on the campaign results, we present lessons learned about our ability to observe and model a potential impactor. Data products derived from astrometric observations were available for inclusion in our risk assessment model almost immediately, allowing real-time updates to the impact probability calculation and possible impact locations. An early NEOWISE diameter measurement provided a significant improvement in the uncertainty on the range of hypothetical impact outcomes. The availability of different characterization methods such as photometry, spectroscopy, and radar provided robustness to our ability to assess the potential impact risk.

We report the discovery of a compact triply eclipsing triple star system in the southern continuous viewing zone of the TESS space telescope. TIC 278825952 is a previously known, but unstudied circular eclipsing binary with a period of 4.781 d with a tertiary component in a wider, circular orbit of 235.55-d period that was found from three sets of third-body eclipses and from light travel-time effect dominated eclipse timing variations. We performed a joint photodynamical analysis of the eclipse timing variation curves, photometric data, and the spectral energy distribution, coupled with the use of PARSEC stellar isochrones. We find that the inner binary consists of slightly evolved, near twin stars of masses of 1.12 and 1.09 M⊙ and radii of 1.40 and 1.31 R⊙. The third, less massive star has a mass of 0.75 M⊙ and radius of 0.70 R⊙. The low mutual inclination and eccentricities of the orbits show that the system is highly coplanar and surprisingly circular.

We present extensive observations of SN 2018zd covering the first $\sim450$\,d after the explosion. This SN shows a possible shock-breakout signal $\sim3.6$\,hr after the explosion in the unfiltered light curve, and prominent flash-ionisation spectral features within the first week. The unusual photospheric temperature rise (rapidly from $\sim 12,000$\,K to above 18,000\,K) within the earliest few days suggests that the ejecta were continuously heated. Both the significant temperature rise and the flash spectral features can be explained with the interaction of the SN ejecta with the massive stellar wind ($0.18^{+0.05}_{-0.10}\, \rm M_{\odot}$), which accounts for the luminous peak ($L_{\rm max} = [1.36\pm 0.63] \times 10^{43}\, \rm erg\,s^{-1}$) of SN 2018zd. The luminous peak and low expansion velocity ($v \approx 3300$ km s$^{-1}$) make SN 2018zd to be like a member of the LLEV (luminous SNe II with low expansion velocities) events originated due to circumstellar interaction. The relatively fast post-peak decline allows a classification of SN 2018zd as a transition event morphologically linking SNe~IIP and SNe~IIL. In the radioactive-decay phase, SN 2018zd experienced a significant flux drop and behaved more like a low-luminosity SN~IIP both spectroscopically and photometrically. This contrast indicates that circumstellar interaction plays a vital role in modifying the observed light curves of SNe~II. Comparing nebular-phase spectra with model predictions suggests that SN 2018zd arose from a star of $\sim 12\,\rm M_{\odot}$. Given the relatively small amount of $^{56}$Ni ($0.013 - 0.035 \rm M_{\odot}$), the massive stellar wind, and the faint X-ray radiation, the progenitor of SN 2018zd could be a massive asymptotic giant branch star which collapsed owing to electron capture., Accepted for publication in MNRAS, 20 pages, 11 figures

Extreme debris discs are characterized by unusually strong mid-infrared excess emission, which often proves to be variable. The warm dust in these discs is of transient nature and is likely related to a recent giant collision occurring close to the star in the terrestrial region. Here we present the results of a 877 days long, gap-free photometric monitoring performed by the Spitzer Space Telescope of the recently discovered extreme debris disc around TYC 4209-1322-1. By combining these observations with other time-domain optical and mid-infrared data, we explore the disc variability of the last four decades with particular emphasis on the last 12 yr. During the latter interval the disc showed substantial changes, the most significant was the brightening and subsequent fading between 2014 and 2018 as outlined in WISE data. The Spitzer light curves outline the fading phase and a subsequent new brightening of the disc after 2018, revealing an additional flux modulation with a period of ~39 days on top of the long-term trend. We found that all these variations can be interpreted as the outcome of a giant collision that happened at an orbital radius of ~0.3 au sometime in 2014. Our analysis implies that a collision on a similar scale could have taken place around 2010, too. The fact that the disc was already peculiarly dust rich 40 yr ago, as implied by IRAS data, suggests that these dust production events belong to a chain of large impacts triggered by an earlier even more catastrophic collision., Comment: 18 pages, 12 figures, published in MNRAS

DI Cha A is K0-type pre-main sequence star, the brightest component of a quadruple stellar system. Here we report on a detailed study of this star based on archival VLTI/MIDI and VLTI/PIONIER infrared interferometric observations, as well as optical--infrared photometric monitoring from ground-based and space-born instruments. We determined the structure of the circumstellar disk by fitting simultaneously the interferometric visibilities and the spectral energy distribution, using both analytical models and the radiative transfer code RADMC-3D. The modeling revealed that the radial density distribution of the disk appears to have a gap between 0.21 and 3.0 au. The inner ring, whose inner size coincides with the sublimation radius, is devoid of small, submicrometer-sized dust grains. The inner edge of the outer disk features a puffed-up rim, typically seen in intermediate-mass stars. Grain growth, although less progressed, was also detected in the outer disk. The inner ring is variable at mid-infrared wavelengths on both daily and annual timescales, while the star stays remarkably constant in the optical, pointing to geometrical or accretion changes in the disk as possible explanation for the flux variations., 19 pages, 18 figures

The satellite of (225088) 2007 OR10 was discovered on archival Hubble Space Telescope images and along with new observations with the WFC3 camera in late 2017 we have been able to determine the orbit. The orbit's notable eccentricity, e$\approx$0.3, may be a consequence of an intrinsically eccentric orbit and slow tidal evolution, but may also be caused by the Kozai mechanism. Dynamical considerations also suggest that the moon is small, D$_{eff}$ $, Comment: Accepted for publication in Icarus

Pronounced global cooling around the Eocene–Oligocene transition (EOT) was a pivotal event in Earth’s climate history, controversially associated with the opening of the Drake Passage. Using a physical laboratory model we revisit the fluid dynamics of this marked reorganization of ocean circulation. Here we show, seemingly contradicting paleoclimate records, that in our experiments opening the pathway yields higher values of mean water surface temperature than the “closed” configuration. This mismatch points to the importance of the role ice albedo feedback plays in the investigated EOT-like transition, a component that is not captured in the laboratory model. Our conclusion is supported by numerical simulations performed in a global climate model (GCM) of intermediate complexity, where both “closed” and “open” configurations were explored, with and without active sea ice dynamics. The GCM results indicate that sea surface temperatures would change in the opposite direction following an opening event in the two sea ice dynamics settings, and the results are therefore consistent both with the laboratory experiment (slight warming after opening) and the paleoclimatic data (pronounced cooling after opening). It follows that in the hypothetical case of an initially ice-free Antarctica the continent could have become even warmer after the opening, a scenario not indicated by paleotemperature reconstructions.

V899 Mon is an eruptive young star showing characteristics of both FUors and EXors. It reached a peak brightness in 2010, then briefly faded in 2011, followed by a second outburst. We conducted multi-filter optical photometric monitoring, as well as optical and near-infrared spectroscopic observations of V899 Mon. The light curves and color-magnitude diagrams show that V899 Mon has been gradually fading after its second outburst peak in 2018, but smaller accretion bursts are still happening. Our spectroscopic observations taken with Gemini/IGRINS and VLT/MUSE show a number of emission lines, unlike during the outbursting stage. We used the emission line fluxes to estimate the accretion rate and found that it has significantly decreased compared to the outbursting stage. The mass loss rate is also weakening. Our 2D spectro-astrometric analysis of emission lines recovered jet and disk emission of V899 Mon. We found the emission from permitted metallic lines and the CO bandheads can be modeled well with a disk in Keplerian rotation, which also gives a tight constraint for the dynamical stellar mass of 2 ${M_{\odot}}$. After a discussion of the physical changes that led to the changes in the observed properties of V899 Mon, we suggest this object is finishing its second outburst., 31 pages, 26 figures, accepted for publication in ApJ

We describe the target selection procedure by which stars are selected for 2-minute and 20-second observations by TESS. We first list the technical requirements of the TESS instrument and ground systems processing that limit the total number of target slots. We then describe algorithms used by the TESS Payload Operation Center (POC) to merge candidate targets requested by the various TESS mission elements (the Target Selection Working Group, TESS Asteroseismic Science Consortium, and Guest Investigator office). Lastly, we summarize the properties of the observed TESS targets over the two-year primary TESS mission. We find that the POC target selection algorithm results in 2.1 to 3.4 times as many observed targets as target slots allocated for each mission element. We also find that the sky distribution of observed targets is different from the sky distributions of candidate targets due to technical constraints that require a relatively even distribution of targets across the TESS fields of view. We caution researchers exploring statistical analyses of TESS planet-host stars that the population of observed targets cannot be characterized by any simple set of criteria applied to the properties of the input Candidate Target Lists., 15 pages, 6 figures, accepted for publication in PASP

Light curves of young star systems show photometric variability due to different kinematic and physical processes. One of the main contributors to the photometric variability is the changing mass accretion rate, which regulates the interplay between the forming young star and the protoplanetary disk. We collected high-resolution spectroscopy in eight different epochs, as well as ground-based and space-borne multiepoch optical and infrared photometry of WX Cha, an M0 binary system, with an almost edge-on disk (i = 87°) in the Chamaeleon I star-forming region. Spectroscopic observations cover 72 days, the ground-based optical monitoring covers 42 days while space-borne TESS photometry extends for 56 days. The multiwavelength light curves exhibit quasi-periodic variability of 0.35–0.53 mag in the near-infrared, and of 1.3 mag in the g band. We studied the variability of selected emission lines that trace the accretion, computed the accretion luminosity and the mass accretion rate using empirical relations, and obtained values between L acc ∼ 1.6 L ⊙ − 3.2 L ⊙ and M ̇ acc ∼ 3.31 × 10 − 7 M ⊙ yr − 1 − 7.76 × 10 − 7 M ⊙ yr−1. Our results show that WX Cha is accreting at a rate larger than what is typical for T Tauri stars in the same star-forming region with the same stellar parameters. We theorize that this is due to the higher disk mass of WX Cha than what is usual for stars with similar stellar mass and to the binary nature of the system. Daily changes in the accretion luminosity and in the extinction can explain the photometric variability.

There are thousands of previously known and newly identified eclipsing binaries observed almost continuously for nearly a year long during Year 2 of the TESS mission. These precise and unique data sets give us an opportunity to measure a significant amount of their eclipse times on a time-scale comparable to the usual outer orbital periods of compact hierarchical triple stellar systems (below 1-2 years). In our project, we obtain the light curves from TESS Full-Frame Images of more than 4,000 eclipsing binaries observed at least during 8 Sectors and determine their eclipse timing variations in order to search for the signals of third-body perturbations and/or light travel time effect caused by an unseen, more distant component. We present orbital solutions for the identified triple systems and statistical considerations will also be discussed., {"references":["Borkovits T., Hajdu T., Sztakovics J. et al., 2016, MNRAS, 455, 4136","Hippke M., David T. J., Mulders G. D., Heller R., 2019, AJ, 158, 143","Pál A., 2012, MNRAS, 421, 1825"]}

The Herschel Space Observatory had two imaging instruments, working in the far-infrared and submillimetre regimes: the PACS cameras at 70/100 and 160 μm and the SPIRE photometers at 250, 350 and 500 μm. Small solar system bodies, especially main belt asteroids were serendipitously present in the field of view mainly in large scan maps. We identified these objects with the original aim to mark the affected sources in the Herschel PACS and SPIRE Point Source Catalogues. In our present study we are extracting flux densities in the PACS bands for asteroids above the detection limit, either using existing standard level 2 data products from the Herschel Science Archive, or re-reducing the PACS maps in the co-moving frame of the target. We are performing aperture photometry on the PACS maps with the HIPE software, because it has all the necessary built-in functions and values for the aperture correction. To obtain high S/N flux densities we are using somewhat smaller aperture radii, namely 4, 5 and 8 arcsecs in blue, green and red bands, see Figure 1. For the faintest objects, where the S/N is below 10, we are using larger apertures of 6, 7, and 10 arcsec in radius, to get a higher encircled energy fraction. We are making a 5' cutout around the asteroid, which is being stored as a png file to allow us to make visual inspections of the source and its surroundings. Using HIPE's annularSkyAperturePhotometry function we are performing the photometry, using the above mentioned aperture sizes and the built-in centroid algorithm. The initial coordinates are derived from RA and DEC values from JPL/Horizons. On the raw photometry results we perform the aperture correction with the photApertureCorrectionPointSource task, taking into account the band, the observation mode (parallel or not) and the scan speed. The centered aperture is also plotted on the cutout. One of the most important and difficult task is to give a proper error estimation for the extracted fluxes. First, we are using the photometry error given by the HIPE tasks, and the coverage, to determine how much time the asteroid spent on the detector and assess the quality of the photometry. Finally we are adding 5 percent absolute flux calibration error. Figure 1. Left: (16) Psyche in red band. The aperture for photometry, and the annulus for background estimation is plotted on the image, as well as the apertures for an alternative error estimation. Right: Four serendipitously observed asteroids on blue maps with no shape model and very limited or no thermal measurements. Left upper: (475) Ocllo, Left bottom: (11014) Svatopluk, Right upper: (634) Ute, Right bottom: (13029) 1989 HA With this method we extracted 633 flux densities of 275 asteroids, on 272 maps, 229 in the blue band, 179 in the green and 225 in the red band.Five of the serendipitously observed asteroids are mission targets, and are well known, e.g. (1) Ceres, (4) Vesta. 59 of them have shape model and sufficient multi-mission thermal measurements, e.g. (13) Egeria, (43) Ariadne, 12 have shape model, but very limited or no thermal measurements, e.g. (212) Medea, or (675) Ludmilla. In these cases the new flux densities can help confirming the existing shape model and to develop thermophysical models for the asteroids. 48 have sufficient multi-mission thermal measurements, but no shape models, e.g. (58) Concordia or (128) Nemesis. Finally, 19 asteroid from our work don't have shape model and has very limited or no thermal measurements, e.g. (475) Ocllo or (634) Ute, see Figure 1. We are planning to publish the new fluxes and these new flux densities will be included in the Small Bodies: Near and Far (SBNAF) Infrared Database (Szakáts et al., 2020). The fluxes obtained from Herschel are excellent for radiometric studies to get the object's size, albedo and maybe also thermal properties, when combined with other measurements (Alí-Lagoa et al., 2020). A natural continuation of our current work is the extension to SPIRE maps to find serendipitous asteroids, extract new submm fluxes and flag sources in the SPIRE Point Source Catalog for possible contamination. For those asteroids which have thermal emission data available at shorter wavelengths the Herschel PACS and SPIRE measurements allow us to determine the far infrared and submmm emissivities.

Despite that they are the brightest representatives of trans-Neptunian objects, many dwarf planets still miss a well-determined light curve and hence rotation period, in many cases due to the biases and aliases of ground-based observations. As it was previously shown in the case of e.g. Gonggong (2007 OR10), one of the largest Kuiper belt objects known, a well determined rotation period can be deduced from light curve data obtained by space surveys with original goals very different from the observations of Solar system objects (Pál et al., 2016; Kiss et al., 2017). In the case of Gonggong the observations were taken in the framework of the K2 mission of the Kepler Space Telescope (Howell et al., 2014) which was initially dedicated to measure exoplanet transits. However, in its K2 mission the telescope was pointed to the ecliptic, and allowed the long, in some cases up to ~80 day measurements of small solar system bodies, resulting in a previously undetected overabundance of long rotation periods in many of the small body populations, including the Hildas, Jovian Trojans, Centaurs, as well as trans-Neptunian objects (Kalup et al., 2021; Kiss et al., 2020; Marton et al., 2020; Szabó et al., 2017, 2020). The Transiting Exoplanet Survey Satellite (TESS, Ricker et al., 2015) was also designed originally to find exoplanets and avoid the vicinity of the ecliptic. However, it was still able to observe many asteroid at intermediate ecliptic latitudes. While the maximum length of the TESS observations are shorter than those in the K2 mission (28 days vs. 80 days) these are still much longer than that a typical long observing campaign can provide from the ground. Recently, a large number of main belt asteroid light curves was obtained by TESS, leading to the discovery of a previously undetected group of slow rotators (Pál et al., 2020). Due to the smaller light collecting area most trans-Neptunian objects are too faint to be observed with TESS, but the largest ones -- the dwarf planets -- are bright enough to provide a meaningful light curve. In this presentation we show results from the analysis of TESS light curves of the Kuiper belt dwarf planets Eris, Makemake and Orcus. These data are, in some cases, combined with ground-based data, and also with thermal emission measurements to constrain the spin axis orientation and determine how the light curve can be caused by surface features. As all the three targets have companions, the rotational information can also be incorporated into tidal evolution models and can be used to constrain the formation scenarios of the system, and, in some cases, give hints on the internal structure and material properties. Figure 1: Right panels: The apparent path of Eris during the first year of the TESS observations just in between the boundaries of TESS Sectors #3 (red) and #4 (blue). During its retrograde motion, Eris just misses Camera #1, CCD #3 of Sector #3. The thin line shows the full path of Eris while the thick red and blue sections correspond to the observing runs of Sectors #3 and #4, respectively. Left panels: The same plot for Orcus and for TESS Sectors #8 (red) and #9 (blue). Orcus has been observed during Sector #9 by CCD #1 of Camera #1, pretty close to the CCD edges. Bottom panel: The apparent track of Orcus within the field of view of TESS Camera #1 CCD #1, section [44:144]x[1675:1705]. The semi-transparent green curve shows the full track where observations were available while the thin red curve shows the track from which data points were involved in the photometric analysis. The direction of this section with respect to the sky is mirrored (east is nearly to the right and north points nearly to up): Orcus moved from the left to right due to its retrograde motion at its opposition. References: - Howell, S.B., et al., 2014, PASP, 126, 398- Kalup, C., et al., 2021, ApJS, 254, 7- Kiss C., et al., 2017, ApJL, 838, L1- Kiss, C., et al., 2020, EPSC2020-536 - Marton, G., et al., 2020, Icarus, 345, 113721- Pál, A., et al., 2020, ApJS, 247, 26- Ricker, G.R., et al. 2015, J. Astron. Telesc. Instrum. Syst. Vol. 1,id. 014003- Szabó et al., 2017, A&A, 599, 44- Szabó et al., 2020, ApJS, 247, 26

To study astrophysical transit phenomena we follow an alternative strategy for getting high-cadence observations of the field. This can be achieved with our new Fly’s Eye Camera System that monitors the entire sky above 30◦ horizontal altitude. With this instrument one can observe all phenomena brighter than ∼ 15m in Sloan r-band (u’, g’,i’ and z’ filters are also available). If we stack together a few hour of images we canobserve ∼ 17 m faint sources. This small-sized instrument is designed for time-domain astronomy with its 150 sec cadence. Due to the hexapod-based motion control, the instrument can be installed anywhere without any modifications, it can accomplish sky tracking automatically. These parallel robots have 6 degrees of freedom (DoF), but since any kind of rotation can be done by using only 3 DoF, the tracking with hexapods is independent from the geographical coordinates. Even polar alignment is not required, because Fly’s Eye can calibrate itself based on its own observed data. The system is optimal for time-domain astronomy: detecting novae, supernovae, optical afterglows of gamma-ray bursts and other bright, fast transients, and, from the observation database such data can be obtained - even from before the discovery of the transient event. In the future when the direction of the gravitational waves will be defined precisely we will be able to detect their first multiwaveband counterparts. In addition the Fly’s Eye will support the “Transient Astrophysical Objects” project which will use two new 80 cm robotic telescopes for follow-up observations of transients.

We are designing a fleet of nanosatellites to perform accurate position determinations of short-duration gamma-ray bursts by measuring arrival time differences. To achieve sufficient photon statistics to measure the arrival times precisely under the severe limitations of size, mass and power consumption, we propose the use of large-area CsI(Tl) scintillators which provide a high light output and the use of a small-sized multipixel photon counters (MPPC), a kind of SiPM, that have low power consumption. We plan to use one of the latest-model MPPCs provided by Hamamatsu Photonics, namely, S13360-6050CS, which has an active area of 6 x 6 mm(2). We have compared the performance of two scintillators of different sizes (150x75x5 mm(3) and 100x75x5 mm(3)); the bigger one is the maximum size that can be mounted on a three-unit satellite, according to CubeSat standards. We have found that the light yield of the bigger scintillator is lower than that of the smaller one, but the difference is only approximately 13%, and each has an energy threshold of similar to 10 keV at 25 degrees C. We have tested one-and two-MPPC readout and confirmed the same energy thresholds. The light yield for the two-MPPC readout is 1.4 times as high as that for the one-MPPC readout. We have also examined the position dependence of the light yield by using radiation from a Am-241 (59.5 keV) source, and confirmed that uniformity was improved, and energy resolution got better by similar to 7% for two-MPPC readout.

We report on a laboratory model of the large-scale flow phenomena in the Southern Ocean with a closed Drake Passage, imitating the situation before the Eocene-Oligocene transition (EOT) ca. 34 million years ago. In a differentially heated rotating annular wave tank an insulating ‘meridional’ barrier is installed to block zonal circulation. The obtained temperature fields and time series are compared to the ones from control runs with partially blocked and fully opened ‘passages’. In the ‘closed’ case a persistent azimuthal temperature gradient emerges whose magnitude scales linearly with the ‘meridional’ temperature contrast. The anomalous temperature distribution is accompanied with perturbations of the background flow yielding significantly larger low frequency variability than in the ‘opened’ configuration. The experimental findings have implications to the EOT but at present validation of results is hindered by the lack of deep-sea drilling data from the Southeast Pacific and adjacent parts of the Southern Ocean. Our model may inspire further paleoceanographic research in this largely unexplored region.

During World War 2 Szeged served as the key cultural, economic, military and trade centre of the Southern Great Plain in Hungary, also being an important water transport and railway junction. No wonder that Allied intelligence services took notice of the town early on, starting to collect reconnaissance information about Szeged’s most important establishments and objects in order to prepare for their destruction. In my study I dissect the town’s Second World War military geography, and look into how strategic targets in Szeged were located.

Primarily thanks to recently discovered American archives documents, the comprehensive history of the American air raids against Hungary during the Second World War is now ready for review. On the 75th anniversary of the events I briefly summarise and analyse the in-depth data of the American air raids against Budapest, and I investigate whether, in comparison to rural cities and considering its significance, Budapest was an overrepresented target during the 1944–1945 air war. Based on my conclusions and in view of the whole picture, there is no doubt about Budapest having been more targeted and overrepresented than the rural cities. One obvious reason for that was the centralised industrial transporting lines and railway. It is certain that Budapest was the most heavily bombed location in Hungary by the U.S. Air Force during World War 2, at least the statistics support this fact. However, if upon investigating individual air raids we examine the target groups and further break them down to targets, we find that one by one these had to withstand almost the same strain as a rural city.

In recent years the number of CubeSats (U-class spacecrafts) launched into space has increased exponentially marking the dawn of the nanosatellite technology. In general these satellites have a much smaller mass budget compared to conventional scientific satellites which limits shielding of scientific instruments against direct and indirect radiation in space. In this paper we present a simulation framework to quantify the signal in large field-of-view gamma-ray scintillation detectors of satellites induced by X-ray/gamma-ray transients, by taking into account the response of the detector. Furthermore, we quantify the signal induced by X-ray and particle background sources at a Low-Earth Orbit outside South Atlantic Anomaly and polar regions. Finally, we calculate the signal-to-noise ratio taking into account different energy threshold levels. Our simulation can be used to optimize material composition and predict detectability of various astrophysical sources by CubeSats. We apply the developed simulation to a satellite belonging to a planned CAMELOT CubeSat constellation. This project mainly aims to detect short and long gamma-ray bursts (GRBs) and as a secondary science objective, to detect soft gamma-ray repeaters (SGRs) and terrestrial gamma-ray flashes (TGFs). The simulation includes a detailed computer-aided design (CAD) model of the satellite to take into account the interaction of particles with the material of the satellite as accurately as possible. Results of our simulations predict that CubeSats can complement the large space observatories in high-energy astrophysics for observations of GRBs, SGRs and TGFs. For the detectors planned to be on board of the CAMELOT CubeSats the simulations show that detections with signal-to-noise ratio of at least 9 for median GRB and SGR fluxes are achievable., 86 pages, 30 figures, 8 tables, accepted for publication in Journal of Astronomical Telescopes, Instruments, and Systems

Pronounced global cooling around the Eocene-Oligocene transition (EOT) was a pivotal event in Earth’s climate history, controversially associated with the opening of the Drake Passage. Using a physical laboratory model we revisit the fluid dynamics of this marked reorganization of ocean circulation. Our differentially heated rotating annulus is a widely studied experimental set-up designed to model mid-latitude circulation in the atmosphere and the ocean, as well. Here we show, seemingly contradicting paleoclimate records that in our experiments opening the pathway yields higher values of mean water surface temperature than the “closed” configuration. This mismatch points to the importance of the crucial role ice albedo feedback plays in the investigated EOT-like transition, a component that is not captured in the laboratory model. Our conclusion is supported by numerical simulations performed in a global climate model (GCM) of intermediate complexity, where both “closed” and “open” configurations were explored, with and without active ice-albedo feedback. The GCM results indicate that sea surface temperatures would change in the opposite direction following an opening event in the two ice dynamics settings, and the results are therefore consistent both with the laboratory experiment (slight warming after opening) and the paleoclimatic data (pronounced cooling after opening). It follows that in the hypothetical case of an initially ice-free Antarctica the continent could have become even warmer after the opening, a scenario not indicated by paleotemperature reconstruction. These results provide circumstantial evidence supporting a particular EOT scenario in which Antarctica had already been – at least partially – covered with ice when the Drake Passage fully opened.

As more exoplanets are being discovered around ultracool dwarfs, understanding their magnetic activity -- and the implications for habitability -- is of prime importance. To find stellar flares and photometric signatures related to starspots, continuous monitoring is necessary, which can be achieved with spaceborn observatories like the Transiting Exoplanet Survey Satellite (TESS). We present an analysis of TRAPPIST-1 like ultracool dwarfs with TESS full-frame image photometry from the first two years of the primary mission. A volume-limited sample up to 50 pc is constructed consisting of 339 stars closer than 0.5 mag to TRAPPIST-1 on the Gaia colour-magnitude diagram. The 30-min cadence TESS light curves of 248 stars were analysed, searching for flares and rotational modulation caused by starspots. The composite flare frequency distribution of the 94 identified flares shows a power law index similar to TRAPPIST-1, and contains flares up to $E_\mathrm{TESS} = 3 \times 10^{33}$ erg. Rotational periods shorter than 5 days were determined for 42 stars, sampling the regime of fast rotators. The ages of 88 stars from the sample were estimated using kinematic information. A weak correlation between rotational period and age is observed, which is consistent with magnetic braking., 19 pages, 23 figures, accepted by A&A

Due to the advancement of nano-satellite technology, CubeSats and fleets of CubeSats can form an alternative to high-cost large-size satellite missions with the advantage of extended spatial coverage. One of these initiatives is the Cubesats Applied for MEasuring and LOcalising Transients (CAMELOT) mission concept, aimed at detecting and localizing gamma-ray bursts with an efficiency and accuracy comparable to large gamma-ray space observatories. While precise attitude control is not necessary for such a mission, attitude determination is an important issue in the interpretation of scintillator detector data as well as optimizing downlink telemetry. The employment of star trackers is not always a viable option for such small satellites, hence another alternative is necessary. A new method is proposed in this series of papers, utilizing thermal imaging sensors to provide simultaneous measurement of the attitude of the Sun and the horizon by employing a homogeneous array of such detectors. The combination with Sun and horizon detection w.r.t. the spacecraft would allow the full 3-DoF recovery of its attitude. In this paper we determine the spherical projection function of the MLX90640 infrasensors planned to be used for this purpose. We apply a polynomial transformation with radial corrections to map the spatial coordinates to the sensor plane. With the determined projection function the location of an infrared point source can be determined with an accuracy of $\sim 40^{\prime }$ , well below the design goals of a nano-satellite designed for gamma-ray detection.

Extreme debris disks (EDDs) are rare systems with peculiarly large amounts of warm dust that may stem from recent giant impacts between planetary embryos during the final phases of terrestrial planet growth. Here we report on the identification and characterization of six new EDDs. These disks surround F5-G9 type main-sequence stars with ages >100 Myr, have dust temperatures higher than 300K and fractional luminosities between 0.01 and 0.07. Using time-domain photometric data at 3.4 and 4.6$��$m from the WISE all sky surveys, we conclude that four of these disks exhibited variable mid-infrared emission between 2010 and 2019. Analyzing the sample of all known EDDs, now expanded to 17 objects, we find that 14 of them showed changes at 3-5$��$m over the past decade suggesting that mid-infrared variability is an inherent characteristic of EDDs. We also report that wide-orbit pairs are significantly more common in EDD systems than in the normal stellar population. While current models of rocky planet formation predict that the majority of giant collisions occur in the first 100 Myr, we find that the sample of EDDs is dominated by systems older than this age. This raises the possibility that the era of giant impacts may be longer than we think, or that some other mechanism(s) can also produce EDDs. We examine a scenario where the observed warm dust stems from the disruption and/or collisions of comets delivered from an outer reservoir into the inner regions, and explore what role the wide companions could play in this process., 39 pages, 11 figures, accepted for publication in ApJ

We report the discovery of a compact, coplanar, quadruply-lined, eclipsing quadruple star system from TESS data, TIC 454140642, also known as TYC 0074-01254-1. The target was first detected in Sector 5 with 30-min cadence in Full-Frame Images and then observed in Sector 32 with 2-min cadence. The light curve exhibits two sets of primary and secondary eclipses with periods of PA = 13.624 days (binary A) and PB = 10.393 days (binary B). Analysis of archival and follow-up data shows clear eclipse-timing variations and divergent radial velocities, indicating dynamical interactions between the two binaries and confirming that they form a gravitationally-bound quadruple system with a 2+2 hierarchy. The Aa+Ab binary, Ba+Bb binary, and A-B system are aligned with respect to each other within a fraction of a degree: the respective mutual orbital inclinations are 0.25 degrees (A vs B), 0.37 degrees (A vs A-B), and 0.47 degrees (B vs A-B). The A-B system has an orbital period of 432 days - the second shortest amongst confirmed quadruple systems - and an orbital eccentricity of 0.3., Comment: 25 pages, 17 figures, 9 tables; accepted for publication in the Astrophysical Journal

Various properties of Jovian trojan asteroids such as composition, rotation periods, and photometric amplitudes, or the rate of binarity in the population can provide information and constraints on the evolution of the group and of the Solar System itself. Here we present new photometric properties of 45 Jovian trojans from the K2 mission of the Kepler space telescope, and present phase-folded light curves for 44 targets, including (11351) Leucus, one of the targets of the Lucy mission. We extend our sample to 101 asteroids with previous K2 Trojan measurements, then compare their combined amplitude- and frequency distributions to other ground-based and space data. We show that there is a dichotomy in the periods of Trojans with a separation at $\sim 100$ hr. We find that 25% of the sample are slow rotators (P$\geq$30 hr), which excess can be attributed to binary objects. We also show that 32 systems can be classified as potential detached binary systems. Finally, we calculate density and rotation constraints for the asteroids. Both the spin barrier and fits to strengthless ellipsoid models indicate low densities and thus compositions similar to cometary and TNO populations throughout the sample. This supports the scenario of outer Solar System origin for Jovian trojans., Comment: 25 pages, 19 figures, accepted for publication in ApJS

We present our photometric and spectroscopic observations on the peculiar transient AT2018cow. The multi-band photometry covers from peak to $\sim$70 days and the spectroscopy ranges from 5 to $\sim$50 days. The rapid rise ($t_{\mathrm{r}}$$\lesssim$2.9 days), high luminosity ($M_{V,\mathrm{peak}}\sim-$20.8 mag) and fast decline after peak make AT2018cow stand out of any other optical transients. While we find that its light curves show high resemblance to those of type Ibn supernovae. Moreover, the spectral energy distribution remains high temperature of $\sim$14,000 K after $\sim$15 days since discovery. The spectra are featureless in the first 10 days, while some broad emission lines due to H, He, C and O emerge later, with velocity declining from $\sim$14,000 km s$^{-1}$ to $\sim$3000 km s$^{-1}$ at the end of our observations. Narrow and weak He I emission lines emerge in the spectra at $t>$20 days since discovery. These emission lines are reminiscent of the features seen in interacting supernovae like type Ibn and IIn subclasses. We fit the bolometric light curves with a model of circumstellar interaction (CSI) and radioactive decay (RD) of \Ni and find a good fit with ejecta mass $M_{\mathrm{ej}}\sim$3.16 M$_{\odot}$, circumstellar material mass $M_{\mathrm{CSM}}\sim$0.04 M$_{\odot}$, and ejected \Ni mass $M_{^{56}\mathrm{Ni}}\sim$0.23 M$_{\odot}$. The CSM shell might be formed in an eruptive mass ejection of the progenitor star. Furthermore, host environment of AT2018cow implies connection of AT2018cow with massive stars. Combining observational properties and the light curve fitting results, we conclude that AT2018cow might be a peculiar interacting supernova originated from a massive star., Accepted for publication in ApJ