Papers for Wednesday, Feb 01 2023

The role of magnetic fields in galaxy evolution is still an unsolved question in astrophysics. We have previously shown that magnetic fields play a crucial role in major mergers between disc galaxies; in hydrodynamic simulations of such mergers, the Auriga model produces compact remnants with a distinctive bar and ring morphology. In contrast, in magnetohydrodynamic (MHD) simulations, remnants form radially-extended discs with prominent spiral arm structure. In this paper, we analyse a series of cosmological "zoom-in" simulations of major mergers and identify exactly $\textit{how}$ magnetic fields are able to alter the outcome of the merger. We find that magnetic fields modify the transport of angular momentum, systematically hastening the merger progress. The impact of this altered transport depends on the orientation of the field, with a predominantly azimuthal (non-azimuthal) orientation providing support against collapse (increasing the central baryonic concentration). Both effects act to suppress an otherwise existent bar-instability, which in turn leads to a fundamentally different morphology and manifestation of feedback. We note, in particular, that stellar feedback is substantially less influential in MHD simulations, which allows for the later accretion of higher angular momentum gas and the subsequent rapid radial growth of the remnant disc. A corollary of the increased baryonic concentration in MHD simulations is that black holes are able to grow twice as large, although this turns out to have little impact on the remnant's development. Our results show that galaxy evolution cannot be modelled correctly without including magnetic fields.

Frampton PLB 835 (2022) put forward the idea that charged, primordial extremely massive black holes between $10^{11}$ and $10^{22}$ solar masses could exist and may explain the observed accelerated expansion of the universe in lieu of dark energy. The extreme charge turns these black holes into naked singularities, making it difficult to derive observable signatures in their proximity. Here, we derive the electro-magnetic and gravitational lensing effects caused by such extreme objects at distances much larger than their extent and discuss the most promising observables to search for their existence in the least cosmology-dependent way. Restricting searches to black holes between $10^{12}$ to $10^{14}$ solar masses, we show that such objects do not cause any totally disruptive catastrophes, like dissociation of neutral hydrogen clouds or proton decay induced by strong electro-magnetic fields. Einstein rings of the order of 10 arcsec and rotation measures of plasma clouds subject to the magnetic fields induced by the black holes are identified as optimum observable signatures of these black holes in current sky surveys.

Five stars in the extreme outskirts (from $\sim5$ to $\sim12$ elliptical half-light radii, r$_h$) of the Ursa Minor (UMi) dwarf galaxy have been identified as potential new members using a Bayesian algorithm applied to \textit{Gaia} EDR3 data. These targets were observed with the GRACES spectrograph, resulting in precise radial velocities and metallicities that confirm their association with UMi. For the brightest and outermost star (Target~1, at $\sim12$ r$_h$), the chemical abundances of $\alpha$- (Mg, Ca, Ti), odd-Z (Na, K, Sc), Fe-peak (Fe, Ni, Cr), and neutron-capture process (Ba) elements have also been determined. We also discuss data from the literature and from APOGEE DR17. We find the chemical patterns in UMi are consistent with a star formation history that includes yields from core collapse supernovae, asymptotic giant branch stars, and supernovae Ia. Evidence for a knee in the [$\alpha$/Fe] ratios near [Fe/H] $\sim-2.1$ indicates a low star formation efficiency similar to that in other dwarf galaxies. Given the distance of Target~1 from the centre of UMi (R$\sim$4.5 kpc), we show that UMi has a more extended structure than previously thought. This "stellar halo" around UMi could be a secondary feature resulting from tidal stripping after multiple orbits around the Galaxy, or maybe a primary UMi feature due to early hierarchical accretion activity or to strong gravitational fluctuations prompted by feedback in the early star formation phase. Also consistent with observations is a late-time merger-free scenario where outside-in star formation is accompanied by gradual supernovae Ia enrichment.

Black hole evaporation is generally considered inevitable for low-mass black holes, yet there is no confirmation of this remarkable hypothesis. Here, we propose a phenomenological model that appeals to the possible survival of light quasi-extremal primordial black holes as a significant dark matter component and show that the related cosmological and astrophysical constraints disappear for reasonable degrees of quasi-extremality. The results obtained are general, conservative and should be taken as a proof of principle for future, model-specific analyses.

Galaxies and their dark-matter halos are commonly presupposed to spin. But it is an open question how this spin manifests in halos and soliton cores made of scalar dark matter (SDM, including fuzzy/wave/ultralight-axion dark matter). One way spin could manifest in a necessarily irrotational SDM velocity field is with a vortex. But recent results have cast doubt on this scenario, finding that vortices are generally unstable except with substantial repulsive self-interaction. In this paper, we introduce an alternative route to stability: in both (non-relativistic) analytic calculations and simulations, a black hole or other central mass at least as massive as a soliton can stabilize a vortex within it. This conclusion may also apply to stellar-scale Bose stars.

While interstellar gas is known to be supersonically turbulent, the injection processes of this turbulence are still unclear. Many studies suggest a dominant role of gravitational instabilities. However, their effect on galaxy morphology and large-scale dynamics vary across cosmic times, in particular due to the evolution of the gas fraction of galaxies. In this paper, we propose numerical simulations to follow the isothermal turbulent cascade of purely gravitationally-driven turbulence from its injection scale down to 0.095 pc for a gas-poor spiral disk and a gas-rich clumpy disk. To this purpose, and to lift the memory-footprint technical lock of sufficiently resolving the interstellar medium of a galaxy, we developed an encapsulated zoom method that allows us to probe self-consistently the self-generated turbulence cascade over three orders of magnitude on spatial scales. We follow this cascade for 10 Myrs. We find that the turbulent cascade follows the same scaling laws in both setups. Namely, in both cases the turbulence is close to equipartition between its compressive and solenoidal modes, the velocity power spectrum follows the Burgers' scaling and the density power spectrum is rather shallow, with a power-law slope of -0.7. Last, gravitationally-bound substructures follow a mass distribution with a -1.8 slope, similar to that of CO clumps. These simulations thus suggest a universality of gravity-driven isothermal turbulent cascade in disk galaxies across cosmic time.

The presence of an obscuring torus at pc-scale distances from the central black hole is the main ingredient for the Unified Model of Active Galactic Nuclei (AGN), as obscured sources are thought to be seen through this structure. However, the Unified Model fails to describe a class of sources that undergo dramatic spectral changes, transitioning from obscured to unobscured and vice-versa through time. The variability in such sources, so-called Changing Look AGN (CLAGN), is thought to be produced by a clumpy medium at much smaller distances than the conventional obscuring torus. ESO 323-G77 is a CLAGN that was observed in various states through the years with Chandra, Suzaku, Swift-XRT and XMM-Newton, from unobscured ($N_{\rm H}<3\times10^{22}$ cm$^{-2}$) to Compton-thin ($N_{\rm H}\sim1-6\times10^{23}$ cm$^{-2}$) and even Compton-thick ($N_{\rm H}>1\times10^{24}$ cm$^{-2}$), with timescales as short as one month. We present the analysis of the first NuSTAR monitoring of ESO 323-G77, consisting of 5 observations taken at different timescales (1, 2, 4 and 8 weeks from the first one) in 2016-2017, in which the AGN was caught in a persistent Compton-thin obscured state ($N_{\rm H}\sim2-4\times10^{23}$ cm$^{-2}$). We find that a Compton-thick reflector is present ($N_{\rm H,refl}=5\times10^{24}$ cm$^{-2}$), most likely associated with the presence of the putative torus. Two ionized absorbers are unequivocally present, located within maximum radii of $r_{\rm max,1}=1.5$ pc and $r_{\rm max,2}=0.01$ pc. In one of the observations, the inner ionized absorber is blueshifted, indicating the presence of a possible faster ($v_{\rm out}=0.2c$) ionized absorber, marginally detected at $3\sigma$. Finally, we are able to constrain the coronal temperature and the optical depth of ESO 323-G77, obtaining $kT_e=38$ keV or $kT_e=36$ keV, and $\tau=1.4$ or $\tau=2.8$, depending on the coronal geometry assumed.

Context. Strong lenses are a biased subset of the general population of galaxies. Aims. The goal of this work is to quantify how lens galaxies and lensed sources differ from their parent distribution, namely the strong lensing bias. Methods. We first studied how the strong lensing cross-section varies as a function of lens and source properties. Then, we simulated strong lensing surveys with data similar to that expected for Euclid and measured the strong lensing bias in different scenarios. We focused particularly on two quantities: the stellar population synthesis mismatch parameter, $\alpha_{sps}$, defined as the ratio between the true stellar mass of a galaxy and the stellar mass obtained from photometry, and the central dark matter mass at fixed stellar mass and size. Results. Strong lens galaxies are biased towards larger stellar masses, smaller half-mass radii and larger dark matter masses. The amplitude of the bias depends on the intrinsic scatter in the mass-related parameters of the galaxy population and on the completeness in Einstein radius of the lens sample. For values of the scatter that are consistent with observed scaling relations and a minimum detectable Einstein radius of $0.5''$, the strong lensing bias in $\alpha_{sps}$ is $10\%$, while that in the central dark matter mass is $5\%$. The bias has little dependence on the properties of the source population: samples of galaxy-galaxy lenses and galaxy-quasar lenses that probe the same Einstein radius distribution are biased in a very similar way. Quadruply imaged quasar lenses, however, are biased towards higher ellipticity galaxies. Conclusions. Given current uncertainties, strong lensing observations can be used directly to improve our current knowledge of the inner structure of galaxies, without the need to correct for selection effects.

X-ray radiation, in particular radiation between 0.1 keV and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the {\sc TreeRay} reverse raytrace algorithm implemented in the {\sc Flash} magneto-hydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV $< E_{\gamma} <$ 100 keV, dubbed {\sc XrayTheSpot}. {\sc XrayTheSpot} allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multi-bin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disk and simulations of molecular clouds with active chemistry, radiation pressure, protostellar radiation feedback from infrared to X-ray radiation.

We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC-pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion ($\sigma_v$) versus gas surface density ($\Sigma_{\mathrm{mol}}$) to the one observed in local spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For GMC-pixels in simulated mergers, we see a significant increase of factor of 5 - 10 in both $\Sigma_{\mathrm{mol}}$ and $\sigma_v$, which puts these pixels above the trend of PHANGS galaxies in the $\sigma_v$ vs $\Sigma_{\mathrm{mol}}$ diagram. This deviation may indicate that GMCs in the simulated mergers are much less gravitationally bound compared with simulated control galaxies with virial parameter ($\alpha_{\mathrm{vir}}$) reaching 10 - 100. Furthermore, we find that the increase in $\alpha_{\mathrm{vir}}$ happens at the same time as the increase in global star formation rate (SFR), which suggests stellar feedback is responsible for dispersing the gas. We also find that the gas depletion time is significantly lower for high $\alpha_{\mathrm{vir}}$ GMCs during a starburst event. This is in contrast to the simple physical picture that low $\alpha_{\mathrm{vir}}$ GMCs are easier to collapse and form stars on shorter depletion times. This might suggest that some other physical mechanisms besides self-gravity are helping the GMCs in starbursting mergers collapse and form stars.

This conference proceeding presents an overview of the modern approaches in the study of baryonic matter at high densities, focusing on the use of online repositories such as CompOSE and MUSES for the calculation of neutron star properties. In this context, relevant astrophysical constraints for the equations of state (mass-radius relation, speed of sound, tidal deformability) are discussed.

The phenomenon of pulsar nulling, observed as the temporary inactivity of a pulsar, remains poorly understood both observationally and theoretically. Most observational studies that quantify nulling employ a variant of Ritchings (1976)'s algorithm which can suffer significant biases for pulsars where the emission is weak. Using a more robust mixture model method, we study pulsar nulling in a sample of 22 recently discovered pulsars, for which we publish the nulling fractions for the first time. These data clearly demonstrate biases of the former approach and show how an otherwise non-nulling pulsar can be classified as having significant nulls. We show that the population-wide studies that find a positive correlation of nulling with pulsar period/characteristic age can similarly be biased because of the bias in estimating the nulling fraction. We use our probabilistic approach to find the evidence for periodicity in the nulls in a subset of three pulsars in our sample. In addition, we also provide improved timing parameters for 17 of the 22 pulsars that had no prior follow-up.

ALMA observations of the disk around HD 163296 have resolved a crescent-shape substructure at around 55 au, inside and off-center from a gap in the dust that extends from 38 au to 62 au. In this work we propose that both the crescent and the dust rings are caused by a compact pair (period ratio $\simeq 4:3$) of sub-Saturn-mass planets inside the gap, with the crescent corresponding to dust trapped at the $L_5$ Lagrange point of the outer planet. This interpretation also reproduces well the gap in the gas recently measured from the CO observations, which is shallower than what is expected in a model where the gap is carved by a single planet. Building on previous works arguing for outer planets at $\approx 86$ and $\approx 137$ au, we provide with a global model of the disk that best reproduces the data and show that all four planets may fall into a long resonant chain, with the outer three planets in a 1:2:4 Laplace resonance. We show that this configuration is not only an expected outcome from disk-planet interaction in this system, but it can also help constraining the radial and angular position of the planet candidates using three-body resonances.

We present a high time resolution, multi-frequency linear polarization analysis of Very Large Array (VLA) radio observations during some of the brightest radio flaring (~1 Jy) activity of the 2015 outburst of V404 Cygni. The VLA simultaneously captured the radio evolution in two bands (each with two 1 GHz base-bands), recorded at 5/7 GHz and 21/26 GHz, allowing for a broadband polarimetric analysis. Given the source's high flux densities, we were able to measure polarization on timescales of ~13 minutes, constituting one of the highest temporal resolution radio polarimetric studies of a black hole X-ray binary (BHXB) outburst to date. Across all base-bands, we detect variable, weakly linearly polarized emission (<1%) with a single, bright peak in the time-resolved polarization fraction, consistent with an origin in an evolving, dynamic jet component. We applied two independent polarimetric methods to extract the intrinsic electric vector position angles and rotation measures from the 5 and 7 GHz base-band data and detected a variable intrinsic polarization angle, indicative of a rapidly evolving local environment or a complex magnetic field geometry. Comparisons to the simultaneous, spatially-resolved observations taken with the Very Long Baseline Array at 15.6 GHz, do not show a significant connection between the jet ejections and the polarization state.

The initial mass function (IMF) of planetesimals is of key importance for understanding the initial stages of planet formation, yet theoretical predictions so far have been insufficient in explaining the variety of IMFs found in simulations. Here, we connect diffusion-tidal-shear limited planetesimal formation within the framework of a Toomre-like instability in the particle mid-plane of a protoplanetary disk to an analytic prediction for the planetesimal IMF. The shape of the IMF is set by the stability parameter $Q_\mathrm{p}$, which in turn depends on the particle Stokes number, the Toomre $Q$ value of the gas, the local dust concentration and the local diffusivity. We compare our prediction to high-resolution numerical simulations of the streaming instability and planetesimal formation via gravitational collapse. We find that our IMF prediction agrees with numerical results, and is consistent with both the `planetesimals are born big' paradigm and the power law description commonly found in simulations.

Blazar emission is dominated by nonthermal radiation processes that are highly variable across the entire electromagnetic spectrum. Turbulence, which can be a major source of nonthermal particle acceleration, can widely exist in the blazar emission region. The Turbulent Extreme Multi-Zone (TEMZ) model has been widely used to describe turbulent radiation signatures. Recent particle-in-cell (PIC) simulations have also revealed the stochastic nature of the turbulent emission region and particle acceleration therein. However, radiation signatures have not been systematically studied via first-principle-integrated simulations. In this paper, we perform combined PIC and polarized radiative transfer simulations to study synchrotron emission from magnetic turbulence in the blazar emission region. We find that the multi-wavelength flux and polarization are generally characterized by stochastic patterns. Specifically, the variability time scale and average polarization degree (PD) are governed by the correlation length of the turbulence. Interestingly, magnetic turbulence can result in polarization angle (PA) swings with arbitrary amplitudes and duration, in either directions, that are not associated to changes in flux or PD. Surprisingly, these swings, which are of stochastic nature, can appear either bumpy or smooth, although large amplitude swings ($>180^{\circ}$) are very rare as expected. Our radiation and polarization signatures from first-principle-integrated simulations are consistent with the TEMZ model.

Despite the importance of AGN in galaxy evolution, accurate AGN identification is often challenging, as common AGN diagnostics can be confused by contributions from star formation and other effects (e.g., Baldwin-Phillips-Terlevich diagrams). However, one promising avenue for identifying AGNs are ``coronal emission lines" (``CLs"), which are highly ionized species of gas with ionization potentials $\ge$ 100 eV. These CLs may serve as excellent signatures for the strong ionizing continuum of AGN. To determine if CLs are in fact strong AGN tracers, we assemble and analyze the largest catalog of optical CL galaxies using the Sloan Digital Sky Survey's Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) catalog. We detect CL emission in 71 MaNGA galaxies, out of the 10,010 unique galaxies from the final MaNGA catalog, with $\ge$ 5$\sigma$ confidence. In our sample, we measure [NeV]$\lambda$3347, $\lambda$3427, [FeVII]$\lambda$3586, $\lambda$3760, $\lambda$6086, and [FeX]$\lambda$6374 emission and crossmatch the CL galaxies with a catalog of AGNs that were confirmed with broad line, X-ray, IR, and radio observations. We find that [NeV] emission, compared to [FeVII] and [FeX] emission, is best at identifying high luminosity AGN. Moreover, we find that the CL galaxies with the least dust extinction yield the most iron CL detections. We posit that the bulk of the iron CLs are destroyed by dust grains in the galaxies with the highest [OIII] luminosities in our sample, and that AGN in the galaxies with low [OIII] luminosities are possibly too weak to be detected using traditional techniques.

Microquasar SS 433 located at the geometric center of radio nebula W50 is a suitable source for investigating the physical process of how galactic jets affect the surrounding interstellar medium (ISM). Previous studies have searched for evidence of the interaction between the SS 433 jet and ISM, such as neutral hydrogen gas and molecular clouds; however, it is still unclear which ISM interacts with the jet. We looked for new molecular clouds that possibly interact at the terminal of the SS 433 eastern jet using the Nobeyama 45-m telescope and the Atacama Submillimeter Telescope Experiment (ASTE). We identified two molecular clouds, comprising many small clumps, in the velocity range of 30.1--36.5 km s$^{-1}$ for the first time. These clouds have complex velocity structures, and one of them has a density gradient toward SS 433. Although it is difficult to conclude the relation between the molecular clouds and the SS 433/W50 system, there is a possibility that the eastern structure of W50 constructed by the SS 433 jet swept up tiny molecular clumps drifting in the surroundings and formed the molecular clouds that we identified in this study.

Alfven wave is the single most important physical phenomenon of magneto-hydrodynamic turbulence and has far-reaching impact to almost all studies related to astrophysical magnetic field. Yet the restoration of the Alfven wave fluctuations from a given magnetic field, aka the local Alfven wave problem, is never properly addressed in literature albeit its importance. Previous works model the Alfven wave fluctuation as the perturbation along a straight-line, constant magnetic field. However, Lazarian & Pogosyan (2012) suggested that the decomposition of Alfven wave along a straight line, aka. the global frame decomposition, has a factor of discrepancy to the true local Alfven wave fluctuation. Here we provide a geometric interpretation on how the local Alfven wave is related to the global frame through the use of vector frame formulation. We prove both analytically and numerically that the local frame Alfven wave is an orthogonal transformation of that of the global frame and related by the local Alfvenic Mach number. In other words, when we observe Alfven wave in the global frame of reference, some of the Alfven wave will be mistaken as compressible waves. The importance of frame choices have a far-reaching impact to the analytical studies of MHD turbulence. Combining the frame formalism and the new techniques we can have accurate measurement to some of the fundamental turbulence properties like the inclination angle of mean magnetic field relative to the line of sight.

We have carried out the first spectro-polarimetric study of the bright NS-LMXB GX 9+9 using IXPE and AstroSat observations. We report a significant detection of polarization of $1.7\pm 0.4\%$ over the $2-8$ keV energy band, with a polarization angle of $63^{\circ}\pm 7^{\circ}$. The polarization is found to be energy-dependent, with a $3\sigma$ polarization degree consistent with null polarization in $2-4$ keV, and $3.2\%$ in $4-8$ keV. Typical of the spectra seen in NS-LMXBs, we find that a combination of soft thermal emission from the accretion disc and Comptonized component from the optically thick corona produces a good fit to the spectra. We also attempt to infer the individual polarization of these components, and obtain a $3\sigma$ upper limit of $\sim 11\%$ on the polarization degree of the thermal component, and constrain that of the Comptonized component to $\sim 3\%$. We comment on the possible corona geometry of the system based on our results.

The shear measurement from DECaLS (Dark Energy Camera Legacy Survey) provides an excellent opportunity for galaxy-galaxy lensing study with DESI (Dark Energy Spectroscopic Instrument) galaxies, given the large ($\sim 9000$ deg$^2$) sky overlap. We explore this potential by combining the DESI 1\% survey and DECaLS DR8. With $\sim 106$ deg$^2$ sky overlap, we achieve significant detection of galaxy-galaxy lensing for BGS and LRG as lenses. Scaled to the full BGS sample, we expect the statistical errors to improve from $18(12)\%$ to a promising level of $2(1.3)\%$ at $\theta>8^{'}(<8^{'})$. This brings stronger requirements for future systematics control. To fully realize such potential, we need to control the residual multiplicative shear bias $|m|<0.01$ and the bias in the mean redshift $|\Delta z|<0.015$. We also expect significant detection of galaxy-galaxy lensing with DESI LRG/ELG full samples as lenses, and cosmic magnification of ELG through cross-correlation with low-redshift DECaLS shear. {If such systematical error control can be achieved,} we find the advantages of DECaLS, comparing with KiDS (Kilo Degree Survey) and HSC (Hyper-Suprime Cam), are at low redshift, large-scale, and in measuring the shear-ratio (to $\sigma_R\sim 0.04$) and cosmic magnification.

Galaxy shear - cosmic microwave background (CMB) lensing convergence cross-correlations contain additional information on cosmology to auto-correlations. While being immune to certain systematic effects, they are affected by the galaxy intrinsic alignments (IA). This may be responsible for the reported low lensing amplitude of the galaxy shear $\times$ CMB convergence cross-correlations, compared to the standard Planck $\Lambda$CDM (cosmological constant and cold dark matter) cosmology prediction. In this work, we investigate how IA affects the Kilo-Degree Survey (KiDS) galaxy lensing shear - Planck CMB lensing convergence cross-correlation and compare it to previous treatments with or without IA taken into consideration. More specifically, we compare marginalization over IA parameters and the IA self-calibration (SC) method (with additional observables defined only from the source galaxies) and prove that SC can efficiently break the degeneracy between the CMB lensing amplitude $A_{\rm lens}$ and the IA amplitude $A_{\rm IA}$. We further investigate how different systematics affect the resulting $A_{\rm IA}$ and $A_{\rm lens}$, and validate our results with the MICE2 simulation. We find that by including the SC method to constrain IA, the information loss due to the degeneracy between CMB lensing and IA is strongly reduced. The best-fit values are $A_{\rm lens}=0.84^{+0.22}_{-0.22}$ and $A_{\rm IA}=0.60^{+1.03}_{-1.03}$, while different angular scale cuts can affect $A_{\rm lens}$ by $\sim10\%$. We show that appropriate treatment of the boost factor, cosmic magnification, and photometric redshift modeling is important for obtaining the correct IA and cosmological results.

We present a simple phenomenological model of hadronic interaction between protons accelerated in an old supernova remnant (SNR) and cold protons situated within the associated molecular clouds (MCs). The accelerated protons from the old SNR escaped the SNR shock front, and got injected into the MCs at an earlier time, producing ultra high energy gamma-rays and neutrinos through inelastic proton-proton interaction. We also take into account the acceleration and subsequent escape of electrons from the SNR shock front. The escaped electrons produce gamma-rays through various radiative cooling mechanisms, after getting injected into the MCs. We use the model discussed in this letter to explain the multiwavelength (MWL) spectral energy distribution (SED) of unidentified Galactic ultra high energy gamma-ray source LHAASO J2108+5157. We also discuss the feasibility of applying this model in other cases as well. Future observations can test the viability of the model discussed in this letter, which will in turn confirm that the SNRs can, in fact, accelerate particles up to PeV energies.

Stellar coronal mass ejections (CMEs) have recently received much attention for their impacts on exoplanets and stellar evolution. Detecting prominence eruptions, the initial phase of CMEs, as the blue-shifted excess component of Balmer lines is a technique to capture stellar CMEs. However, most of prominence eruptions identified thus far have been slow and less than the surface escape velocity. Therefore, whether these eruptions were developing into CMEs remained unknown. In this study, we conducted simultaneous optical photometric observations with Transiting Exoplanet Survey Satellite and optical spectroscopic observations with the 3.8m Seimei Telescope for the RS CVn-type star V1355 Orionis that frequently produces large-scale superflares. We detected a superflare releasing $7.0 \times 10^{35} \: \mathrm{erg}$. In the early stage of this flare, a blue-shifted excess component of $\mathrm{H \alpha}$ extending its velocity up to $760-1690 \: \mathrm{km \: s^{-1}}$ was observed and thought to originate from prominence eruptions. The velocity greatly exceeds the escape velocity (i.e., $\sim 350 \: \mathrm{km \: s^{-1}}$), which provides important evidence that stellar prominence eruptions can develop into CMEs. Furthermore, we found that the prominence is very massive ($9.5 \times 10^{18} \: \mathrm{g} < M < 1.4 \times 10^{21} \: \mathrm{g}$). These data will clarify whether such events follow existing theories and scaling laws on solar flares and CMEs even when the energy scale far exceeds solar cases.

Recent development of the velocity gradient technique shows the {\toreferee capability} of the technique in the way of tracing magnetic fields morphology in diffuse interstellar gas and molecular clouds. In this paper, we perform the numerical systemic study of the performance of velocity and synchrotron gradient for a wide range of magnetization in the sub-sonic environment. Addressing the studies of magnetic field in atomic hydrogen, we also study the formation of velocity caustics in the spectroscopic channel maps in the presence of the thermal broadening. We show that the velocity caustics can be recovered when applied to the Cold Neutral Medium (CNM) and the Gradient Technique (GT) can reliably trace magnetic fields there. Finally, we discuss the changes of the anisotropy of observed structure functions when we apply to the analysis the procedures developed within the framework of GT studies.

Copious star formation occurs in the dense Central Molecular Zone (CMZ) of our Galaxy, but at a much smaller rate than occurs in a comparable mass of molecular gas in the Galactic disk. The combination of large turbulent velocity dispersions, a relatively strong magnetic field, and a strong tidal field all contribute to inhibiting star formation (SF) in different ways in different CMZ locations. Nonetheless, there are spectacular displays of recent and ongoing SF in the CMZ, including massive young stellar clusters, sites of abundant SF in progress, and numerous spots of protostellar or YSO activity. The presence of giant molecular clouds in the CMZ that are almost entirely devoid of SF indicates that SF requires a trigger that is not present everywhere. The dominant provocation of SF is likely to be cloud compression, either by large-scale shocks or by orbital motion of clouds into a region of enhanced tidal compression and/or enhanced external pressure. Recent hypotheses for where and how SF takes place in the CMZ are constrained by the recent orbital determinations of the massive Arches and Quintuplet clusters. Star formation in the central parsec is subject to a very different set of physical conditions, and is less well understood, but is important for the co-evolution of the central black hole and the nuclear star cluster.

AR Scorpii (AR Sco), the only-known radio-pulsing white dwarf binary, shows unusual pulsating emission at the radio, infrared, optical and ultraviolet bands. To determine its astrometric parameters at the radio band independently, we conducted multi-epoch Very Long Baseline Interferometry (VLBI) phase-referencing observations with the European VLBI Network (EVN) at 5 GHz and the Chinese VLBI Network (CVN) plus the Warkworth 30-metre telescope (New Zealand) at 8.6 GHz. By using the differential VLBI astrometry, we provide high-precision astrometric measurements on the parallax ($\pi=8.52_{-0.07}^{+0.04}$ mas), and proper motion ($\mu_{\alpha}=9.48_{-0.07}^{+0.04}$ mas yr$^{-1}$, $\mu_{\delta}=-51.32_{-0.38}^{+0.22}$ mas yr$^{-1}$). The new VLBI results agree with the optical Gaia astrometry. Our kinematic analysis reveals that the Galactic space velocities of AR Sco are quite consistent with that of both intermediate polars (IPs) and polars. Combined with the previous tightest VLBI constraint on the size, our parallax distance suggests that the radio emission of AR Sco should be located within the light cylinder of its white dwarf.

In the past decade, numerical simulations started to reveal the possible existence of planet 9 in our solar system. The planet 9 scenario can provide an excellent explanation to the clustering in orbital elements for Kuiper Belt objects. However, no optical counterpart has been observed so far to verify the planet 9 scenario. Therefore, some recent studies suggest that planet 9 could be a dark object, such as a primordial black hole. In this article, we show that the probability of capturing large trans-Neptunian objects (TNOs) by planet 9 to form a satellite system in the scattered disk region (between the inner Oort cloud and Kuiper Belt) is large. By adopting a benchmark model of planet 9, we show that the tidal effect can heat up the satellites significantly, which can give sufficient thermal radio flux for observations, even if planet 9 is a dark object. This provides a new indirect way for examining the planet 9 hypothesis and revealing the basic properties of planet 9.

We present in this reference paper an instrumental project dedicated to the monitoring of solar activity during solar cycle 25. It concerns the survey of fast evolving chromospheric events implied in Space Weather, such as flares, coronal mass ejections, filament instabilities and Moreton waves. Coronal waves are produced by large flares around the solar maximum and propagate with chromospheric counterparts; they are rare, faint, difficult to observe, and for that reason, challenging. They require systematic observations with automatic, fast and multi-channel optical instruments. MeteoSpace is a high cadence telescope assembly specially designed for that purpose. The large amount of data will be freely available to the solar community. We describe in details the optical design, the qualification tests and capabilities of the telescopes, and show how waves can be detected. MeteoSpace will be installed at Calern observatory (C{\^o}te d'Azur, 1270 m) and will be in full operation in 2023.

Fast radio bursts (FRBs) are being detected with increasing regularity. However, their spontaneous and often once-off nature makes high-precision burst position and frequency-time structure measurements difficult without specialised real-time detection techniques and instrumentation. The Australian Square Kilometre Array Pathfinder (ASKAP) has been enabled by the Commensal Real-time ASKAP Fast Transients Collaboration (CRAFT) to detect FRBs in real-time and save raw antenna voltages containing FRB detections. We present the CRAFT Effortless Localisation and Enhanced Burst Inspection pipeline (CELEBI), an automated software pipeline that extends CRAFT's existing software to process ASKAP voltages in order to produce sub-arcsecond precision localisations and polarimetric data at time resolutions as fine as 3 ns of FRB events. We use Nextflow to link together Bash and Python code that performs software correlation, interferometric imaging, and beamforming, making use of common astronomical software packages.

In evolved and dusty circumstellar discs, two planets with masses comparable to Jupiter and Saturn that migrate outwards while maintaining an orbital resonance can produce distinctive features in the dust distribution. Dust accumulates at the outer edge of the common gas gap, which behaves as a dust trap, where the local dust concentration is significantly enhanced by the planets outward motion. Concurrently, an expanding cavity forms in the dust distribution inside the planets orbits, because dust does not filter through the common gaseous gap and grain depletion in the region continues via inward drifting. There is no cavity in the gas distribution because gas can filter through the gap, although ongoing gas accretion on the planets can reduce the gas density in the inner disc. Such behaviour was demonstrated by means of simulations neglecting the effects of dust diffusion due to turbulence and of dust backreaction on the gas. Both effects may alter the formation of the dust peak at the gap outer edge and of the inner dust cavity, by letting grains filter through the dust trap. We performed high resolution hydrodynamical simulations of the coupled evolution of gas and dust species, the latter treated as pressureless fluids, in the presence of two giant planets. We show that diffusion and backreaction can change some morphological aspects of the dust distribution but do not alter some main features, such as the outer peak and the expanding inner cavity. These findings are confirmed for different parametrizations of gas viscosity.

M-dwarf stars are the dominant stellar population in the MilkyWay and they are important for a wide variety of astrophysical topics. The Gaia mission has delivered a superb collection of data, nevertheless, ground-based photometric surveys are still needed to study faint objects. Therefore, the present work aims to identify and characterise M-dwarf stars in the direction of the Galactic bulge using photometric data and with the help of Virtual Observatory tools. Using parallax measurements and proper motions from Gaia Data Release 3, in addition to different colour-cuts based on VISTA filters, we identify and characterise 7 925 M-dwarf stars in the b294 field from the Vista Variables in the V\'ia L\'actea (VVV) survey. We performed a spectral energy distribution fitting to obtain the effective temperature for all objects using photometric information available at Virtual Observatory archives. The objects in our sample have temperatures varying from 2800 to 3900 K. We also search for periodic signals in VVV light curves with up to 300 epochs, approximately. As a secondary outcome, we obtain periods for 82 M dwarfs by applying two methods: the Lomb-Scargle and Phase Dispersion Minimization methods, independently. These objects, with periods ranging from 0.14 to 34 d, are good candidates for future ground-based follow up. Our sample has increased significantly the number of known M dwarfs in the direction of the Galactic bulge and within 500 pc, showing the importance of ground-based photometric surveys in the near-infrared.

The last decade has seen tremendous developments in gamma-ray astronomy with the extragalactic sky becoming highly populated by Active Galactic Nuclei (AGNs). This brief review highlights some of the progress in AGN research achieved over the years, and discusses exemplary advances in the theory and physics of gamma-ray loud AGNs, including black-hole magnetospheric processes, the physics of pc-scales jets, as well as particle acceleration and high-energy emission in the large-scale jets of AGNs.

One of the few firm predictions of string theory is the existence of a massless scalar field coupled to gravity, the dilaton. In its presence, the value of the fundamental constants of the universe, such as the fine-structure constant, will vary with the time-dependent vacuum expectation value of this field, in direct violation of the Einstein Equivalence Principle. The \emph{runaway dilaton} proposed by Damour, Piazza, and Veneziano provides a physically motivated cosmological scenario which reconciles the existence of a massless dilaton with observations, while still providing non-standard and testable predictions. Furthermore, the field can provide a natural candidate for dynamical dark energy. While this model has been previously constrained from local laboratory experiments and low-redshift observations, we provide here the first full self-consistent constraints, also including high redshift data, in particular from the cosmic microwave background. We consider various possible scenarios in which the field could act as quintessence. Despite the wider parameter space, we make use of recent observational progress to significantly improve constraints on the model's coupling parameters, showing that order unity couplings (which would be natural in string theory) are ruled out.

We report the discovery with $TESS$ of a third set of eclipses from V994 Herculis (TIC 424508303), previously only known as a doubly-eclipsing system. The key implication of this discovery and our analyses is that V994 Her is the second fully-characterized (2+2) + 2 sextuple system, in which all three binaries eclipse. In this work, we use a combination of ground-based observations and $TESS$ data to analyze the eclipses of binaries A and B in order to update the parameters of the inner quadruple's orbit (with a derived period of 1062 $\pm$ 2d). The eclipses of binary C that were detected in the $TESS$ data were also found in older ground-based observations, as well as in more recently obtained observations. The eclipse timing variations of all three pairs were studied in order to detect the mutual perturbations of their constituent stars, as well as those of the inner pairs in the (2+2) core. At the longest periods they arise from apsidal motion, which may help constraining parameters of the component stars' internal structure. We also discuss the relative proximity of the periods of binaries A and B to a 3:2 mean motion resonance. This work represents a step forward in the development of techniques to better understand and characterize multiple star systems, especially those with multiple eclipsing components.

We present new neutral hydrogen (HI) observations of the nearby galaxy NGC 2403 to determine the nature of a low-column density cloud that was detected earlier by the Green Bank Telescope. We find that this cloud is the tip of a complex of filaments of extraplanar gas that is coincident with the thin disk. The total HI mass of the complex is $2\times10^{7}\text{ M}_\odot$ or 0.6% of the total HI mass of the galaxy. The main structure, previously referred to as the 8-kpc filament, is now seen to be even more extended, along a 20 kpc stream. The kinematics and morphological properties of the filaments are unlikely to be the result of outflows related to galactic fountains. It is more likely that the 20 kpc filament is related to a recent galaxy interaction. In this context, a $\sim$ 50 kpc long stellar stream has been recently detected connecting NGC 2403 with the nearby dwarf satellite DDO 44. Intriguingly, the southern tip of this stream overlaps with that of 20 kpc HI filament. We conclude that the HII anomalies in NGC 2403 are the result of a recent ($\sim2\text{ Gyr}$) interaction with DDO 44 leading to the observed filamentary complex.

In a study to provide ground truth data for mid infrared observations of the surface of Mercury with the MERTIS (Mercury Radiometer and Thermal Infrared Spectrometer) instrument onboard the ESA/JAXA BepiColombo mission, we have studied 17 synthetic glasses. These samples have the chemical compositions of characteristic Hermean surface areas based on MESSENGER data. The samples have been characterized using optical microscopy, EMPA and Raman spectroscopy. Mid infrared spectra have been obtained from polished thin sections using Micro FTIR, and of powdered size fractions of bulk material (0-25, 25-63, 93-125 and 125-250 micron) in the 2.5 to 18 micron range. The synthetic glasses display mostly spectra typical for amorphous materials with a dominating, single Reststrahlen Band (RB) at 9.5 micron to 10.7 micron. RB Features of crystalline forsterite are found in some cases at 9.5 to 10.2 micron, 10.4 to 11.2 micron, and at 11.9 micron. Dendritic crystallization starts at a MgO content higher than 23 wt.% MgO. The Reststrahlen Bands, Christiansen Features (CF), and Transparency Features (TF) shift depending on the SiO2 and MgO contents. Also a shift of the Christiansen Feature of the glasses compared with the SCFM (SiO2/(SiO2+CaO+FeO+MgO)) index is observed. This shift could potentially help distinguish crystalline and amorphous material in remote sensing data. A comparison between the degree of polymerization of the glass and the width of the characteristic strong silicate feature shows a weak positive correlation. A comparison with a high-quality mid-IR spectrum of Mercury shows some moderate similarity to the results of this study, but does not explain all features.

We investigate the structure and dynamics of the magnetospheric accretion region and associated outflows on a scale smaller than 0.1 au around the young transitional disk system GM Aur. We monitored the variability of the system on timescales ranging from days to months, using high-resolution optical and near-infrared spectroscopy, multiwavelength photometry, and low-resolution near-infrared spectroscopy, over a total duration of six months (30 rotational cycles). We analyzed the photometric and line profile variability to characterize the accretion and ejection processes. The luminosity of the system is modulated by surface spots at the stellar rotation period of 6.04 days. The Balmer, Paschen, and Brackett hydrogen lines as well as the HeI 5876 A and HeI 10830 A line profiles are modulated on the same period. The PaB line flux correlates with the photometric excess in the u' band, which suggests that most of the line emission originates from the accretion process. High-velocity redshifted absorptions reaching below the continuum periodically appear in the near-infrared line profiles at the rotational phase in which the veiling and line fluxes are the largest. These are signatures of a stable accretion funnel flow and associated accretion shock at the stellar surface. This large-scale magnetospheric accretion structure appears fairly stable over at least 15 and possibly up to 30 rotational periods. In contrast, outflow signatures randomly appear as blueshifted absorption components in the Balmer and HeI 10830 A line profiles and disappear on a timescale of a few days. The coexistence of a stable, large-scale accretion pattern and episodic outflows supports magnetospheric ejections as the main process occurring at the star-disk interface. Stable magnetospheric accretion and episodic outflows appear to be physically linked on a scale of a few stellar radii in this system.

Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting `Bringing Stellar Evolution and Feedback Together' in April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in.

PSR B1706$-$44 is an energetic gamma-ray pulsar located inside supernova remnant (SNR) G343.1$-$2.3 and it powers a compact pulsar wind nebula (PWN) that shows torus and jet structure in X-rays. We present a radio study of the PWN using Australia Telescope Compact Array (ATCA) observations at 3, 6, 13, and 21\,cm. We found an overall arc-like morphology at 3 and 6\,cm, and the ``arc" shows two distinct peaks at 6\,cm. The radio emission is faint inside the X-ray PWN and only brightens beyond that. We develop a thick torus model with Doppler boosting effect to explain the radio PWN structure. The model suggests a bulk flow speed of $\sim 0.2c$, which could indicate significant deceleration of the flow from the X-ray emitting region. Our polarization result reveals a highly ordered toroidal $B$-field in the PWN. Its origin is unclear given that the supernova reverse shock should have interacted with the PWN. At a larger scale, the 13 and 21\,cm radio images detected a semi-circular rim and an east-west ridge of G343.1$-$2.3. We argue that the latter could possibly be a pulsar tail rather than a filament of the SNR, as supported by the flat radio spectrum and the alignment between the magnetic field and its elongation.

Low-frequency radio observations have been expected to serve as a powerful tool for Space Weather (SW) observations for decades. Radio observations are sensitive to a wide range of SW-related observations ranging from emissions from coronal mass ejections (CMEs) to the solar wind. Ground-based radio observatories allow one gathering of high-sensitivity data at high time and spectral resolution for an extended period, which remains a challenge for most space-based observatories. While radio techniques like Interplanetary Scintillation (IPS) are well established, radio imaging studies have remained technically challenging. This is now changing with the confluence of data from instruments, like the Murchison Widefield Array (MWA), and robust unsupervised analysis pipelines. This pipeline delivers full Stokes radio images with unprecedented fidelity and dynamic range. This will serve as a powerful tool for coronal and heliospheric studies. We present the recent developments and achievements to measure the magnetic fields of the CME plasma and shock front at coronal heights and also share the current status of the objective to measure the heliospheric Faraday rotation towards numerous background linearly polarised radio sources with the Sun in the field of view. We envision that in the coming years, the availability of new-generation radio instruments combined with the Aditya-L1 and PUNCH mission will mark the start of a new era in Space Weather modeling and prediction.

The tomography of the polarized Sunyaev-Zeldvich effect due to free electrons of galaxy clusters can be used to constrain the nature of dark energy because CMB quadrupoles at different redshifts as the polarization source are sensitive to the integrated Sachs-Wolfe effect. Here we show that the low multipoles of the temperature and E-mode polarization anisotropies from the all-sky CMB can improve the constraint further through the correlation between them and the CMB quadrupoles viewed from the galaxy clusters. Using a Monte-Carlo simulation, we find that low multipoles of the temperature and E-mode polarization anisotropies potentially improve the constraint on the equation of state of dark energy parameter by $\sim 17$ percent.

Solar flares are efficient particle accelerators with a large fraction of released magnetic energy (10-50%) converted into energetic particles such as hard X-ray producing electrons. This energy transfer process is not well constrained, with competing theories regarding the acceleration mechanism(s), including MHD turbulence. We perform a detailed parameter study examining how various properties of the acceleration region, including its spatial extent and the spatial distribution of turbulence, affect the observed electron properties, such as those routinely determined from X-ray imaging and spectroscopy. Here, a time-independent Fokker-Planck equation is used to describe the acceleration and transport of flare electrons through a coronal plasma of finite temperature. Motivated by recent non-thermal line broadening observations that suggested extended regions of turbulence in coronal loops, an extended turbulent acceleration region is incorporated into the model. We produce outputs for the density weighted electron flux, a quantity directly related to observed X-rays, modelled in energy and space from the corona to chromosphere. We find that by combining several spectral and imaging diagnostics (such as spectral index differences or ratios, energy or spatial-dependent flux ratios, and electron depths into the chromosphere) the acceleration properties, including the timescale and velocity dependence, can be constrained alongside the spatial properties. Our diagnostics provide a foundation for constraining the properties of acceleration in an individual flare from X-ray imaging spectroscopy alone, and can be applied to past, current and future observations including those from RHESSI and Solar Orbiter.

Active galactic nuclei (AGNi) can launch and sustain powerful winds featuring mildly relativistic velocity and wide opening angle. Such winds, known as ultra-fast outflows (UFOs), can develop a bubble structure characterized by a forward shock expanding in the host galaxy and a wind termination shock separating the fast cool wind from the hot shocked wind. In this work we explore whether diffusive shock acceleration can take place efficiently at the wind termination shock of UFOs. We calculate the spectrum of accelerated particles and find that protons can be energized up to the EeV range promoting UFOs to promising candidates for accelerating ultra-high energy cosmic rays (UHECRs). We also compute the associated gamma-ray and neutrino fluxes and compare them with available data in the literature. We observe that high-energy (HE) neutrinos are efficiently produced up to hundreds of PeV while the associated gamma rays are efficiently absorbed beyond a few tens of GeV. By assuming a typical source density of non-jetted AGNi we expect that UFO could play a dominant role as diffuse sources of UHECRs and HE neutrinos. We finally apply our model to the recently observed NGC1068 and we find out that an obscured UFO could provide a sizeable contribution to the observed gamma-ray flux while only contributing up to $\sim 10\%$ to the associated neutrino flux.

In this paper, we present YunMa, a model which enables the study of cloud microphysics and radiative properties in exoplanetary atmospheres. YunMa simulates the vertical distribution and sizes of cloud particles and their corresponding scattering signature in transit spectra. We validated YunMa against results from the literature. When coupled to the TauREx 3 platform, an open Bayesian framework for spectral retrievals, YunMa enables the retrieval of the cloud properties and parameters from transit spectra of exoplanets. The sedimentation efficiency ($f_{sed}$), which controls the cloud microphysics, is set as a free parameter in retrievals. We assess the retrieval performances of YunMa through 28 instances of a cloudy super-Earth's atmosphere. This work also highlights the need for cloud radiative transfer and microphysics modelling to retrieve next-generation data of exoplanets.

A comprehensive wideband spectral analysis of the brightest black hole X-ray binary 4U $1543-47$ during its 2021 outburst is carried out for the first time using NICER, NuSTAR, and AstroSat observations by phenomenological and reflection modelling. The source attains a super-Eddington peak luminosity and remains in the soft state, with a small fraction ($< 3\%$) of the inverse-Comptonized photons. The spectral modelling reveals a steep photon index ($\Gamma \sim 2-2.6$) and relatively high inner disk temperature ($T_{in}\sim 0.9-1.27$ keV). The line-of-sight column density varies between ($0.45-0.54$)$\times10^{22}$ cm$^{-2}$. Reflection modelling using the RELXILL model suggests that 4U $1543-47$ is a low-inclination system ($\theta \sim 32^\circ - 40^\circ$). The accretion disk is highly ionized (log $\xi$ > 3) and has super solar abundance (3.6$-$10 $A_{Fe,\odot}$) over the entire period of study. We detected a prominent dynamic absorption feature between $\sim 8-11$ keV in the spectra throughout the outburst. This detection is the first of its kind for X-ray binaries. We infer that the absorption of the primary X-ray photons by the highly ionized, fast-moving disk-winds can produce the observed absorption feature. The phenomenological spectral modelling also shows the presence of a neutral absorption feature $\sim 7.1 - 7.4$ keV, and both ionized and neutral absorption components follow each other with a delay of a typical viscous timescale of $10-15$ days.

The inverse-Compton effect (IC) is a widely recognized cooling mechanism for both relativistic and thermal electrons in various astrophysical environments, including the intergalactic medium and X-ray emitting plasmas. Its effect on thermal electrons is however frequently overlooked in theoretical and numerical models of colliding-wind binaries (CWB). In this article, we provide a comprehensive investigation of the impact of IC cooling in CWBs, presenting general results for when the photon fields of the stars dominate the cooling of the thermal plasma and when shocks at the stagnation point are expected to be radiative. Our analysis shows that IC cooling is the primary cooling process for the shocked-wind layer over a significant portion of the relevant parameter space, particularly in eccentric systems with large wind-momentum ratios, e.g., those containing a Wolf-Rayet and O-type star. Using the binary system WR140 as a case study, we demonstrate that IC cooling leads to a strongly radiative shocked wind near periastron, which may otherwise remain adiabatic if only collisional cooling was considered. Our results are further supported by 2D and 3D simulations, which illustrate the impact of including or neglecting IC cooling. Specifically, 3D magnetohydrodynamic simulations of WR140 show a significant decrease in hard-X-ray emission around periastron, in agreement with observations but in contrast to equivalent simulations that omit IC cooling. A novel method is proposed for constraining mass-loss rates of both stars in eccentric binaries where the wind-collision zone switches from adiabatic to radiative approaching periastron. This study highlights the importance of including the IC process for thermal electrons in theoretical and numerical models of colliding-wind binaries.

We present the combination of the Canada-France-Hawaii Telescope (CHFT) Large Area $U$-bands Deep Survey (CLAUDS) and the Hyper-Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP) data over their four deep fields. We provide photometric catalogs for $u$, $u^*$ (CFHT--MegaCam), $g$, $r$, $i$, $z$, and $y$ (Subaru--HSC) bands over $\sim 20~{\rm deg}^2$, complemented in two fields by data from the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey and the UltraVISTA survey, thus extending the wavelength coverage toward near-infrared with VIRCAM $Y$, $J$, $H$, and $K_s$ observations over $5.5~{\rm deg}^2$. The extraction of the photometry was performed with two different softwares: the HSC pipeline hscPipe and the standard and robust SExtractor software. Photometric redshifts were computed with template-fitting methods using the new Phosphoros code for the hscPipe photometry and the well-known Le Phare code for the SExtractor photometry. The products of these methods were compared with each other in detail. We assessed their quality using the large spectroscopic sample available in those regions, together with photometry and photometric redshifts from COSMOS2020, the latest version of the Cosmic Evolution Survey catalogs. We find that both photometric data sets are in good agreement in $Ugrizy$ down to magnitude$\sim26$, and to magnitude$\sim24.5$ in the $YJHK_s$ bands. We achieve good performance for the photometric redshifts, reaching precisions of $\sigma_{NMAD} \lesssim 0.04$ down to ${m}_i\sim25$, even using only the CLAUDS and HSC bands. At the same magnitude limit, we measured an outlier fraction of $\eta \lesssim 10\%$ when using the $Ugrizy$ bands, and down to $\eta \lesssim 6\%$ when considering near-infrared data. [abridged]

We present the photometric analysis of \textit{BVR} and \textit{TESS} light curves of three eclipsing binaries (RU~UMi and purely studied VY~UMi, GSC 04364-00648), together with their period changes considering archival data and new minima times from our and \textit{TESS} observations. For the first time we detected wave-like variations with low-amplitude in $O-C$ residua of RU UMi, which can be interpreted as a consequence of the light-time effect caused by the 3rd invisible component with period 7370 days. Period increase with rate 2.56(9)$\times10^{-7}$~d/yr$^{-1}$ detected in the VY UMi system corresponds to mass transfer from the secondary to the primary component. For GSC 04364-00648 binary system we find some quadratic changes on the $O-C$ diagram, which corresponds to a period decrease with a high rate of $-2.26(5)\times10^{-5}$~d/yr$^{-1}$. We cannot assumptions about their nature, mainly due to short time of observation and uneven coverage of $O-C$ diagram. We also determined the absolute parameter of their components using the photometric solution and \textit{GAIA} distances.

The Near-Infrared Spectrograph (NIRSpec) is one of the four focal plane instruments on the James Webb Space Telescope. In this paper, we summarize the in-orbit performance of NIRSpec, as derived from data collected during its commissioning campaign and the first few months of nominal science operations. More specifically, we discuss the performance of some critical hardware components such as the two NIRSpec Hawaii-2RG (H2RG) detectors, wheel mechanisms, and the micro-shutter array. We also summarize the accuracy of the two target acquisition procedures used to accurately place science targets into the slit apertures, discuss the current status of the spectro-photometric and wavelength calibration of NIRSpec spectra, and provide the as measured sensitivity in all NIRSpec science modes. Finally, we point out a few important considerations for the preparation of NIRSpec science programs.

The Bayesian evidence is a key tool in model selection, allowing a comparison of models with different numbers of parameters. Its use in analysis of cosmological models has been limited by difficulties in calculating it, with current numerical algorithms requiring supercomputers. In this paper we give exact formulae for the Bayesian evidence in the case of Gaussian likelihoods with arbitrary correlations and top-hat priors, and approximate formulae for the case of likelihood distributions with leading non-Gaussianities (skewness and kurtosis). We apply these formulae to cosmological models with and without isocurvature components, and compare with results we previously obtained using numerical thermodynamic integration. We find that the results are of lower precision than the thermodynamic integration, while still being good enough to be useful.

In this work we report the independent discovery and analysis of nine new compact triply eclipsing triple star systems found with the TESS mission: TICs 47151245, 81525800, 99013269, 229785001, 276162169, 280883908, 294803663, 332521671, and 356324779. Each of these nine systems exhibits distinct third-body eclipses where the third (`tertiary') star occults the inner eclipsing binary (EB), or vice versa. We utilize a photodynamical analysis of the TESS photometry, archival photometric data, TESS eclipse timing variations of the EBs, available archival spectral energy distribution curves (SED), and, in some cases, newly acquired radial velocity observations, to solve for the parameters of all three stars, as well as most of the orbital elements. From these analyses we find that the outer orbits of all nine systems are viewed nearly edge on (i.e., within $\lesssim 4^\circ$), and 6 of the systems are coplanar to within $5^\circ$; the others have mutual inclination angles of $20^\circ$, $41^\circ$, and possibly $179^\circ$ (i.e., a retrograde outer orbit). The outer orbital periods range from 47.8 days to 604 days, with eccentricities spanning 0.004 to 0.61. The masses of all 18 EB stars are in the range of 0.9-2.6 M$_\odot$ and are mostly situated near the main sequence. By contrast, the masses and radii of the tertiary stars range from 1.4-2.8 M$_\odot$ and 1.5-13 R$_\odot$, respectively. We make use of the system parameters from these 9 systems, plus those from a comparable number of compact triply eclipsing triples published previously, to gain some statistical insight into their properties.

Numerical simulations of neutron star mergers represent an essential step toward interpreting the full complexity of multimessenger observations and constraining the properties of supranuclear matter. Currently, simulations are limited by an array of factors, including computational performance and input physics uncertainties, such as the neutron star equation of state. In this work, we expand the range of nuclear phenomenology efficiently available to simulations by introducing a new analytic parametrization of cold, beta-equilibrated matter that is based on the relativistic enthalpy. We show that the new $\textit{enthalpy parametrization}$ can capture a range of nuclear behavior, including strong phase transitions. We implement the enthalpy parametrization in the $\texttt{SpECTRE}$, code, simulate isolated neutron stars, and compare performance to the commonly used spectral and polytropic parametrizations. We find comparable computational performance for nuclear models that are well represented by either parametrization, such as simple hadronic EoSs. We show that the enthalpy parametrization further allows us to simulate more complicated hadronic models or models with phase transitions that are inaccessible to current parametrizations.

NectarCAM is a Cherenkov camera which is going to equip the Medium-Sized Telescopes (MST) of the northern site of the Cherenkov Telescope Array Observatory (CTAO). NectarCAM is equipped with 265 modules, each consisting of 7 photo-multiplier tubes (PMTs), a Front-End Board and a local camera trigger system used for data acquisition. This paper addresses the timing performances of NectarCAM which are crucial to reduce the noise in shower images and improve image cleaning as well as to discriminate between gamma-ray photons and cosmic-ray background and finally to allow coincidence identification with neighbouring telescopes for stereoscopic operations. Verification tests of the system have been performed in a dark room using various light sources to illuminate the first NectarCAM unit. The resulting timing precision and accuracy of the trigger arrival relative to a laser source, of individual and multiple pixel signals have been studied and are shown to comply to CTAO requirements.

It may be possible to detect biosignatures of photosynthesis in an exoplanet's atmosphere. However, such a detection would likely require a dedicated study, occupying a large amount of telescope time. It is therefore prudent, while searching for signs of life that we may recognise, to pick the best target possible. In this work, we present a new region, the ``photosynthetic habitable zone'' - the distance from a star where both liquid water and oxygenic photosynthesis can occur. It is therefore the region where detectable biosignatures of oxygenic photosynthesis are most likely to occur. Our analysis indicates that in the most ideal conditions for life and no atmospheric and greenhouse effects, the photosynthetic habitable zone is almost as broad as the habitable zone. On the other hand, if conditions for life are anything less than excellent and atmospheric attenuation and greenhouse effects are even moderate, the photosynthetic habitable zone is concentrated at larger separations around more massive stars. Such cases are also not tidally locked to their host star, which could result in planetary rotation periods similar to the Earth's. We identify five planets, Kepler-452 b, Kepler-1638 b, Kepler-1544 b and Kepler-62 e and Kepler-62 f, that are consistently in the photosynthetic habitable zone for a variety of conditions, and we predict their day lengths to be between 9 and 11 hours. We conclude that the parameter space in which we should search for signs of life is much narrower than the standard habitable zone.

Metric type II radio bursts are usually early indicators of CME-driven shocks and other space weather phenomena in the solar corona. This paper presents a detailed investigation of the spectral properties of band-splitting type II radio bursts and their association with sunspot number. Using type II radio bursts in a frequency range 20 to 200MHz observed by CALLISTO from 2010 to 2017, it was discovered that the analyzed type II shock height, magnetic field strength, CME shock speed, and Alfven speed synchronize with the trajectory of the solar cycle 24. Also, the study revealed that the onset of the declining phase of solar cycle 24 has the highest electron density. The analysis ascertained that the frequency of type II bursts depicts a bimodal distribution during the study period, with peaks in 2012 and 2015. Further, a good correlation(with correlation factor R=0.87) was obtained between the estimated CME shock speeds from the dynamic spectra, and the associated CME speeds from SOHO-LASCO. Moreover, the study confirmed a significant correlation(R=0.8) between the absolute drift rates and the plasma frequency. Additionally, the study explored that 60% of the type II radio bursts considered in this study emanated from the western longitudes. Hence, these findings emphasize that the temporal dynamics of the physical conditions of band-splitting type II radio are essential parameters in space weather monitoring and forecasting.

NANOGrav and other Pulsar Timing Arrays (PTAs) have discovered a common-spectrum process in the nHz range that may be due to gravitational waves (GWs): if so, they are likely to have been generated by black hole (BH) binaries with total masses $> 10^9 M_{\odot}$. Using the Extended Press-Schechter formalism to model the galactic halo mass function and a simple relation between the halo and BH masses suggests that these binaries have redshifts $z = {\cal O}(1)$ and mass ratios $\gtrsim 10$, and that the GW signal at frequencies above ${\cal O}(10)$~nHz may be dominated by relatively few binaries that could be distinguished experimentally and would yield observable circular polarization. Extrapolating the model to higher frequencies indicates that future GW detectors such as LISA and AEDGE could extend the PTA observations to lower BH masses $\in (10^6, 10^9 ) M_{\odot}$ and $\in (10^3, 10^9) M_{\odot}$.

We calculate deviations in cosmological observables as a function of parameters in a class of connection-based models of quantum gravity. In this theory non-trivial modifications to the background cosmology can occur due to a distortion of the wave function of the Universe at the transition from matter to dark energy domination (which acts as a "reflection" in connection space). We are able to exclude some regions of parameter space and show with projected constraints that future experiments like DESI will be able to further constrain these models. An interesting feature of this theory is that there exists a region of parameter space that could naturally alleviate the $S_8$ tension.

In this work, we reconstruct the cosmological unified dark fluid model proposed previously by Elkhateeb \cite{Elkhateeb:2017oqy} in the framework of $f(R)$ gravity. Utilizing the equivalence between the scalar-tensor theory and the $f(R)$ gravity theory, the scalar field for the dark fluid is obtained, whence the $f(R)$ function is extracted and its viability is discussed. The linear growth of matter perturbations on the low redshifts is studied in our constructed $f(R)$ function, and the results are in good agreement with those from the $\Lambda$CDM model for the dark universe. The ability of our function to describe early time inflation is also tested. The early time scalar field potential is extracted and the slow roll inflation parameters are derived. Results of the tensor-to-scalar ratio $r$ and the scalar spectral index $n_s$ are in good agreement with results from Planck-2018 TT+lowE data for the model parameter $m \ge 2$.

In the past years, evidence has started piling up that some high-energy cosmic neutrinos can be associated with blazars in flaring states. On February 26, 2022, a new blazar-neutrino coincidence has been reported: the track-like neutrino event IC220225A detected by IceCube is spatially coincident with the flat-spectrum radio quasar PKS 0215+015. Like previous associations, this source was found to be in a high optical and ${\gamma}$-ray state. Moreover, the source showed a bright radio outburst, which substantially increases the probability of a true physical association. We have performed six observations with the VLBA shortly after the neutrino event with a monthly cadence and are monitoring the source with the Effelsberg 100m-Telescope, and with the Australia Compact Telescope Array. Here, we present first results on the contemporary parsec-scale jet structure of PKS 0215+015 in total intensity and polarization to constrain possible physical processes leading to neutrino emission in blazars.

The discovery of transitional millisecond pulsars (tMSPs) provided conclusive proof that neutron star (NS) low-mass X-ray binaries (LMXBs) comprise part of the evolutionary pathway towards binary millisecond pulsars (MSPs). Redback and black widow `spider' pulsars are a sub-category of binary MSPs that `devour' their companions through ablation - the process through which material is lifted from the stellar surface by a pulsar wind. In addition to reducing the companion star's mass, ablation introduces observable characteristics like extended, energy-dependent and asymmetric eclipse profiles in systems observed at a sufficiently high inclination. Here, we present a detailed study and comparison of the X-ray eclipses of two NS LMXBs; $\textit{Swift}$ J1858.6$-$0814 and EXO 0748$-$676. Some of the X-ray eclipse characteristics observed in these two LMXBs are similar to the radio eclipse characteristics of eclipsing redback and black widow pulsars, suggesting that they may also host ablated companion stars. X-ray irradiation or a pulsar wind could drive the ablation. We conduct orbital phase-resolved spectroscopy for both LMXBs to map the column density, ionization and covering fraction of the material outflow. From this, we infer the presence of highly ionized and clumpy ablated material around the companion star in both systems. We term LMXBs undergoing ablation, $\textit{false widows}$, and speculate that they may be the progenitors of redback pulsars under the assumption that ablation begins in the LMXB stage. Therefore, the false widows could provide a link between LMXBs and spider pulsars. The detection of radio pulsations during non-accreting states can support this hypothesis.

Exploration of plasma dynamics in space, including turbulence, is entering a new era of multi-satellite constellation measurements that will determine fundamental properties with unprecedented precision. Familiar but imprecise approximations will need to be abandoned and replaced with more advanced approaches. We present a preparatory study of the evaluation of second- and third-order statistics, using simultaneous measurements at many points. Here, for specificity, the orbital configuration of the NASA Helioswarm mission is employed in conjunction with three-dimensional magnetohydrodynamics numerical simulations of turbulence. The Helioswarm 9-spacecraft constellation flies virtually through the turbulence to compare results with the exact numerical statistics. We demonstrate novel increment-based techniques for the computation of (1) the multidimensional spectra and (2) the turbulent energy flux. This latter increment-space estimate of the cascade rate, based on the third-order Yaglom-Politano-Pouquet theory, uses numerous increment-space tetrahedra. Our investigation reveals that Helioswarm will provide crucial information on the nature of astrophysical turbulence.

Superradiant instability of rotating black holes (BHs) leads to the formation of a cloud of ultralight bosons, such as axions. When the BH with the cloud belongs to a binary system and is in an inspiraling orbit, the resonant transition between the axion's bound states can occur. We study the history of the evolution of the binary system accompanying the cloud composed of the fastest growing mode, and its impact on the observational signatures, especially for small mass ratio cases. In this case, the hyperfine resonance, which has a very small resonance frequency, is relevant. Therefore, due to the long timescale, we should take into account the decaying process of axions in the transition destination mode, the backreaction to the orbital motion and the central BH, and gravitational emission from the cloud. We present a formulation to examine the evolution of the system around the resonance and useful expressions for the analysis. As a result, we found the mass of the cloud that can remain after the resonance is, at most, about $10^{-5}$ of the central BH. The maximum remaining cloud mass is achieved when the mass ratio of the binary is $q\sim10^{-3}$. In addition, we show that the resonant transition hardly changes the BH mass and spin distribution, while the associated modification of the gravitational wave frequency evolution when the binary pass through the resonance can be a signature of the presence of the cloud.

Besides the transient effect, the passage of a gravitational wave also causes a persistent displacement in the relative position of an interferometer's test masses through the nonlinear memory effect. This effect is generated by the gravitational backreaction of the waves themselves, and encodes additional information about the source. In this work, we explore the implications of using this information for the parameter estimation of massive binary black holes with LISA. Based on a Fisher analysis, our results show that the memory can help to reduce the degeneracy between the luminosity distance and the inclination for binaries observed only for a short time (~ few hours) before merger. To assess how many such short signals will be detected, we utilized state-of-the-art predictions for the population of massive black hole binaries and models for the gaps expected in the LISA data. We forecast from tens to few hundreds of binaries with observable memory, but only ~ O(0.1) events in 4 years for which the memory helps to reduce the degeneracy between distance and inclination. Based on this, we conclude that the new information from the non-linear memory, while promising for testing general relativity in the strong field regime, has probably a limited impact on further constraining the uncertainty on massive black hole binary parameters with LISA.

We investigate the possibility of the enhancement of parity-violation signal in bouncing cosmology. Specifically, we are interested in deciding which phase should generate the most significant parity-violation signals. We find that the dominant contribution comes from the bouncing phase, while the contraction phase has a smaller contribution. Therefore, bouncing cosmology can enhance the parity-violation signals during the bouncing phase. Moreover, since the bouncing phase has the highest energy scale in bouncing cosmology, we can also probe new physics at this scale by studying the parity-violation effect.

It has been argued that there exist extragalactic magnetic fields in the range from $10^{-17}$ to $10^{-9}$ Gauss with a cosmological coherence length. One plausible explanation for the origin of the extragalactic magnetic fields would be quantum fluctuations of the electromagnetic fields during inflation. It is also believed that primordial gravitational waves (PGWs) arise out of quantum fluctuations during inflation. We study the graviton-photon conversion process in the presence of background magnetic fields and find that the process induces the tachyonic instability of the PGWs. As a consequence, a peak appears in the power spectrum of PGWs. It turns out that the peak height depends on the direction of observation. The peak frequency could be in the range from $10^{-5}$ to $10^{3}$ Hertz for GUT scale inflation. Hence, the observation of PGWs could provide a new window for probing primordial magnetic fields.

The toolbox to study the Universe grew on 14 September 2015 when the LIGO-Virgo collaboration heard a signal from two colliding black holes between 30-250 Hz. Since then, many more gravitational waves have been detected as detectors increased sensitivity. However, the current detector design sensitivity curves still have a lower cut-off of 10 Hz. To detect even lower-frequency gravitational-wave signals, the Lunar Gravitational-wave Antenna will use an array of seismic stations in a permanently shadowed crater. It aims to detect the differential between the elastic response of the Moon and the suspended inertial sensor proof mass motion induced by gravitational waves. A cryogenic superconducting inertial sensor is under development that aims for fm/rtHz sensitivity or better down to 1 Hz and is planned to be deployed in seismic stations. Here, we describe the current state of research towards the inertial sensor, its applications and additional auxiliary technologies in the payload of the lunar gravitational-wave detection mission.

Any representation of data involves arbitrary investigator choices. Because those choices are external to the data-generating process, each choice leads to an exact symmetry, corresponding to the group of transformations that takes one possible representation to another. These are the passive symmetries; they include coordinate freedom, gauge symmetry and units covariance, all of which have led to important results in physics. Our goal is to understand the implications of passive symmetries for machine learning: Which passive symmetries play a role (e.g., permutation symmetry in graph neural networks)? What are dos and don'ts in machine learning practice? We assay conditions under which passive symmetries can be implemented as group equivariances. We also discuss links to causal modeling, and argue that the implementation of passive symmetries is particularly valuable when the goal of the learning problem is to generalize out of sample. While this paper is purely conceptual, we believe that it can have a significant impact on helping machine learning make the transition that took place for modern physics in the first half of the Twentieth century.

We study null and timelike constant radii geodesics in the environment of an over-spinning Kerr-type naked singularity. We are particularly interested in two topics: first, the differences of the shadows of the naked rotating singularity and the Kerr black hole; and second, the spinning down effect of the particles falling from the accretion disk. Our findings are as follows: around the naked singularity, the non-equatorial prograde orbits in the in Kerr black hole remain intact up to a critical rotation parameter ($\alpha=\frac{4\sqrt{2}}{3\sqrt{3}}$) and cease to exist above this value. This has an important consequence in the shadow of the naked singularity if the shadow is registered by an observer on the rotation plane or close to it as the shadow cannot be distinguished from that of a Kerr black hole viewed from the same angle. We also show that the timelike retrograde orbits in the equatorial plane immediately (after about an 8% increase in mass) reduce the spin parameter of the naked singularity from larger values to $\alpha=1$ at which an event horizon appears. This happens because the retrograde orbits have a larger capture cross-section than the prograde ones. So if a naked singularity happens to have an accretion disk, it will not remain naked for long, an event horizon forms.

The basic idea of this work is to achieve the observed relic density of a non-thermal dark matter(DM) and its connection with Cosmic Microwave Background (CMB) via additional relativistic degrees of freedom which are simultaneously generated during the period $T_{\rm BBN}~{\rm to}~T_{\rm CMB}$ from a long-lived dark sector particle. To realize this phenomena we minimally extend the type-I seesaw scenario with a Dirac fermion singlet($\chi$) and a complex scalar singlet ($\phi$) which transform non-trivially under an unbroken symmetry $\mathcal{Z}_3$. $\chi$ being the lightest particle in the dark sector acts as a stable dark matter candidate while the next to lightest state $\phi$ operates like a long lived dark scalar particle. The initial density of $\phi$ can be thermally produced through either self-interacting number changing processes ($3 \phi \to 2 \phi$) within dark sector or the standard annihilation to SM particles ($2 \phi \to 2~ {\rm SM}$). The late time (after neutrino decoupling) non-thermal decay of $\phi$ can produce dark matter in association with active neutrinos. The presence of extra relativistic neutrino degrees of freedom at the time of CMB can have a significant impact on $\Delta \rm N_{eff}$. Thus the precise measurement of $\Delta \rm N_{ eff}$ by current PLANCK 2018 collaboration and future experiments like SPT-3G and CMB-IV can indirectly probe this non-thermal dark matter scenario which is otherwise completely secluded due its tiny coupling with the standard model.