Papers for Thursday, Jul 22 2021

Interacting binaries are of general interest as laboratories for investigating the physics of accretion, which gives rise to the bulk of high-energy radiation in the Galaxy. They allow us to probe stellar evolution processes that cannot be studied in single stars. Understanding the orbital evolution of binaries is essential in order to model the formation of compact binaries. Here we focus our attention on studying orbital evolution driven by angular momentum loss through stellar winds in massive binaries. We run a suite of hydrodynamical simulations of binary stars hosting one mass losing star with varying wind velocity, mass ratio, wind velocity profile and adiabatic index, and compare our results to analytic estimates for drag and angular momentum loss. We find that, at leading order, orbital evolution is determined by the wind velocity and the binary mass ratio. Small ratios of wind to orbital velocities and large accreting companion masses result in high angular momentum loss and a shrinking of the orbit. For wider binaries and binaries hosting lighter mass-capturing companions, the wind mass-loss becomes more symmetric, which results in a widening of the orbit. We present a simple analytic formula that can accurately account for angular momentum losses and changes in the orbit, which depends on the wind velocity and mass ratio. As an example of our formalism, we compare the effects of tides and winds in driving the orbital evolution of high mass X-ray binaries, focusing on Vela X-1 and Cygnus X-1 as examples.

We use machine learning techniques to investigate their performance in classifying active galactic nuclei (AGNs), including X-ray selected AGNs (XAGNs), infrared selected AGNs (IRAGNs), and radio selected AGNs (RAGNs). Using known physical parameters in the Cosmic Evolution Survey (COSMOS) field, we are able to well-established training samples in the region of Hyper Suprime-Cam (HSC) survey. We compare several Python packages (e.g., scikit-learn, Keras, and XGBoost), and use XGBoost to identify AGNs and show the performance (e.g., accuracy, precision, recall, F1 score, and AUROC). Our results indicate that the performance is high for bright XAGN and IRAGN host galaxies. The combination of the HSC (optical) information with the Wide-field Infrared Survey Explorer (WISE) band-1 and WISE band-2 (near-infrared) information perform well to identify AGN hosts. For both type-1 (broad-line) XAGNs and type-1 (unobscured) IRAGNs, the performance is very good by using optical to infrared information. These results can apply to the five-band data from the wide regions of the HSC survey, and future all-sky surveys.

Extragalactic planetary nebulae (PNe) offer a way to determine the distance to their host galaxies thanks to the nearly universal shape of the planetary nebulae luminosity function (PNLF). Accurate PNe distance measurements rely on obtaining well-sampled PNLFs and the number of observed PNe scales with the encompassed stellar mass. This means either disposing of wide-field observations or focusing on the bright central regions of galaxies. In this work we take this second approach and conduct a census of the PNe population in the central regions of galaxies in the Fornax cluster, using VLT/MUSE data for the early-type galaxies observed over the course of the Fornax3D survey. Using such integral-field spectroscopic observations to carefully separate the nebular emission from the stellar continuum, we isolated [OIII] 5007 {\AA} sources of interest, filtered out unresolved impostor sources or kinematic outliers, and present a catalogue of 1350 unique PNe sources across 21 early-type galaxies, which includes their positions, [OIII] 5007 {\AA} line magnitudes, and line-of-sight velocities. Using the PNe catalogued within each galaxy, we present independently derived distance estimates based on the fit to the entire observed PNLF observed while carefully accounting for the PNe detection incompleteness. With these individual measurements, we arrive at an average distance to the Fornax cluster itself of 19.86 $\pm$ 0.32 Mpc ($\mu_{PNLF}$ = 31.49 $\pm$ 0.04 mag). Our PNLF distance measurements agree well with previous distances based on surface brightness fluctuations, finding no significant systematic offsets between the two methods as otherwise reported in previous studies.

We present and analyze ALMA submillimeter observations from a multi-wavelength campaign of Sgr A* during 18 July 2019. In addition to the submillimeter, we utilize concurrent mid-IR (Spitzer) and X-ray (Chandra) observations. The submillimeter emission lags $\delta t=21.48^{+3.44}_{-3.57}$ minutes behind the mid-IR data. The entire submillimeter flare was not observed, raising the possibility that the time delay is a consequence of incomplete sampling of the light curve. The decay of the submillimeter emission is not consistent with synchrotron cooling. Therefore, we analyze these data adopting an adiabatically expanding synchrotron source that is initially optically thick or thin in the submillimeter, yielding time-delayed or synchronous flaring with the IR, respectively. The time-delayed model is consistent with a plasma blob of radius $0.8~R_{\text{S}}$ (Schwarzschild radius), electron power-law index $p=3.5$ ($N(E)\propto E^{-p}$), equipartition magnetic field of $B_{\text{eq}}\approx90$ Gauss, and expansion velocity $v_{\text{exp}}\approx0.004c$. The simultaneous emission is fit by a plasma blob of radius $2~R_{\text{S}}$, $p=2.5$, $B_{\text{eq}}\approx27$ Gauss, and $v_{\text{exp}}\approx0.014c$. Since the submillimeter time delay is not completely unambiguous, we cannot definitely conclude which model better represents the data. This observation presents the best evidence for a unified flaring mechanism between submillimeter and X-ray wavelengths and places significant constraints on the source size and magnetic field strength. We show that concurrent observations at lower frequencies would be able to determine if the flaring emission is initially optically thick or thin in the submillimeter.

To date, studies of $\textit{Laplace Surface}$ dynamics have concerned themselves with test particle orbits of fixed shape and orientation in the combined field of an oblate central body (to which the particle is bound) and a distant, inclined, companion which is captured to quadrupolar order. While amply sufficient for satellites around planets on near-circular orbits, the quadrupolar approximation fails to capture essential dynamical features induced by a wide binary companion (be it a star, a planet or a black hole) on a fairly eccentric orbit. With similar such astronomical settings in mind, we extend the classical Laplace framework to higher multipoles, and map out the backbone of stationary orbits, now complexified by the broken axial symmetry. Eccentric and inclined Laplace equilibria, which had been presaged in systems of large enough mutual inclination, are here delineated over a broad range of mutually inclined perturbations. We recover them for test particles in the field of a hot Jupiter and a wide eccentric stellar binary, highlighting their relevance for the architecture of multi-planet systems in binaries. We then extend and deploy our machinery closer to home, as we consider the secular dynamics of Trans-Neptunian Objects (TNOs) in the presence of a putative ninth planet. We show how generalized Laplace equilibria seed islands for Trans-Neptunian objects to be sheltered around, islands within chaotic seas which we capture via Poincar$\'{e}$ sections, while highlighting a beautiful interplay between Laplace and Kozai-Lidov secular dynamical structures. An eminently classical tale revived for the exo-planetary 21st century!

The possible origin of millisecond bursts from the giant elliptical galaxy M87 has been scrutinized since the earliest searches for extragalactic fast radio transients undertaken in the late 1970s. Motivated by rapid technological advancements in recent years, we conducted $\rm \simeq 10~hours$ of L-band ($\rm 1.15-1.75~GHz$) observations of the core of M87 with the Arecibo radio telescope in 2019. Adopting a matched filtering approach, we searched our data for single pulses using trial dispersion measures up to $\rm 5500~pc~cm^{-3}$ and burst durations between $\rm 0.3-123~ms$. We find no evidence of astrophysical bursts in our data above a 7$\sigma$ detection threshold. Our observations thus constrain the burst rate from M87 to $\rm \lesssim 0.1~bursts~hr^{-1}$ above $\rm 1.4~Jy~ms$, the most stringent upper limit obtained to date. Our non-detection of radio bursts is consistent with expectations of giant pulse emission from a Crab-like young neutron star population in M87. However, the dense, strongly magnetized interstellar medium surrounding the central $\sim 10^9 \ M_{\odot}$ supermassive black hole of M87 may potentially harbor magnetars that can emit detectable radio bursts during their flaring states.

Binary black hole spins are among the key observables for gravitational wave astronomy. Among the spin parameters, their orientations within the orbital plane, $\phi_1$, $\phi_2$ and $\Delta \phi=\phi_1-\phi_2$, are critical for understanding the prevalence of the spin-orbit resonances and merger recoils in binary black holes. Unfortunately, these angles are particularly hard to measure using current detectors, LIGO and Virgo. Because the spin directions are not constant for precessing binaries, the traditional approach is to measure the spin components at some reference stage in the waveform evolution, typically the point at which the frequency of the detected signal reaches 20 Hz. However, we find that this is a poor choice for the orbital-plane spin angle measurements. Instead, we propose measuring the spins at a fixed \emph{dimensionless} time or frequency near the merger. This leads to significantly improved measurements for $\phi_1$ and $\phi_2$ for several gravitational wave events. Furthermore, using numerical relativity injections, we demonstrate that $\Delta \phi$ will also be better measured near the merger for louder signals expected in the future. Finally, we show that numerical relativity surrogate models are key for reliably measuring the orbital-plane spin orientations, even at moderate signal-to-noise ratios like $\sim 30-45$.

Binary black hole spin measurements from gravitational wave observations can reveal the binary's evolutionary history. In particular, the spin orientations of the component BHs within the orbital plane, $\phi_1$ and $\phi_2$, can be used to identify binaries caught in the so-called spin-orbit resonances. In a companion paper, we demonstrate that $\phi_1$ and $\phi_2$ are best measured near the merger of the two black holes. In this work, we use these spin measurements to constrain the distribution of $\phi_1$ and $\Delta \phi=\phi_1 - \phi_2$ over the population of merging binary black holes. We find that there is a preference for $\Delta \phi \sim \pm \pi$ in the population, which can be a signature of spin-orbit resonances. We also find a preference for $\phi_1 \sim -\pi/4$ with respect to the line of separation near merger, which has not been predicted for any astrophysical formation channel. However, we are unable to constrain the widths of the $\phi_1$ and $\Delta \phi$ distributions, and more observations may be necessary to confirm these trends. Finally, we derive constraints on the distribution of recoil kicks in the population, and use this to estimate the fraction of merger remnants retained by globular and nuclear star clusters.

We consider topological configurations of the magnetically coupled spinning stellar binaries (e.g., merging neutron stars or interacting star-planet systems). We discuss conditions when the stellar spins and the orbital motion nearly `compensate' each other, leading to very {\it slow} overall winding of the coupled magnetic fields; slowly winding configurations allow gradual accumulation of magnetic energy that is eventually released in a flare when the instability threshold is reached. We find that this slow winding can be global and/or local. We describe the topology of the relevant space $\mathbb{F}=T^1S^2$ as the unit tangent bundle of the two-sphere and find conditions for slowly winding configurations in terms of magnetic moments, spins and orbital momentum. These conditions become ambiguous near the topological bifurcation points; in certain cases they also depend on the relative phases of the spin and orbital motions. In the case of merging magnetized neutron stars, if one of the stars is a millisecond pulsar, spinning at $\sim$ 10 msec, the global resonance $\omega_1+\omega_2= 2 \Omega$ (spin-plus beat is two times the orbital period) occurs approximately a second before the merger; the total energy of the flare can be as large as $10\%$ of the total magnetic energy, producing bursts of luminosity $\sim 10^{44}$ erg s$^{-1}$.

We present results of an optical spectroscopic survey using SALT and NOT to build a WISE mid-infrared color based, dust-unbiased sample of powerful radio-bright ($>$200 mJy at 1.4 GHz) AGN for the MeerKAT Absorption Line Survey (MALS). Our sample has 250 AGN (median $z=1.8$) showing emission lines, 26 with no emission lines, and 27 without optical counterparts. M1312-2026, the highest redshift object ($z=5.068$) in our sample, is the most radio-loud ($R$=$1.4\times 10^4$) AGN known at $z > 5$. Overall, our sample is fainter ($\Delta i$=0.6 mag) and redder ($\Delta(g-i)$=0.2 mag) than radio-selected quasars, and representative of fainter quasar population detected in optical surveys. About 20% of the sources are narrow line AGN (NLAGN) $-$ 65% of these, at $z < 0.5$ are galaxies without strong nuclear emission, and 10% at $z>1.9$, have emission line ratios similar to radio galaxies. The farthest NLAGN in our sample is M1513$-$2524 ($z_{em}=3.132$), and the largest (size$\sim$330 kpc) is M0909$-$3133 ($z_{em}=0.884$). We discuss in detail 110 AGN at $1.9 < z < 3.5$. Despite representing the radio loudest quasars (median $R$=3685), their Eddington ratios are similar to the SDSS quasars having lower $R$. We detect 4 CIV BALQSOs, all among AGN with least $R$, and highest black hole masses and Eddington ratios. The BAL detection rate ($4^{+3}_{-2}$%) is consistent with that seen in extremely powerful ($L_{1.4GHz}>10^{25}$ WHz$^{-1}$) quasars. Using optical light-curves, radio polarization and $\gamma$-ray detections, we identify 7 high-probability BL Lacs. We also summarize the full MALS footprint to search for HI 21-cm and OH 18-cm lines at $z<2$.

We examine nonrelativistic particles that decay into relativistic products in big rip, little rip, and pseudo-rip models for the future evolution of the universe. In contrast to decays that occur in standard $\Lambda$CDM, the evolution of the ratio $r$ of the energy density of the relativistic decay products to the energy density of the initially decaying particles can decrease with time in all of these models. In big rip and little rip models, $r$ always goes to zero asymptotically, while this ratio evolves to infinity or a constant in pseudo-rip models.

We conduct searches for continuous gravitational waves from seven pulsars, that have not been targeted in continuous wave searches of Advanced LIGO data before. We target emission at exactly twice the rotation frequency of the pulsars and in a small band around such frequency. The former search assumes that the gravitational wave quadrupole is changing phase-locked with the rotation of the pulsar. The search over a range of frequencies allows for differential rotation between the component emitting the radio signal and the component emitting the gravitational waves, for example the crust or magnetosphere versus the core. Timing solutions derived from the Arecibo 327-MHz Drift-Scan Pulsar Survey (AO327) observations are used. No evidence of a signal is found and upper limits are set on the gravitational wave amplitude. For one of the pulsars we probe gravitational wave intrinsic amplitudes just a factor of 3.8 higher than the spin-down limit, assuming a canonical moment of inertia of $10^{38}$ kg m$^2$. Our tightest ellipticity is $1.7 \times 10^{-8}$, which is a value well within the range of what a neutron star crust could support.

The merger remnant NGC 34 is a local luminous infrared galaxy (LIRG) hosting a nuclear starburst and a hard X-ray source associated with a putative, obscured Seyfert~2 nucleus. In this work, we use adaptive optics assisted near infrared (NIR) integral field unit observations of this galaxy to map the distribution and kinematics of the ionized and molecular gas in its inner $\mathrm{1.2\,kpc \times 1.2\,kpc}$, with a spatial resolution of 70~pc. The molecular and ionized gas kinematics is consistent with a disc with projected major axis along a mean PA~=~$\mathrm{-9^{\circ}.2 \pm 0^{\circ}.9}$. Our main findings are that NGC~34 hosts an AGN and that the nuclear starburst is distributed in a circumnuclear star-formation ring with inner and outer radii of $\approx$~60 and 180~pc, respectively, as revealed by maps of the $\mathrm{[Fe II] / Pa\beta}$ and $\mathrm{H_{2} / Br\gamma}$ emission-line ratios, and corroborated by PCA Tomography analysis. The spatially resolved NIR diagnostic diagram of NGC~34 also identifies a circumnuclear structure dominated by processes related to the stellar radiation field and a nuclear region where $[Fe II]$ and H$_2$ emissions are enhanced relative to the hydrogen recombination lines. We estimate that the nuclear X-ray source can account for the central H$_2$ enhancement and conclude that $[Fe II]$ and H$_2$ emissions are due to a combination of photo-ionization by young stars, excitation by X-rays produced by the AGN and shocks. These emission lines show nuclear, broad, blue-shifted components that can be interpreted as nuclear outflows driven by the AGN.

A test particle in a noncoplanar orbit about a member of a binary system can undergo Kozai-Lidov oscillations in which tilt and eccentricity are exchanged. An initially circular highly inclined particle orbit can reach high eccentricity. We consider the nonlinear secular evolution equations previously obtained in the quadrupole approximation. For the important case that the initial eccentricity of the particle orbit is zero, we derive an analytic solution for the particle orbital elements as a function of time that is exact within the quadrupole approximation. The solution involves only simple trigonometric and hyperbolic functions. It simplifies in the case that the initial particle orbit is close to being perpendicular to the binary orbital plane. The solution also provides an accurate description of particle orbits with nonzero but sufficiently small initial eccentricity. It is accurate over a range of initial eccentricity that broadens at higher initial inclinations. In the case of an initial inclination of pi/3, an error of 1% at maximum eccentricity occurs for initial eccentricities of about 0.1.

Large-scale surveys over the last years have revealed about 300 QSOs at redshift above 6. Follow-up observations identified surprising properties, such as the very high black hole (BH) masses, spatial correlations with surrounding cold gas of the host galaxy, or high CIV-MgII velocity shifts. In particular, the discovery of luminous high-redshift quasars suggests that at least some black holes likely have large masses at birth and grow efficiently. We aim at quantifying quasar pairs at high redshift for a large sample of objects. This provides a new key constraint on a combination of parameters related to the origin and assembly for the most massive black holes: BH formation efficiency and clustering, growth efficiency and relative contribution of BH mergers. We observed 116 spectroscopically confirmed QSOs around redshift 6 with the simultaneous 7-channel imager GROND in order to search for companions. Applying identical colour-colour cuts as for those which led to the spectroscopically confirmed QSO, we perform LePHARE fits to the 26 best QSO pair candidates, and obtained spectroscopic observations for 11 of those. e do not find any QSO pair with a companion brighter than M1450(AB) < -26 mag within our 0.1-3.3 h^-1 cMpc search radius, in contrast to the serendipitous findings in the redshift range 4--5. However, a low fraction of such pairs at this luminosity and redshift is consistent with indications from present-day cosmological-scale galaxy evolution models. In turn, the incidence of L- and T-type brown dwarfs which occupy a similar colour space as z ~ 6 QSOs, is higher than expected, by a factor of 5 and 20, respectively.

As material from an infalling protostellar envelope hits the forming disk, an accretion shock may develop which could (partially) alter the envelope material entering the disk. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) indicate that emission originating from warm SO and SO$_2$ might be good tracers of such accretion shocks. The goal of this work is to test under what shock conditions the abundances of gas-phase SO and SO$_2$ increase in an accretion shock at the disk-envelope interface. Detailed shock models including gas dynamics are computed using the Paris-Durham shock code for non-magnetized J-type accretion shocks in typical inner envelope conditions. The effect of pre-shock density, shock velocity, and strength of the ultraviolet (UV) radiation field on the abundance of warm SO and SO$_2$ is explored. Warm gas-phase chemistry is efficient in forming SO under most J-type shock conditions considered. In lower-velocity (~3 km/s) shocks, the abundance of SO is increased through subsequent reactions starting from thermally desorbed CH$_4$ toward H$_2$CO and finally SO. In higher velocity (>4 km/s) shocks, both SO and SO$_2$ are formed through reactions of OH and atomic S. The strength of the UV radiation field is crucial for SO and in particular SO$_2$ formation through the photodissociation of H$_2$O. Thermal desorption of SO and SO$_2$ ice is only relevant in high-velocity (>5 km/s) shocks at high densities ($>10^7$ cm$^{-3}$). Warm emission from SO and SO$_2$ is a possible tracer of accretion shocks at the disk-envelope interface as long as a local UV field is present. Additional observations with ALMA at high-angular resolution could provide further constraints. Moreover, the James Webb Space Telescope will give access to other possible slow, dense shock tracers such as H$_2$, H$_2$O, and [S I] 25$\mu$m.

We report on two Chandra observations of the quasar PSO J0309+27, the most distant blazar observed so far (z=6.1), performed eight months apart, in March and November 2020. Previous Swift-XRT observation showed that this object is one of the brightest X-ray sources beyond redshift 6.0 ever observed so far. This new data-set confirmed the high flux level and unveiled a spectral change occurred on a very short timescale (250s rest-frame), caused by a significant softening of the emission spectrum. This kind of spectral variability, on a such short interval, has never been reported in the X-ray emission of a flat spectrum radio quasar. A possible explanation is given by the emission produced by the inverse Compton scatter of the quasar UV photons by the cold electrons present in a fast shell moving along the jet. Although this bulk comptonization emission should be an unavoidable consequence of the standard leptonic jet model, this would be the first time that it is observed.

We compute the spectrum of tensor perturbations in warm inflation. We find that the spectrum, besides the standard component $\propto {H^2}/{M_P^2}$ associated to the amplification of the tensor vacuum fluctuations, acquires a component $\propto {\ell_{\rm mfp}\,T^5}/{M_P^4}$, where $\ell_{\rm mfp}$ and $T$ are respectively the mean free path and the temperature of the thermal degrees of freedom. The new contribution is due to the direct production of gravitational waves by the thermal bath, and can exceed the standard one in a viable region of parameter space. This contribution is dominated by thermal fluctuations at scales longer than $\ell_{\rm mfp}$.

Whether studying neutrinos for their own sake or as a messenger particle, neutrino cross sections are critically important for numerous analyses. On the low energy side, measurements from accelerator experiments reach up to a few 100s of GeV. On the high energy side, neutrino-earth absorption measurements extend down to a few TeV. The intermediate energy range has yet to be measured experimentally. This work is made possible by the linear relationship between the event rate and cross section, and will utilize IceCube muon neutrino data collected between 2010 and 2018. An advanced energy reconstruction, tailored to the unique properties of the energy range and using the full description of photon propagation in ice, is applied to an event sample of neutrino-induced through-going muons to perform a forward folding analysis.

Star-forming galaxies with stellar masses below $10^{10}\,\mathrm{M}_\odot$ make up the bulk of the galaxy population at $z>2$. The properties of the cold gas in these galaxies can only be probed in very deep ALMA observations or by targeting strongly lensed galaxies. Here we report the results of a pilot survey using the Atacama Compact Array (ACA) of molecular gas in the most strongly magnified galaxies selected as giant arcs in optical data. The selection in rest-frame UV wavelengths ensures that sources are regular star forming galaxies, without a priori indications of intense dusty starburst activity. We conducted Band 4 and Band 7 observations to detect mid-J CO, [C I] and thermal continuum as molecular gas tracers from three strongly lensed systems at $z\approx 2-3$: our targets are SGAS J1226651.3+215220, SGAS J003341.5+024217 and the Sunburst Arc. The measured molecular mass is then projected onto the source plane with detailed lens models developed from high resolution HST observations. Multiwavelength photometry is then used to obtain the intrinsic stellar mass and star formation rate via SED fitting. In only one of the four sources are the three tracers robustly detected, while in the other three sources they are either undetected or detected in continuum only. The implied molecular gass masses range from $4\times 10^9\,\mathrm{M}_\odot$ in the detected source to an upper limit of $\sim 10^8\,\mathrm{M}_\odot$ in the most magnified source. The inferred gas fraction and gas depletion timescale are found to lie approximately 0.5 to 1.0 dex below the established scaling relations based on previous studies of unlensed massive galaxies. Our results indicate that the cold gas content of intermediate to low mass galaxies should not be extrapolated from the trends seen in more massive high-$z$ galaxies. (Abridged abstract)

We present a model for the formation of the Magellanic Stream (MS) due to ram pressure stripping. We model the history of the Small and Large Magellanic Clouds in the recent cosmological past in a static Milky Way potential with diffuse halo gas, using observationally motivated orbits for the Magellanic Clouds derived from HST proper motions within the potential of the Milky Way. This model is able to reproduce the trailing arm but does not reproduce the leading arm feature, which is common for models of the stream formation that include ram pressure stripping effects. Our model produces a good match to observations (including the densities and line-of-sight velocities of the stream, as well as the positions and velocities of the satellites at present day) when we include a diffuse halo component for the Milky Way. From analyzing our grid of models, we find that there is a direct correlation between the observed stream length in our simulations and the mass of the Milky Way. For the observed MS length, the inferred Milky Way mass is $1.5 \pm 0.3 \times 10^{12}$ $M_\odot$, which agrees closely with other independent measures of the Milky Way mass. We also discuss the MS in the context of HI streams in galaxy clusters, and find that the MS lies on the low-mass end of a continuum from Hickson groups to the Virgo cluster. As a tracer of the dynamical mass in the outer halo, the MS is a particularly valuable probe of the Milky Way's potential.

The IceCube Neutrino Observatory is a cubic kilometer-sized detector designed to detect neutrinos of astrophysical origin. However, muons created by cosmic rays interacting in the atmosphere pose a significant background for these astrophysical neutrinos particularly in the southern equatorial sky. Identifying neutrino events that start in the detector allows us to reduce the atmospheric muon component while retaining a high rate of starting neutrino events. The method presented today also rejects atmospheric neutrinos if they are accompanied by muons from the same cosmic ray shower, lowering the 50$\%$ purity threshold for astrophysical-to-atmospheric neutrinos from 100 TeV to ~10 TeV at declinations less than -25{\deg}. We use 10$\%$ (burn sample) of 9.5 years IceCube data to demonstrate the status of this dataset. We outline a planned measurement of the diffuse neutrino flux inclusive of theoretical and detector systematic uncertainties. In addition, we discuss searches for neutrino point sources and diffuse galactic plane neutrino emission in the Southern sky and plans to release high astrophysical-purity real-time alerts to the multi-messenger community.

We performed a time-resolved spectral analysis of 53 bright gamma-ray bursts (GRBs) observed by \textit{Fermi}/GBM. Our sample consists of 908 individual spectra extracted from the finest time slices in each GRB. We fitted them with the synchrotron radiation model by considering the electron distributions in five different cases: mono-energetic, single power-law, Maxwellian, traditional fast cooling, and broken power-law. Our results were further qualified through Bayesian Information Criterion (BIC) by comparing with the fit by empirical models, namely the so-called Band function and cut-off power-law models. Our study showed that the synchrotron models, except for the fast-cooling case, can successfully fit most observed spectra, with the single power-law case being the most preferred. We also found that the electron distribution indices for the single power-law synchrotron fit in more than half of our spectra exhibits flux-tracking behavior, i.e., the index increases/decreases with the flux increasing/decreasing, implying that the distribution of the radiating electrons is increasingly narrower with time before the flux peaks and becomes more spreading afterward. Our results indicate that the synchrotron radiation is still feasible as a radiation mechanism of the GRB prompt emission phase.

Atmospheric gravity (buoyancy) waves (GWs) are of great importance for the energy and momentum budget of all planetary atmospheres. Propagating upward waves carry energy and momentum from the lower atmosphere to thermospheric altitudes and re-distribute them there. On Mars, GWs dominate the variability of the thermosphere and ionosphere. We provide a comprehensive climatology of Martian thermospheric GW activity at solar minimum (end of Solar Cycle 24) inferred from measurements by Neutral Gas and Ions Mass Spectrometer on board Mars Atmosphere and Volatile EvolutioN (NGIMS/MAVEN). The results are compared and interpreted using a one-dimensional spectral nonlinear GW model. Monthly mean GW activity varies strongly as a function of altitude (150-230 km) between 6-25%, reaching a maximum at $\sim$170 km. GW activity systematically exhibits a local time variability with nighttime values exceeding those during daytime, in accordance with previous studies. The analysis suggests that the day-night difference is primarily caused by a competition between dissipation due to molecular diffusion and wave growth due to decreasing background density. Thus, convective instability mechanism is likely to play a less important role in limiting GW amplitudes in the upper thermosphere, which explains their local time behavior.

We found that ZTF J185139.81+171430.3 = ZTF18abnbzvx shows a very short [0.00858995(3) d = 12.37 min] and large-amplitude (0.8 mag) coherent variations using Public Data Release of Zwicky Transient Facility observations. The only known object that shows similar very short period, large-amplitude and coherent variations is the unique white dwarf pulsar AR Sco. The variations in ZTF J185139.81+171430.3 may arise from the mechanism as in AR Sco and should deserve attention.

A key problem in protoplanetary disc evolution is understanding the efficiency of dust radial drift. This process makes the observed dust disc sizes shrink on relatively short timescales, implying that discs started much larger than what we see now. In this paper we use an independent constraint, the gas radius (as probed by CO rotational emission), to test disc evolution models. In particular, we consider the ratio between the dust and gas radius, $R_{\rm CO}/R_{\rm dust}$. We model the time evolution of protoplanetary discs under the influence of viscous evolution, grain growth, and radial drift. Then, using the radiative transfer code RADMC with approximate chemistry, we compute the dust and gas radii of the models and investigate how $R_{\rm CO}/R_{\rm dust}$ evolves. Our main finding is that, for a broad range of values of disc mass, initial radius, and viscosity, $R_{\rm CO}/R_{\rm dust}$ becomes large (>5) after only a short time (<1 Myr) due to radial drift. This is at odds with measurements in young star forming regions such as Lupus, which find much smaller values, implying that dust radial drift is too efficient in these models. Substructures, commonly invoked to stop radial drift in large, bright discs, must then be present, although currently unresolved, in most discs.

Blazars are Active Galactic Nuclei characterized by relativistic jets launched in the vicinity of the central engine (i.e. a supermassive black hole), that are oriented close to our line of sight. Their peculiar orientation makes them very efficient tracers of the overall jetted population, and due to their brightness they can be visible up to very high redshifts. A deep knowledge of these objects can provide fundamental clues to the models of formation and growth of the first supermassive black holes, but their search in the early Universe must be careful and follow a systematic approach. The discovery in the last $\sim15$ years of extremely massive blazars at very high redshifts ($M_{\rm BH}>10^9M_\odot, \, z>4$) revolutionized our perception of their earliest evolution: there seem to be different formation epochs for extremely massive black holes hosted in jetted ($z\sim4$) and non-jetted systems ($z\sim2.5$). This is not easy to explain, since one would expect that jetted sources accrete less efficiently. Small differences in the population are also derived from the search of such high-$z$ sources: we will go through the open questions in order to understand where the common knowledge stands and which steps must be undertaken to better understand the formation and common evolution of supermassive black holes and jets in the early Universe.

The IceCube Upgrade is an extension of the IceCube detector at the geographic South Pole. It consists of seven new strings with novel instrumentation. More than 430 multi-PMT optical modules called "mDOMs", housing 24 3-inch PMTs each, will be produced for the Upgrade. This will require testing and pre-calibration on a short timescale of more than 10,000 PMTs prior to assembly and deployment. We present the design of a PMT testing facility that enables simultaneous testing of roughly 100 PMTs per day at temperatures down to -20{\deg}C. The design is implemented at RWTH Aachen University and TU Dortmund University in parallel to achieve a throughput of up to 1,000 PMTs per week. This will enable a steady supply of tested PMTs to the production sites, which is critical for the Upgrade, as well as the future IceCube-Gen2 project.

We analyze radio bursts observed in events with interacting/non-interacting CMEs that produced major SEPs (Ip $>$ 10 MeV) fromApril 1997 to December 2014.We compare properties of meter (m), deca-hectometer (DH) type II as well as DH type III bursts, and time lags for interacting-CME-associated (IC) events and non-interacting-CME-associated (NIC) events. About 70\% of radio emissions were observed in events of both types from meters to kilometers. We found high correlations between the drift rates and mid-frequencies of type II radio bursts calculated as the mean geometric between their starting and ending frequencies for both NIC and IC-associated events (Correlation coefficient \textit{R}$^{2}$ = 0.98, power-law index $\varepsilon$ = 1.68 $\pm $ 0.16 and \textit{R}$^{2}$ = 0.93, $\varepsilon$ = 1.64 $\pm $ 0.19 respectively).We also found a correlation between the frequency drift rates of DH type II bursts and space speeds of CMEs in NIC-associated events. The absence of such correlation for IC-associated events confirms that the shock speeds changed in CME--CME interactions. For the events with western source locations, the mean peak intensity of SEPs in IC-associated events is four times larger than that in NIC-associated SEP events. From the mean time lags between the start times of SEP events and the start of m, DH type II, and DH type III radio bursts, we inferred that particle enhancements in NIC-associated SEP events occurred earlier than in IC-associated SEP events. The difference between NIC events and IC events in the mean values of parameters of type II and type III bursts is statistically insignificant.

While an axion-clouded black hole (BH) encounters a pulsar (PSR) or has a PSR companion, a "gravitational molecule" can be formed. In such a system, the axion cloud evolves at the binary hybrid orbitals, as it happens at microscopic level to electron cloud in a chemical molecule. To demonstrate this picture, we develop a semi-analytical formalism using the method of linear combination of atomic orbitals with an adiabatic approximation. An oscillating axion-cloud profile and a perturbed binary rotation, together with unique and novel detection signals, are then predicted. Remarkably, the proposed PSR timing and polarization observables, namely the oscillation of periastron time shift and the birefringence with multiple modulations, correlate in pattern, and thus can be properly combined to strengthen the detection.

At the IceCube Neutrino Observatory, a Surface Array Enhancement is planned, consisting of 32 hybrid stations, placed within the current IceTop footprint. This surface enhancement will considerably increase the detection sensitivity to cosmic rays in the 100 TeV to 1 EeV primary energy range, measure the effects of snow accumulation on the existing IceTop tanks and serve as R&D for the possible future large-scale surface array of IceCube-Gen2. Each station has one central hybrid DAQ, which reads out 8 scintillation detectors and 3 radio antennas. The radio antenna SKALA-2 is used in this array due to its low-noise, high amplification and sensitivity in the 70-350 MHz frequency band. Every scintillation detector has an active area of 1.5 m$^2$ organic plastic scintillators connected by wavelength-shifting fibers, which are connected to a silicon photomultiplier. The signals from the scintillation detectors are integrated and digitized by a local custom electronics board and transferred to the central DAQ. When triggered by the scintillation detectors, the filtered and amplified analog waveforms from the radio antennas are read out and digitized by the central DAQ. A full prototype station has been developed and built and was installed at the South Pole in January 2020. It is planned to install the full array by 2026. In this contribution the hardware design of the array as well as the installation plans will be presented.

Globular clusters (GCs) are important tools to rebuild the accretion history of a galaxy. There are newly discovered GCs in the Sagittarius (Sgr) dwarf galaxy, that can be used as probes of the accretion event onto the Milky Way (MW). Our main aim is to characterize the Sgr GC system by measuring its main physical parameters. We build the optical and near-IR color-magnitude diagrams (CMDs) for 21 new Sgr GCs using the VISTA Variables in the Via Lactea Extended Survey (VVVX) near-IR database combined with the Gaia EDR3 optical database. We derive metallicities and ages for all targets, using the isochrone-fitting method and the RGB-slope and metallicity relation. The total luminosities are calculated both in the near-IR and in the optical. We construct the metallicity distribution (MD), the globular cluster luminosity function (GCLF), and the age-metallicity relation for the Sgr GC system. We find 17 metal-rich GCs with -0.9<[Fe/H]<-0.3, 4 metal-poor GCs with -2.0<[Fe/H]<-1.1 in the new Sgr GC sample. Even though our age estimates are rough, we find that the metal-poor GCs are consistent with an old population with an average age of ~13 Gyr, while the metal-rich GCs show a wider age range, between 6-8 Gyr and 10-13 Gyr. We compare the MD and the GCLF for the Sgr GC system with those of the MW, M31 and Large Magellanic Cloud (LMC) galaxies. We conclude that the majority of the metal-rich GCs are located within the main body of Sgr galaxy. We confirm that the GCLF is not a universal distribution, since the Sgr GCLF peaks at fainter luminosities than the GCLFs of the MW, M31 and LMC. The MD shows a double-peaked distribution, and we note that the metal-rich population looks like the MW bulge GCs. We compared our results with the literature concluding that the Sgr progenitor could have been a reasonably large galaxy able to retain the SNe ejecta, thus enriching its ISM.

The IceCube Neutrino Observatory first observed a diffuse flux of high energy astrophysical neutrinos in 2013. Since then, this observation has been confirmed in multiple detection channels such as high energy starting events, cascades, and through-going muon tracks. Combining these event selections into a high statistics global fit of 10 years of IceCube's neutrino data could strongly improve the understanding of the diffuse astrophysical neutrino flux: challenging or confirming the simple unbroken power-law flux model as well as the astrophysical neutrino flux composition. One key component of such a combined analysis is the consistent modelling of systematic uncertainties of different event selections. This can be achieved using the novel SnowStorm Monte Carlo method which allows constraints to be placed on multiple systematic parameters from a single simulation set. We will report on the status of a new combined analysis of through-going muon tracks and cascades. It is based on a consistent all flavor neutrino signal and background simulation using, for the first time, the SnowStorm method to analyze IceCube's high-energy neutrino data. Estimated sensitivities for the energy spectrum of the diffuse astrophysical neutrino flux will be shown.

Context. Intermediate- to high-mass stars are the least numerous types of stars and they are less well understood than their more numerous low-mass counterparts in terms of their internal physical processes. Modelling the photometric variability of a large sample of main-sequence intermediate- to high-mass stars in eclipsing binary systems will help to improve the models for such stars. Aims. Our goal is to compose a homogeneously compiled sample of main-sequence intermediate- to high-mass OBA-type dwarfs in eclipsing binary systems from TESS photometry. We search for binaries with and without pulsations and determine their approximate ephemerides. Methods. Our selection starts from a catalogue of dwarfs with colours corresponding to those of OBA-type dwarfs in the TESS Input Catalog. We develop a new automated method aimed at detecting eclipsing binaries in the presence of strong pulsational and/or rotational signal relative to the eclipse depths and apply it to publicly available 30-min cadence TESS light curves. Results. Using targets with TESS magnitudes below 15 and cuts in the 2MASS magnitude bands of $J - H < 0.045$ and $J - K < 0.06$ as most stringent criteria, we arrive at a total of 189 981 intermediate- to high-mass candidates, 91193 of which have light curves from at least one of two data reduction pipelines. The eclipsing binary detection and subsequent manual check for false positives resulted in 3155 unique OBA-type eclipsing binary candidates. Conclusions. Our sample of eclipsing binary stars in the intermediate- to high-mass regime allows for future binary (and asteroseismic) modelling with the aim to better understand the internal physical processes in this hot part of the main sequence.

In this third paper of a series describing direction dependent corrections for polarimetric radio imaging, we present the the A-to-Z solver methodology to model the full Jones antenna aperture illumination pattern (AIP) with Zernike polynomials. In order to achieve thermal noise limited imaging with modern radio interferometers, it is necessary to correct for the instrumental effects of the antenna primary beam (PB) as a function of time, frequency, and polarization. The wideband AW projection algorithm enables those corrections provided an accurate model of the AIP is available. We present the A-to-Z solver as a more versatile algorithm for the modeling of the AIP. It employs the orthonormal circular Zernike polynomial basis to model the measured full Jones AIP. These full Jones models are then used to reconstruct the full Mueller AIP repsonse of an antenna, in principle accounting for all the off-axis leakage effects of the primary beam. The A-to-Z solver is general enough to accomodate any interferometer for which holographic measurements exist, we have successfully modelled the AIP of VLA, MeerKAT and ALMA as a demonstration of its versatility. We show that our models capture the PB morphology to high accuracy within the first 1-2 sidelobes, and show the viability of full Mueller gridding and deconvolution for any telescope given high quality holographic measurements.

Galaxy clusters provide an excellent probe in various research fields in astrophysics and cosmology. However, the number of galaxy clusters detected so far in the $AKARI$ North Ecliptic Pole (NEP) field is limited. In this work, we provide galaxy cluster candidates in the $AKARI$ NEP field with the minimum requisites based only on coordinates and photometric redshift (photo-$z$) of galaxies. We used galaxies detected in 5 optical bands ($g$, $r$, $i$, $z$, and $Y$) by the Subaru Hyper Suprime-Cam (HSC), assisted with $u$-band from Canada-France-Hawaii Telescope (CFHT) MegaPrime/MegaCam, and IRAC1 and IRAC2 bands from the $Spitzer$ space telescope for photo-$z$ estimation. We calculated the local density around every galaxy using the 10$^{th}$-nearest neighbourhood. Cluster candidates were determined by applying the friends-of-friends algorithm to over-densities. 88 cluster candidates containing 4390 member galaxies below redshift 1.1 in 5.4 deg$^2$ have been detected. The reliability of our method was examined through false detection tests, redshift uncertainty tests, and applications on the COSMOS data, giving false detection rates of 0.01 to 0.05 and recovery rate of 0.9 at high richness. 3 X-ray clusters previously observed by $ROSAT$ and $Chandra$ were recovered. The cluster galaxies show higher stellar mass and lower star formation rate (SFR) compared to the field galaxies in two-sample Z-tests. These cluster candidates are useful for environmental studies of galaxy evolution and future astronomical surveys in the NEP, where $AKARI$ has performed unique 9-band mid-infrared photometry for tens of thousands of galaxies.

We present $V$-band photometry of the 20,000 brightest asteroids using data from the All-Sky Automated Survey for Supernovae (ASAS-SN) between 2012 and 2018. We were able to apply the convex inversion method to more than 5,000 asteroids with more than 60 good measurements in order to derive their sidereal rotation periods, spin axis orientations, and shape models. We derive unique spin state and shape solutions for 760 asteroids, including 163 new determinations. This corresponds to a success rate of about 15%, which is significantly higher than the success rate previously achieved using photometry from surveys. We derive the first sidereal rotation periods for additional 69 asteroids. We find good agreement in spin periods and pole orientations for objects with prior solutions. We obtain a statistical sample of asteroid physical properties that is sufficient for the detection of several previously known trends, such as the underrepresentation of slow rotators in current databases, and the anisotropic distribution of spin orientations driven by the nongravitational forces. We also investigate the dependence of spin orientations on the rotation period. Since 2018, ASAS-SN has been observing the sky with higher cadence and deeper limiting magnitude, which will lead to many more new solutions in just a few years.

Black hole low-mass X-ray binaries (BH LMXBs) evolve in a similar way during outburst. Based on the X-ray spectrum and variability, this evolution can be divided into three canonical states: low/hard, intermediate and high/soft state. BH LMXBs evolve from the low/hard to the high/soft state through the intermediate state in some outbursts (here called "full outbursts"). However, in other cases, BH LMXBs undergo outbursts in which the source never reaches the high/soft state, here called "Failed-Transition outburst" (FT outbursts). From a sample of 56 BH LMXBs undergoing 128 outbursts, we find that $\sim$36% of these BH LMXBs experienced at least one FT outburst, and that FT outbursts represent $\sim$33% of the outbursts of the sample, showing that these are common events. We compare all the available X-ray data of full and FT outbursts of BH LMXBs from RXTE/PCA, Swift/BAT and MAXI and find that FT and full outbursts cannot be distinguished from their X-ray light curves, HIDs or X-ray variability during the initial 10-60 days after the outburst onset. This suggests that both types of outbursts are driven by the same physical process. We also compare the optical and infrared (O/IR) data of FT and full outbursts of GX 339-4. We found that this system is generally brighter in O/IR bands before an FT outburst, suggesting that the O/IR flux points to the physical process that later leads to a full or an FT outburst. We discuss our results in the context of models that describe the onset and evolution of outbursts in accreting X-ray binaries.

The tau lepton's supersymmetric partner, the stau, appears in some models as the next-to-lightest supersymmetric particle. Their decay process into the lightest superpartner is usually suppressed by supersymmetry breaking, which makes it a long-lived particle. In this scenario, its signature is a long, minimally ionizing track when traveling through the IceCube detector. Independent of their primary energy, the stau tracks appear like low-energy muons in the detector. A potential signal of staus would thus be an excess over muon tracks induced by atmospheric muon neutrinos. Our analysis focuses on the region around the horizon as here the ratio between stau signal and atmospheric background is largest. We will present the first sensitivity to constrain the stau mass using IceCube and demonstrate the potential of this analysis with future improvements.

Miras are fascinating stars. A kappa-mechanism in their atmosphere drives pulsations which produce changes in their photometric brightness, apparent spectral type and effective temperature. These pulsations also drive the formation of Balmer emission lines in the spectrum. This behaviour can be observed and investigated with small telescopes. We report on a three-year project combining spectroscopy and photometry to analyse the behaviour of Mira stars SU Cam and RY Cep, and describe how their brightness, colour, spectral type, effective temperature and Balmer emission vary over four pulsation cycles.

Magnetic reconnection associated with the tearing instability occurring in double-current sheet systems is investigated within the framework of reduced resistive magnetohydrodynamics (MHD) in a two-dimensional Cartesian geometry. The explosive non linear phase is particularly explored using the adaptive finite-element FINMHD code. The critical aspect ratio, that is defined as the minimum $L/x_s$ ratio (with $L$ and $x_s$ being the system length and half-distance between the two current layers respectively) necessary for non linear destabilization after the linear and early non linear saturation phases, is obtained. The latter threshold is independent of the details of the chosen initial equilibrium (double Harris-like magnetic profile) and of the resistivity. Its value is shown to be $4.7$, that is close and slightly smaller than the value of order $5$ deduced using a more particular equilibrium configuration in previous studies. The time dependence of the kinetic energy ($E_K$) is shown to follow a double exponential law, $E_K \propto \exp \ [e^{(\gamma^* t)} ]$, with a pseudo-growth rate $\gamma^* \simeq 0.1 \ t_A^ {-1}$ ($t_A$ being the characteristic Alfv\'en time) that is again independent of the configuration and resistivity. The mechanism offers a possible explanation for the sudden onset of explosive magnetic energy release occurring on the fast Alfv\'en time scale in disruptive events of astrophysical plasmas with pre-existing double current sheets like in the solar corona.

Clusters of galaxies -- with their turbulent magnetic fields and abundant matter content -- are a promising class of potential neutrino sources. Cosmic rays accelerated within the large-scale shocks,Active GalacticNuclei (AGN), or both can be confined in galaxy clusters over cosmological timescales and produce a steady flux of neutrinos in secondary interactions. The IceCube Neutrino Observatory has detected a diffuse flux of high-energy astrophysical neutrinos. After ten years of operations, however, the origin of this flux remains largely unconstrained. In this work, we perform a stacked search for neutrinos, using a population of over one thousand galaxy clusters detected by the Planck Satellite via the Sunyaev-Zeldovich (SZ) effect up to a redshift $z = 1$. We present the first results on the contribution of galaxy clusters to the diffuse neutrino flux and discuss the implications for various models of cosmic-ray acceleration in large-scale structures.

Since the discovery of the first fast radio burst (FRB) in 2007, and their confirmation as an abundant extragalactic population in 2013, the study of these sources has expanded at an incredible rate. In our 2019 review on the subject we presented a growing, but still mysterious, population of FRBs -- 60 unique sources, 2 repeating FRBs, and only 1 identified host galaxy. However, in only a few short years new observations and discoveries have given us a wealth of information about these sources. The total FRB population now stands at over 600 published sources, 24 repeaters, and 14 host galaxies. Higher time resolution data, sustained monitoring, and better localisation tools have given us insight into repeaters, host galaxies, burst morphology, source activity, progenitor models, and the use of FRBs as cosmological probes. The recent detection of a bright FRB-like burst from Galactic magnetar SGR 1935+2154 provides an important link between FRBs and magnetars. There also continue to be surprising discoveries, like periodic activity from repeaters and the localisation of one FRB source to a relatively nearby globular cluster associated with the M81 galaxy. In this review, we summarise the exciting observational results from the past few years and their impact on our understanding of the FRB population and proposed progenitor models. We build on the introduction to FRBs in our earlier review, update our readers on recent results, and discuss interesting avenues for exploration as the field enters a new regime where hundreds to thousands of new FRBs will be discovered and reported each year.

HiPERCAM is a portable, quintuple-beam optical imager that saw first light on the 10.4-m Gran Telescopio Canarias (GTC) in 2018. The instrument uses re-imaging optics and 4 dichroic beamsplitters to record $u_s g_s r_s i_s z_s$ ($320-1060$ nm) images simultaneously on its five CCD cameras, each of 3.1 arcmin (diagonal) field of view. The detectors in HiPERCAM are frame-transfer devices cooled thermo-electrically to 183 K, thereby allowing both long-exposure, deep imaging of faint targets, as well as high-speed (over 1000 windowed frames per second) imaging of rapidly varying targets. A comparison-star pick-off system in the telescope focal plane increases the effective field of view to 6.7 arcmin for differential photometry. Combining HiPERCAM with the world's largest optical telescope enables the detection of astronomical sources to $g_s \sim 23$ in 1 s and $g_s \sim 28$ in 1 h. In this paper we describe the scientific motivation behind HiPERCAM, present its design, report on its measured performance, and outline some planned enhancements.

We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study.

For a track based polarimeter, such as the Imaging X-ray Polarimetry Explorer (IXPE), the sensitivity to polarization depends on the modulation factor, which is a strong function of energy. In previous work, a likelihood method was developed that would account for this variation in order to estimate the minimum detectable polarization (MDP). That method essentially required that the position angles of individual events should be known precisely. In a separate work, however, it was shown that using a machine learning method for measuring event tracks can generate track angle uncertainties, which can be used in the analysis. Here, the maximum likelihood method is used as a basis for revising the estimate of the MDP in a general way that can include uncertainties in event track position angles. The resultant MDP depends solely upon the distribution of track angle uncertainties present in the input data. Due to the physics of the IXPE detectors, it is possible to derive a simple relationship between these angular uncertainties and the energy-dependent modulation function as a step in the process.

Numerous magnetohydrodynamic oscillations have been reported within solar pores over the past decades, including in line-of-sight (LOS) velocities, intensities, and magnetic field strengths. Our aim is to identify whether high-amplitude oscillations in the LOS magnetic field strength can be detected within a pore located in Active Region 12748 and to investigate which physical mechanisms could be responsible for them. A solar pore was observed on the 1st September 2019 using the GREGOR Infrared Spectrograph for around one hour. Full-Stokes vectors were sampled in a 37 A window containing the Fe I 15648.52 A line (effective Lande g-factor of 3). The LOS magnetic field strength is inferred using the strong-field approximation. The Stokes Inversion based on Response functions code is used to gain a more complete understanding the properties of the solar atmosphere at the locations of these oscillations. Oscillations of more than 100 G are observed in the LOS magnetic field in the period window 600-1272 s at three localised (>1"^2) regions. These oscillations have coherence across individual regions indicating that jitter cannot account for their occurrence. Longer-period amplitude variations, amplitudes over 200 G, are also detected but these have periods outside of the cone-of-influence. Numerical inversions confirm both oscillations in the LOS magnetic field strength at optical depths of around log-tau_5000=-0.5 (potentially caused by compression) and other effects (e.g., changes in the optical depth or the inclination of the field) may account for these changes. The oscillations in the separations of the Stokes-V lobes of the 15648.52 A line appear to be solar in nature. Future work will be required to understand whether these are truly oscillations in the magnetic field strength at a specific depth in the solar atmosphere or whether other effects are responsible for these signatures.

Redundant calibration is a technique in radio astronomy that allows calibration of radio arrays, whose antennas lie on a lattice by exploiting the fact that redundant baselines should see the same sky signal. Because the number of measured visibilities scales quadratically with the number of antennas, but the number of unknowns describing the individual antenna responses and the available information about the sky scale only linearly with the array size, the problem is always over-constrained as long as the array is big and dense enough. This is true even for non-lattice array configurations. In this work we study a generalized algorithm in which a per-antenna gain is replaced with a number of gains. We show that it can successfully describe data from an approximately redundant array on square lattice with pointing and geometry errors. We discuss the parameterization, limitations and possible extensions of this algorithm.

The Wavelength-shifting Optical Module (WOM) is a novel optical sensor that uses wavelength shifting and light guiding to substantially enhance the photosensitive area of UV optical modules. It has been designed for the IceCube Upgrade, a seven-string extension of the IceCube detector planned for the 2022/2023 South Pole deployment season. The WOM consists of a hollow quartz cylinder coated in wavelength shifting paint which serves as detection area and has two photomultipliers (PMTs) attached to the end faces. The light-collecting tube increases the effective photocathode area of the PMTs without producing additional dark current, making it suitable for low-signal, low-noise applications. We report on the design and performance of the WOM with a focus on the 12 modules in production for deployment in the IceCube Upgrade. While the WOM will be deployed in IceCube, its design is applicable to any large-volume particle detector based on the detection of Cherenkov light.

We present in this paper one of the largest galaxy morphological classification catalogues to date, including over 20 million of galaxies, using the Dark Energy Survey (DES) Year 3 data based on Convolutional Neural Networks (CNN). Monochromatic $i$-band DES images with linear, logarithmic, and gradient scales, matched with debiased visual classifications from the Galaxy Zoo 1 (GZ1) catalogue, are used to train our CNN models. With a training set including bright galaxies ($16\le{i}<18$) at low redshift ($z<0.25$), we furthermore investigate the limit of the accuracy of our predictions applied to galaxies at fainter magnitude and at higher redshifts. Our final catalogue covers magnitudes $16\le{i}<21$, and redshifts $z<1.0$, and provides predicted probabilities to two galaxy types -- Ellipticals and Spirals (disk galaxies). Our CNN classifications reveal an accuracy of over 99\% for bright galaxies when comparing with the GZ1 classifications ($i<18$). For fainter galaxies, the visual classification carried out by three of the co-authors shows that the CNN classifier correctly categorises disky galaxies with rounder and blurred features, which humans often incorrectly visually classify as Ellipticals. As a part of the validation, we carry out one of the largest examination of non-parametric methods, including $\sim$100,000 galaxies with the same coverage of magnitude and redshift as the training set from our catalogue. We find that the Gini coefficient is the best single parameter discriminator between Ellipticals and Spirals for this data set.

The paper presents an analysis of multi-wavelength data of two Lynds Bright Nebulae (LBN), LBN 140.07+01.64 and LBN 140.77$-$1.42. The 1420 MHz continuum map reveals an extended Y-shaped feature (linear extent ~3.7 deg), which consists of a linear part and a V-like structure. The sites LBN 140.07+01.64 and AFGL 437 are located toward the opposite sides of the V-like structure, and LBN 140.77$-$1.42 is spatially seen toward the linear part. Infrared-excess sources are traced toward the entire Y-feature, suggesting star formation activities. Infrared and sub-millimeter images show the presence of at least two large-scale dust filaments extended toward the LBN sources. The Herschel maps, which are available only toward the northern and central parts of the Y-feature, display the presence of higher column density (> 2.4 X 10^{21} cm^{-2}) of materials toward the filaments. Using the 12CO(1-0) line data, the distribution of molecular gas at [-42.7, -34.4] km/s traces the cloud associated with the Y-feature, and confirms the existence of filaments. The large-scale filaments appear to be possibly spatially twisted. There is a hint of an oscillatory-like velocity pattern along both the filaments, favouring their proposed twisted nature. It is the first study showing the possible twisting of filaments, which is more prominent in the northern and central parts of the Y-feature. This possible twisting/coupling of the large-scale filaments appears to be responsible for the observed star formation (including known OB-stars). The proposed physical process and the energetics of OB-stars together seem to explain the origin of the ionized Y-feature.

We forecast the number of galaxy clusters that can be detected via the thermal Sunyaev-Zeldovich (tSZ) signals by future cosmic microwave background (CMB) experiments, primarily the wide area survey of the CMB-S4 experiment but also CMB-S4's smaller delensing survey and the proposed CMB-HD experiment. We predict that CMB-S4 will detect 75,000 clusters with its wide survey of $f_{\rm sky}$ = 50% and 14,000 clusters with its deep survey of $f_{\rm sky}$ = 3%. Of these, approximately 1350 clusters will be at $z \ge 2$, a regime that is difficult to probe by optical or X-ray surveys. We assume CMB-HD will survey the same sky as the S4-Wide{}, and find that CMB-HD will detect $\times3$ more overall and an order of magnitude more $z \ge 2$ clusters than CMB-S4. These results include galactic and extragalactic foregrounds along with atmospheric and instrumental noise. Using CMB-cluster lensing to calibrate cluster tSZ-mass scaling relation, we combine cluster counts with primary CMB to obtain cosmological constraints for a two parameter extension of the standard model ($\Lambda CDM+\sum m_{\nu}+w_{0}$). Besides constraining $\sigma(w_{0})$ to $\lesssim 1\%$, we find that both surveys can enable a $\sim 2.5-4.5\sigma$ detection of $\sum m_{\nu}$, substantially strengthening CMB-only constraints. We also study the evolution of intracluster medium by modelling the cluster virialization ${\rm v}(z)$ and find tight constraints from CMB-S4, with further factors of 3-4 improvement for CMB-HD.

Using 3+1 classical lattice simulations, we follow the symmetry breaking pattern and subsequent non-linear evolution of a spectator field non-minimally coupled to gravity when the post-inflationary dynamics is given in terms of a stiff equation-of-state parameter. We find that the gradient energy density immediately after the transition represents a non-negligible fraction of the total energy budget, steadily growing to equal the kinetic counterpart. This behaviour is reflected on the evolution of the associated equation-of-state parameter, which approaches a universal value $1/3$, independently of the shape of non-linear interactions. Combined with kination, this observation allows for the generic onset of radiation domination for arbitrary self-interacting potentials, significantly extending previous results in the literature. The produced spectrum at that time is, however, non-thermal, precluding the naive extraction of thermodynamical quantities like temperature. Potential identifications of the spectator field with the Standard Model Higgs are also discussed.

The search for extraterrestrial intelligence (SETI) is a scientific endeavor which struggles with unique issues -- a strong indeterminacy in what data to look for and when to do so. This has led to attempts at finding both fundamental limits of the communication between extraterrestrial intelligence and human civilizations, as well as benchmarks so as to predict what kinds of signals we might most expect. Previous work has been formulated in terms of the information-theoretic task of communication, but we instead argue it should be viewed as a detection problem, specifically one-shot (asymmetric) hypothesis testing. With this new interpretation, we develop fundamental limits as well as provide simple examples of how to use this framework to analyze and benchmark different possible signals from extraterrestrial civilizations. We show that electromagnetic signaling for detection requires much less power than for communication, that detection as a function of power can be non-linear, and that much of the analysis in this framework may be addressed using computationally efficient optimization problems, thereby demonstrating tools for further inquiry.

We present a code to simulate the propagation of GWs in a potential well in the time domain. Our code uses the finite element method (FEM) based on the publicly available code {\it deal.ii}. We test our code using a point source monochromatic spherical wave. We examine not only the waveform observed by a local observer but also the global energy conservation of the waves. We find that our numerical results agree with the analytical predictions very well. Based on our code, we study the propagation of the leading wavefront of GWs in a potential well. We find that our numerical results agree with the results obtained from tracing null geodesics very well. Based on our simulations, we also test the accuracy of the thin-lens model in predicting the positions of the wavefront. We find that the analytical formula of the Shapiro-time delay is only accurate in regimes that are far away from the center of the potential well. However, near the optic axis, the analytical formula shows significant differences from the simulated ones. Besides these results, we find that unlike the conventional images in geometric optics, GWs can not be sheltered by the scatterer due to wave effects. The signals of GWs can circle around the scatterer and travel along the optic axis until arrive at a distant observer, which is an important observational consequence in such a system.

Random, multifield functions can set generic expectations for landscape-style cosmologies. We consider the inflationary implications of a landscape defined by a Gaussian random function, which is perhaps the simplest such scenario. Many key properties of this landscape, including the distribution of saddles as a function of height in the potential, depend only on its dimensionality, N, and a single parameter, {\gamma}, which is set by the power spectrum of the random function. We show that for saddles with a single downhill direction the negative mass term grows smaller, relative to the average mass, as N increases, a result with potential implications for the {\eta}-problem in landscape scenarios. For some power spectra Planck-scale saddles have {\eta} ~ 1 and eternal, topological inflation would be common in these scenarios. Lower-lying saddles typically have large {\eta}, but the fraction of these saddles which would support inflation is computable, allowing us to identify which scenarios can deliver a universe that resembles ours. Finally, by drawing inferences about the relative viability of different multiverse proposals we also illustrate ways in which quantitative analyses of multiverse scenarios are feasible.

We study gamma-ray line signatures from electroweakly interacting non-abelian spin-1 dark matter (DM). In this model, $Z_2$-odd spin-1 particles including a DM candidate have the SU(2)$_L$ triplet-like features, and the Sommerfeld enhancement is relevant in the annihilation processes. We derive the annihilation cross sections contributing to the photon emission and compare with the SU(2)$_L$ triplet fermions, such as Wino DM in the supersymmetric Standard Model. The Sommerfeld enhancement factor is approximately the same in both systems, while our spin-1 DM predicts the larger annihilation cross sections into $\gamma \gamma/ Z \gamma$ modes than those of the Wino by $\frac{38}{9}$. This is because a spin-1 DM pair forms not only $J=0$ but also $J=2$ partial wave states where $J$ denotes the total spin angular momentum. Our spin-1 DM also has a new annihilation mode into $Z_2$-even extra heavy vector and photon, $Z' \gamma$. For this mode, the photon energy depends on the masses of DM and the heavy vector, and thus we have a chance to probe the mass spectrum. The latest gamma-ray line search in the Galactic Center region gives a strong constraint on our spin-$1$ DM. We can probe the DM mass for $\lesssim 25.3~$TeV by the Cherenkov Telescope Array experiment even if we assume a conservative DM density profile.