Papers for Tuesday, Feb 20 2024

The Sun and solar-type stars exhibit irregular cyclic variations in their magnetic activity over long time scales. To understand this irregularity, we employed the flux transport dynamo models to investigate the behavior of one solar mass star at various rotation rates. To achieve this, we have utilized a mean-field hydrodynamic model to specify differential rotation and meridional circulation, and we have incorporated stochastic fluctuations in the Babcock-Leighton source of the poloidal field to capture inherent fluctuations in the stellar convection. Our simulations successfully demonstrated consistency with the observational data, revealing that rapidly rotating stars exhibit highly irregular cycles with strong magnetic fields and no Maunder-like grand minima. On the other hand, slow rotators produce smoother cycles with weaker magnetic fields, long-term amplitude modulation, and occasional extended grand minima. We observed that the frequency and duration of grand minima increase with the decreasing rotation rate. These results can be understood as the tendency of a less supercritical dynamo in slow rotators to be more prone to produce extended grand minima. We further explore the possible existence of the dynamo in the subcritical regime in a Babcock-Leighton-type framework and in the presence of a small-scale dynamo.

These notes summarise the contents of the lectures I delivered at the International School of Physics "Enrico Fermi" on "Foundations of Cosmic Ray Astrophysics". The lectures were dealing with the physics of Pulsars and Pulsar Wind Nebulae (PWNe) in the Cosmic Ray (CR) perspective. It has become now clear that the processes taking place in the environment of fast rotating, highly magnetized neutron stars, often detected as pulsars, play a crucial role in the formation of the CR spectrum detected at the Earth. These lectures discuss the main aspects of this connection. Pulsars are likely contributors of the CR lepton flux at the Earth thanks to their nature of electron-positron factories. Pulsars and their nebulae are the best potential leptonic PeVatron in the Galaxy, and the Crab Nebula, the prototype of the Pulsar Wind Nebula class is the only established PeVatron in the Galaxy. Pulsars are however also potential sources of high energy hadrons, up to the energies relevant for UHECRs. Pulsars and their nebulae are the best potential leptonic PeVatrons in the Galaxy, and the Crab Nebula, the prototype of the Pulsar Wind Nebula class, is the only established PeVatron in the Galaxy. Finally, regions of suppressed particle diffusion have been observed around evolved pulsars, the so-called TeV halos, which could have an impact on galactic CR transport. These lectures discuss the physics of pulsars and PWNe, summarising what we know about these systems and what pieces of information are still missing to fully assess their role in all the above mentioned Cosmic Ray connected aspects.

Magnetic bright points on the solar photosphere mark the footpoints of kilogauss magnetic flux tubes extending toward the corona. Convective buffeting of these tubes is believed to excite magnetohydrodynamic waves, which can propagate to the corona and there deposit heat. Measuring wave excitation via bright-point motion can thus constrain coronal and heliospheric models, and this has been done extensively with centroid tracking, which can estimate kink-mode wave excitation. DKIST is the first telescope to provide well-resolved observations of bright points, allowing shape and size measurements to probe the excitation of other wave modes that have been difficult, if not impossible, to study to date. In this work, we demonstrate a method of automatic bright-point tracking that robustly identifies the shapes of bright points, and we develop a technique for interpreting measured bright-point shape changes as the driving of a range of thin-tube wave modes. We demonstrate these techniques on a MURaM simulation of DKIST-like resolution. These initial results suggest that modes other than the long-studied kink mode could increase the total available energy budget for wave-heating by 50%. Pending observational verification as well as modeling of the propagation and dissipation of these additional wave modes, this could represent a significant increase in the potency of wave-turbulence heating models.

The Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory is poised to observe thousands of quasars using the Deep Drilling Fields (DDF) across six broadband filters over a decade. Understanding quasar accretion disc (AD) time delays is pivotal for probing the physics of these distant objects. Pozo Nu\~nez et al. (2023) has recently demonstrated the feasibility of recovering AD time delays with accuracies ranging from 5\% to 20\%, depending on the quasar's redshift and time sampling intervals. Here we reassess the potential for measuring AD time delays under the current DDF observing cadence, which is placeholder until a final cadence is decided. We find that contrary to prior expectations, achieving reliable AD time delay measurements for quasars is significantly more challenging, if not unfeasible, due to the limitations imposed by the current observational strategies.

In 1953 I heard of an experiment in 1925 by Bengt Str\"omgren where he observed transit times with the meridian circle at the Copenhagen University Observatory measuring the current in a photocell behind slits when a star was crossing. In 1954 just 22 years old I was given the task as a student to make first test observations with a new meridian circle of the observatory. I became fascinated by the instrument and by the importance of astrometry for astronomy. Work at four meridian circles, two in Denmark, one in Hamburg, one in Lund, and Pierre Lacroute's vision of space astrometry in France had by 1973 created the foundation for development of the Hipparcos satellite, and Gaia followed. In 2013 I proposed a successor satellite which has gained momentum especially thanks to the efforts of David Hobbs and it has a good chance to be launched by ESA about 2045. But 70 years ago, optical astrometry was considered a dying branch of astronomy, unattractive compared with astrophysics. The following growth built on the still active interest in astrometry in Europe in those years and it was supported by ESA, the European Space Agency. This review is only about astrometry where I was personally involved.

Energy and momentum feedback from stars is a key element of models for galaxy formation and interstellar medium dynamics, but resolving the relevant length scales to directly include this feedback remain out of reach of current-generation simulations. We aim to constrain the energy feedback by winds, photoionisation and supernovae (SNe) from massive stars. We measure the thermal and kinetic energy imparted to the interstellar medium on various length scales, calculated from high-resolution 1D radiation-hydrodynamics simulations. Our grid of simulations covers a broad range of densities, metallicities, and state-of-the-art evolutionary models of single and binary stars. We find that a single star or binary system can carve a cavity of tens-of-pc size into the surrounding medium. During the pre-SN phase, post-main-sequence stellar winds and photoionisation dominate. While SN explosions dominate the total energy budget, the pre-SN feedback is of great importance by reducing the circumstellar gas density and delaying the onset of radiative losses in the SN remnant. Contrary to expectations, the metallicity dependence of the stellar wind has little effect on the cumulative energy imparted by feedback to the ISM; the only requirement is the existence of a sufficient level of pre-SN radiative and mechanical feedback. The ambient medium density determines how much and when feedback energy reaches to distance $\gtrsim 10-20$ pc and affects the division between kinetic and thermal feedback. Our results can be used as a sub-grid model for feedback in large-scale simulations of galaxies. The results reinforce that the uncertain mapping of stellar evolution sequences to SN explosion energy is very important to determining the overall feedback energy from a stellar population.

Using a set of general-relativistic magnetohydrodynamics simulations that include proper neutrino transfer, we assess for the first time the role played by the fallback accretion onto the remnant from a binary neutron-star merger over a timescale of hundreds of seconds. In particular, we find that, independently of the equation of state, the properties of the binary, and the fate of the remnant, the fallback material reaches a total mass of $\gtrsim 10^{-3}\,M_\odot$, i.e. about $50\%$ of the unbound matter, and that the fallback accretion rate follows a power-law in time with slope $\sim t^{-5/3}$. Interestingly, the timescale of the fallback and the corresponding accretion luminosity are in good agreement with the so-called ``extended emission'' observed in short gamma-ray bursts (GRBs). Using a simple electromagnetic emission model based on the self-consistent thermodynamical state of the fallback material heated by r-process nucleosynthesis, we show that this fallback material can shine in the gamma- and X-rays with luminosities $\gtrsim \,10^{48}\,{\rm erg/s}$ for hundreds of seconds, thus making it a good and natural candidate to explain the extended emission in short GRBs. In addition, our model for the emission by the fallback material reproduces well and rather naturally some of the phenomenological traits of the extended emission, such as its softer spectra with respect to the prompt emission and the presence of exponential cutoffs in time. Our results clearly highlight that fallback flows onto merger remnants cannot be neglected and the corresponding emission represents a very promising and largely unexplored avenue to explain the complex phenomenology of GRBs.

Theoretical and observational approaches to settle the important questions of the progenitor systems and the explosion mechanism of normal Type Ia supernovae thus far have failed. With its unique capability to obtain continuous spectra through the near- and mid-Infrared, JWST now offers completely new insights into Type Ia supernovae. In particular, observing them in the nebular phase allows us to directly see the central ejecta and thereby constrain the explosion mechanism. We aim to understand and quantify differences in the structure and composition of the central ejecta of various Type Ia supernova explosion models. We cover the currently most popular explosion scenarios using self-consistent multi-dimensional explosion simulations of delayed-detonation and pulsationally-assisted gravitationally-confined delayed detonation Chankdrasekhar-mass models, and double detonation sub-Chandrasekhar-mass and violent merger models. We focus on the distribution of radioactive and stable nickel in the final ejecta in density and velocity space. Of all models we study, we find that only in the double detonation scenario with a surviving companion the inner ejecta are anywhere close to spherical symmetry. It is thus crucial to simulate in particular Chandrasekhar-mass explosion models in 3D, because the lack of buoyancy in 1D models leads to an unphysical inverted stratification of radioactive and stable nickel in 1D models. The ejecta of Chandrasekhar-mass and sub-Chandrasekhar-mass merger models with an exploding secondary white dwarf are inherently 3D. In all of them nickel is distributed over a wide range of densities even at fixed velocity. Therefore, it is crucial to compute synthetic observables in the nebular phase, when the ejecta are optically thin, in 3D from the full 3D ejecta, because spherical averaging instead leads to unphysical ejecta properties. (abridged)

Electromagnetic fundamental and harmonic emission is ubiquitously observed throughout the heliosphere, and in particular it is commonly associated with the occurrence of Type II and III solar radio bursts. Classical analytic calculations for the plasma-emission process, though useful, are limited to idealized situations; a conclusive numerical verification of this theory is still lacking, with earlier studies often providing contradicting results on e.g. the precise parameter space in which fundamental and harmonic emission can be produced. To accurately capture the chain of mechanisms underlying plasma emission - from precursor plasma processes to the generation of electromagnetic waves over long times - we perform large-scale, first-principles simulations of beam-plasma instabilities. By employing a very large number of computational particles we achieve very low numerical noise, and explore (with an array of simulations) a wide parameter space determined by the beam-plasma density ratio and the ion-to-electron temperature ratio. In particular, we observe direct evidence of both fundamental and harmonic plasma emission when the beam-to-background density ratio $\le$0.005 (with beam-to-background energy ratio ~0.5), tightly constraining this threshold. We observe that, asymptotically, in this regime ~0.1% of the initial beam energy is converted into harmonic emission, and ~0.001% into fundamental emission. In contrast with previous studies, we also find that this emission is independent of the ion-to-electron temperature ratio. In addition, we report the direct detection of third-harmonic emission in all of our simulations, at power levels compatible with observations. Our findings have important consequences for understanding the viable conditions leading to plasma emission in space systems, and for the interpretation of observed electromagnetic signals throughout the heliosphere.

We carry out an in-depth analysis of the prompt-collapse behaviour of binary neutron star (BNS) mergers. To this end, we perform more than $80$ general relativistic BNS merger simulations using a family of realistic Equations of State (EOS) with different stiffness, which feature a first order deconfinement phase transition between hadronic and quark matter. From these simulations we infer the critical binary mass $M_{\rm crit}$ that separates the prompt from the non-prompt collapse regime. We show that the critical mass increases with the stiffness of the EOS and obeys a tight quasi-universal relation, $M_{\rm crit}/M_{\rm TOV}\approx 1.41\pm 0.06$, which links it to the maximum mass $M_{\rm TOV}$ of static neutron stars, and therefore provides a straightforward estimate for the total binary mass beyond which prompt collapse becomes inevitable. In addition, we introduce a novel gauge independent definition for a one-parameter family of threshold masses in terms of curvature invariants of the Riemann tensor which characterizes the development toward a more rapid collapse with increasing binary mass. Using these diagnostics, we find that the amount of matter remaining outside the black hole sharply drops in supercritical mass mergers compared to subcritical ones and is further reduced in mergers where the black hole collapse is induced by the formation of a quark matter core. This implies that $M_{\rm crit}$, particularly for merger remnants featuring quark matter cores, imposes a strict upper limit on the emission of any detectable electromagnetic counterpart in BNS mergers.

Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846,000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We identified 13 candidate white dwarfs with a counterpart in VLASS within 2". Five of those were found not to be white dwarfs in follow-up or archival spectroscopy, whereas seven others were found to be chance alignments with a background source in higher-resolution optical or radio images. The remaining source, WDJ204259.71+152108.06, is found to be a white dwarf and M-dwarf binary with an orbital period of 4.1 days and long-term stochastic optical variability, as well as luminous radio and X-ray emission. For this binary, we find no direct evidence of a background contaminant, and a chance alignment probability of only ~2 per cent. However, other evidence points to the possibility of an unfortunate chance alignment with a background radio and X-ray emitting quasar, including an unusually poor Gaia DR3 astrometric solution for this source. With at most one possible radio emitting white dwarf found, we conclude that strong (> 1-3 mJy) radio emission from white dwarfs in the 3 GHz band is virtually nonexistent outside of interacting binaries.

Hyper-luminous X-ray sources (HLXs) are extragalactic off-nuclear X-ray sources with luminosities exceeding the theoretical limit for accretion onto stellar-mass compact objects. Many HLXs may represent intermediate-mass black holes (IMBHs) deposited in galaxy halos through mergers, and properties of the stellar cores surrounding HLXs provide powerful constraints on this scenario. Therefore, we have systematically built the largest sample of HLX candidates with archival Hubble Space Telescope (HST) imaging (24) for the first uniform population study of HLX stellar cores down to low masses. Based on their host galaxy redshifts, at least 21 (88%) have stellar core masses >=10^7 Msun and hence are consistent with accretion onto massive black holes from external galaxies. In 50% of the sample, the HST imaging reveals features connecting the HLXs with their host galaxies, strongly suggesting against the background/foreground contaminant possibility in these cases. Assuming a mass scaling relation for active galactic nuclei and accounting for an estimated contamination fraction of 29%, up to ~60% of our sample may be associated with IMBHs. Similar to previously known HLXs, the X-ray luminosities are systematically elevated relative to their stellar core masses, possibly from merger-driven accretion rate enhancements. The least massive stellar cores are preferentially found at larger nuclear offsets and are more likely to remain wandering in their host galaxy halos. The HLX galaxy occupation fraction is ~10^-2 and has a strong inverse mass dependence. Up to three of the HLX candidates (12%) are potentially consistent with formation within globular clusters or with exceptionally luminous X-ray binaries.

Long-period radio transients have unusual properties that challenge their interpretation as pulsars or magnetars. We examine whether they might instead be powered by primordial black holes (PBHs) making repeated passages through a host star, thereby providing a signature of elusive dark-matter candidates. We demonstrate that constraints derived from the transients' period and period derivative alone already rule out this scenario for most potential host stars. While white dwarfs may satisfy these constraints, they are unlikely to capture PBHs in the required mass range.

Aims. We present a new sub-grid model, HYACINTH -- HYdrogen And Carbon chemistry in the INTerstellar medium in Hydro simulations, for computing the non-equilibrium abundances of ${\rm H_2}$ and its carbon-based tracers, namely ${\rm CO}$, ${\rm C}$, and ${\rm C^+}$, in cosmological simulations of galaxy formation. Methods. The model accounts for the unresolved density structure in simulations using a variable probability distribution function of sub-grid densities and a temperature-density relation. Included is a simplified chemical network tailored for hydrogen and carbon chemistry within molecular clouds and easily integrated into large-scale simulations with minimal computational overhead. As an example, we apply HYACINTH to a simulated galaxy at redshift $z\sim2.5$ in post-processing and compare the resulting abundances with observations. Results. The chemical predictions from HYACINTH show a good agreement with high-resolution molecular-cloud simulations. We reproduce the $\rm H\,I - {\rm H_2}$ transition in the $f_{\rm H_2} - N_{\rm H}$ plane for both Milky-Way and LMC-like conditions. We also match the $N_{\rm CO} - N_{\rm H_2}$ values inferred from absorption measurements towards Milky-Way molecular clouds. Column density maps reveal that ${\rm CO}$ is concentrated in the peaks of the $\rm H_2$ distribution, while atomic carbon more broadly traces the bulk of ${\rm H_2}$ in our post-processed galaxy. Based on surface density profiles of stars and different gas species in the post-processed galaxy, we find that ${\rm C^+}$ extends farther than all other components and maintains a substantially high surface density out to $\sim 10 \, \rm kpc$. This is similar to the $[\rm C\,II]$ halos found in some recent observations at high redshifts.

We show that the abundance of primordial black holes, if formed through the collapse of large fluctuations generated during inflation and unless the power spectrum of the curvature perturbation is very peaked, is always dominated by the broadest profile of the compaction function, even though statistically it is not the most frequent. The corresponding threshold is therefore 2/5. This result exacerbates the tension when combining the primordial black hole abundance with the signal seen by pulsar timing arrays and originated from gravitational waves induced by the same large primordial perturbations.

The SETI Ellipsoid is a strategy for technosignature candidate selection which assumes that extraterrestrial civilizations who have observed a galactic-scale event -- such as supernova 1987A -- may use it as a Schelling point to broadcast synchronized signals indicating their presence. Continuous wide-field surveys of the sky offer a powerful new opportunity to look for these signals, compensating for the uncertainty in their estimated time of arrival. We explore sources in the TESS continuous viewing zone, which corresponds to 5% of all TESS data, observed during the first three years of the mission. Using improved 3D locations for stars from Gaia Early Data Release 3, we identified 32 SN 1987A SETI Ellipsoid targets in the TESS continuous viewing zone with uncertainties better than 0.5 ly. We examined the TESS light curves of these stars during the Ellipsoid crossing event and found no anomalous signatures. We discuss ways to expand this methodology to other surveys, more targets, and different potential signal types.

Using the upGREAT instrument on SOFIA, we have imaged [C II] 157.74 and [O I] 63.18 micron line emission from a bright photodissociation region (PDR) associated with an ionized ``bubble'' located in the Nessie Nebula, a filamentary infrared dark cloud. A comparison with ATCA data reveals a classic PDR structure, with a uniform progression from ionized gas, to photodissociated gas, and on to molecular gas from the bubble's interior to its exterior. [O I] line emission from the bubble's PDR reveals self-absorption features. Toward a FIR-bright protostar, both [O I] and [C II] show an absorption feature at a velocity of $-18$ km/s, the same velocity as an unrelated foreground molecular cloud. Since the gas density in typical molecular clouds is well below the [O I] and [C II] critical densities, the excitation temperatures for both lines are low (~20 K). The Meudon models demonstrate that the surface of a molecular cloud, externally illuminated by a standard G_0 = 1 interstellar radiation field, can produce absorption features in both transitions. Thus, the commonly observed [O I] and [C II] self-absorption and absorption features plausibly arise from the subthermally excited, externally illuminated, photodissociated envelopes of molecular clouds. The luminous young stellar object AGAL337.916-00.477, located precisely where the expanding bubble strikes the Nessie filament, is associated with two shock tracers: NH3 (3,3) maser emission and SiO 2-1 emission, indicating interaction between the bubble and the filament. The interaction of the expanding bubble with its parental dense filament has triggered star formation.

Context. Small bodies change their brightness due to different motives: Rotation along their axis or axes, combined with irregular shapes and/or changing surface properties, or changes in the geometry of observations. In this work, we tackle the problem of Phase curves, which show the change in brightness due to changes in the fraction of illuminated surface as seen by the observer. Aims. We aim to study the effect of the phase curves in the five wavelengths of the Sloan Digital Sky Survey in scores of objects (several tens of thousands), focusing particularly on the spectral slopes and the colors and their changes with phase angle. Methods. We used a Bayesian inference method and Monte Carlo techniques to retrieve the absolute magnitudes in five wavelengths, using the results to study the phase coloring effect in different bins of the semi-major axis. Results. We obtained absolute magnitudes in the five filters for over 40 000 objects. Although some outliers are identified, most of the usual color-color space is recovered by the data presented. We also detect a dual behavior in the spectral slopes, with a change at ${\alpha\approx}$ 5 deg.

Dense core collisions, previously regarded as minor in star formation, are proposed to play a significant role in structure formation around protostellar envelopes and binary formation. Using archival data of nearby star-forming regions, we determine the frequencies of core collisions. Our calculations reveal that a typical core is likely to undergo multiple interactions with other cores throughout its lifetime. To further investigate the core collision process, we employ adaptive mesh refinement hydrodynamic simulations with sink particles. Our simulations demonstrate that following the formation of a protostar within a gravitationally-unstable core, the merging core's accreting gas gives rise to a rotationally-supported circumstellar disk. Meanwhile, the region compressed by the shock between the cores develops into asymmetric arms that connect with the disk. Gas along these arms tends to migrate inward, ultimately falling toward the protostar. One of the arms, a remnant of the shock-compressed region, dominates over the second core gas, potentially exhibiting a distinct chemical composition. This is consistent with recent findings of large-scale streamers around protostars. Additionally, we found that collisions with velocities of $\sim$ 1.5 km s$^{-1}$ result in the formation of a binary system, as evidenced by the emergence of a sink particle within the dense section of the shocked layer. Overall, dense core collisions are highlighted as a critical process in creating $10^3$ au-scale streamers around protostellar systems and binary stars.

There is a heated debate regarding the specific roles played by ideal magnetohydrodynamic (MHD) instability and magnetic reconnection in the causes of solar eruptions. In the context with a pre-existing magnetic flux rope (MFR) before an eruption, it is widely believed that an ideal MHD instability, in particular, the torus instability, is responsible for triggering and driving the eruption, while reconnection, as invoked in the wake of the erupting MFR, plays a secondary role. Here we present a new numerical MHD model in which the eruption of a pre-existing MFR is primarily triggered and driven by reconnection. In this model, a stable MFR embedded in a strapping field is set as the initial condition. A surface converging flow is then applied at the lower boundary, pushing magnetic flux towards to the main polarity inversion line. It drives a quasi-static evolution of the system, during which a current layer is built up below the MFR with decreasing thickness. Once reconnection starts in the current sheet, the eruption commences, which indicates that the reconnection plays a determining role in triggers the eruption. By further analyzing the works done by in the magnetic flux of the pre-existing MFR and the newly reconnected flux during the acceleration stage of the eruption, we find that the latter plays a major role in driving the eruption. Such a model may explain observed eruptions in which the pre-eruption MFR has not reached the conditions for ideal instability.

Kink oscillations are ubiquitously observed in solar coronal loops, their understanding being crucial in the contexts of coronal seismology and atmospheric heating. We study kink modes supported by a straight coronal loop embeded in an asymmetric environment using three-dimensional magnetohydrodynamic (MHD) simulations. We implement the asymmetric effect by setting different exterior densities below and above the loop interior, and initiate the simulation using a kink-like velocity perturbation perpendicular to the loop plane, mimicking the frequently measured horizontally polarized kink modes. We find that the external velocity fields show fan blade structures propagating in the azimuthal direction as a result of the successive excitation of higher azimuthal Fourier modes. Resonant absorption and phase mixing can still occur despite an asymmetric environment, leading to the development of small scales at loop boundaries. These small scales nonetheless develop asymmetrically at the upper and lower boundaries due to the different gradients of the Alfven speed. These findings enrich our understanding of kink modes in coronal loops embedded within an asymmetric environment, providing insights helpful for future high-resolution observations.

In this work, we present the results from a study using the Giant Meterwave Radio Telescope (GMRT) to search for radio {emission} from planets around three evolved stars namely $\alpha$~Tau, $\beta$~UMi, and $\beta$~Gem. Both $\alpha$~Tau and $\beta$~UMi host massive $\sim$ 6 $M_J$ mass planets at about $\sim$1.4 au from the central star, while $\beta$~Gem is host to a 2.9 $M_J$ mass planet at 1.7 au from the host star. We observe $\alpha$~Tau and $\beta$~ UMi at two u(upgraded)GMRT bands; band~3 (250-500~MHz) and band~4 (550-900~MHz). We also analyzed the archival observations from $\beta$ Gem at 150~MHz from GMRT. We did not detect any radio signals from these systems. At 400~MHz, the 3$\sigma$ upper limit is 87 $\mu$Jy/beam for $\alpha$~Tau~{b} and 77.4 $\mu$Jy/beam for $\beta$~UMi~{b}. From our observations at 650~MHz, we place a 3$\sigma$ upper limit of 28.2 $\mu$Jy/beam for $\alpha$~Tau~b and 33.6 $\mu$Jy/beam for $\beta$~UMi~b. For $\beta$ Gem b, at 150~MHz, we place an upper limit of 2.5 mJy. At 400~MHz and 650~MHz, our observations are the deepest radio images for any exoplanetary system.

Herein, we present the statistical properties of active galactic nuclei (AGNs) for approximately 1 million member galaxies of galaxy groups and clusters, with 0.1 $<$ cluster redshift ($z_{\rm cl}$) $<$ 1.4, selected using Subaru Hyper Suprime-Cam, the so-called CAMIRA clusters. In this research, we focused on the AGN power fraction ($f_{\rm AGN}$), which is defined as the proportion of the contribution of AGNs to the total infrared (IR) luminosity, $L_{\rm IR}$ (AGN)/$L_{\rm IR}$, and examined how $f_{\rm AGN}$ depends on (i) $z_{\rm cl}$ and (ii) the distance from the cluster center. We compiled multiwavelength data using the ultraviolet--mid-IR range. Moreover, we performed spectral energy distribution fits to determine $f_{\rm AGN}$ using the CIGALE code with the SKIRTOR AGN model. We found that (i) the value of $f_{\rm AGN}$ in the CAMIRA clusters is positively correlated with $z_{\rm cl}$, with the correlation slope being steeper than that for field galaxies, and (ii) $f_{\rm AGN}$ exhibits a high value at the cluster outskirts. These results indicate that the emergence of AGN population depends on the redshift and environment and that galaxy groups and clusters at high redshifts are important in AGN evolution. Additionally, we demonstrated that cluster--cluster mergers may enhance AGN activity at the outskirts of particularly massive galaxy clusters. Our findings are consistent with a related study on the CAMIRA clusters that was based on the AGN number fraction.

The Mg II k \& h line intensity ratios can be used to probe the characteristics of the plasma in the solar atmosphere. In this study, using the observations recorded by the Interface Region Imaging Spectrometer (IRIS), we study the variation of the Mg II k \& h intensity ratio for three flares belonging to X-class, M-class, and C-class, throughout their evolution. We also study the k-to-h intensity ratio as a function of magnetic flux density obtained from the line-of-sight magnetograms recorded by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). Our results reveal that while the intensity ratios are independent of magnetic flux density, they show significant changes during the evolution of the C-class and M-class flares. The intensity ratios start to increase at the start of the flare and peak during the impulsive phase before the flare peak and decrease rapidly thereafter. The values of the ratios fall even below the pre-flare level during the peak and decline phases of the flare. These results are important in light of the heating and cooling of localized plasma and provide further constraints on the understanding of flare physics.

The resonant scattering interaction between Ly$\alpha$ photons and neutral hydrogen implies that a partially neutral IGM can significantly impact the detectability of Ly$\alpha$ emission in galaxies. The redshift evolution of the Ly$\alpha$ equivalent width distribution of galaxies thus offers a key probe of the degree of ionization during the Epoch of Reionization (EoR). Previous in-depth investigations at $z$ $\geq$ 7 were limited by ground-based instrument capabilities. We present an extensive study of Ly$\alpha$ emission from galaxies at 4 < $z$ < 8.5, observed from the CEERS and JADES surveys in the JWST NIRSpec/PRISM configuration. The sample consists of 235 galaxies, among which we identify 65 as Ly$\alpha$ emitters. We first measure Ly$\alpha$ escape fractions from Balmer lines, and explore the correlations with the inferred galaxies' physical properties, which are similar to those found at lower redshift. We also investigate the possible connection between the escape of Ly$\alpha$ photons and the inferred escape fractions of LyC photons obtained from indirect indicators. We then analyze the redshift evolution of the Ly$\alpha$ emitter fraction, finding lower average values at $z$ = 5 and 6 compared to ground-based observations. At $z$ = 7 we find a very large difference in Ly$\alpha$ visibility between the EGS and GOODS-South fields, possibly due to the presence of early reionized regions in the EGS. Such large variance is also expected in the Cosmic Dawn II radiation-hydrodynamical simulation. Our findings suggest a scenario in which the ending phase of the EoR is characterized by $\sim$ 1 pMpc ionized bubbles around a high fraction of moderately bright galaxies. Finally, we characterize such two ionized regions found in the EGS at $z$ = 7.18 and $z$ = 7.49 by estimating the radius of the ionized bubble that each of the spectroscopically-confirmed members could have created.

The magnetic field in the O9.7 V star \hd{} has been monitored for almost a decade. Spectropolarimetric observations reveal a rather strong mean longitudinal magnetic field that varies with a period of about 7.41 yr. Observations in the \halp{} line show a variation with the same period, while the \hbeta{} line shows only little variation. Assuming the periodic variation to be caused by a slow rotation and a dipolar magnetic field, we find a magnetic field strength of $\ge$2 kG at the magnetic poles. With the relatively low mass loss rate of $10^{-9} M_\odot$ yr$^{-1}$, this star is a case of extremely strong magnetic confinement. Both theoretical arguments and numerical simulations indicate the presence of an extended disk of increased gas density in the equatorial plane of the magnetic field, where gas from the line-driven stellar wind is trapped. This disk is likely to be the origin of the observed H$\alpha$ emission, which peaks together with the strongest line-of-sight magnetic field. The profile of the H$\alpha$ line is resolved in several components and shows a remarkable variability with the rotation period.

One of the major challenges we face is how to quickly and accurately create the three-dimensional (3D) density distributions of interstellar dust in the Milky Way using extinction and distance measurements of large samples of stars. In this study, we introduce a novel machine-learning approach that utilizes a convolution neural network, specifically a V-net, to infer the 3D distribution of dust density. Experiments are performed within two regions located towards the Galactic anti-center. The neural network is trained and tested using 10,000 simulations of dust density and line-of-sight extinction maps. Evaluation of the test sample confirms the successful generation of dust density maps from extinction maps by our model. Additionally, the performance of the trained network is evaluated using data from the literature. Our results demonstrate that our model is capable of capturing detailed dust density variations and can recover dust density maps while reducing the ``fingers of god" effect. Moving forward, we plan to apply this model to real observational data to obtain the fine distribution of dust at large and small scales in the Milky Way.

In this paper we proposed a joint reconstruction of \gray events using both extensive air array (EAS) and Imaging air Cherenkov Telescope array (IACT). We considered eight Cherenkov telescopes to be built on the LHAASO (Large High Altitude Air Shower Observatory) site and investigate the improvement in differential sensitivity when combining the information from both IACT and Moun detectors of LHAASO-KM2A. We found that due to the higher cosmic ray background rejection power and higher gamma ray retention ratio provided by muon detectors of LHAASO, such a joint reconstruction can significantly improve the sensitivity of IACTs, especially for extended sources and long exposure time. In this article, we have shown the performance of an eight-telescopes mini array, and our results indicate that above $10~\rm TeV$, the sensitivity can be improved by muon detector from $25\% - 60\%$ in different energy ranges.

Wavelet analysis, in addition to power density spectra, is another method to study the quasi-periodic signals in the light curves, but has been rarely used in black hole X-ray transients. We performed wavelet analysis of X-ray timing features and quasi-periodic oscillations (QPOs) based on NICER observations of the black hole candidate MAXI J1535-571 in this paper. Separating the light curves by the confidence level of wavelet results, we find significant differences exist in the PDS, hardness ratio and mean count between light curve segments above and below the confidence level. The S-factor, which is defined as the ratio of the effective oscillation time and the total time, demonstrates distinct values between type-C and type-B QPOs. Based on our results, the S-factor for type-B QPO is very close or equal to 0, no matter the confidence level is set as 95\% or 68\%, while the S-factor of type-C QPO is significantly higher, especially in the 68\% confidence level case. We discuss the implications of the wavelet results on resolving type-B and type-C QPOs in black hole X-ray binaries.

The second brightest gamma-ray burst, GRB 230307A (with a duration $T_{90}$ ~ 40 s), exhibited characteristics indicative of a magnetar engine during the prompt emission phase. Notably, a suspected kilonova was identified in its follow-up optical and infrared observations. Here we propose that the origin of GRB 230307A is a neutron star-white dwarf (NS-WD) merger, as this could naturally interpret the long duration and the large physical offset from the center of its host galaxy. In the framework of such a NS-WD merger event, the late-time kilonova-like emission is very likely to be powered by the spin-down of the magnetar and the radioactive decay of $^{56}$Ni, rather than by the decay of r-process elements as these heavy elements may not be easy to be synthesized in a NS-WD merger. It is demonstrated that the above scenario can be supported by our fit to the late-time observational data, where a mass of ~ $10^{-3} \ \rm M_{\odot}$ $^{56}$Ni is involved in the ejecta of a mass of ~ $0.1 \ \rm M_{\odot}$. Particularly, the magnetar parameters required by the fit are consistent with those derived from the early X-ray observation.

During recombination, the cosmic background radiation is disturbed, in particular, by Lyman-alpha emissions from neutral hydrogen. It is proposed to account for the subsequent time-dependent partial thermalization of the Lyman-alpha energy content in an analytically solvable nonlinear diffusion model. The amplitude of the partially thermalized and redshifted Ly-$\alpha$ line is found to be too low to be visible in the cosmic microwave spectrum, in accordance with previous numerical models and Planck observations.

We analyse the X-ray spectrum of the black hole X-ray binary MAXI J1820+070 using observations from XMM-Newton and NuSTAR during 'hard' states of its 2018-2019 outburst. We take a fully Bayesian approach, and this is one of the first papers to present a fully Bayesian workflow for the analysis of an X-ray binary X-ray spectrum. This allows us to leverage the relatively well-understood distance and binary system properties (like inclination and black hole mass), as well as information from the XMM-Newton RGS data to assess the foreground X-ray absorption. We employ a spectral model for a `vanilla' disc-corona system: the disc is flat and in the plane perpendicular to the axis of the jet and the black hole spin, the disc extends inwards to the innermost stable circular orbit around the black hole, and the (non-thermal) hard X-ray photons are up-scattered soft X-ray photons originating from the disc thermal emission. Together, these provide tight constraints on the spectral model and, in combination with the strong prior information about the system, mean we can then constrain other parameters that are poorly understood such as the disc colour correction factor. By marginalising over all the parameters, we calculate a posterior density for the black hole spin parameter, $a$. Our modelling suggests a preference for low or negative spin values, although this could plausibly be reproduced by higher spins and a modest degree of disc truncation. This approach demonstrates the efficacy and some of the complexities of Bayesian methods for X-ray spectral analysis.

The Parkes 20 cm Multibeam pulsar surveys have discovered nearly half of the known pulsars and revealed many distant pulsars with high dispersion measures. Using a sample of 1,301 pulsars from these surveys, we have explored the spatial distribution and birth rate of normal pulsars. The pulsar distances used to calculate the pulsar surface density are estimated from the YMW16 electron-density model. When estimating the impact of the Galactic background radiation on our survey, we projected pulsars in the Galaxy onto the Galactic plane, assuming that the flux density distribution of pulsars is uniform in all directions, and utilized the most up-to-date background temperature map. We also used an up-to-date version of the ATNF Pulsar Catalogue to model the distribution of pulsar flux densities at 1400 MHz. We derive an improved radial distribution for the pulsar surface density projected on to the Galactic plane, which has a maximum value at $\sim$4 kpc from the Galactic Centre. We also derive the local surface density and birthrate of pulsars, obtaining 47 $\pm$ 5 $\mathrm{kpc^{-2}}$ and $\sim$ 4.7 $\pm$ 0.5 $\mathrm{kpc^{-2}\ Myr^{-1}}$, respectively. For the total number of potentially detectable pulsars in the Galaxy, we obtain (1.1 $\pm$ 0.2) $\times$ $10^{4}$ and (1.1 $\pm$ 0.2) $\times$ $10^{5}$ before and after applying the TM98 beaming correction model. The radial distribution function is used to estimate the proportion of pulsars in each spiral arm and the Galactic centre.

The ring-like structures in protoplanetary discs that are observed in the cold dust emission by ALMA, might be explained by dust aggregates trapped aerodynamically in pressure maxima. The effect of a transient pressure maximum is investigated that develops between two regimes with different turbulent levels. We study how such a pressure maximum collects dust aggregates and transforms them into large planetesimals and Moon-mass cores that can further grow to a few Earth-mass planets by pebble accretion, and eventually to giant planets, by considering the accretion of a gaseous envelope. A numerical model is developed, incorporating the evolution of gaseous disc, growth and transport of pebbles, N-body interactions of growing planetary cores and their backreaction to gas disc by opening a partial gap. Planetesimal formation by streaming instability is parametrized in our model. A transient pressure maximum efficiently accumulates dust particles that can grow larger than mm-size. If this happens, dust aggregates can be transformed by the streaming instability process into such large planetesimals, which can grow further by pebble accretion, according to our assumptions. As the gas evolves to its steady state, the pressure maximum vanishes, and the concentrated pebbles that are not transformed to planetesimals and accreted by the growing planet, drift inward. During this inward drift, if the conditions of the streaming instability are met, planetesimals are formed in a wide radial range of the disc. Conclusions. A transient pressure maximum is a favourable place for planetesimal and planet formation during its lifetime and the concentration of pebbles induces continuous formation of planetesimals even after its disappearance. Besides, the formation of a planet can trigger the formation of planetesimals over a wide area of the protoplanetary disc.

We report the results of first ever spectro-polarimetric analyses of the Galactic ultra-luminous X-ray pulsar Swift J0243.6$+6124$ during the 2023 outburst using quasi-simultaneous {\it IXPE}, {\it NICER} and {\it NuSTAR} observations. A pulsation of period $\sim 9.79$ s is detected in {\it IXPE} and {\it NuSTAR} observations with pulse fractions (PFs) $\sim 18\%$ ($2-8$ keV) and $\sim 28\%$ ($3-78$ keV), respectively. Energy-dependent study of the pulse profiles with {\it NuSTAR} indicates an increase in PF from $\sim 27\%$ ($3-10$ keV) to $\sim 50\%$ ($40-78$ keV). Further, epoch-dependent polarimetric measurements during decay phase of the outburst confirm the detection of significant polarization degree (PD) varying as $\sim 2-3.1\%$ with polarization angle $\sim 8.6^{\circ}-10.8^{\circ}$ in $2-8$ keV energy range. We also observe that PD increases up to $\sim 4.8\%$ at higher energies ($\gtrsim 5$ keV) with dominating \texttt{bbodyrad} flux contribution ($1.5 \lesssim F_{\rm BB}/F_{\rm PL} \lesssim 3.4$) in {\it IXPE} spectra. The phase-resolved polarimetric study yields PD as $\sim 1.7-3.1\%$ suggesting a marginal correlation with the pulse profiles. Moreover, the broad-band ($0.6-70$ keV) energy spectrum of combined {\it NICER} and {\it NuSTAR} observations is well described by the combination of \texttt{bbodyrad} and \texttt{cutoffpl} components with seed photon temperature ($kT_{\rm bb}$) $\sim 0.86 \pm 0.03$ keV and photon index ($\Gamma$) $\sim 0.98 \pm 0.01$. With the above findings, we infer that the observed `low' PD in Swift J0243.6$+6124$ is resulted possibly due to `vacuum resonance' effect occurred between the overheated and relatively cooler regions of the neutron star boundary layer.

While stellar processes are believed to be the main source of feedback in dwarf galaxies, the accumulating discoveries of AGN in dwarf galaxies over recent years arouse the interest to also consider AGN feedback in them. Fast, AGN-driven outflows, a major mechanism of AGN feedback, have indeed been discovered in dwarf galaxies and may be powerful enough to provide feedback to their dwarf hosts. In this paper, we search for outflows traced by the blueshifted ultraviolet absorption features in three dwarf galaxies with AGN from the sample examined in our previous ground-based study. We confirm outflows traced by blueshifted absorption features in two objects and tentatively detect an outflow in the third object. In one object where the outflow is clearly detected in multiple species, photoionization modeling suggests that this outflow is located $\sim$0.5 kpc from the AGN, implying a galactic-scale impact. This outflow is much faster and possesses higher kinetic energy outflow rate than starburst-driven outflows in sources with similar star formation rates, and is likely energetic enough to provide negative feedback to its host galaxy as predicted by simulations. Much broader ($\sim$4000 km s$^{-1}$) absorption features are also discovered in this object which may have the same origin as that of broad absorption lines in quasars. Additionally, strong He II $\lambda$1640 emission is detected in both objects where the transition falls in the wavelength coverage, and is consistent with an AGN origin. In one of these two objects, blueshifted He II emission line is clearly detected, likely tracing a highly-ionized AGN wind.

For decades, the boundary of cosmic filaments have been a subject of debate. In this work, we determine the physically-motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at z = 0 in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from z = 4 to z = 1, and then slightly increases. Moreover, we find that the filament radius depends on the filament length, the distance from connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experience rapid contraction before z = 1, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves.

Observations of high-redshift quasars hosting billion solar mass black holes at $z\gtrsim6$ challenge our understanding of early supermassive black hole (SMBH) growth. In this work, we conduct a near-infrared spectroscopic study of $19$ quasars at $6.2\lesssim z\lesssim 7.5$, using the Folded-port InfraRed Echellette (FIRE) instrument on the $6.5$-meter Magellan/Baade Telescope. We estimate the single-epoch masses of the quasars' SMBHs by means of the MgII emission line and find black hole masses of $M_{\text{BH}} \approx(0.2-4.8)\,\times\,10^9\,M_\odot$. Furthermore, we measure the sizes of the quasars' proximity zones, which are regions of enhanced transmitted flux bluewards of the Ly$\alpha\,$ emission line, ionized by the quasars' radiation itself. While it has been shown that the proximity zone sizes correlate with the quasars' lifetimes due to the finite response time of the intergalactic medium to the quasars' radiation, we do not find any correlation between the proximity zone sizes and the black hole mass, which suggests that quasar activity and the concomitant black hole growth are intermittent and episodic.

We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe (EMU) early science observations, and Wide-field Infrared Survey Explorer (WISE), we compare this technique with a selection of different AGN diagnostics to explore the similarities and differences in AGN classification. We find that diagnostics based on optical and near-infrared criteria (the standard BPT diagram, the WISE colour criterion, and the mass-excitation, or MEx diagram) tend to favour detection of AGN in high mass, high luminosity systems, while the ``ProSpect'' SED fitting tool can identify AGN efficiently in low mass systems. We investigate an explanation for this result in the context of proportionally lower mass black holes in lower mass galaxies compared to higher mass galaxies and differing proportions of emission from AGN and star formation dominating the light at optical and infrared wavelengths as a function of galaxy stellar mass. We conclude that SED-derived AGN classification is an efficient approach to identify low mass hosts with low radio luminosity AGN.

We perform a comprehensive spectral study of a carefully selected sample (total exposure $\sim 50.5$ ks) of NICER observations of the atoll neutron star low mass X-ray binary 4U 1702-429. Our sample encompasses nearly all classical spectral states found within the NICER dataset. We require two thermal emission components, originating from the accretion disc and the boundary layer, to describe the soft state spectra in the energy band 0.3-10.0 keV. In contrast, in our model, only the disc component directly contributes to the intermediate/hard state. Additionally, we use a thermally Comptonised component (or a power-law with pegged normalisation) to represent the hard coronal emission in the soft and intermediate/hard state spectra. The boundary layer emerges as the principal source providing soft seed photons for Comptonisation across all spectral states. In contrast to a previously held assertion regarding this source, our analyses reveal a decrease in the inner disc temperature coupled with the retreat of the inner disc from the NS surface as the source evolves from the soft to the intermediate/hard state. The reflection features are either absent or weak ($\sim 3-4\sigma$) in all these observations. Further investigation using broad-band NuSTAR (3.0-50.0 keV) and AstroSat spectra (1.3-25.0 keV) shows a slightly stronger iron emission line ($\sim 5.8\sigma$) in the NuSTAR spectra. However, this feature is not significantly detected in the AstroSat observation. The AstroSat data suggests a highly ionised disc, explaining the absence of reflection features. In the case of NuSTAR, a truncated disc is likely responsible for the weak reflection features.

Thermal microwave emissions detected from stellar atmospheres contain information on stellar activity. However, even for the Sun, the relationship between multifrequency microwave data and other activity indices remains unclear. We investigated the relationships among the thermal microwave fluxes with 1, 2, 3.75 and 9.4 GHz, their circular polarizations, and several activity indices recorded during recent solar cycles and observed that these relationships can be categorized into two groups. In the first group, the relationship between the microwave fluxes and solar indices, which are strongly related to the active regions, can be well-fitted by using a linear function. In the second group, the fitting function is dependent on frequency. Specifically, the microwave fluxes at 1 and 2 GHz can be well-fitted to the total unsigned magnetic and extreme ultraviolet fluxes by employing a power-law function. The trend changes around 3.75 GHz, and of the trend for the 9.4 GHz fluxes can be fitted by using a linear function. For the first time, we present the relationship between circular polarization and solar indices. Moreover, we extrapolated these relationships of the solar microwave fluxes to higher values and compared them with the solar-type stars. We found that epsilon Eri, whose microwave emission originates from thermal plasma, follows the extrapolated relationship. However, to date, only one star's emission at 1--10 GHz has been confirmed as thermal emission. More solar-type stars should be observed with future radio interferometers to confirm that relationships based on solar data can be applied to stellar microwave data.

We present Submillimeter Array (SMA) mapping of $^{12}$CO $J=2\rightarrow 1$, $^{13}$CO $J=2\rightarrow 1$, and CN $N=2\rightarrow 1$ emission from the Ring-like planetary nebula (PN) NGC 3132, one of the subjects of JWST Early Release Observation (ERO) near-infrared imaging. The $\sim$5$''$ resolution SMA data demonstrate that the Southern Ring's main, bright, molecule-rich ring is indeed an expanding ring, as opposed to a limb-brightened shell, in terms of its intrinsic (physical) structure. This suggests that NGC 3132 is a bipolar nebula viewed more or less pole-on (inclination $\sim$15--30$^\circ$). The SMA data furthermore reveal that the nebula harbors a second expanding molecular ring that is aligned almost orthogonally to the main, bright molecular ring. We propose that this two-ring structure is the remnant of an ellipsoidal molecular envelope of ejecta that terminated the progenitor star's asymptotic giant branch evolution and was subsequently disrupted by a series of misaligned fast, collimated outflows or jets resulting from interactions between the progenitor and one or more companions.

Gamma-ray bursts (GRBs), especially short GRBs, are often considered potential candidates for exhibiting kilohertz quasi-periodic oscillations (QPOs) due to their origin from binary mergers. It has already been discovered that two bursts exhibit QPOs. While systematic searches for QPOs in GRBs typically concentrate on the kilohertz range, there has been no comprehensive exploration in the hundred-hertz range. In this study, we systematically conducted QPO searches on all BATSE short burst data within the 0-1000 Hz range. Using nested significance tests, we observed that the reference distributions for different GRBs are quite similar. This observation prompted us to analyze the data by selectively focusing on those with larger statistical values, obviating the need to iterate through all the data and significantly reducing computational workload. Ultimately, our findings did not reveal any compelling evidence for QPOs, which may suggest that the GRB jet has lost the early merging memory.

HESS J1813-178 is one of the brightest and most compact TeV $\gamma$-ray sources, and whether its $\gamma$-ray emission is associated with supernova remnant (SNR), pulsar wind nebula (PWN) or young stellar cluster (YSC) is still under debate. By analysing the GeV $\gamma$-ray data in the field of HESS J1813-178 using 14 years of PASS 8 data recorded by the Fermi Large Area Telescope (Fermi-LAT), we report the discovery of three different sources with different spectra in this region. The hard source with a power law spectral index of 2.11 $\pm$ 0.08 has a small size extension, which is spatially and spectrally coincident with the TeV $\gamma$-ray emission from HESS J1813-178. CO observations display the dense molecular clouds surrounding HESS J1813-178 in the velocity range of 45-60 km s$^{\rm -1}$. The possible origins of the $\gamma$-ray emission from HESS J1813-178 are discussed, including SNR G12.82-0.02, the PWN driven by the energetic X-ray pulsar PSR J1813-1749, and YSC Cl 1813-178. However, none of them can be ruled out clearly. Note that the maximum energy of protons in the hadronic model should exceed a few hundred TeV, which makes HESS J1813-178 to be a promising PeVatron. The detailed LHAASO data analysis about the morphology and spectrum would be helpful to investigate the origin of the $\gamma$-ray emission in this region and test its PeVatron nature.

We present the results of a spectral analysis using simultaneous multifrequency (22, 43, 86, and 129 GHz) very long baseline interferometry (VLBI) observations of the Korean VLBI Network (KVN) on BL Lac object, Markarian 421 (Mrk 421). The data we used was obtained from January 2013 to June 2018. The light curves showed several flux enhancements with global decreases. To separate the variable and quiescent components in the multifrequency light curves for milliarcsecond-scale emission regions, we assumed that the quiescent radiation comes from the emission regions radiating constant optically-thin synchrotron emissions (i.e., a minimum flux density with an optically thin spectral index). The quiescent spectrum determined from the multifrequency light curves was subtracted from the total CLEAN flux density, yielding a variable component in the flux that produces the time-dependent spectrum. We found that the observed spectra were flat at 22-43 GHz, and relatively steep at 43-86 GHz, whereas the quiescent-corrected spectra are sometimes quite different from the observed spectra (e.g., sometimes inverted at 22-43 GHz ). The quiescent-corrected spectral indices were much more variable than the observed spectral indices. This spectral investigation implies that the quiescent-spectrum correction can significantly affect the multifrequency spectral index of variable compact radio sources such as blazars. Therefore, the synchrotron self-absorption B-field strength (B_SSA) can be significantly affected because B_SSA is proportional to the fifth power of turnover frequency.

Mrk 501 is a prototypical high-synchrotron-peaked blazar (HBL) and serves as one of the primary targets for the Imaging X-ray Polarimetry Explorer (IXPE). In this study, we report X-ray polarization measurements of Mrk 501 based on six IXPE observations. The $>99\%$ confidence detection of X-ray polarization is obtained in four out of six observations in the entire energy band (2--8 keV) of IXPE, with a maximum polarization degree ($\Pi_{\rm X}$) of $15.8\%\pm2.8\%$ and a polarization angle ($\psi_{\rm X}$) of $98.0^{\circ}\pm5.1^{\circ}$ at a confidence level of $5.2\sigma$. During the other two observations, the source only exhibits the $>99\%$ confidence detection of X-ray polarization within an energy bin of 1 keV. Interestingly, both the second and sixth IXPE observations present energy-dependent variability in X-ray polarization, which is also observed in our optical spectropolarimetry. The chromatic behavior of $\Pi$ and the consistent values of $\psi$ across different frequencies, along with the agreement between $\psi$ and jet position angle, strongly support the interpretation of the energy-stratified model with shock-accelerated particles in the jet of Mrk 501. Additionally, the possibility of the presence of a global helical magnetic field in the jet of Mrk 501 is discussed.

The magnetic dipole moment $\mu$ of an accretion-powered pulsar in magnetocentrifugal equilibrium cannot be inferred uniquely from time-averaged pulse period and aperiodic X-ray flux data, because the radiative efficiency $\eta_0$ of the accretion is unknown, as are the mass, radius, and distance of the star. The degeneracy associated with the radiative efficiency is circumvented, if fluctuations of the pulse period and aperiodic X-ray flux are tracked with a Kalman filter, whereupon $\mu$ can be measured uniquely up to the uncertainties in the mass, radius, and distance. Here the Kalman filter analysis is demonstrated successfully in practice for the first time on Rossi X-ray Timing Explorer observations of the X-ray transient SXP 18.3 in the Small Magellanic Cloud, which is monitored regularly. The analysis yields $\mu = 8.0^{+1.3}_{-1.2} \, \times \, 10^{30} \, {\rm G \, cm^3}$ and $\eta_0 = 0.04^{+0.02}_{-0.01}$, compared to $\mu = 5.0^{+1.0}_{-1.0} \times 10^{30} \, {\rm G \, cm^3}$ as inferred traditionally from time-averaged data assuming $\eta_0=1$. The analysis also yields time-resolved estimates of two hidden state variables, the mass accretion rate and the Maxwell stress at the disk-magnetosphere boundary. The success of the demonstration confirms that the Kalman filter analysis can be applied in the future to study the magnetic moments and disk-magnetosphere physics of accretion-powered pulsar populations in the Small Magellanic Cloud and elsewhere.

Our analysis focus on the dragging effects on the accretion flows and jet emission in Kerr super-spinars. These attractors are characterized by peculiar accretion structures as double tori, or special dragged tori in the ergoregion, produced by the balance of the hydrodynamic and centrifugal forces and also effects of super-spinars repulsive gravity. We investigate the accretion flows, constituted by particles and photons, from toroids orbiting a central Kerr super-spinar. As results of our analysis, in both accretion and jet flows, properties characterizing these geometries, that constitute possible strong observational signatures or these attractors, are highlighted. We found that the flow is characterized by closed surfaces, defining inversion coronas (spherical shell), with null the particles flow toroidal velocity ($u^\phi=0$) embedding the central singularity. We proved that this region distinguishes proto-jets and accretion driven flows, co-rotating and counter-rotating flows. Therefore in both cases the flow carries information about the accretion structures around the central attractor, demonstrating that inversion points can constitute an observational aspect capable of distinguishing the super-spinars.

Recent studies reveal a radial acceleration relation (RAR) in galaxies, which illustrates a tight empirical correlation connecting the observational acceleration and the baryonic acceleration with a characteristic acceleration scale. However, a distinct RAR has been revealed on BCG-cluster scales with a seventeen times larger acceleration scale by the gravitational lensing effect. In this work, we systematically explored the acceleration and mass correlations between dynamical and baryonic components in 50 Brightest Cluster Galaxies (BCGs). To investigate the dynamical RAR in BCGs, we derived their dynamical accelerations from the stellar kinematics using the Jeans equation through Abel inversion and adopted the baryonic mass from the SDSS photometry. We explored the spatially resolved kinematic profiles with the largest integral field spectroscopy (IFS) data mounted by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. Our results demonstrate that the dynamical RAR in BCGs is consistent with the lensing RAR on BCG-cluster scales as well as a larger acceleration scale. This finding may imply that BCGs and galaxy clusters have fundamental differences from field galaxies. We also find a mass correlation, but it is less tight than the acceleration correlation.

We use spectral energy distribution (SED) fitting to place constraints on the stellar populations of 59 ultra-diffuse galaxies (UDGs) in the low-to-moderate density fields of the MATLAS survey. We use the routine PROSPECTOR, coupled with archival data in the optical from DECaLS, and near- and mid-infrared imaging from WISE, to recover the stellar masses, ages, metallicities and star formation timescales of the UDGs. We find that a subsample of the UDGs lies within the scatter of the mass-metallicity relation (MZR) for local classical dwarfs. However, another subsample is more metal-poor, being consistent with the evolving MZR at high-redshift. We investigate UDG positioning trends in the mass-metallicity plane as a function of surface brightness, effective radius, axis ratio, local volume density, mass-weighted age, star formation timescale, globular cluster (GC) counts and GC specific frequency. We find that our sample of UDGs can be separated into two main classes. Class A: Comprised of UDGs with lower stellar masses, prolonged star formation histories (SFHs), more elongated, inhabiting less dense environments, hosting fewer GCs, younger, consistent with the classical dwarf MZR, and fainter. Class B: UDGs with higher stellar masses, rapid SFHs, rounder, inhabiting the densest of our probed environments, hosting on average the most numerous GC systems, older, consistent with the high-redshift MZR (i.e., consistent with early-quenching), and brighter. The combination of these properties suggests that UDGs of Class A are consistent with a `puffed-up dwarf' formation scenario, while UDGs of Class B seem to be better explained by `failed galaxy' scenarios.

Disc winds and planet-disc interactions are two crucial mechanisms that define the structure, evolution and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms. We study the interplay between an X-ray photoevaporative disc wind and the substructures generated due to planet-disc interactions to determine how their mutual interactions affect the disc's and the planet's evolution. We perform three-dimensional hydrodynamic simulations of viscous ($\alpha = 6.9\cdot10^{-4}$) discs that host a Jupiter-like planet and undergo X-ray photoevaporation. We trace the gas flows within the disc and wind and measure the accretion rate onto the planet, as well as the gravitational torque that is acting on it. Our results show that the planetary gap takes away the wind's pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower, and the flow of mass across the gap in both directions is significantly increased, as well as the planet's mass-accretion rate (by factors $\approx 5$ and $\approx 2$, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes of the inner disc, contrary to the expectation in a planet-induced photoevaporation (PIPE) scenario that has been proposed in the past.

NGC 1566 is known for exhibiting recurrent outbursts, which are accompanied by changes in spectral type. The most recent transient event occurred from 2017 to 2019 and was reported to be accompanied by a change in Seyfert classification from Seyfert 1.8 to Seyfert 1.2. We analyze data from an optical spectroscopic variability campaign of NGC 1566 taken with the 9.2m SALT between July 2018 and October 2019 and supplement our data set with optical to near-infrared spectroscopic archival data taken by VLT/MUSE in September 2015 and October 2017. We observe the emergence and fading of a strong power-law-like blue continuum as well as strong variations in the Balmer, HeI, HeII lines and the coronal lines [FeVII], [FeX] and [FeXI]. Moreover, we detect broad double-peaked emission line profiles of OI 8446 and the CaII 8498,8542,8662 triplet. This is the first time that genuine double-peaked OI 8446 and CaII 8498,8542,8662 emission in AGN is reported in the literature. All broad lines show a clear redward asymmetry with respect to their central wavelength and we find indications for a significant blueward drift of the total line profiles during the transient event. We show that the double-peaked emission line profiles are well approximated by emission from a low-inclination, relativistic eccentric accretion disk, and that single-peaked profiles can be obtained by broadening due to scale-height dependent turbulence. Small-scale features in the OI and CaII lines suggest the presence of inhomogeneities in the broad-line region. We conclude that the broad-line region in NGC 1566 is dominated by the kinematics of a relativistic eccentric accretion disk. The broad-line region can be modeled to be vertically stratified with respect to scale-height turbulence. The observed blueward drift might be attributed to a low-optical-depth wind launched during the transient event.

Sizes of narrow emission line regions (NLRs) of AGN could be estimated by [O~{\sc iii}] line luminosity $L_{O3}$ through the known $R_{NLRs}-L_{O3}$ empirical relations. Unfortunately, it is not convenient to test the $R_{NLRs}-L_{O3}$ empirical relations through structure properties of spatially resolved NLRs of large samples of AGN. In this manuscript, a method is proposed to test the $R_{NLRs}-L_{O3}^{\sim0.25}$ empirical relations for AGN NLRs through SDSS Type-2 AGN having few orientation effects on NLRs sizes expected by AGN unified model, after considering sizes $R_{fib}$ of SDSS fiber covered regions. Comparing $R_{fib}$ and $R_{NLRs}$ estimated by $L_{O3}$, Type-2 AGN with $R_{fib}>R_{NLRs}$ (Sample-II) and with $R_{fib}<R_{NLRs}$ (Sample-I) should have different physical properties of NLRs. Accepted electron density gradients in AGN NLRs, statistically higher electron densities (traced by lower flux ratio $R_{S2}$ of [S~{\sc ii}]$\lambda6717$\AA~ to [S~{\sc ii}]$\lambda6731$\AA) could be expected for the Type-2 AGN in the Sample-I. Then, through the collected 1062 SDSS Type-2 AGN in the Sample-I and 3658 SDSS Type-2 AGN in the Sample-II, statistically lower $R_{S2}$ for the Type-2 AGN in the Sample-I can be confirmed with confidence level higher than 5$\sigma$, even after considering necessary effects. Therefore, the results in this manuscript can provide strong clues to support that the reported $R_{NLRs}~\propto~L_{O3}^{0.25}$ empirical relation is preferred to estimate NLRs sizes of SDSS AGN through SDSS fiber spectroscopic results, and also to support the commonly expected electron density gradients in AGN NLRs.

The fourth-DR3 version (4FGL-DR3) of the Fermi/LAT catalogue of $\gamma$-ray sources contains $\sim$ 1000 objects at a galactic latitude |b| > 10$^{\circ}$ which are not identified with an optical counterpart (UGS). We performed a systematic study of these sources, focusing on 190 objects that have a unique X-ray counterpart in the available Swift/XRT observations. Optical counterparts are then selected, and for 33 sources optical spectra were found in the literature. We found that 21 can be classified as BL Lac objects. Among these we were able to provide the redshift for 8 of them while for 2 others we established a lower limit to the redshift by detecting intervening absorption. The other 12 objects display optical spectra with prominent emission lines (0.036<z<1.65). These spectra are characterized by both broad and narrow emission lines with the exception of 3 sources. One of them displays only broad emission lines, while the other two exclusively exhibit narrow lines. On the basis of the radio/optical flux ratio, all BL Lac objects in this study are radio loud. Four sources out of the 12 with prominent emission lines can be classified as radio loud, while at least 5 of the 12 sources with prominent lines are radio quiet. This is somewhat unexpected comparing with the radio-loudness distribution of the 4FGL-associated blazars.

Fast radio bursts (FRBs) are bright, millisecond-duration radio bursts of cosmic origin. There have been several dozen FRBs found to repeat. Among them, those precisely localized provide the best opportunity to probe their multi-wavelength counterparts, local environment, and host galaxy that would reveal their origins. Here we report our X-ray, ultraviolet (UV) and optical observations with the $Swift$ satellite that were performed simultaneously in the radio band with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) observations of the repeating FRB 20190520B, aiming at detection of possible multi-wavelength bursts in association with radio bursts and multi-wavelength counterpart of the persistent radio source (PRS). While a total of 10 radio bursts were detected by FAST at the same time of $Swift$ observations, we detected neither X-ray, UV or optical bursts in accompany of the radio bursts, nor persistent multi-wavelength counterpart of the PRS. We obtained the energy upper limits ($3\sigma$) on any multi-wavelength bursts as $5.03 \times 10^{47}$ erg in the hard X-ray band (15-150 keV), $7.98 \times 10^{45}$ erg in the soft X-ray band (0.3-10 keV), and $4.51 \times 10^{44}$ erg in the U band, respectively. The energy ratio between soft X-ray (0.3-10 keV) and radio emission of the bursts is constrained as $<6\times10^{7}$, and the ratio between optical (U band) and radio as $<1.19\times10^{6}$. The 3$\sigma$ luminosity upper limits at the position of PRS are 1.04$\times10^{47}$ (15-150 keV), 8.81$\times10^{42}$ (0.3-10 keV), 9.26$\times10^{42}$ (UVW1), and 2.54$\times10^{42}$ erg s$^{-1}$ (U), respectively. We show that the PRS is much more radio loud than representative pulsar wind nebulae, supernova remnants, extended jet of Galactic X-ray binaries and ultraluminous X-ray sources, suggestive of boosted radio emission of the PRS.

The discovery of the first transiting hot Jupiters (HJs; giant planets on orbital periods shorter than $P\sim10$ days) was announced more than twenty years ago. As both ground- and space-based follow-up observations are piling up, we are approaching the temporal baseline required to detect secular variations in their orbital parameters. In particular, several recent studies focused on constraining the efficiency of the tidal decay mechanism to better understand the evolutionary time scales of HJ migration and engulfment. This can be achieved by measuring a monotonic decrease of orbital period $\mathrm{d}P/\mathrm{d}t<0$ due to mechanical energy being dissipated by tidal friction. WASP-12b was the first HJ for which a tidal decay scenario appeared convincing, even though alternative explanations have been hypothesized. Here we present a new analysis based on 28 unpublished high-precision transit light curves gathered over a twelve-year baseline and combined with all the available archival data, and an updated set of stellar parameters from HARPS-N high-resolution spectra, which are consistent with a main sequence scenario, close to the hydrogen exhaustion in the core. Our values of $\mathrm{d}P/\mathrm{d}t$ = $-30.72 \pm 2.67$ and $Q_{\ast}^{'}$ = $(2.13 \pm 0.18) \times 10^{5}$ are statistically consistent with previous studies, and indicate that WASP-12b is undergoing fast tidal dissipation. We additionally report the presence of an excess scatter in the timing data and discuss its possible origin.

Unraveling the internal kinematics of open clusters is crucial for understanding their formation and evolution. However, there is a dearth of research on this topic, primarily due to the lack of high-quality kinematic data. Using the exquisite-precision astrometric parameters and radial velocities provided by Gaia data release 3, we investigate the internal rotation in three of the most nearby and best-studied open clusters, namely the Pleiades, Alpha Persei, and Hyades clusters. Statistical analyses of the residual motions of the member stars clearly indicate the presence of three-dimensional rotation in the three clusters. The mean rotation velocities of the Pleiades, Alpha Persei, and Hyades clusters within their tidal radii are estimated to be 0.24 (0.04), 0.43 (0.08), and 0.09 (0.03) km s-1, respectively. Similar to the Praesepe cluster that we have studied before, the rotation of the member stars within the tidal radii of these three open clusters can be well interpreted by Newton's theorem. No expansion or contraction is detected in the three clusters either. Furthermore, we find that the mean rotation velocity of open clusters may be positively correlated with the cluster mass, and the rotation is likely to diminish as open clusters age.

Anisotropies in the distance-redshift relation of cosmological sources are expected due to large-scale inhomogeneities in the local Universe. When the observed sources are tracing a large-scale matter flow in a general spacetime geometry, the distance-redshift relation with its anisotropies can be described with a geometrical prediction that generalises the well-known Friedmann-Lema\^itre-Robertson-Walker result. Furthermore, it turns out that a finite set of multipole coefficients contain the full information about a finite-order truncation of the distance-redshift relation of a given observer. The multipoles of the distance-redshift relation are interesting new cosmological observables that have a direct physical interpretation in terms of kinematical quantities of the underlying matter flow. Using light cones extracted from $N$-body simulations we quantify the anisotropies expected in a $\Lambda$ cold dark matter cosmology by running a Markov chain Monte Carlo analysis on the observed data. In this observational approach the survey selection implements an implicit smoothing scale over which the effective rest frame of matter is fitted. The perceived anisotropy therefore depends significantly on the redshift range and distribution of sources. We find that the multipoles of the expansion rate, as well as the observer's velocity with respect to the large-scale matter flow, can be determined robustly with our approach.

In recent years the gas kinematics probed by molecular lines detected with ALMA has opened a new window to study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure and correctly interpreting these data allow us to gain a better comprehension of the planet formation process. We investigate the impact of thermal stratification on the azimuthal velocity of protoplanetary disks. High resolution gas observations are showing velocity differences between CO isotopologues, which cannot be adequately explained with vertically isothermal models. The aim of this work is to determine whether a stratified model can explain this discrepancy. We analytically solve the hydrostatic equilibrium for a stratified disk and we derive the azimuthal velocity. We test the model with SPH numerical simulations and then we use it to fit for star mass, disk mass and scale radius of the sources in the MAPS sample. In particular, we use 12CO and 13CO datacubes.

We have shown in a recent study, using 3D climate simulations, that dayside land cover has a substantial impact on the climate of a synchronously rotating temperate rocky planet such as Proxima Centauri b. Building on that result, we generate synthetic transit spectra from our simulations to assess the impact of these land-induced climate uncertainties on water vapour transit signals. We find that distinct climate regimes will likely be very difficult to differentiate in transit spectra, even under the more favourable conditions of smaller planets orbiting ultracool dwarfs. Further, we show that additional climate ambiguities arise when both land cover and atmosphere mass are unknown, as is likely to be the case for transiting planets. While water vapour may be detectable under favourable conditions, it may be nearly impossible to infer a rocky planet's surface conditions or climate state from its transit spectrum due to the interdependent effects of land cover and atmosphere mass on surface temperature, humidity, and terminator cloud cover.

In order to improve our understanding on the pre-requisites of eruptive solar flares, we study and compare different measures that characterize the eruptive potential of solar active regions - the critical height for torus instability as a local measure and the helicity ratio as a global measure - with the structural properties of the underlying magnetic field, namely the altitude of the center of the current-carrying magnetic structure. Using time series of 3D optimization-based nonlinear force-free magnetic field models for 10 different active regions (ARs) around the time of large solar flares, we determine the altitudes of the current-weighted centers of the non-potential model structures. Based on the potential magnetic field, we inspect the decay index, $n$, in multiple vertical planes oriented along of or perpendicular to the flare-relevant polarity inversion line, and estimate the critical height ($h_{\mathrm{crit}}$) for torus instability (TI) using different thresholds of $n$. The critical heights are interpreted with respect to the altitudes of the current-weighted centers of the associated non-potential structures, as well as the eruptive character of the associated flares, and the eruptive potential of the host AR, as characterized by the helicity ratio. Our most important findings are that (i) $h_{\mathrm{crit}}$ is more segregated in terms of flare type than the helicity ratio, and that (ii) coronal field configurations with a higher eruptive potential (in terms of the helicity ratio) also appear to be more prone to TI. Furthermore, we find no pronounced differences in the altitudes of the non-potential structures prior to confined and eruptive flares.

We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young sub-solar mass star Tau 042021, acquired as part of the Cycle 1 GO program "Mapping Inclined Disk Astrochemical Signatures (MIDAS)." These data resolve the mid-IR spatial distributions of H$_2$, revealing X-shaped emission extending to ~200 au above the disk midplane with a semi-opening angle of $35 \pm 5$ degrees. We do not velocity-resolve the gas in the spectral images, but the measured semi-opening angle of the H$_2$ is consistent with an MHD wind origin. A collimated, bipolar jet is seen in forbidden emission lines from [Ne II], [Ne III], [Ni II], [Fe II], [Ar II], and [S III]. Extended H$_2$O and CO emission lines are also detected, reaching diameters between ~90 and 190 au, respectively. Hot molecular emission is not expected at such radii, and we interpret its extended spatial distribution as scattering of inner disk molecular emission by dust grains in the outer disk surface. H I recombination lines, characteristic of inner disk accretion shocks, are similarly extended, and are likely also scattered light from the innermost star-disk interface. Finally, we detect extended PAH emission at 11.3 microns co-spatial with the scattered light continuum, making this the first low-mass T Tauri star around which extended PAHs have been confirmed, to our knowledge. MIRI MRS line images of edge-on disks provide an unprecedented window into the outflow, accretion, and scattering processes within protoplanetary disks, allowing us to constrain the disk lifetimes and accretion and mass loss mechanisms.

In their early, formative stages star clusters can undergo rapid dynamical evolution leading to strong gravitational interactions and ejection of ``runaway'' stars at high velocities. While O/B runaway stars have been well studied, lower-mass runaways are so far very poorly characterised, even though they are expected to be much more common. We carried out spectroscopic observations with MAG2-MIKE to follow-up 27 high priority candidate runaways consistent with having been ejected from the Orion Nebula Cluster (ONC) $>2.5$ Myr ago, based on Gaia astrometry. We derive spectroscopic youth indicators (Li \& H$\alpha$) and radial velocities, enabling detection of bona fide runaway stars via signatures of youth and 3D traceback. We successfully confirmed 10 of the candidates as low-mass Young Stellar Objects (YSOs) on the basis of our spectroscopic criteria and derived radial velocities (RVs) with which we performed 3D traceback analysis. Three of these confirmed YSOs have kinematic ejection ages $>4\:$Myr, with the oldest being 4.7~Myr. This yields an estimate for the overall formation time of the ONC to be at least $\sim 5\:$Myr, i.e., about 10 free-fall times, and with a mean star formation efficiency per free-fall time of $\bar{\epsilon}_{\rm ff}\lesssim0.05$. These results favor a scenario of slow, quasi-equilibrium star cluster formation, regulated by magnetic fields and/or protostellar outflow feedback.

Stellar coronal mass ejections (CMEs) from host stars are an important factor that affects the habitability of exoplanets. Although their solar counterparts have been well observed for decades, it is still very difficult to find solid evidence of stellar CMEs. Using the spectral line profile asymmetry caused by the Doppler shift of erupting plasma, several stellar CME candidates have been identified from spectral lines formed at chromospheric or transition region temperatures of the stars. However, a successful detection of stellar CME signals based on the profile asymmetries of coronal lines is still lacking. It is unclear whether we can detect such signals. Here we construct an analytical model for CMEs on solar-type stars, and derive an expression of stellar extreme-ultraviolet (EUV) line profiles during CME eruptions. For different instrumental parameters, exposure times, CME propagation directions and stellar activity levels, we synthesized the corresponding line profiles of Fe IX 171.07 \AA\ and Fe XV 284.16 \AA. Further investigations provide constraints on the instrumental requirements for successful detection and characterization of stellar CMEs. Our results show that it is possible to detect stellar CME signals and infer their velocities based on spectral profile asymmetries using an EUV spectrometer with a moderate spectral resolution and signal-to-noise ratio. Our work provides important references for the design of future EUV spectrometers for stellar CME detection and the development of observation strategies.

Stellar multiplicity is a key aspect of exoplanet diversity, as the presence of more than one star in a planetary system can have both devastating and positive effects on its formation and evolution. In this paper, we present the results of a Lucky Imaging survey of 212 exoplanet host stars performed with AstraLux at CAHA 2.2 m. The survey includes data from seven observing epochs between August 2015 and September 2020, and data for individual targets from four earlier observing epochs. The targets of this survey are nearby, bright, solar-like stars with high proper motions. In total, we detected 46 co-moving companions of 43 exoplanet host stars. Accordingly, this survey shows that the minimum multiplicity rate of exoplanet host stars is 20 $\pm$ 3 %. In total, 33 binary and ten hierarchical triple star systems with exoplanets have been identified. All companions were found to have a common proper motion with the observed exoplanet host stars, and with our astrometry we even find evidence of orbital motion for 28 companions. For all targets, we determined the detection limit and explore the detection space for possible additional companions of these stars. Based on the reached detection limit, additional co-moving companions beyond the detected ones can be excluded around all observed exoplanet host stars. The increasing number of exoplanets discovered in multiple stellar systems suggests that the formation of planets in such systems is by no means rare, but common. Therefore, our study highlights the need to consider stellar multiplicity in future studies of exoplanet habitability.

Earth's surface is deficient in available forms of many elements considered limiting for prebiotic chemistry. In contrast, many extraterrestrial rocky objects are rich in these same elements. Limiting prebiotic ingredients may, therefore, have been delivered by exogenous material; however, the mechanisms by which exogeneous material may be reliably and non-destructively supplied to a planetary surface remains unclear. Today, the flux of extraterrestrial matter to Earth is dominated by fine-grained cosmic dust. Although this material is rarely discussed in a prebiotic context due to its delivery over a large surface area, concentrated cosmic dust deposits are known to form on Earth today due to the action of sedimentary processes. Here we combine empirical constraints on dust sedimentation with dynamical simulations of dust formation and planetary accretion to show that localized sedimentary deposits of cosmic dust could have accumulated in arid environments on early Earth, in particular glacial settings that today produce cryoconite sediments. Our results challenge the widely held assumption that cosmic dust is incapable of fertilizing prebiotic chemistry. Cosmic dust deposits may have plausibly formed on early Earth and acted to fertilize prebiotic chemistry.

This study investigates accelerated cosmic expansion using various cosmological models within Horava-Lifshitz gravity, including BBCCFHKO, Seljak, ASSS, PADE-I, PADE-II, and BAZS, utilizing Equation of State Parametrization. To constrain the cosmological parameters of each model, we incorporate 24 Baryon Acoustic Oscillation points, 30 Cosmic Chronometer points, 40 Type Ia Supernovae points, 24 quasar Hubble diagram points, and 162 Gamma Ray Bursts points, along with the latest Hubble constant measurement (R22). We treat $r_{d}$ as a free parameter, aiming to extract $H_{0}$ and $r_{d}$ using late-time datasets and obtain optimal fitting values for each parameter in every model. The benefits of treating $r_{d}$ as a free parameter include reduced bias, improved precision, and enhanced dataset compatibility. The resulting values of $H_{0}$ and $r_{d}$ are relative to the $\Lambda$CDM model, showcasing alignment with early Planck and SDSS estimations. We additionally minimize errors by simulating random correlations in the covariance matrix. Furthermore, the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) support all models, with the $\Lambda$CDM model exhibiting the lowest AIC. Evaluation via reduced chi-square statistic confirms reasonable fits for all models. While $\Lambda$CDM remains favored, extensions warrant further investigation. This study underscores the importance of exploring alternative cosmological models to deepen our understanding of the universe's fundamental properties and evolution. Continued refinement of models is essential for advancing cosmological research.

We perform a detailed study of the black hole spin of Cyg X-1, using accurate broad-band X-ray data obtained in the soft spectral state by simultaneous NICER and NuSTAR observations, supplemented at high energies by INTEGRAL data. We use the relativistic disk model kerrbb together with different models of the Comptonization high energy tail and the relativistically-broadened reflection features. Unlike most previous studies, we tie the spin parameters of the disk and relativistic broadening models, thus combining the continuum and reflection methods of spin determination. We also consider a likely increase of the disk color correction due to a partial support of the disk by large scale magnetic fields. We find that such models yield the spin parameter of $a_*\approx 0.88^{+0.04}_{-0.01}$ if the disk inclination is allowed to be free, with $i\approx 40^{+1}_{-4}$ deg. Assuming $i=27.5$ deg, as determined by optical studies of the binary, worsens the fit but only marginally changes the spin, to $a_*\approx 0.92^{+0.07}_{-0.05}$. In addition, we consider the presence of a warm Comptonization layer on top of the disk, motivated by successful modeling of soft X-ray excesses in other sources with such a model. This dramatically lowers the spin, to $a_*\leq 0.3$, consistent with the spin measurements from black-hole mergers. On the other hand, if the natal spin of Cyg X-1 was low but now $a_*\approx 0.9$, a period of effective super-critical accretion had to take place in the past. Such accretion could be facilitated by photon advection, as proposed for ultraluminous X-ray sources.

Recent studies have suggested the possibility of Hycean worlds, characterised by deep liquid water oceans beneath H$_2$-rich atmospheres. These planets significantly widen the range of planetary properties over which habitable conditions could exist. We conduct internal structure modelling of Hycean worlds to investigate the range of interior compositions, ocean depths and atmospheric mass fractions possible. Our investigation explicitly considers habitable oceans, where the surface conditions are limited to those that can support potential life. The ocean depths depend on the surface gravity and temperature, confirming previous studies, and span 10s to $\sim$1000 km for Hycean conditions, reaching ocean base pressures up to $\sim$6$\times$10$^4$ bar before transitioning to high-pressure ice. We explore in detail test cases of five Hycean candidates, placing constraints on their possible ocean depths and interior compositions based on their bulk properties. We report limits on their atmospheric mass fractions admissible for Hycean conditions, as well as those allowed for other possible interior compositions. For the Hycean conditions considered, across these candidates we find the admissible mass fractions of the H/He envelopes to be $\lesssim$10$^{-3}$. At the other extreme, the maximum H/He mass fractions allowed for these planets can be up to $\sim$4-8$\%$, representing purely rocky interiors with no H$_2$O layer. These results highlight the diverse conditions possible among these planets and demonstrate their potential to host habitable conditions under vastly different circumstances to the Earth. Upcoming JWST observations of candidate Hycean worlds will allow for improved constraints on the nature of their atmospheres and interiors.

In this work, we consider the dynamics of the self-induced collapse of the tachyon wave function in inflationary scenarios. We analyze the modifications on the power spectrum by considering the \beta-exponential potential, whose parameters have updated constraints by the Planck 2018 baseline data. Moreover, we show that for this kind of potential, just for a narrow range of \beta-parameter, there is agreement between the theoretical predictions and the current observational data. Considering the proposal for a collapse scheme that leads to the modification of Schr\"odinger evolution of the inflaton wave function from the employment of a Continuous Spontaneous Localization (CSL) approach, we derive the scalar spectral index and tensor-to-scalar ratio. We then obtained the constraints on both collapse and \beta-parameters that, in turn, yield deviations in the n_{s} vs. r plane when compared to the \beta-exponential potential standard estimate. The CSL scheme applied to tachyonic inflation driven by a \beta-potential offers an adequate description of the recent data.

With the new generation of space telescopes such as the James Webb Space Telescope (JWST), it is possible to better characterize the atmospheres of exoplanets. The atmospheres of Hot and Ultra Hot Jupiters are highly heterogeneous and asymmetrical. The difference between the temperatures on the day-side and the night-side is especially extreme in the case of Ultra Hot Jupiters. We introduce a new tool to compute synthetic lightcurves from 3D GCM simulations, developed in the Pytmosph3R framework. We show how rotation induces a variation of the flux during the transit that is a source of information on the chemical and thermal distribution of the atmosphere. We find that the day-night gradient linked to Ultra Hot Jupiters has an effect close to the stellar limb-darkening, but opposite to tidal deformation. We confirm the impact of the atmospheric and chemical distribution on variations of the central transit time, though the variations found are smaller than that of available observational data, which could indicate that the east-west asymmetries are underestimated, due to the chemistry or clouds.

We present a study of X-ray normal galaxies using data from the first all-sky scan of the eROSITA X-ray survey. eRASS1 provides the first unbiased X-ray census of normal galaxies allowing us to study the X-ray emission from XRBs and the hot ISM in the full range of stellar population parameters present in the local Universe. By combining the HECATE value-added galaxy catalogue with the eRASS1, we study the X-ray emission from normal galaxies as a function of their SFR, M$_{*}$, Metallicity, and stellar population age. After applying optical and mid-IR activity classification criteria, we constructed a sample of 18790 star-forming galaxies with measurements of their L$_{X}$. By stacking the X-ray data in SFR-M$_{*}$-distance bins we study the correlation between the average L$_{X}$ and stellar population parameters. We also present updated L$_{\rm{X}}$-SFR and L$_{\rm{X}}$/SFR-Metallicity scaling relations accounting for the scatter dependence on the SFR. We find that the integrated L$_{X}$ of the HEC-eR1 star-forming galaxies is significantly elevated with respect to that expected from the current scaling relations. The observed scatter is also significantly larger. This excess persists even when we measure the average L$_{X}$ of galaxies in SFR-M$_{*}$-distance and metallicity bins and it is stronger in lower SFRs. The excess is not the result of hot gas, LMXBs, background AGN, LLAGN (including TDEs), or stochastic sampling of the XRB XLF. We find that while the excess correlates with lower metallicity, its primary driver is the age of the stellar populations. Our analysis reveals a sub-population of X-ray luminous starbursts with high sSFRs, low metallicities, and young stellar populations. This population drives upwards the X-ray scaling relations for star-forming galaxies, and has important implications for understanding the population of XRBs in the local and high-z Universe.

The Transiting Exoplanet Survey Satellite (TESS) mission measured light from stars in ~85% of the sky throughout its two-year primary mission, resulting in millions of TESS 30-minute cadence light curves to analyze in the search for transiting exoplanets. To search this vast dataset, we aim to provide an approach that is both computationally efficient, produces highly performant predictions, and minimizes the required human search effort. We present a convolutional neural network that we train to identify short period variables. To make a prediction for a given light curve, our network requires no prior target parameters identified using other methods. Our network performs inference on a TESS 30-minute cadence light curve in ~5ms on a single GPU, enabling large scale archival searches. We present a collection of 14156 short-period variables identified by our network. The majority of our identified variables fall into two prominent populations, one of short-period main sequence binaries and another of Delta Scuti stars. Our neural network model and related code is additionally provided as open-source code for public use and extension.

Recent observations have revealed detailed structures of radio relics in a wide range of frequencies. In this work, we perform three-dimensional magnetohydrodynamical simulations of merger shocks propagating through a turbulent magnetized intracluster medium, and employ on-the-fly Lagrangian particles to explore the physical processes originating radio substructures and their appearances in high and low-frequency observations. We employ two cosmic-ray (CR) electron acceleration models: the fresh injection of electrons from the thermal pool and the re-acceleration of mildly relativistic electrons. We use the relative surface brightness fluctuations, $\delta S_{\nu}$, to define a "degree of patchiness''. We find that: 1) Patchiness is produced if the shock's surface has a distribution of Mach numbers, rather than a single Mach number; 2) Radio relics appear patchier if the Mach number distribution consists of a large percentage of low Mach numbers ($\mathcal{M}\lesssim2.5$); 3) As the frequency increases, the patchiness also becomes larger. Nevertheless, if radio relics are patchy at high frequencies (e.g., 18.6 GHz), they are necessarily also at low frequencies (e.g., 150 MHz); 4) To produce noticeable differences in the patchiness at low and high frequencies, the shock front should have a Mach number spread of $\sigma_{\mathcal{M}}\gtrsim0.3$-0.4; 5) The amount of the patchiness depends on the Mach number distribution as well as the CR acceleration model. We propose $\delta S_{\nu}$ as a potential tool for extracting merger shock properties and information about particle acceleration processes at shocks in radio observations.

The holographic principle suggests that regions of space contain fewer physical degrees of freedom than would be implied by conventional quantum field theory. Meanwhile, in Hilbert spaces of large dimension $2^n$, it is possible to define $N \gg n$ Pauli algebras that are nearly anti-commuting (but not quite) and which can be thought of as "overlapping degrees of freedom". We propose to model the phenomenology of holographic theories by allowing field-theory modes to be overlapping, and derive potential observational consequences. In particular, we build a Fermionic quantum field whose effective degrees of freedom approximately obey area scaling and satisfy a cosmic Bekenstein bound, and compare predictions of that model to cosmic neutrino observations. Our implementation of holography implies a finite lifetime of plane waves, which depends on the overall UV cutoff of the theory. To allow for neutrino flux from blazar TXS 0506+056 to be observable, our model needs to have a cutoff $k_{\mathrm{UV}} \lesssim 500\, k_{\mathrm{LHC}}\,$. This is broadly consistent with current bounds on the energy spectrum of cosmic neutrinos from IceCube, but high energy neutrinos are a potential challenge for our model of holography. We motivate our construction via quantum mereology, \ie using the idea that EFT degrees of freedom should emerge from an abstract theory of quantum gravity by finding quasi-classical Hilbert space decompositions. We also discuss how to extend the framework to Bosons. Finally, using results from random matrix theory we derive an analytical understanding of the energy spectrum of our theory. The numerical tools used in this work are publicly available within the GPUniverse package, https://github.com/OliverFHD/GPUniverse .

We prove the non-linear stability of a large class of spherically symmetric equilibrium solutions of both the collisonless Boltzmann equation and of the Euler equations in MOND. This is the first such stability result that is proven with mathematical rigour in MOND. While we strive to prove our stability theorems, we develop new, genuinely Mondian ideas how arising mathematical difficulties can be solved. At some points it was necessary to restrict our analysis to spherical symmetry. We discuss every point where this extra assumption was necessary and outline which efforts must be undertaken to get along without it in future works. In the end we show on the example of a polytropic model how our stability result can be applied.

We apply the generalized Lomb-Scargle periodogram to 22 years data of solar $^{8}$B neutrino fluxes detected by Super-Kamiokande. The primary motivation of this work was to check if the sinusoidal modulation at a frequency of 9.43/year (with a period of 38 days), which we had found to be marginally significant with the first five years of Super-K data, persists, with the accumulated data. We use four different metrics for the calculation of significance. We do not find any evidence for periodicity at the aforementioned frequency or any other frequency with the updated data. Therefore the marginally detected periodicity at 9.43/year with the first five years of data was only a statistical fluctuation.

This review provides a conceptual and technical survey of methods for parameter estimation of gravitational wave signals in ground-based interferometers such as LIGO and Virgo. We introduce the framework of Bayesian inference and provide an overview of models for the generation and detection of gravitational waves from compact binary mergers, focusing on the essential features that are observable in the signals. Within the traditional likelihood-based paradigm, we describe various approaches for enhancing the efficiency and robustness of parameter inference. This includes techniques for accelerating likelihood evaluations, such as heterodyne/relative binning, reduced-order quadrature, multibanding and interpolation. We also cover methods to simplify the analysis to improve convergence, via reparametrization, importance sampling and marginalization. We end with a discussion of recent developments in the application of likelihood-free (simulation-based) inference methods to gravitational wave data analysis.

Evidence has accumulated that there are supermassive black holes (SMBHs) in the centers of most galaxies, and that these were formed in the very early universe by some as yet unknown process. In particular, there is evidence [15] that at least some galaxies formed as early as $10^8$ to $10^9$ years after the Big Bang host SMBHs. We suggest that the holographic model of inflation, whose dark matter candidates are primordial black holes carrying a discrete gauge charge, which originated as a small subset of the inflationary horizon volumes in the very early universe, can provide the seeds for this early structure formation. Aspects of the model pointed out long ago suggested an early era of structure formation, with structures dominated by dark matter. The additional assumption that the dark matter consists of discretely charged black holes implies black hole dominance of early structures, which seems to be implied by JWST data.

We consider the classic question posed by Pardo and Spergel about the price of abandoning dark matter in the context of an invariant, metric-based theory of gravity. Our answer is that the price is nonlocality. This has been known for some time in the context of the quasi-static regime. We show that it also applies for cosmology and we exhibit a model which reproduces standard CDM successes such as perturbations in the cosmic microwave background, baryon acoustic oscillations and structure formation.

Kinetic inductance traveling-wave parametric amplifiers (KI-TWPA) have a wide instantaneous bandwidth with near quantum-limited sensitivity and a relatively high dynamic range. Because of this, they are suitable readout devices for cryogenic detectors and superconducting qubits and have a variety of applications in quantum sensing. This work discusses the design, fabrication, and performance of a KI-TWPA based on four-wave mixing in a NbTiN microstrip transmission line. This device amplifies a signal band from 4 to 8~GHz without contamination from image tones, which are produced in a separate higher frequency band. The 4 - 8~GHz band is commonly used to read out cryogenic detectors, such as microwave kinetic inductance detectors (MKIDs) and Josephson junction-based qubits. We report a measured maximum gain of over 20 dB using four-wave mixing with a 1-dB gain compression point of -58 dBm at 15 dB of gain over that band. The bandwidth and peak gain are tunable by adjusting the pump-tone frequency and power. Using a Y-factor method, we measure an amplifier-added noise of $ 0.5 \leq N_{added} \leq 1.5$ photons from 4.5 - 8 GHz.

Sterile neutrino is a fascinating candidate for dark matter. In this paper, we examine the Affleck-Dine (AD) leptogenesis scenario generating a large lepton asymmetry, which can induce the resonant production of sterile neutrino dark matter via the Shi-Fuller (SF) mechanism. We also revisit the numerical calculation of the SF mechanism and the constraints from current X-ray and Lyman-$\alpha$ forest observations. We find that the AD leptogenesis scenario can explain the production of sterile neutrino dark matter by incorporating a non-topological soliton with a lepton charge called L-ball. Finally, we discuss an enhancement of second-order gravitational waves at the L-ball decay and investigate the testability of our scenario with future gravitational wave observations.

Nested sampling is widely used in astrophysics for reliably inferring model parameters and comparing models within a Bayesian framework. To address models with many parameters, Markov Chain Monte Carlo (MCMC) random walks are incorporated within nested sampling to advance a live point population. Diagnostic tools for nested sampling are crucial to ensure the reliability of astrophysical conclusions. We develop a diagnostic to identify problematic random walks that fail to meet the requirements of nested sampling. The distance from the start to the end of the random walk, the jump distance, is divided by the typical neighbor distance between live points, computed robustly with the MLFriends algorithm, to obtain a relative jump distance (RJD). We propose the geometric mean RJD and the fraction of RJD>1 as new summary diagnostics. Relative jump distances are investigated with mock and real-world inference applications, including inferring the distance to gravitational wave event GW170817. Problematic nested sampling runs are identified based on significant differences to reruns with much longer MCMC chains. These consistently exhibit low average RJDs and f(RJD>1) values below 50 percent. The RJD is more sensitive than previous tests based on the live point insertion order. The RJD diagnostic is proposed as a widely applicable diagnostic to verify inference with nested sampling. It is implemented in the UltraNest package in version 4.1.

Using the effective theory of dark energy (EFT) we show that the time evolution of non minimal gravity coupling can provide a natural explanation to the apparent Hubble tension. The non minimal coupling induces a modification of the Einstein frame Friedman equation, which can explain the difference between low and high red-shift estimations of $H_0$. Since the EFT predicts the non minimal coupling to be only a function of time, tests of the equivalence principle are insensitive to it, if experiments are performed in regions of space-time where the time scale of the coupling time evolution is much larger than the experiment time scale. The effects of a time varying non minimal gravity coupling only manifest on sufficiently long time scales, such as in cosmological observations at different red-shift, and if ignored lead to apparent tensions in the values of cosmological parameters estimated with observations from different epochs of the Universe history.