Papers for Friday, May 19 2023

At 7.7 pc, the A-type star Fomalhaut hosts a bright debris disk with multiple radial components. The disk is eccentric and misaligned, strongly suggesting that it is sculpted by interaction with one or more planets. Compact sources are now being detected with JWST, suggesting that new planet detections may be imminent. However, to confirm such sources as companions, common proper motion with the star must be established, as with unprecedented sensitivity comes a high probability that planet candidates are actually background objects. Here, ALMA and Keck observations of Fomalhaut are found to show significant emission at the same sky location as multiple compact sources in JWST MIRI coronagraphic observations, one of which has been dubbed the "Great Dust Cloud" because it lies within the outer belt. Since the ground-based data were obtained between 6 to 18 years prior to the JWST observations, these compact sources are unlikely to be common proper motion companions to Fomalhaut. More generally, this work illustrates that images collected at a range of wavelengths can be valuable for rejecting planet candidates uncovered via direct imaging with JWST.

Interloper contamination due to line misidentification is an important issue in the future low-resolution spectroscopic surveys. We realize that the algorithm previously used for photometric redshift self-calibration, with minor modifications, can be particularly applicable to calibrate the interloper bias. In order to explore the robustness of the modified self-calibration algorithm, we construct the mock catalogues based on China Space Station Telescope (CSST), taking two main target emission lines, H$\alpha$ and [O III]. The self-calibration algorithm is tested in cases with different interloper fractions at 1 per cent, 5 per cent and 10 per cent. We find that the interloper fraction and mean redshift in each redshift bin can be successfully reconstructed at the level of ~ 0.002 and ~ 0.001(1+z), respectively. We also find the impact of the cosmic magnification can be significant, which is usually ignored in previous works, and therefore propose a convenient and efficient method to eliminate it. Using the elimination method, we show that the calibration accuracy can be effectively recovered with slightly larger uncertainty.

Gravitational wave (GW) sources at cosmological distances can be used to probe the expansion rate of the Universe. GWs directly provide a distance estimation of the source but no direct information on its redshift. The optimal scenario to obtain a redshift is through the direct identification of an electromagnetic (EM) counterpart and its host galaxy. With almost 100 GW sources detected without EM counterparts (dark sirens), it is becoming crucial to have statistical techniques able to perform cosmological studies in the absence of EM emission. Currently, only two techniques for dark sirens are used on GW observations: the spectral siren method, which is based on the source-frame mass distribution to estimate conjointly cosmology and the source's merger rate, and the galaxy survey method, which uses galaxy surveys to assign a probabilistic redshift to the source while fitting cosmology. It has been recognized, however, that these two methods are two sides of the same coin. In this paper, we present a novel approach to unify these two methods. We apply this approach to several observed GW events using the \textsc{glade+} galaxy catalog discussing limiting cases. We provide estimates of the Hubble constant, modified gravity propagation effects, and population properties for binary black holes. We also estimate the binary black hole merger rate per galaxy to be $10^{-6}-10^{-5} {\rm yr^{-1}}$ depending on the galaxy catalog hypotheses.

We use hybrid (kinetic ions -- fluid electrons) kinetic simulations to investigate particle acceleration and magnetic field amplification at non-relativistic, weakly magnetized, quasi-perpendicular shocks. Unlike 2D simulations, 3D runs show that protons develop a non-thermal tail spontaneously (i.e., from the thermal bath and without pre-existing magnetic turbulence). They are rapidly accelerated via shock drift acceleration up to a maximum energy determined by their escape upstream.

We present the results of deep, ground based U-band imaging with the Large Binocular Telescope of the Cosmic Evolution Survey (COSMOS) field as part of the near-UV imaging program, UVCANDELS. We utilize a seeing sorted stacking method along with night-to-night relative transparency corrections to create optimal depth and optimal resolution mosaics in the U-band, which are capable of reaching point source magnitudes of AB 26.5 mag at 3 sigma. These ground based mosaics bridge the wavelength gap between the HST WFC3 F27W and ACS F435W images and are necessary to understand galaxy assembly in the last 9-10 Gyr. We use the depth of these mosaics to search for the presence of U-band intragroup light (IGrL) beyond the local Universe. Regardless of how groups are scaled and stacked, we do not detect any U-band IGrL to unprecedented U-band depths of 29.1-29.6 mag/arcsec2, which corresponds to an IGrL fraction of less than 1% of the total group light. This stringent upper limit suggests that IGrL does not contribute significantly to the Extragalactic Background Light at short wavelengths. Furthermore, the lack of UV IGrL observed in these stacks suggests that the atomic gas observed in the intragroup medium (IGrM) is likely not dense enough to trigger star formation on large scales. Future studies may detect IGrL by creating similar stacks at longer wavelengths or by pre-selecting groups which are older and/or more dynamically evolved similar to past IGrL observations of compact groups and loose groups with signs of gravitational interactions.

The recent detection of 12 gamma-ray Galactic sources well above E > 100 TeV by the LHAASO observatory has been a breakthrough in the context of Cosmic Ray (CR) origin search. Although most of these sources are unidentified, they are often spatially correlated with leptonic accelerators, like pulsar and pulsar wind nebulae (PWNe). This dramatically affects the paradigm for which a gamma-ray detection at E > 100 TeV implies the presence of a hadronic accelerator of PeV particles (PeVatron). Moreover, the LHAASO results support the idea that sources other than the standard candidates, Supernova Remnants, can accelerate Galactic CRs. In this context, the good angular resolution of future Cherenkov telescopes, such as the ASTRI Mini-Array and CTA, and the higher sensitivity of future neutrino detectors, such as KM3NeT and IceCube-Gen2, will be of crucial importance. In this brief review, we want to summarize the efforts done up to now, from both theoretical and experimental points of view, to fully understand the LHAASO results in the context of the CR acceleration issue.

Galaxies exhibiting a specific large-scale extended radio emission, such as X-shaped radio galaxies, belong to a rare class of winged radio galaxies. The morphological evolution of these radio sources is explained using several theoretical models, including galaxy mergers. However, such a direct link between a perturbed radio morphology and a galaxy merger remains observationally sparse. Here we investigate a unique radio galaxy J1159+5820, whose host CGCG 292-057 displays the optical signature of a post-merger system with a distinct tidal tail feature, and an X-shaped radio morphology accompanied by an additional pair of inner lobes. We observed the target on a wide range of radio frequencies ranging from 147 MHz to 4959 MHz, using dedicated GMRT and VLA observations, and supplemented it with publicly available survey data for broadband radio analysis. Particle injection models were fitted to radio spectra of lobes and different parts of the wings. Spectral ageing analysis performed on the lobes and the wings favors a fast jet realignment model with a reorientation timescale of a few million years. We present our results and discuss the possible mechanisms for the formation of the radio morphology.

Global stellar oscillations probe the internal structure of stars. In low- to intermediate-mass red giants, these oscillations provide signatures from both the outer regions of the star as well as from the core. These signatures are imprinted in e.g. the frequency of maximum oscillation power, and in the differences in periods of non-radial oscillations (period spacings), respectively. In core helium burning giants with masses below about 1.7 solar masses, i.e. stars that have gone through a helium flash, the asymptotic period spacings take values of about 220 -350 s at frequency of maximum oscillation power of $\sim$30-50 $\mu$Hz. A set of stars with asymptotic period spacings lower than about 200 s at similar frequencies separations has recently been discovered by Elsworth and collaborators. In this work, we present a hypothesis for the formation scenario of these stars. We find that these stars can be the result of a mass-loss event at the end of the red-giant branch phase of stars massive enough to not have a degenerate core, i.e. one of the scenarios to form hot subdwarf stars. Therefore, these stars can be classified as `hot subdwarf analogues'. Interestingly, if mass loss continues gradually during the core helium burning phase, these stars turn hotter and denser, and could, therefore, be hot subdwarf progenitors as they shed more of their envelope.

The environment of gamma-ray burst (GRB) has an important influence on the evolution of jet dynamics and of its afterglow. Here we investigate the afterglow polarizations in a stratified medium with the equal arrival time surface (EATS) effect. Polarizations of multi-band afterglows are predicted. The effects of the parameters of the stratified medium on the afterglow polarizations are also investigated. We found the influences of the EATS effect on the afterglow polarizations become important for off-axis detections and PD bumps move to later times with the EATS effect. Even the magnetic field configurations, jet structure and observational angles are fixed, polarization properties of the jet emission could still evolve. Here, we assume a large-scale ordered magnetic field in the reverse-shock region and a two-dimensional random field in the forward-shock region. Then PD evolution is mainly determined by the evolution of $f_{32}$ parameter (the flux ratio between the reverse-shock region and forward-shock region) at early stage and by the evolution of the bulk Lorentz factor $\gamma$ at late stage. Through the influences on the $f_{32}$ or $\gamma$, the observational energy band, observational angles, and the parameters of the stratified medium will finally affect the afterglow polarizations.

We construct a new formulation that allows efficient exploration of steady-state accretion processes onto compact objects. Accretion onto compact objects is a common scenario in astronomy. These systems serve as laboratories to probe the nuclear burning of the accreted matter. Conventional stellar evolution codes have been developed to simulate in detail the nuclear reactions on the compact objects. In order to follow the case of steady burning, however, using these codes can be very expensive as they are designed to follow a time-dependent problem. Here we introduce our new code $\textsc{StarShot}$, which resolves the structure of the compact objects for the case of stable thermonuclear burning, and is able to follow all nuclear species using an adaptive nuclear reaction network and adaptive zoning. Compared to dynamical codes, the governing equations can be reduced to time-independent forms under the assumption of steady-state accretion. We show an application to accreting low mass X-ray binaries (LMXBs) with accretion onto a neutron-star as compact object. The computational efficiency of $\textsc{StarShot}$ allows to to explore the parameter space for stable burning regimes, and can be used to generate initial conditions for time-dependent evolution models.

How the environment influences the most massive galaxies is still unclear. To explore the environmental effects on morphology and star formation in the most massive galaxies at high redshift, we select galaxies with stellar mass $\log(M_{\star}/M_{\odot})>11$ at $0.5<z<2.5$ in the COSMOS-DASH field, which is the largest field with near-infrared photometrical observations using HST/WFC3 to date. Combining with the newly published COSMOS2020 catalog, we estimate the localized galaxy overdensity using a density estimator within the Bayesian probability framework. With the overdensity map, no significant environmental dependence is found in the distributions of S\'{e}rsic index and effective radius. When we consider the star formation state, galaxies in lower density are found to have higher median specific star formation rate (sSFR) at $0.5<z<1.5$. But for star-forming galaxies only, sSFR is independent of the environment within the whole redshift range, indicating that the primary effect of the environment might be to control the quiescent fraction. Based on these observations, the possible environmental quenching process for these massive galaxies might be mergers.

The equilibrium configuration of a solid strange star in the final inspiral phase with another compact object is generally discussed, and the starquake-related issue is revisited, for a special purpose to understand the precursor emission of binary compact star merger events (e.g., that of GRB211211A). As the binary system inspirals inward due to gravitational wave radiation, the ellipticity of the solid strangeon star increases due to the growing tidal field of its compact companion. Elastic energy is hence accumulated during the inspiral stage which might trigger a starquake before the merger when exceeds a critical value. The energy released during such starquakes is calculated and compared to the precursor observation of GRB211211A. The result shows that the energy might be insufficient for binary strangeon-star case unless the entire solid strangeon star shatters, and hence favors a black hole-strangeon star scenario for GRB211211A. The timescale of the precursor as well as the frequency of the observed quasi-periodic-oscillation have also been discussed in the starquake model.

Well-studied very metal-poor (VMP, [Fe/H] < -2 ) stars in the inner Galaxy are few in number, and they are of special interest because they are expected to be among the oldest stars in the MilkyWay. We present high-resolution spectroscopic follow-up of the carbon-enhanced metal-poor (CEMP) star Pristine_184237.56-260624.5 (hereafter Pr184237) identified in the Pristine Inner Galaxy Survey. This star has an apocentre of about 2 kpc. Its atmospheric parameters (Teff = 5100 K, log g = 2.0, [Fe/H] = -2.60) were derived based on the non-local thermodynamic equilibrium (NLTE) line formation. We determined abundances for 32 elements, including 15 heavy elements beyond the iron group. The NLTE abundances were calculated for 13 elements from Na to Pb. Pr184237 is strongly enhanced in C, N, O, and both s- and r-process elements from Ba to Pb; it reveals a low carbon isotope ratio of 12C/13C = 7. The element abundance pattern in the Na-Zn range is typical of halo stars. With [Ba/Eu] = 0.32, Pr184237 is the first star of the CEMP-r/s subclass identified in the inner Galaxy. Variations in radial velocity suggest binarity. We tested whether a pollution by the s- or i-process material produced in the more massive and evolved companion can form the observed abundance pattern and find that an i-process in the asymptotic giant branch star with a progenitor mass of 1.0-2.0 Msun can be the solution.

Very high energy (VHE; 100 GeV $<$ E $\leq$ 100 TeV) and high energy (HE; 100 MeV $<$ E $\leq$ 100 GeV) gamma-rays were observed from the symbiotic recurrent nova RS Ophiuchi (RS Oph) during its outburst in August 2021, by various observatories such as High Energy Stereoscopic System (H.E.S.S.), Major Atmospheric Gamma Imaging Cherenkov (MAGIC), and {\it Fermi}-Large Area Telescope (LAT). The models explored so far tend to favor a hadronic scenario of particle acceleration over an alternative leptonic scenario. This paper explores a time-dependent lepto-hadronic scenario to explain the emission from the RS Oph source region. We have used simultaneous low frequency radio data observed by various observatories, along with the data provided by H.E.S.S., MAGIC, and \textit{Fermi}-LAT, to explain the multi-wavelength (MWL) spectral energy distributions (SEDs) corresponding to 4 days after the outburst. Our results show that a lepto-hadronic interpretation of the source not only explains the observed HE-VHE gamma-ray data but the corresponding model synchrotron component is also consistent with the first 4 days of low radio frequency data, indicating the presence of non-thermal radio emission at the initial stage of nova outburst. We have also calculated the expected neutrino flux from the source region and discussed the possibility of detecting neutrinos.

SMSS\,J114447.77-430859.3 ($z=0.83$) has been identified in the SkyMapper Southern Survey as the most luminous quasar in the last $\sim 9\,\rm Gyr$. In this paper, we report on the eROSITA/Spectrum-Roentgen-Gamma (SRG) observations of the source from the eROSITA All Sky Survey, along with presenting results from recent monitoring performed using Swift, XMM-Newton, and NuSTAR. The source shows a clear variability by factors of $\sim 10$ and $\sim 2.7$ over timescales of a year and of a few days, respectively. When fit with an absorbed power law plus high-energy cutoff, the X-ray spectra reveal a $\Gamma=2.2 \pm 0.2$ and $E_{\rm cut}=23^{+26}_{-5}\,\rm keV$. Assuming Comptonisation, we estimate a coronal optical depth and electron temperature of $\tau=2.5-5.3\, (5.2-8)$ and $kT=8-18\, (7.5-14)\,\rm keV$, respectively, for a slab (spherical) geometry. The broadband SED is successfully modelled by assuming either a standard accretion disc illuminated by a central X-ray source, or a thin disc with a slim disc emissivity profile. The former model results in a black hole mass estimate of the order of $10^{10}\,M_\odot$, slightly higher than prior optical estimates; meanwhile, the latter model suggests a lower mass. Both models suggest sub-Eddington accretion when assuming a spinning black hole, and a compact ($\sim 10\,r_{\rm g}$) X-ray corona. The measured intrinsic column density and the Eddington ratio strongly suggest the presence of an outflow driven by radiation pressure. This is also supported by variation of absorption by an order of magnitude over the period of $\sim 900\,\rm days$.

We study the latitudinal distribution and temporal evolution of the sunspot penumbra-umbra ratio (q) for the even and odd Solar Cycles 12-24 of RGO sunspot groups, SC21-SC24 of Debrecen sunspot groups and Kodaikanal sunspot dataset for SC16-SC24. We find that RGO even (odd) Cycles have q-values 5.20 (4.75), Kodaikanal even (odd) cycles have q-values 5.27 (5.43), and Debrecen cycles has q-value 5.74 on the average. We also show that q is at lowest around the Equator of the Sun and increases towards higher latitudes having maximum values at about 10-25 degrees. This is understandable, because smaller sunspots and groups locate nearer to Equator and have smaller q-values than larger sunspots and groups, which maximize at about 10-20 degrees at both hemispheres. The error limits are very wide and thus the confidence of this result is somewhat vague. For Debrecen dataset we find a deep valley in the temporal q-values before the middle of the cycle. We show that this exists simultaneously with the Gnevyshev gap (GG) in the graph of the total and umbral areas of the large sunspot groups. Other databases do not show GG in their q-graphs, although GG exists in their temporal total area and umbral area.

Accreting low mass X-ray binaries (LMXBs) are capable of launching powerful outflows such as accretion disc winds. In disc winds, vast amounts of material can be carried away, potentially greatly impacting the binary and its environment. Previous studies have uncovered signatures of disc winds in the X-ray, optical, near-infrared, and recently even the UV band, predominantly in LMXBs with large discs ($P_{orb}{\geq}20$ hrs). Here, we present the discovery of transient UV outflow features in UW CrB, a high-inclination ($i{\geq}77$\deg) neutron star LMXB with an orbital period of only $P_{orb}{\approx}111$ min. We present P-Cygni profiles detected for Si iv 1400\r{A} and tentatively for N v 1240\r{A} in one 15 min exposure, which is the only exposure covering orbital phase $\phi{\approx}0.7{-}0.8$, with a velocity of ${\approx}1500$ km/s. We show that due to the presence of black body emission from the neutron star surface and/or boundary layer, a thermal disc wind can be driven despite the short $P_{orb}$, but explore alternative scenarios as well. The discovery that thermal disc winds may occur in NS-LMXBs with $P_{orb}$ as small as ${\approx}111$ min, and can potentially be transient on time scales as short as ${\approx}15$ min, warrants further observational and theoretical work.

Using the HI self-absorption data from the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we perform a study of the cold atomic gas in the Cygnus-X North region. The most remarkable HI cloud is characterized by a filamentary structure, associated in space and in velocity with the principle molecular filament in the Cygnus-X North region. We investigate the transition from the atomic filament to the molecular filament. We find that the HII regions Cygnus OB2 and G081.920+00.138 play a critical role in compressing and shaping the atomic Cygnus-X North filament, where the molecular filament subsequently forms. The cold HI in the DR21 filament has a much larger column density (N(HI) $\sim$ 1 $\times$ 10$^{20}$ cm$^{-2}$) than the theoretical value of the residual atomic gas ($\sim$ 1 $\times$ 10$^{19}$ cm$^{-2}$), suggesting that the HI-to-H$_2$ transition is still in progress. The timescale of the HI-to-H$_2$ transition is estimated to be 3 $\times$ 10$^{5}$ yr, which approximates the ages of massive protostars in the Cygnus-X North region. This implies that the formation of molecular clouds and massive stars may occur almost simultaneously in the DR21 filament, in accord with a picture of rapid and dynamic cloud evolution.

We describe the optical alignment method for the Prime-focus Infrared Microlensing Experiment (PRIME) telescope which is a prime-focus near-infrared (NIR) telescope with a wide field of view for the microlensing planet survey toward the Galactic center that is the major task for the PRIME project. There are three steps for the optical alignment: preliminary alignment by a laser tracker, fine alignment by intra- and extra-focal (IFEF) image analysis technique, and complementary and fine alignment by the Hartmann test. We demonstrated that the first two steps work well by the test conducted in the laboratory in Japan. The telescope was installed at the Sutherland Observatory of South African Astronomical Observatory in August, 2022. At the final stage of the installation, we demonstrated that the third method works well and the optical system satisfies the operational requirement.

Gaia Data Release 3 contains accurate photometric observations of more than 150,000 asteroids covering a time interval of 34 months. With a total of about 3,000,000 measurements, a typical number of observations per asteroid ranges from a few to several tens. We aimed to reconstruct the spin states and shapes of asteroids from this dataset. We computed the viewing and illumination geometry for each individual observation and used the light curve inversion method to find the best-fit asteroid model, which was parameterized by the sidereal rotation period, the spin axis direction, and a low-resolution convex shape. To find the best-fit model, we ran the inversion for tens of thousands of trial periods on interval 2-10,000 h, with tens of initial pole directions. To find the correct rotation period, we also used a triaxial ellipsoid model for the shape approximation. In most cases the number of data points was insufficient to uniquely determine the rotation period. However, for about 8600 asteroids we were able to determine the spin state uniquely together with a low-resolution convex shape model. This large sample of new asteroid models enables us to study the spin distribution in the asteroid population. The distribution of spins confirms previous findings that (i) small asteroids have poles clustered toward ecliptic poles, likely because of the YORP-induced spin evolution, (ii) asteroid migration due to the Yarkovsky effect depends on the spin orientation, and (iii) members of asteroid families have the sense of rotation correlated with their proper semimajor axis: over the age of the family, orbits of prograde rotators evolved, due to the Yarkovsky effect, to larger semimajor axes, while those of retrograde rotators drifted in the opposite direction.

We present the high-precision optical polarimetric observations of black hole X-ray binary Cyg X-1, spanning several cycles of its 5.6 day orbital period. Week-long observations on two telescopes located in opposite hemispheres allowed us to track the evolution of the polarization within one orbital cycle with the highest temporal resolution to date. Using the field stars, we determine the interstellar polarization in the source direction and subsequently its intrinsic polarization. The optical polarization angle is aligned with that in the X-rays as recently obtained with the Imaging X-ray Polarimetry Explorer. Furthermore, it is consistent, within the uncertainties, with the position angle of the radio ejections. We show that the intrinsic PD is variable with the orbital period with the amplitude of $\sim$0.2% and discuss various sites of its production. Assuming the polarization arises from a single Thomson scattering of the primary star radiation by the matter that follows the black hole in its orbital motion, we constrain the inclination of the binary orbit $i>120^\circ$ and its eccentricity $e<0.08$. The asymmetric shape of the orbital profiles of Stokes parameters implies also the asymmetry of the scattering matter distribution about the orbital plane, which may arise from the tilted accretion disk. We compare our data to the polarimetric observations made over 1975-1987 and find good, within $1^\circ$, agreement between the intrinsic polarization angles. On the other hand, the PD decreased by 0.4% over half a century, suggesting the presence of secular changes in the geometry of accreting matter.

The intracluster light (ICL) fraction, measured at certain specific wavelengths, has been shown to provide a good marker for determining the dynamical stage of galaxy clusters, i.e., merging versus relaxed, for small to intermediate redshifts. Here, we apply it for the first time to a high-redshift system, SPT-CLJ0615-5746 at z=0.97, using its RELICS (Reionization Lensing Cluster Survey) observations in the optical and infrared. We find the ICL fraction signature of merging, with values ranging from 16 to 37%. A careful re-analysis of the X-ray data available for this cluster points to the presence of at least one current merger, and plausibly a second merger. These two results are in contradiction with previous works based on X-ray data, which claimed the relaxed state of SPT-CLJ0615-5746, and confirmed the evidences presented by kinematic analyses. We also found an abnormally high ICL fraction in the rest-frame near ultraviolet wavelengths, which may be attributed to the combination of several phenomena such as an ICL injection during recent mergers of stars with average early-type spectra, the reversed star formation-density relation found at this high redshift in comparison with lower-redshift clusters, and projection effects.

Characterizing the structural properties of galaxies in high-redshift protoclusters is key to our understanding of the environmental effects on galaxy evolution in the early stages of galaxy and structure formation. In this study, we assess the structural properties of 85 and 87 Halpha emission-line candidates (HAEs) in the densest regions of two massive protoclusters, BOSS1244 and BOSS1542, respectively, using HST H-band imaging data. Our results show a true pair fraction of 22+-5 (33+-6) percent in BOSS1244 (BOSS1542), which yields a merger rate of 0.41+-0.09 (0.52+-0.04) per Gyr for massive HAEs with log (M_*/M_sun) > 10.3. This rate is 1.8 (2.8) times higher than that of the general fields at the same epoch. Our sample of HAEs exhibits half-light radii and Sersic indices that cover a broader range than field star-forming galaxies. Additionally, about 15 percent of the HAEs are as compact as the most massive (log(M_*/M_sun) > 11) spheroid-dominated population. These results suggest that the high galaxy density and cold dynamical state (i.e., velocity dispersion of <400 km/s) are key factors that drive galaxy mergers and promote structural evolution in the two protoclusters. Our findings also indicate that both the local environment (on group scales) and the global environment play essential roles in shaping galaxy morphologies in protoclusters. This is evident in the systematic differences observed in the structural properties of galaxies between BOSS1244 and BOSS1542.

We study the morphological properties of mid-infrared selected galaxies at $1.0<z<1.7$ in the SMACS J0723.3-7327 cluster field, to investigate the mechanisms of galaxy mass assembly and structural formation at cosmic noon. We develop a new algorithm to decompose the dust and stellar components of individual galaxies by utilizing high-resolution images in the MIRI F770W and NIRCam F200W bands. Our analyses reveal that most galaxies in the stellar mass range ${\rm 10^{9.5}<M_*/M_\odot<10^{10.5}}$ have dust cores relatively compact compared to their stellar cores, whereas the most massive ($\rm{M_* \sim 10^{10.9}\,M_\odot}$) galaxy in our sample displays a comparably compact stellar core as to dust. The observed compactness of the dust component is potentially attributed to the presence of a (rapidly growing) massive bulge, in some cases associated with elevated star formation. Expanding the sample size through a joint analysis of multiple Cycle~1 deep-imaging programs can help to confirm the inferred picture. Our pilot study highlights that MIRI offers an efficient approach to studying the structural formation of galaxies from cosmic noon to the modern universe.

The Large Magellanic Cloud (LMC) has a complex dynamics driven by both internal and external processes. The external forces are due to tidal interactions with the Small Magellanic Cloud and the Milky Way, while internally its dynamics mainly depends on the stellar, gas, and dark matter mass distributions. Despite the overall complexity of the system, very often simple physical models can give us important insights about the main driving factors. Here we focus on the internal forces and attempt to model the proper motions of $\sim10^6$ stars in the LMC as measured by Gaia Data Release 3 with an axisymmetric dynamical model, based on the Jeans equations. We test both cored and cusped spherical Navarro-Frenk-White dark matter halos to fit the LMC gravitational potential. We find that this simple model is very successful at selecting a clean sample of genuine LMC member stars and predicts the geometry and orientation of the LMC with respect to the observer within the constraint of axisymmetry. Our Jeans dynamical models describe well the rotation profile and the velocity dispersion of the LMC stellar disc, however they fail to describe the motions of the LMC bar, which is a non-axisymmetric feature dominating the central region. We plan a triaxial Schwarzschild approach as a next step for the dynamical modelling of the LMC.

We present a new nucleosynthesis process that may take place on neutron-rich ejecta experiencing an intensive neutrino flux. The nucleosynthesis proceeds similarly to the standard $r$-process, a sequence of neutron-captures and beta-decays, however with charged-current neutrino absorption reactions on nuclei operating much faster than beta-decays. Once neutron capture reactions freeze-out the produced $r$-process neutron-rich nuclei undergo a fast conversion of neutrons into protons and are pushed even beyond the $\beta$-stability line producing the neutron-deficient $p$-nuclei. This scenario, which we denote as the $\nu r$-process, provides an alternative channel for the production of $p$-nuclei and the short-lived nucleus $^{92}$Nb. We discuss the necessary conditions posed on the astrophysical site for the $\nu r$-process to be realized in nature. While these conditions are not fulfilled by current neutrino-hydrodynamic models of $r$-process sites, future models, including more complex physics and a larger variety of outflow conditions, may achieve the necessary conditions in some regions of the ejecta.

Observations from the NASA Solar Dynamic Observatory Atmospheric Imaging Assembly were employed to investigate targeted physical properties of coronal active region structures across the entirety of Solar Cycle 24 (dates). This is the largest consistent study to date which analyses emergent trends in structural width, location, and occurrence rate by performing an automatic and long-term examination of observable coronal limb features within equatorial active region belts across four extreme ultraviolet wavelengths (171, 193, 211, and 304 angstroms). This has resulted in over thirty thousand observed coronal structures and hence allows for the production of spatial and temporal distributions focused upon the rise, peak and decay activity phases of Solar Cycle 24. Employing a self-organized-criticality approach as a descriptor of coronal structure formation, power law slopes of structural widths versus frequency are determined, ranging from -1.6 to -3.3 with variations of up to 0.7 found between differing periods of the solar cycle, compared to a predicted Fractal Diffusive Self Organized Criticality (FD-SOC) value of -1.5. The North-South hemispheric asymmetry of these structures was also examined with the northern hemisphere exhibiting activity that is peaking earlier and decaying slower than the southern hemisphere, with a characteristic "butterfly" pattern of coronal structures detected. This represents the first survey of coronal structures performed across an entire solar cycle, demonstrating new techniques available to examine the composition of the corona by latitude in varying wavelengths at selected altitudes.

Individual vibrational band spectroscopy presents an opportunity to examine exoplanet atmospheres in detail by distinguishing where the vibrational state populations of molecules differ from the current assumption of a Boltzmann distribution. Here, retrieving vibrational bands of OH in exoplanet atmospheres is explored using the hot Jupiter WASP-33b as an example. We simulate low-resolution spectroscopic data for observations with the JWST's NIRSpec instrument and use high resolution observational data obtained from the Subaru InfraRed Doppler instrument (IRD). Vibrational band-specific OH cross section sets are constructed and used in retrievals on the (simulated) low and (real) high resolution data. Low resolution observations are simulated for two WASP-33b emission scenarios: under the assumption of local thermal equilibrium (LTE) and a toy non-LTE model for vibrational excitation of selected bands. We show that mixing ratios for individual bands can be retrieved with sufficient precision to allow the vibrational population distributions of the forward models to be reconstructed. A simple fit for the Boltzmann distribution in the LTE case shows that the vibrational temperature is recoverable in this manner. For high resolution, cross-correlation applications, we apply the individual vibrational band analysis to an IRD spectrum of WASP-33b, applying an 'un-peeling' technique. Individual detection significances for the two strongest bands are shown to be in line with Boltzmann distributed vibrational state populations consistent with the effective temperature of the WASP-33b atmosphere reported previously. We show the viability of this approach for analysing the individual vibrational state populations behind observed and simulated spectra including reconstructing state population distributions.

The volume of data from current and future observatories has motivated the increased development and application of automated machine learning methodologies for astronomy. However, less attention has been given to the production of standardised datasets for assessing the performance of different machine learning algorithms within astronomy and astrophysics. Here we describe in detail the MiraBest dataset, a publicly available batched dataset of 1256 radio-loud AGN from NVSS and FIRST, filtered to $0.03 < z < 0.1$, manually labelled by Miraghaei and Best (2017) according to the Fanaroff-Riley morphological classification, created for machine learning applications and compatible for use with standard deep learning libraries. We outline the principles underlying the construction of the dataset, the sample selection and pre-processing methodology, dataset structure and composition, as well as a comparison of MiraBest to other datasets used in the literature. Existing applications that utilise the MiraBest dataset are reviewed, and an extended dataset of 2100 sources is created by cross-matching MiraBest with other catalogues of radio-loud AGN that have been used more widely in the literature for machine learning applications.

We develop a tool, which we name Protoplanetary Disk Operator Network (PPDONet), that can predict the solution of disk-planet interactions in protoplanetary disks in real-time. We base our tool on Deep Operator Networks (DeepONets), a class of neural networks capable of learning non-linear operators to represent deterministic and stochastic differential equations. With PPDONet we map three scalar parameters in a disk-planet system -- the Shakura \& Sunyaev viscosity $\alpha$, the disk aspect ratio $h_\mathrm{0}$, and the planet-star mass ratio $q$ -- to steady-state solutions of the disk surface density, radial velocity, and azimuthal velocity. We demonstrate the accuracy of the PPDONet solutions using a comprehensive set of tests. Our tool is able to predict the outcome of disk-planet interaction for one system in less than a second on a laptop. A public implementation of PPDONet is available at \url{https://github.com/smao-astro/PPDONet}.

The apparent brightness of satellites is calculated as a function of satellite position as seen by a ground-based observer in darkness. Both direct illumination of the satellite by the sun as well as indirect illumination due to reflection from the Earth are included. The reflecting properties of each satellite component and of the Earth must first be estimated (the Bidirectional Reflectance Distribution Function, BRDF). Integrating over all scattering surfaces leads to the angular pattern of the flux reflected from the satellite. Finally, the apparent brightness of the satellite as seen by an observer at a given location is calculated as a function of satellite position. We validate our calculations by comparing to observations of selected Starlink satellites and show significant improvement on previous satellite brightness models. With multiple observations of a satellite at various solar angles and with minimal assumptions regarding the satellite, BRDF model coefficients for each satellite component can be accurately inferred, obviating the need to import direct BRDF lab measurements. This widens the effectiveness of this model approach to virtually all satellites. This work finds application in satellite design and operations, and in planning observatory data acquisition and analysis. Similar methodology for predicting satellite brightness has already informed mitigation strategies for next generation Starlink satellites.

We present detailed optical photometry and spectroscopy of the Type IIn supernova (SN) 2021qqp. Its unusual light curve is marked by a long gradual brightening (i.e., precursor) for about 300 days, a rapid increase in brightness for about 60 days, and then a sharp increase of about 1.6 mag in only a few days to a first peak of $M_r\approx -19.5$ mag. The light curve then turns over and declines rapidly, until it re-brightens to a second distinct and sharp peak with $M_r\approx -17.3$ mag centered at about 335 days after the first peak. The spectra are dominated by Balmer-series lines with a complex morphology that includes a narrow component with a width of $\approx 1300$ km s$^{-1}$ (first peak) and $\approx 2500$ km s$^{-1}$ (second peak) that we associate with the circumstellar medium (CSM), and a P Cygni component with an absorption velocity of $\approx 8500$ km s$^{-1}$ (first peak) and $\approx 5600$ km s$^{-1}$ (second peak) that we associate with the SN-CSM interaction shell. Using the bolometric light curve and velocity evolution, we construct an analytical model to extract the CSM profile and SN properties. We find two significant mass-loss episodes with peak mass loss rates of $\approx 10$ M$_\odot$ yr$^{-1}$ and $\approx 5$ M$_\odot$ yr$^{-1}$ about 0.8 and 2 years before explosion, and a total CSM mass of $\approx 2-4\,M_\odot$. We show that the most recent mass-loss episode can explain the precursor for the year preceding the explosion. The SN ejecta mass is constrained to be $M_{\rm SN}\approx 5-30\,M_\odot$ for an explosion energy of $E_{\rm SN}\approx (3-10)\times10^{51}\,{\rm erg}$. We discuss eruptive massive stars (luminous blue variable, pulsational pair instability) and an extreme stellar merger with a compact object as possible progenitor channels for generating the energetic explosion in the complex CSM environment.

A gravitational machine is defined as an arrangement of gravitating masses from which useful energy can be extracted. It is shown that such machines may exist if the masses are of normal astronomical size. A simple example of a gravitational machine, consisting of a double star with smaller masses orbiting around it, is described. It is shown that an efficient gravitational machine will also be an emitter of gravitational radiation. The emitted radiation sets a limit on the possible performance of gravitational machines, and also provides us with a possible means for detecting such machines if they exist.

Gravitational parity violation arises in a variety of theories beyond general relativity. Gravitational waves in such theories have their propagation altered, leading to birefringence effects in both the amplitude and speed of the wave. In this work, we introduce a generalized, theory-motivated parametrization scheme to study parity violation in gravitational wave propagation. This parametrization maps to parity-violating gravity theories in a straightforward way. We find that the amplitude and velocity birefringence effects scale with an effective distance measure that depends on how the dispersion relation is modified. Furthermore, we show that this generic parametrization can be mapped to the parametrized-post-Einsteinian (ppE) formalism with convenient applications to gravitational wave observations and model-agnostic tests of general relativity. We derive a mapping to the standard ppE waveform of the gravitational wave response function, and also find a ppE waveform mapping at the level of the polarization modes, $h_+$ and $h_\times$. Finally, we show how existing constraints in the literature translate to bounds on our new parity-violating parameters and discuss avenues for future analysis.

We use numerical simulations of scalar field dark matter evolving on a moving black hole background to confirm the regime of validity of (semi-)analytic expressions derived from first principles for both dynamical friction and momentum accretion in the relativistic regime. We cover both small and large clouds (relative to the de Broglie wavelength of the scalars), and light and heavy particle masses (relative to the BH size). In the case of a small dark matter cloud, the effect of accretion is a non-negligible contribution to the total force on the black hole, even for small scalar masses. We confirm that this momentum accretion transitions between two regimes (wave- and particle-like) and we identify the mass of the scalar at which the transition between regimes occurs.

In this work, we complete our CT18qed study with the neutron's photon parton distribution function (PDF), which is essential for the nucleus scattering phenomenology. Two methods, CT18lux and CT18qed, based on the LUXqed formalism and the DGLAP evolution, respectively, to determine the neutron's photon PDF have been presented. Various low-$Q^2$ non-perturbative variations have been carefully examined, which are treated as additional uncertainties on top of those induced by quark and gluon PDFs. The impacts of the momentum sum rule as well as isospin symmetry violation have been explored, and turn out to be negligible. A detailed comparison with other neutron's photon PDF sets has been performed, which shows a great improvement in the precision and a reasonable uncertainty estimation in our results. Finally, two phenomenological implications are demonstrated with photon-initiated processes: neutrino-nucleus $W$-boson production, which is important for the near-future TeV--PeV neutrino observations, and the axion-like particle production at a high-energy muon beam-dump experiment.

Post-inflationary reheating phase is usually said to be solely governed by the decay of coherently oscillating inflaton into radiation. In this submission, we explore a new avenue toward reheating through the evaporation of primordial black holes (PBHs). After the inflation, if PBHs form, depending on its initial mass, abundance, and inflaton coupling with the radiation, we found two physically distinct possibilities of reheating the universe. In one possibility, the thermal bath is solely obtained from the decay of PBHs while inflaton plays the role of dominant energy component in the entire process. In the other possibility, we found that PBHs itself dominate the total energy budget of the Universe during the course of evolution, and then its subsequent evaporation leads to radiation dominated universe. Furthermore, we analyze the impact of both monochromatic and extended PBH mass functions and estimate the detailed parameter ranges for which those distinct reheating histories are realized.

General relativity is one of the pillars of modern physics. For decades, the theory has been mainly tested in the weak field regime with experiments in the Solar System and radio observations of binary pulsars. Until 2015, the strong field regime was almost completely unexplored. Thanks to new observational facilities, the situation has dramatically changed in the last few years. Today we have gravitational wave data of the coalesce of stellar-mass compact objects from the LIGO-Virgo-KAGRA Collaboration, images at mm wavelengths of the supermassive black holes in M87$^*$ and SgrA$^*$ from the Event Horizon Telescope Collaboration, and X-ray data of accreting compact objects from a number of X-ray missions. Gravitational wave tests and black hole imaging tests are certainly more popular and are discussed in other articles of this Special Issue. The aim of the present manuscript is to provide a pedagogical review on X-ray tests of general relativity with black holes and to compare this kind of tests with those possible with gravitational wave data and black hole imaging.

The relativistic spin-precession equations for black-hole binaries have four different equilibrium solutions that correspond to systems where the two individual black hole spins are either aligned or anti-aligned with the orbital angular momentum. Surprisingly, it was demonstrated that only three of these equilibrium solutions are stable. Binary systems in the up-down configuration, where the spin of the heavier (lighter) black hole is co- (counter-) aligned with the orbital angular momentum, might be unstable to small perturbations of the spin directions. After the onset of the up-down instability, that occurs after a specific critical orbital separation $r_\mathrm{UD+}$, the binary becomes unstable to spin precession leading to large misalignment of the spins. In this work, we present a Bayesian procedure based on the Savage-Dickey density ratio to test the up-down origin of gravitational-wave events. We apply this procedure to look for promising candidates among the events detected so far during the first three observing runs performed by LIGO/Virgo.

Warm inflation has normalized two ideas in cosmology, that in the early universe the initial primordial density perturbations generally could be of classical rather than quantum origin and that during inflation, particle production from interactions amongst quantum field, and its backreaction effects, can occur concurrent with inflationary expansion. When we first introduced these ideas, both were met with resistance, but today they are widely accepted as possibilities with many models and applications based on them, which is an indication of the widespread influence of warm inflation. Open quantum field theory, which has been utilized in studies of warm inflation, is by now a relevant subject in cosmology, in part due to this early work. In this review I first discuss the basic warm inflation dynamics. I then outline how to compute warm inflation dynamics from first principles quantum field theory (QFT) and in particular how a dissipative term arises. Warm inflation models can have an inflaton mass bigger than the Hubble scale and the inflaton field excursion can remain sub-Planckian, thus overcoming the most prohibitive problems of inflation model building. I discuss the early period of my work in developing warm inflation that helped me arrive at these important features of its dynamics. Inflationary cosmology today is immersed in hypothetical models, which by now are acting as a diversion from reaching any endgame in this field. I discuss better ways to approach model selection and give necessary requirements for a well constrained and predictive inflation model. I point out a few warm inflation models that could be developed to this extent. I discuss how at this stage more progress would be made in this subject by taking a broader view on the possible early universe solutions that include not just inflation but the diverse range of options.

The violation of the null energy condition (NEC) is closely related to potential solutions for the cosmological singularity problem and may therefore play a crucial role in the very early universe. We explore a novel approach to generate primordial black holes (PBHs) via the violation of the NEC in a single-field inflationary scenario. In our scenario, the universe transitions from a first slow-roll inflation stage with a Hubble parameter H = Hinf1 to a second slow-roll inflation stage with H = Hinf2 > Hinf1, passing through an intermediate stage of NEC violation. The resulting primordial scalar power spectrum is naturally enhanced by the NEC violation at a certain wavelength. As a result, PBHs with masses and abundances of observational interest can be produced in our scenario. We also examine the phenomenological signatures of scalar-induced gravitational waves (SIGWs). Our work highlights the significance of utilizing a combination of PBHs, SIGWs, and primordial gravitational waves as a powerful probe for exploring the NEC violation during inflation.

Single-photon avalanche detectors (SPADs) are well-suited for satellite-based quantum communication because of their advantageous operating characteristics as well as their relatively straightforward and robust integration into satellite payloads. However, space-borne SPADs will encounter damage from space radiation, which usually manifests itself in the form of elevated dark counts. Methods for mitigating this radiation damage have been previously explored, such as thermal and optical (laser) annealing. Here we investigate in a lab, using a CubeSat payload, laser annealing protocols in terms of annealing laser power and annealing duration, for their possible later use in orbit. Four Si SPADs (Excelitas SLiK) irradiated to an equivalent of 10 years in low Earth orbit exhibit very high dark count rates (>300 kcps at -22 C operating temperature) and significant saturation effects. We show that annealing them with optical power between 1 and 2 W yields reduction in dark count rate by a factor of up to 48, as well as regaining SPAD sensitivity to a very faint optical signal (on the order of single photon) and alleviation of saturation effects. Our results suggest that an annealing duration as short as 10 seconds can reduce dark counts, which can be beneficial for power-limited small-satellite quantum communication missions. Overall, annealing power appears to be more critical than annealing duration and number of annealing exposures.

Background: The equations of state (EoSs) which determine the properties of neutron stars (NSs) are often characterized by the iso-scalar and iso-vector nuclear matter parameters (NMPs). Recent attempts to relate the radius and tidal deformability of a NS to the individual NMPs have been inconclusive. These properties display strong correlations with the pressure of NS matter which depends on several NMPs. The knowledge of minimal NMPs that determine the NS properties will be necessary to address any connection between NS properties (e.g., tidal deformability) and that of finite nuclei. Purpose: To identify the important NMPs required to describe the tidal deformability of neutron star for astrophysically relevant range of their gravitational masses (1.2 -- 1.8 M$_\odot$) as encountered in the binary neutron star merger events. Method: We construct a large set of EoSs using four iso-scalar and five iso-vector NMPs. These EOSs are employed to perform a systematic analysis to isolate the NMPs that predominantly determine the tidal deformability, over a wide range of NS mass. The tidal deformability is then directly mapped to these NMPs. Results: The tidal deformability of the NS with mass 1.2-1.8 M$_\odot$ can be determined within 10$\%$ directly in terms of four nuclear matter parameters, namely, the incompressibility $K_0$ and skewness $Q_0$ of symmetric nuclear matter, and the slope $L_0$ and curvature parameter $K_{\rm sym,0}$ of symmetry energy. Conclusion: A function that quickly estimates the value of tidal deformability in terms of minimal nuclear matter parameters is developed. Our method can also be extended to other NS observables.

We propose a simple model in the type-III seesaw framework to explain the recently reported W-mass anomaly by CDF-II collaboration, neutrino mass, asymmetric dark matter, and baryon asymmetry of the Universe. We extend the standard model with a vector-like singlet lepton ($\chi$) and a hypercharge zero scalar triplet ($\Delta$) in addition to three hypercharge zero triplet fermions($\Sigma_i~,i=1,2,3$). A $Z_2$ symmetry is imposed under which $\chi$ and $\Delta$ are odd, while all other particles are even. As a result, the lightest $Z_2$ odd particle $\chi$ behaves as a candidate of dark matter. In the early Universe, the CP-violating out-of-equilibrium decay of heavy triplet fermions to the Standard Model lepton ($L$) and Higgs ($H$) generate a net lepton asymmetry, while that of triplet fermions to $\chi$ and $\Delta$ generate a net asymmetric dark matter. The lepton asymmetry is converted to the required baryon asymmetry of the Universe via the electroweak sphalerons, while the asymmetry in $\chi$ remains as a dark matter relic that we observe today. We introduce a singlet scalar $\phi$, with mass $m_\phi < m_\chi$, which not only assists to deplete the symmetric component of $\chi$ through the annihilation process: $\bar{\chi} \chi \to \phi \phi$ but also paves a path to detect dark matter $\chi$ at direct search experiments through $\phi-H$ mixing. The $Z_2$ symmetry is broken softly resulting in an unstable asymmetric dark matter with mass ranging from a few MeV to a few tens of GeV. The softly broken $Z_2$ symmetry also induces a vacuum expectation value (vev) of $\Delta$ due to which the asymmetry in $\Delta$ disappears. Moreover, the vev of $\Delta$ enhances the W-boson mass as reported by CDF-II collaboration with $7\sigma$ statistical significance, while keeping the $Z$-boson mass intact.