Papers for Wednesday, Oct 27 2021

The emergence of the first, so-called Population III (Pop III), stars shaped early cosmic history in ways that crucially depends on their initial mass function (IMF). However, because of the absence of direct observational constraints, the detailed IMF remains elusive. Nevertheless, numerical simulations agree in broad terms that the first stars were typically massive and should often end their lives in violent, explosive deaths. These fates include extremely luminous pair-instability supernovae (PISNe) and bright gamma-ray bursts (GRBs), the latter arising from the collapse of rapidly rotating progenitor stars into black holes. These high-redshift transients are expected to be within the detection limits of upcoming space telescope missions, allowing to place effective constraints on the shape of the primordial IMF that is not easily accessible with other probes. This paper presents a framework to probe the Pop III IMF, utilizing the cosmological source densities of high-redshift PISNe and GRBs. Considering these transients separately could provide useful constraints on the Pop III IMF, but tighter bounds are obtainable by combining PISN and GRB counts. This combined diagnostic is more robust as it is independent of the underlying Pop III star formation rate density, an unknown prior. Future surveys promise to capture most high-redshift GRBs across the entire sky, but high-redshift PISN searches with future telescopes, e.g. Roman Space Telescope, will likely be substantially incomplete. Nevertheless, we demonstrate that even such lower bounds on the PISN count will be able to provide key constraints on the primordial IMF, in particular, if it is top-heavy or not.

TESS will target the ecliptic plane in Sectors 42 - 46. These sectors overlap with campaigns from the K2 mission, providing a unique opportunity for multi-mission light curve analysis. This data release note describes the combined analysis of K2 and TESS light curves as part of the Quick-Look Pipeline (QLP) procedure, which processes light curves for all targets in TESS Full-Frame Images (FFIs) down to TESS magnitude T = 13.5. We describe updates to our codebase, and the planet transit search, candidate triage, and report generation that are affected by this combined analysis.

The ionizing photon escape fraction (LyC $f_{\rm{esc}}$) of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshift. Here we circumvent this uncertainty by inferring $f_{\rm{esc}}$ with resolved (R>4000) LyA profiles from the X-SHOOTER LyA survey at z=2 (XLS-z2). We select leakers ($f_{\rm{esc}}>20$%) and non-leakers ($f_{\rm{esc}}<5$%) from a representative sample of $>0.2 L^{*}$ LyA emitters (LAEs). With median stacked spectra of these subsets covering 1000-8000 {\AA} (rest-frame) we investigate the conditions for LyC $f_{\rm{esc}}$. We find the following differences between leakers vs. non-leakers: (i) strong nebular CIV and HeII emission vs. non-detections, (ii) O32~8.5 vs. ~3, (iii) Ha/Hb indicating no dust vs. E(B-V)~0.3, (iv) MgII emission close to the systemic velocity vs. redshifted, optically thick MgII, (v) LyA $f_{\rm{esc}}$ of ~50% vs. ~10%. The extreme EWs in leakers (O3+Hb~1100 {\AA}) constrain the characteristic timescale of LyC escape to ~3-10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers -- extremely ionizing stellar populations, low column densities, a dust-free, high ionization state ISM -- occur simultaneously in the $f_{\rm{esc}}>20\%$ stack, suggesting they are causally connected, and motivating why indicators like O32 may suffice to constrain $f_{\rm{esc}}$ at z>6 with JWST. The leakers comprise half our sample, have a median LyC $f_{\rm{esc}}$~50%, and an ionising production efficiency $\log({\xi_{\rm{ion}}/\rm{Hz\ erg^{-1}}})$~25.9. These results show LAEs -- the type of galaxies rare at z=2, but that become the norm at higher redshift -- are highly efficient ionizers, with extreme $\xi_{\rm{ion}}$ and prolific $f_{\rm{esc}}$ occurring in sync. (ABRIDGED)

We study galactic enrichment with rapid neutron capture (r-process) elements in cosmological, magnetohydrodynamical simulations of a Milky Way-mass galaxy. We include a variety of enrichment models, based on either neutron star mergers or a rare class of core-collapse supernova as sole r-process sources. For the first time in cosmological simulations, we implement neutron star natal kicks on-the-fly to study their impact. With kicks, neutron star mergers are more likely to occur outside the galaxy disc, but how far the binaries travel before merging also depends on the kick velocity distribution and shape of the delay time distribution for neutron star mergers. In our fiducial model, the median r-process abundance ratio is somewhat lower and the trend with metallicity is slightly steeper when kicks are included. In a model 'optimized' to better match observations, with a higher rate of early neutron star mergers, the median r-process abundances are fairly unaffected by kicks. In both models, the scatter in r-process abundances is much larger with natal kicks, especially at low metallicity, giving rise to more r-process enhanced stars. We experimented with a range of kick velocities and find that with lower velocities, the scatter is reduced, but still larger than without natal kicks. We discuss the possibility that the observed scatter in r-process abundances is predominantly caused by natal kicks removing the r-process sources far from their birth sites, making enrichment more inhomogeneous, rather than the usual interpretation that the scatter is set by the rarity of its production source.

Dissipative dark matter predicts rich observable phenomena that can be tested with future large-scale structure surveys. As a specific example, we study atomic dark matter, consisting of a heavy particle and a light particle charged under a dark electromagnetism. In particular, we calculate the cosmological evolution of atomic dark matter focusing on dark recombination and dark-molecule formation. We have obtained the relevant interaction-rate coefficients by re-scaling the rates for normal hydrogen, and evolved the abundances for ionized, atomic, and molecular states using a modified version of RecFAST++. We also provide an analytical approximation for the final abundances. We then calculate the effects of the atomic dark matter on the linear power spectrum, which enter through a dark-photon diffusion and dark acoustic oscillations. At the formation time, the atomic dark matter model suppresses halo abundances on scales smaller than the diffusion scale, just like the warm dark matter models suppress the abundance below the free-streaming scale. The subsequent evolution with radiative cooling, however, will alter the halo mass function further.

Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can help understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing observations. To minimise the chances of biases, we focus on galaxy properties that can be determined robustly: the surface brightness profile and the colour. Methods. We selected 2439 central massive elliptical galaxies from the SDSS spectroscopic sample. We first measured their surface brightness profile and colours by fitting Sersic models to photometric data from the Kilo-Degree Survey (KiDS). We fitted their halo mass distribution as a function of redshift, rest-frame $r-$band luminosity, half-light radius and rest-frame $u-g$ colour, using KiDS weak lensing data and a Bayesian hierarchical approach. For the sake of robustness to assumptions on the large-radii behaviour of the surface brightness, we repeated the analysis replacing total luminosity and half-light radius with the luminosity within a 10~kpc aperture, $L_{r,10}$, and the light-weighted surface brightness slope, $\Gamma_{10}$. Results. We did not detect any correlation between halo mass and either half-light radius or colour, at fixed redshift and luminosity. Conclusions. Our results indicate that the average star formation efficiency of massive elliptical galaxies has little dependence on their final size or colour. This suggests that the origin of the diversity in the size and colour distribution of these objects lies with properties other than halo mass.

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$\alpha$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$\alpha$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Ly\alpha}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$\alpha$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$\alpha$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $\xi_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), we find that a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. We propose that rapid flashes of LyC from such rare galaxies produce a "disco" ionizing background that may result in significant fluctuations at $z>5$. We conclude with proposed observational tests to further develop our suggested Ly$\alpha$-anchored formalism.

We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of binary neutron stars in quasi-circular orbit that merge and undergo delayed or prompt collapse to a black hole (BH). The stars are irrotational and modeled using an SLy or an H4 nuclear equation of state. To assess the impact of the initial magnetic field configuration on jet launching, we endow the stars with a purely poloidal magnetic field that is initially unimportant dynamically and is either confined to the stellar interior or extends from the interior into the exterior as in typical pulsars. Consistent with our previous results, we find that only the BH + disk remnants originating from binaries that form hypermassive neutron stars (HMNSs) and undergo delayed collapse can drive magnetically-powered jets. We find that the closer the total mass of the binary is to the threshold value for prompt collapse, the shorter is the time delay between the gravitational wave peak amplitude and jet launching. This time delay also strongly depends on the initial magnetic field configuration. We also find that seed magnetic fields confined to the stellar interior can launch a jet over~$\sim 25\,\rm ms$ later than those with pulsar-like magnetic fields. The lifetime of the jet [$\Delta t\lesssim 150\,\rm ms$] and its outgoing Poynting luminosity [$L_{\rm EM}\sim 10^{52\pm 1}\rm erg/s$] are consistent with typical short gamma-ray burst central engine lifetimes, as well as with the Blandford--Znajek mechanism for launching jets and their associated Poynting luminosities. Our numerical results also suggest that the dynamical ejection of matter can be enhanced by the magnetic field. Therefore, GRMHD studies are required to fully understand kilonova signals from GW170818-like events.

Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields.

Dark matter that is dissipative may cool sufficiently to form compact objects, including black holes. Determining the abundance and mass spectrum of those objects requires an accurate model of the chemistry relevant for the cooling of the dark matter gas. Here we introduce a chemistry tool for dark matter, DarkKROME, an extension of the KROME software package. DarkKROME is designed to include all atomic and molecular processes relevant for dark matter with two unequal-mass fundamental fermions, interacting via a massless-photon mediated $U(1)$ force. We use DarkKROME to perform one-zone collapse simulations and study the evolution of temperature-density phase diagrams for various dark-sector parameters.

We present a search for galaxy-scale strong gravitational lenses in the initial 2 500 square degrees of the Canada-France Imaging Survey (CFIS). We design a convolutional neural network (CNN) committee that we apply on a selection of 2 344 002 exquisite-seeing r-band images of color-selected luminous red galaxies (LRGs). Our training set is particularly realistic, since the deflector and source images of our mock lensing systems are taken from real CFIS r-band and Hubble Space Telescope (HST) images. Only the lensing effect is simulated. A total of 9 460 candidates obtain a score above 0.5 with the CNN committee. After a visual inspection of the candidates, we find a total of 133 lens candidates, among which 104 are completely new. The set of false positives mainly contains ring, spiral and merger galaxies and to a smaller extent galaxies with nearby companions. We classify 32 of the lens candidates as secure lenses and 101 as maybe lenses. For the 32 best-quality lenses, we also fit a singular isothermal ellipsoid mass profile with external shear along with an elliptical Sersic profile for the lens and source light. This modeling step is fully automated and provides distributions of properties for both sources and lenses. We also use auto-encoders to provide a lens/source deblended image of the best lens candidates.

We present the discovery of a new double detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with an P=76.34179(2) min orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits we find a low mass sdB star, $M_{\rm sdB}=0.383\pm0.028$ M$_\odot$ with a massive white dwarf companion, $M_{\rm WD}=0.725\pm0.026$ M$_\odot$. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of $\approx$25 Myrs whereas our MESA model predicts an sdB age of $\approx$170 Myrs. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion $\approx$25 Myrs ago. Using the MESA stellar evolutionary code we find that the sdB star will start mass transfer in $\approx$6 Myrs and in $\approx$60 Myrs the white dwarf will reach a total mass of $0.92$ M$_\odot$ with a thick helium layer of $0.17$ M$_\odot$. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 2238+7430 is only the second confirmed candidate for a double detonation thermonuclear supernova. Using both systems we estimate that at least $\approx$1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.

We compute the isotropic radiation background due to Hawking emission from primordial black holes (PBHs), and examine if this background is a viable option in explaining the excess radiowave background observed by the ARCADE2 and LWA1 experiments at $\lesssim 1\,$GHz. We find that even under the extreme assumption that all of the dark matter is in the form of PBHs, the radio brightness temperature induced by Hawking evaporation of PBHs is $\mathcal{O}(10^{-46})\,$K, highly subdominant compared to the cosmic microwave background. The main reason for this is that for PBHs in the mass range $\sim10^{12}$-$10^{14}\,$kg, which can be constrained by Hawking emission, the spectrum peaks at $10^7$ to $10^5\,$eV. As the Hawking spectrum is power law suppressed towards lower energies, negligible flux of $\mu$eV photons is obtained. The peak of the Hawking spectrum shifts to lower energies for higher masses, but the number density is low and so is the specific intensity. Because Hawking emission from PBHs is thus unable to explain the observed excess radio background, we also consider the alternative possibility of radio emission from gas accretion onto supermassive PBHs. These PBHs can readily produce strong radio emission that could easily explain the ARCADE2/LWA1 excess.

During its departure from Pluto, New Horizons used its LORRI camera to image a portion of Pluto's southern hemisphere that was in a decades-long seasonal winter darkness, but still very faintly illuminated by sunlight reflected by Charon. Recovery of this faint signal was technically challenging. The bright ring of sunlight forward-scattered by haze in the Plutonian atmosphere encircling the nightside hemisphere was severely overexposed, defeating the standard smeared-charge removal required for LORRI images. Reconstruction of the overexposed portions of the raw images, however, allowed adequate corrections to be accomplished. The small solar elongation of Pluto during the departure phase also generated a complex scattered-sunlight background in the images that was three orders of magnitude stronger than the estimated Charon-light flux (the Charon-light flux is similar to the flux of moonlight on Earth a few days before first quarter). A model background image was constructed for each Pluto image based on principal component analysis (PCA) applied to an ensemble of scattered-sunlight images taken at identical Sun-spacecraft geometry to the Pluto images. The recovered Charon-light image revealed a high-albedo region in the southern hemisphere. We argue that this may be a regional deposit of N_2 or CH_4 ice. The Charon-light image also shows that the south polar region currently has markedly lower albedo than the north polar region of Pluto, which may reflect the sublimation of N_2 ice or the deposition of haze particulates during the recent southern summer.

We conducted a deep spectroscopic survey, named SSA22-HIT, in the SSA22 field with DEIMOS on the Keck telescope, designed to tomographically map of high-$z$ HI gas through analysis of Ly$\alpha$ absorption in background galaxies' spectra. In total, 203 galaxies were spectroscopically confirmed at $2.5 < z < 6$ in the $26 \times 15$ arcmin$^2$ area, of which 150 were newly determined in this study. Our redshift measurements were merged with previously confirmed redshifts available in the $34 \times 27$ arcmin$^2$ area of the SSA22 field. This compiled catalog containing $\sim 740$ galaxies of various types at $z > 2$ is useful for various applications, and it is made publicly available. Our SSA22-HIT survey has increased by approximately twice the number of spectroscopic redshifts of sources at $z > 3.2$ in the observed field. This greatly benefits the HI tomographic mapping of the proto-cluster region at $z = 3.1$, which we present in a parallel study. From a comparison with publicly available redshift catalogs, we show that our compiled redshift catalog in the SSA22 field is comparable to those among major extra-galactic survey fields in terms of a combination of wide area and high surface number density of objects at $z > 2$.

The theory of substellar evolution predicts that there is a sharp mass boundary between lithium and non-lithium brown dwarfs, not far below the substellar-mass limit. The imprint of thermonuclear burning is carved on the surface lithium abundance of substellar-mass objects during the first few hundred million years of their evolution, leading to a sharp boundary between lithium and non-lithium brown dwarfs, so-called, the lithium test. New optical spectroscopic observations of the binaries DENIS+J063001.4-184014 and DENIS+J225210.7-173013 obtained using the 10.4-m Gran Telescopio de Canarias are reported here. They allow us to re-determine their combined optical spectral types (M9.5 and L6.5, respectively) and to search for the presence of the LiI resonance doublet. The non detection of the LiI feature in the combined spectrum of DENIS\,J063001.4$-$184014AB is converted into estimates for the depletion of lithium in the individual components of this binary system. In DENIS\,J225210.7$-$173013AB we report the detection of a weak LiI feature which we tentatively ascribe as arising from the contribution of the T3.5-type secondary. Combining our results with data for seven other brown dwarf binaries in the literature treated in a self-consistent way, we confirm that there is indeed a sharp transition in mass for lithium depletion in brown dwarfs, as expected from theoretical calculations. We estimate such mass boundary is observationally located at 51.48$^{+0.22}_{-4.00}$ $M_\mathrm{Jup}$, which is lower than the theoretical determinations.

[abridged] We study how mock-observed stellar morphological and structural properties of massive galaxies are built up between $z=0.5$ and $z=3$ in the TNG50 cosmological simulation. We generate mock images with the properties of the CANDELS survey and derive Sersic parameters and optical rest-frame morphologies as usually done in the observations. Overall, the simulation reproduces the observed evolution of the abundances of different galaxy morphological types of star-forming and quiescent galaxies. The $\log{M_*}-\log R_e$ and $\log{M_*}-\log\Sigma_1$ relations of the simulated star-forming and quenched galaxies also match the observed slopes and zeropoints to within 1-$\sigma$. In the simulation, galaxies increase their observed central stellar mass density ($\Sigma_1$) and transform in morphology from irregular/clumpy systems to normal Hubble-type systems in the Star Formation Main Sequence at a characteristic stellar mass of $\sim 10^{10.5}~M_\odot$. This morphological transformation is connected to the activity of the central Super Massive Black Holes (SMBHs). At low stellar masses ($10^9$ < $M_*/M_\odot$ < $10^{10}$) SMBHs grow rapidly, while at higher mass SMBHs switch into the kinetic feedback mode and grow more slowly. During this low-accretion phase, SMBH feedback leads to the quenching of star-formation, along with a simultaneous growth in $\Sigma_1$. More compact massive galaxies grow their SMBHs faster than extended ones of the same mass and end up quenching earlier. In the TNG50 simulation, SMBHs predominantly grow via gas accretion before galaxies quench, and $\Sigma_1$ increases substantially after SMBH growth slows down. The simulation predicts therefore that quiescent galaxies have higher $\Sigma_1$ values than star-forming galaxies for the same SMBH mass, which disagrees with alternative models, and may potentially be in tension with some observations.

Context. The correct modeling of the scattering polarization signals observed in several strong resonance lines requires taking partial frequency redistribution (PRD) phenomena into account. Aims. This work aims at assessing the impact and the range of validity of the angle-averaged AA approximation with respect to the general angle-dependent (AD) treatment of PRD effects in the modeling of scattering polarization in strong resonance lines, with focus on the solar Ca i 4227 {\AA} line. Methods. Spectral line polarization is modeled by solving the radiative transfer problem for polarized radiation, under nonlocal thermodynamic equilibrium conditions, taking PRD effects into account, in static one-dimensional semi-empirical atmospheric models presenting arbitrary magnetic fields. The problem is solved through a two-step approach. In step 1, the problem is solved for intensity only, considering a multi-level atom. In step 2, the problem is solved including polarization, considering a two-level atom with an unpolarized and infinitely sharp lower level, and fixing the lower level population calculated at step 1. Results. The results for the Ca i 4227 {\AA} line show a good agreement between the AA and AD calculations for the Q/I and U/I wings signals. However, AA calculations reveal an artificial trough in the line-core peak of the linear polarization profiles, whereas AD calculations show a sharper peak in agreement with observations. Conclusions. An AD treatment of PRD effects is essential to correctly model the line-core peak of the scattering polarization signal of the Ca i 4227 {\AA} line. By contrast, in the considered static case, the AA approximation seems to be suitable to model the wing scattering polarization lobes and their magnetic sensitivity through magneto-optical effects.

Neptune's moon Triton shares many similarities with Pluto, including volatile cycles of N2, CH4 and CO, and represents a benchmark case for the study of surface-atmosphere interactions on volatile-rich KBOs. Within the context of New Horizons observations of Pluto as well as recent Earth-based observations of Triton, we adapt a Plutonian VTM to Triton, and test its ability to simulate its volatile cycles, thereby aiding our understanding of its climate. We present VTM simulations exploring the volatile cycles on Triton over long-term and seasonal timescales for varying model parameters. We explore what scenarios and model parameters allow for a best match of the available observations. In particular, our set of observational constraints include Voyager 2 observations, ground-based NIR (0.8 to 2.4 {\mu}m) disk-integrated spectra and the evolution of surface pressure as retrieved from stellar occultations. Our results show that Triton's poles act as cold traps for volatile ices and favor the formation of polar caps extending to lower latitudes through glacial flow. As previously evidenced by other VTMs, North-South asymmetries in surface properties can favor the development of one cap over the other. Our best-case simulations are obtained for a global reservoir of N2 ice thicker than 200 m and a bedrock thermal inertia larger than 500 SI. The large N2 ice reservoir implies a permanent N2 southern cap extending to the equator. Our results also suggest that a small permanent polar cap exists in the northern (currently winter) hemisphere if the internal heat flux remains radiogenic (< 3 mW m-2). Finally, we provide predictions for the evolution of ice distribution, surface pressure, CO and CH4 atmospheric mixing ratios in the next decades. We also model the thermal lightcurves of Triton in 2022, which serve as predictions for future JWST observations.

We introduce CALAMITY, a precision bandpass calibration method for radio interferometry. CALAMITY can solve for direction independent gains with arbitrary frequency structure to the high precision required for 21 cm cosmology with minimal knowledge of foregrounds or antenna beams and does not require any degree of redundancy (repeated identical measurements of the same baseline). We have achieved this through two key innovations. Firstly, we model the foregrounds on each baseline independently using a flexible and highly efficient set of basis functions that have minimal overlap with 21 cm modes and enforce spectral smoothness in the calibrated foregrounds. Secondly, we use an off-the-shelf GPU accelerated API (TENSORFLOW) to solve for per-baseline foregrounds simultaneously with per-frequency antenna gains in a single optimization loop. GPU acceleration is critical for our technique to be able to solve for the large numbers of foreground and gain parameters simultaneously across all frequencies for an interferometer with $\gtrsim 10$ antennas in a reasonable amount of time. In this paper, we give an overview of our technique and using realistic simulations and demonstrate its performance in solving for and removing pathological gain structures to the level necessary to measure fluctuations in the 21 cm emission field from Hydrogen gas during the Cosmic Dawn and Reionization. If you want to start using CALAMITY now, you can find a tutorial notebook at https://github.com/aewallwi/calamity/blob/main/examples/Calamity_Tutorial.ipynb .

We show that the galaxy 4-Point Correlation Function (4PCF) can test for cosmological parity violation. The detection of cosmological parity violation would reflect previously unknown forces present at the earliest moments of the Universe. Recent developments both in rapidly evaluating galaxy $N$-Point Correlation Functions (NPCFs) and in determining the corresponding covariance matrices make the search for parity violation in the 4PCF possible in current and upcoming surveys such as those undertaken by Dark Energy Spectroscopic Instrument (DESI), the $Euclid$ satellite, and the Vera C. Rubin Observatory (VRO).

The broad high-energy spectral component in blazars is usually attributed to various inverse Compton scattering processes in the relativistic jet, but has not been clearly identified in most cases due to degeneracies in physical models. AP Librae, a low-synchrotron-peaking BL Lac object (LBL) detected in 2015 by H.E.S.S. at very high energies (VHE; $>$ 0.5 TeV), has an extremely broad high-energy spectrum, covering $\sim$ 9 decades in energy. Standard synchrotron self-Compton models generally fail to reproduce the VHE emission, which has led to the suggestion that it might arise not from the blazar core, but on kiloparsec scales from inverse Compton scattering of cosmic microwave background (CMB) photons by a still-relativistic jet (IC/CMB). IC/CMB models for the TeV emission of AP Librae in prior works have implied a high level of infrared emission from the kpc-scale jet. With newly obtained Hubble Space Telescope imaging, we obtain a deep upper limit on the kpc-scale jet emission at 1.6 $\mu$m, well below the expected level. High-resolution ALMA imaging in bands 3-9 reveals a residual dust disk signature after core subtraction, with a clearly thermal spectrum, and an extent ($\sim$500 pc) which matches with a non-jet residual emission seen after PSF subtraction in our 1.6 $\mu$m HST imaging. We find that the unusually broad GeV and VHE emission in AP Librae can be reproduced through the combined IC scattering of photons from the CMB and the dust disk, respectively, by electrons in both the blazar core and sub-kpc jet.

We measure the absolute proper motion of Leo I using a WFPC2/HST data set that spans up to 10 years, to date the longest time baseline utilized for this satellite. The measurement relies on ~ 2300 Leo I stars located near the center of light of the galaxy; the correction to absolute proper motion is based on 174 Gaia EDR3 stars and 10 galaxies. Having generated highly-precise, relative proper motions for all Gaia EDR3 stars in our WFPC2 field of study, our correction to the absolute EDR3 system does not rely on these Gaia stars being Leo I members. This new determination also benefits from a recently improved astrometric calibration of WFPC2. The resulting proper-motion value, (mu_alpha, mu_delta) = (-0.007 +- 0.035, -0.119 +-0.026) mas/yr is in agreement with recent, large-area, Gaia EDR3-based determinations. We discuss all the recent measurements of Leo I's proper motion and adopt a combined, multi-study average of (mu_alpha_3meas, mu_delta_3meas) = (-0.036 +- 0.016, -0.130 +- 0.010) mas/yr. This value of absolute proper motion for Leo I indicates its orbital pole is well aligned with that of the Vast Polar Structure, defined by the majority of the brightest dwarf-spheroidal satellites of the Milky Way.

Coronal mass ejections (CMEs) are huge expulsions of magnetized matter from the Sun and stars, traversing space with speeds of millions of kilometers per hour. Solar CMEs can cause severe space weather disturbances and consumer power outages on Earth, whereas stellar CMEs may even pose a hazard to the habitability of exoplanets. While CMEs ejected by our Sun can be directly imaged by white-light coronagraphs, for stars this is not possible. So far, only a few candidates for stellar CME detections are reported. Here we demonstrate a different approach, based on sudden dimmings in the extreme-ultraviolet (EUV) and X-ray emission caused by the CME mass loss. We report dimming detections associated with flares on cool stars, indicative of stellar CMEs and benchmarked by Sun-as-a-star EUV measurements. This study paves the way for comprehensive detections and characterizations of CMEs on stars, important for planetary habitability and stellar evolution.

Tidal evolution of eccentric binary systems containing at least one massive main-sequence (MS) star plays an important role in the formation scenarios of merging compact-object binaries. The dominant dissipation mechanism in such systems involves tidal excitation of outgoing internal gravity waves at the convective-radiative boundary and dissipation of the waves at the stellar envelope/surface. We have derived analytical expressions for the tidal torque and tidal energy transfer rate in such binaries for arbitrary orbital eccentricities and stellar rotation rates. These expressions can be used to study the spin and orbital evolution of eccentric binaries containing massive MS stars, such as the progenitors of merging neutron star binaries. Applying our results to the PSR J0045-7319 system, which has a massive B-star companion and an observed, rapidly decaying orbit, we find that for the standard radius of convective core based on non-rotating stellar models, the B-star must have a significant retrograde and differential rotation in order to explain the observed orbital decay rate. Alternatively, we suggest that the convective core may be larger as a result of rapid stellar rotation and/or mass transfer to the B-star in the recent past during the post-MS evolution of the pulsar progenitor.

We report the results of our follow-up campaign for the neutron star - black hole (NSBH) merger GW200115 detected during the O3 run of the Advanced LIGO and Advanced Virgo detectors. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering ~20% of the total probability area down to a limiting magnitude of $w$=20.5 AB at ~23 h after the merger. Our search for counterparts returns a single candidate (AT2020aeo), likely not associate to the merger. In total, only 25 sources of interest were identified by the community and later discarded as unrelated to the GW event. We compare our upper limits with the emission predicted by state-of-the-art kilonova simulations and disfavor high mass ejecta (>0.1$M_{\odot}$), indicating that the spin of the system is not particularly high. By combining our optical limits with gamma-ray constraints from $Swift$ and $Fermi$, we disfavor the presence of a standard short duration burst for viewing angles $\lesssim$15 deg from the jet axis. Our conclusions are however limited by the large localization region of this GW event, and accurate prompt positions remain crucial to improving the efficiency of follow-up efforts.

Luminous Blue Variables (LBVs) are massive stars that are likely to be a transitionary phase between O stars and hydrogen-free classical Wolf-Rayet stars. The variability of these stars has been an area of study for both professional and amateur astronomers for more than a century. In this paper, we present five years of precision photometry of the classical LBV P Cygni taken with the BRITE-Constellation nanosatellites. We have analyzed these data with Fourier analysis to search for periodicities that could elucidate the drivers of variability for these stars. These data show some long-timescale variability over the course of all six calendar years of observations, but the frequencies needed to reproduce the individual light curves are not consistent from one year to the next. These results likely show that there is no periodic phenomenon present for P Cygni, meaning that the variability is largely stochastic. We interpret the data as being caused by internal gravity waves similar to those seen in other massive stars, with P Cygni exhibiting a larger amplitude and lower characteristic frequency than the main-sequence or blue supergiant stars previously studied. These results show evidence that LBVs may be an extrapolation of the blue supergiants, which have previously been shown to be an extension of main-sequence stars in the context of the stochastic low-frequency photometric variability.

We present 5--14~$\mu$m spectra at two different positions across the Orion Bar photodissociation region (PDR) obtained with the Infrared Spectrograph onboard the Spitzer Space Telescope and 3.3~$\mu$m PAH observations obtained with the Stratospheric Observatory for Infrared Astronomy (SOFIA). We aim to characterize emission from Polycyclic Aromatic Hydrocarbon (PAH), dust, atomic and molecular hydrogen, argon, sulfur, and neon as a function of distance from the primary illuminating source. We find that all the major PAH bands peak between the ionization front and the PDR front, as traced by H$_{2}$, while variations between these bands become more pronounced moving away from this peak into the face-on PDRs behind the PDR front and at the backside of the \HII\, region. While the relative PAH intensities are consistent with established PAH characteristics, we report unusual behaviours and attribute these to the PDR viewing angle and the strength of the FUV radiation field impinging on the PDRs. We determine the average PAH size which varies across the Orion Bar. We discuss subtle differences seen between the cationic PAH bands and highlight the photo-chemical evolution of carbonaceous species in this PDR environment. We find that PAHs are a good tracer of environmental properties such as the strength of the FUV radiation field and the PAH ionization parameter.

A model of the optical light detected following the merger of two neutron stars reveals polarization to be a unique probe of the geometry of the kilonova explosion that accompanied the gravitational waves.

We study the sizes of galaxies in the Epoch of Reionization using a sample of ~100,000 galaxies from the BlueTides cosmological hydrodynamical simulation from z=7 to 11. We find an inverse relationship between stellar mass and the size measured from stellar mass maps, suggesting that the most massive galaxies are more compact and dense than lower mass galaxies, which have flatter mass distributions. We find a mildly negative relation between intrinsic far-ultraviolet luminosity and size, while we find a positive size--luminosity relation when measured from dust-attenuated images. This suggests that dust is the predominant cause of the observed positive size--luminosity relation, with dust preferentially attenuating bright sight lines resulting in a flatter emission profile and thus larger measured sizes. We study the size--luminosity relation across the rest-frame ultraviolet and optical, and find that the slope decreases at longer wavelengths; this is a consequence of the relation being caused by dust, which produces less attenuation at longer wavelengths. We find that the far-ultraviolet size--luminosity relation shows mild evolution from z=7 to 11, and galaxy size evolves with redshift as $R\propto(1+z)^{-m}$, where $m=0.559\pm0.008$. Finally, we investigate the sizes of z=7 quasar host galaxies, and find that while the intrinsic sizes of quasar hosts are small relative to the overall galaxy sample, they have comparable sizes when measured from dust-attenuated images.

We present 75 near-infrared (NIR; 0.8$-$2.5 $\mu$m) spectra of 34 stripped-envelope core-collapse supernovae (SESNe) obtained by the Carnegie Supernova Project-II (CSP-II), encompassing optical spectroscopic Types IIb, Ib, Ic, and Ic-BL. The spectra range in phase from pre-maximum to 80 days past maximum. This unique data set constitutes the largest NIR spectroscopic sample of SESNe to date. NIR spectroscopy provides observables with additional information that is not available in the optical. Specifically, the NIR contains the resonance lines of He I and allows a more detailed look at whether Type Ic supernovae are completely stripped of their outer He layer. The NIR spectra of SESNe have broad similarities, but closer examination through statistical means reveals a strong dichotomy between NIR "He-rich" and "He-poor" SNe. These NIR subgroups correspond almost perfectly to the optical IIb/Ib and Ic/Ic-BL types, respectively. The largest difference between the two groups is observed in the 2 $\mu$m region, near the He I $\lambda$2.0581 $\mu$m line. The division between the two groups is not an arbitrary one along a continuous sequence. Early spectra of He-rich SESNe show much stronger He I $\lambda$2.0581 $\mu$m absorption compared to the He-poor group, but with a wide range of profile shapes. The same line also provides evidence for trace amounts of He in half of our SNe in the He-poor group.

We analyzed 82,548 carefully curated observations of 4,420 asteroids with Wide-field Infrared Survey Explorer (WISE) 4-band data to produce estimates of diameters and infrared emissivities. We also used these diameter values in conjunction with absolute visual magnitudes to infer estimates of visible-band geometric albedos. We provide solutions to 131 asteroids not analyzed by the NEOWISE team and to 1,778 asteroids not analyzed with 4-band data by the NEOWISE team. Our process differs from the NEOWISE analysis in that it uses an accurate solar flux, integrates the flux with actual bandpass responses, obeys Kirchhoff's law, and does not force emissivity values in all four bands to an arbitrary value of 0.9. We used a regularized model fitting algorithm that yields improved fits to the data. Our results more closely match stellar occultation diameter estimates than the NEOWISE results by a factor of ~2. Using 24 high-quality stellar occultation results as a benchmark, we found that the median error of 4-infrared-band diameter estimates in a carefully curated data set is 9.3%. Our results also suggest the presence of a size-dependent bias in the NEOWISE diameter estimates, which may pollute estimates of asteroid size distributions and slightly inflate impact hazard risk calculations. For more than 90% of asteroids in this sample, the primary source of error on the albedo estimate is the error on absolute visual magnitude.

Silicon abundances were determined by applying the spectrum-fitting technic to the Si II doublet lines at 6347 and 6371A for a sample of 120 main-sequence stars in the T_eff range of ~7000-14000K (comprising not only normal stars but also non-magnetic chemically peculiar stars) with an aim of investigating their behaviors (e.g., correlation with stellar parameters and abundances of other elements such as Fe or C) and the background physical mechanisms involved therein, where attention was paid to taking into account of the non-LTE effect and to assigning a reasonable value of microturbulence. The following trends were revealed from the analysis: (i) The resulting [Si/H] values, mostly ranging from ~-0.5 to ~+0.3, show a positive correlation with [Fe/H]. (ii) A kind of anti-correlation exists between Si and C as seen from the tendency of [C/Si] steeply decreasing with [Si/H]. (iii) Si abundances do not show any clear dependence upon T_eff or vsini, while Am and HgMn stars appear to show comparatively higher [Si/H] than normal stars. Although it is not straightforward to explain these observational facts, different physical processes (gas-dust separation and atomic diffusion) are likely to be intricately involved in producing these characteristic behaviors of Si composition in the surface of late A through late B dwarfs.

Tidal interactions play an important role in many astrophysical systems, but uncertainties regarding the tides of rapidly rotating, centrifugally distorted stars and gaseous planets remain. We have developed a precise method for computing the dynamical, non-dissipative tidal response of rotating planets and stars, based on summation over contributions from normal modes driven by the tidal potential. We calculate the normal modes of isentropic polytropes rotating at up to ~90% of their critical breakup rotation rates, and tabulate fits to mode frequencies and tidal overlap coefficients that can be used to compute the frequency-dependent tidal response (via potential Love numbers). Although fundamental modes (f-modes) dominate the tides at all rotation rates, we find that the strong coupling of retrograde inertial modes (i-modes) to tesseral (l>|m|) components of the tidal potential produces resonances that may be relevant to gas giants like Jupiter and Saturn. The coupling of f-modes in rapid rotators to multiple components of both the driving tidal potential and the induced gravitational field also affect the tesseral response, leading to significant deviations from treatments of rotation that neglect centrifugal distortion and high-order corrections. For very rapid rotation rates (>70% of breakup), mixing between prograde f-modes and i-modes significantly enhances the sectoral (l=|m|) tidal overlap of the latter. The tidal response of very rapidly rotating, centrifugally distorted planets or stars can also be modified by resonant sectoral f-modes that are secularly unstable via the Chandrasekhar-Friedman-Schutz (CFS) mechanism.

Accretion-induced collapse (AIC) from oxygen/neon/magnesium composition white dwarf (ONeMg WD) + stripped helium (He) star binaries is one promising channel to form peculiar neutron star objects. It has been discussed that the WD's magnetic field may alter the accretion phase in the WD binary evolution. By considering non-magnetic and sufficiently magnetized WDs, we investigate the evolution of ONeMg WD + He star binaries with detailed stellar evolution and binary population synthesis simulations. The role of the magnetically confined accretion in the possible formation pathway for like millisecond pulsars (MSPs) and magnetars is also studied. Comparing with the case of spherically symmetric accretion, the mass accumulation efficiency of the WDs is enhanced at low mass transfer rate under the magnetic confinement model. The initial parameter space of the potential AIC progenitor systems moves toward shorter orbital period and lower donor mass (but not so significantly) due to the effect of the magnetic confinement. This also allows final MSPs to have lower-mass WD companions and shorter orbital periods. There is no significant difference between the Galactic birthrates of the AIC derived with and without the magnetic confinement, which implies that the magnetic field of the WD does not dramatically change the number of ONeMg WD + He star binaries which can produce AIC. It is worth noting that these conclusions can be applied for the CO (carbon/oxgen) WD + He star binaries as progenitors of type Ia supernovae, because the accretion phases of ONeMg WDs and CO WDs are similar. The Galactic rate of magnetars possibly formed via AIC of highly magnetized WDs is $0.34\times10^{-4}\,{\rm yr}^{-1}$.

Earth's surface environment is largely influenced by its budget of major volatile elements: carbon (C), nitrogen (N), and hydrogen (H). Although the volatiles on Earth are thought to have been delivered by chondritic materials, the elemental composition of the bulk silicate Earth (BSE) shows depletion in the order of N, C, and H. Previous studies have concluded that non-chondritic materials are needed for this depletion pattern. Here, we model the evolution of the volatile abundances in the atmosphere, oceans, crust, mantle, and core through the accretion history by considering elemental partitioning and impact erosion. We show that the BSE depletion pattern can be reproduced from continuous accretion of chondritic bodies by the partitioning of C into the core and H storage in the magma ocean in the main accretion stage and atmospheric erosion of N in the late accretion stage. This scenario requires a relatively oxidized magma ocean ($\log_{10} f_{\rm O_2}$ $\gtrsim$ $\rm{IW}$$-2$, where $f_{\rm O_2}$ is the oxygen fugacity, ${\rm IW}$ is $\log_{10} f_{\rm O_2}^{\rm IW}$, and $f_{\rm O_2}^{\rm IW}$ is $f_{\rm O_2}$ at the iron-w\"{u}stite buffer), the dominance of small impactors in the late accretion, and the storage of H and C in oceanic water and carbonates in the late accretion stage, all of which are naturally expected from the formation of an Earth-sized planet in the habitable zone.

A three-component Stackel model of the Galaxy, including the bulge, disk, and halo, is constructed. Parameter estimates of the potential are obtained as a result of fitting the model rotation curve to azimuthal velocities found from data on trigonometric parallaxes and spatial velocities of masers. The fitting method takes into account the measurement and natural dispersions of azimuthal velocities and uses an algorithm for excluding objects with excessive residuals. In order to obtain more uniform samples, the objects were divided into two groups: masers associated with high-mass star forming regions and masers of other types. A significant kinematic inhomogeneity of these groups was identified and taken into account: the azimuthal velocity dispersion is $\sigma_{0,1}=4.3\pm 0.4$~km\,s$^{-1}$, in the first group and $\sigma_{0,2}=15.2\pm1.3$~km\,s$^{-1}$ in the second. After constructing the model of the Galactic-plane potential, it was generalized to the entire space under the assumption of the existence of a third quadratic integral of motion. When reconstructing the Galactic rotation curve in detail, the used algorithm gives an analytical expression for the Stackel potential, which significantly simplifies the task of constructing the Galaxy's phase density model in the Stackel approximation. In order to make the Stackel model more realistic, one needs to develop methods of direct account of data on the vertical distribution of density in the Galaxy.

Context. Puppis A is a medium-age supernova remnant (SNR), which is visible as a very bright extended X-ray source. While numerous studies have investigated individual features of the SNR, at this time, no comprehensive study of the entirety of its X-ray emission exists. Aims. Using field-scan data acquired by the SRG/eROSITA telescope during its calibration and performance verification phase, we aim to investigate the physical conditions of shocked plasma and the distribution of elements throughout Puppis A. Methods. Using broad- and narrow-band imaging, we investigate the large-scale distribution of absorption and plasma temperature as well as typical emission lines. This approach is complemented by spatially resolved spectral analysis of the shocked plasma in Puppis A, for which we divide the SNR into around 700 distinct regions, resulting in maps of key physical quantities over its extent. Results. We find a strong peak of foreground absorption in the southwest quadrant, which in conjunction with high temperatures at the northeast rim creates the well-known strip of hard emission crossing Puppis A. Furthermore, using the observed distribution of ionization ages, we attempt to reconstruct the age of the shock in the individual regions. We find a quite recent shock interaction for the prominent northeast filament and ejecta knot, as well as for the outer edge of the bright eastern knot. Finally, elemental abundance maps reveal only a single clear enhancement of the plasma with ejecta material, consistent with a previously reported knot, and no obvious ejecta enrichment in the remainder of the SNR. Within this region, we confirm the spatial separation of silicon-rich ejecta from those dominated by lighter elements.

The redshift drift (also known as the Sandage Test) is a model-independent probe of fundamental cosmology, enabling us to watch the universe expand in real time, and thereby to confirm (or not) the recent acceleration of the universe without any model-dependent assumptions. On the other hand, by choosing a fiducial model one can also use it to constrain the model parameters, thereby providing a consistency test for results obtained with other probes. The drift can be measured by the Extremely Large Telescope and also by the full SKA. Recently two alternative measurement methods have been proposed: the cosmic accelerometer, and the differential redshift drift. Here we summarize a comparative analysis of the various methods and their possible outcomes, using both Fisher Matrix and MCMC techniques. We find that no single method is uniformly better than the others. Instead, their comparative performance depends both on experimental parameters (including the experiment time and redshift at which the measurement is made) and also on the scientific goal (e.g., detecting the drift signal with high statistical significance, constraining the matter density, or constraining the dark energy properties). In other words, the experiment should be optimized for the preferred scientific goal.

In this paper we present PyCF3, a python client for the cosmological distance-velocity calculator CosmicFlow-3. The project has a cache and retry system designed with the objective of reducing the stress on the server and mitigating the waiting times of the users in the calculations. We also address Quality Assurance code standards and availability of the code.

This note collects together useful unit conversions and numerical values from early universe cosmology. It is a quick reference that can be used to make easy order-of-magnitude estimates. Included are tables for unit conversions, the thermal history of the universe, and collected properties of astronomical objects. The note also introduces a modifiable Mathematica package NaturalUnits, which makes it easy to convert between natural and physical units.

The COSmic Monopole Observer (COSMO) is an experiment to measure low-level spectral distortions in the isotropic component of the Cosmic Microwave Background (CMB). Deviations from a pure blackbody spectrum are expected at low level ($<$ 1 ppm) due to several astrophysical and cosmological phenomena, and promise to provide important independent information on the early and late phases of the universe. They have not been detected yet, due to the extreme accuracy required, the best upper limits being still those from the COBE-FIRAS mission. COSMO is based on a cryogenic differential Fourier Transform Spectrometer, measuring the spectral brightness difference between the sky and an accurate cryogenic blackbody. The first implementation of COSMO, funded by the Italian PRIN and PNRA programs, will operate from the Concordia station at Dome-C, in Antarctica, and will take advantage of a fast sky-dip technique to get rid of atmospheric emission and its fluctuations, separating them from the monopole component of the sky brightness. Here we describe the instrument design, its capabilities, the current status. We also discuss its subsequent implementation in a balloon-flight, which has been studied within the COSMOS program of the Italian Space Agency.

Investigations of some type Iax supernovae have led to the suggestion that their ejecta must be layered to some degree. Such an ejecta structure has been argued as inconsistent with the well-mixed composition predicted by pure deflagrations. Based on explosion models, we create toy models in which the ejecta are artificially stratified and progressively mixed until a uniform composition is obtained. We find that models that are heavily mixed, containing burned and unburned material at all velocities, produce reasonably good agreement with SN 2012Z, for which a layered structure has been suggested. We also discuss how existing ejecta compositions determined for type Iax supernovae do not necessarily contradict pure deflagration models and may be consistent with a steeper density profile. We investigate previous claims that differences in line profile shapes may be due to strong blending, by presenting a series of models with different plasma states. These models indicate that blending could indeed explain differences in the observed profiles. Alternatively, stratification could also explain such differences, however all of our models indicate that this does not necessarily require stratification in abundance. Sufficient stratification in ionisation state can be achieved even for a well-mixed model. Based on our analysis, we demonstrate that there is insufficient evidence to suggest the ejecta of type Iax supernovae must be layered and therefore argue the pure deflagration scenario is not ruled out, even for the brightest type Iax supernovae. Our analysis does not indicate the ejecta cannot be layered to some degree, but observations within days of explosion are necessary to determine the extent to which the outer ejecta could be layered.

Here, we report the new detection of three shock fronts using archival \textit{Chandra} X-ray observations in Abell 1914, which also hosts a radio halo, a radio phoenix, and a head-tail galaxy. In this study, we report the X-ray shock front at the position of the radio phoenix, which further strengthens the scenario that radio phoenix traces old plasma that gets lit up when compressed by shock passage. We further analyze the thermodynamic structure of the cluster in detail. We create temperature maps of A1914 using three different techniques, Adaptive Circular Binning (ACB), Weighted Voronoi Tessellations (WVT), and Contour binning (Contbin) method. These thermodynamic maps, along with the pseudo pressure and the pseudo entropy maps for the cluster, are evidence of disturbed morphology produced by multiple merger events. These merger events create cluster-wide turbulence, which may re-accelerate the relativistic particles and result in a radio halo within the cluster. Further, comparing X-ray and radio images reveals that the radio halo is contained within two X-ray shock fronts. Our analysis suggests that A1914 has both equatorial shock and axial shock within the cluster's ICM. We proposed a dual merging scenario based on the shock position and analysis of the thermodynamic maps obtained from the deeper \textit{Chandra} X-ray observations.

Probing the existence of hypothetical particles beyond the Standard model often deals with extreme parameters: large energies, tiny cross-sections, large time scales, etc. Sometimes laboratory experiments can test required regions of parameter space, but more often natural limitations leads to poorly restrictive upper limits. In such cases astrophysical studies can help to expand the range of values significantly. Among astronomical sources, used in interests of fundamental physics, compact objects -- neutron stars and white dwarfs, -- play a leading role. We review several aspects of astroparticle physics studies related to observations and properties of these celestial bodies. Dark matter particles can be collected inside compact objects resulting in additional heating or collapse. We summarize regimes and rates of particle capturing as well as possible astrophysical consequences. Then we focus on a particular type of hypothetical particles -- axions. Their existence can be uncovered due to observations of emission originated due to Primakoff process in magnetospheres of neutron stars or white dwarfs. Alternatively, they can contribute to cooling of these compact objects. We present results in these areas, including upper limits based on recent observations.

We report the discovery of a close quasar pair candidate at $z=5.66$, J2037--4537. J2037--4537 is resolved into two quasar images at the same redshift in ground-based observations. Followup spectroscopy shows significant differences in both the continuum slopes and emission line properties of the two images. The two quasar images have a projected separation of $1\farcs24$ ($7.3\text{~kpc}$ at $z=5.66$) and a redshift difference of $\Delta z\lesssim0.01$. High-resolution images taken by {\em Hubble Space Telescope} do not detect the foreground lensing galaxy. The observational features of J2037--4537 strongly disfavor the lensing hypothesis. If J2037--4537 is a physical quasar pair, it indicates a quasar clustering signal of $\sim10^5$ at a separation of $\sim10$ proper kpc (pkpc), and gives the first observational constraint on the pair fraction of $z>5$ quasars, $f_\text{pair}(r<30\text{~pkpc})>0.3\%$. The properties of J2037--4537 are consistent with those of merger-triggered quasar pairs in hydrodynamical simulations of galaxy mergers.

Exploring Cosmic History and Origins (ECHO), popularly known as `CMB-Bh$\overline{a}$rat', is a space mission that has been proposed to the Indian Space Research Organisation (ISRO) for the scientific exploitation of the Cosmic Microwave Background (CMB) at the next level of precision and accuracy. The quest for the CMB polarization $B$-mode signals, generated by inflationary gravitational waves in the very early universe, is one of the key scientific goals of its experimental design. This work studies the potential of the proposed ECHO instrumental configuration to detect the target tensor-to-scalar ratio $r \sim 10^{-3}$ at $3\sigma$ significance level, which covers the predictions of a large class of inflationary models. We investigate the performance of two different component separation pipelines, NILC and Commander, for the measurement of $r$ in presence of different physically motivated models of astrophysical foregrounds. For a simplistic foreground model (only polarized dust and synchrotron), both component separation pipelines can achieve the desired sensitivity of ECHO, i.e. $\sigma (r =0) \sim (0.4 - 0.7)\times 10^{-3}$. NILC performs better than Commander in terms of bias on recovered $r$ for complex spectral models (power-law and curved power-law) of the synchrotron emission and complex dust models (dust decorrelation). Assuming 84 % delensing, we can achieve an improvement of $\sigma (r = 0)$ by approximately 50 % as compared to the results obtained for the same configuration without any lensing correction.

We attempt to constrain the efficiency of additional transport or mixing processes that reduce the effect of atomic diffusion in stellar atmospheres. We apply spectral synthesis methods to spectra observed with the GIRAFFE spectrograph on the VLT to estimate abundances of Mg, Ti, Fe, and Ba in stars in the metal-poor globular cluster M30. To the abundances we fit trends of abundances predicted by stellar evolution models assuming different efficiencies of additional transport or mixing processes. The fitting procedure attempts to take into account the effects of parameter-dependent systematic errors in the derived abundances. We find that the parameter $T_0$, which describes the efficiency of additional transport or mixing processes, can almost certainly be constrained to the narrow range $\log_{10}{\left( T_0 / \left[ \mathrm{K} \right] \right)}$ between $6.09$ and $6.2$. This corresponds to decreased abundances for stars at the main sequence turn-off point compared to the red giant branch by $0.2\,\mathrm{dex}$ for Mg, $0.1\,\mathrm{dex}$ for Fe, and $0.07\,\mathrm{dex}$ for Ti. We also find that while our estimates do have non-negligible systematic errors stemming from the continuum placement and the assumed microturbulence, our method can take them into account. Our results partly amend the results of an earlier paper in this article series, that tentatively used a value of $\log_{10}{\left( T_0 / \left[ \mathrm{K} \right] \right)} = 6.0$ when modelling the Spite plateau of lithium. To more easily distinguish physical effects from systematic errors, we recommend that studies of this kind focus on elements for which the expected surface abundances as functions of effective temperature have a distinct structure and cover a wide range.

Context. High-resolution observations of the solar photosphere reveal the presence of fine structures, in particular the so-called magnetic bright points (MBPs), which are small-scale features associated with strong magnetic field regions of the order of kilogauss (kG). It is especially relevant to study these magnetic elements, which are extensively detected at all moments of the solar cycle, in order to establish their contribution to the behaviour of the solar atmosphere, and ultimately a plausible role within the coronal heating problem. Aims. We aim to characterise the size and velocity distributions of MBPs in the solar photosphere in two different datasets of quiet Sun images acquired with the Solar Optical Telescope SOT/Hinode and the High-resolution Fast Imager HiFI/GREGOR, in the G-band (4308 angstroms). Methods. In order to detect the MBPs, an automatic segmentation and identification algorithm was used. Next, the identified features were tracked to measure their proper motions. Finally, a statistical analysis of hundreds of MBPs was carried out, generating histograms for areas, diameters, and horizontal velocities. Results. This work establishes that areas and diameters of MBPs display log-normal distributions that are well fitted by two different components, whereas the velocity vector components follow Gaussians, and the vector magnitude follows a Rayleigh distribution again revealing a two-component composition for all vector elements. Conclusions. The results can be interpreted as due to the presence of two different populations of MBPs in the solar photosphere, one likely related to stronger network magnetic flux elements and the other one to weaker intranetwork flux elemens. In particular, this work concludes on the effect of the different spatial resolutions of the GREGOR and Hinode telescopes, affecting detections and average values.

We address the physical origin of the ultrarelativistic prompt emission (UPE) phase of GRB 190114C observed in the interval 1.9-3.99 s, by the Fermi-GBM in 10 keV-10 MeV . Thanks to high S/N ratio of Fermi-GBM data, a time resolved spectral analysis has evidenced a sequence of similar blackbody plus cutoff power-law spectra, on ever decreasing time intervals during the entire UPE phase. We assume that during the UPE phase, the inner engine of the GRB, composed of a Kerr black hole and a uniform test magnetic field B0, aligned with the BH rotation axis, operates in an overcritical field. We infer an $e^+e^-$ pair electromagnetic plasma in presence of a baryon load, a PEMB pulse, originating from a vacuum polarization quantum process in the inner engine. This initially optically thick plasma self-accelerates, giving rise at the transparency radius to the MeV radiation observed by Ferm-GBM. At trf > 3.99 s, the electric field becomes undercritical, and the inner engine operates in the classical electrodynamics regime and generate the GeV emission. During both the quantum and the classical electrodynamics processes, we determine the time varying mass and spin of the Kerr BH in the inner engine, fulfilling the Christodoulou-Hawking-Ruffini mass-energy formula. For the first time, we quantitatively show how the inner engine, by extracting the rotational energy of the Kerr BH, produces a series of PEMB pulses. We follow the quantum vacuum polarization process in sequences with decreasing time bins. We compute the Lorentz factors, the baryon loads and the radii at transparency, as well as the value of the magnetic field, assumed to be constant in each sequence. The fundamental hierarchical structure, linking the quantum electrodynamics regime to the classical electrodynamics regime, is characterized by the emission of blackholic quanta with a timescale $t=10^{-9}$s, and energy $E=10^{45}$ erg.

Using fully-kinetic plasma simulations, we study the non-resonant (Bell) streaming instability driven by energetic leptons. We identify the necessary conditions to drive it and the differences from the standard proton-driven case in both linear and saturated stages. A simple analytic theory is presented to explain simulations. Our findings are crucial for understanding the phenomenology of astrophysical environments where only electrons may be accelerated (e.g., oblique shocks) or where relativistic pairs are produced (e.g., around pulsar wind nebulae).

The theory of magnetohydrodynamic (MHD) turbulence predicts that Alfv\'enic and slow-mode-like compressive fluctuations are energetically decoupled at small scales in the inertial range. However, when the magnetorotational instability (MRI) drives the turbulence, it is difficult to resolve numerically the scale at which both types of fluctuations start to be decoupled because the MRI energy injection occurs in a broad range of wavenumbers, and both types of fluctuations are usually expected to be coupled even at relatively small scales. In this study, we focus on MRI turbulence threaded by a near-azimuthal mean magnetic field, which is naturally produced by the differential rotation of a disc. We show that the decoupling scales are reachable using a reduced MHD model that includes differential-rotation effects. In this model, the Alfv\'enic and compressive fluctuations are coupled only through the linear terms that are proportional to the angular velocity of the accretion disc. We numerically solve for the turbulence in this model and show that the Alfv\'enic and compressive fluctuations are decoupled at the small scales of our simulations as the nonlinear energy transfer dominates the linear coupling below the MRI-injection scale. We show that the energy flux of compressive fluctuations contained in the small scales is almost double that of Alfv\'enic fluctuations. Combining this finding with our recent results found on hybrid-gyrokinetic simulations leads to an ion-to-electron heating prescription that takes into account both the driving of turbulence via MRI at MHD scales and the dissipation at kinetic scales. This prescription suggests that ions are heated at least twice as strongly as the electrons in the portion of accretion disc with near-azimuthal magnetic field, which can be a useful model in interpreting observations of hot accretion discs by the Event Horizon Telescope.

Numerical simulations with varying realism indicate an emergent principle -- multiphase condensation and large cavity power occur when the ratio of the cooling time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) falls below a threshold value close to 10. Observations indeed show cool-core signatures when this ratio falls below 20-30, but the prevalence of cores with \tctf~ratio below 10 is rare as compared to simulations. In X-ray observations, we obtain projected spectra from which we have to infer radial gas density and temperature profiles. Using idealized models of X-ray cavities and multiphase gas in the core and 3-D hydro jet-ICM simulations, we quantify the biases introduced by deprojection based on the assumption of spherical symmetry in determining $t_{\rm cool}/t_{\rm ff}$. We show that while the used methods are able to recover the $t_{\rm cool}/t_{\rm ff}$ ratio for relaxed clusters, they have an uncertainty of a factor of $2-3$ in systems containing large cavities ($\gtrsim 20$ kpc). We also show that the mass estimates from these methods, in absence of X-ray spectra close to the virial radius, suffer from a degeneracy between the virial mass ($M_{200}$) and the concentration parameter ($c$) in the form of $M_{200}\: c^2 \approx$ constant. Additionally, lack of soft-X-ray ($\lesssim 0.5$ keV) coverage and poor spatial resolution make us overestimate min($t_{\rm cool}/t_{\rm ff}$) by a factor of few in clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$. This bias can largely explain the lack of cool-core clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$.

The radiation formation length for relativistic particles, $l_c \sim \gamma^2 \lambda$ ($\gamma$ is the Lorentz factor, $\lambda$ is the emitted wavelength), is much lager than the inter-particle distances in many astrophysical applications. This leads to the importance of plasma effects even for the high energy emission. The consequences are nontrivial: (i) averaging of the phases of the emitting particles reduces the power (a.k.a., a circle current does not emit); (ii) density fluctuations may lead to the sporadic production of coherent emission; (iii) plasma effects during assembly of a photon may lead to the suppression of the emission (Razin-Tsytovich effect for the superluminal modes), or, in the opposite limit of subluminous normal modes, to the newly discussed synchrotron super-radiance. For synchrotron emission the radiation formation length is the same for all emitted waves, $\sim c/\omega_B$ (non-relativistic Larmor length); for curvature emission it is $R/\gamma$ - macroscopically long in pulsar magnetospheres (e.g., kilometers for radio). The popular model of coherent curvature emission by bunches", with kilometers-long radiation formation length, particles swinging-out in a rotating magnetospheres before they finish emitting a wave, extreme requirements on the momentum spreads, and demands on the electric energy needed to keep the electrostatically repulsing charges together, all make that model internally inconsistent. Long radiation formation lengths affect how emission from PIC simulations should be interpreted: phases of the emitted wave should be added over the radiation formation length, not just the powers from the instantaneous acceleration of each particle.

Kinematics of the molecular clouds in the star forming complex M17 is studied using the high-resolution CO-line mapping data at resolution ($20" \sim 0.2$ pc) with the Nobeyama 45-m telescope. The northern molecular cloud of M17, which we call the molecular "lobe", is shown to have an elongated shell structure around a top-covered cylindrical cavity. The lobe is expanding at $\sim 5$ \kms in the minor axis direction, and at $ \sim 3/\cos \ i$ km s$^{-1}$ in the major axis direction, where $i$ is the inclination of the major axis. The kinetic energy of the expanding motion is on the order of $\sim 3\times 10^{49}$ ergs. We show that the lobe is a backyard structure having the common origin to the denser molecular "horn" flowing out from NGC 6618 toward the south, so that the lobe and horn compose a bipolar outflow. Intensity distributions across the lobe and horn show a double-peak profile typical for a cylinder around a cavity. Position-velocity diagrams (PVD) across the lobe and horn exhibit open ring structure with the higher- and/or lower-velocity side(s) being lacking or faded. This particular PVD behavior can be attributed to outflow in a conical cylinder with the flow velocity increasing toward the lobe and horn axes.

Several recent investigations indicate the existence of gender-related systematic trends in the peer review of proposals for observations on astronomical facilities. This includes the National Radio Astronomy Observatory (NRAO) where there is evidence of a gender imbalance in the rank of proposals with male principal investigators (PIs) favored over female PIs. Since semester 2017A (17A), the NRAO has taken the following steps: (1) inform science review panels (SRPs) and the telescope time allocation committee (TAC) about the gender imbalance; and (2) increase the female representation on SRPs and the TAC to reflect the community demographics. Here we analyze SRP normalized rank-ordered scores, or linear ranks, by PI gender for NRAO observing proposals from semesters 12A-21A. We use bootstrap resampling to generate modeled distributions and the Anderson-Darling (AD) test to evaluate the probability that the linear rank distributions for male and female PIs are drawn from the same parent sample. We find that between semesters 12A-17A that male PIs are favored over female PIs (AD p-value 0.0084), whereas between semesters 17B-21A female PIs are favored over male PIs, but at a lower significance (AD p-value 0.11). Therefore the gender imbalance is currently being ameliorated, but this imbalance may have been reversed. Regardless, we plan to adopt a dual-anonymous approach to proposal review to reduce the possibility of bias to occur.

Solar flare forecasting can be realized by means of the analysis of magnetic data through artificial intelligence techniques. The aim is to predict whether a magnetic active region (AR) will originate solar flares above a certain class within a certain amount of time. A crucial issue is concerned with the way the adopted machine learning method is implemented, since forecasting results strongly depend on the criterion with which training, validation, and test sets are populated. In this paper we propose a general paradigm to generate these sets in such a way that they are independent from each other and internally well-balanced in terms of AR flaring effectiveness. This set generation process provides a ground for comparison for the performance assessment of machine learning algorithms. Finally, we use this implementation paradigm in the case of a deep neural network, which takes as input videos of magnetograms recorded by the Helioseismic and Magnetic Imager on-board the Solar Dynamics Observatory (SDO/HMI). To our knowledge, this is the first time that the solar flare forecasting problem is addressed by means of a deep neural network for video classification, which does not require any a priori extraction of features from the HMI magnetograms.

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $\mu$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust.

The bulk chemical composition and interior structure of rocky exoplanets are of fundamental importance to understanding their long-term evolution and potential habitability. Observations of the chemical compositions of the solar system rocky bodies and of other planetary systems have increasingly shown a concordant picture that the chemical composition of rocky planets reflects that of their host stars for refractory elements, whereas this expression breaks down for volatiles. This behavior is explained by devolatilization during planetary formation and early evolution. Here, we apply a devolatilization model calibrated with the solar system bodies to the chemical composition of our nearest Sun-like stars -- $\alpha$ Centauri A and B -- to estimate the bulk composition of any habitable-zone rocky planet in this binary system ("$\alpha$-Cen-Earth"). Through further modeling of likely planetary interiors and early atmospheres, we find that compared to Earth, such a planet is expected to have (i) a reduced (primitive) mantle that is similarly dominated by silicates albeit enriched in carbon-bearing species (graphite/diamond); (ii) a slightly larger iron core, with a core mass fraction of $38.4_{-5.1}^{+4.7}$ wt% (cf. Earth's 32.5 $\pm$ 0.3 wt%); (iii) an equivalent water-storage capacity; and (iv) a CO$_2$-CH$_4$-H$_2$O-dominated early atmosphere that resembles that of Archean Earth. Further taking into account its $\sim$ 25% lower intrinsic radiogenic heating from long-lived radionuclides, an ancient $\alpha$-Cen-Earth ($\sim$ 1.5-2.5 Gyr older than Earth) is expected to have less efficient mantle convection and planetary resurfacing, with a potentially prolonged history of stagnant-lid regimes.

The cosmic microwave background (CMB) observation by the Planck satellite precisely determines primordial curvature fluctuations on larger scales than $\mathcal O(1)\,\mathrm{Mpc}$, while the small-scale curvature fluctuation is still less constrained. The constraint on small-scale fluctuations is highly improved if we assume the standard thermal relic dark matter scenario. When small-scale fluctuations are large enough, dense regions collapse to form small halos even in a redshift $z\gtrsim10^3$, which is called ultracompact minihalos. These minihalos enhance the annihilation of the dark matter and it is constrained by observations such as extragalactic gamma rays and the CMB. We revisit the effect of minihalos formed by the small-scale density fluctuations and calculate the ionization history modified by the dark matter annihilation. We perform the Markov Chain Monte Carlo method to constrain the size of small-scale curvature fluctuations by the CMB power spectrum. It is found that the constraint from the CMB power spectrum is comparable to that from the extragalactic gamma rays. We confirm that our constraint mainly comes from the energy injection in early time ($z\gtrsim 100$) and hence it is independent of the uncertainty of minihalo properties in the late time.

White dwarfs (WDs) often show metal lines in their spectra, indicating accretion of asteroidal material. Our Sun is to become a WD in several Gyr. Here, we examine how the solar WD accretes from the three major small body populations: the main belt asteroids (MBAs), Jovian trojan asteroids (JTAs), and trans-Neptunian objects (TNOs). Owing to the solar mass loss during the giant branch, 40\% of the JTAs are lost but the vast majority of MBAs and TNOs survive. During the WD phase, objects from all three populations are sporadically scattered onto the WD, implying ongoing accretion. For young cooling ages $\lesssim 100$ Myr, accretion of MBAs predominates; our predicted accretion rate $\sim10^6$ g/s falls short of observations by two orders of magnitude. On Gyr timescales, thanks to the consumption of the TNOs that kicks in $\gtrsim 100$ Myr, the rate oscillates around $10^6-10^7$ g/s until several Gyr and drops to $\sim10^5$ g/s at 10 Gyr. Our solar WD accretion rate from 1 Gyr and beyond agrees well with those of the extrasolar WDs. We show that for the solar WD, the accretion source region evolves in an inside-out pattern. Moreover, in a realistic small body population with individual sizes covering a wide range as WD pollutants, the accretion is dictated by the largest objects. As a consequence, the accretion rate is lower by an order of magnitude than that from a population of bodies of a uniform size and the same total mass and shows greater scatter.

Neptune's incomplete ring arcs have been stable since their discovery in 1984 by stellar occultation. Although these structures should be destroyed in a few months through differential Keplerian motion, imaging data over the past couple of decades has shown that these structures are persistent. We present here the first SPHERE near-infrared observations of Neptune's ring arcs taken at 2.2 $\mu$m (BB-Ks) with the IRDIS camera at the Very Large Telescope in August 2016. The images were aligned using the ephemerides of the satellite Proteus and were suitably co-added to enhance ring and satellite signals. We analyse high-angular resolution near-infrared images of Neptune's ring arcs obtained in 2016 at the ESO VLT-UT3 with the adaptive-optics fed camera SPHERE-IRDIS. We derive here accurate mean motion values for the arcs and the nearby satellite Galatea. The trailing arcs Fraternit\'e and Egalit\'e are stable since they were last observed in 2007. Furthermore, we confirm the fading away of the leading arcs Courage and Libert\'e. Finally, we confirm the mismatch between the arcs' position and 42:43 inclined and eccentric corotation resonances with Galatea; thus demonstrating that no 42:43 corotation model works to explain the azimuthal confinement of the arcs' materiel.

The subject of this study is oscillations in the lower atmosphere in coronal-hole regions, where the conditions are favorable for propagation between the atmospheric layers. Based on spectroscopic observations in photospheric and chromospheric lines, we analysed the features of the oscillations that show signs of propagation between layers of the solar atmosphere. Using the cross-spectrum wavelet algorithm, we found that both chromospheric and photospheric signals under coronal holes share a range of significant oscillations of periods around five minutes, while the signals outside of coronal holes show no mutual oscillations in the photosphere and chromosphere. The phase shift between the layers indicates a predominantly upward propagation with partial presence of standing waves. We also tested the assumption that torsional Alfv\'en waves propagating in the corona originate in the lower atmosphere. However, the observed line-width oscillations, although similar in period to the Alfv\'en waves observed earlier in the corona of open-field regions, seem to be associated with other MHD modes. If we assume that the oscillations that we observed are related to Alfv\'en waves, then perhaps this is only through the mechanisms of the slow MHD wave transformation.

The GOLF instrument on board SoHO has been in operation for almost 25 years but aging of the instrument has now strongly affected its performance, especially in the low-frequency p-mode region. At the end of the SoHO mission, the ground-based network BiSON will remain the only facility able to perform Sun-integrated helioseismic observations. Therefore, we want to assess the helioseismic performances of an \'echelle spectrograph like SONG. Indeed, the high precision of such an instrument and the quality of the data acquired for asteroseismic purpose calls for an evaluation of the instrument ability to perform global radial-velocity measurements of the solar disk. Data acquired during the Solar-SONG 2018 observation campaign at the Teide Observatory are used to study mid- and low-frequency p modes. A Solar-SONG time series of 30-day duration is reduced with a combination of the traditional IDL iSONG pipeline and a new Python pipeline described in this paper. A mode fitting method built around a Bayesian approach is then performed on the Solar-SONG and contemporaneous GOLF, BiSON, and HMI data. For this contemporaneous time series, Solar-SONG is able to characterise p modes at a lower frequency than BiSON and GOLF (1750{\mu}Hz against 1946 and 2157 {\mu}Hz respectively), while for HMI it is possible to characterise a mode at 1686 {\mu}Hz. The decrease of GOLF sensitivity is then evaluated through the evolution of its low-frequency p-mode characterisation abilities over the years. [abridged]

The orbit of the outer satellite Alexhelios of (216) Kleopatra is already constrained by adaptive-optics astrometry, obtained with the VLT/SPHERE instrument. However, there is also a preceding occultation event in 1980 attributed to this satellite. Hereinafter, we try to link all observations, spanning 1980--2018. We find the nominal orbit exhibits an unexplained shift by $+60^\circ$ in the true longitude. Using both periodogram analysis and an $\ell = 10$ multipole model suitable for the motion of mutually interacting moons about the irregular body, we confirmed that it is not possible to adjust the respective osculating period $P_2$. Instead, we were forced to use a model with tidal dissipation (and increasing orbital periods) to explain the shift. We also analyzed light curves, spanning 1977--2021, and searched for the expected spin deceleration of Kleopatra. According to our best-fit model, the observed period rate is $\dot P_2 = (1.8\pm 0.1)\cdot 10^{-8}\,{\rm d}\,{\rm d}^{-1}$ and the corresponding time lag $\Delta t_2 = 42\,{\rm s}$ of tides, for the assumed value of the Love number $k_2 = 0.3$. It is the first detection of tidal evolution for moons orbiting 100-km asteroids. The corresponding dissipation factor $Q$ is comparable with other terrestrial bodies, albeit at a higher loading frequency $2|\omega-n|$. We also predict a secular evolution of the inner moon, $\dot P_1 = 5.0\cdot 10^{-8}$, as well as a spin deceleration of Kleopatra, $\dot P_0 = 1.9\cdot 10^{-12}$. In alternative models, with moons captured in the 3:2 mean-motion resonance or more massive moons, the respective values of $\Delta t_2$ are a factor of 2--3 lower. Future astrometric observations by direct imaging or occultations should allow to distinguish between these models, which is important for the internal structure and mechanical properties of (216) Kleopatra.

To figure out the effect of stellar mass and local environment on morphological transformation and star formation quenching in galaxies, we use the massive ($M_* \geq 10^{10} M_{\odot}$) galaxies at $0.5 \leq z \leq 2.5$ in five fields of 3D-HST/CANDELS. Based on the {\it UVJ} diagnosis and the possibility of possessing spheroid, our sample of massive galaxies are classified into four populations: quiescent early-type galaxies (qEs), quiescent late-type galaxies (qLs), star-forming early-type galaxies (sEs), and star-forming late-type galaxies (sLs). It is found that the quiescent fraction is significantly elevated at the high ends of mass and local environmental overdensity, which suggests a clear dependence of quenching on both mass and local environment. Over cosmic time, the mass dependence of galaxy quiescence decreases while the local environment dependence increases. The early-type fraction is found to be larger only at high-mass end, indicating a evident mass dependence of morphological transformation. This mass dependence becomes more significant at lower redshifts. Among the four populations, the fraction of active galactic nucleus (AGN) in the qLs peaks at $2<z \leq 2.5$, and rapidly declines with cosmic time. The sEs are found to have higher AGN fractions of $20-30\%$ at $0.5\leq z<2$ . The redshift evolution of AGN fractions in the qLs and sEs suggests that the AGN feedback could have played important roles in the formation of the qLs and sEs.

The Saturn Kilometric Radiation (SKR) was observed for the first time during the flyby of Saturn by the Voyager spacecraft in 1980. These radio emissions, in the range of a few kHz to 1 MHz, are emitted by electrons travelling around auroral magnetic field lines. Their study is useful to understand the variability of a magnetosphere and its coupling with the solar wind. Previous studies have shown a strong correlation between the solar wind pressure and the SKR intensity. However, up to now, the effect of an Interplanetary Coronal Mass Ejection (ICME) has never been examined in detail, due to the lack of SKR observations at the time when an ICME can be tracked and its different parts be clearly identified. In this study, we take advantage of a large ICME that reached Saturn mid-November 2014 (Witasse et al., 2017). At that time, the Cassini spacecraft was fortunately travelling within the solar wind for a few days, and provided a very accurate timing of the ICME structure. A survey of the Cassini data for the same period indicated a significant increase in the SKR emissions, showing a good correlation after the passage of the ICME shock with a delay of 13 hours and after the magnetic cloud passage with a delay of 25-42 hours. In between, a smaller SKR burst could be correlated with a proton flux peak occurring during the passage of the ejecta of the ICME.

We present a search for disturbed, candidate ram pressure stripping galaxies across more than 50 spectroscopically selected SDSS groups and clusters. Forty-eight ram pressure candidates are visually identified in these systems using high quality UNIONS imaging from the Canada-France Hawaii Telescope, covering ~6200 and ~2800 square degrees in the u- and r-bands respectively. Ram pressure candidates are found in groups and clusters spanning a wide range in halo mass and include ~30 ram pressure candidates in the group regime ($M_h < 10^{14}$). The observed frequency of ram pressure candidates shows substantial scatter with group/cluster mass, but on average is larger in clusters ($M_h > 10^{14}\,M_\odot$) than groups ($M_h < 10^{14}\,M_\odot$) by a factor of ~2. We find that ram pressure candidates are most commonly low-mass galaxies and have enhanced star formation rates relative to star-forming field galaxies. The enhancement in star formation is largely independent of galaxy mass and strongest for galaxies in clusters. As a result of the large survey footprint and excellent image quality from UNIONS, we are able to identify disturbed galaxies, potentially affected by ram pressure stripping, across a wide range of host environment.

Ram pressure stripping is a crucial evolutionary driver for cluster galaxies. It is thought to be able to accelerate the evolution of their star formation, trigger the activity of their central active galactic nucleus (AGN) and the interplay between the galactic and environmental gas, and eventually dissipate their gas reservoir. We explored the outcomes of ram pressure stripping by studying the non-thermal radio emission of the jellyfish galaxy JW100 in the cluster Abell 2626 ($z=0.055$) by combining LOFAR, MeerKAT, and VLA observations from 0.144 to 5.5 GHz. We studied the integrated spectra of the stellar disk, the stripped tail and the AGN, mapped the spectral index over the galaxy, and constrained the magnetic field intensity to be between 11 and 18 $\mu$G in the disk and $<10$ $\mu$G in the tail. The stellar disk radio emission is dominated by a radiatively old plasma, likely related to an older phase of high star formation rate. This suggests that the star formation was quickly quenched by a factor of 4 in a few $10^7$ yr. The radio emission in the tail is consistent with the stripping scenario, where the radio plasma originally accelerated in the disk is then displaced in the tail. The morphology of the radio and X-ray emissions supports the scenario of accretion of the magnetized environmental plasma onto the galaxy. The AGN non-thermal spectrum indicates that the relativistic electron acceleration may have occurred simultaneously with a central ionized gas outflow, thus suggesting a physical connection between the two processes.

Investigating the spatial correlation between different emissions in an extended astrophysical source can provide crucial insights into their physical connection, hence it can be the key to understand the nature of the system. The point-to-point analysis of surface brightness is a reliable method to do such an analysis. In this work, we present PT-REX, a software to carry out these studies between radio and X-ray emission in extended sources. We discuss how to reliably carry out this analysis and its limitation and we introduce the Monte Carlo point-to-point analysis, which allows to extend this approach to poorly-resolved sources. Finally we present and discuss the application of our tool to study the diffuse radio emission in a galaxy cluster.

We propose a method for diagnosing the physical conditions in the solar atmosphere using an increase in oscillation amplitudes resulting from minuscule solar flares. As an example, we consider a B2 flare, which caused a sharp short-lived increase in the amplitude of three- and five-minute oscillations in the lower layers of the solar atmosphere. Enhanced three- and five-minute oscillations propagated from the lower layers of the atmosphere into the corona. Such short oscillation trains made it possible to remove the uncertainties arising in the measurements of the phase and group lags between the layers. In addition, the amplification of the oscillations that reach the corona may add to the likelihood of a repeated flare. Studying oscillations in small flare events has the advantage of exploring the atmosphere in its quasi-quiet condition as opposed to powerful flares, which cause substantial and prolonged disturbance of the environment. In addition, small flares are much more common than powerful flares, which allows one to choose from a larger sample of observational material.

Astronomers have typically set out to solve supervised machine learning problems by creating their own representations from scratch. We show that deep learning models trained to answer every Galaxy Zoo DECaLS question learn meaningful semantic representations of galaxies that are useful for new tasks on which the models were never trained. We exploit these representations to outperform existing approaches at several practical tasks crucial for investigating large galaxy samples. The first task is identifying galaxies of similar morphology to a query galaxy. Given a single galaxy assigned a free text tag by humans (e.g. `#diffuse'), we can find galaxies matching that tag for most tags. The second task is identifying the most interesting anomalies to a particular researcher. Our approach is 100\% accurate at identifying the most interesting 100 anomalies (as judged by Galaxy Zoo 2 volunteers). The third task is adapting a model to solve a new task using only a small number of newly-labelled galaxies. Models fine-tuned from our representation are better able to identify ring galaxies than models fine-tuned from terrestrial images (ImageNet) or trained from scratch. We solve each task with very few new labels; either one (for the similarity search) or several hundred (for anomaly detection or fine-tuning). This challenges the longstanding view that deep supervised methods require new large labelled datasets for practical use in astronomy. To help the community benefit from our pretrained models, we release our fine-tuning code zoobot. Zoobot is accessible to researchers with no prior experience in deep learning.

We developed a GPU based single-pulse search pipeline (GSP) with candidate-archiving database. Largely based upon the infrastructure of Open source pulsar search and analysis toolkit (PRESTO), GSP implements GPU acceleration of the de-dispersion and integrates a candidate-archiving database. We applied GSP to the data streams from the commensal radio astronomy FAST survey (CRAFTS), which resulted in a quasi-real-time processing. The integrated candidate database facilitates synergistic usage of multiple machine-learning tools and thus improves efficient identification of radio pulsars such as rotating radio transients (RRATs) and Fast Radio Bursts (FRBs). We first tested GSP on pilot CRAFTS observations with the FAST Ultra-Wide Band (UWB) receiver. GSP detected all pulsars known from the the Parkes multibeam pulsar survey in the respective sky area covered by the FAST-UWB. GSP also discovered 13 new pulsars. We measured the computational efficiency of GSP to be ~120 times faster than the original PRESTO and ~60 times faster than a MPI-parallelized version of PRESTO.

Recent observations from Parker Solar Probe have revealed that the solar wind has a highly variable structure. How this complex behaviour is formed in the solar corona is not yet known, since it requires omnipresent fluctuations, which constantly emit material to feed the wind. In this article we analysed 14 upflow regions in the solar corona to find potential sources for plasma flow. The upflow regions were derived from spectroscopic data from the EUV Imaging Spectrometer (EIS) onboard Hinode determining their Doppler velocity and defining regions which have blueshifts stronger than $-6\,km\,s^{-1}$. To identify the sources of this blueshift data from the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI), onboard the Solar Dynamics Observatory (SDO), and the X-ray Telescope (XRT), onboard Hinode, were used. The analysis revealed that only 5 out of 14 of the upflows were associated with frequent transients, like obvious jets or bright points. In contrast to that, seven events were associated with small-scale features, which show a large variety of dynamics. Some resemble small bright points, while others show an eruptive nature, all of which are faint and only live for a few minutes; we can not rule out that several of these sources may be fainter and, hence, less obvious jets. Since the complex structure of the solar wind is known, this suggests that new sources have to be considered or better methods used to analyse the known sources. This work shows that small and frequent features, which were previously neglected, can cause strong upflows in the solar corona. These results emphasise the importance of the first observations from the Extreme-Ultraviolet Imager (EUI) onboard Solar Orbiter, which revealed complex small-scale coronal structures.

In 2002, V838 Monocerotis (V838 Mon) erupted in a red novae event which has been interpreted to be a stellar merger. Soon after reaching peak luminosity, it began to cool, and its spectrum evolved to later spectral types. Dust was also formed in the post-merger remnant, making it bright in the mid-infrared. Interferometric studies at these wavelengths have suggested the presence of a flattened, elongated structure. We investigate, for the first time, the structure and orientation of the dusty envelope surrounding V838 Mon in the $L$(2.8-4.2 $\mu$m) band using recent observations with the MATISSE instrument at the VLTI. We perform simple geometrical modelling of the interferometric observables using basic models (disks, Gaussians, point sources, along with their combinations). We also reconstructed an image and analyzed the corresponding $L$-band spectrum. This study indicates the presence of an elongated, disk-like structure near 3.5$\mu$m, similar to what has been observed in other wavelength regimes. In particular, the orientation at a position angle of -40 degrees agrees with prior measurements in other bands. The dusty elongated structure surrounding V838 Mon appears to be a stable and long lived feature that has been present in the system for over a decade. Its substructure and origin remain unclear, but may be related to mass loss phenomena that took place in the orbital plane of the merged binary.

The purpose of this work is to characterize the diffuse Galactic polarized synchrotron, which is the dominant CMB foreground emission at low frequency. We present EE, BB, and EB power spectra estimated from polarization frequency maps at 23 and 30 GHz as observed respectively by the WMAP K-band and the Planck lowest frequency channel, for a set of six sky regions covering from 30% to 94% of the sky. We study the synchrotron polarization angular distribution and spectral energy distribution (SED) by means of the so-called pseudo-$C_\ell$ formalism, provided by the NaMaster package, in the multipole interval 30 $\leq$ $\ell$ $\leq$ 300. Best results are obtained cross-correlating Planck and WMAP data. The EE and BB angular power spectra show a steep decay of the spectral amplitude as a function of multipole, approximated by a power law $C^{EE,BB} \propto \ell^{\alpha_{EE,BB}}$, with $\alpha_{EE}= -2.79\pm0.05$ and $\alpha_{BB}=-2.77\pm0.15$. The B/E power asymmetry is proved with a B-to-E ratio, computed as the amplitude ratio at the pivot multipole $\ell = 80$, of 0.22$\pm$0.02. The EB cross-component is compatible with zero at 1$\sigma$, with upper constraint on the EB/EE ratio of 1.8%. We show that the EE and BB power-law model with null EB cross-correlation describes reasonably well the diffuse synchrotron polarization emission for the full sky if the bright Galactic center and point sources are masked. The recovered SED shows power-law spectral indices $\beta_{EE}$ and $\beta_{BB}$ compatible between themselves, with a weighted average value, in the frequency range 23-30 GHz, of -2.96$\pm$0.09. Results also seem to indicate that the SED gets steeper from low to high Galactic latitude.

In this paper, the recently proposed mass dimension one fermionic field is supposed to be responsible for the dark matter halo around galactic nuclei, through the quantum degeneracy pressure effect of the field. It will be showed that the mass-ratio relation for dwarf galaxies can be well explained for a particle dark matter mass of about $100 - 200$eV. For a large galaxy, as Milky Way, the observational data for rotation curve can be well reproduced for a particle mass of about 23eV, with the addition of other substructures.

We derive an exact, time-dependent analytical magnetic field solution for the inner heliosheath, which satisfies both the induction equation of ideal magnetohydrodynamics in the limit of infinite electric conductivity and the magnetic divergence constraint. To this end, we assume that the magnetic field is frozen into a plasma flow resembling the characteristic interaction of the solar wind with the local interstellar medium. Furthermore, we make use of the ideal Ohm's law for the magnetic vector potential and the electric scalar potential. By employing a suitable gauge condition that relates the potentials and working with a characteristic coordinate representation, we thus obtain an inhomogeneous first-order system of ordinary differential equations for the magnetic vector potential. Then, using the general solution of this system, we compute the magnetic field via the magnetic curl relation. Finally, we analyze the well-posedness of the corresponding Dirichlet boundary value problem, specify compatibility conditions for the boundary values, and outline the implementation of boundary conditions.

Studying lunar dust is vital to the exploration of the Moon and other airless planetary bodies. The Ultraviolet and Visible Spectrometer (UVS) on board the Lunar Atmosphere and Dust Environment Explorer (LADEE) spacecraft conducted a series of Almost Limb activities to look for dust near the dawn terminator region. During these activities the instrument stared at a fixed point in the zodiacal background off the Moon's limb while the spacecraft moved in retrograde orbit from the sunlit to the unlit side of the Moon. The spectra obtained from these activities probe altitudes within a few kilometers of the Moon's surface, a region whose dust populations were not well constrained by previous remote-sensing observations from orbiting spacecraft. Filtering these spectra to remove a varying instrumental signal enables constraints to be placed on potential signals from a dust atmosphere. These filtered spectra are compared with those predicted for dust atmospheres with various exponential scale heights and particle size distributions to yield upper limits on the dust number density for these potential populations. For a differential size distribution proportional to $s^{-3}$ (where $s$ is the particle size) and a scale height of 1 km, we obtain an upper limit on the number density of dust particles at the Moon's surface of 142 $m^{-3}$.

Small reconnection events in the lower solar atmosphere can lead to its heating, but whether such heating can propagate into higher atmospheric layers and potentially contribute to coronal heating is an open question. We carry out a large statistical analysis of all IRIS observations from 2013 and 2014. We identified "IRIS burst" (IB) spectra via a k-means analysis by classifying and selecting Si IV spectra with superimposed blend lines on top of bursts, which indicate low atmospheric heating. We found that ~8% of all observations show IBs with about 0.01% of all recorded IRIS spectra being IB spectra. We found varying blend absorption levels, which may indicate different depths of the reconnection event and heating. IBs are statistically visible with similar properties and timings in the spectral lines Mg II, C II, and Si IV, but invisible in Fe XXI. By statistically analyzing co-spatial AIA lightcurves, we found systematic enhancements in AIA 1600 and AIA 1700, but no clear response to bursts in all other AIA wavelengths (94, 131, 171, 193, 211, 304, 335) in a timeframe of $\pm 6$ minutes around the burst. This may indicate that heating due to IBs is confined within the lower atmosphere and dissipates before reaching temperatures or formation heights covered by the hotter AIA lines. Our developed methods are applicable for statistical analyses of any co-observed data sets and allow us to efficiently analyze millions of spectra and lightcurves simultaneously.

We investigate neutron star-black hole (NS-BH) merger candidates as a test for compact exotic objects. Using the events GW190814, GW200105 and GW200115 measured by the LIGO-Virgo collabration, which represent a broad profile of the masses in the NS mass spectrum, we demonstrate the constraining power for the parameter spaces of compact stars consisting of dark matter for future measurements. We consider three possible cases of dark matter stars: self-interacting, purely bosonic or fermionic dark matter stars, stars consisting of a mixture of interacting bosonic and fermionic matter, as well as the limiting case of selfbound stars. We find that the scale of those hypothetical objects are dominated by the one of the strong interaction. The presence of fermionic dark matter requires a dark matter particle of the GeV mass scale, while the bosonic dark matter particle mass can be arbitrarily large or small. In the limiting case of a selfbound linear equation of state, we find that the vacuum energy of those configurations has to be similar to the one of QCD.

The coevolution between supermassive black holes (SMBHs) and their host galaxies has been proposed for more than a decade, albeit with little direct evidence about black hole accretion activities regulating galaxy star formation at $z>1$. In this paper, we study the lifetimes of X-ray active galactic nuclei (AGNs) in $UV$-selected red sequence (RS), blue cloud (BC) and green valley (GV) galaxies, finding that AGN accretion activities are most prominent in GV galaxies at $z\sim1.5-2$, compared with RS and BC galaxies. We also compare AGN accretion timescales with typical color transition timescales of $UV$-selected galaxies. We find that the lifetime of GV galaxies at $z\sim1.5-2$ is very close to the typical timescale when the AGNs residing in them stay in the high-accretion-rate mode at these redshifts; for BC galaxies, the consistency between the color transition timescale and the black hole strong accretion lifetime is more likely to happen at lower redshifts ($z<1$). Our results support the scenario where AGN accretion activities govern $UV$ color transitions of host galaxies, making galaxies and their central SMBHs coevolve with each other.

Apart from regular, low-level stochastic variability, some AGN occasionally show exceptionally large changes in the luminosity, spectral shape, and/or X-ray absorption. The most notable are the changes of the spectral type when the source classified as a Seyfert 1 becomes a Seyfert 2 galaxy or vice versa. Thus a name was coined of 'Changing-Look AGN' (CL AGN). The origin of this phenomenon is still unknown, but for most of the sources, there are strong arguments in favor of the intrinsic changes. Understanding the nature of such rapid changes is a challenge to the models of black hole accretion flows since the timescales of the changes are much shorter than the standard disk viscous timescales. We aim to model the CL AGN phenomenon assuming that the underlying mechanism is the time-dependent evolution of a black hole accretion disk unstable due to the dominant radiation pressure. We use a 1-dimensional, vertically integrated disk model, but we allow for the presence of the hot coronal layer above the disk and the presence of the inner purely hot flow. We focus on the variability timescales and amplitudes, which can be regulated by the action of large-scale magnetic fields, the description of the disk-corona coupling, and the presence of an inner optically thin flow, like Advection-Dominated Accretion Flow (ADAF). We compare model predictions for the accretion disk around black hole mass 10$^7$M$_{\odot}$.

A fully-sampled and hitherto highest resolution and sensitivity observation of neutral hydrogen (HI) in the Leo Triplet (NGC 3628, M 65/NGC 3623, and M 66/NGC 3627) reveals six HI structures beyond the three galaxies. We present detailed results of the morphologies and kinematics of these structures, which can be used for future simulations. In particular, we detect a two-arm structure in the plume of NGC 3628 for the first time, which can be explained by a tidal interaction model. The optical counterpart of the plume is mainly associated with the southern arm. The connecting part (base) of the plume (directed eastwards) with NGC 3628 is located at the blueshifted (western) side of NGC 3628. Two bases appear to be associated with the two arms of the plume. A clump with reversed velocity gradient (relative to the velocity gradient of M 66) and a newly detected tail, i.e. M 66SE, is found in the southeast of M 66. We suspect that M 66SE represents gas from NGC 3628 which was captured by M 66 in the recent interaction between the two galaxies. Meanwhile gas is falling toward M 66, resulting in features already previously observed in the southeastern part of M 66, e.g. large line widths and double peaks. An upside-down `Y'-shaped HI gas component (M 65S) is detected in the south of M 65 which suggests that M 65 may also have been involved in the interaction. We strongly encourage modern hydrodynamical simulations of this interacting group of galaxies to reveal the origin of the gaseous debris surrounding all three galaxies.

The study of Gamma Ray Bursts (GRBs) has the potential to improve our understanding of high energy astrophysical phenomena. In order to reliably use GRBs to this end, we first need to have a well-developed grasp of the mechanism that produces the radiation within GRB jets and how that relates to their structure. One model for the emission mechanism of GRBs invokes radiation produced deep in the jet which eventually escapes the jet at its photosphere. While this model has been able to explain a number of observed GRB characteristics, it is currently lacking in predictive power and in ability to fully reproduce GRB spectra. In order to address these shortcomings of the model, we have expanded the capabilities of the MCRaT code, a state of the art radiative transfer code that can now simulate optical to gamma ray radiation propagating in a hydrodynamically simulated GRB jet. Using the MCRaT code, we have constructed mock observed light curves, spectra, and polarization from optical to gamma ray energies for the simulated GRBs. Using these mock observables, we have compared our simulations of photospheric emission to observations and found much agreement between the two. Furthermore, the MCRaT calculations combined with the hydrodynamical simulations allow us to connect the mock observables to the structure of the simulated GRB jet in a way that was not previously possible. While there are a number of improvements that can be made to the analyses, the steps taken here begin to pave the way for us to fully understand the connection between the structure of a given GRB jet and the radiation that would be expected from it.

The demographics of dark matter substructure depend sensitively on the nature of dark matter. Optimally leveraging this probe requires accurate theoretical predictions regarding the abundance of subhaloes. These predictions are hampered by artificial disruption in numerical simulations, by large halo-to-halo variance, and by the fact that the results depend on the baryonic physics of galaxy formation. In particular, numerical simulations have shown that the formation of a central disc can drastically reduce the abundance of substructure compared to a dark matter-only simulation, which has been attributed to enhanced destruction of substructure due to disc shocking. We examine the impact of discs on substructure using the semi-analytical subhalo model SatGen, which accurately models the tidal evolution of substructure free of the numerical disruption that still hampers $N$-body simulations. Using a sample of 10,000 merger trees of Milky Way-like haloes, we study the demographics of subhaloes that are evolved under a range of composite halo-disc potentials with unprecedented statistical power. We find that the overall subhalo abundance is relatively insensitive to properties of the disc aside from its total mass. For a disc that contains $5\%$ of $M_\mathrm{vir}$, the mean subhalo abundance within $r_\mathrm{vir}$ is suppressed by ${\lesssim}10\%$ relative to the no-disc case, a difference that is dwarfed by halo-to-halo variance. For the same disc mass, the abundance of subhaloes within 50 kpc is reduced by ${\sim}30\%$. We argue that the disc mainly drives excess mass loss for subhaloes with small pericentric radii and that the impact of disc shocking is negligible.

We present near infrared spectroscopy of Nova Herculis 2021 (V1674 Her), obtained over the first 70 days of its evolution. This fastest nova on record displays a rich emission line spectrum, including strong coronal line emission with complex structures. The hydrogen line fluxes, combined with a distance of 4.7 (+1.3 / -1.0) kpc, give an upper limit to the hydrogen ejected mass of 1.4 (+0.8 / -1.2) 10^{-3} solar masses. The coronal lines appeared at day 11.5, the earliest onset yet observed for any classical nova, before there was an obvious source of ionizing radiation. We argue that the gas cannot be photoionized, at least in the earliest phase, and must be shocked. Its temperature is estimated to be 10^{5.57 +/- 0.05} K on day 11.5. Tentative analysis indicates a solar abundance of aluminum and an underabundance of calcium, relative to silicon, with respect to solar values in the ejecta. Further, we show that the vexing problem of whether collisional or photoionization is responsible for coronal emission in classical novae can be resolved by correlating the temporal sequence in which the X-ray supersoft phase and the near-infrared coronal line emission appear.

We present the discovery and characterisation of two transiting planets discovered by TESS in the light curves of the young and bright (V=9.67) star HD 73583 (TOI-560). We perform an intensive spectroscopic and photometric space- and ground-based follow-up in order to confirm and characterise the system. We found that HD 73583 is a young (750 Myr) active star with a rotational period of $12.1$ d and a mass and radius of $ 0.71 M_\odot$ and $0.66 R_\odot$, respectively. HD 73583 b ($P_{\rm b}=$ $6.4$d) has a mass and radius of $10 M_\oplus$ and $2.83 \pm 0.10 R_\oplus$, respectively, that gives a density of $2.43 {\rm g\,cm^{-3}}$. HD 73583 c ($P_{\rm c}= 18.9$ d) has a mass and radius of $9.6 M_\oplus$ and $2.37 R_\oplus$, respectively, this translates to a density of $3.97 {\rm g\,cm^{-3}}$. Both planets are consistent with worlds made of a solid core surrounded by a thick envelope. Because of their youth and host star brightness, they both are excellent candidates to perform transmission spectroscopy studies and search for mass-loss signatures. We expect atmospheric mass-loss rates of $2.4 \times 10^{10}\,{\rm g\,s^{-1}}$ and $5.4 \times 10^{9}\,{\rm g\,s^{-1}}$ for HD73583 b and c, respectively. We expect that the detection of evaporating signatures on H and He would be challenging, but doable with present and future instruments.

Our ability to explore the cosmos by direct contact has been limited to a small number of lunar and interplanetary missions. However, the NASA Starlight program points a path forward to send small, relativistic spacecraft far outside our solar system via standoff directed-energy propulsion. These miniaturized spacecraft are capable of robotic exploration but can also transport seeds and organisms, marking a profound change in our ability to both characterize and expand the reach of known life. Here we explore the biological and technological challenges of interstellar space biology, focusing on radiation-tolerant microorganisms capable of cryptobiosis. Additionally, we discuss planetary protection concerns and other ethical considerations of sending life to the stars.

We present light curve analyses of two newly identified detached eclipsing binaries, HD 96609 and HD 303734, in the region of the richly populated open cluster NGC 3532. HD 96609 is composed of two main sequence stars (B9-A0V + A2V) with masses and radii of $M_{1}=2.66\pm0.02$$M_{\odot}$, $M_{2}=1.84\pm0.01$$M_{\odot}$, $R_{1}=2.740\pm0.006$$R_{\odot}$, $R_{2}=1.697\pm0.005$$R_{\odot}$. The positions of the components on $log~M-log~R$ plane suggests log(age/yr) 8.55, corresponding $350\pm40$ Myr of age, which agrees with the $300\pm100$ Myr age of NGC 3532 estimated in previous studies. We find the distance of HD 96609 as $460\pm17$ pc, which is consistent with the $484^{+35}_{-30}$ pc distance of NGC 3532, estimated from GAIA parallaxes. HD 303734 is an interesting totally eclipsing binary with a quite shallow secondary eclipse. Using photometric properties of the system in conjunction with theoretical calibrations, we estimate that HD 303734 consists of A6V + K3V components. HD 96609 and HD 303734 are the second and third eclipsing binaries discovered in the region of NGC 3532, after the first one, GV Car.

In this paper we present a sample of variability-selected AGN from a parent sample of dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward-modeled mid-IR photometry of time-resolved Wide-field Infrared Survey Explorer (WISE) coadded images. We found that 44 out of 25,714 dwarf galaxies had optically variable AGN candidates, and 158 out of 79,879 dwarf galaxies had mid-IR variable AGN candidates, corresponding to active fractions of $0.17\pm0.03$% and $0.20\pm0.02$% respectively. Only two objects, NSA164884 and NSA451469, broad line AGN with virial masses $M_{\text{BH}}=10^{6.9}M_\odot$ and $M_{\text{BH}}=10^{6.3}M_\odot$ respectively, were found in both optical and mid-IR searches. We find that spectroscopic approaches to AGN identification would have missed 81% of our ZTF IMBH candidates and 69% of our WISE IMBH candidates. Only $9$ candidates have been detected previously in radio, X-ray, and variability searches for dwarf galaxy AGN. The IMBHs with broad Balmer lines have virial masses of $10^5M_\odot<M_{\text{BH}}<10^7M_\odot$ but for the rest of the sample, BH masses predicted from host galaxy mass range between $10^{4.8}M_\odot<M_{\text{BH}}<10^7M_\odot$. We found that only 5 of 152 previously reported variability-selected AGN candidates from the Palomar Transient Factory in common with our parent sample were variable in ZTF. We also determined a nuclear supernova fraction of $0.05\pm0.01$% year$^{-1}$ for dwarf galaxies in ZTF. Our ZTF and WISE IMBH candidates show the promise of variability searches for discovery of otherwise hidden low mass AGN in preparation for the Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory.

Searching for gravitational waves from compact binary coalescences (CBC) is performed by matched filtering the observed strain data from gravitational-wave observatories against a discrete set of waveform templates designed to accurately approximate the expected gravitational-wave signal, and are chosen to efficiently cover a target search region. The computational cost of matched filtering scales with both the number of templates required to cover a parameter space and the in-band duration of the waveform. Both of these factors increase in difficulty as the current observatories improve in sensitivity, especially at low frequencies, and may pose challenges for third-generation observatories. Reducing the cost of matched filtering would make searches of future detector's data more tractable. In addition, it would be easier to conduct searches that incorporate the effects of eccentricity, precession or target light sources (e.g. subsolar). We present a hierarchical scheme based on a reduced bases method to decrease the computational cost of conducting a matched-filter based search. Compared to the current methods, we estimate without any loss in sensitivity, a speedup by a factor of $\sim$ 18 for sources with signal-to-noise ratio (SNR) of at least $= 6.0$, and a factor of $8$ for SNR of at least 5. Our method is dominated by linear operations which are highly parallelizable. Therefore, we implement our algorithm using graphical processing units (GPUs) and evaluate commercially motivated metrics to demonstrate the efficiency of GPUs in CBC searches. Our scheme can be extended to generic CBC searches and allows for efficient matched filtering using GPUs.

Recent estimates of the characteristics of Planet Nine have suggested that it could be closer than originally assumed. Such a Planet Nine would also be brighter than originally assumed, suggesting the possibility that it has already been observed in wide-field moderate-depth surveys. We search for Planet Nine in the Zwicky Transient Facility public archive and find no candidates. Using known asteroids to calculate the magnitude limit of the survey, we find that we should have detected Planet Nine throughout most of the northern portion of its predicted orbit -- including within the galactic plane -- to a 95% detection efficiency of approximately $V=20.5$. To aid in understanding detection limits for this and future analyses, we present a full-sky synthetic Planet Nine population drawn from a statistical sampling of predicted Planet Nine orbits. We use this reference population to estimate that this survey rules out 56% of predicted Planet Nine phase space, and we demonstrate how future analyses can use the same synthetic population to continue to constrain the amount of parameter space effectively searched for Planet Nine.

Extreme TeV blazars (ETBs) are active galactic nuclei with jets presumably pointing towards the observer having their intrinsic spectral energy distributions (SEDs) peaked at an energy in excess of 1 TeV. These sources typically reveal relatively weak and slow variability as well as an extremely high frequency of the low-energy SED peak compared to other classes of blazars. It proved to be exceedingly hard to incorporate all these peculiar properties of ETBs into the framework of a reasonable $\gamma$-ray emission model. ETB physics have recently attracted great attention in the astrophysical community, underlying the importance of the development of self-consistent ETB emission model(s). We propose a new scenario for the formation of X-ray and $\gamma$-ray spectra of ETBs assuming that electromagnetic cascades develop in the infrared photon field surrounding the central blazar engine. This scenario does not invoke compact fast-moving sources of radiation (so-called "blobs"), in agreement with the apparent absence of fast and strong variability of ETBs. For the case of the extreme TeV blazar 1ES 0229+200 we propose a specific emission model in the framework of the considered scenario. We demonstrate that this model allows to obtain a good fit to the measured SED of 1ES 0229+200.

We present a closed-form normalization method suitable for the study of the secular dynamics of small bodies inside the trajectory of Jupiter. The method is based on a convenient use of a book-keeping parameter introduced not only in the Lie series organization but also in the Poisson bracket structure employed in all perturbative steps. In particular, we show how the above scheme leads to a redefinition of the remainder of the normal form at every step of the formal solution of the homological equation. An application is given for the semi-analytical representation of the orbits of main-belt asteroids.

A unique characteristic of exponentially growing scattering amplitude arises in an anomalous Abelian effective field theory when an extremely light Dirac neutrino mass is introduced to break the symmetry. We show that the low energy effective Lagrangian can be made explicitly gauge invariant with the help of a nonlinear representation of the Goldstone or Stueckelberg field. We study the peculiar feature of exponential growth in the ultra-high-energy neutrino-nucleon inelastic scattering. It is found that the inelastic scattering cross section is highly sensitive to the ratio of gauge coupling to the gauge boson mass, $g_X/m_X$. When the IceCube measurement of ultra-high-energy neutrinos, which is consistent with the standard model prediction up to $E_\nu \sim 6$ PeV, is taken into account, the inferred constraint on $g_X/m_X$ is more severe than that obtained from the events of mono-lepton$+$missing transverse energy at the LHC. A muon collider with a collision energy of $10$ TeV can be a good environment other than hadron colliders to probe the novel effect.

Recent observational data on transiently-accreting neutron stars has unequivocally shown fast-cooling sources, such as in the case of neutron star MXB 1659-29. Previous calculations have estimated its total neutrino luminosity and heat capacity, as well as suggested that direct Urca reactions take place in $1 \%$ of the volume of the core. In this paper, we reproduce the inferred luminosity of this source with detailed models of equations of state (EOS) and nuclear pairing gaps. We show that three superfluidity gap models are inconsistent with data for all EOS and another three are disfavoured because of fine tuning arguments. We also calculate the total heat capacity for all constructed stars and show that independent observations of mass and luminosity could set constraints on the core superfluidity of a source as well as the density slope of the symmetry energy, L. This is an important step towards defining a universal equation of state for neutron stars and therefore, towards a better understanding of the phase diagram of asymmetric matter at high densities.

The cosmological constant $\Lambda$ is a measure of the energy density of the vacuum. Therefore properties of the energy of the system in the metastable vacuum state reflect properties of $\Lambda = \Lambda(t)$. We analyze properties of the energy, $E(t)$, of a general quantum system in the metastable state in various phases of the decay process: In the exponential phase, in the transition phase between the exponential decay and the later phase, where decay law as a function of time $t$ is in the form of powers of $1/t$, and also in this last phase. We found that this energy having an approximate value resulting from the Weisskopf--Wigner theory in the exponential decay phase is reduced very fast in the transition phase to its asymptotic value $E(t) \simeq E_{min} + \alpha_{2}/t^{2}+\ldots$ in the late last phase of the decay process. (Here $E_{min}$ is the minimal energy of the system). This quantum mechanism reduces the energy of the system in the unstable state by a dozen or even several dozen orders or more. We show that if to assume that a universe was born in metastable false vacuum state then according to this quantum mechanism the cosmological constant $\Lambda$ can have a very great value resulting from the quantum field theory calculations in the early universe in the inflationary era, $\Lambda \simeq \Lambda_{qft}$, and then it can later be quickly reduced to the very, very small values.

Axion-like particles (ALPs) are undiscovered pseudo-scalar particles that are candidates for ultralight dark matter. ALPs interact with photons slightly and cause the rotational oscillation of linearly polarized light. Dark matter Axion search with riNg Cavity Experiment (DANCE) searches for ALP dark matter by amplifying the rotational oscillation with a bow-tie ring cavity. Simultaneous resonance of linear polarizations is necessary to amplify both the carrier field and the ALP signal, and to achieve the design sensitivity. The sensitivity of the current prototype experiment DANCE Act-1 is less than expectation by around three orders of magnitude due to the resonant frequency difference between s- and p-polarization in the bow-tie ring cavity. In order to tune the resonant frequency difference, the method of introducing an auxiliary cavity was proposed. We designed an auxiliary cavity that can cancel out the resonant frequency difference and realize simultaneous resonance, considering optical loss. We also confirmed that the sensitivity of DANCE Act-1 with the auxiliary cavity can reach the original sensitivity.

This paper investigates the application of the Koopman Operator theory to the motion of a satellite about a libration point in the Circular Restricted Three-Body Problem. Recently, the Koopman Operator has emerged as a promising alternative to the geometric perspective for dynamical systems, where the Koopman Operator formulates the analysis and dynamical systems in terms of observables. This paper explores the use of the Koopman Operator for computing both 2D and 3D periodic orbits near libration points. Further, simulation results show that the Koopman Operator provides analytical solutions with high accuracy for both Lyapunov and Halo orbits, which are then applied to a station-keeping application.

We carefully perform a Hamiltonian Dirac's constraint analysis of $\omega=-\frac{3}{2}$ Brans-Dicke theory with Gibbons-Hawking-York (GHY) boundary term. The Poisson brackets are computed via functional derivatives. After a brief summary of the results for $\omega\neq-\frac{3}{2}$ case, we derive all Hamiltonian Dirac's constraints and constraint algebra both in the Jordan and Einstein frames. Confronting and contrasting Dirac's constraint algebra in both frames, it is shown that they are not equivalent. This highlights the transformations from the Jordan to the Einstein frames are not Hamiltonian canonical transformations.

Molecular chirality is inherent to biology and cellular chemistry. In this report, the origin of enantiomeric selectivity is analyzed from the viewpoint of the "RNA World" model, based on the autocatalytic self-replication of glyceraldehyde as a precursor for simple sugars, and in particular ribose, as promoted by the formose reaction. Autocatalytic coupling of formaldehyde and glycolaldehyde produces glyceraldehyde, which contains a chiral carbon center that is carried through in formation of the ribose ring. The parity non-conserving weak interaction is the only inherently handed property in nature and is herein shown to be sufficient to differentiate between two enantiomeric forms in an autocatalytic reaction performed over geologically relevant time scales, but only in the presence of a catalytic metal ion such as divalent calcium or higher Z alkaline earth elements. This work details calculations of the magnitude of the effect, the impact of various geologically-available metal ions, and the influence on evolution and dominance of chirality in the molecules of life.

An alternative scenario about the phenomenology of primordial Universe is k-inflation. According to this concept, inflation can be achieved by nonstandard kinetic term of scalar field, namely the inflaton. In this project we focus on k-essence models in the presence of a higher order and a linear kinetic term. Furthermore, the inflationary phenomenology with a Dirac-Born-Infeld scalar field is briefly examined, which arises from quantum theories of gravity such as superstring theory. Our approach about the inflationary era is that it can be described in the context of Einstein's gravity involving quantum corrections such as the Chern-Simons string inspired parity violating gravitational term. The equations of motion namely, the Friedmann equation, the Raychadhuri equation and the Klein-Gordon equation for an expanding background are extracted from the gravitational action utilizing the variational principle. The consequential system of differential equations with respect to Hubble's parameter and the inflaton field was quite perplexed in order to be solved with an analytic way. Therefore, the slow-roll conditions during inflationary era were imposed and terms with minor numerically contribution were neglected. From the overall phenomenological analysis it is proved that, models with exotic kinetic terms can generate viable results in consistency with the latest Planck data. Finally, the presence of Chern-Simons quantum corrections shifts the primordial spectral tensor index to blue. Even though blue gravitational waves have yet to be observed, if detected, compatibility with the aforementioned theory can be achieved.

The developments of sounding rocket and CubeSat are a total game changer to the space program and it allows building space instruments to be more achievable and affordable. Therefore, it gives us a good opportunity to build a small cosmic ray detector which has capabilities to measure the flux, direction, and even energy of cosmic rays at the height above the limitation of balloon experiments, and it may open a new door for the cosmic ray physics. Compact Scintillator Array Detector (ComSAD) is dedicated for the sounding rocket mission of Taiwan's National Space Organization. In paper, we present the idea, design, and performance of ComSAD which is also suitable for CubeSat missions in the future.

Emission of anisotropic gravitational radiation from compact binary system leads to a flux of linear momentum. This results in the recoil of the system. We investigate the rate of loss of Linear momentum flux in the far zone of the source using various mass type and current type multipole moments for inspiralling compact binary mergers in quasi-elliptical orbits at 2.5 Post Newtonian order. We compute the linear momentum flux accurate up to $\mathcal{O}(e_t)$ in harmonic coordinate. A 2.5 Post Newtonian Quasi-Keplarian representation of the parametric solution to the Post Newtonian equation of motion for the compact binary system has been adopted here. We also provide a closed-form expression for the accumulated linear momentum from the remote past through the binary evolution.

Core-collapse supernovae (CCSN) are a prime source of gravitational waves. Estimations of their typical frequencies make them perfect targets for the current network of advanced, ground-based detectors. A successful detection could potentially reveal the underlying explosion mechanism through the analysis of the waveform. This has been illustrated using the SupernovaModel Evidence Extractor (SMEE; Logue et al. (2012)), an algorithm based on principal-component analysis and Bayesian model selection. Here, we present a complementary approach to SMEE based on (supervised) dictionary-learning and show that it is able to reconstruct and classify CCSN signals according to their morphology. Our waveform signals are obtained from (a) two publicly available catalogs built from numerical simulations of neutrino-driven (Mur) and magneto-rotational (Dim) CCSN explosions and (b) from a third 'mock' catalog of simulated sine-Gaussian (SG) waveforms. Those signals are injected into coloured Gaussian noise to simulate the background noise of Advanced LIGO in its broadband configuration and scaled to a freely-specifiable signal-to-noise ratio (SNR). We show that our approach correctly classifies signals from all three dictionaries. In particular, for SNR=15-20, we obtain perfect matches for both Dim and SG signals and about 85% true classifications for Mur signals. These results are comparable to those reported by SMEE for the same CCSN signals when those are injected in only one LIGO detector. We discuss the main limitations of our approach as well as possible improvements.