Mitigation of the threat from airbursting asteroids requires an understanding of the potential risk they pose for the ground. How asteroids release their kinetic energy in the atmosphere is not well understood due to the rarity of significant impacts. Ordinary chondrites, in particular L chondrites, represent a frequent type of Earth-impacting asteroids. Here, we present the first comprehensive, space-to-lab characterization of an L chondrite impact. Small asteroid 2023 CX1 was detected in space and predicted to impact over Normandy, France, on 13 February 2023. Observations from multiple independent sensors and reduction techniques revealed an unusual but potentially high-risk fragmentation behavior. The nearly spherical 650 $\pm$ 160 kg (72 $\pm$ 6 cm diameter) asteroid catastrophically fragmented around 28 km altitude, releasing 98% of its total energy in a concentrated region of the atmosphere. The resulting shockwave was spherical, not cylindrical, and released more energy closer to the ground. This type of fragmentation increases the risk of significant damage at ground level. These results warrant consideration for a planetary defense strategy for cases where a >3-4 MPa dynamic pressure is expected, including planning for evacuation of areas beneath anticipated disruption locations.

We report the discovery and characterization of three transiting giant planets in the TIC118798035 system. The three planets were identified as transiting candidates from data of the TESS mission, and confirmed with ground-based photometric transit observations along with radial velocity variations obtained with FEROS, HARPS and ESPRESSO. The three planets present transit timing variations (TTVs). We performed a N-body orbital fitting to the TTVs and radial velocities finding that TIC118798035 b is as warm low-density Neptune with a mass of 0.0250$\pm$0.0023 $M_J$, a radius of 0.655$\pm$0.018 $R_J$, and an orbital period of 11.507 d; TIC118798035 c is a warm Saturn with a mass of 0.403$\pm$0.024 $M_J$, a radius of 0.973$\pm$0.023 $R_J$, and an orbital period of 22.564 d; and TIC118798035 d is a warm Jupiter with a mass of 0.773$\pm$0.052 $M_J$, a radius of 0.923$\pm$0.044 $R_J$, and an orbital period of 48.925 d. The bulk metallicities of the three planets don't fully follow the mass-metallicity correlation found for the giant planets of the solar system, which hints at a somewhat different formation history for the planets of the TIC118798035 system. TIC118798035 is the only system having more than two transiting planets larger than 0.5 $R_J$ with a precise orbital and physical characterization, amenable for future atmospheric studies.

The energy equilibrium between the corona and the underlying disk in a two-phase accretion flow sets a lower limit on the achievable photon index. A slab corona may not explain the hard state observations of X-ray binaries (XRBs). We incorporate energy feedback to the accretion disk resulting from illumination by an extended corona, and vice versa. The interaction between these two components allows for the possibility of finding an energetically self-consistent equilibrium solution for a given disk-corona system. We have upgraded the existing Monte Carlo radiative transfer code, MONK, to incorporate the interaction between the disk and the extended corona within the general relativistic framework. We introduce an albedo parameter to specify the fraction of the incident flux that is reflected by the disk, while the remainder is absorbed and added to the intrinsic dissipation. Reflection is modeled assuming a semi-infinite electron atmosphere. We find global equilibrium solutions by iterating interaction between disk and extended slab corona. A higher black hole spin, higher coronal temperature, and higher albedo all lead to harder spectra. For typical coronal temperatures and disk albedo, the lowest achievable photon index with a static slab corona fully covering the disk is approximately 1.7-1.8. With the upgraded version of MONK, we are now able to achieve global energy equilibrium for a given disk-corona system. This approach holds significant potential for constraining the coronal geometry using not only the observed flux but also polarization. A static slab does not appear to be a favorable coronal geometry for the hard state of XRBs, even when global energy balance is taken into account. In future work, we will explore truncated disk geometries and outflowing coronae as potential alternatives. (shortened)

The forthcoming space-based gravitational-wave observatory Laser Interferometer Space Antenna (LISA) should enable the detection of Extreme Mass Ratio Inspirals (EMRIs), in which a stellar-mass compact object gradually inspirals into a supermassive black hole while emitting gravitational waves. Modeling the waveforms of such systems is a challenging task, requiring precise computation of energy and angular momentum fluxes as well as proper treatment of orbital resonances, during which two fundamental orbital frequencies become commensurate. In this work, we perform a systematic comparison of fluxes derived from three approaches: the quadrupole formula, post-Newtonian approximations, and time-domain solutions of the Teukolsky equation. We show that quadrupole-based fluxes remain in good agreement with Teukolsky results across a broad range of orbital configurations, including perturbed orbits. Building on these insights, we explore the dynamical impact of resonance crossings within the adiabatic approximation. By introducing novel numerical methods, we reduce computational costs and uncover diverse resonance-crossing behaviors. These results contribute to the effort to understand theoretically and model adequately resonance crossings during an EMRI.

We report the discovery of X-ray polarization from the X-ray-bright filament. G0.13-0.11 in the Galactic center (GC) region. This filament features a bright, hard X-ray source that is most plausibly a pulsar wind nebula (PWN) and an extended and structured diffuse component. Combining the polarization signal from IXPE with the imaging/spectroscopic data from Chandra, we find that X-ray emission of G0.13-0.11 is highly polarized PD=$57(\pm18)$% in the 3-6 keV band, while the polarization angle is PA=$21^\circ(\pm9^\circ)$. This high degree of polarization proves the synchrotron origin of the X-ray emission from G0.13-0.11. In turn, the measured polarization angle implies that the X-ray emission is polarized approximately perpendicular to a sequence of nonthermal radio filaments that may be part of the GC Radio Arc. The magnetic field on the order of $100\,{\rm\mu G}$ appears to be preferentially ordered along the filaments. The above field strength is the fiducial value that makes our model self-consistent, while the other conclusions are largely model independent.

Based on the Gaia DR2 catalogue of hot subdwarf star candidates, we identified 1587 hot subdwarf stars with spectra in LAMOST DR7. We present atmospheric parameters for these stars by fitting the LAMOST spectra with {\sc Tlusty/Synspec} non-LTE synthetic spectra. Combining LAMOST radial velocities and Gaia Early Data Release 3 (EDR3) parallaxes and proper motions, we also present the Galactic space positions, velocity vectors, orbital parameters and the Galactic population memberships of the stars. With our He classification scheme, we identify four groups of He rich hot subdwarf stars in the $T_{\rm eff}-\log\,g$and$T_{\rm eff}-\log{(n{\rm He}/n{\rm H})}$ diagrams. We find two extreme He-rich groups ($e$He-1 and $e$He-2) for stars with $\log{(n{\rm He}/n{\rm H})}\geq0$and two intermediate He-rich groups ($i$He-1 and$i$He-2) for stars with -1\le\log{(n{\rm He}/n{\rm H})}&lt;0. We also find that over half of the stars in Group $e$He-1 are thick disk stars, while over half of the stars in Group $e$He-2 correspond to thin disk stars. The disk population fractions of Group $i$He-1 are between those of Group $e$He-1 and $e$He-2. Almost all stars in Group $i$He-2 belong to the thin disk. These differences indicate that the four groups probably have very different origins. Comparisons between hot subdwarf stars in the halo and in the Galactic globular cluster $\omega$ Cen show that only He-deficient stars with $-2.2\le\log{(n{\rm He}/n{\rm H})}<-1$have similar fractions. Hot subdwarfs with$\log{(n{\rm He}/n{\rm H})}\ge 0$in$\omega$ Cen have no counterparts in the thick disk and halo populations, but they appear in the thin disk.

We consider a spinning test particle around a rotating black hole and compare the Mathisson-Papapetrou-Dixon (MPD) formalism under the Tulczyjew-Dixon spin supplementary condition to the test-mass limit of the effective-one-body (EOB) Hamiltonian of [Phys. Rev. D.90, 044018(2014)], with enhanced spin-orbit sector. We focus on circular equatorial orbits: we first compare the constants of motion at their linear in secondary spin approximation and then we compute the gravitational-wave (GW) fluxes using a frequency domain Teukolsky equation solver. We find no difference between the EOB and MPD fluxes when the background spacetime is Schwarzschild, while the difference for a Kerr background is maximum for large, positive spins. Our work could be considered as a first step to improve the radiation reaction of the EOB model, in view of the needs of the next-generation of GW detectors.

The Laser Interferometer Space Antenna (LISA) mission, scheduled for launch in the mid-2030s, is a gravitational wave observatory in space designed to detect sources emitting in the millihertz band. LISA is an ESA flagship mission, currently entering the Phase B development phase. It is expected to help us improve our understanding about our Universe by measuring gravitational wave sources of different types, with some of the sources being at very high redshifts $z\sim 20$. On the 23rd of February 2022 we organized the 1$^\mathrm{st}$ {\it LISA in Greece Workshop}. This workshop aimed to inform the Greek scientific and tech industry community about the possibilities of participating in LISA science and LISA mission, with the support of the Hellenic Space Center (HSC). In this white paper, we summarize the outcome of the workshop, the most important aspect of it being the inclusion of $15$ Greek researchers to the LISA Consortium, raising our total number to $22$. At the same time, we present a road-map with the future steps and actions of the Greek Gravitational Wave community with respect to the future LISA mission.

The superb image quality, stability and sensitivity of the JWST permit deconvolution techniques to be pursued with a fidelity unavailable to ground-based observations. We present an assessment of several deconvolution approaches to improve image quality and mitigate effects of the complex JWST point spread function (PSF). The optimal deconvolution method is determined by using WebbPSF to simulate JWST's complex PSF and MIRISim to simulate multi-band JWST/Mid-Infrared Imager Module (MIRIM) observations of a toy model of an active galactic nucleus (AGN). Five different deconvolution algorithms are tested: (1) Kraken deconvolution, (2) Richardson-Lucy, (3) Adaptive Imaging Deconvolution Algorithm, (4) Sparse regularization with the Condat-V\~u algorithm, and (5) Iterative Wiener Filtering and Thresholding. We find that Kraken affords the greatest FWHM reduction of the nuclear source of our MIRISim observations for the toy AGN model while retaining good photometric integrity across all simulated wavebands. Applying Kraken to Galactic Activity, Torus, and Outflow Survey (GATOS) multi-band JWST/MIRIM observations of the Seyfert 2 galaxy NGC 5728, we find that the algorithm reduces the FWHM of the nuclear source by a factor of 1.6-2.2 across all five filters. Kraken images facilitate detection of a SE to NW $\thicksim$2".5 ($\thicksim$470 pc, PA $\simeq$115\deg) extended nuclear emission, especially in the longest wavelengths. We demonstrate that Kraken is a powerful tool to enhance faint features otherwise hidden in the complex JWST PSF.

The vertical (Lense-Thirring) precession of the innermost accretion flows has been discussed as a sensitive indicator of the rotational properties of neutron stars (NSs) and their equation of state because it vanishes for a non-rotating star. In this work, we apply the Hartle-Thorne spacetimes to study the frequencies of the precession for both geodesic and non-geodesic (fluid) flows. We build on previous findings on the effect of the NS quadrupole moment, which revealed the importance of the interplay between the relativistic and classical precession. Because of this interplay, the widely used Lense-Thirring metric, linear in the NS angular momentum, is insufficient to calculate the behaviour of the precession frequency across an astrophysically relevant range of NS angular momentum values. We find that even for a moderately oblate NSs, the dependencies of the precession frequency on the NS angular momentum at radii within the innermost accretion region have maxima that occur at relatively low values of the NS angular momentum. We conclude that very different groups of accreting NSs, slow and fast rotators, can display the same precession frequencies. This may explain the lack of evidence for a correlation between the frequencies of the observed low-frequency quasiperiodic oscillations and the NS spin. In our work, we provide a full, general description of precession behaviour, and also examples that assume specific NS and quark star (MIT bag) equation of state. Our calculations are reproducible using the associated Wolfram Mathematica notebook.

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

We develop a fully analytical waveform model for precessing binaries with arbitrary spin vectors using post-Newtonian~(PN) theory in the extreme mass-ratio limit and a hierarchical multi-scale analysis. The analytical model incorporates leading PN order spin precession dynamics from spin-orbit, spin-spin, and quadrupole-monopole couplings, and 2PN order dissipative dynamics truncated to first order in the mass ratio $q \ll 1$. Due to the pure analytic nature of the model, the framework developed herein can readily be extended to both higher PN and higher-$q$ order. Although the PN series is asymptotic to this limit, our results can be used to estimate how precession affects the measurability of certain binary parameters, and to inform and compare with other waveform approximants, such as effective-one-body models, hybrid waveforms, and self-force calculations.

We present a study of the high-temperature spectral component in meteor fireballs, with a particular focus on neutral hydrogen at 656.28 nm and ionised silicon doublet at 634.71 nm and 637.14 nm. By analysing spectra from the European Fireball Network (EN) that exhibit H$\alpha$ and Si~II emissions, we investigated the relationship between H and Si abundances across different meteoroid types. The plasma temperature of the high-temperature component remains independent of meteor velocity. This allows us to directly compare relative intensities of volatile hydrogen with less volatile silicon in bodies with different velocities. Our results confirmed that the H/Si value remains largely independent of meteor velocity. We show a positive correlation with photometric mass for cometary meteoroids, suggesting that larger bodies better preserve their volatile content, namely hydrogen. This correlation persists across the meteor showers, showing a physical process related to volatile preservation rather than specific parent body composition. Our data suggest that the abundance of hydrogen in large cometary meteoroids is not only higher than in CI chondrites, but is also comparable to or higher than the measured abundances in small particles of dust from Halley's comet, depending on the assumed plasma conditions. This work brought new constraints on the distribution and preservation of volatile elements in Solar System bodies and new insights into the potential delivery mechanisms of water to Earth. The prevalence of hydrogen in larger cometary meteoroids supports models where comets could be significant contributors to Earth's volatile inventory.

We report on the Atmospheric Imaging Assembly (AIA) observations of plasma outflows originating in a coronal dimming during the 2015 April 28th filament eruption. After the filament started to erupt, two flare ribbons formed, one of which had a well-visible hook enclosing a core (twin) dimming region. Along multiple funnels located in this dimming, a motion of plasma directed outwards started to be visible in the 171 and 193 filter channels of the instrument. In time-distance diagrams, this motion generated a strip-like pattern, which lasted for more than five hours and which characteristics did not change along the funnel. We therefore suggest the motion to be a signature of outflows corresponding to velocities ranging between $\approx70$ and 140 km s$^{-1}$. Interestingly, the pattern of the outflows as well as their velocities were found to be similar to those we observed in a neighboring ordinary coronal hole. Therefore, the outflows were most likely a signature of a CME-induced slow solar wind flowing along the open-field structures rooted in the dimming region. Further, the evolution of the hook encircling the dimming region was examined in the context of the latest predictions imposed for the three-dimensional magnetic reconnection. The observations indicate that the filament's footpoints were, during their transformation to the dimming region, reconnecting with surrounding canopies. To our knowledge, our observations present the first imaging evidence for outflows of plasma from a dimming region.

This work provides gravitational wave energy and angular momentum asymptotic fluxes from a spinning body moving on generic orbits in a Kerr spacetime up to linear in spin approximation. To achieve this, we have developed a new frequency domain Teukolsky equation solver that calculates asymptotic amplitudes from generic orbits of spinning bodies with their spin aligned with the total orbital angular momentum. However, the energy and angular momentum fluxes from these orbits in the linear in spin approximation are appropriate for adiabatic models of extreme mass ratio inspirals even for spins non-aligned to the orbital angular momentum. To check the newly obtained fluxes, they were compared with already known frequency domain results for equatorial orbits and with results from a time domain Teukolsky equation solver called Teukode for off-equatorial orbits. The spinning body framework of our work is based on the Mathisson-Papapetrou-Dixon equations under the Tulczyjew-Dixon spin supplementary condition.

This handbook chapter gives a modern introduction to spectral analysis of celestial X-ray sources. Concepts presented include the instrumentation response, the linear modelling approximation, Poisson count statistics and the Gaussian approximation, data re-binning, visualisation techniques, handling backgrounds, Bayesian and frequentist viewpoints, uncertainty quantification of model parameters, model checking, model comparison, and inferring population distributions. Realistic hands-on example exercises with accompanying data files and code are included to apply the theoretical concepts in practice.

GRB 221009A is the brightest gamma-ray burst (GRB) observed to date. Extensive observations of its afterglow emission across the electromagnetic spectrum were performed, providing the first strong evidence of a jet with a nontrivial angular structure in a long GRB. We carried out an extensive observation campaign in very-high-energy (VHE) gamma rays with the first Large-Sized Telescope (LST-1) of the future Cherenkov Telescope Array Observatory (CTAO), starting on 2022 October 10, about one day after the burst. A dedicated analysis of the GRB 221009A data is performed to account for the different moonlight conditions under which data were recorded. We find an excess of gamma-like events with a statistical significance of 4.1$\sigma$ during the observations taken 1.33 days after the burst, followed by background-compatible results for the later days. The results are compared with various models of afterglows from structured jets that are consistent with the published multiwavelength data, but entail significant quantitative and qualitative differences in the VHE emission after one day. We disfavor models that imply VHE flux at one day considerably above $10^{-11}$ erg cm$^{-2}$ s$^{-1}$. Our late-time VHE observations can help disentangle the degeneracy among the models and provide valuable new insight into the structure of GRB jets.

The LAMOST DR1 survey contains about two million of spectra labelled by its pipeline as stellar objects of common spectral classes. There is, however, a lot of spectra corrupted in some way by both instrumental and processing artifacts, which may mimic spectral properties of interesting celestial objects, namely emission lines of Be stars and quasars. We have tested several clustering methods as well as outliers analysis on a sample of one hundred thousand spectra using Spark scripts running on Hadoop cluster consisting of twenty-four sixteen-core nodes. This experiment was motivated by an attempt to find rare objects with interesting spectra as outliers most dissimilar from all common spectra. The result of this time-consuming procedure is a list of several hundred candidates where different artifacts are prominent, but also tens of very interesting emission-line spectra requiring further detailed examination. Many of them may be quasars or even blazars as well as yet unknown Be-stars. It deserves mentioning that most of the work benefitted considerably from technologies of Virtual Observatory.

The width of a resonance in a nearly integrable system, i.e. in a non-integrable system where chaotic motion is still not prominent, can tell us how a perturbation parameter is driving the system away from integrability. Although the tool that we are presenting here can be used is quite generic and can be used in a variety of systems, our particular interest lies in binary compact object systems known as extreme mass ratio inspirals (EMRIs). In an EMRI a lighter compact object, like a black hole or a neutron star, inspirals into a supermassive black hole due to gravitational radiation reaction. During this inspiral the lighter object crosses resonances, which are still not very well modeled. Measuring the width of resonances in EMRI models allows us to estimate the importance of each perturbation parameter able to drive the system away from resonances and decide whether its impact should be included in EMRI waveform modeling or not. To tackle this issue in our study we show first that recurrence quantifiers of orbits carry imprints of resonant behavior, regardless of the system's dimensionality. As a next step, we apply a long short-term memory machine learning architecture to automate the resonance detection procedure. Our analysis is developed on a simple standard map and gradually we extend it to more complicated systems until finally we employ it in a generic deformed Kerr spacetime known in the literature as the Johannsen-Psaltis spacetime.

Creating a website to promote one's scientific work has become commonplace in many scientific disciplines. A plethora of options exist for framework to generate your website content, hosting it, and registering a domain name. The goal of this document is to provide early career scientists (1) an overview of the current options for creating a website to promote their professional persona, and (2) general advice concerning website written content and one's web presence. To get a sense of how other scientists created their websites, I created a survey asking colleagues about the services they used to create their websites and advice they have for someone creating a website. I received 54 responses from 53 astronomers and one computer scientist of which 23 were in an academic position beyond postdoc (faculty, scientist, etc.), 1 was an individual research fellow, 4 were in their third postdoc, 4 were in their second postdoc, 16 were in their first postdoc, and 6 were graduate students. I report the results of this survey here.