Papers for Thursday, Sep 19 2024

We show that some holographic dark energy models can lead to a future evolution of the universe in which the scale factor $a$ is asymptotically constant, while $\dot a \rightarrow 0$ and the corresponding energy and pressure densities also vanish. We provide specific examples of such models and general conditions that can lead to an asymptotically static universe, which we have called the ``long freeze." We show that in some cases, such evolution can follow an arbitrarily long exponential expansion essentially identical to the asymptotic evolution of $\Lambda$CDM.

According to the lambda CDM scenario, galaxies are formed through the hierarchical accretion of building blocks. Our Galaxy is a privileged place to look for the remnants of accretion events through the study of the chemical and kinematic properties of its halo stellar populations. Due to its low density, the stellar halo holds the most favorable conditions for chemical tagging. However, chemical tagging alone often yields weak results due to both uncertainties in chemical abundances and to overlapping chemical properties among different populations. To overcome this problem, the use of chemical and kinematic properties can be combined. In this Thesis, we developed a machine learning algorithm, named the CREEK, which combines orbital and chemical properties of halo stars observed by two large public spectroscopic surveys, Gaia-ESO and APOGEE. The CREEK operates as follows: 1)Data selection: We selected halo stars from the APOGEE and Gaia-ESO surveys based both on their velocity and metallicity and we computed their orbital parameters. 2)Using kinematics: The selected data were passed to a Siamese Neural Network that established links between stars based on their kinematic similarities. 3)Using chemistry: The graph was passed through a Graph Neural Network (GNN) auto-encoder that took as input the selected abundances. The abundances were chosen to maximize homogeneity within stars from the same cluster while ensuring distinctiveness between stars from different clusters. Additionally, we prioritised elements with smallest errors. The GNN auto-encoder computed a mean of the abundances of all connected stars, weighted on the number of links of each star and mapped the chemical space into a more efficient representation in the latent space. 4)Recovering structures: Finally, OPTICS was applied to the latent space, providing groups based on the chemical similarities of the stars.

Soon after the release of the WISE all-sky catalogue of 500 million mid-infrared (IR) objects, suggestions were made that it could be used to search for extrasolar devices constructed by an advanced civilization to convert a significant fraction of their host star's luminosity into useful work: "technostructures", "megastructures" or "Dyson spheres/structures", hereafter DSMs, whose inevitable waste heat would be seen by WISE at mid-IR wavelengths. However, a trawl of several million potentially-habitable Gaia-detected stars for mid-IR-excess signatures is fraught with danger, due to both noise from such a large sample and, more importantly, confusion with the emission from dusty background galaxies. In light of a recent claim of seven potential DSMs in MNRAS, a brief rebuttal appeared on arXiv. Further to this response, the relevance of WISE-detected galaxies is discussed in more detail, leading to a seemingly tight limit on the number and lifetime of DSMs, and indeed intelligent worlds, in the ~600-pc-radius region patrolled by Gaia. However, the detectability of DSMs is questioned: a DSM might extinguish its star at optical/near-IR wavelengths, and thus either not appear or appear anomalously faint in a stellar catalogue. Moreover, a civilization advanced enough to construct a DSM is likely to be advanced enough to use countermeasures to mask its presence from us.

Autonomous precision navigation to land onto the Moon relies on vision sensors. Computer vision algorithms are designed, trained and tested using synthetic simulations. High quality terrain models have been produced by Moon orbiters developed by several nations, with resolutions ranging from tens or hundreds of meters globally down to few meters locally. The SurRender software is a powerful simulator able to exploit the full potential of these datasets in raytracing. New interfaces include tools to fuse multi-resolution DEMs and procedural texture generation. A global model of the Moon at 20m resolution was integrated representing several terabytes of data which SurRender can render continuously and in real-time. This simulator will be a precious asset for the development of future missions.

We identify a point-symmetric morphology of three pairs of ears/clumps in the core-collapse supernova (CCSN) remnant (CCSNR) Puppis A, supporting the jittering jets explosion mechanism (JJEM). In the JJEM, the three pairs of jets that shaped the three pairs of ears/clumps in Puppis A are part of a large, about 10 to 30 pairs of jets that exploded Puppis A. Some similarities in morphological features between CCSNR Puppis A and three multipolar planetary nebulae considered to have been shaped by jets solidify the claim for shaping by jets. Puppis A has a prominent dipole structure, where one side is bright with a well-defined boundary, while the other is faint and defused. The neutron star (NS) remnant of Puppis A has a proper velocity, its natal kick velocity, in the opposite direction to the denser part of the dipole structure. We propose a new mechanism in the frame of the JJEM that imparts a natal kick to the NS, the kick-by-early asymmetrical pair (kick-BEAP) mechanism. At the early phase of the explosion process, the NS launches a pair of jets where one jet is much more energetic than the counter jet. The more energetic jet compresses a dense side to the CCSNR, and, by momentum conservation, the NS recoils in the opposite direction. Our study supports the JJEM as the primary explosion mechanism of CCSNe and enriches this explosion mechanism by introducing the novel kick-BEAP mechanism.

We analyze the dark matter (DM) halos of a sample of dwarf Ellitpicals (dE) and discuss cosmological and evolutionary implications. Using orbit modeling we recover their density slopes and, for the first time, the halo flattening. We find the `cusp-core' tension is mild, on average dEs have central slopes slightly below the Navarro Frenk White (NFW) predictions. However, the measured flattenings are still more spherical than cosmological simulations predict. Unlike brighter ETGs the total density slopes of dEs are shallower, and their average DM density does not follow their scaling relation with luminosity. Conversely, dE halos are denser and the densities steeper than in LTGs. We find average DM density and slope are strongly correlated with the environment and moderately with the angular momentum. Central, non-rotating dEs have dense and cuspy halos, whereas rotating dEs in Virgo's outskirts are more cored and less dense. This can be explained by a delayed formation of the dEs in the cluster outskirts, or alternatively, by the accumulated baryonic feedback the dEs in the outskirts have experienced during their very different star formation history. Our results suggest halo profiles are not universal (they depend on assembly conditions) and they evolve only mildly due to internal feedback. We conclude dEs in the local Universe have assembled at a higher redshift than local spirals. In these extreme conditions (e.g. star-formation, halo assembly) were very different, suggesting no new dEs are formed at present.

Rocky planets are thought to form with a magma ocean that quickly solidifies. The horizontal and vertical extent of this magma ocean depends on the interior thermal evolution of the planet, and possibly exogenous processes such as planet migration. We present a model for simulating the thermal history of tidally locked lava planets. We initiate the model with a completely molten mantle and evolve it for ten billion years. We adopt a fixed surface temperature of 3000 K for the irradiated day-side, but allow the night-side temperature to evolve along with the underlying layers. We simulate planets of radius 1.0$R_{\oplus}$ and 1.5$R_{\oplus}$ with different core mass fractions, although the latter does not significantly impact the thermal evolution. We confirm that the day-side magma ocean on these planets has a depth that depends on the planetary radius. The night-side, on the other hand, begins crystallizing within a few thousand years and completely solidifies within 800 million years in the absence of substantial tidal heating or day-night heat transport. We show that a magma ocean can be sustained on the night-side of a lava planet if at least 20 per cent of absorbed stellar power is transmitted from the day-side to the night-side via magma currents. Such day-night transport could be sustained if the magma has a viscosity of $10^{-3}$ Pa s, which is plausible at these temperatures. Alternatively, the night-side could remain molten if the mush layer is tidally heated at the rate of $8 \times 10^{-4}$ W/kg of mush, which is plausible for orbital eccentricities of $e > 7 \times 10^{-3}$. Night-side cooling is a runaway process, however: the magma becomes more viscous and the mush solidifies, reducing both day-night heat transport and tidal heating. Measurements of the night-sides of lava planets are therefore a sensitive probe of the thermal history of these planets.

We performed spectroscopic analyses of five local compact star-forming galaxies (CSFGs) with extremely high [OIII]/[OII] (O$_{32}$) ratios ($>20$). These targets remarkably share similar properties with high-redshift CIV emitters at $z>6$: high H$\beta$ equivalent widths (EWs $>200$Å), extreme O$_{32}$ ratios, low metallicities (12+log(O/H) $\lesssim7.8$), low C/O abundances (log(C/O) $<-0.6$), and high ionization conditions (log$U>-2$). The UV spectra were acquired using the Hubble Space Telescope's (HST) Cosmic Origins Spectrograph (COS) and Space Telescope Imaging Spectrograph (STIS). We have identified a wealth of rest-frame UV emission lines (CIV, HeII, OIII], CIII]) in the HST spectra. Notably, all our targets show intense CIV emission lines with rest-frame EWs $>10$Å, indicative of hard ionizing radiation. The rest-frame UV emission line diagnostics disfavor an AGN and could be consistent with significant shock contributions to the source of ionizing radiation. Four of our targets show high CIV/CIII] ratios ($\geq1.4$), suggestive of strong Lyman-continuum leakage (LyC escape fraction, $f_{\rm esc,LyC}>10$%) from these sources. This is consistent with their Ly$\alpha$-inferred LyC escape fractions ($f_{\rm esc,LyC}=$ 9 - 31%). We derive relative C/O abundances from our sources, showing log(C/O) values from $-1.12$ to $-0.61$, comparable to those of reionization-era galaxies at $z\gtrsim6$. The properties of the CSFGs, particularly their intense CIV emission and high O$_{32}$ ratios, which suggest significant LyC escape fractions, are similar to those of the reionization-era CIV emitters. These similarities reinforce the hypothesis that these CSFGs are the closest analogs of significant contributors to the reionization of the intergalactic medium.

We use the latest dataset of supernova (SN) host galaxies to investigate how the host properties -- stellar mass, star formation rate, metallicity, absolute magnitude, and colour -- differ across SN types, with redshift-driven selection effects controlled. SN Ib and Ic host galaxies, on average, are more massive, metal-rich, and redder than SN II hosts. For subtypes, SN Ibn and Ic-BL have bluer hosts than their normal SN Ib and Ic siblings; SN IIb has consistent host properties with SN Ib, while hosts of SN IIn are more metal-rich than those of SN II. Hydrogen-deficient superluminous supernovae feature bluer and lower luminosity hosts than most subtypes of core-collapse supernova (CC SN). Assuming simple proportionality of CC SN rates and host star formation rates (SFRs) does not recover the observed mean host properties; either a population of long-lived progenitors or a metallicity-dependent SN production efficiency better reproduces the observed host properties. Assuming the latter case, the rates of SN II are insensitive to host metallicity, but the rates of SN Ib and Ic are substantially enhanced in metal-rich hosts by a factor of ~10 per dex increase in metallicity. Hosts of SN Ia are diverse in their observed properties; subtypes including SN Ia-91T, Ia-02cx, and Ia-CSM prefer star-forming hosts, while subtypes like SN Ia-91bg and Ca-rich prefer quiescent hosts. The rates of SN Ia-91T, Ia-02cx, and Ia-CSM are closely dependent on, or even proportional to, their host SFRs, indicating relatively short-lived progenitors. Conversely, the rates of SN Ia-91bg and Ca-rich transients are proportional to the total stellar mass, favoring long-lived progenitors.

We present a new method to determine the star formation rate (SFR) density of the Universe at $z \gtrsim 5$ that includes the contribution of dust-obscured star formation. For this purpose, we use a [CII] (158 $\mu$m) selected sample of galaxies serendipitously identified in the fields of known $z\gtrsim 4.5$ objects to characterize the fraction of obscured SFR. The advantage of a [CII] selection is that our sample is SFR-selected, in contrast to a UV-selection that would be biased towards unobscured star formation. We obtain a sample of 23 [CII] emitters near star-forming (SF) galaxies and QSOs -- three of which we identify for the first time -- using previous literature and archival ALMA data. 18 of these serendipitously identified galaxies have sufficiently deep rest-UV data and are used to characterize the obscured fraction of the star formation in galaxies with SFRs $\gtrsim 30\ \text{M}_{\odot} \ \text{yr}^{-1}$. We find that [CII] emitters identified around SF galaxies have $\approx$63\% of their SFR obscured, while [CII] emitters around QSOs have $\approx$93\% of their SFR obscured. By forward modeling existing wide-area UV luminosity function (LF) determinations, we derive the intrinsic UV LF using our characterization of the obscured SFR. Integrating the intrinsic LF to $M_{UV}$ = $-$20 we find that the obscured SFRD contributes to $>3\%$ and $>10\%$ of the total SFRD at $z \sim 5$ and $z \sim 6$ based on our sample of companions galaxies near SFGs and QSOs, respectively. Our results suggest that dust obscuration is not negligible at $z\gtrsim 5$, further underlining the importance of far-IR observations of the $z\gtrsim 5$ Universe.

For the first time, we systematically search for galaxies with extended emission line and potential outflows features using medium-band images in the GOODS-S field by comparing the morphology in medium-band images to adjacent continuum and UV bands. We look for galaxies that have a maximum extent 50\% larger, an excess area 30\% greater, or an axis ratio difference of more than 0.3 in the medium band compared to the reference bands. After visual inspection, we find 326 candidate galaxies at $1 < z < 6$, with a peak in the population near cosmic noon, benefiting from the good coverage of the medium-band filters. By examining their SEDs, we find that the candidate galaxies are at least 20\% more bursty in their star-forming activity and have 60\% more young stellar populations compared to a control sample selected based on the continuum band flux. Additionally, these candidates exhibit a significantly higher production rate of ionizing photons. We further find that candidates hosting known AGN produce extended emission that is more anisotropic compared to non-AGN candidates. A few of our candidates have been spectroscopically confirmed to have prominent outflow signatures through NIRSpec observations, showcasing the robustness of the photometric selection. Future spectroscopic follow-up will better help verify and characterize the kinematics and chemical properties of these systems.

We present an analysis of 152,355 radio sources identified in the second data release of the LOFAR Two Metre Sky Survey (LoTSS-DR2) with Sloan Digital Sky Survey (SDSS) spectroscopic redshifts in the range 0.00 < z < 0.57. Using Monte Carlo simulations we determine the reliability of each source exhibiting an excess in radio luminosity relative to that predicted from their Ha emission, and, for a subset of 124,023 sources we combine this measurement with a full BPT analysis. Using these two independent diagnostics we determine the reliability of each source hosting a supermassive black hole of high or low Eddington-scaled accretion rate, and combine the measurements to determine the reliability of sources belonging to each of four physical classes of objects: star forming galaxies (SFGs), radio-quiet active galactic nuclei (RQAGN), and high- or low-excitation radio galaxies (HERGs or emission-line LERGs). The result is a catalogue which enables user-defined samples of radio sources with a reliability threshold suited to their science goal e.g. prioritising purity or completeness. Here we select high-confidence samples of radio sources (>90% reliability) to report: 38,588 radio-excess AGN in the LoTSS DR2 sample (362 HERGs, and 12,648 emission-line LERGs), together with 38,729 SFGs, and 18,726 RQAGN. We validate these results through comparison to literature using independent emission-line measurements, and to widely-adopted WISE photometric selection techniques. While our use of SDSS spectroscopy limits our current analysis to ~4 percent of the LoTSS-DR2 catalogue, our method is directly applicable to data from the forthcoming WEAVE-LOFAR survey which will obtain over a million spectra of 144 MHz selected sources.

The reduced speed of light approximation has been employed to speed up radiative transfer simulations of reionization by a factor of $\gtrsim 5-10$. However, it has been shown to cause significant errors in the HI-ionizing background near reionization's end in simulations of representative cosmological volumes. This can bias inferences on the galaxy ionizing emissivity required to match observables, such as the Ly$\alpha$ forest. In this work, we show that using a reduced speed of light is, to a good approximation, equivalent to re-scaling the global ionizing emissivity in a redshift-dependent way. We derive this re-scaling and show that it can be used to ``correct'' the emissivity in reduced speed of light simulations. This approach of re-scaling the emissivity after the simulation has been run is useful in contexts where the emissivity is a free parameter. We test our method by running full speed of light simulations using these re-scaled emissivities and comparing them with their reduced speed of light counterparts. We find that for reduced speeds of light $\tilde{c} \geq 0.2$, the 21 cm power spectrum at $0.1 \leq k /[h{\rm Mpc}^{-1}] \leq 0.2$ and key Ly$\alpha$ forest observables agree to within $20\%$ throughout reionization, and often better than $10\%$. Position-dependent time-delay effects cause inaccuracies in reionization's morphology on large scales that produce errors up to a factor of $2$ for $\tilde{c} \leq 0.1$. Our method enables a factor of $5$ speedup of radiative transfer simulations of reionization in situations where the emissivity can be treated as a free parameter.

The recent observations by the James Webb Space Telescope (JWST) have revealed a larger number of bright galaxies at $z\gtrsim10$ than was expected. The origin of this excess is still under debate, although several possibilities have been presented. We propose that gamma-ray bursts (GRBs) are a powerful probe to explore the origin of the excess and, hence, the star and galaxy formation histories in the early universe. Focusing on the recently launched mission, Einstein Probe (EP), we find that EP can detect several GRBs annually at $z\gtrsim10$, assuming the GRB formation rate calibrated by events at $z\lesssim6$ can be extrapolated. Interestingly, depending on the excess scenarios, the GRB event rate may also show an excess at $z\simeq10$, and its detection will help to discriminate between the scenarios that are otherwise difficult to distinguish. Additionally, we discuss that the puzzling, red-color, compact galaxies discovered by JWST, the so-called ``little red dots'', could host dark GRBs if they are dust-obscured star forming galaxies. We are eager for unbiased follow-up of GRBs and encourage future missions such as high-z GUNDAM to explore the early universe.

We study the implications of relaxing the requirement for ultralight axions to account for all dark matter in the Universe by examining mixed dark matter (MDM) cosmologies with axion fractions $f \leq 0.3$ within the fuzzy dark matter (FDM) window $10^{-25}$ eV $\lesssim m \lesssim 10^{-23}$ eV. Our simulations, using a new MDM gravity solver implemented in AxiREPO, capture wave dynamics across various scales with high accuracy down to redshifts $z\approx 1$. We identify halos with Rockstar using the CDM component and find good agreement of inferred halo mass functions (HMFs) and concentration-mass relations with theoretical models across redshifts $z=1-10$. This justifies our halo finder approach a posteriori as well as the assumptions underlying the MDM halo model AxionHMcode. Using the inferred axion halo mass - cold halo mass relation $M_{\text{a}}(M_{\text{c}})$ and calibrating a generalised smoothing parameter $\alpha$ to our MDM simulations, we present a new version of AxionHMcode. The code exhibits excellent agreement with simulations on scales $k< 20 \ h$ cMpc$^{-1}$ at redshifts $z=1-3.5$ for $f\leq 0.1$ around the fiducial axion mass $m = 10^{-24.5}$ eV $ = 3.16\times 10^{-25}$ eV, with maximum deviations remaining below 10%. For axion fractions $f\leq 0.3$, the model maintains accuracy with deviations under 20% at redshifts $z\approx 1$ and scales $k< 10 \ h$ cMpc$^{-1}$, though deviations can reach up to 30% for higher redshifts when $f=0.3$. Reducing the run-time for a single evaluation of AxionHMcode to below $1$ minute, these results highlight the potential of AxionHMcode to provide a robust framework for parameter sampling across MDM cosmologies in Bayesian constraint and forecast analyses.

We introduce adaptive particle refinement for compressible smoothed particle hydrodynamics (SPH). SPH calculations have the natural advantage that resolution follows mass, but this is not always optimal. Our implementation allows the user to specify local regions of the simulation that can be more highly resolved. We test our implementation on practical applications including a circumbinary disc, a planet embedded in a disc and a flyby. By comparing with equivalent globally high resolution calculations we show that our method is accurate and fast, with errors in the mass accreted onto sinks of less than 9 percent and speed ups of 1.07-6.62 times for the examples shown. Our method is adaptable and easily extendable, for example with multiple refinement regions or derefinement.

The rise of time-domain astronomy including electromagnetic counterparts to gravitational waves, gravitational microlensing, explosive phenomena, and even astrometry with Gaia, are showing the power and need for surveys with high-cadence, large area, and long time baselines to study the transient universe. A constellation of SmallSats or CubeSats providing wide, instantaneous sky coverage down to 21 Vega mag at optical wavelengths would be ideal for addressing this need. We are assembling CuRIOS-ED (CubeSats for Rapid Infrared and Optical Survey--Exploration Demo), an optical telescope payload which will act as a technology demonstrator for a larger constellation of several hundred 16U CubeSats known as CuRIOS. In preparation for CuRIOS, CuRIOS-ED will launch in late 2025 as part of the 12U Starspec InspireSat MVP payload. CuRIOS-ED will be used to demonstrate the StarSpec ADCS pointing capabilities to <1" and to space-qualify a commercial camera package for use on the full CuRIOS payload. The CuRIOS-ED camera system will utilize a Sony IMX455 CMOS detector delivered in an off-the-shelf Atik apx60 package which we modified to be compatible with operations in vacuum as well as the CubeSat form factor, power, and thermal constraints. By qualifying this commercial camera solution, the cost of each CuRIOS satellite will be greatly decreased (~100x) when compared with current space-qualified cameras with IMX455 detectors. We discuss the CuRIOS-ED mission design with an emphasis on the disassembly, repackaging, and testing of the Atik apx60 for space-based missions. Characterization of the apx60's read noise, dark current, patterned noise, and thermal behavior are reported for a range of temperatures (-35 C to 40 C) and exposure times (0.001s to 30 s). Additionally, we comment on preliminary environmental testing results from a successful thermal vacuum test.

Statistical properties of LSS serve as powerful tools to constrain the cosmological properties of our Universe. Tracing the gas pressure, the tSZ effect is a biased probe of mass distribution and can be used to test the physics of feedback or cosmological models. Therefore, it is crucial to develop robust modeling of hot gas pressure for applications to tSZ surveys. Since gas collapses into bound structures, it is expected that most of the tSZ signal is within halos produced by cosmic accretion shocks. Hence, simple empirical halo models can be used to predict the tSZ power spectra. In this study, we employed the HMx halo model to compare the tSZ power spectra with those of several hydrodynamical simulations: the Horizon suite and the Magneticum simulation. We examined various contributions to the tSZ power spectrum across different redshifts, including the one- and two-halo term decomposition, the amount of bound gas, the importance of different masses and the electron pressure profiles. Our comparison of the tSZ power spectrum reveals discrepancies that increase with redshift. We find a 20% to 50% difference between the measured and predicted tSZ angular power spectrum over the multipole range $\ell=10^3-10^4$. Our analysis reveals that these differences are driven by the excess of power in the predicted two-halo term at low k and in the one-halo term at high k. At higher redshifts (z~3), simulations indicate that more power comes from outside the virial radius than from inside suggesting a limitation in the applicability of the halo model. We observe differences in the pressure profiles, despite the fair level of agreement on the tSZ power spectrum at low redshift with the default calibration of the halo model. In conclusion, our study suggests that the properties of the halo model need to be carefully controlled against real or mock data to be proven useful for cosmological purposes.

The Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst (CHIME/FRB) Project has a new VLBI Outrigger at the Green Bank Observatory (GBO), which forms a 3300km baseline with CHIME operating at 400-800MHz. Using 100ms long full-array baseband "snapshots" collected commensally during FRB and pulsar triggers, we perform a shallow, wide-area VLBI survey covering a significant fraction of the Northern sky targeted at the positions of compact sources from the Radio Fundamental Catalog. In addition, our survey contains calibrators detected from two 1s long trial baseband snapshots for a deeper survey with CHIME and GBO. In this paper, we present the largest catalog of compact calibrators suitable for 30-milliarcsecond-scale VLBI observations at sub-GHz frequencies to date. Our catalog consists of 200 total calibrators in the Northern Hemisphere that are compact on 30-milliarcsecond scales with fluxes above 100mJy. This calibrator grid will enable the precise localization of hundreds of FRBs a year with CHIME/FRB-Outriggers.

The Mittelman-di Cicco-Walker (MDW) H$\alpha$ Sky Survey is an autonomously-operated and ongoing all-sky imaging survey in the narrowband H$\alpha$ wavelength. The survey was founded by amateur astronomers, and is presented here in its first stage of refinement for rigorous scientific use. Each field is exposed through an H$\alpha$ filter with a 3nm bandwidth for a total of four hours, with a pixel scale of 3.2 arcsec. Here, we introduce the first Data Release of the MDW H$\alpha$ Survey (Data Release 0, or DR0), spanning 238 fields in the region of Orion (~3100 deg$^2$). DR0 includes: calibrated mean fields, star-removed mean fields, a point source catalog matched to Data Release 1 of the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS1) and the INT Galactic Plane Survey (IGAPS), and mosaics.

Aims. The goal of this work is to comprehensively characterize the impact of stellar multiplicity on Class II disks in the L1641 and L1647 regions of Orion A (~1-3 Myr), part of the Survey of Orion Disks with ALMA (SODA). We characterize the protostellar multiplicity using the Atacama Large Millimeter/submillimeter Array (ALMA), the ESO-VISTA, and Hubble Space telescopes. The resulting sample of 65 multiple systems represents the largest catalogue of wide binary systems to date (projected separation >1000 AU), allowing a more robust statistical characterization of the evolution and properties of protoplanetary disks. Methods. The disk population was observed in continuum with ALMA at 225 GHz, with a median rms of 1.5 Mearth. Combining these data (resolution ~1.1arcsec ) with the ESO-VISTA near-infrared survey of the Orion A cloud (resolution ~0.7arcsec ), multiple systems are assembled and selected by an iterative inside-out search in projected separation (>1000 AU). Results. We identify 61 binary systems, 3 triple systems, and one quadruple system. The separation range is between 1000 and 10^4 AU. The dust mass distributions inferred with the Kaplan-Meier estimator yield a median mass of 3.23+0.6-0.4 Mearth for primary disks and 3.88+0.3-0.3 Mearth for secondary disks.

Many accretion discs have been found to be distorted: either warped due a misalignment in the system, or non-circular as a result of orbital eccentricity or tidal deformation by a binary companion. Warped, eccentric, and tidally distorted discs are not in vertical hydrostatic equilibrium, and thus exhibit vertical oscillations in the direction perpendicular to the disc, a phenomenon that is absent in circular and flat discs. In extreme cases, this vertical motion is manifested as a vertical `bouncing' of the gas, potentially leading to shocks and heating, as observed in recent global numerical simulations. In this paper we isolate the mechanics of vertical disc oscillations by means of quasi-2D and fully 3D hydrodynamic local (shearing-box) models. To determine the numerical and physical dissipation mechanisms at work during an oscillation we start by investigating unforced oscillations, examining the effect of initial oscillation amplitude, as well as resolution, boundary conditions, and vertical box size on the dissipation and energetics of the oscillations. We then drive the oscillations by introducing a time-dependent gravitational potential. A key result is that even a purely vertically oscillating disc is (parametrically) unstable to developing inertial waves, as we confirm through a linear stability analysis. The most important of these has the character of a bending wave, whose radial wavelength depends on the frequency of the vertical oscillation. The nonlinear phase of the instability exhibits shocks, which dampen the oscillations, although energy can also flow from the bending wave back to the vertical oscillation.

The largest geomagnetic storm in two decades occurred in 2024 May with a minimum $D_{\rm st}$ of $-412$ nT. We examine its solar and interplanetary origins by combining multipoint imaging and in situ observations. The source active region, NOAA AR 13664, exhibited extraordinary activity and produced successive halo eruptions, which were responsible for two complex ejecta observed at the Earth. In situ measurements from STEREO A, which was $12.6^{\circ}$ apart, allow us to compare the ``geo-effectiveness" at the Earth and STEREO A. We obtain key findings concerning the formation of solar superstorms and how mesoscale variations of coronal mass ejections affect geo-effectiveness: (1) the 2024 May storm supports the hypothesis that solar superstorms are ``perfect storms" in nature, i.e., a combination of circumstances resulting in an event of an unusual magnitude; (2) the first complex ejecta, which caused the geomagnetic superstorm, shows considerable differences in the magnetic field and associated ``geo-effectiveness" between the Earth and STEREO A, despite a mesoscale separation; and (3) two contrasting cases of complex ejecta are found in terms of the geo-effectiveness at the Earth, which is largely due to different magnetic field configurations within the same active region.

Magnetars are neutron stars that host huge, complex magnetic fields which require supporting currents to flow along the closed field lines. This makes magnetar atmospheres different from those of passively cooling neutron stars because of the heat deposited by backflowing charges impinging on the star surface layers. This particle bombardment is expected to imprint the spectral and, even more, the polarisation properties of the emitted thermal radiation. We present solutions for the radiative transfer problem for bombarded plane-parallel atmospheres in the high magnetic field regime. The temperature profile is assumed a priori, and selected in such a way to reflect the varying rate of energy deposition in the slab (from the impinging currents and/or from the cooling crust). We find that thermal X-ray emission powered entirely by the energy released in the atmosphere by the magnetospheric back-bombardment is linearly polarised and X-mode dominated, but its polarisation degree is significantly reduced (down to $10-50\%$) when compared with that expected from a standard atmosphere heated only from the cooling crust below. By increasing the fraction of heat flowing in from the crust the polarisation degree of the emergent radiation increases, first at higher energies ($\sim 10$ keV) and then in the entire soft X-ray band. We use our models inside a ray-tracing code to derive the expected emission properties as measured by a distant observer and compare our results with recent IXPE observations of magnetar sources.

This Chapter outlines the basic properties of waves in solar partially ionized plasmas. It provides a summary of the main sets of equations, from the single-fluid formalism, to the multi-fluid one, giving examples for purely hydrogen, and for hydrogen-helium plasmas. It then discusses the solutions for waves under the single-fluid frame: the influence of the ambipolar diffusion, diamagnetic effect, and the Hall effect on the propagation, dissipation, and mode conversion of the magnetohydrodynamic waves. The Chapter continues by outlining the wave solutions in the multi-fluid formalism: the influence of the elastic inter-particle collisions into the propagation, damping and dissipation of different magnetohydrodynamic modes. Both parts discuss linear and non-linear wave solutions, and the effects of the gravitational stratification of the solar atmosphere.

The identification of persistent radio sources (PRSs) coincident with two repeating fast radio bursts (FRBs) supports FRB theories requiring a compact central engine. However, deep non-detections in other cases highlight the diversity of repeating FRBs and their local environments. Here, we perform a systematic search for radio sources towards 37 CHIME/FRB repeaters using their arcminute localizations and a combination of archival surveys and targeted observations. Through multi-wavelength analysis of individual radio sources, we identify two (20181030A-S1 and 20190417A-S1) for which we disfavor an origin of either star formation or an active galactic nucleus in their host galaxies and thus consider them candidate PRSs. We do not find any associated PRSs for the majority of the repeating FRBs in our sample. For 8 FRB fields with Very Large Array imaging, we provide deep limits on the presence of PRSs that are 2--4 orders of magnitude fainter than the PRS associated with FRB\,20121102A. Using Very Large Array Sky Survey imaging of all 37 fields, we constrain the rate of luminous ($\gtrsim$10$^{40}$ erg s$^{-1}$) PRSs associated with repeating FRBs to be low. Within the context of FRB-PRS models, we find that 20181030A-S1 and 20190417A-S1 can be reasonably explained within the context of magnetar, hypernebulae, gamma-ray burst afterglow, or supernova ejecta models -- although we note that both sources follow the radio luminosity versus rotation measure relationship predicted in the nebula model framework. Future observations will be required to both further characterize and confirm the association of these PRS candidates with the FRBs.

Context. Magnetic bright points (MBPs) are one of the smallest manifestations of the magnetic field in the solar atmosphere and are observed to extend from the photosphere up to the chromosphere. As such, they represent an excellent feature to use in searches for types of magnetohydrodynamic (MHD) waves and mode coupling in the solar atmosphere. Aims. In this work, we aim to study wave propagation in the lower solar atmosphere by comparing intensity oscillations in the photosphere with the chromosphere via a search for possible mode coupling, in order to establish the importance of these types of waves in the solar atmosphere, and their contribution to heating the chromosphere. Methods. Observations were conducted in July 2011 with the ROSA and the HARDCam instruments at the Dunn Solar Telescope. We used wavelet analysis to identify traveling MHD waves and derive frequencies in the G-band and H$\alpha$wave bands. We isolated a large sample of MBPs using an automated tracking algorithm throughout our observations. Two dozen of the brightest MBPs were selected from the sample for further study. Results. We find oscillations in the G-band MBPs, with frequencies between 1.5 and 3.6 mHz. Corresponding MBPs in the lower solar chromosphere observed in H$\alpha$ show a frequency range of 1.4 to 4.3 mHz. In about 38\% of the MBPs, the ratio of H$\alpha$ to G-band frequencies was near two. Thus, these oscillations show a form of mode coupling where the transverse waves in the photosphere are converted into longitudinal waves in the chromosphere. Conclusions. From simple estimates we find an energy flux of $\approx$45 $\times 10^{3}$ W m$^{-2}$ and show that the energy flowing through MBPs is enough to heat the chromosphere, and mode coupling is important in helping us understand the types of MHD waves in the lower solar atmosphere and the overall energy budget.

We present a new algorithm for identifying superbubbles in HI column density maps of both observed and simulated galaxies that has only a single adjustable parameter. The algorithm includes an automated galaxy-background separation step to focus the analysis on the galactic disk. To test the algorithm, we compare the superbubbles it finds in a simulated galactic disk with the ones it finds in 21cm observations of a similar galactic disk. The sizes and radial distribution of those superbubbles are indeed qualitatively similar. However, superbubbles in the simulated galactic disk have lower central HI column densities. The HI superbubbles in the simulated disk are spatially associated with pockets of hot gas. We conclude that the algorithm is a promising method for systematically identifying and characterizing superbubbles using only HI column density maps that will enable standardized tests of stellar feedback models used in galaxy simulations.

Pulsars are rotating neutron stars that are observed to be slowing down, implying a loss of their kinetic energy. There can be several different physical mechanisms involved in their spin-down process. The properties of fast-rotating pulsars depend on the nature of the neutron star matter, which can also affect the spin-down mechanisms. In this work, we examine three different physical phenomena contributing to the spin-down: magnetic dipole radiation, gravitational mass quadrupole radiation due to the ``mountain" formation, gravitational mass current quadrupole radiation or the r-modes, and calculate the expressions for the braking indices due to all of them. We have also considered jointly the implications of the uncertainties of the equation of the state of neutron star matter and rapid rotation on the braking indices corresponding to the aforementioned processes and their combinations. In all cases, the rapid rotation results in a departure from the standard values in the literature for the braking index when the rotational effects are ignored. If generated with a saturation amplitude within the range of $10^{-4} - 10^{-1}$, the r-mode oscillations dominate the spin-down of millisecond pulsars. We also explore the braking index in the context of millisecond magnetars. This study examines two braking index measurements in the context of newly born millisecond magnetars from two observed short $\gamma$-ray bursts. The measured braking indices for these objects are consistent with our estimation, which allows us to conclude that the spin frequency of the remnants is within the range of $\sim 550-850$ Hz.

orbitize! is a package for Bayesian modeling of the orbital parameters of resolved binary objects from time series measurements. It was developed with the needs of the high-contrast imaging community in mind, and has since also become widely used in the binary star community. A generic orbitize! use case involves translating relative astrometric time series, optionally combined with radial velocity or astrometric time series, into a set of derived orbital posteriors. This paper is published alongside the release of orbitize! version 3.0, which has seen significant enhancements in functionality and accessibility since the release of version 1.0 (Blunt et al., 2020).

We explore the Schmidt-Kennicutt (SK) relations and the star formation efficiency per free-fall time ($\eff$), mirroring observational studies, in numerical simulations of filamentary molecular clouds undergoing gravitational contraction. We find that {\it a)} collapsing clouds accurately replicate the observed SK relations for galactic clouds and {\it b)} the so-called efficiency per free-fall time ($\eff$) is small and constant in space and in time, with values similar to those found in local clouds. This constancy is a consequence of the similar radial scaling of the free-fall time and the internal mass in density structures with spherically-averaged density profiles near $r^{-2}$. We additionally show that {\it c)} the star formation rate (SFR) increases rapidly in time; {\it d)} the low values of $\eff$ are due to the different time periods over which $\tauff$ and $\tausf$ are evaluated, together with the fast increasing SFR, and {\it e)} the fact that star clusters are significantly denser than the gas clumps from which they form is a natural consequence of the fast increasing SFR, the continuous replenishment of the star-forming gas by the accretion flow, and the near $r^{-2}$ density profile generated by the collapse Finally, we argue that the interpretation of $\eff$ as an efficiency is problematic because its maximum value is not bounded by unity, and because the total gas mass in the clouds is not fixed, but rather depends on the environment where clouds are embedded. In summary, our results show that the SK relation, the typical observed values of $\eff$, and the mass density of clusters arise as a natural consequence of gravitational contraction.

We have investigated the evolutionary connections of the isolated neutron star (NS) populations including radio pulsars (RPs), anomalous X-ray pulsars (AXPs), soft gamma repeaters (SGRs), dim isolated NSs (XDINs), ``high-magnetic-field'' RPs (``HBRPs''), central compact objects (CCOs), rotating radio transients (RRATs), and long-period pulsars (LPPs) in the fallback disc model. The model can reproduce these NS families as a natural outcome of different initial conditions (initial period, disc mass, and dipole moment, $\mu$) with a continuous $\mu$ distribution in the $\sim 10^{27} - 5 \times 10^{30}$ G cm$^3$ range. Results of our simulations can be summarised as follows: (1) A fraction of ``HBRPs'' with relatively high $\mu$ evolve into the persistent AXP/SGR properties, and subsequently become LPPs. (2) Persistent AXP/SGRs do not have evolutionary links with CCOs, XDINs, and RRATs. (3) For a wide range of $\mu$, most RRATs evolve passing through RP or ``HBRP'' properties during their early evolutionary phases. (4) A fraction of RRATs which have the highest estimated birth rate seem to be the progenitors of XDINs. (5) LPPs, whose existence was predicted by the fallback disc model, are the sources evolving in the late stage of evolution before the discs become inactive. These results provide concrete support to the ideas proposing evolutionary connections between the NS families to account for the ``birth-rate problem'', the discrepancy between the cumulative birth rate estimated for these systems and the core-collapse supernova rate.

The tumultuous effects of ultraviolet photons that source cosmic reionization, the subsequent compression and shock-heating of low-density regions, and the modulation of baryons in shallow potential wells induced by the passage of ionization fronts, collectively introduce perturbations to the evolution of the intergalactic medium in the post-reionization era. These enduring fluctuations persist deep into the post-reionization era, casting a challenge upon precision cosmology endeavors targeting tracers in this cosmic era. Simultaneously, these relics from reionization also present a unique opportunity to glean insights into the astrophysics that govern the epoch of reionization. In this work, we propose a first study of the cross-correlation of \lya forest and 21 cm intensity mapping, accounting for the repercussions of inhomogeneous reionization in the post-reionization era. We investigate the ability of SKA $\times$ DESI-like, SKA $\times$ MUST-like, and PUMA $\times$ MUST-like instrumental setups to achieve a high signal-to-noise ratio (SNR) in the redshift range $3.5 \leq z \leq 4$. Moreover, we assess how alterations in integration time, survey area, and reionization scenarios impact the SNR. Furthermore, we forecast the cross-correlation's potential to constrain cosmological parameters under varying assumptions: considering or disregarding reionization relics, marginalizing over reionization astrophysics, and assuming perfect knowledge of reionization. Notably, our findings underscore the remarkable capability of a futuristic PUMA $\times$ MUST-like setup, with a modest 100-hour integration time over a 100 sq. deg. survey, to constrain the ionization efficiency error to $\sigma_\zeta = 3.42 $.

The TRAPPIST-1 system has been extensively observed with JWST in the near-infrared with the goal of measuring atmospheric transit transmission spectra of these temperate, Earth-sized exoplanets. A byproduct of these observations has been much more precise times of transit compared with prior available data from Spitzer, HST, or ground-based telescopes. In this note we use 23 new timing measurements of all seven planets in the near-infrared from five JWST observing programs to better forecast and constrain the future times of transit in this system. In particular, we note that the transit times of TRAPPIST-1h have drifted significantly from a prior published analysis by up to tens of minutes. Our newer forecast has a higher precision, with median statistical uncertainties ranging from 7-105 seconds during JWST Cycles 4 and 5. Our expectation is that this forecast will help to improve planning of future observations of the TRAPPIST-1 planets, whereas we postpone a full dynamical analysis to future work.

We trained denoiser autoencoding neural networks on medium resolution simulated optical spectra of late-type stars to demonstrate that the reconstruction of the original flux is possible at a typical relative error of a fraction of a percent down to a typical signal-to-noise ratio of 10 per pixel. We show that relatively simple networks are capable of learning the characteristics of stellar spectra while still flexible enough to adapt to different values of extinction and fluxing imperfections that modifies the overall shape of the continuum, as well as to different values of Doppler shift. Denoised spectra can be used to find initial values for traditional stellar template fitting algorithms and - since evaluation of pre-trained neural networks is significantly faster than traditional template fitting - denoiser networks can be useful when a fast analysis of the noisy spectrum is necessary, for example during observations, between individual exposures.

The understanding of the formation and evolution of the solar system still has many unanswered questions. Formation of solids in the solar system, mineral and organic mixing, and planetary body creation are all topics of interest to the community. Studying these phenomena is often performed through observations, remote sensing, and in-situ analysis, but there are limitations to the methods. Limitations such as IR diffraction limits, spatial resolution issues, and spectral resolution issues can prevent detection of organics, detection and identification of cellular structures, and the disentangling of granular mixtures. Optical-PhotoThermal InfraRed (O-PTIR) spectroscopy is a relatively new method of spectroscopy currently used in fields other than planetary sciences. O-PTIR is a non-destructive, highly repeatable, and fast form of measurement capable of reducing these limitations. Using a dual laser system with an IR source tuned to the mid-IR wavelength we performed laboratory O-PTIR measurements to compare O-PTIR data to existing IR absorption data and laboratory FTIR measurements for planetary materials. We do this for the purpose of introducing O-PTIR to the planetary science community. The technique featured here would serve to better measurements of planetary bodies during in-situ analysis. We find that, unlike other fields where O-PTIR produces almost one-to-one measurements with IR absorption measurements of the same material, granular materials relevant to planetary science do not. However, we do find that the materials compared were significantly close and O-PTIR was still capable of identifying materials relevant to planetary science.

Twenty-one-centimetre signals from the Epoch of Reionization (EoR) are expected to be detected in the low-frequency radio window by the next-generation interferometers, particularly the Square Kilometre Array (SKA). However, precision data analysis pipelines are required to minimize the systematics within an infinitesimal error budget. Consequently, there is a growing need to characterize the sources of errors in EoR analysis. In this study, we identify one such error origin, namely source blending, which is introduced by the overlap of objects in the densely populated observing sky under SKA1-Low's unprecedented sensitivity and resolution, and evaluate its two-fold impact in both the spatial and frequency domains using a novel hybrid evaluation (HEVAL) pipeline combining end-to-end simulation with an analytic method to mimic EoR analysis pipelines. Sky models corrupted by source blending induce small but severe frequency-dependent calibration errors when coupled with astronomical foregrounds, impeding EoR parameter inference with strong additive residuals in the two-dimensional power spectrum space. We report that additive residuals from poor calibration against sky models with blending ratios of 5 and 0.5 per cent significantly contaminate the EoR window. In contrast, the sky model with a 0.05 per cent blending ratio leaves little residual imprint within the EoR window, therefore identifying a blending tolerance at approximately 0.05 per cent. Given that the SKA observing sky is estimated to suffer from an extended level of blending, strategies involving de-blending, frequency-dependent error mitigation, or a combination of both, are required to effectively attenuate the calibration impact of source-blending defects.

Cloud-cloud collisions in splash bridges produced in gas-rich disk galaxy collisions offer a brief but interesting environment to study the effects of shocks and turbulence on star formation rates in the diffuse IGM, far from the significant feedback effects of massive star formation and AGN. Expanding on our earlier work, we describe simulated collisions between counter-rotating disk galaxies of relatively similar mass, focusing on the thermal and kinematic effects of relative inclination and disk offset at the closest approach. This includes essential heating and cooling signatures, which go some way towards explaining the luminous power in H$_2$ and [CII] emission in the Taffy bridge, as well as providing a partial explanation of the turbulent nature of the recently observed compact CO-emitting clouds observed in Taffy by ALMA. The models show counter-rotating disk collisions result in swirling, shearing kinematics for the gas in much of the post-collision bridge. Gas with little specific angular momentum due to collisions between counter-rotating streams accumulates near the center of mass. The disturbances and mixing in the bridge drive continuing cloud collisions, differential shock heating, and cooling throughout. A wide range of relative gas phases and line-of-sight velocity distributions are found in the bridges, depending sensitively on initial disk orientations and the resulting variety of cloud collision histories. Most cloud collisions can occur promptly or persist for quite a long duration. Cold and hot phases can largely overlap throughout the bridge or can be separated into different parts of the bridge.

Accretion disc outbursts are re-occurring events observed in various astrophysical systems, including X-ray binaries and cataclysmic variables. These outbursts are characterized by a sudden increase in luminosity due to various instabilities in the accretion disc. We need to investigate the time-dependent accretion flow models to understand the mechanisms driving these outbursts. Time-dependent models incorporate the disc's time evolution and can capture the build-up of instabilities. This review aims to give a basic overview of accretion disc outburst and stability analysis. The paper highlights the necessity of considering the hierarchy of different timescales, dynamical, viscous, and thermal, when investigating the instabilities occurring in the accretion disc. The importance and observational implications of studying these accretion disc outbursts are also discussed.

Reliable photo-nuclear reaction rates at the stellar conditions are essential to understand the origin of the heavy stable neutron-deficient isotopes between $^{74}$Se and $^{196}$Hg-p-nuclei, however, many reaction rates of relevance still have to rely on the Hauser-Feshbach model due to rare experimental progress. One such case is in the mass range of 160 for Dy, Er, Ho and Tm isotopes. In this work we attempt to constrain the Hauser-Feshbach model in the TALYS package by reproducing the available experimental data of $^{160}$Dy($p,\gamma$)$^{161}$Ho and $^{162}$Er($p,\gamma$)$^{163}$Tm in the $A\sim 160$ mass region, and examine the effects of level density, gamma strength function and the optical model potential. The constrained model then allows us to calculate the reaction rates of $^{157, 159}$Ho($\gamma$, $p$) and $^{163,165}$Tm($\gamma$, $p$) for the $\gamma$-process nucleosynthesis in carbon-deflagration SNe Ia model. Our recommended rates differ from the JINA REACLIB by more than 1 order of magnitude in the temperature range of 2-3 GK. This results in the changes of final abundance of $p$-nuclei in the $A\sim 160$ mass range by -5.5-3\% from those with JINA, which means that the ($\gamma$, $p$) reactions uncertainty is not predominant for the synthesis of these nuclei.

CXOU J110926.4-650224 is a candidate transitional millisecond pulsar (tMSP) with X-ray and radio emission properties reminiscent of those observed in confirmed tMSPs in their X-ray 'subluminous' disc state. We present the results of observing campaigns that, for the first time, characterise the optical and near-infrared variability of this source and establish a connection with the mode-switching phenomenon observed in X-rays. The optical emission exhibited flickering activity, frequent dipping episodes where it appeared redder, and a multi-peaked flare where it was bluer. The variability pattern was strongly correlated with that of the X-ray emission. Each dip matched an X-ray low-mode episode, indicating that a significant portion of the optical emission originates from nearly the same region as the X-ray emission. The near-infrared emission also displayed remarkable variability, including a dip of 20 min in length during which it nearly vanished. Time-resolved optical spectroscopic observations reveal significant changes in the properties of emission lines from the disc and help infer the spectral type of the companion star to be between K0 and K5. We compare the properties of CXOU J110926.4-650224 with those of other tMSPs in the X-ray subluminous disc state and discuss our findings within the context of a recently proposed scenario that explains the phenomenology exhibited by the prototypical tMSP PSR J1023+0038.

In this work, we performed a spectro-temporal investigation of the low-mass X-ray binary GX 9+9 using the Large Area X-ray Proportional Counter (LAXPC) and Soft X- ray Telescope (SXT) observation on board AstroSat. The source was detected in the soft state during the observation, which results in a disk dominating energy spectrum within the energy range of 0.7-25.0 keV. We carried out the analysis at different flux levels. In the temporal analysis, LAXPC data in all flux levels showed the presence of noise components, describing broad Lorentzian components. We modeled the energy-dependent temporal properties of the source in order to identify the radiative origin of the observed variability. This source is not a well-studied source; hence we attempt to estimate various source characteristics like inner-disk radius, flux, and inner-disk temperature.

One way to understand planet formation is through studying the correlations between planet occurrence rates and stellar mass. However, measuring stellar mass in the red giant regime is very difficult. In particular, the spectroscopic masses of certain evolved stars, often referred to as "retired A-stars", have been questioned in the literature. Efforts to resolve this mass controversy using spectroscopy, interferometry and asteroseismology have so far been inconclusive. A recent ensemble study found a mass-dependent mass offset, but the result was based on only 16 stars. With NASA's Transiting Exoplanet Survey Satellite (TESS), we expand the investigation of the mass discrepancy to a total of 92 low-luminosity stars, synonymous with the retired A-stars. We measure their characteristic oscillation frequency, $\mathrm{\nu}_{\mathrm{max}}$, and the large frequency separation, $\mathrm{\Delta\nu}$, from their TESS photometric time series. Using these measurements and asteroseismic scaling relations, we derive asteroseismic masses and compare them with spectroscopic masses from five surveys, to comprehensively study the alleged mass-dependent mass offset. We find a mass offset between spectroscopy and seismology that increases with stellar mass. However, we note that adopting the seismic mass scale does not have a significant effect on the planet occurrence-mass-metallicity correlation for the so-called retired A-stars. We also report seismic measurements and masses for 157 higher luminosity giants (mostly helium-core-burning) from the spectroscopic surveys.

We report on the detection of unintended electromagnetic radiation (UEMR) from the second-generation of Starlink satellites. Observations with the LOFAR radio telescope between 10 to 88MHz and 110 to 188MHz show broadband emission covering the frequency ranges from 40 to 70MHz and 110 to 188MHz from the v2-Mini and v2-Mini Direct-to-Cell Starlink satellites. The spectral power flux density of this broadband UEMR varies from satellite to satellite, with values ranging from 15Jy to 1300Jy, between 56 and 66MHz, and from 2 to 100Jy over two distinct 8MHz frequency ranges centered at 120 and 161MHz. We compared the detected power flux densities of this UEMR to that emitted by the first generation v1.0 and v1.5 Starlink satellites. When correcting for the observed satellite distances, we find that the second-generation satellites emit UEMR that is up to a factor of 32 stronger compared to the first generation. The calculated electric field strengths of the detected UEMR exceed typical electromagnetic compatibility standards used for commercial electronic devices as well as recommended emission thresholds from the Radiocommunication Sector of the International Telecommunications Union (ITU-R) aimed at protecting the 150.05-153MHz frequency range allocated to radio astronomy. We characterize the properties of the detected UEMR with the aim of assisting the satellite operator with the identification of the cause of the UEMR.

This brief review is based on a lecture given by one of the authors at the international youth conference AYSS-2023. It is devoted to multimessenger astronomy, which studies astrophysical objects and phenomena using various particles and waves that bring information from space. The messengers include electromagnetic and gravitational waves, neutrinos, and cosmic rays. We discuss new opportunities that open up with the combined use of several carriers of information. Combination of data obtained through various observation channels allows one to obtain more complete and accurate information about the processes occurring in the Universe, and even to use it for studying fundamental physics.

We present a new constraint on the Galactic $^{12}$C/$^{13}$C gradient with sensitive HCO$^+$ absorption observations against strong continuum sources. The new measurements suffer less from beam dilution, optical depths, and chemical fractionation, allowing us to derive the isotopic ratios precisely. The measured $^{12}$C/$^{13}$C ratio in the Solar neighborhood (66$\pm$5) is consistent with those obtained from CH$^+$. Two measurements toward the Galactic Center are 42.2$\pm$1.7 and 37.5$\pm$6.5. Though the values are a factor of 2$\sim$3 higher than those derived from dense gas tracers (e.g., H$_2$CO, complex organic molecules) toward Sagittarius (Sgr) B2 regions, our results are consistent with the absorption measurements from c-C$_3$H$_2$ toward Sgr B2 ($\sim$40), and those from CH$^+$ toward Sgr A$^*$ and Sgr B2(N) ($>$30). We calculate a new Galactic $^{12}$C/$^{13}$C gradient of (6.4$\pm$1.9)$R_{\rm GC}$/kpc+(25.9$\pm$10.5), and find an increasing trend of $^{12}$C/$^{13}$C gradient obtained from high-density to low-density gas tracers, suggesting opacity effects and chemical fractionation may have a strong impact on the isotopic ratios observed at high-density regions.

Context: Mutual neutralization between cations and anions play an important role in determining the charge-balance in certain astrophysical environments. However, empirical data for such reactions involving complex molecular species has been lacking due to challenges in performing experimental studies, leaving the astronomical community to rely on decades old models with large uncertainties for describing these processes in the interstellar medium. Aims: To investigate the mutual neutralization (MN) reaction, C$_{60}^+$ + C$_{60}^-$ $\rightarrow$ C$_{60}^*$ + C$_{60}$, for collisions at interstellar-like conditions. Methods: The mutual neutralization reaction between C$_{60}^+$ and C$_{60}^-$ at collision energies of 100\,meV was studied using the Double ElectroStatic Ion Ring ExpEriment, DESIREE, and its merged-beam capabilities. To aid in the interpretation of the experimental results, semi-classical modeling based on the Landau-Zener approach was performed for the studied reaction. Results: We experimentally identify a narrow range of kinetic energies for the neutral reaction products. Modeling was used to calculate the quantum state-selective reaction probabilities, absolute cross sections, and rate coefficients of these MN reactions, using the experimental results as a benchmark. The MN cross sections are compared with model results for electron attachment to C$_{60}$ and electron recombination with C$_{60}^+$. Conclusions: The present results show that it is crucial to take mutual polarization effects, the finite sizes, and the final quantum states of both molecular ions into account for reliable predictions of MN rates expected to strongly influence the charge-balance and chemistry in, e.g., dense molecular clouds.

We present a spectro-timing analysis of the black hole X-ray transient GX 339-4 using simultaneous observations from Astrosat and NICER during the 2021 outburst period. The combined spectrum obtained from NICER, LAXPC, and SXT data is effectively described by a model comprising a thermal disk component, hard Comptonization component, and reflection component with an edge. Our analysis of the Astrosat and NICER spectra indicates the source to be in a low/hard state, with a photon index of ~1.64. The Power Density Spectra (PDS) obtained from both Astrosat and NICER observations exhibit two prominent broad features at 0.22 Hz and 2.94 Hz. We generated energy-dependent time lag and fractional root mean square (frms) at both frequencies in a broad energy range of 0.5-30 keV and found the presence of hard lags along with a decrease in variability at higher energy levels. Additionally, we discovered that the correlated variations in accretion rate, inner disc radius, coronal heating rate, and the scattering fraction, along with a delay between them, can explain the observed frms and lag spectra for both features.

Ultra-high-energy neutrinos and cosmic rays are excellent probes of astroparticle physics phenomena. For astroparticle physics analyses, robust and accurate reconstruction of signal parameters like arrival direction and energy is essential. Current reconstruction methods ignore bin-to-bin noise correlations, which limits reconstruction resolution and so far has prevented calculations of event-by-event uncertainties. In this work, we present a likelihood description of neutrino or cosmic-ray signals in a radio detector with correlated noise, as present in all neutrino and cosmic-ray radio detectors. We demonstrate with a toy-model reconstruction that signal parameters such as energy and direction, including event-by-event uncertainties with correct coverage, can be obtained. Additionally, by correctly accounting for correlations, the likelihood description constrains the best-fit parameters better than alternative methods and thus improves experimental reconstruction capabilities.

GRS 1716-249 is a stellar-mass black hole in a low-mass X-ray binary that underwent a gaint outburst in 2016/17. In this paper we use simultaneous observations of Insight-HXMT and NuSTAR to determine its basic parameters. The observations were performed during the softest part of the outburst, and the spectra show clear thermal disk emission and reflection features. We have fitted the X-ray energy spectra using the joint fitting method of the continuum and reflection components with the kerrbb2+ relxill model. Since there is a possibility that the distance to this source was previously underestimated, we use the latest distance parameter of 6.9 kpc in our study, in contrast to previous work in which the distance was set at 2.4 kpc. Through spectral fitting of fixing black hole mass at 6.4 $M_{\rm \odot}$, we observe a strong dependence of the derived spin on the distance: $a_{*}=0.972_{-0.005}^{+0.004}$ at an assumed distance of 2.4 kpc and $a_{*}=0.464_{-0.007}^{+0.016}$ at an assumed distance of 6.9 kpc, at a confidence level of 90%. If considering the uncertainties in the distance and black hole mass, there will be a wider range of spin with $a_{*}$ < 0.78. The fitting results with the new distance indicate that GRS 1716-249 harbors a moderate spin black hole with an inclined ($i\sim 40-50^{\circ}$) accretion disk around it. Additionally, we have also found that solely using the method of the reflection component fitting but ignoring the constraints on the spin from the accretion disk component will result in an extremely high spin.

SWEET-Cat (Stars With ExoplanETs Catalogue) was originally introduced in 2013, and since then, the number of confirmed exoplanets has increased significantly. A crucial step for a comprehensive understanding of these new worlds is the precise and homogeneous characterization of their host stars. We used a large number of high-resolution spectra to continue the addition of new stellar parameters for planet-host stars in SWEET-Cat following the new detection of exoplanets listed both at the Extrasolar Planets Encyclopedia and at the NASA exoplanet archive. We obtained high-resolution spectra for a significant number of these planet-host stars, either observed by our team or collected through public archives. For FGK stars, the spectroscopic stellar parameters were derived for the spectra following the same homogeneous process using ARES+MOOG as for the previous SWEET-Cat releases. The stellar properties are combined with the planet properties to study possible correlations that could shed more light into the star-planet connection studies. We increase the number of stars with homogeneous parameters by 232 ($\sim$ 25\% - from 959 to 1191). We then focus on the exoplanets with both mass and radius determined to review the mass-radius relation where we find consistent results with the ones previously reported in the literature. For the massive planets we also revisit the radius anomaly where we confirm a metallicity correlation for the radius anomaly already hinted in previous results.

We propose a novel approach for determining the orbital inclination of low-mass X-ray binary systems by modelling the H$\alpha$ and H$\beta$ line profiles emitted by the accretion disc, with a Newtonian version of diskline. We applied the model to two sample sources, Swift J1357.2-0933 and MAXI J1305-704, which are both transient black hole systems, and analyse two observations that were collected during a quiescent state and one observation of an outburst. The line profile is well described by the diskline model, although we had to add a Gaussian line to describe the deep inner core of the double-peaked profile, which the diskline model was unable to reproduce. The H$\beta$ emission lines in the spectrum of Swift J1357.2-0933 and the H$\alpha$ emission lines in that of MAXI J1305-704 during the quiescent state are consistent with a scenario in which these lines originate from a disc ring between $(9.6-57) \times 10^{3}, \rm{R_{g}}$ and $(1.94-20) \times 10^{4}, \rm{R_{g}}$, respectively. We estimate an inclination angle of $81 \pm 5$ degrees for Swift J1357.2-0933 and an angle of $73 \pm 4$ degrees for MAXI J1305-704. This is entirely consistent with the values reported in the literature. In agreement with the recent literature, our analysis of the outburst spectrum of MAXI J1305-704 revealed that the radius of the emission region deviates from expected values. This outcome implies several potential scenarios, including alternative disc configuration or even a circumbinary disc. We caution that these results were derived from a simplistic model that may not fully describe the complicated physics of accretion discs. Despite these limitations, our results for the inclination angles are remarkably consistent with recent complementary studies, and the proposed description of the emitting region remains entirely plausible.

We review the advantages of fitting with a Two Component Advective Flow (TCAF) which uses only four physical parameters. We then present the results of hydrodynamic simulations to highlight the fact that the primary component of a black hole accretion remains the sub-Keplerian or the low angular momentum flow independent of whether we have a high, intermediate or low mass X-ray binary. Every aspect of spectral and timing properties, including the disk-jet connection could be understood well only if such a component is present along with a Keplerian component of variable size and accretion rate.

We confirm for the first time the existence of distinctive halo bias associated with the quadrupolar type of statistical anisotropy (SA) of the linear matter density field using cosmological $N$-body simulations. We find that the coefficient of the SA-induced bias for cluster-sized halos takes negative values and exhibits a decreasing trend with increasing halo mass. This results in the quadrupole halo power spectra in a statistically anisotropic universe being less amplified compared to the monopole spectra. The anisotropic feature in halo bias that we found presents a promising new tool for testing the hypothesis of a statistically anisotropic universe, with significant implications for the precise verification of anisotropic inflation scenarios and vector dark matter and dark energy models.

Accurately characterizing the true redshift (true-$z$) distribution of a photometric redshift (photo-$z$) sample is critical for cosmological analyses in imaging surveys. Clustering-based techniques, which include clustering-redshift (CZ) and self-calibration (SC) methods--depending on whether external spectroscopic data are used--offer powerful tools for this purpose. In this study, we explore the joint inference of the true-$z$ distribution by combining SC and CZ (denoted as SC+CZ). We derive simple multiplicative update rules to perform the joint inference. By incorporating appropriate error weighting and an additional weighting function, our method shows significant improvement over previous algorithms. We validate our approach using a DES Y3 mock catalog. The true-$z$ distribution estimated through the combined SC+CZ method is generally more accurate than using SC or CZ alone. To account for the different constraining powers of these methods, we assign distinct weights to the SC and CZ contributions. The optimal weights, which minimize the distribution error, depend on the relative constraining strength of the SC and CZ data. Specifically, for a spectroscopic redshift sample that represents 1% of the photo-$z$ sample, the optimal combination reduces the total error by 20% (40%) compared to using CZ (SC) alone, and it keeps the bias in mean redshift [$\Delta \bar{z} / (1 + z) $] at the level of 0.3%. Furthermore, when CZ data is only available in the low-$z$ range and the high-$z$ range relies solely on SC data, SC+CZ enables consistent estimation of the true-$z$ distribution across the entire redshift range. Our findings demonstrate that SC+CZ is an effective tool for constraining the true-$z$ distribution, paving the way for clustering-based methods to be applied at $z\gtrsim 1$.

Faraday rotation is an important probe of the magnetic fields and magnetized plasma around active galactic nuclei (AGN) jets. We present a Faraday rotation measure image of the M87 jet between 85.2 GHz and 101.3 GHz with a resolution of ~2" with the Atacama Large Millimeter/submillimeter Array (ALMA). We found that the rotation measure (RM) of the M87 core is $\rm (4.5\pm 0.4)\times10^{4}\ rad\ m^{-2}$ with a low linear polarization fraction of $\rm (0.88\pm 0.08)\%$. The spatial RM gradient in the M87 jet spans a wide range from $\sim -2\times10^4\rm~rad\ m^{-2}$ to $\sim 3\times10^4\rm~rad\ m^{-2}$ with a typical uncertainty of $0.3\times10^4\rm~rad\ m^{-2}$. A comparison with previous RM measurements of the core suggests that the Faraday rotation of the core may originate very close to the super massive black hole (SMBH). Both an internal origin and an external screen with a rapidly varying emitting source could be possible. As for the jet, the RM gradient indicates a helical configuration of the magnetic field that persists up to kpc scale. Combined with the kpc-scale RM measurements at lower frequencies, we found that RM is frequency-dependent in the jet. One possible scenario to explain this dependence is that the kpc-scale jet has a trumpet-like shape and the jet coil unwinds near its end.

We characterize blue straggler stars (BSS) and yellow straggler stars (YSS) of an open cluster (OC) Berkeley 39 using multi-wavelength observations including Swift/UVOT. Our analysis also makes use of ultraviolet (UV) data from GALEX, optical data from Gaia DR3 and Pan-STARRS, and infrared data from 2MASS, Spitzer/IRAC, and WISE. Berkeley 39 is a ~6 Gyr old Galactic OC located at a distance of ~4200 pc. We identify 729 sources as cluster members utilizing a machine learning algorithm, ML-MOC, on Gaia DR3 data. Of these, 17 sources are classified as BSS candidates and four as YSS candidates. We construct multi-wavelength spectral energy distributions (SEDs) of 16 BSS and 2 YSS candidates, within the Swift/UVOT field, to analyze their properties. Out of these, 8 BSS candidates and both the YSS candidates are successfully fitted with single-component SEDs. Five BSS candidates show marginal excess in the near-UV (fractional residual < 0.3 in all but one UVOT filter), whereas three BSS candidates show moderate to significant excess in the near-UV (fractional residual > 0.3 in at least two UVOT filters). We present the properties of the BSS and YSS candidates, estimated based on the SED fits.

Accurate modeling of tidal interactions is crucial for interpreting recent JWST observations of the thermal emissions of TRAPPIST-1~b and c and for characterizing the surface conditions and potential habitability of the other planets in the system. Indeed, the rotation state of the planets, driven by tidal forces, significantly influences the heat redistribution regime. Due to their proximity to their host star and the estimated age of the system, the TRAPPIST-1 planets are commonly assumed to be in a synchronization state. In this work, we present the recent implementation of the co-planar tidal torque and forces equations within the formalism of Kaula in the N-body code Posidonius. This enables us to explore the hypothesis of synchronization using a tidal model well suited to rocky planets. We studied the rotational state of each planet by taking into account their multi-layer internal structure computed with the code Burnman. Simulations show that the TRAPPIST-1 planets are not perfectly synchronized but oscillate around the synchronization state. Planet-planet interactions lead to strong variations on the mean motion and tides fail to keep the spin synchronized with respect to the mean motion. As a result, the sub-stellar point of each planet experiences short oscillations and long-timescale drifts that lead the planets to achieve a synodic day with periods varying from $55$~years to $290$~years depending on the planet.

Nowadays, transit timing variations (TTVs) are proving to be a very valuable tool in exoplanetary science to detect exoplanets by observing variations in transit times. To study the transit timing variation of the hot Jupiter, TrES-2b, we have combined 64 high-quality transit light curves from all seven sectors of NASA's Transiting Exoplanet Survey Satellite (TESS) along with 60 best-quality light curves from the ground-based facility Exoplanet Transit Database (ETD) and 106 mid-transit times from the previous works. From the precise transit timing analysis, we have observed a significant improvement in the orbital ephemerides, but we did not detect any short period TTVs that might result from an additional body. The inability to detect short-term TTVs further motivates us to investigate long-term TTVs, which might be caused by orbital decay, apsidal precession, Applegate mechanism, and $R{\phi}$mer effect and the orbital decay appeared to be a better explanation for the observed TTV with $\Delta BIC$ = 4.32. The orbital period of the hot Jupiter TrES-2b appears to be shrinking at a rate of $-5.58 \pm 1.81$ ms/yr. Assuming this decay is primarily caused by tidal dissipation within the host star, we have subsequently calculated the stellar tidal quality factor value to be 9900, which is 2 to 3 orders of magnitude smaller than the theoretically predicted values for other hot-Jupiter systems and its low value indicates more efficient tidal dissipation within the host star. Additional precise photometric and radial velocity observations are required to pinpoint the cause of the change in the orbital period.

Context. Recent observations of decayless transverse oscillations have shown two branches in the relationship between periods and loop lengths. One is a linear relationship, interpreted as a standing mode. The other shows almost no correlation and has not yet been interpreted conclusively. Aims. We investigated the undersampling effect on observed periods of decayless oscillations. Methods. We considered oscillating coronal loops that closely follow the observed loop length distribution. Assuming that all oscillations are standing waves, we modeled a signal that represents decayless oscillations where the period is proportional to the loop length and the amplitude and phase are randomly drawn. A downsampled signal was generated from the original signal by considering different sample rates that mimic temporal cadences of telescopes, and periods for sampled signals were analysed using the fast Fourier transform. Results. When the sampling cadence is getting closer to the actual oscillation period, a tendency for overestimating periods in short loops is enhanced. The relationship between loop lengths and periods of the sampled signals shows the two branches as in the observation. Conclusions. We find that long periods of decayless oscillations occurring in short loops could be the result of undersampling.

Spectral distortions of the cosmic microwave background (CMB) provide stringent constraints on energy and entropy production in the post-BBN (Big Bang Nucleosynthesis) era. This has been used to constrain dark photon models with COBE/FIRAS and forecast the potential gains with future CMB spectrometers. Here, we revisit these constraints by carefully considering the photon to dark photon conversion process and evolution of the distortion signal. Previous works only included the effect of CMB energy density changes but neglected the change to the photon number density. We clearly define the dark photon distortion signal and show that in contrast to previous analytic estimates the distortion has an opposite sign and a $\simeq 1.5$ times larger amplitude. We furthermore extend the treatment into the large distortion regime to also cover the redshift range $\simeq 2\times 10^6-4\times 10^7$ between the $\mu$-era and the end of BBN using CosmoTherm. This shows that the CMB distortion constraints for dark photon masses in the range $10^{-4}\,{\rm eV}\lesssim m_{\rm dp}\lesssim 10^{-3}\,{\rm eV}$ were significantly underestimated. We demonstrate that in the small distortion regime the distortion caused by photon to dark photon conversion is extremely close to a $\mu$-type distortion independent of the conversion redshift. This opens the possibility to study dark photon models using CMB distortion anisotropies and the correlations with CMB temperature anisotropies as we highlight here.

One of the fundamental questions of cosmology is the origin and mechanism(s) responsible for the reionization of the Universe beyond $z\sim6$. To address this question, many studies over the past decade have focused on local ($z\sim0.3$) galaxies which leak ionizing radiation (Lyman continuum or LyC). However, line-of-sight effects and data quality have prohibited deeper insight into the nature of LyC escape. To circumvent these limitations, we analyze stacks of a consolidated sample of {\it HST}/COS observations of the LyC in 89 galaxies at $z\sim0.3$. From fitting of the continuum, we obtain information about the underlying stellar populations and neutral ISM geometry. We find that most LyC non-detections are not leaking appreciable LyC ($f_{esc}^{\rm LyC}<1$\%) but also that exceptional cases point to spatial variations in the LyC escape fraction $f_{esc}^{\rm LyC}$. Stellar populations younger than 3 Myr lead to an increase in ionizing feedback, which in turn increases the isotropy of LyC escape. Moreover, mechanical feedback from supernovae in 8-10 Myr stellar populations is important for anisotropic gas distributions needed for LyC escape. While mechanical feedback is necessary for any LyC escape, high $f_{esc}^{\rm LyC}$ ($>5$\%) also requires a confluence of young stars and ionizing feedback. A two-stage burst of star formation could facilitate this optimal LyC escape scenario.

We present a weak lensing analysis of the galaxy cluster Abell 2390 at z = 0.23 using second moment shape measurements made in 411 short 60s exposures. The exposures are obtained in three broadband photometric filters (g, r, i) using WIYN-ODI. Shape measurement in individual exposures is done using a moment matching algorithm. Forced measurement is used when the moment matching algorithm fails to converge at low signal to noise ratio (SNR). The measurements made in individual images are combined using inverse error weight to obtain accurate shape of sources and hence recover shear. We use PhoSim simulations to validate shear measurements recovered by our pipeline. We find the mass of Abell 2390 is in agreement with previously published results. We also find the E-Mode maps show filamentary structures consistent with baryonic structures and recovers most clusters/groups of galaxies found using Optical and X-Ray data. Thus we demonstrate the feasibility of using Weak Lensing to map large scale structure of the universe. We also find the central portion of the cluster has a bimodal mass distribution and the relative orientation of the peaks are similar to X-Ray. We discuss earlier research on this galaxy cluster and show that a late stage merger accounts for all the observed data.

Oscillatory reconnection is a specific type of time-dependent reconnection which involves periodic changes in the magnetic topology of a null point. The mechanism has been reported for a variety of magnetic field strengths and configurations, background temperatures and densities. All these studies report an oscillation in the current density at the null point, but also report a variety of periods, amplitudes and overall behaviors. We conduct a parametric study for equilibrium magnetic field strength and initial background temperature, solving 2D resistive MHD equations around a magnetic X-point. We introduce a parameter space for the ratio of internal-to-magnetic energy and find self-similar solutions for simulations where this ratio is below 0.1 (which represents a magnetically-dominated environment or, equivalently, a low-beta plasma). Self-similarity can be seen in oscillations in the current density at the null (including amplitude and period), Ohmic heating and the temperature generated via reconnection jets. The parameter space of energy ratios also allows us to contextualize previous studies of the oscillatory reconnection mechanism and bring those different studies together into a single unified understanding.

Ultra-luminous X-ray sources (ULXs) are high-mass X-ray binaries with an X-ray luminosity above $10^{39}$ erg s$^{-1}$. These ULXs can be powered by black holes that are more massive than $20M_\odot$, accreting in a standard regime, or lighter compact objects accreting supercritically. There are only a few ULXs with known optical or UV counterparts, and their nature is debated. Determining whether optical/UV radiation is produced by the donor star or by the accretion disc is crucial for understanding ULX physics and testing massive binary evolution. We conduct, for the first time, a fully consistent multi-wavelength spectral analysis of a ULX and its circumstellar nebula. We aim to establish the donor star type and test the presence of strong disc winds in the prototypical ULX Holmberg II X-1 (Ho II X-1). We intent to obtain a realistic spectral energy distribution of the ionising source, which is needed for robust nebula analysis. We acquired new UV spectra of Ho II X-1 with the HST and complemented them with archival optical and X-ray data. We explored the spectral energy distribution of the source and analysed the spectra using the stellar atmosphere code PoWR and the photoionisation code Cloudy. Our analysis of the X-ray, UV, and optical spectra of Ho II X-1 and its nebula consistently explains the observations. We do not find traces of disc wind signatures in the UV and the optical, rejecting previous claims of the ULX being a supercritical accretor. The optical/UV counterpart of HoII X-1 is explained by a B-type supergiant donor star. Thus, the observations are fully compatible with Ho II X-1 being a close binary consisting of an $\gtrsim 66\,M_\odot$ black hole accreting matter from an $\simeq 22 M_\odot$ B-supergiant companion. Also, we propose a possible evolution scenario for the system, suggesting that Ho II X-1 is a potential gravitational wave source progenitor.

We introduce $\texttt{A}$strophysical $\texttt{H}$ybrid-$\texttt{K}$inetic simulations with the $\texttt{flASH}$ code ($\texttt{AHKASH}$) -- a new Hybrid particle-in-cell (PIC) code developed within the framework of the multi-physics code $\texttt{FLASH}$. The new code uses a second-order accurate Boris integrator and a predictor-predictor-corrector algorithm for advancing the Hybrid-kinetic equations, using the constraint transport method to ensure that magnetic fields are divergence-free. The code supports various interpolation schemes between the particles and grid cells, with post-interpolation smoothing to reduce finite particle noise. We further implement a $\delta f$ method to study instabilities in weakly collisional plasmas. The new code is tested on standard physical problems such as the motion of charged particles in uniform and spatially varying magnetic fields, the propagation of Alfvén and whistler waves, and Landau damping of ion acoustic waves. We test different interpolation kernels and demonstrate the necessity of performing post-interpolation smoothing. We couple the $\texttt{TurbGen}$ turbulence driving module to the new Hybrid PIC code, allowing us to test the code on the highly complex physical problem of the turbulent dynamo. To investigate steady-state turbulence with a fixed sonic Mach number, it is important to maintain isothermal plasma conditions. Therefore, we introduce a novel cooling method for Hybrid PIC codes and provide tests and calibrations of this method to keep the plasma isothermal. We describe and test the `hybrid precision' method, which significantly reduces (by a factor $\sim1.5$) the computational cost, without compromising the accuracy of the numerical solutions. Finally, we test the parallel scalability of the new code, showing excellent scaling up to 10,000~cores.

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide an unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium. We present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of 1 < p < 2 and 2 $\leq$ p. We show that the SSC scenario plays a relevant role in the radiative parameter $\epsilon$, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies $\geq$ 100 MeV in a sample of nine bursts from the second Fermi/LAT GRB catalog that exhibited temporal and spectral indices with $\geq$ 1.5 and $\approx$ 2, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modeled the multi-wavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor and the mechanism responsible for explaining the energetic GeV photons.

The Milky Way's Central Molecular Zone (CMZ) differs dramatically from our local solar neighbourhood, both in the extreme interstellar medium conditions it exhibits (e.g. high gas, stellar, and feedback density) and in the strong dynamics at play (e.g. due to shear and gas influx along the bar). Consequently, it is likely that there are large-scale physical structures within the CMZ that cannot form elsewhere in the Milky Way. In this paper, we present new results from the Atacama Large Millimeter/submillimeter Array (ALMA) large programme ACES (ALMA CMZ Exploration Survey) and conduct a multi-wavelength and kinematic analysis to determine the origin of the M0.8$-$0.2 ring, a molecular cloud with a distinct ring-like morphology. We estimate the projected inner and outer radii of the M0.8$-$0.2 ring to be 79" and 154", respectively (3.1 pc and 6.1 pc at an assumed Galactic Centre distance of 8.2 kpc) and calculate a mean gas density $> 10^{4}$ cm$^{-3}$, a mass of $\sim$ $10^6$ M$_\odot$, and an expansion speed of $\sim$ 20 km s$^{-1}$, resulting in a high estimated kinetic energy ($> 10^{51}$ erg) and momentum ($> 10^7$ M$_\odot$ km s$^{-1}$). We discuss several possible causes for the existence and expansion of the structure, including stellar feedback and large-scale dynamics. We propose that the most likely cause of the M0.8$-$0.2 ring is a single high-energy hypernova explosion. To viably explain the observed morphology and kinematics, such an explosion would need to have taken place inside a dense, very massive molecular cloud, the remnants of which we now see as the M0.8$-$0.2 ring. In this case, the structure provides an extreme example of how supernovae can affect molecular clouds.

Accretion among planets is a poorly understood phenomenon, due to lack of both observational and theoretical studies. Detection of emission lines from accreting gas giants facilitate detailed investigations into this process. This work presents a detailed analysis of Balmer lines from one of the few known young, planetary-mass objects with observed emission, the isolated L2 dwarf 2MASS J11151597+1937266 with a mass 7-21 Mj and age 5-45 Myr, located at 45+-2 pc. We obtained the first high-resolution (R~50,000) spectrum of the target with VLT/UVES, a spectrograph in the near-UV to visible wavelengths (3200-6800 AA). We report resolved H3-H6 and He I (5875.6 AA) emission in the spectrum. Based on the asymmetric line profiles of H3 and H4, 10% width of H3 (199+-1 km/s), tentative He I 6678 AA emission and indications of a disk from MIR excess, we confirm ongoing accretion at this object. Using the Gaia update of the parallax, we revise its temperature to 1816+-63 K and radius to 1.5+-0.1 Rj. Analysis of observed H I profiles using 1D planet-surface shock model implies a pre-shock gas velocity of v0=120(+80,-40) km/s and a pre-shock density of log(n0/cm^-3)=14(+0,-5). Pre-shock velocity points to a mass of 6(+8,-4) Mj for the target. Combining the H I line luminosities and planetary Lline-Lacc scaling relations, we derive a mass accretion rate of 1.4(+2.8,-0.9)x10^-8 Mj/yr.

Deviations of the cosmic microwave background (CMB) energy spectrum from a perfect blackbody uniquely probe a wide range of physics, ranging from fundamental physics in the primordial Universe ($\mu$-distortion) to late-time baryonic feedback processes (y-distortion). While the y-distortion can be detected with a moderate increase in sensitivity over that of COBE/FIRAS, the $\Lambda$CDM-predicted $\mu$-distortion is roughly two orders of magnitude smaller and requires substantial improvements, with foregrounds presenting a serious obstacle. Within the standard model, the dominant contribution to $\mu$ arises from energy injected via Silk damping, yielding sensitivity to the primordial power spectrum at wavenumbers $k \approx 1-10^{4}$ Mpc$^{-1}$. Here, we present a new instrument concept, SPECTER, with the goal of robustly detecting $\mu$. The instrument technology is similar to that of LiteBIRD, but with an absolute temperature calibration system. Using a Fisher approach, we optimize the instrument's configuration to target $\mu$ while robustly marginalizing over foreground contaminants. Unlike Fourier-transform-spectrometer-based designs, the specific bands and their individual sensitivities can be independently set in this instrument, allowing significant flexibility. We forecast SPECTER to observe the $\Lambda$CDM-predicted $\mu$-distortion at $\approx 5\sigma$ (10$\sigma$) assuming an observation time of 1 (4) year(s) (corresponding to mission duration of 2 (8) years), after foreground marginalization. Our optimized configuration includes 16 bands spanning 1-2000 GHz with degree-scale angular resolution at 150 GHz and 1046 total detectors. SPECTER will additionally measure the y-distortion at sub-percent precision and its relativistic correction at percent-level precision, yielding tight constraints on the total thermal energy and mean temperature of ionized gas.