Papers for Thursday, Oct 05 2023

Recent space-based missions have ushered in a new era of observational astronomy, where high-cadence photometric light curves for thousands to millions of stars in the solar neighborhood can be used to test and apply stellar age-dating methods, including gyrochronology. Combined with precise kinematics, these data allow for powerful new insights into our understanding of the Milky Way's dynamical history. Using TESS data, we build a series of rotation period measurement and confirmation pipelines and test them on 1,560 stars across five benchmark samples: the Pleiades, Pisces--Eridanus, Praesepe, the Hyades, and field stars from the MEarth Project. Our pipelines' recovery rates across these groups are on average 89\%. We then apply these pipelines to 4,085 likely single stars with TESS light curves in two interesting regions of Galactic action space. We identify 141 unique, rapidly rotating stars in highly eccentric orbits in the disk, some of which appear as rotationally young as the 120-Myr-old Pleiades. Pending spectroscopic analysis to confirm their youth, this indicates these stars were subject to fast-acting dynamical phenomena, the origin of which will be investigated in later papers in this series.

We consider collective Thomson scattering of an incident X-mode wave (with the electric vector perpendicular to the background magnetic field) in magnetized electron and positron pair plasma. The collective effects do not exactly cancel out in contrast to the non-magnetized case. Still, the cross-section is comparable to the non-collective one, with the same suppression by the square of the cyclotron frequency in a strong magnetic field. The comparable cross-section holds even though the net current is nearly zero from the drift motion of electrons and positrons. The plasma response does not also affect the cross-section so much. The spectrum of the scattered wave in finite temperature plasma peaks at cyclotron overtones. Based on these results, we also estimate induced Compton scattering in strongly magnetized pair plasma. Implications for pulsars and fast radio bursts are discussed.

The advent of JWST (the James Webb Space Telescope) now allows entire star cluster populations to be imaged in galaxies at cosmologically significant redshifts, bringing with it the need to apply K-corrections to their magnitudes and colour indices. Since the stellar populations within star clusters can be well approximated by a single age and metallicity, their spectral energy distributions are very different from those of galaxies or supernovae, and their K-corrections behave differently. We derive the photometric K-corrections versus redshift for model star clusters that cover a wide range of ages and metallicities, illustrating the results particularly for the broadband filters on the HST/ACS and the JWST/NIRCam cameras that are most commonly being used for imaging of populations of star clusters in distant galaxies. In an Appendix, we introduce a simple webtool called RESCUER that can generate K-values for any user-defined combination of cluster properties.

We aim at a direct measurement of the compactness of three galaxy-scale lenses in massive clusters, testing the accuracy of the scaling laws that describe the members in strong lensing (SL) models of galaxy clusters. We selected the multiply imaged sources MACS J0416.1$-$2403 ID14 ($z=3.221$), MACS J0416.1$-$2403 ID16 ($z=2.095$), and MACS J1206.2$-$0847 ID14 ($z=3.753$). Eight images were observed for the first SL system, and six for the latter two. We focused on the main deflector of each galaxy-scale SL system (identified as members 8971, 8785, and 3910, respectively), and modelled its total mass distribution with a truncated isothermal sphere. We accounted for the lensing effects of the remaining cluster components, and included the uncertainty on the cluster-scale mass distribution through a bootstrapping procedure. We measured a truncation radius value of $6.1^{+2.3}_{-1.1} \, \mathrm{kpc}$, $4.0^{+0.6}_{-0.4} \, \mathrm{kpc}$, and $5.2^{+1.3}_{-1.1} \, \mathrm{kpc}$ for members 8971, 8785, and 3910, respectively. Alternative non-truncated models with a higher number of free parameters do not lead to an improved description of the SL system. We measured the stellar-to-total mass fraction within the effective radius $R_e$ for the three members, finding $0.51\pm0.21$, $1.0\pm0.4$, and $0.39\pm0.16$, respectively. We find that a parameterisation of the properties of cluster galaxies in SL models based on power-law scaling relations with respect to the total luminosity cannot accurately describe their compactness over their full total mass range. Our results agree with modelling of the cluster members based on the Fundamental Plane relation. Finally, we report good agreement between our values of the stellar-to-total mass fraction within $R_e$ and those of early-type galaxies from the SLACS Survey. Our work significantly extends the regime of the current samples of lens galaxies.

Theoretical models of spiral arms suggest that the spiral arms provoke a vertical bulk motion in disc stars. By analysing the breathing motion, a coherent asymmetric vertical motion around the mid-plane of the Milky Way disc, with $\textit{Gaia}$ DR3, we found that a compressing breathing motion presents along the Local arm. On the other hand, with an $N$-body simulation of an isolated Milky Way-like disc galaxy, we found that the transient and dynamic spiral arms induce compressing breathing motions when the arms are in the growth phase, while the expanding breathing motion appears in the disruption phase. The observed clear alignment of the compressing breathing motion with the Local arm is similar to what is seen in the growth phase of the simulated spiral arms. Hence, we suggest that the Local arm's compressing breathing motion can be explained by the Local arm being in the growth phase of a transient and dynamic spiral arm. We also identified the tentative signatures of the expanding breathing motion associated with the Perseus arm and also the Outer arm coinciding with the compressing breathing motion. This may infer that the Perseus and Outer arms are in the disruption and growth phases, respectively.

We use CANUCS JWST/NIRCam imaging of galaxies behind the gravitationally-lensing cluster MACS J0417.5-1154 to investigate star formation burstiness in low-mass ($M_\star\sim10^8\ M_\odot$) galaxies at $z\sim4.7-6.5$. Our sample of 123 galaxies is selected using the Lyman break selection and photometric emission-line excess methods. Sixty percent of the 123 galaxies in this sample have H$\alpha$-to-UV flux ratios that deviate significantly from the range of $\eta_{1500}$ values consistent with smooth and steady star formation histories. This large fraction indicates that the majority of low-mass galaxies is experiencing bursty star formation histories at high redshift. We also searched for interacting galaxies in our sample and found that they are remarkably common ($\sim40\ \%$ of the sample). Compared to non-interacting galaxies, interacting galaxies are more likely to have very low H$\alpha$-to-UV ratios, suggesting that galaxy-galaxy interactions enhance star formation burstiness and enable faster quenching (with timescales of $\lesssim100$ Myr) that follows the rapid rise of star formation activity. Given the high frequency of galaxy-galaxy interactions and the rapid SFR fluctuations they appear to cause, we conclude that galaxy-galaxy interactions could be a leading cause of bursty star formation in low-mass, high-$z$ galaxies. They could thus play a significant role in the evolution of the galaxy population at early cosmological times.

When measuring the observed pressure, density or temperature profiles of the intracluster gas, and hence the mass of clusters of galaxies, projection effects or departures from the spherical symmetry hypothesis may induce biases. In order to estimate how strongly the cluster observed properties depend on the direction of observation, we use a constrained hydrodynamical simulation of Virgo cluster that replicates the actual cluster of galaxies. In this case study, we analyse Virgo properties when they are projected along different directions including along the Milky Way-Virgo axis which mimics our observation direction. We compare the hydrostatic mass and the hydrostatic mass bias from the projection along the different observation directions to that derived from the 3D simulation. We show that projection effects impact the determination of Virgo mass. We particularly demonstrate that the mass and pressure along the line of sight correlates with the 2D and 3D-deprojected electron density and pressure profiles intensity and thus impacts the derived hydrostatic mass. We also show that the deviations to the hydrostatic equilibrium induced by pressure discontinuities within the cluster are emphasized by the deprojection process and thus makes the hydrostatic mass estimation invalid at these radii.

Multi-planet systems provide important laboratories for exploring dynamical interactions within the range of known exoplanetary system architectures. One such system is GJ 357, consisting of a low-mass host star and three orbiting planets, the outermost (planet d) of which does not transit but lies within the Habitable Zone (HZ) of the host star. The minimum mass of planet d causes its nature to be unknown, both in terms of whether it is truly terrestrial and if it is a candidate for harboring surface liquid water. Here, we use three sectors of photometry from the Transiting Exoplanet Survey Satellite (TESS) to show that planets c and d do not transit the host star, and therefore may have masses higher than the derived minimum masses. We present the results for a suite of dynamical simulations that inject an Earth-mass planet within the HZ of the system for three different orbital and mass configurations of planet d. These results show that planet d, rather than being a potentially habitable planet, is likely a source of significant orbital instability for other potential terrestrial planets within the HZ. We find that relatively small eccentricities of planet d cause a majority of the HZ to be unstable for an Earth-mass planet. These results highlight the importance of dynamical stability for systems that are prioritized in the context of planetary habitability.

The arrival of gravitational wave astronomy and a growing number of time-domain focused observatories are set to lead to a increasing number of detections of short gamma-ray bursts (GRBs) launched with a moderate inclination to Earth. Being nearby events, these are also prime candidates for very long-term follow-up campaigns and very-long-baseline interferometry (VLBI), which has implications for multi-messenger modelling, data analysis, and statistical inference methods applied to these sources. Here we present a comprehensive modelling update that directly incorporates into afterglowpy astrometric observations of the GRB position, Poissonian statistics for faint sources, and modelling of a trans-relativistic population of electrons. We use the revolutionary event GW170817 to demonstrate the impact of these extensions both for the best-fit physics parameters and model selection methods that assess the statistical significance of additional late-time emission components. By including in our analysis the latest Chandra X-ray observations of GRB 170817A, we find only weak evidence (less than two sigma) for a new emission component at late times, which makes for a slightly more natural fit to the centroid evolution and prediction for the external medium density.

We investigate the potential of a pure linear expansion for the rest-frame flux of a type Ia supernova light curve fitter based on the well known Spectral Adaptive Light Curve Template 2 (SALT2). We generate the expansion components by performing Principal Component Analysis (PCA) and Factor Analysis (FA) onto a representative training set. Then, we derive a Tripp-like expression for the distance modulus and fit the $\Lambda$CDM cosmological model on the Pantheon sample. The constraining power of the model, dubbed Pure Expansion Template for Supernovae (PETS), and SALT2 is evaluated and we found compatible results for $\Omega_{m0}$ and $\Omega_{\Lambda 0}$ within $68\%$ uncertainty between the two models, with PETS' fit parameters exhibiting non negligible linear correlations with SALT2' parameters. We find non negligible correlations between PETS's fit parameters and the supernovae host galaxies masses, while the Hubble Diagram residues show no correlation with fit parameters, redshift or host galaxy mass. The model nuisance parameters, $\alpha$ and $\beta$, are slighted correlated and we find evidence for redshift evolution for $\beta$. The intrinsic scatter, $\sigma_{\textrm{int}}$, shows a subtle redshift evolution that should be further investigated increasing the number of high redshift supernovae in the cosmology sample.

In the standard quasar model, the accretion disk obscuration is due to the canonical dusty torus. Here, we argue that a substantial part of the quasar obscuration can come from the interstellar medium (ISM) when the quasars are embedded in compact starbursts. We use an obscuration-unbiased sample of 578 infrared (IR) quasars at $z\approx 1-3$ and archival ALMA submillimeter host galaxy sizes to investigate the ISM contribution to the quasar obscuration. We calculate SFR and ISM column densities for the IR quasars and a control sample of submillimeter galaxies (SMGs) not hosting quasar activity and show that: (1) the quasar obscured fraction is constant up to $\rm SFR\approx 300 \: M_{\odot} \: yr^{-1}$, and then increases towards higher SFR, suggesting that the ISM obscuration plays a significant role in starburst host galaxies, and (2) at $\rm SFR\gtrsim 300 \: M_{\odot} \: yr^{-1}$, the SMGs and IR quasars have similarly compact submillimeter sizes ($R_{\rm e}\approx 0.5-3\rm \: kpc$) and, consequently, the ISM can heavily obscure the quasar, even reaching Compton-thick ($N_{\rm H}>10^{24} \rm \: cm^{-2}$) levels in extreme cases. Based on our results, we infer that $\approx 10-30\%$ of the IR quasars with $\rm SFR\gtrsim 300 \: M_{\odot} \: yr^{-1}$ are obscured solely by the ISM.

We have developed a new computer code, RELDAFNA, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite volume Godunov scheme using the full characteristic decomposition of the SRHD equations, to achieve second and third order accuracy in space (both PLM and PPM reconstruction). For time integration, we use the method of directional splitting with symmetrization, which is second order accurate in time. We have also implemented second and third order Runge-Kutta time integration scheme for comparison. In addition to the hydrodynamics solvers we have implemented approximate Riemann solvers along with an exact Riemann solver. We examine the ability of RELDAFNA to accurately simulate special relativistic flows efficiently in number of processors, computer memory and over all integration time. We show that a wide variety of test problems can be solved as accurately as solved by higher order programs, such as RAM, GENESIS, or PLUTO, but with a less number of variables kept in memory and computer calculations than most schemes, an advantage which is crucial for 3D high resolution simulations to be of practical use for scientific research in computational astrophysics. RELDAFNA has been tested in one, two and three dimensions and in Cartesian, cylindrical and spherical geometries. We present the ability of RELDAFNA to assist with the understanding of open questions in high energy astrophysics which involve relativistic flows.

For temperate exoplanets orbiting M dwarf hosts, the proximity of the habitable zone to the star necessitates careful consideration of tidal effects. Spin synchronization of the planetary orbital period and rotation period is a common assumption for habitable zone planets across the entire M spectral type. This predicted tidal locking has important implications for conditions on the surface. In this manuscript, we investigate the plausibility of capture into Cassini State 2 for a known sample of 280 multiplanet systems orbiting M dwarf hosts. This resonance of the spin precession and orbital precession frequencies is capable of exciting planets into stable nonzero rotational obliquities. A large axial tilt resulting from this predicted resonance capture can preclude synchronous rotation, inducing some version of "day" and "night." Considering each planetary pair and estimating the spin and orbital precession frequencies, we report that 75% of detected planets orbiting M dwarfs may be plausibly excited to a high obliquity over long timescales. This is consistent with similar findings for planets orbiting close-in to FGK dwarfs. However, it is only for M dwarf planets that the parameter space relevant for capture into Cassini State 2 overlaps with the stellar habitable zone. This effect is strongest for host stars with effective temperatures $T_{eff}<3000$ K, where more than half of planets with $T_{eq}<400$ K could possess non-zero obliquity due to residence in Cassini State 2. This overlap renders the potential capture into Cassini States extremely relevant to understanding the galaxy's most common temperate planets.

We present a quantitative analysis of the properties of galaxies and structures evolving in universes dominated by different modified gravitational models, including two variants of the f(R)-gravity (F) and two of the Dvali-Gabdadze-Poratti (N) braneworld model, which respectively feature the chameleon and Vainshtein screening mechanisms. Using the Simulation HYdrodynamics BeyONd Einstein (SHYBONE) cosmological hydrodynamical full-physics simulations suite, we study the departures in the properties of galaxies residing in different environments with respect to the standard model (GR). Using two different criteria to compare, we find that structures formed within modified gravity tend to show a denser gas density profile than their GR counterparts. Within the different modified gravity models, N1 and F5 gravity models show greater departures from the standard model, with gas density profiles $\rho_{\rm IGM} \geq 30\%$ denser in the outskirts for the N1 model, and in the inner parts for the F5 model. Additionally, we find that haloes evolving in MG universes show, in general, larger quenched fractions than GR, reaching up to $20\%$ larger quenching fractions in F5 regardless of the stellar mass of the galaxy. With respect to the other models, F6, N1 and N5 show slightly larger quenched fractions, but no strong differences can be found. These results directly impact the colour distribution of galaxies, making them in MG models redder and older than their GR counterparts. Like GR, once the environment starts to play a role, galaxies rapidly get quenched and the differences between models vanish.

Primordial non-Gaussianities (PNGs) are signatures in the density field that encode particle physics processes from the inflationary epoch. Such signatures have been extensively studied using the Cosmic Microwave Background, through constraining the amplitudes, $f^{X}_{\rm NL}$, with future improvements expected from large-scale structure surveys; specifically, the galaxy correlation functions. We show that weak lensing fields can be used to achieve competitive and complementary constraints. This is shown via the new Ulagam suite of N-body simulations, a subset of which evolves primordial fields with four types of PNGs. We create full-sky lensing maps and estimate the Fisher information from three summary statistics measured on the maps: the moments, the cumulative distribution function, and the 3-point correlation function. We find that the year 10 sample from the Rubin Observatory Legacy Survey of Space and Time (LSST) can constrain PNGs to $\sigma(f^{\rm\,eq}_{\rm NL}) \approx 110$, $\sigma(f^{\rm\,or,lss}_{\rm NL}) \approx 120$, $\sigma(f^{\rm\,loc}_{\rm NL}) \approx 40$. For the former two, this is better than or comparable to expected galaxy clustering-based constraints from the Dark Energy Spectroscopic Instrument (DESI). The PNG information in lensing fields is on non-linear scales and at low redshifts ($z \lesssim 1.25$), with a clear origin in the evolution history of massive halos. The constraining power degrades by $\sim\!\!60\%$ under scale cuts of $\gtrsim 20{\,\rm Mpc}$, showing there is still significant information on scales mostly insensitive to small-scale systematic effects (e.g. baryons). We publicly release the Ulagam suite to enable more survey-focused analyses.

The processes responsible for the assembly of cold and warm gas in early-type galaxies (ETGs) are not well-understood. We report on the multiwavelength properties of 15 non-central, nearby ($z \leq$ 0.00889) ETGs primarily through Multi-Unit Spectroscopic Explorer (MUSE) and Chandra X-ray observations, to address the origin of their multiphase gas. The MUSE data reveals 8/15 sources contain warm ionized gas traced by the H$\alpha$ emission line. The morphology of this gas is found to be filamentary in 3/8 sources, similar to that observed in many group and cluster-centered galaxies. All H$\alpha$ filamentary sources have X-ray luminosities exceeding the expected emission from the stellar population, suggesting the presence of diffuse hot gas which likely cooled to form the cooler phases. The morphology of the remaining 5/8 sources are rotating gas disks, not as commonly observed in higher mass systems. Chandra X-ray observations (when available) of the ETGs with rotating H$\alpha$ disks indicate that they are nearly void of hot gas. A mixture of stellar mass loss and external accretion was likely the dominant channel for the cool gas in NGC 4526 and NGC 4710. These ETGs show full kinematic alignment between their stars and gas, and are fast rotators. The H$\alpha$ features within NGC 4191 (clumpy, potentially star-forming ring), NGC 4643 and NGC 5507 (extended structures) along with loosely overlapping stellar and gas populations allow us to attribute external accretion to be the primary formation channel of the cool gas in these systems.

We present SOFIA observations with HAWC+ and FIFI-LS of the peculiar galaxy Arp 25, also known as NGC 2276 or UGC 3740, whose morphology is deformed by its impact with the intra-group medium of the NGC 2300 galaxy group. These observations show the first direct proof of the enhancement of [CII] emission due to shocks caused by ram pressure in a group of galaxies. By comparing the [CII] emission to UV attenuation, dust emission, PAH, and CO emission in different regions of the galaxy, we find a clear excess of [CII] emission along the impact front with the intra-group medium. We estimate that the shock due to the impact with the intra-group medium increases the [CII] emission along the shock front by 60% and the global [CII] emission by approximately 25% with respect to the predicted [CII] emission assuming only excitation caused by stellar radiation. This result shows the danger of interpreting [CII] emission as directly related to star formation since shocks and other mechanisms can significantly contribute to the total [CII] emission from galaxies in groups and clusters.

Dense stellar clusters are expected to house the ideal conditions for binary black hole (BBH) formation, both through binary stellar evolution and through dynamical encounters. We use theoretical arguments as well as $N$-body simulations to make predictions for the evolution of BBHs formed through stellar evolution inside clusters from the cluster birth (which we term $\textbf{primordial binaries}$), and for the sub-population of merging BBHs. We identify three key populations: (i) BBHs that form in the cluster, and merge before experiencing any $\textit{strong}$ dynamical interaction; (ii) binaries that are ejected from the cluster after only one dynamical interaction; and, (iii) BBHs that experience more than one strong interaction inside the cluster. We find that populations (i) and (ii) are the dominant source of all BBH mergers formed in clusters with escape velocity $v_{\mathrm{esc}}\leq 30$ $\mathrm{km\,s^{-1}}$. At higher escape velocities, dynamics are predicted to play a major role both for the formation and subsequent evolution of BBHs. Finally, we argue that for sub-Solar metallicity clusters with $v_{\mathrm{esc}}\lesssim100$ $\mathrm{km\,s^{-1}}$, the dominant form of interaction experienced by primordial BBHs (BBHs formed from primordial binaries) within the cluster is with other BBHs. The complexity of these binary-binary interactions will complicate the future evolution of the BBH and influence the total number of mergers produced.

Liquid water oceans are thought to underlie the ice shells of Europa and Enceladus. However, ocean properties can be challenging to measure due to the overlying ice cover. Here, we explore how measurements of ice shell thickness, which may be easier to obtain, could be used to infer information about the subsurface ocean. In particular, we consider lateral gravity-driven flow of the ice shells of icy satellites and relate this to ocean freeze and melt rates. We employ a first-principles approach applicable to conductive ice shells. We derive a scaling law under which ocean freeze and melt rates can be estimated from thickness measurements of a shell with a vertically-varying temperature-dependent viscosity. Under a steady-state assumption, ocean freeze and melt rates can be inferred from measurements of ice thickness; however, these rates depend on the basal viscosity, a key unknown. Depending on a characteristic thickness scale and basal viscosity, the characteristic freeze and melt rates range from about O(10$^{-1}$) to O(10$^{-5}$) mm/year. We validate our scaling in an Earth environment with ice-penetrating radar measurements of ice thickness and modelled snow accumulation for Roosevelt Island, Antarctica. Our model, coupled with the forthcoming observations of shell thickness from upcoming missions, may help bound the magnitudes of estimated ocean freeze and melt rates on icy satellites and shed light on potential ocean stratifications.

We conduct one-dimensional stellar evolutionary numerical simulations under the assumption that an efficient dynamo operates in the core of massive stars years to months before core collapse, and find that the magnetic activity enhances mass loss rate and might trigger binary interaction that leads to outbursts. We assume that the magnetic flux tubes that the dynamo forms in the inner core buoy out to the outer core where there is a steep entropy rise and a molecular weight drop. There the magnetic fields turn to thermal energy, i.e., by reconnection. We simulate this energy deposition where the entropy steeply rises and find that for our simulated cases the envelope radius increases by a factor of ~1.2-2 and luminosity by about an order of magnitude. These changes enhance mass loss rate. The envelope expansion can trigger a binary interaction that power an outburst. Because magnetic field amplification depends positively on the core rotation rate and operates in cycles, not in all cases the magnetic activity will be powerful enough to change envelope properties. Namely, only a fraction of core collapse supernovae experience pre-explosion outbursts.

The trajectories of coronal mass ejections (CMEs) are often seen to substantially deviate from a purely radial propagation direction. Such deviations occur predominantly in the corona and have been attributed to "channeling" or deflection of the eruptive flux by asymmetric ambient magnetic fields. Here, we investigate an additional mechanism that does not require any asymmetry of the pre-eruptive ambient field. Using magnetohydrodynamic numerical simulations, we show that the trajectory of CMEs through the solar corona can significantly deviate from a radial direction when propagation takes place in a unipolar radial field. We demonstrate that the deviation is most prominent below ~15 solar radii and can be attributed to an "effective IxB force" that arises from the intrusion of a magnetic flux rope with a net axial electric current into a unipolar background field. These results are important for predictions of CME trajectories in the context of space weather forecasts, as well as for reaching a deeper understanding of the fundamental physics underlying CME interactions with the ambient fields in the extended solar corona.

Context. Stars tend to form in clusters, but many escape their birth clusters very early. Identifying the escaped members of clusters can inform us about the dissolution of star clusters, but also about the stellar dynamics in the galaxy. Methods capable of finding escaped stars from many clusters are required to fully exploit the large amounts of data in the Gaia era. Aims. We present a new method of identifying escaped members of nearby clusters and apply it to ten young clusters. Methods. We assumed the escaped stars were close to the cluster in the past and performed traceback computations based on the Gaia DR3 radial velocity subsample. For each individual star, our method produces a probability estimate that it is an escaped member of a cluster, and for each cluster it also estimates the field star contamination rate of the identified fugitives. Results. Our method is capable of finding fugitives that have escaped from their cluster in the last few ten million years. In many cases the fugitives form an elongated structure that covers a large volume. Conclusions. The results presented here show that traceback computations using Gaia DR3 data can identify stars that have recently escaped their cluster. Our method will be even more useful when applied to future Gaia data releases that contain more radial velocity measurements.

We conducted the time-resolved simultaneous optical spectroscopic and photometric observations of mid M dwarf flare stars YZ CMi, EV Lac, and AD Leo. Spectroscopic observations were obtained using Apache Point Observatory 3.5m and Small \& Moderate Aperture Research Telescope System 1.5m telescopes during 31 nights. Among the 41 detected flares, seven flares showed clear blue wing asymmetries in the H$\alpha$ line, with various correspondences in flare properties. The duration of the blue wing asymmetries range from 20 min to 2.5 hours, including a flare showing the shift from blue to red wing asymmetry. Blue wing asymmetries can be observed during both white-light and candidate non white-light flares. All of the seven flares showed blue wing asymmetries also in the H$\beta$ line, but there are large varieties on which other chromospheric lines showed blue wing asymmetries. One among the 7 flares was also observed with soft X-ray spectroscopy, which enabled us to estimate the flare magnetic field and length of the flare loop. The line-of-sight velocities of the blue-shifted components range from -73 to -122 km s$^{-1}$. Assuming that the blue-shifts were caused by prominence eruptions, the mass of upward moving plasma was estimated to be 10$^{15}$ -- 10$^{19}$ g, which are roughly on the relation between flare energy and erupting mass expected from solar coronal mass ejections (CMEs). Although further investigations are necessary for understanding the observed various properties, these possible prominence eruptions on M-dwarfs could evolve into CMEs, assuming the similar acceleration mechanism with solar eruptions.

Considering subglacial liquid water, a significant extension of the classical Habitable Zone is obtained. Elaborating on the model of Wandel (2023) it is shown how an atmosphere and liquid water could survive on tidally locked planets closely orbiting an M-dwarf host, extending the Habitable Zone boundary inwards. In addition, subglacial liquid water could extend the Habitable Zone beyond the outer boundary of the conservative Habitable Zone as well. These two results enhance the circumstellar region with a potential for liquid water well beyond the conservative boundaries of the classical Habitable Zone. It is argued that the probable recent JWST detection of atmospheric water vapor on the rocky Earth-sized exoplanet GJ 486 b, along with earlier detections of water on other planets orbiting M-dwarf stars gives an empirical answer to the much-argued question, of whether such planets can support liquid water, organic chemistry and eventually life. It is shown how water on terrestrial planets closely orbiting M-dwarf stars may sustain in a subglacial melting layer. Finally, the model is applied to a few exoplanets demonstrating how water detection may constrain their atmospheric properties.

The Local Group (LG) of galaxies is dominated by the M31 and Milky Way (MW) pair, a configuration which suggests that the mass distribution in the LG and its surroundings should be highly anisotropic. We use the APOSTLE cosmological simulations to examine how this anisotropy manifests on the spatial distribution and kinematics of dwarf galaxies out to a distance of 3 Mpc from the MW. The simulations indicate a clear preference for dwarfs to be located close to the axis defined by the MW-M31 direction, even for dwarfs in the LG periphery (LGP; i.e., those at distances 1.25<d/Mpc<3). The LGP "Hubble flow" is also affected; at fixed distance from the MW the mean recession speed, <V_{rad} >, varies with angular distance to M31, peaking in the anti-M31 direction and reaching a minimum behind M31. The combined M31-MW mass decelerates the local expansion; the LG "turnaround radius" (where <V_{rad}>=0) in APOSTLE is located at r ~ 1.25 Mpc from the LG barycentre and the pure Hubble flow (where <V_{rad}> ~ H_0*d) is not reached out to at least 3 Mpc. The predicted flow is very cold, with a barycentric dispersion of <40 km/s. A comparison of these predicted features with existing observations gives mixed results. There is clear observational evidence for an angular anisotropy in V_{rad} around the LGP, but little evidence for a preferred spatial distribution of LGP dwarfs along the MW-M31 axis. The observed local Hubble flow is also peculiar. Although the coldness of the flow is consistent with the simulations, it is significantly less decelerated: relative to the MW, on average, all galaxies beyond d~1.25 Mpc seem to be on a pure Hubble flow. We argue that these oddities may result from incompleteness and inhomogeneous sky coverage, but a full explanation may need to await the completion of deep all-sky surveys able to fill the gaps in our current inventory of nearby dwarfs.

We follow up our photometric study of the post-outburst evolution of the FU Ori object V960 Mon with a complementary spectroscopic study at high dispersion that uses time series spectra from Keck/HIRES. Consistent with the photometric results reported in Carvalho et al. 2023, we find that the spectral evolution of V960 Mon corresponds to a decrease in the temperature of the inner disk, driven by a combination of decreasing accretion rate and increasing inner disk radius. We also find that although the majority of the absorption lines are well-matched by our accretion disk model spectrum, there are several strong absorption line families and a few emission lines that are not captured by the model. By subtracting the accretion disk model from the data at each epoch, we isolate the wind/outflow components of the system. The residuals show both broad and highly blueshifted profiles, as well as narrow and only slightly blueshifted profiles, with some lines displaying both types of features.

Blazars are the most abundant type of extragalactic gamma-ray source, usually presenting high variability across the electromagnetic spectrum. Their Very High Energy (VHE, above 0.1 TeV) emission has been studied in detail using Air Cherenkov Telescopes, with observations biased to flaring periods while their average activity has not been properly characterized. In this work, we report the results of 2090 days of quasi-continuous observations of the blazars 1ES 1215+303 and VER J0521+211, carried out with the High Altitude Water Cherenkov (HAWC) gamma-ray observatory. Fitting a power-law attenuated by photon-photon interaction with the extragalactic background light, we obtained a 6.2 $\sigma$ level detection for 1ES 1215+303 and a 4.3 $\sigma$ marginal detection for VER J0521+211. With the inclusion of the HAWC TeV spectrum, we built quasi-simultaneous multiwavelength spectral distributions and fitted a leptonic emission model to the observed data.

We present an update of the survey of Active Galaxies with the High Altitude Water Cherenkov (HAWC) gamma-ray observatory. This work adds 567 days of HAWC data to the previously published survey, providing a refined analysis of an updated total exposure of 2090 days. The sample includes 138 nearby AGNs from the 3FHL catalog. We fit a modified power-law to their very high energy spectra, including the exponential attenuation caused by the Extragalactic Background Light. We found four sources with significant detections (above 5$\sigma$): the radio galaxy M87 and the BL Lac objects Mkn 421, Mkn 501 and 1ES 1215+303. We also report eight sources with a marginal detection (between 3$\sigma$ and 5$\sigma$) of which seven are classified as BL Lac objects and one as a radio galaxy.

Although neural-network-based emulators enable efficient parameter estimation in 21-cm cosmology, the accuracy of such constraints is poorly understood. We employ nested sampling to fit mock data of the global 21-cm signal and high-$z$ galaxy ultraviolet luminosity function (UVLF) and compare for the first time the emulated posteriors obtained using the global signal emulator ${\tt globalemu}$ to the `true' posteriors obtained using the full model on which the emulator is trained using ${\tt ARES}$. Of the eight model parameters we employ, four control the star formation efficiency (SFE), and thus can be constrained by UVLF data, while the remaining four control UV and X-ray photon production, and the minimum virial temperature of star-forming halos ($T_{\rm min}$), and thus are uniquely probed by reionization and 21-cm measurements. For noise levels of 50 and 250 mK in the 21-cm data being jointly-fit, the emulated and `true' posteriors are consistent to within $1\sigma$. However, at lower noise levels of 10 and 25 mK, ${\tt globalemu}$ overpredicts $T_{\rm min}$ and underpredicts $\gamma_{\rm lo}$, an SFE parameter, by $\approx3-4\sigma$, while the `true' ${\tt ARES}$ posteriors capture their fiducial values within $1\sigma$. We find that jointly-fitting the mock UVLF and 21-cm data significantly improves constraints on the SFE parameters by breaking degeneracies in the ${\tt ARES}$ parameter space. Our results demonstrate the astrophysical constraints that can be expected for global 21-cm experiments for a range of noise levels from pessimistic to optimistic, and also the potential for probing redshift evolution of SFE parameters by including UVLF data.

For future extremely large telescopes, error in extreme adaptive optics systems at small angular separations will be highly impacted by the lag time of the correction, which is typically on millisecond timescales; one solution is to apply a predictive correction to catch up with the system delay. Predictive control leads to significant RMS error reductions in simulation (on the order of 5-10x improvement in RMS error compared with a standard integral controller), but shows only modest improvement on-sky (less than 2x in RMS error). This performance limitation is likely impacted by elements of pseudo open loop (POL) reconstruction, which requires assumptions about the response of the deformable mirror and accuracy of the wavefront measurements that are difficult to verify in practice. In this work, we explore a closed-loop method for data-driven prediction using a reformulated empirical orthogonal functions (EOF). We examine the performance of the open and closed-loop methods in simulation on perfect systems and systems with an inaccurate understanding of the DM response.

The inner $500\rm pc$ in the galactic center is dense with stars and a dynamically hot environment. Here, we focus on wide binaries as a source of tidally or collisional interactions. These binaries were previously ignored as sources of binary interaction because they are too wide to have a close pericenter passage, or they will quickly become unbound (ionized) due to gravitational interactions with passing neighbors. However, we show that wide binaries tend to interact more frequently with neighboring stars due to their larger cross-section for gravitational impulse interactions. These interactions mainly torque the wide system, causing it to change its eccentricity. As a result, the eccentricity might be excited to sufficiently high values, causing the binary to interact at the pericenter. As a proof of concept, we present four channels of such interactions: binary main-sequence (MS), white-dwarf (WD) - MS, black hole - MS, and lastly, WD-WD. During Galaxy's lifetime, we predict tens of thousands of MS-MS interacting binaries that may form G2-like objects later appear younger than their environment. X-ray signatures and, perhaps, supernovae may result from thousands of WD-MS and WD-WD interacting binaries from this channel. Lastly, we predict a few hundred BH-MS interacting binaries at the inner $500$~pc.

Context. Multi-wavelength studies of galaxies and galactic nuclei allow us to build a relatively more complete picture of the interstellar medium (ISM), especially in the dusty regions of starburst galaxies. An understanding of the physical processes in nearby galaxies can assist in the study of more distant sources at higher redshifts, which cannot be resolved. Aims. We aimed to use observations presented in the first part of this series of papers to model the physical conditions of the ISM in the nuclear region of NGC 253, in order to obtain primary parameters such as gas densities and metallicities. From the created model we further calculated secondary parameters such as gas masses of the different phases, and estimated the fraction of [C II] 158 um from the different phases, which allowed us to probe the nuclear star-formation rate. Methods. To compare theory with our observations we used MULTIGRIS, a probabilistic tool that determines probabilities for certain ISM parameters from a grid of Cloudy models together with a set of spectroscopic lines. Results. We find that the hypothetical active galactic nucleus within NGC 253 has only a minor impact compared to the starburst on the heating of the ISM as probed by the observed lines. We characterise the ISM and obtain parameters such as a solar metallicity, a mean density of ~230cm-3 , an ionisation parameter of log U = -3, and an age of the nuclear cluster of ~2 Myr. Furthermore, we estimate the masses of the ionised (3.8 x 10^6 M_sol ), neutral atomic (9.1 x 10^6 M_sol ), and molecular (2.0 x 10^8 M_sol ) gas phases as well as the dust mass (1.8 x 10^6 M_sol ) in the nucleus of NGC 253.

According to the $\Lambda$ cold dark matter model of galaxy formation, the hierarchical assembly process is scale-free and interactions between galaxies in all mass ranges are expected. The effects of interactions between dwarf galaxies on their evolution are not well understood. In this study, we aim to understand the effect of low-mass galaxy interactions on their star formation rate (SFR). We estimated the SFR of 22 interacting and 36 single gas-rich dwarf galaxies in the Lynx-Cancer void region using their far-ultraviolet (FUV) images from the GALEX mission. We find an enhancement in SFR by a factor of 3.4$\pm$1.2 for interacting systems compared to single dwarf galaxies in the stellar mass range of 10$^{7}$ - 10$^{8}$ M$\odot$. Our results indicate that dwarf - dwarf galaxy interactions can lead to an enhancement in their SFR. These observations are similar to the predictions based on the simulations of dwarf galaxies at lower redshifts. Future deeper and higher-spatial-resolution UV studies will help us to understand the effect of dwarf galaxy interactions on the spatial distribution of star forming clumps and to identify star formation in tidal tails.

Gamma-ray bursts (GRB) are powerful transient events that emit a large output of gamma rays within a few seconds. Studying these short bursts is vital for cosmological research since they originate from sources observed at large redshifts. To effectively carry out these studies, it is crucial to establish a correlation between the observable features of GRBs while reducing their uncertainty. For these reasons, a comprehensive description of the general GRB light curve (LC) would be crucial for the studies. However, unevenly spaced observations and significant gaps in the LC, which are primarily unavoidable for various reasons, make it difficult to characterize GRBs. Therefore, the general classification of GRB LCs remains challenging. In this study, we present a novel approach to reconstruct gamma-ray burst (GRB) light curves using bidirectional Long Short-Term Memory (BiLSTM). Experimental results show that the BiLSTM approach performs better than traditional methods and produces smoother and more convincing reconstructions for GRBs.

We report results of a spectropolarimetric monitoring of the young Sun-like star V1298~Tau based on data collected with the near-infrared spectropolarimeter SPIRou at the Canada-France-Hawaii Telescope between late 2019 and early 2023. Using Zeeman-Doppler Imaging and the Time-dependent Imaging of Magnetic Stars methods on circularly polarized spectra, we reconstructed the large-scale magnetic topology of the star (and its temporal evolution), found to be mainly poloidal and axisymmetric with an average strength varying from 90 to 170 G over the ~3.5 years of monitoring. The magnetic field features a dipole whose strength evolves from 85 to 245 G, and whose inclination with respect to the stellar rotation axis remains stable until 2023 where we observe a sudden change, suggesting that the field may undergo a polarity reversal, potentially similar to those periodically experienced by the Sun. Our data suggest that the differential rotation shearing the surface of V1298 Tau is about 1.5 times stronger than that of the Sun. When coupling our data with previous photometric results from K2 and TESS and assuming circular orbits for all four planets, we report a $3.9\sigma$ detection of the radial velocity signature of the outermost planet (e), associated with a most probable mass, density and orbital period of $M_e=0.95^{+0.33}_{-0.24} \ \rm M_{\rm jup}$, $\rho_e=1.66^{+0.61}_{-0.48}$ $\rm g\,cm^{-3}$ and $P_e=53.0039\pm0.0001 \ \rm d$, respectively. For the 3 inner planets, we only derive 99\% confidence upper limits on their mass of $0.44\ \rm M_{\rm jup}$, $0.22\ \rm M_{\rm jup}$ and $0.25\ \rm M_{\rm jup}$, for b, c and d, respectively.

Mini-EUSO (Multiwavelength Imaging New Instrument for the Extreme Universe Space Observatory) is a telescope observing the Earth in the ultraviolet band (290-430 nm) from the Russian Zvezda module of the International Space Station since 2019. The telescope is capable of observing UV emissions of cosmic, atmospheric, and terrestrial origin on different time scales. Among the atmospheric phenomena that can be studied, ELVES (Emission of Light and Very low-frequency perturbations due to Electromagnetic pulse Sources) have been photographed by Mini-EUSO with a time resolution of 2.5 us. ELVES are rapidly expanding rings of optical and ultraviolet emissions, 75-95 km in height, resulting from the de-excitation of molecular nitrogen and oxygen in the lower ionosphere following a lightning-associated ElectroMagnetic wave Pulse (EMP). A detailed study of their characteristics, such as radius, speed, and energy, is required for the understanding of these phenomena. In this work, results from the observation of about 30 ELVES with Mini-EUSO will be presented. Using dedicated algorithms, their electro-optical dynamics and morphological characteristics have been thoroughly investigated.

Ultra High Energy Cosmic Rays (UHECRs) offer a unique chance to study the universe at energies inaccessible by man-made accelerators. Observations by ground based observatories have clarified several characteristics of these particles, but their origin, nature, and acceleration mechanisms are still unclear, mostly due to their extremely low flux. Space-based observatories have the potential for an increase in statistics, up to several orders of magnitude, and would be able to cover the whole sky, allowing for a direct comparison of spectra and direction of arrival, but the detector design poses several formidable technical challenges. The JEM-EUSO program has been addressing this problem, trying to open the road of space-based UHECR observations. Several missions have already been completed (on the ground: EUSO-TA; with stratospheric ballons: EUSO-Balloon, EUSO-SPB1 and EUSO-SPB2; in space: TUS\cite{Klimov2017}, MINI-EUSO). Others are under study (K-EUSO) or proposed for the next decade (POEMMA)\cite{Olinto_2021}. In this work we report on the status of the JEM-EUSO program and the technology developed so far.

We report the orbital decay rate of the high mass X-ray binary GX 301$-$2 from an analysis of its long-term X-ray light curves and pulsed flux histories from CGRO/BATSE, RXTE/ASM, Swift/BAT, Fermi/GBM and MAXI by timing the pre-periastron flares over a span of almost 30 years. The time of arrival of the pre-periastron flares exhibits an energy dependence (hard lag) and the orbital period decay was estimated after correcting for it. This method of orbital decay estimation is unaffected by the fluctuations in the spin rate of the X-ray pulsar associated with variations in the mass accretion rate. The resulting $\dot P_\textrm{orb}$ $=-(1.98\pm0.28)\times10^{-6}$ s s$^{-1}$ indicates a rapid evolution timescale of $|P_\textrm{orb}/\dot P_\textrm{orb}|\sim 0.6\times10^{5}$ yr, making it the high mass X-ray binary with the fastest orbital decay. Our estimate of $\dot P_\mathrm{orb}$ is off by a factor of $\sim2$ from the previously reported value of $-(3.7\pm0.5)\times10^{-6}$ s s$^{-1}$ estimated from pulsar TOA analysis. We discuss various possible mechanisms that could drive this rapid orbital decay and also suggest that GX 301$-$2 is a prospective Thorne-\.{Z}ytkow candidate.

In this study, we treat Earth as an exoplanet and investigate our home planet by means of a potential future mid-infrared (MIR) space mission called the Large Interferometer For Exoplanets (LIFE). We combine thermal spectra from an empirical dataset of disk-integrated Earth observations with a noise model for LIFE to create mock observations. We apply a state-of-the-art atmospheric retrieval framework to characterize the planet, assess the potential for detecting the known bioindicators, and investigate the impact of viewing geometry, seasonality, and patchy clouds on the characterization. Key findings include: (1) we are observing a temperate habitable planet with significant abundances of CO2, H2O, O3, and CH4; (2) seasonal variations in the surface and equilibrium temperature, and in the Bond albedo are detectable; (3) the viewing geometry and the spatially and temporally unresolved nature of our observations only have a minor impact on the characterization; (4) Earth's variable H2O profile and patchy cloud coverage lead to biased retrieval results for the atmospheric structure and trace gas abundances; (5) the limited extent of Earth's seasonal variations in biosignature abundances makes the direct detection of its biosphere through atmospheric seasonality unlikely. Our results suggest that LIFE could correctly identify Earth as a planet where life could thrive, with detectable levels of bioindicators, a temperate climate, and surface conditions allowing liquid surface water. Even if atmospheric seasonality is not easily observed, our study demonstrates that next generation, optimized space missions can assess whether nearby temperate terrestrial exoplanets are habitable or even inhabited.

The generation mechanism of compressible fluid turbulence at kiloparsec scales in the Interstellar Medium (ISM) is a long-lasting puzzle. In this work, we explore the nature of large-scale turbulence in the external spiral galaxy NGC~6946. We use the Visibility Moment Estimator (VME) to measure the \HI column density and line of sight turbulent velocity power spectra combining the new observations of A array configuration of Karl G. Jansky Very Large Array (VLA) with the VLA B, C, D array observations from The \HI Nearby Galaxy Survey (THINGS). The estimated power spectra are obeying a power law with a slope of $-0.96\pm0.05$ in column density and $-1.81\pm0.07$ in line of sight velocity in length scales ranging from $6$ kpc to $170$ pc. This points towards a forward energy cascade in the plane of the disc with a driving scale at least as large as $6$ kpc. The values of the power law indices indicate a combination of solenoidal and compressive force responsible in driving the measured turbulence. The presence of strong regular magnetic fields from the magnetic spiral arms in the galaxy is possibly contributing to the solenoidal part, while self-gravity or gravitational instability can mostly be the input for the compressive part of the forcing in the driving mechanism.

This study aims to investigate the ambiguous source and the underlying physical processes of the solar type III radio bursts that occurred on April 3, 2019, through the utilization of multiwavelength observations from the LOFAR radio telescope and the PSP space mission, as well as incorporating results from a PFSS and MHD models. The primary goal is to identify the spatial and temporal characteristics of the radio sources, as well as the plasma conditions along their trajectory. Data preprocessing techniques are applied to combine high- and low-frequency observations from LOFAR and PSP between 2.6 kHz and 80 MHz. We then extract information on the frequency drift and speed of the accelerated electron beams from the dynamic spectra. Additionally, we use LOFAR interferometric observations to image the sources of the radio emission at multiple frequencies and determine their locations and kinematics in the corona. Lastly, we analyze the plasma parameters and magnetic field along the trajectories of the radio sources using PFSS and MHD model results. We present several notable findings related to type III radio bursts. Firstly, through our automated implementation, we were able to effectively identify and characterize 9 type III radio bursts in the LOFAR-PSP combined dynamic spectrum and 16 type III bursts in the LOFAR dynamic spectrum. Secondly, our imaging observations show that the electrons responsible for these bursts originate from the same source and within a short time frame of fewer than 30 minutes. Finally, our analysis provides informative insights into the physical conditions along the path of the electron beams. For instance, we found that the plasma density obtained from the MAS model is significantly lower than the expected theoretical density.

We examine the long-term spectroscopic and photometric variability of EX~Lupi and TW~Hya, studying the presence of stable accretion and the role it plays in the observed variability. Analysing the velocity modulations of the emission lines with STAR-MELT, we obtain information on the structure of the accretion columns and the disk-star connection. The emission line radial velocities reveal that TW Hya, like EX Lupi, has a remarkably stable or slow-varying accretion column footprint, locked to the star for several years. The line-emitting regions are non-polar for both EX Lupi and TW Hya, and species with different energies differ in position. In contrast, the continuum emission as observed in the photometry is very variable and can be modelled by hot spot(s) that change over time in phase, shape, temperature, size, and location with respect to the emission line region. The continuum emission region may not be limited to the stellar surface, especially during episodes of high accretion. The broad line emission observed in EX Lupi during episodes of increased accretion reveals a further structure, which can be fitted by non-axisymmetric disk in Keplerian rotation inwards of the corotation radius. Since the radial velocity modulation due to accretion footprints is so stable, we used it to search for further velocity modulations. While no residual modulation (other than caused by stellar rotation) is found in these objects, a similar analysis could help to identify young planets/companions. Therefore, determining whether stable accretion footprints are common among young stars is a key to detect young planets.

The Square Kilometer array is expected to measure the 21cm signal from the Epoch of Reionization in the coming decade, and its pathfinders may provide a statistical detection even earlier. Current reported upper limits start putting constraints on the astrophysical parameters of the models of the Epoch of Reionization (EoR). In order to interpret such data with 3D radiative hydrodynamics simulations using Bayesian inference, we present the latest developments of the Licorice code. Relying mainly on an implementation of the halo Conditional Mass Function to account for unresolved star formation, they allow for accurate simulations of the EoR at $256^3$ resolution. We use this version of Licorice to produce the first iteration of LoReLi, a public dataset now containing hundreds of 21cm signals computed from radiative hydrodynamics simulations. We train a neural network on LoReLi to provide a fast emulator of the Licorice powerspectra, LorEMU, that has $\sim 5\%$ rms error relative to the simulated signals. LorEMU is used in an MCMC framework to perform Bayesian inference, first on a mock observation composed of a simulated signal and a thermal noise corresponding to 100h observations with the SKA. We then apply our inference pipeline to the latest measurements from the HERA interferometer. We report constraints on the X-ray emissivity, and confirm that cold reionization scenarios are unlikely to accurately represent our Universe.

We have created an up-to-date catalogue of 214 brown dwarfs (BDs) in binaries with $P < 10^4$ d. This allows us to examine the population statistics of the brown dwarf desert. We searched $\textit{Gaia}$ DR3 NSS results for orbital inclinations of BD candidates, deriving 12 new masses. Three remain as desert BDs whereas nine candidates are found to be low-mass stars. We improved the RV solutions for three previously studied BD candidates. A further 19 BD masses with periods less than $\sim$ 1200 d were identified in the DR3 $\texttt{binary_masses}$ database. We confirm a valley in the mass distribution with a minimum around 30-35 M$_{\textrm{jup}}$, and find that periods $<$ 100 d are still under-populated in comparison to longer periods. The updated mass and eccentricity distribution of BDs still shows a marginally statistically significant split into high- and low-mass BD populations. This hints at two different parent distributions, and two potential origins $-$ either akin to planetary formation, or stellar. There are no low eccentricity BDs at periods around 100 d. The mass-metallicity distribution of BDs indicates that core accretion is not the dominant formation mechanism for BDs as they do not follow the same trends that giant exoplanets do with metallicity. We identify a diagonal envelope bounding the Gaia BDs in the mass-period plane due to the detection thresholds of the currently available NSS solutions from 34 months of data.

Significant nonradial, nongravitational accelerations with magnitudes incompatible with radiation-driven effects have been reported in seven small, photometrically inactive near-Earth objects. Two of these objects exhibit large transverse accelerations (i.e., within the orbital plane but orthogonal to the radial direction), and six exhibit significant out-of-plane accelerations. Here, we find that anisotropic outgassing resulting from differential heating on a nucleus with nonzero spin-pole obliquity, averaged over an eccentric orbit, can explain these accelerations for most of the objects. This balanced outgassing model depends on three parameters -- the spin pole orientation (R.A. and Dec.) and an acceleration magnitude. For these "dark comets" (excepting 2003 RM), we obtain parameter values that reproduce the observed nongravitational accelerations. We derive formulae for the component accelerations under certain assumptions for the acceleration scaling over heliocentric distance. Although we lack estimates of these objects' spin axes to confirm our values, this mechanism is nevertheless a plausible explanation for the observed accelerations, and produces accurate perturbations to the heliocentric motions of most of these objects. This model may also be applied to active objects outside of the dark comets group.

Context. The observed scarcity of brown dwarfs in close orbits (within 10 au) around solar-type stars poses significant questions about the origins of these substellar companions. These questions impact our broader understanding of planetary formation processes. However, to resolve these formation mechanisms, accurate observational constraints are essential. Most of the brown dwarfs have been discovered by radial velocity surveys, but this method introduces uncertainties due to its inability to determine the orbital inclination, leaving the true mass-and thus their true nature-unresolved. This highlights the crucial role of astrometric data, helping us distinguish between genuine brown dwarfs and stars. Aims. We aim to refine the mass estimates of massive companions to solar-type stars, mostly discovered through radial velocity measurements and subsequently validated using Gaia DR3 astrometry, to gain a clearer understanding of their true mass and occurrence rates. Methods. We selected a sample of 31 sources with substellar companion candidates validated by Gaia DR3 and with available radial velocities. Using the Gaia DR3 solutions as prior information, we performed an MCMC fit with the available radial velocity measurements to integrate these two sources of data and thus obtain an estimate of their true mass. Results. Combining radial velocity measurements with Gaia DR3 data led to more precise mass estimations, leading us to reclassify several systems initially labeled as brown dwarfs as low-mass stars. Out of the 32 analyzed companions, 13 are determined to be stars, 17 are sub-stellar, and 2 have inconclusive results with the current data. Importantly, using these updated masses, we reevaluated the occurrence rate of brown dwarf companions (13-80 M$_{jup}$ ) on close orbits (<10 au) in the CORALIE sample, determining a tentative occurrence rate of $0.8^{+0.3}_{-0.2}\%$

Galactic winds play a crucial role in the ejection of the interstellar medium (ISM) into intergalactic space. This study presents a model that classifies possible transonic solutions of galactic winds in the gravitational potential of the dark matter halo and stellar component under spherically symmetric and steady assumptions. Our model includes injections of mass and energy resulting from supernovae feedback along a flow line. The mass flux in galactic winds is a critical factor in determining the acceleration process of the flow and revealing the impact of galactic winds on galaxy evolution. We apply the transonic galactic wind model to the observed outflow velocities of star-forming galaxies to estimate the mass flux. Dividing the mass flux by the star formation rate (SFR) yields the mass loading rate (and mass loading factor), which indicates the entrainment effect of the ISM by the hot gas flow. Our results demonstrate that the mass loading rate is inversely correlated with galaxy mass and SFR. In less massive galaxies (stellar mass $\sim 10^{7-8} M_\odot$), the mass loading rate exceeds unity, indicating effective ejection of the ISM into intergalactic space. However, in massive galaxies (stellar mass $\sim 10^{10-11} M_\odot$), the mass loading rate falls below unity, meaning that the mass flux cannot exceed the injected mass by supernovae, thus resulting in the ineffective ejection of the ISM.

Nearby associations are ideal regions to study coeval samples of protoplanetary and debris disks down to late M-type stars. Those aged 5-10,Myrs, where most of the disk should have already dissipated forming planets, are of particular interest. We present the first complete study of both protoplanetary and debris disks in a young region, using the $\eta$ Chamaeleontis association as a test bench to study the cold disk content. We obtained sub-millimeter data for the entire core population down to late M-type stars, plus a few halo members. We performed a continuum sub-millimeter survey with APEX/LABOCA of all the core populations of $\eta$ Cha association. Disk properties have been derived by modeling protoplanetary and debris disks using RADMC 2D and DMS, respectively. We find that protoplanetary disks in $\eta$ Cha typically have holes with radii of the order of 0.01 to 0.03 AU, while ring-like emission from the debris disks is located between 20 au and 650 au from the central star. The parallaxes and Gaia eDR3 photometry, in combination with the PARSEC and COLIBRI isochrones, enable us to confirm an age of $\eta$ Cha between 7 and 9 Myrs. In general, the disk mass seems insufficient to support accretion over a long time, even for the lowest mass accretors, a clear difference compared with other regions and also a sign that the mass budget is further underestimated. We do not find a correlation between the stellar masses, accretion rates, and disk masses, although this could be due to sample issues. We confirm that the presence of inner holes is not enough to stop accretion unless accompanied by dramatic changes to the total disk mass content. Comparing $\eta$ Cha with other regions at different ages, we find that the physical processes responsible for debris disks (e.g., dust growth, dust trapping) efficiently act in less than 5 Myrs.

When stars depart from the main-sequence, various changes occur including the loss of angular momentum owing to changes in the stellar interior and the impact of stellar winds. These processes affect the amount of outer atmospheric heating and emission as revealed by observations in the UV and X-ray spectral regimes. From a theoretical perspective, both magnetic and acoustic energy generation are affected as indicated by detailed theoretical simulations. Here, I will summarize selected observational and theoretical results, including recent work for Beta Hydri (G2~IV), a star constituting a prime example and proxy for the future Sun.

The galaxy total mass inside the effective radius encode important information on the dark matter and galaxy evolution model. Total "central" masses can be inferred via galaxy dynamics or with gravitational lensing, but these methods have limitations. We propose a novel approach, based on Random Forest, to make predictions on the total and dark matter content of galaxies using simple observables from imaging and spectroscopic surveys. We use catalogs of multi-band photometry, sizes, stellar mass, kinematic "measurements" (features) and dark matter (targets) of simulated galaxies, from Illustris-TNG100 hydrodynamical simulation, to train a Mass Estimate machine Learning Algorithm (Mela). We separate the simulated sample in passive early-type galaxies (ETGs), both "normal" and "dwarf", and active late-type galaxies (LTGs) and show that the mass estimator can accurately predict the galaxy dark masses inside the effective radius in all samples. We finally test the mass estimator against the central mass estimates of a series of low redshift (z$\leq$0.1) datasets, including SPIDER, MaNGA/DynPop and SAMI dwarf galaxies, derived with standard dynamical methods based on Jeans equations. Dynamical masses are reproduced within 0.30 dex ($\sim2\sigma$), with a limited fraction of outliers and almost no bias. This is independent of the sophistication of the kinematical data collected (fiber vs. 3D spectroscopy) and the dynamical analysis adopted (radial vs. axisymmetric Jeans equations, virial theorem). This makes Mela a powerful alternative to predict the mass of galaxies of massive stage-IV surveys' datasets.

Neutrino losses play a crucial role in the evolution of massive stars. We study the neutrino luminosity of stars ranging from 20 to 90 M_{\odot} from Zero Age Main Sequence (ZAMS) to Fe Core Collapse (FeCC) with different rotation and metallicity in a neutrino Hertzsprung-Russell diagram. In our simulations, we consider {\omega}/{\omega}crit = 0 and 0.7 to represent non-rotation and high rotation, respectively, and set the metallicities to 0.014, 0.001, and 0.0001. During hydrogen burning stages, neutrino luminosity primarily originates from CNO cycle, and increases with higher stellar mass while decreasing with increasing metallicity. For the high metallicity models (Z = 0.014) during the helium burning stage, the reduction of the hydrogen envelope caused by a larger mass loss rate leads to a gradual decrease in neutrino luminosity. The rapid rotation results in extra mixing inside massive stars, which increases the neutrino luminosity during main sequence (MS), while decreases the neutrino luminosity during helium burning phase. Simultaneously, the rapid rotation also increases CO core mass, which enhances the neutrino luminosity during C and O burning phase. We also investigate the effect of neutrino magnetic moment (NMM) on the massive stars. We find that the energy loss caused by the NMM does not have effects on the evolutionary destiny of massive stars, and it does not significant change the compactness at the time of Fe core collapse.

Approximately one hundred sources of very-high-energy (VHE) gamma rays are known in the Milky Way. A survey of the entire Galactic Plane in the energy range from a few tens of GeV to a few hundred TeV has been proposed as a Key Science Project for the upcoming Cherenkov Telescope Array Observatory (CTAO). This article presents the status of the studies towards the Galactic Plane Survey (GPS). We build and make publicly available a sky model that combines data from observations of known gamma-ray emitters with state-of-the-art physically-driven models of synthetic populations of the main classes of established Galactic VHE sources, as well as of interstellar emission from cosmic-ray interactions in the Milky Way. We also perform an optimisation of the observation strategy. We use the improved sky model and observation strategy to simulate GPS data that are analysed using the methods and software tools under development for real data. We show that the GPS has the potential to increase the number of known Galactic VHE emitters by almost a factor of five. This corresponds to the detection of more than two hundred pulsar wind nebulae and a few tens of supernova remnants at average integral fluxes one order of magnitude lower than in the existing sample above 1 TeV, therefore opening the possibility to perform unprecedented population studies. The GPS also has the potential to provide new VHE detections of binary systems and pulsars, and to identify any bright PeVatrons. Furthermore, the GPS will constitute a pathfinder for deeper follow-up observations of these source classes. Finally, we show that we can extract from GPS data an estimate of the contribution to diffuse emission from unresolved sources, and that there are good prospects of detecting interstellar emission and statistically distinguishing different scenarios. (Abridged)

Class II CH3OH masers are used as a convenient tracer of disc-like structures in high-mass star formation. However, more than half of them show a complex distribution in Very Long Baseline Interferometry (VLBI) maps. The origin of such a complex distribution is still unknown. We conducted VLBI monitoring observations to unveil the origin of a complex class II CH3OH maser in the high-mass star-forming region G59.783+0.065. We observed the CH3OH maser at 6.7 GHz and the H2O maser at 22 GHz to probe detailed circumstellar kinematics and structures by the Japanese VLBI network and the VLBI Exploration of Radio Astrometry. We found similar bipolar distributions in both masers, specifically two clusters located 2000 au apart along the East-West direction. We detected a linear distribution of CH3OH masers in the Western cluster. A position-velocity diagram shows that the Western CH3OH masers trace a rotating disc-wind or infalling component inside an edge-on disc-like structure. In contrast to the simple bipolar expanding motions of the H2O masers, the CH3OH masers exhibited complex motions despite their spatial coincidence. Some of the Eastern CH3OH masers showed bipolar expansions similar to the H2O masers, while others displayed random or even inward motions. Such complex kinematics and their close association with the H2O maser could occur at the boundary between outflow and inflow. We suggest that the complex distribution of class II CH3OH masers, like G59.783+0.065 arises from several distinct circumstellar structures that simultaneously achieve maser excitation.

We studied the timing and spectral properties of the Be/X-ray pulsar Swift J1626.6-5156 using the \emph{NICER} observations of its 2021 outburst. The most important observation is the positive correlation of the centroid energy of the fundamental cyclotron line with the luminosity. This observation agrees with the usual positive correlation of the centroid energy cyclotron line with luminosity in the sub-critical regime. The correlation between the two quantities is verified using two different continuum models. The photon index decreases with an increase in flux. Thus, the spectrum is softer when the flux is low, which may be due to a decrease in the optical depth of the accretion column with a decrease in the flux.

Solar flares sometimes lead to coronal mass ejections that directly affect the Earth's environment. However, a large fraction of flares, including on solar-type stars, are confined flares. What are the differences in physical properties between confined and eruptive flares? For the first time, we quantify thermodynamic and magnetic properties of hundreds of confined and eruptive flares of GOES class C5.0 and above, 480 flares total. We first analyze large flares of GOES class M1.0 and above observed by the Solar Dynamics Observatory (SDO): 216 flares total, including 103 eruptive and 113 confined flares, from 2010 until 2016 April, we then look at the entire dataset above C5.0 of 480 flares. We compare GOES X-ray thermodynamic flare properties, including peak temperature and emission measure, and active-region and flare-ribbon magnetic field properties, including reconnected magnetic flux and peak reconnection rate. We find that for fixed peak X-ray flux, confined and eruptive flares have similar reconnection fluxes; however, for fixed peak X-ray flux confined flares have on average larger peak magnetic reconnection rates, are more compact, and occur in larger active regions than eruptive flares. These findings suggest that confined flares are caused by reconnection between more compact, stronger, lower lying magnetic-fields in larger active regions that reorganizes smaller fraction of these regions' fields. This reconnection proceeds at faster rates and ends earlier, potentially leading to more efficient flare particle acceleration in confined flares.

One of the current challenges of planet formation theory is to explain the enrichment of observed exoplanetary atmospheres. Past studies have focused on scenarios where either pebbles or planetesimals were the heavy element enrichment's drivers, we combine here both approaches to understand whether the composition of a planet can constrain its formation pathway. We study three different formation scenarios: pebble accretion, pebble accretion with planetesimal formation, combined pebble and planetesimal accretion. We use the chemcomp code to perform semi-analytical 1D simulations of protoplanetary discs, including viscous evolution, pebble drift, and simple chemistry to simulate the growth of planets from planetary embryos to gas giants as they migrate through the disc, while tracking their composition. Our simulations confirm that the composition of the planetary atmosphere is dominated by the accretion of gas enriched by inward drifting and evaporating pebbles. Including planetesimal formation hinders the enrichment, because the pebbles locked into planetesimals cannot evaporate and enrich the disc. This results in a big drop of the accreted heavy elements both in the planetesimal formation and accretion case, proving that planetesimal formation needs to be inefficient in order to explain planets with high heavy element content. Accretion of planetesimals enhances the refractory component of the atmosphere, leading to low volatile to refractory ratios, contrary to the pure pebble scenario. Such low volatile to refractory ratios can also be achieved by planets migrating in the inner disc in pure pebble scenario. Distinguishing these two scenarios requires knowledge about the planet's atmospheric C/H and O/H ratios, which are higher for pure pebble accretion. Therefore, a detailed knowledge of the composition of planetary atmospheres could help to constrain the planet's formation pathway.

(abridged) The overwhelming majority of diagnostic tools for galactic activity are focused on active galaxies. Passive or dormant galaxies are often excluded from these diagnostics which usually employ emission line features. In this work, we use infrared and optical colors in order to build an all-inclusive galactic activity diagnostic tool that can discriminate between star-forming, AGN, LINER, composite, and passive galaxies, and which can be used in local and low-redshift galaxies. We explore classification criteria based on infrared colors from the 3 WISE bands supplemented with optical colors from the u, g, and r SDSS bands. From these we aim to find the minimal combination of colors for optimal results. Furthermore, to mitigate biases related to aperture effects, we introduce a new WISE photometric scheme combing different sized apertures. We develop a diagnostic tool using machine learning methods that includes both active and passive galaxies under one unified scheme using 3 colors. We find that the combination of W1-W2, W2-W3, and g-r colors offers good performance while the broad availability of these colors for a large number of galaxies ensures wide applicability on large galaxy samples. The overall accuracy is $\sim$81% while the achieved completeness for each class is $\sim$81% for star-forming, $\sim$56% for AGN, $\sim$68% for LINER, $\sim$65% for composite, and $\sim$85% for passive galaxies. Our diagnostic provides a significant improvement over existing IR diagnostics by including all types of active, as well as passive galaxies, and extending them to the local Universe. The inclusion of the optical colors improves their performance in identifying low-luminosity AGN which are generally confused with star-forming galaxies, and helps to identify cases of starbursts with extreme mid-IR colors which mimic obscured AGN galaxies, a well-known problem for most IR diagnostics.

Using deep archival observations from the Chandra X-ray Observatory, we present an analysis of linear X-ray-emitting features located within the southern portion of the Galactic center chimney, and oriented orthogonal to the Galactic plane, centered at coordinates l = 0.08 deg, b = -1.42 deg. The surface brightness and hardness ratio patterns are suggestive of a cylindrical morphology which may have been produced by a plasma outflow channel extending from the Galactic center. Our fits of the feature's spectra favor a complex two-component model consisting of thermal and recombining plasma components, possibly a sign of shock compression or heating of the interstellar medium by outflowing material. Assuming a recombining plasma scenario, we further estimate the cooling timescale of this plasma to be on the order of a few hundred to thousands of years, leading us to speculate that a sequence of accretion events onto the Galactic Black Hole may be a plausible quasi-continuous energy source to sustain the observed morphology.

We present Rainbow, a physically motivated framework which enables simultaneous multi-band light curve fitting. It allows the user to construct a 2-dimensional continuous surface across wavelength and time, even in situations where the number of observations in each filter is significantly limited. Assuming the electromagnetic radiation emission from the transient can be approximated by a black-body, we combined an expected temperature evolution and a parametric function describing its bolometric light curve. These three ingredients allow the information available in one passband to guide the reconstruction in the others, thus enabling a proper use of multi-survey data. We demonstrate the effectiveness of our method by applying it to simulated data from the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) as well as real data from the Young Supernova Experiment (YSE DR1).We evaluate the quality of the estimated light curves according to three different tests: goodness of fit, time of peak prediction and ability to transfer information to machine learning (ML) based classifiers. Results confirm that Rainbow leads to equivalent (SNII) or up to 75% better (SN Ibc) goodness of fit when compared to the Monochromatic approach. Similarly, accuracy when using Rainbow best-fit values as a parameter space in multi-class ML classification improves for all classes in our sample. An efficient implementation of Rainbow has been publicly released as part of the light curve package at https://github.com/light-curve/light-curve. Our approach enables straight forward light curve estimation for objects with observations in multiple filters and from multiple experiments. It is particularly well suited for situations where light curve sampling is sparse.

We study the angular distribution of the cosmic microwave background (CMB) temperature fluctuations to probe the statistical isotropy of the universe by using precise full-sky CMB data with a model-independent approach. We investigate the temperature-temperature angular correlations in the four Planck foreground-cleaned CMB maps, recently released. This study performs a directional analysis on the CMB sphere looking for directions where temperature-temperature angular correlations are extremes. Our analyses confirm a preferred axis in the CMB sphere, pointing in the direction $(l,b) \simeq (250^{\circ}, 130^{\circ})$, at $98\%$-$99\%$ confidence level, direction where the CMB angular correlations show the maximum excess over the antipodal direction. This preferred direction is unexpected in the $\Lambda$CDM cosmological model, and represents a significant deviation compared to results obtained applying the same procedure to simulated statistically isotropic CMB maps. This result confirms the North-South asymmetry in the most recent Planck data, a phenomenon that is one of the previously reported CMB anomalies. We perform a robust detection of the North-South asymmetry in the temperature-temperature angular correlations -- with slightly different statistical significance -- in the four Planck foreground-cleaned CMB maps. Moreover, we perform consistency tests by adding foregrounds and noise, both Planck data products, to the CMB map in study, and also investigate and discard possible bias in our methodology. After these detailed analyses we conclude that the North-South asymmetry phenomenon is present, with high statistical significance, in the Planck CMB maps studied, result that confirms previous reports in the literature in the last 20 years.

We explore the potential of the Gaussian Mixture Model (GMM), an unsupervised machine learning method, to identify coherent physical structures in the ISM. The implementation we present can be used on any kind of spatially and spectrally resolved data set. We provide a step-by-step guide to use these models on different sources and data sets. Following the guide, we run the models on NGC 1977, RCW 120 and RCW 49 using the [CII] 158 $\mu$m mapping observations from the SOFIA telescope. We find that the models identified 6, 4 and 5 velocity coherent physical structures in NGC 1977, RCW 120 and RCW 49, respectively, which are validated by analysing the observed spectra towards these structures and by comparison to earlier findings. In this work we demonstrate that GMM is a powerful tool that can better automate the process of spatial and spectral analysis to interpret mapping observations.

The luminosity of X-ray pulsars is their key parameter determining the geometry and physical conditions of the accretion flow both on the spatial scales of a binary system and on much smaller scales of emitting regions located close to the stellar surface. Traditionally, the luminosity of X-ray pulsars is estimated out of the X-ray energy flux averaged over the pulsed period and the estimated distance to the source. Due to the anisotropy of X-ray emission, the luminosity estimated on the base of the observed pulse profile can differ from the actual one. Super-critical X-ray pulsars with accretion columns are of particular interest because the X-ray flux from columns is a matter of strong gravitational lensing by a neutron star. Using toy model of an accretion column, we simulate beam patterns in super-critical X-ray pulsars, construct theoretical pulse profiles for different geometries and mutual orientations of pulsars and distant observers and show that despite strong light bending, the typical deviation of the apparent luminosity from the actual one is $\sim 20\%$ only, and in $\sim 90\%$ of cases, the apparent luminosity $0.8 L\lesssim L_{\rm app}\lesssim 1.25 L$. However, the shape of the pulse profiles is strongly affected by the geometry of the emitting region. We show that the appearance and growth of accretion columns tend to be accompanied by an increase of observed pulsed fraction, which is in agreement with the recent observations of bright X-ray transients.

One of the primary outstanding questions in extragalactic astronomy is the formation and early evolution of the supermassive black holes that are seen in nearly every massive galaxy. Low metallicity dwarf galaxies may offer the most representative local analogs to pristine early galaxies, making them a vital tool in probing black hole seed models through the study of the intermediate mass black holes (IMBHs) possibly hosted therein, though these dwarf galaxies, and the IMBHs they may host, are typically not as well-studied in this context as their higher metallicity and higher mass counterparts. In this paper, we explore the X-ray properties of a sample of 37 low metallicity dwarf galaxies using archival XMM observations, and we compare the properties of this population against a representative sample of higher metallicity counterparts. We report the detection of ten sources with 0.3-10 keV luminosity in excess of $10^{40}$~erg~s$^{-1}$ within the low metallicity sample, which we highlight for follow-up as potential intermediate mass black hole candidates. Finally, we discuss the differing multi-wavelength scaling relations (e.g., $L_X - L_{W2}$, $L_X-SFR$) between the two galaxy populations, as well as the sample's $L_X$ as a function of metallicity.

EUSO-SPB2 (Extreme Universe Space Observatory on a Super Pressure Balloon II) is a precursor mission for a future space observatory for multi-messenger astrophysics, planned to be launched in Spring 2023 with a flight duration target of 100 days. The Fluorescence Telescope (FT) hosted on board is designed to detect Ultra High Energy Cosmic Rays via the UV fluorescence emission of the Extensive Air Showers in the atmosphere. The Data Processor (DP) of the FT is the component of the electronics system that performs data management and instrument control for the telescope. The DP controls front-end electronics, tags events with arrival time and payload position through a GPS system, provides signals for time synchronization of the event and measures the live and dead time of the telescope. Furthermore, it manages mass memory for data storage, performs housekeeping monitoring, and controls the power-on and power-off sequences. Finally, the data processor combines the data from the PDMs and onboard differential GPS and prioritizes data for download. The long duration of the flight poses strict requirements on electronics and data handling. The operations at high altitude in an unpressurized environment represent a technological challenge for heat dissipation. This contribution will provide an overview of the innovative elements developed and the results of the integration and field test campaigns. We will also present some preliminary analysis of the performance during the flight.

Using the Very Long Baseline Array, we observed the active galactic nucleus (AGN) in NGC 3079 over a span of six months to test for variability in the two main parsec-scale radio components, $A$ and $B$, which lie on either side of the AGN. We found evidence for positional differences in the positions of $A$ and $B$ over the six months consistent with the apparent motion of these components extrapolated from older archival data, finding that their projected rate of separation, $(0.040\pm0.003)$ c, has remained constant since $\sim2004$ when a slowdown concurrent with a dramatic brightening of source $A$ occurred. This behavior is consistent with an interaction of source $A$ with the interstellar medium (ISM), as has previously been suggested in the literature. We calculated the amount of mechanical feedback on the ISM for both the scenario in which $A$ is an expulsion of material from the central engine and the scenario in which $A$ is a shock front produced by a relativistic jet, the latter of which is favored by several lines of evidence we discuss. We find that the cumulative mechanical feedback on the ISM is between $2 \times 10^{44}$ erg to $1 \times 10^{48}$ erg for the expulsion scenario or between $3\times 10^{50}$ erg to $1 \times 10^{52}$ erg for the jet scenario. Integrated over the volume-complete FRAMEx sample, our results imply that jet-mode mechanical feedback plays a negligible role in the energetics of AGNs in the local universe.

Different star-formation models at Cosmic Dawn produce detectable signatures in the observables of upcoming 21-cm experiments. In this work, we consider the physical scenario of feedback-free starbursts (FFB), according to which the star-formation efficiency (SFE) is enhanced in sufficiently massive halos at early enough times, thus explaining the indication from the James Webb Space Telescope for an excess of bright galaxies at $z \geq 10$. We model the contribution of FFBs to popII SFE and compute the impact these have on the 21-cm global signal and power spectrum. We show that FFBs affect the evolution of the brightness temperature and the 21-cm power spectrum, but they only have a limited effect on the neutral hydrogen fraction. We investigate how the observables are affected by changes in the underlying star formation model and by contribution from popIII stars. Finally, we forecast the capability of next-generation Hydrogen Epoch of Reionization Array (HERA) to detect the existence of FFB galaxies via power spectrum measurements. Our results show the possibility of a significant detection, provided that popII stars are the main drivers of lowering the spin temperature. Efficient popIII star formation will make the detection more challenging.

Multi-scale modeling of expanding plasmas is crucial for understanding the dynamics and evolution of various astrophysical plasma systems such as the solar and stellar winds. In this context, the Expanding Box Model (EBM) provides a valuable framework to mimic plasma expansion in a non-inertial reference frame, co-moving with the expansion but in a box with a fixed volume, which is especially useful for numerical simulations. Here, fundamentally based on the Vlasov equation for magnetized plasmas and the EBM formalism for coordinates transformations, for the first time we develop a first-principles description of radially expanding plasmas in the EB frame. From this approach, we aim to fill the gap between simulations and theory at microscopic scales to model plasma expansion at the kinetic level. Our results show that expansion introduces non-trivial changes in the Vlasov equation (in the EB frame), especially affecting its conservative form through non-inertial forces purely related to the expansion. In order to test the consistency of the equations, we also provide integral moments of the modified Vlasov equation, obtaining the related expanding moments (i.e., continuity, momentum, and energy equations). Comparing our results with the literature, we obtain the same fluids equations (ideal-MHD), but starting from a first principles approach. We also obtained the tensorial form of the energy/pressure equation in the EB frame. These results show the consistency between the kinetic and MHD descriptions. Thus, the expanding Vlasov kinetic theory provides a novel framework to explore plasma physics at both micro and macroscopic scales in complex astrophysical scenarios.

We perform a state-of-the-art global study of the cosmological thermal histories of a simple Yukawa model, and find higher perturbative orders to be important for determining both the presence and strength of strong first-order phase transitions. Using high-temperature effective field theory, we calculate the free energy density of the model up to $\mathcal{O}(y^5T^4)$, where $y$ is the Yukawa coupling and $T$ is the temperature. The locations of phase transitions are found using the results of lattice Monte-Carlo simulations, and the strength of first-order transitions are evaluated within perturbation theory, to 3-loop order. This is the first global study of any model at this order. Compared to a vanilla 1-loop analysis, accurate to $\mathcal{O}(y^2 T^4)$, reaching such accuracy enables on average a five-fold reduction in the relative error in the predicted critical temperature $T_\text{c}$, and an additional $\sim50\%$ strong first-order transitions with latent heat $L/T_\text{c}^4 > 0.1$ to be identified in our scan.

Models for the observational appearance of astrophysical black holes rely critically on accurate general-relativistic ray tracing and radiation transport to compute the intensity measured by a distant observer. In this paper, we illustrate how the choice of coordinates and initial conditions affect this process. In particular, we show that propagating rays from the camera to the source leads to different solutions if the spatial part of the momentum of the photon points towards the horizon or away from it. In doing this, we also show that coordinates that are well suited for numerical General-Relativistic MagnetoHydroDynamic (GRMHD) simulations are typically not optimal for generic ray tracing. We discuss the implications for black-hole images and show that radiation transport in optimal and non-optimal spacetime coordinates lead to the same images up to numerical errors and algorithmic choices.

Accurate determination of Cassini's spacecraft potential in Titan's ionosphere is important for interpreting measurements by its low energy plasma instruments. Estimates of the floating potential varied significantly, however, between the various different plasma instruments. In this study we utilize 3-D particle-in-cell simulations to understand the key features of Cassini's plasma interaction in Titan's ionosphere. The spacecraft is observed to charge to negative potentials for all scenarios considered, and close agreement is found between the current onto the simulated Langmuir Probe and that observed in Titan's ionosphere. These simulations are therefore shown to provide a viable technique for modeling spacecraft interacting with Titan's dusty ionosphere.

The next generation of ground-based gravitational-wave detectors will look much deeper into the Universe and have unprecedented sensitivities and low-frequency capabilities. Especially alluring is the possibility of detecting an early-Universe cosmological stochastic background that could provide important insights into the beginnings of our Universe and fundamental physics at extremely high energies. However, even if next-generation detectors are sensitive to cosmological stochastic backgrounds, they will be masked by more dominant astrophysical backgrounds, namely the residual background from the imperfect subtraction of resolvable compact binary coalescences (CBCs) as well as the CBC background from individually unresolvable CBCs. Using our latest knowledge of masses, rates, and delay time distributions, we present a data-driven estimate of the unresolvable CBC background that will be seen by next-generation detectors. Accounting for statistical and systematic errors, this estimate quantifies an important piece in the CBC noise budget for next-generation detectors and can help inform detector design and subtraction algorithms. We compare our results with predictions for backgrounds from several cosmological sources in the literature, finding that the unresolvable background will likely be a significant impediment for many models. This motivates the need for simultaneous inference methods or other statistical techniques to detect early-Universe cosmological backgrounds.

This research focuses on scalar field cosmologies with a generalized harmonic potential. Our attention is centred on the anisotropic LRS Bianchi I and III metrics, Bianchi V metrics, and their isotropic limits. We provide a comprehensive overview of the first two metrics classes and offer new findings for Bianchi V metrics. We show that the Hubble parameter is a time-dependent perturbation parameter that controls the magnitude of the error between full-system and time-averaged solutions as it decreases, such that those complete and time-averaged systems have the same asymptotic behaviour. Therefore, oscillations entering the system can be controlled and smoothed out, which simplifies the problem at hand.

Assuming that neutrinos are spacelike (tachyonic) fermions, we calculate width for the kinematically allowed, lepton number conserving, three-body decay $\nu_{\alpha}\rightarrow \nu_{\alpha} \; \nu_{\beta} \bar{\nu}_{\beta}$ in the Standard Model. Decays of tachyonic neutrinos over cosmological distances can lead to a reduction of the neutrino flux in the high-energy end of the spectrum. We estimate upper limits on the spacelike neutrino mass based on the PeV-energy cosmological neutrino events observed in the IceCube experiment. These limits are close to those deduced from the measurements of $m_{\nu}^2$ in the tritium-decay experiment KATRIN.

The post-Newtonian orbital effects induced by the mass quadrupole and spin octupole moments of an isolated, oblate spheroid of constant density that is rigidly and uniformly rotating on the motion of a test particle are analytically worked out for an arbitrary orbital configuration and without any preferred orientation of the body's spin axis. The resulting expressions are specialized to the cases of a) equatorial and b) polar orbits. The opportunity offered by a hypothetical new spacecraft moving around Jupiter along a Juno-like highly elliptical, polar orbit to measure them is preliminarily studied. Although more difficult to be practically implemented, also the case of a less elliptical orbit is considered since it yields much larger figures for the relativistic effects of interest. The possibility of using the S-stars orbiting the supermassive black hole in Sgr A$^\ast$ at the Galactic Center as probes to potentially constrain some parameters of the predicted extended mass distribution surrounding the hole by means of the aforementioned orbital effects is briefly examined.

When a test particle moves about an oblate spheroid, it is acted upon, among other things, by two standard perturbing accelerations. One, of Newtonian origin, is due to the quadrupole mass moment $J_2$ of the orbited body. The other one, of the order of $\mathcal{O}\left(1/c^2\right)$, is caused by the static, post-Newtonian field arising solely from the mass of the central object. Both of them concur to induce \textrm{indirect}, \textrm{mixed} orbital effects of the order of $\mathcal{O}\left(J_2/c^2\right)$. They are of the same order of magnitude of the \textrm{direct} ones induced by the post-Newtonian acceleration arising in presence of an oblate source, not treated here. We calculate these less known features of motion in their full generality in terms of the osculating Keplerian orbital elements. Subtleties pertaining the correct calculation of their mixed net \textrm{precessions} per orbit to the full order of $\mathcal{O}\left(J_2/c^2\right)$ are elucidated. The obtained results hold for arbitrary orbital geometries and for any orientation of the body's spin axis $\boldsymbol{\hat{k}}$ in space. The method presented is completely general, and can be extended to any pair of post-Keplerian accelerations entering the equations of motion of the satellite, irrespectively of their physical nature.