Papers for Thursday, Apr 27 2023

We present multi-epoch spectropolarimetry and spectra for a sample of 14 Type IIn supernovae (SNe IIn). We find that after correcting for likely interstellar polarization, SNe IIn commonly show intrinsic continuum polarization of 1--3% at the time of peak optical luminosity, although a few show weaker or negligible polarization. While some SNe IIn have even stronger polarization at early times, their polarization tends to drop smoothly over several hundred days after peak. We find a tendency for the intrinsic polarization to be stronger at bluer wavelengths, especially at early times. While polarization from an electron scattering region is expected to be grey, scattering of SN light by dusty circumstellar material (CSM) may induce such a wavelength-dependent polarization. For most SNe IIn, changes in polarization degree and wavelength dependence are not accompanied by changes in the position angle, requiring that asymmetric pre-SN mass loss had a persistent geometry. While 2--3% polarization is typical, about 30% of SNe IIn have very low or undetected polarization. Under the simplifying assumption that all SN IIn progenitors have axisymmetric CSM (i.e. disk/torus/bipolar), then the distribution of polarization values we observe is consistent with similarly asymmetric CSM seen from a distribution of random viewing angles. This asymmetry has very important implications for understanding the origin of pre-SN mass loss in SNe IIn, suggesting that it was shaped by binary interaction.

Gravitational redshift is a classical effect of Einstein's General Relativity, already measured in stars, quasars and clusters of galaxies. We here identify the signature of gravitational redshift in the emission lines of active galaxies due to supermassive black holes and discuss their impact on cosmological inference from type Ia supernovae. Firstly, from the full width at half maximum of $H_{\beta}$ lines of 75 Seyfert type I galaxies of the AGN Black Hole Mass Database, we derive a gravitational redshift $z_g = (2.4 \pm 0.9) \times 10^{-4}$. Expanding this analysis to 86755 quasars from DR14 of SDSS we have a mean value $z_g \approx 2.7 \times 10^{-4}$. Then, by comparing the redshifts of 34 lines measured at the central and outer regions of LINER galaxies in the SAMI survey we obtain $z_g = (0.68 \pm 0.09) \times 10^{-4}$. These numbers are compatible with central black holes of $\approx 10^9$ solar masses and broad line regions of $\approx 1$~pc. For non-AGN galaxies the gravitational redshift is compatible with zero and, as they constitute most of SNe Ia host galaxies, the impact on the cosmological parameters is negligible.

The stochastic gravitational-wave background is imprinted on the times of arrival of radio pulses from millisecond pulsars. Traditional pulsar timing analyses fit a timing model to each pulsar and search the residuals of the fit for a stationary time correlation. This method breaks down at gravitational-wave frequencies below the inverse observation time of the array; therefore, existing analyses restrict their searches to frequencies above 1 nHz. An effective method to overcome this challenge is to study the correlation of secular drifts of parameters in the pulsar timing model itself. In this paper, we show that timing model correlations are sensitive to sub-nanohertz stochastic gravitational waves and perform a search using existing measurements of binary spin-down rates and pulsar spin-decelerations. We do not observe a signal at our present sensitivity, constraining the stochastic gravitational-wave relic energy density to $\Omega_\text{GW} ( f ) < 3.9 \times 10 ^{ - 9} $ at 450 pHz with sensitivity which scales as the frequency squared until approximately 10 pHz. We place additional limits on the amplitude of a power-law spectrum of $A_\star \lesssim 8\times10^{-15}$ for the spectral index expected from supermassive black hole binaries, $\gamma = 13/3$. If a detection of a supermassive black hole binary signal above 1 nHz is confirmed, this search method will serve as a critical complementary probe of the dynamics of galaxy evolution.

Fermi-LAT observations provide continuous and regularly-sampled measurements of gamma-ray photon flux for hundreds of blazars. Many of these light curves, spanning almost 15 years, have been thoroughly examined for periodicity in multiple studies. However, the possibility that blazars may exhibit irregularly repeating flaring patterns in their gamma-ray light curves has not been systematically explored. In this study, we aim to find repeating episodes of flaring activity in the 100 brightest blazars using Fermi-LAT light curves with various integration times. We use a Bayesian Blocks representation to convert the time series into strings of symbols and search for repeating sub-strings using a fuzzy search algorithm. As a result, we identify 27 repeated episodes in the gamma-ray light curves of 10 blazars. We find that the patterns are most likely produced in structured jets composed of a fast spine and a slower sheath. When individual emission features propagate in the spine, they scatter seed photons produced in the non-uniform sheath through the inverse Compton mechanism, resulting in a set of gamma-ray flares with a similar profile every such passage. Additionally, we explore the theoretically-predicted possibility that the spine-sheath structure facilitates the production of high-energy neutrinos in blazar jets. Using the catalogue of track-like events detected by the IceCube neutrino telescope, we find evidence supporting this hypothesis at a $3.5\sigma$ significance level.

We present a new spectroscopic study of the dwarf galaxy Hercules (d ~ 132 kpc) with data from the Anglo-Australian Telescope and its AAOmega spectrograph together with the Two Degree Field multi-object system to solve the conundrum that whether Hercules is tidally disrupting. We combine broadband photometry, proper motions from Gaia, and our Pristine narrow-band and metallicity-sensitive photometry to efficiently weed out the Milky Way contamination. Such cleaning is particularly critical in this kinematic regime, as both the transverse and heliocentric velocities of Milky Way populations overlap with Hercules. Thanks to this method, three new member stars are identified, including one at almost 10rh of the satellite. All three have velocities and metallicities consistent with that of the main body. Combining this new dataset with the entire literature cleaned out from contamination shows that Hercules does not exhibit a velocity gradient (d<v>/dX = 0.1+0.4/-0.2 km s-1 arcmin-1) and, as such, does not show evidence to undergo tidal disruption.

Stars far beyond the half-light radius of a galaxy suggest the existence of a mechanism able to move stars out of the region where most star formation has taken place. The formation of these "stellar halos" are usually ascribed to the effects of early mergers or Galactic tides, although fluctuations in the gravitational potential due to stellar feedback is also a possible candidate mechanism. A Bayesian algorithm is used to find new candidate members in the extreme outskirts of the Sculptor dwarf galaxy. Precise metallicities and radial velocities for two distant stars are measured from their spectra taken with the Gemini South GMOS spectrograph. The radial velocity, proper motion and metallicity of these targets are consistent with Sculptor membership. As a result, the known boundary of the Sculptor dwarf extends now out to an elliptical distance of $\sim10$ half-light radii, which corresponds to a projected physical distance of $\sim3$ kpc. As reported in earlier work, the overall distribution of radial velocities and metallicities indicate the presence of a more spatially and kinematically dispersed metal-poor population that surrounds the more concentrated and colder metal-rich stars. Sculptor's density profile shows a "kink" in its logarithmic slope at a projected distance of $\sim25$ arcmin (620 pc), which we interpret as evidence that Galactic tides have helped to populate the distant outskirts of the dwarf. We discuss further ways to test and validate this tidal interpretation for the origin of these distant stars.

Recent body of work points out that the mean-field approximation, widely employed to mimic the neutrino field within a neutrino-dense source, might give different results in terms of flavor evolution with respect to the correspondent many-body treatment. In this paper, we investigate whether such conclusions derived within a constrained framework should hold in an astrophysical context. We show that the plane waves, commonly adopted in the many-body literature to model the neutrino field, provide results that are crucially different with respect to the ones obtained using wavepackets of finite size streaming with a non-zero velocity. The many-body approach intrinsically includes coherent and incoherent scatterings. The mean-field approximation, on the other hand, only takes into account the coherent scattering in the absence of the collision term. Even if incoherent scatterings are included in the mean-field approach, the nature of the collision term is different from that in the many-body approach. Because of this, if only a finite number of neutrinos is considered, as often assumed, the two approaches naturally lead to different flavor outcomes. These differences are further exacerbated by vacuum mixing. We conclude that existing many-body literature, based on closed neutrino systems with a finite number of particles, is neither able to rule out nor assess the validity of the mean-field approach adopted to simulate the evolution of the neutrino field in dense astrophysical sources, which are open systems.

We present the first measurement of the HI mass function (HIMF) using data from MeerKAT, based on 276 direct detections from the MIGHTEE Survey Early Science data covering a period of approximately a billion years ($0 \leq z \leq 0.084 $). This is the first HIMF measured using interferometric data over non-group or cluster field, i.e. a deep blank field. We constrain the parameters of the Schechter function which describes the HIMF with two different methods: $1/\rm V_{\rm max}$ and Modified Maximum Likelihood (MML). We find a low-mass slope $\alpha=-1.29^{+0.37}_{-0.26}$, `knee' mass $\log_{10}(M_{*}/{\rm M_{\odot}}) = 10.07^{+0.24}_{-0.24}$ and normalisation $\log_{10}(\phi_{*}/\rm Mpc^{-3})=-2.34^{+0.32}_{-0.36}$ ($H_0 = 67.4$ kms$^{-1}$ Mpc$^{-1}$) for $1/\rm V_{\rm max}$ and $\alpha=-1.44^{+0.13}_{-0.10}$, `knee' mass $\log_{10}(M_{*}/{\rm M_{\odot}}) = 10.22^{+0.10}_{-0.13}$ and normalisation $\log_{10}(\phi_{*}/\rm Mpc^{-3})=-2.52^{+0.19}_{-0.14}$ for MML. When using $1/\rm V_{\rm max}$ we find both the low-mass slope and `knee' mass to be consistent within $1\sigma$ with previous studies based on single-dish surveys. The cosmological mass density of HI is found to be slightly larger than previously reported: $\Omega_{\rm HI}=5.46^{+0.94}_{-0.99} \times 10^{-4}h^{-1}_{67.4}$ from $1/\rm V_{\rm max}$ and $\Omega_{\rm HI}=6.31^{+0.31}_{-0.31} \times 10^{-4}h^{-1}_{67.4}$ from MML but consistent within the uncertainties. We find no evidence for evolution of the HIMF over the last billion years.

If long-range attractive forces exist and are stronger than gravity then cosmic halo formation can begin in the radiation-dominated era. We study a simple realization of this effect in a system where dark matter fermions have Yukawa interactions mediated by scalar particles, analogous to the Higgs boson in the standard model. We develop a self-consistent description of the system including exact background dynamics of the scalar field, and precise modelling of the fermion density fluctuations. For the latter, we provide accurate approximations for the linear growth as well as quantitative modelling of the nonlinear evolution using N-body simulations. We find that halo formation occurs exponentially fast and on scales substantially larger than simple estimates predict. The final fate of these halos remains uncertain, but could be annihilation, dark stars, primordial black holes, or even the existence of galaxy-sized halos at matter-radiation equality. More generally, our results demonstrate the importance of mapping scalar-mediated interactions onto structure formation outcomes and constraints for beyond the standard model theories.

Boson stars are gravitationally bound objects that arise in ultralight dark matter models and form in the centers of galactic halos or axion miniclusters. We systematically study the excitations of a boson star, taking into account the mixing between positive and negative frequencies introduced by gravity. We show that the spectrum contains zero-energy modes in the monopole and dipole sectors resulting from spontaneous symmetry breaking by the boson star background. We analyze the general properties of the eigenmodes and derive their orthogonality and completeness conditions which have non-standard form due to the positive-negative frequency mixing. The eigenvalue problem is solved numerically for the first few energy levels in different multipole sectors and the results are compared to the solutions of the Schr\"odinger equation in fixed boson star gravitational potential. The two solutions differ significantly for the lowest modes, but get close for higher levels. We further confirm the normal mode spectrum in 3D wave simulations where we inject perturbations with different multipoles. As an application of the normal mode solutions, we compute the matrix element entering the evaporation rate of a boson star immersed in a hot axion gas. The computation combines the use of exact wavefunctions for the low-lying bound states and of the Schr\"odinger approximation for the high-energy excitations.

With the precision now afforded by modern space-based photometric observations from the retired K2 and current TESS missions, the effects of general relativity (GR) may be detectable in the light curves of pulsating white dwarfs (WDs). Almost all WD models are calculated using a Newtonian description of gravity and hydrodynamics. To determine if inclusion of GR leads to observable effects, we used idealized models of compact stars and made side-by-side comparison of mode periods computed using a (i) Newtonian and (ii) GR description of the equilibrium structure and nonradial pulsations. For application to white dwarfs, it is only necessary to include the first post-Newtonian (1PN) approximation to GR. The mathematical nature of the linear nonradial pulsation problem is then qualitatively unchanged and the GR corrections can be written as extensions of the classic Dziembowski equations. As such, GR effects might easily be included in existing asteroseismology codes. The idealized stellar models are (i) \pn1 relativistic polytropes and (ii) stars with cold degenerate-electron equation of state featuring a near-surface chemical transition from $\mu_e = 2$ to $\mu_e = 1$, simulating a surface hydrogen layer. Comparison of Newtonian and 1PN normal mode periods reveals fractional differences on the order of the surface gravitational redshift $z$. For a typical WD, this fractional difference is $\sim 10^{-4}$ and is greater than the period uncertainty $\sigma_{\Pi}/\Pi$ of many white dwarf pulsation modes observed by TESS. A consistent theoretical modeling of periods observed in these stars should in principle include GR effects to 1PN order.

Lyman $\alpha$ (Ly$\alpha$) is now routinely used as a tool for studying high-redshift galaxies and its resonant nature means it can trace neutral hydrogen around star-forming galaxies. Integral field spectrograph measurements of high-redshift Ly$\alpha$ emitters indicate that significant extended Ly$\alpha$ halo emission is ubiquitous around such objects. We present a sample of redshift 0.23 to 0.31 galaxies observed with the Hubble Space Telescope selected to match the star formation properties of high-$z$ samples while optimizing the observations for detection of low surface brightness Ly$\alpha$ emission. The Ly$\alpha$ escape fractions range between 0.7\% and 37\%, and we detect extended Ly$\alpha$ emission around six out of seven targets. We find Ly$\alpha$ halo to UV scale length ratios around 6:1 which is marginally lower than high-redshift observations, and halo flux fractions between 60\% and 85\% -- consistent with high-redshift observations -- when using comparable methods. However, our targets show additional extended stellar UV emission: we parametrize this with a new double exponential model. We find that this parametrization does not strongly affect the observed Ly$\alpha$ halo fractions. We find that deeper H$\alpha$ data would be required to firmly determine the origin of Ly$\alpha$ halo emission, however, there are indications that H$\alpha$ is more extended than the central FUV profile, potentially indicating conditions favorable for the escape of ionizing radiation. We discuss our results in the context of high-redshift galaxies, cosmological simulations, evolutionary studies of the circumgalactic medium in emission, and the emission of ionizing radiation.

Extreme mass ratio inspirals (EMRIs) occur when stellar-mass compact objects begin a gravitational wave (GW) driven inspiral into massive black holes. EMRI waveforms can precisely map the surrounding spacetime, making them a key target for future space-based GW interferometers such as {\it LISA}, but their event rates and parameters are massively uncertain. One of the largest uncertainties is the ratio of true EMRIs (which spend at least thousands of orbits in the {\it LISA} band) and direct plunges, which are in-band for at most a handful of orbits and are not detectable in practice. In this paper, we show that the traditional dichotomy between EMRIs and plunges -- EMRIs originate from small semimajor axes, plunges from large -- does not hold for intermediate-mass black holes with masses $M_\bullet \lesssim 10^5 M_\odot$. In this low-mass regime, a plunge always has an $\mathcal{O}(1)$ probability of failing and transitioning into a novel ``cliffhanger'' EMRI. Cliffhanger EMRIs are more easily produced for larger stellar-mass compact objects, and are less likely for smaller ones. This new EMRI production channel can dominate volumetric EMRI rates $\dot{n}_{\rm EMRI}$ if intermediate-mass black holes are common in dwarf galactic nuclei, potentially increasing $\dot{n}_{\rm EMRI}$ by an order of magnitude.

Angular power spectra are central to the study of our Universe. In this paper, I develop a new method for the numeric evaluation and analytic estimation of the angular cross-power spectrum of two random fields using complex analysis and Picard- Lefschetz theory. The proposed continuous deformation of the integration domain resums the highly oscillatory integral into a convex integral whose integrand decays exponentially. This deformed integral can be quickly evaluated with conventional integration techniques. These methods can be used to quickly evaluate and estimate the angular power spectrum from the three-dimensional power spectrum for all angles (or multipole moments). This method is especially useful for narrow redshift bins, or samples with small redshift overlap, for which the Limber approximation has a large error.

We report LOFAR and VLA detections of extended radio emission in the spiral galaxy HCG 97b, hosted by an X-ray bright galaxy group. The extended radio emission detected at 144 MHz, 1.4 GHz and 4.8 GHz is elongated along the optical disk and has a tail that extends 27 kpc in projection towards the centre of the group at GHz frequencies or 60 kpc at 144 MHz. Chandra X-ray data show two off-nuclear ultra-luminous X-ray sources (ULXs) with the more distant one being a suitable candidate for an accreting intermediate-mass black hole (IMBH) embedded in an environment with an increased density of molecular gas. Given the observed morphology in optical, CO, and radio continuum, we propose that the galaxy is undergoing ram-pressure stripping and the relativistic plasma accelerated in star-forming regions is transported from the galactic disc by galaxy-intragroup medium interaction. Finally, we also demonstrate that the formation of the radio tail could, in principle, be the result of putative IMBH-induced activity, which could facilitate the stripping or inject the radio plasma via jets.

We present the analysis performed using the pyPipe3D pipeline for the 895 galaxies that comprises the eCALIFA data release Sanchez et al. submitted, data with a significantly improved spatial resolution (1.0-1.5"/FWHM). We include a description of (i) the analysis performed by the pipeline, (ii) the adopted datamodel for the derived spatially resolved properties and (iii) the catalog of integrated, characteristics and slope of the radial gradients for a set of observational and physical parameters derived for each galaxy. We illustrate the results of the analysis (i) using the NGC\,2906 as an archetypal galaxy, showing the spatial distribution of the different derived parameters and exploring in detail the properties of the ionized gas, and (ii) showing distribution of the spatial resolved ionized gas across the classical [OIII]/H$\beta$ vs. [NII]/H$\alpha$ for the whole galaxy sample. In general our results agree with previous published ones, however, tracing radial patterns and segregating individual ionized structures is improved when using the current dataset. All the individual galaxy dataproducts and the catalog discussed along this article are distributed as part of the eCALIFA data release this http URL

The wideband timing technique enables the high-precision simultaneous estimation of Times of Arrival (ToAs) and Dispersion Measures (DMs) while effectively modeling frequency-dependent profile evolution. We present two novel independent methods that extend the standard wideband technique to handle simultaneous multi-band pulsar data incorporating profile evolution over a larger frequency span to estimate DMs and ToAs with enhanced precision. We implement the wideband likelihood using the libstempo python interface to perform wideband timing in the tempo2 framework. We present the application of these techniques to the dataset of fourteen millisecond pulsars observed simultaneously in Band 3 (300 - 500 MHz) and Band 5 (1260 - 1460 MHz) of the upgraded Giant Metrewave Radio Telescope (uGMRT) as a part of the Indian Pulsar Timing Array (InPTA) campaign. We achieve increased ToA and DM precision and sub-microsecond root mean square post-fit timing residuals by combining simultaneous multi-band pulsar observations done in non-contiguous bands for the first time using our novel techniques.

In recent years, automatic classifiers of image cutouts (also called "stamps") have shown to be key for fast supernova discovery. The upcoming Vera C. Rubin Observatory will distribute about ten million alerts with their respective stamps each night, which it is expected to enable the discovery of approximately one million supernovae each year. A growing source of confusion for these classifiers is the presence of satellite glints, sequences of point-like-sources produced by rotating satellites or debris. The currently planned Rubin stamps will have a size smaller than the typical separation between these point sources. Thus, a larger field of view image stamp could enable the automatic identification of these sources. However, the distribution of larger field of view stamps would be limited by network bandwidth restrictions. We evaluate the impact of using image stamps of different angular sizes and resolutions for the fast classification of events (AGNs, asteroids, bogus, satellites, SNe, and variable stars), using available data from the Zwicky Transient Facility survey. We compare four scenarios: three with the same number of pixels (small field of view with high resolution, large field of view with low resolution, and a proposed multi-scale strategy) and a scenario with the full ZTF stamp that has a larger field of view and higher resolution. Our multi-scale proposal outperforms all the scenarios, with a macro f1-score of 87.39. We encourage Rubin and its Science Collaborations to consider the benefits of implementing multi-scale stamps as a possible update to the alert specification.

The lateral distribution of particles in extensive air showers from cosmic rays with energy above $10^{17}$ eV registered at the Yakutsk complex array was analyzed. Experimentally measured particle densities were compared to the predictions obtained within frameworks of three ultra-high energy hadron interaction models. The cosmic ray mass composition estimated by the readings of surface-based and underground detectors of the array is consistent with results based on the Cherenkov light lateral distribution data. A comparison was made with the results of direct measurement of the muon component performed at the Pierre Auger Observatory. It is demonstrated that the densities of muon flux measured at Yakutsk array are consistent with results of fluorescent light measurements and disagree with results on muons obtained at the Auger array.

In neutral hydrogen (HI) galaxy survey, a significant challenge is to identify and extract the HI galaxy signal from observational data contaminated by radio frequency interference (RFI). For a drift-scan survey, or more generally a survey of a spatially continuous region, in the time-ordered spectral data, the HI galaxies and RFI all appear as regions which extend an area in the time-frequency waterfall plot, so the extraction of the HI galaxies and RFI from such data can be regarded as an image segmentation problem, and machine learning methods can be applied to solve such problems. In this study, we develop a method to effectively detect and extract signals of HI galaxies based on a Mask R-CNN network combined with the PointRend method. By simulating FAST-observed galaxy signals and potential RFI impacts, we created a realistic data set for the training and testing of our neural network. We compared five different architectures and selected the best-performing one. This architecture successfully performs instance segmentation of HI galaxy signals in the RFI-contaminated time-ordered data (TOD), achieving a precision of 98.64% and a recall of 93.59%.

We present an extensive and well-characterized Chandra X-ray Galaxy Catalog (CGC) of 8557 galaxy candidates in the redshift range z ~ 0.04 - 0.7, optical luminosity 1010 - 1011 Lro, and X-ray luminosity (0.5-7 keV) LX = 2x1040 - 2x1043 erg s-1. We estimate ~5% false match fraction and contamination by QSOs. The CGC was extracted from the Chandra Source Catalog version 2 (CSC2) by cross-correlating with optical and IR all-sky survey data, including SDSS, PanSTARRS, DESI Legacy, and WISE. Our selection makes use of two main criteria that we have tested on the subsample with optical spectroscopical identification. (1) A joint selection based on X-ray luminosity (LX) and X-ray to optical flux ratio (FXO), which recovers 63% of the spectroscopically classified galaxies with a small contamination fraction (7%), a significant improvement over methods using LX or FXO alone (< 50% recovery). (2) A joint W1-W2 (W12) WISE color and LX selection that proves effective in excluding QSOs and improves our selection by recovering 72% of the spectroscopically classified galaxies and reducing the contamination fraction (4%). Of the CGC, 24% was selected by means of optical spectroscopy; 30% on the basis of LX, FXO, and W12; and 46% by using either the LX-FXO or the LX-W12 selection criteria. We have individually examined the data for galaxies with z < 0.1, which may include more than one CSC2 X-ray source, leading to the exclusion of 110 local galaxies. Our catalog also includes near-IR and UV data and galaxy morphological types.

X-ray bright optically normal galaxies (XBONGs) are galaxies with X-ray luminosities consistent with those of active galactic nuclei (AGNs) but no evidence of AGN optical emission lines. Crossmatching the Chandra Source Catalog version 2 (CSC2) with the Sloan Digital Sky Survey (SDSS) sample of spectroscopically classified galaxies, we have identified 817 XBONG candidates with LX > 1042 erg s-1 and X-ray to optical flux ratio FXO > 0.1. Comparisons with WISE colors and NIR, optical, UV, and radio luminosities show that the loci of XBONGs are in-between those of control samples of normal galaxies and quasars and are consistent with low-luminosity quasars. We find that 43% of the XBONG sample have X-ray colors suggesting NH > 1022 cm-2, double the fraction in the QSO sample, suggesting that a large fraction of XBONG are highly obscured AGNs. However, ~50% of the XBONGs are not obscured and have X-ray colors harder than those of normal galaxies. Some of these XBONGs have spatially extended X-ray emission. These characteristics suggest that they may be unidentified galaxy groups and clusters. Using the X-ray luminosity functions of QSOs and galaxies/groups/clusters, we estimate the approximate fraction of extended XBONGs to be < 20%. We also assess the approximate fraction of XBONGs whose AGN signatures are diluted by stellar light of host galaxies to be ~30%, based on their redshift and deviation from the extrapolation of the QSO LX-Lr relation.

The quest to find extraterrestrial life is a critical scientific endeavor with civilization-level implications. Icy moons in our solar system are promising targets for exploration because their liquid oceans make them potential habitats for microscopic life. However, the lack of a precise definition of life poses a fundamental challenge to formulating detection strategies. To increase the chances of unambiguous detection, a suite of complementary instruments must sample multiple independent biosignatures (e.g., composition, motility/behavior, and visible structure). Such an instrument suite could generate 10,000x more raw data than is possible to transmit from distant ocean worlds like Enceladus or Europa. To address this bandwidth limitation, Onboard Science Instrument Autonomy (OSIA) is an emerging discipline of flight systems capable of evaluating, summarizing, and prioritizing observational instrument data to maximize science return. We describe two OSIA implementations developed as part of the Ocean Worlds Life Surveyor (OWLS) prototype instrument suite at the Jet Propulsion Laboratory. The first identifies life-like motion in digital holographic microscopy videos, and the second identifies cellular structure and composition via innate and dye-induced fluorescence. Flight-like requirements and computational constraints were used to lower barriers to infusion, similar to those available on the Mars helicopter, "Ingenuity." We evaluated the OSIA's performance using simulated and laboratory data and conducted a live field test at the hypersaline Mono Lake planetary analog site. Our study demonstrates the potential of OSIA for enabling biosignature detection and provides insights and lessons learned for future mission concepts aimed at exploring the outer solar system.

Plasma lensing displays interesting characteristics that set it apart from gravitational lensing. The magnetised medium induces birefringence in the two polarisation modes. As the lensing deflection grows stronger, e.g. when images form near the critical curve, the geometric delay of the signal can cause rotation in linear polarisation, in addition to Faraday rotation. This rotation has a frequency dependence to the power of four. We study the geometric rotation of the lensed image in a Gaussian density model and find that it is necessary to take into account the geometric rotation when estimating magnetised media, especially in the under-dense lens. At frequencies of $\sim 1$ GHz or lower, the geometric rotation can dominate. We simulate the flux of lensed images and find that when the image forms near the lensing critical curve, the birefringence can convert the linear polarisation and un-polarisation pulse into a circular mode. The lensing magnification has the potential to increase the probability of detecting such events.

The Sun's magnetic fields play an important role in various solar phenomena. Solar pores are regions of intensified magnetic field strength compared to the surrounding photospheric environment, and their study can help us better understand the properties and behaviour of magnetic fields in the Sun. Up to now, there exists only a single study on magnetic field properties at the boundary region of a pore. Therefore, the main goal of this work is to increase the statistics of magnetic properties determining the pore boundary region. We analyse six solar pores using data from the Helioseismic and Magnetic Imager instrument on board the Solar Dynamics Observatory. We apply image processing techniques to extract the relevant features of the solar pores and determine the boundary conditions of the magnetic fields. We find the maximal vertical magnetic field values on the boundaries of the studied solar pores to range from 1400~G to 1600~G, with a standard deviation between 7.8\% and 14.8\%. These values are lower than those reported in the mentioned preceding study. However, this can be explained by differences in spatial resolution as well as the type of data we used. The vertical magnetic field is an important factor in determining the boundary of solar pores, and it plays a more relevant role than the intensity gradient. The obtained information will be useful for future studies on the formation and evolution of magnetic structures of the Sun. Additionally, this study highlights the importance of high spatial resolution data for the purpose of accurately characterising the magnetic properties of solar pores.

We report the detection of ammonia masers in the non-metastable (6, 3), (7, 5) and (6, 5) transitions, the latter is the first unambiguous maser detection of that transition ever made. Our observations include the first VLBI detection of ammonia maser emission, which allowed effective constrain of the (6, 5) maser brightness temperature. The masers were detected towards G358.931-0.030, a site of 6.7-GHz class~II methanol maser emission that was recently reported to be undergoing a period of flaring activity. These ammonia masers appear to be flaring contemporaneously with the class~II methanol masers during the accretion burst event of G358.931-0.030. This newly detected site of ammonia maser emission is only the twelfth such site discovered in the Milky Way. We also report the results of an investigation into the maser pumping conditions, for all three detected masing transitions, through radiative transfer calculations constrained by our observational data. These calculations support the hypothesis that the ammonia (6, 5) maser transition is excited through high colour temperature infrared emission, with the (6, 5) and (7, 5) transition line-ratio implying dust temperatures >400K. Additionally, we detect significant linearly polarised emission from the ammonia (6, 3) maser line. Alongside our observational and radiative transfer calculation results, we also report newly derived rest frequencies for the ammonia (6, 3) and (6, 5) transitions.

We use our photometric time series of more than forty years to analyze the long-term behaviour of EI Eri. Flare activity is investigated using space-borne photometric data obtained with TESS. The MUSICOS campaign aimed to achieve high-resolution spectroscopic observations from many sites around the globe, so that uninterrupted phase coverage of EI Eri became available. We use these data to reconstruct successive surface temperature maps of the star in order to study the changes of starspots on a very short timescale. We use long-term, seasonal period analysis of our photometric time series to study changes in the rotational period. Short-term Fourier-transform is also applied to look for activity cycle-like changes. We also study the phase and frequency distribution of hand-selected flares. We apply our multi-line Doppler imaging code to reconstruct four consecutive Doppler images. These images are also used to measure surface differential rotation by our cross-correlation technique. In addition, we carry out tests to demonstrate how Doppler imaging is affected by the fact that the data came from several different instruments with different spectral resolutions. Seasonal period analysis of the light curve reveals a smooth, significant change in period, possibly indicating the evolution of active latitudes. Temperature curves from $B-V$ and $V-I$ show slight differences, indicating the activity of EI Eri is spot dominated. Short-term Fourier transform reveals smoothly changing cycles between 4.5--5.5 and 8.9--11.6 years. The time-resolved spotted surface of EI Eri from Doppler imaging enabled us to follow the evolution of the different surface features. Cross-correlating the consecutive Doppler maps reveal surface shear of $\alpha=0.036\pm0.007$. Our tests validate our approach and show that the surface temperature distribution is adequately reconstructed by our method.

In this paper we study the initiation mechanism of the first on-disk X-class eruptive flare in solar cycle 25. Coronal magnetic field reconstructions reveal a magnetic flux rope (MFR) with configuration highly consistent with a filament existing for a long period before the flare, and the eruption of the whole filament indicates that the MFR erupted during the flare. However, quantitative analysis shows that the pre-flare MFR resides in a height too low to trigger a torus instability (TI). The filament experienced a slow rise before the flare onset, for which we estimate evolution of the filament height using a triangulation method by combining the SDO and STEREO observations, and find it is also much lower than the critical height for triggering TI. On the other hand, the pre-flare evolution of the current density shows progressive thinning of a vertical current layer on top of the flare PIL, which suggests that a vertical current sheet forms before the eruption. Meanwhile, there is continuously shearing motion along the PIL under the main branch of the filament, which can drive the coronal field to form such a current sheet. As such, we suggest that the event follows a reconnection-based initiation mechanism as recently established using a high-accuracy MHD simulation, in which an eruption is initiated by reconnection in a current sheet that forms gradually within continuously-sheared magnetic arcade. The eruption should be further driven by TI as the filament quickly rises into the TI domain during the eruption.

We report dual-frequency very long baseline interferometry (VLBI) observations at 22 and 43 GHz toward the nucleus of a nearby radio galaxy NGC 4261. In particular, we present a VLBI image of the 22 GHz H2O maser line and its location in the circumnuclear region of NGC 4261. H2O maser emission is marginally detected above the three times the rms level within a velocity range of approximately 2250-2450 km/s, slightly red-shifted with respect to the systemic velocity. H2O maser emission is located approximately 1 milliarcsecond (mas) east of the brightest continuum component at 22 GHz, where the continuum spectrum is optically thick, that is at the free-free absorbed receding jet by ionized gas. A positional coincidence between H2O maser emission and an ionized gas disk implies that the H2O maser emission arises from the near side of the disk, amplifying continuum emission from the background receding jet. If the disk axis is oriented 64 degree relative to the line of sight, the H2O maser emission is expected to be at a mean radius of 0.3 pc in the disk. The broad line width of the H2O maser emission can be attributed to complex kinematics in the immediate vicinity of the supermassive black hole (SMBH), including ongoing gas infall onto the SMBH, turbulence, and outflow. This is analogous to the multi-phase circumnuclear torus model in the nearest radio-loud H2O megamaser source NGC 1052. An alternative explanation for H2O maser association is the shock region between the jet and the ambient molecular clouds. However, this explanation fails to describe the explicit association of H2O maser emission only with the free-free absorbed receding jet.

The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation^{1,2}. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole^3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of 8.4_{-1.1}^{+0.5} Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow.

We investigate the influence of dust particle size evolution on non-ideal magnetohydrodynamic effects during the collapsing phase of star-forming cores, taking both the turbulence intensity in the collapsing cloud core and the fragmentation velocity of dust particles as parameters. When the turbulence intensity is small, the dust particles do not grow significantly, and the non-ideal MHD effects work efficiently in high-density regions. The dust particles rapidly grow in a strongly turbulent environment, while the efficiency of non-ideal MHD effects in such an environment depends on the fragmentation velocity of the dust particles. When the fragmentation velocity is small, turbulence promotes coagulation growth and collisional fragmentation of dust particles, producing small dust particles. In this case, the adsorption of charged particles on the dust particle surfaces becomes efficient and the abundance of charged particles decreases, making non-ideal MHD effects effective at high densities. On the other hand, when the fragmentation velocity is high, dust particles are less likely to fragment, even if the turbulence is strong. In this case, the production of small dust particles become inefficient and non-ideal MHD effects become less effective. We also investigate the effect of the dust composition on the star and disk formation processes. We constrain the turbulence intensity of a collapsing core and the fragmentation velocity of dust for circumstellar disk formation due to the dissipation of the magnetic field.

We present a detailed temporal analysis for the eclipsing high-mass X-ray binary system Cen X-3 using the Insight-HXMT data in 2018 and 2020. Three consecutive and high statistic observations among data are used for the precise timing analysis. The pulse profiles are revealed to vary with energy and time. The pulse profiles for the 2018 observations showed a double peak in the low energy bands below 10 keV and evolved to a single peak in higher energies without the correlation between pulse fraction and flux, and profiles in low energies changed with time. But the pulse profile for the 2020 observation only showed a broad single-peaked pulse in all energy bands with a positive relation between pulse fraction and flux, which may indicate the transition of the emission patterns from a mixture of a pencil and a fan beam to a dominated pencil-like beam. With performing a binary orbital fitting of spin periods, we obtain an accurate value for the spin period and the orbital parameters. The intrinsic spin period of the neutron star is found to be $4.79920 \pm 0.00006$ s at MJD 58852.697, with the orbital period determined at $P_{\rm orb}=2.08695634\pm 0.00000001$ day, and its decay rate of -(1.7832 $\pm$ 0.0001) $\times$ 10$^{-6}$ yr$^{-1}$ for the binary.

We analyzed two observations of the low-mass black hole X-ray binary MAXI J1348-630 from Nuclear Spectroscopic Telescope Array (NuSTAR) during low hard state and hard intermediate state in the 2019 outburst. The reflection components are found in the X-ray spectra, and the spectral fittings give an inclination angle of $\sim 25^\circ-35^\circ$ for the binary system, and there is an absorption line around $\sim 7$ keV coming from highly ionized iron. The photoionization code XSTAR is used to fit the absorption line, which is attributed to outflows with a velocity of $\sim 10^{4}\rm km\ s^{-1}$ along our line of sight and the column density reaching $10^{23}\rm cm^{-2}$ in low hard and hard intermediate states. The physical mechanism launching fast disk winds from the black hole accretion system is still uncertain. These observations strongly support magnetic launching as the dominant mechanism which drives the high velocity, high ionization winds from the inner accretion disk region in hard and hard intermediate states of MAXI J1348-630.

We present the results of the neutron star X-ray binary system Cen X-3 performed by $Insight$-HXMT with two observations during 2017 and 2018. During these two observations, the source reached a X-ray luminosity of $\sim 10^{38}$ erg s$^{-1}$ from 2 -- 105 keV. The analysis of the broadband X-ray spectrum reports the presence of two cyclotron resonance scattering features (CRSFs) with the fundamental line at $\sim$ 28 keV and the harmonic line at $\sim 47 $ keV. The multiple lines exist by fittings with different continuum models, like the absorbed NPEX model and a power-law with high energy exponential cutoff model. This is the first time that both fundamental and harmonic lines are detected in Cen X-3. We also show evidence of two cyclotron lines in the phase-resolved spectrum of Cen X-3. The CRSF and continuum spectral parameters show evolution with the pulse profile, and the two line centroid energy ratio does not change significantly and locates in a narrow value range of $1.6-1.7$ over the pulse phase. The implications of the discovering two cyclotron absorption features and phase-resolved spectral properties are discussed.

MAXI J1348-630 is a low-mass X-ray black hole binary located in the Galaxy and undergone the X-ray outburst in 2019. We analyzed the observation data in very soft state during the outburst between MJD 58588 and MJD 58596 based on the Insight-HXMT observations from 2 -- 20 keV via the continuum fitting method to measure the spin of the stellar-mass black hole in MAXI J1348-630. The inner disk temperature and the apparent inner disk radius were found to be $0.47\pm 0.01 \rm keV$ and $5.33\pm 0.10 \ R_{g}$ from the observation data modeled by the multicolor disc blackbody model. Assuming the distance of the source $D\sim 3.4 \rm kpc$, the mass of the black hole $M\sim 11 \ M_{\odot}$, and the inclination of the system $i\sim 29.2^{\circ}$, the spin is determined to be $a_{\star}=0.41\pm 0.03$ for fixing hardening factor at 1.6 and $n_{H}=8.6\times 10^{21} \rm cm^{-2}$. Besides, considering the uncertainty of the parameters $D, M, i$ of this system, with the Monte Carlo analysis, we still confirm the moderate spin of the black hole as $a_{\star}=0.42^{+0.13}_{-0.50}$. Some spectral parameters (e.g., column density and hardening factor) which could affect the measurements of the BH spin are also briefly discussed.

The integrated perturbation theory (iPT) is a set of methods in nonlinear perturbation theory for the structure formation in the Universe. In Papers~I and II, the basic formalism and technical methods of the iPT for cosmological tensor fields are developed, generalizing the corresponding theory for scalar fields. In previous papers, methods to predict statistical quantities, such as power spectra, correlation functions, etc., of three-dimensional tensor fields are developed based on the iPT. However, observations of tensors, such as angular momenta and shapes of galaxies, etc., are only possible after the three-dimensional tensors are projected onto the two-dimensional sky. In this paper, power spectra and correlation functions of projected two-dimensional tensors are related to those of original three-dimensional tensors, so that one can make predictions for the observable statistics of projected tensor fields from the iPT. The relations are consistently represented on the basis of irreducible decomposition of both two- and three-dimensional tensors.

Supermassive black holes (SMBHs) and nuclear star clusters (NSCs) co-exist in many galaxies. While the formation history of the black hole is essentially lost, NSCs preserve their evolutionary history imprinted onto their stellar populations and kinematics. Studying SMBHs and NSCs in tandem might help us to ultimately reveal the build-up of galaxy centres. In this study, we combine large-scale VLT/MUSE and high-resolution adaptive-optics-assisted VLT/SINFONI observations of the early-type galaxy FCC 47 with the goal being to assess the effect of a spatially (non-)variable initial mass function (IMF) on the determination of the mass of the putative SMBH in this galaxy. We achieve this by performing DYNAMITE Schwarzschild orbit-superposition modelling of the galaxy and its NSC. In order to properly take account of the stellar mass contribution to the galaxy potential, we create mass maps using a varying stellar mass-to-light ratio derived from single stellar population models with fixed and with spatially varying IMFs. Using the two mass maps, we estimate black hole masses of $(7.1^{+0.8}_{-1.1})\times 10^7\,M_{\odot}$ and $(4.4^{+1.2}_{-2.1}) \times 10^7\,M_{\odot}$ at $3\sigma$ signifance, respectively. Compared to models with constant stellar-mass-to-light ratio, the black hole masses decrease by 15% and 48%, respectively. Therefore, a varying IMF, both in its functional form and spatially across the galaxy, has a non-negligible effect on the SMBH mass estimate. Furthermore, we find that the SMBH in FCC 47 has probably not grown over-massive compared to its very over-massive NSC.

We found that SDSS J094002.56+274942.0 underwent a superoutburst in 2019 February based on our observations and Zwicky Transient Facility (ZTF) data. This object showed shallow eclipses during this superoutburst and we established the orbital period to be 0.1635015(1) d in combination with the ZTF and Asteroid Terrestrial-impact Last Alert System (ATLAS) data in quiescence. Superhumps apparently started to develop soon after the object reached the plateau phase and fully grown superhumps were recorded within the initial 6 d of the plateau phase. Using the superhump and orbital periods, we obtained a mass ratio (q) of 0.39(3) and obtained an inclination of 70.5(5) deg by eclipse modeling. These values reproduced the quiescent ellipsoidal variations very well. Using the Gaia parallax and 2MASS observations, we confirmed that the secondary is indistinguishable from an unevolved main-sequence star. The resultant mass ratio and orbital period were the highest among SU UMa stars, and this provided a proof that the 3:1 resonance can develop in less than 6 d even in q=0.39(3). The superoutburst faded relatively rapidly and was followed by a rebrightening, suggesting that the tidal effect in a large-q system was insufficient to maintain a long superoutburst and the remnant matter caused a rebrightening. The presence of such a system among dwarf novae is against the conventional idea that outbursts in dwarf novae are not long enough to develop superhumps, in contrast to novalike variables, under a weak tidal effect. The present observation also supports that the 3:1 resonance is the cause of a long outburst, and not its consequence, even under extreme q. The rapid growth of the 3:1 resonance in a high-q system challenges the generally accepted results of hydrodynamic simulations.

Swift J1357.2$-$0933 is a black-hole candidate X-ray transient, which underwent its third outburst in 2019, during which several multi-wavelength observations were carried out.~Here, we report results from the \emph{Neil Gehrels Swift} and \emph{NICER} observatories and radio data from \emph{AMI}.~For the first time,~millihertz quasi-periodic X-ray oscillations with frequencies varying between ${\sim}$~1--5~$\rm{mHz}$ were found in \emph{NICER} observations and a similar feature was also detected in one \emph{Swift}--\textsc{XRT} dataset.~Our spectral analysis indicate that the maximum value of the measured X-ray flux is much lower compared to the peak values observed during the 2011 and 2017 outbursts.~This value is ${\sim}$~100 times lower than found with \emph{MAXI} on MJD~58558 much ($\sim$~68 days) earlier in the outburst, suggesting that the \emph{Swift} and \emph{NICER} fluxes belong to the declining phase of the 2019 outburst.~An additional soft component was detected in the \textsc{XRT} observation with the highest flux level, but at a relatively low $L_{\rm X}$~$\sim$~$3{\times}10^{34}~(d/{\rm 6~kpc)}^2\rm{erg}~\rm{s}^{-1}$, and which we fitted with a disc component at a temperature of $\sim 0.17$~keV.~The optical/UV magnitudes obtained from \emph{Swift}--\textsc{UVOT} showed a correlation with X-ray observations, indicating X-ray reprocessing to be the plausible origin of the optical and UV emission.~However, the source was not significantly detected in the radio band.~There are currently a number of models that could explain this millihertz-frequency X-ray variability; not least of which involves an X-ray component to the curious dips that, so far, have only been observed in the optical.

Recent infrared and submillimeter observations suggest that the protoplanetary disk lifetime depends on the central stellar mass. The disk dispersal is thought to be driven by viscous accretion, magneto-hydrodynamics (MHD) winds, and photoevaporation by the central star. We perform a set of one-dimensional simulations of long-term disk evolution that include all the three processes. We vary the stellar mass in the range of 0.5-7M$_{\odot}$, and study the mass dependence of the disk evolution. We show that a significant fraction of the disk gas is lost by MHD winds in the early stage, but the later disk evolution is mainly governed by photoevaporation. The disk radius decreases as photoevaporation clears out the gas in the outer disk efficiently. The qualitative evolutionary trends of the disk mass are remarkably similar for the wide range of the central stellar mass we consider, and the time evolution of the disk mass can be well fitted by a simple function. The dispersal time is approximately ten million years for low mass stars with weak mass dependence, but gets as short as two million years around a 7M$_{\odot}$ star. In the latter case, a prominent inner hole is formed by the combined effect of accretion and MHD winds within about one million years. The strength of the MHD wind and viscous accretion controls the overall mass-loss rate, but does not alter the dependence of the dispersal timescale on the central stellar mass.

The solar magnetic structure changes over the solar cycle. It has a dipole structure during solar minimum, where the open flux extends mainly from the polar regions into the interplanetary space. During maximum, a complex structure is formed with low-latitude active regions and weakened polar fields, resulting in spread open field regions. However, the components of the solar magnetic field that is responsible for long-term variations in the interplanetary magnetic field (IMF) are not clear, and the IMF strength estimated based on the solar magnetic field is known to be underestimated by a factor of 3 to 4 against the actual in-situ observations (the open flux problem). To this end, we decomposed the coronal magnetic field into the components of the spherical harmonic function of degree and order $(\ell, m)$ using the potential field source surface model with synoptic maps from SDO/HMI for 2010 to 2021. As a result, we found that the IMF rapidly increased in December 2014 (seven months after the solar maximum), which coincided with the increase in the equatorial dipole, $(\ell, m)=(1, \pm1)$, corresponding to the diffusion of active regions toward the poles and in the longitudinal direction. The IMF gradually decreased until December 2019 (solar minimum) and its variation corresponded to that of the non-dipole component $\ell\geq2$. Our results suggest that the understanding of the open flux problem may be improved by focusing on the equatorial dipole and the non-dipole component and that the influence of the polar magnetic field is less significant.

The observation of varying broad iron lines during the state transition of the black hole X-ray binaries (BHXBs) have been accumulating.In this work, the relation between the normalized intensity and the width of iron lines is investigated, in order to understand better the variation of iron lines and possibly its connection to state transition. Considering the uncertainties due to ionization and illuminating X-rays, only the effects of geometry and gravity are taken into account. Three scenarios were studied, i.e., the continuous disk model, innermost annulus model, and the cloud model. As shown by our calculations, at given iron width, the line flux of the cloud model is smaller than that of the continuous disk model; while for the innermost annulus model, the width is almost unrelated with the flux. The range of the line strength depends on both the BH spin and the inclination of the disk. We then apply to the observation of MAXI J1631-479 by NuSTAR during its decay from the soft state to the intermediate state. We estimated the relative line strength and width according to the spectral fitting results by Xu et al.(2020), and then compared with our theoretical width-flux relation. It was found that the cloud model was more favored. We further modeled the iron line profiles, and found that the cloud model can explain both the line profile and its variation with reasonable parameters.

In this study we present an analysis of the performance and properties of the quasi-periodic (QP) GP kernel, which is the multiplication of the squared-exponential kernel by the exponential-sine-squared kernel, based on an extensive set of synthetic RVs, into which the signature of activity was injected. We find that while the QP-GP rotation parameter matches the simulated rotation period of the star, the length scale cannot be directly connected to the spot lifetimes on the stellar surface. Regarding the setup of the priors for the QP-GP, we find that it can be advantageous to constrain the QP-GP hyperparameters in different ways depending on the application and the goal of the analysis. We find that a constraint on the length scale of the QP-GP can lead to a significant improvement in identifying the correct rotation period of the star, while a constraint on the rotation hyperparameter tends to lead to improved planet detection efficiency and more accurately derived planet parameters. Even though for most of the simulations the Bayesian evidence performed as expected, we identified not far-fetched cases where a blind adoption of this metric would lead to wrong conclusions. We conclude that modeling stellar astrophysical noise by using a QP-GP considerably improves detection efficiencies and leads to precise planet parameters. Nevertheless, there are also cases in which the QP-GP does not perform optimally, for example RV variations dynamically evolving on short timescales or a mixture of a very stable activity component and random variations. Knowledge of these limitations is essential for drawing correct conclusions from observational data.

A striking feature of the Vela pulsar (PSR J0835$-$4510) is that it undergoes sudden increases in its spin frequency, known as glitches, with a fractional amplitude on the order of $10^{-6}$ approximately every 900 days. Glitches of smaller magnitudes are also known to occur in Vela. Their distribution in both time and amplitude is less well constrained but equally important for understanding the physical process underpinning these events. In order to better understand these small glitches in Vela, an analysis of high-cadence observations from the Mount Pleasant Observatory is presented. A hidden Markov model (HMM) is used to search for small, previously undetected glitches across 24 years of observations covering MJD 44929 to MJD 53647. One previously unknown glitch is detected around MJD 48636 (Jan 15 1992), with fractional frequency jump $\Delta f/f = (8.19 \pm 0.04) \times 10^{-10}$ and frequency derivative jump $\Delta\dot{f}/\dot{f} = (2.98 \pm 0.01) \times 10^{-4}$. Two previously reported small glitches are also confidently re-detected, and independent estimates of their parameters are reported. Excluding these events, 90% confidence frequentist upper limits on the sizes of missed glitches are also set, with a median upper limit of $\Delta f^{90\%}/f = 1.35 \times 10^{-9}$. Upper limits of this kind are enabled by the semi-automated and computationally efficient nature of the HMM, and are crucial to informing studies which are sensitive to the lower end of the glitch size distribution.

We have performed experimental investigations of methanol formation via the reactions of low energy $CH_{3}^{+}$ ions with an amorphous solid water (ASW) surface around 10 K. A newly developed experimental apparatus enabled irradiation of the ASW surface by several eV ions and detection of trace amounts of reaction products on the surface. It was found that methanol molecules were produced by low-energy $CH_{3}^{+}$ irradiation of the ASW surface and that hydroxy groups in produced methanol originated from water molecules in ASW, as predicted in a previous theoretical study. Little temperature dependence of observed methanol intensity is apparent in the temperature range 12 - 60 K. Ab-initio molecular dynamics simulations under constant temperature conditions of 10 K suggested that this reaction spontaneously produced a methanol molecule and an $H_{3}O^{+}$ ion, regardless of the contact point of $CH_{3}^{+}$ on the ASW surface. We have performed simulation with an astrochemical model under molecular-cloud conditions, where the reaction between $CH_{3}^{+}$ and $H_{2}O$ ice, leading to methanol formation, was included. We found that the impact of the reaction on methanol abundance was limited only at the edge of the molecular cloud (< 1 mag) because of the low abundance of $CH_{3}^{+}$ in the gas phase, whereas the reaction between the abundant molecular ion $HCO^{+}$ and $H_{2}O$ ice, which has not yet been confirmed experimentally, can considerably affect the abundance of a complex organic molecule. This work sheds light on a new type of reaction between molecular ions and ice surfaces that should be included in astrochemical models.

Context: FUors outbursts are a crucial stage of accretion in young stars. However a complete mechanism at the origin of the outburst still remains missing. Aims: We aim at constraining the instability mechanism in FU Orionis star itself, by directly probing the size and the evolution in time of the outburst region with near-infrared interferometry, and to confront it to physical models of this region. Methods: FU Orionis has been a regular target of near-infrared interferometry. In this paper, we analyze more than 20 years of interferometric observations to perform a temporal monitoring of the region of the outburst, and compare it to the spatial structure deduced from 1D MHD simulations. Results: We measure from the interferometric observations that the size variation of the outburst region is compatible with a constant or slightly decreasing size over time in the H and K band. The temporal variation and the mean sizes are consistently reproduced by our 1D MHD simulations. We find that the most compatible scenario is a model of an outburst occurring in a magnetically layered disk, where a Magneto-Rotational Instability (MRI) is triggered by a Gravitational Instability (GI) at the outer edge of a dead-zone. The scenario of a pure Thermal Instability (TI) fails to reproduce our interferometric sizes since it can only be sustained in a very compact zone of the disk <0.1 AU. The scenario of MRI-GI could be compatible with an external perturbation enhancing the GI, such as tidal interactions with a stellar companion, or a planet at the outer edge of the dead-zone. Conclusions: The layered disk model driven by MRI turbulence is favored to interpret the spatial structure and temporal evolution of FU Orionis outburst region. Understanding this phase gives a crucial link between the early phase of disk evolution and the process of planet formation in the first inner AUs.

We estimated the angular momentum, $J$, of $720$ galaxies in the Local Volume with distances $r < 12$ Mpc. The distribution of the average angular momentum along the Hubble sequence has a maximum at the morphological type $T= 4$ (Sbc), while the dispersion of the $J$-values for galaxies is minimal. Among the Local Volume population, 27 elite spiral galaxies stand out, with an angular momentum greater than 0.15 of the Milky Way, $J > 0.15 J_{MW}$, making the main contribution ($ > 90\%$) to the total angular momentum of galaxies in the considered volume. Using observational data on the kinematics and structure of these galaxies, we determined the direction of their spins. We present the first map of the distribution of the spins of 27 nearby massive spiral galaxies in the sky and note that their pattern does not exhibit significant alignment with respect to the Local Sheet plane. The relationship between the magnitude of the angular momentum and stellar mass of the local galaxies is well represented by a power law with an exponent of ($5/3$) over an interval of $6$ orders of magnitude of the mass of galaxies.

The magneto-rotational instability (MRI) - which is due to an interplay between a sheared background and the magnetic field - is commonly considered a key ingredient for developing and sustaining turbulence in the outer envelope of binary neutron star merger remnants. To assess whether (or not) the instability is active and resolved, criteria originally derived in the accretion disk literature - thus exploiting the symmetries of such systems - are often used. In this paper we discuss the magneto-shear instability as a truly local phenomenon, relaxing common symmetry assumptions on the background on top of which the instability grows. This makes the discussion well-suited for highly dynamical environments such as binary mergers. We find that - although this is somewhat hidden in the usual derivation of the MRI dispersion relation - the instability crucially depends on the assumed symmetries. Relaxing the symmetry assumptions on the background we find that the role of the magnetic field is significantly diminished, as it affects the modes' growth but does not drive it. This suggests that we should not expect the standard instability criteria to provide a faithful indication/diagnostic of what "is actually going on" in mergers. We conclude by making contact with a suitable filtering operation, as this is key to separating background and fluctuations in highly dynamical systems.

In this work, we present a study of the propagation of low-frequency magneto-acoustic waves into the solar chromosphere within small-scale inclined magnetic fields over a quiet-magnetic network region utilizing near-simultaneous photospheric and chromospheric Dopplergrams obtained from the HMI instrument onboard SDO spacecraft and the Multi-Application Solar Telescope (MAST) operational at the Udaipur Solar Observatory, respectively. Acoustic waves are stochastically excited inside the convection zone of the Sun and intermittently interact with the background magnetic fields resulting into episodic signals. In order to detect these episodic signals, we apply the wavelet transform technique to the photospheric and chromospheric velocity oscillations in magnetic network regions. The wavelet power spectrum over photospheric and chromospheric velocity signals show a one-to-one correspondence between the presence of power in the 2.5-4 mHz band. Further, we notice that power in the 2.5-4 mHz band is not consistently present in the chromospheric wavelet power spectrum despite its presence in the photospheric wavelet power spectrum. This indicates that leakage of photospheric oscillations (2.5-4 mHz band) into the higher atmosphere is not a continuous process. The average phase and coherence spectra estimated from these photospheric and chromospheric velocity oscillations illustrate the propagation of photospheric oscillations (2.5-4 mHz) into the solar chromosphere along the inclined magnetic fields. Additionally, chromospheric power maps estimated from the MAST Dopplergrams also show the presence of high-frequency acoustic halos around relatively high magnetic concentrations, depicting the refraction of high-frequency fast mode waves around vA ~ vs layer in the solar atmosphere.

We present near-infrared spectroscopy of the 2022 eruption of the recurrent nova U Sco, over the period from 5.2 to 45.4 days after outburst. This is the most intensive infrared study of this nova. Our observations started early after the outburst and extended almost to the end of the ``Super Soft'' X-ray phase. A major find is the presence of coronal lines from day 9.41, one of the earliest appearances of these in any nova, classical or recurrent. The temperature of the coronal gas is $7\times10^5$ K. There is also evidence for the presence of much cooler ($\lesssim2.5\times10^4$ K) gas. Remarkable changes are seen in the HeI $1.083\mu$m line, the strength of which declines, then recovers, in anti-correlation with the X-ray behaviour. We conclude that shock ionisation is the dominant excitation mechanism for the coronal line emission. There is evidence in the infrared spectra for the presence of black body emission at $\sim20000$ K, which we tentatively identify with the irradiated secondary, and for free-free/free-bound emission. For the previously determined binary inclination of $82.7$ degrees, the implied ejection velocities are as high as 22000 km s$^{-1}$. These velocities appear unprecedented in nova outflows, and are comparable to those seen in supernovae, thereby marking U Sco as a truly remarkable object.

Luminous infrared galaxies at high redshifts ($z$>4) include extreme starbursts that build their stellar mass over short periods of time (>100 Myr). These galaxies are considered to be the progenitors of massive quiescent galaxies at intermediate redshifts ($z\sim$2) but their stellar structure and buildup is unknown. Here, we present the first spatially resolved near-infrared imaging of GN20, one of the most luminous dusty star-forming galaxies known to date, observed at an epoch when the Universe was only 1.5 Gyr old. The 5.6$\mu$m image taken with the JWST Mid-Infrared Instrument (MIRI/JWST) shows that GN20 is a very luminous galaxy (M$_\mathrm{1.1\mu m,AB}$=$-$25.01), with a stellar structure composed of a conspicuous central source and an extended envelope. The central source is an unresolved nucleus that carries 9% of the total flux. The nucleus is co-aligned with the peak of the cold dust emission, and offset by 3.9 kpc from the ultraviolet stellar emission. The diffuse stellar envelope is similar in size to the clumpy CO molecular gas distribution. The centroid of the stellar envelope is offset by 1 kpc from the unresolved nucleus, suggesting GN20 is involved in an interaction or merger event supported by its location as the brightest galaxy in a proto-cluster. The stellar size of GN20 is larger by a factor of about 3-5 than known spheroids, disks, and irregulars at $z\sim$4, while its size and low S\'ersic index are similar to those measured in dusty, infrared luminous galaxies at $z\sim$2 of the same mass. GN20 has all the ingredients necessary for evolving into a massive spheroidal quiescent galaxy at intermediate $z$: it is a large, luminous galaxy at $z$=4.05 involved in a short and massive starburst centred in the stellar nucleus and extended over the entire galaxy, out to radii of 4 kpc, and likely induced by the interaction or merger with a member of the proto-cluster.

Decayless kink oscillations are omnipresent in the solar atmosphere and a viable candidate for coronal heating. Though there have been extensive studies of decayless oscillations in coronal loops with a few hundred Mm lengths, the properties of these oscillations in small-scale ($\sim$10 Mm) loops are yet to be explored. In this study, we present the properties of decayless oscillations in small loops embedded in the quiet corona and coronal holes. We use high resolution observations from the Extreme Ultraviolet Imager onboard Solar Orbiter with pixel scales of 210 km and 5 s cadence or better. We find 42 oscillations in 33 coronal loops with loop lengths varying between 3 to 23 Mm. The average displacement amplitude is found to be 136 km. The oscillations period has a range of 27 to 276 s, and the velocity amplitudes range from 2.2 to 19.3 km s$^{-1}$. The observed kink speeds are lower than those observed in active region coronal loops. The variation of loop length with the period does not indicate a strong correlation. Coronal seismology technique indicated an average magnetic field value of 2.1 G. We estimate the energy flux with a broad range of 0.6-314 W m$^{-2}$. Moreover, we note that the short-period decayless oscillations are not prevalent in the quiet Sun and coronal holes. Therefore, our study suggests that decayless oscillations in small-scale coronal loops are unlikely to provide enough energy to heat the quiet Sun and accelerate solar wind in the coronal holes.

We cross-match and compare characteristics of galaxy clusters identified in observations from two sky surveys using two completely different techniques. One sample is optically selected from the analysis of three years of Dark Energy Survey observations using the redMaPPer cluster detection algorithm. The second is X-ray selected from XMM observations analysed by the XMM Cluster Survey. The samples comprise a total area of 57.4 deg$^2$, bounded by the area of 4 contiguous XMM survey regions that overlap the DES footprint. We find that the X-ray selected sample is fully matched with entries in the redMaPPer catalogue, above $\lambda>$20 and within 0.1$< z <$0.9. Conversely, only 38\% of the redMaPPer catalogue is matched to an X-ray extended source. Next, using 120 optically clusters and 184 X-ray selected clusters, we investigate the form of the X-ray luminosity-temperature ($L_{X}-T_{X}$), luminosity-richness ($L_{X}-\lambda$) and temperature-richness ($T_{X}-\lambda$) scaling relations. We find that the fitted forms of the $L_{X}-T_{X}$ relations are consistent between the two selection methods and also with other studies in the literature. However, we find tentative evidence for a steepening of the slope of the relation for low richness systems in the X-ray selected sample. When considering the scaling of richness with X-ray properties, we again find consistency in the relations (i.e., $L_{X}-\lambda$ and $T_{X}-\lambda$) between the optical and X-ray selected samples. This is contrary to previous similar works that find a significant increase in the scatter of the luminosity scaling relation for X-ray selected samples compared to optically selected samples.

Microscopic horizonless relics could form in the early universe either directly through gravitational collapse or as stable remnants of the Hawking evaporation of primordial black holes. In both cases they completely or partially evade cosmological constraints arising from Hawking evaporation and in certain mass ranges can explain the entirety of the dark matter. We systematically explore the stochastic gravitational-wave background associated with the formation of microscopic dark-matter relics in various scenarios, adopting an agnostic approach and discussing the limitations introduced by existing constraints, possible ways to circumvent the latter, and expected astrophysical foregrounds. Interestingly, this signal is at most marginally detectable with current interferometers but could be detectable by third-generations instruments such as the Einstein Telescope, strengthening their potential as discovery machines.

Since the mid 70ies it is known that the dwarf galaxies around the Milky Way are arranged in a thin, polar structure. The arrangement and motion within this structure has been identified as a severe challenge to the standard model of cosmology, dubbed as the plane of satellites problem. More observational evidence for such structures has been put forward around other galaxies, such as the Andromeda galaxy, Cen\,A or NGC\,253, among others, adding to the previously identified tensions. Solutions to the plane of satellite problem should therefore not only be tailored to the Milky Way, but need to explain all these different observed systems and environments.

This work focuses on the identification of reliable and repeatable spatial (three-dimensional) trajectories that link the Earth and the Moon. For this purpose, this paper aims to extend the 2:1 resonant prograde family and 2:1 resonant retrograde family to three dimensions and to introduce spatial orbits that are not currently present in the literature. These orbits, named the 2:1 resonant spatial family, bifurcate from the two-dimensional families and smoothly transition between them in phase space. The stability properties of this new family of resonant orbits are discussed, and, interestingly, this family includes marginally stable members. Furthermore, this new family of orbits is applied to several engineering problems in the Earth-Moon system. First, this paper selects an appropriate member of 2:1 resonant spatial family on the basis of its stability properties and relationships with other multibody orbits in the regime. Next, this work combines this trajectory with momentum exchange tethers to transit payloads throughout the system in a reliable and repeatable fashion. Finally, this paper studies the process of aborting a catch and related recovery opportunities.

Mergers between galaxy clusters are an important stage in the formation of the large-scale structure of the Universe. Some of the mergers show a spectacular bow shock that formed as a result of recent passage of a smaller cluster through a bigger one, the classic example of this being the so-called bullet cluster. In this paper, I describe ten examples of interacting clusters identified among 200 of the most massive objects, with total masses above $1.4 \times 10^{14}$ M$_{\odot}$, from the IllustrisTNG300 simulation by searching for prominent bow shocks in their temperature maps. Despite different mass ratios of the two merging clusters, the events are remarkably similar in many respects. In all cases, the companion cluster passed close to the main one only once, between 0.9 and 0.3 Gyr ago, with the pericenter distance of 100-530 kpc and a velocity of up to 3400 km s$^{-1}$. The subcluster, typically an order of magnitude smaller in mass than the main cluster before the interaction, loses most of its dark matter and gas in the process. The displacement between the collisionless part of the remnant and the bow shock is such that the remnant typically lags behind the shock or coincides with it, with a single exception of the merger occurring with the largest velocity. Usually about 1% of the gas cells in the merging clusters are shocked, and the median Mach numbers of these gas cells are around two. Due to the relatively small size of the simulation box, no close analog of the bullet cluster was found, but I identified one case that is similar in terms of mass, velocity, and displacement. The presented cases bear more resemblance to less extreme observed interacting clusters such as A520 and Coma.

The X-ray pulsar 1A 0535+262 exhibited a giant outburst in 2020, offering us a unique opportunity for X-ray polarimetry of an accreting pulsar in the supercritical state. Measurement with PolarLight yielded a non-detection in 3-8 keV; the 99% upper limit of the polarization fraction (PF) is found to be 0.34 averaged over spin phases, or 0.51 based on the rotating vector model. No useful constraint can be placed with phase resolved polarimetry. These upper limits are lower than a previous theoretical prediction of 0.6-0.8, but consistent with those found in other accreting pulsars, like Her X-1, Cen X-3, 4U 1626-67, and GRO J1008-57, which were in the subcritical state, or at least not confidently in the supercritical state, during the polarization measurements. Our results suggest that the relatively low PF seen in accreting pulsars cannot be attributed to the source not being in the supercritical state, but could be a general feature.

We carried out a new attempt to check for the presence promethium lines in the spectrum of HD101065 (Przybylski's star). The neutron capture element promethium does not have stable isotopes and the maximum half-life time is about 18 years. Thus its presence in this peculiar star would indicate an ongoing process of irradiation of its surface layers with free neutrons. Unfortunately, almost all promethium lines are heavily blended with lines of other neutron capture elements and other species. We selected and analysed three lines of promethium (Pm I and Pm II) and came to the conclusion that at present it is impossible to definitely claim the presence of this element in Przybylski's star atmosphere.

A conducting cylinder, with a uniform magnetic field along its axis, and radial temperature gradient, is considered. At large temperature gradients the azimuthal Hall electrical current creates the axial magnetic field which strength may be comparable with the original one. It is shown, that the magnetic field generated by the azimuthal Hall current leads to the decrease of magnetic field originated by external sources.

We are in the era of precision cosmology which offers us a unique opportunity to investigate beyond standard model physics. Towards this endeavor, inflaton is assumed to be a perfect new physics candidate. In this submission, we explore the phenomenological impact of the latest observation of PLANCK and BICEP/$Keck$ data on the physics of inflation. We particularly study three different models of inflation, namely $\alpha$-attractor E, T, and the minimal plateau model. We further consider two different post-inflationary reheating dynamics driven by inflaton decaying into Bosons and Fermions. Given the latest data in the inflationary $(n_s-r)$ plane, we derive detailed phenomenological constraints on different inflaton parameters and the associated physical quantities, such as inflationary e-folding number, $N_{ k}$, reheating temperatures $T_{\rm re}$. Apart from considering direct observational data, we further incorporate the bounds from primordial gravitational waves (PGWs) and different theoretical constraints. Rather than in the laboratory, our results illustrate the potential of present and future cosmological observations to look for new physics in the sky.

This paper presents a spectral and photometric study of the poorly studied polar BM CrB. Three states of the polar brightness and signs of transition from one-pole to two-pole accretion mode were found by an analysis of ZTF data. It is shown that the transition from the low state to the high state changes the longitude of the main accretion spot (by $\approx 17^{\circ}$) and increases its elongation (by $\approx 10^{\circ}$). The spectra contain Zeeman absorptions of the H$\alpha$ line which are formed at a magnetic field strength of $15.5\pm1$ MG. These absorptions are likely produced by a cold halo extending from the accretion spot at $\approx {^1/_4}$ of the white dwarf radius. Modeling of the behavior of the H$\alpha$ emission line shows that the main source of emission is the part of the accretion stream near the Lagrangian point L$_1$, which is periodically eclipsed by the donor star. The spectra exhibit a cyclotron component formed in the accretion spot. Its modeling by a simple accretion spot model gives constraints on the magnetic field strength $B=15-40$ MG and the temperature $T_e\gtrsim15$ keV.

Habitable rocky planets around M dwarfs that have H2-dominated atmospheres, if they exist, would permit characterizing habitable exoplanets with detailed spectroscopy using JWST, owing to their extended atmospheres and small stars. However, the H2-dominated atmospheres that are consistent with habitable conditions cannot be too massive, and a moderate-size H2-dominated atmosphere will lose mass to irradiation-driven atmospheric escape on rocky planets around M dwarfs. We evaluate volcanic outgassing and serpentinization as two potential ways to supply H2 and form a steady-state H2-dominated atmosphere. For rocky planets of 1-7 Earth mass and early, mid, and late M dwarfs, the expected volcanic outgassing rates from a reduced mantle fall short of the escape rates by >~1 order of magnitude, and a generous upper limit of the serpentinization rate is still less than the escape rate by a factor of a few. Special mechanisms that may sustain the steady-state H2-dominated atmosphere include direct interaction between liquid water and mantle, heat-pipe volcanism from a reduced mantle, and hydrodynamic escape slowed down by efficient upper-atmospheric cooling. It is thus unlikely to find moderate-size, H2-dominated atmospheres on rocky planets of M dwarfs that would support habitable environments.

Photodissociation regions (PDRs), where the (far-)ultraviolet light from hot young stars interact with the gas in surrounding molecular clouds, provide laboratories for understanding the nature and role of feedback by star formation on the interstellar medium. While the general nature of PDRs is well understood - at least under simplified conditions - the detailed dynamics and chemistry of these regions, including gas clumping, evolution over time etc. can be very complex. We present interferometric observations of the 21 cm atomic hydrogen line, combined with [CII] 158 $\mu$m observations, towards the nearby reflection nebula IC 63. We find a clumpy HI structure in the PDR, and a ring morphology for the HI emission at the tip of IC 63. We further unveil kinematic substructure, of the order of 1~km~s$^{-1}$, in the PDR layers and several legs that will disperse IC 63 in $<$0.5 Myr. We find that the dynamics in the PDR explain the observed clumpy HI distribution and lack of a well-defined HI/H$_{2}$ transition front. However, it is currently not possible to conclude whether HI self-absorption (HISA) and non-equilibrium chemistry also contribute to this clumpy morphology and missing HI/H$_{2}$ transition front.

Aims. We investigate dust attenuation of 122 heavily dust-obscured galaxies detected with the Atacama Large Millimeter Array (ALMA) and Herschel in the COSMOS field. We search for correlations between dust attenuation recipes and the variation of physical parameters, mainly the effective radii of galaxies, their star formation rates (SFR), and stellar masses, and aim to understand which of the commonly used laws best describes dust attenuation in dusty star-forming galaxies at high redshift. Methods. We make use of the extensive photometric coverage of the COSMOS data combined with highly-resolved dust continuum maps from ALMA. We use CIGALE to estimate various physical properties of these dusty objects, mainly their SFR, their stellar masses and their attenuation. We infer galaxy effective radii (Re) using GALFIT in the Y band of HSC and ALMA continuum maps. We use these radii to investigate the relative compactness of the dust continuum and the extension of the rest-frame UV/optical Re(y)/Re(ALMA). Results. We find that the physical parameters calculated from our models strongly depend on the assumption of dust attenuation curve. As expected, the most impacted parameter is the stellar mass, which leads to a change in the "starburstiness" of the objects. We find that taking into account the relative compactness of star-to-dust emission prior to SED fitting is crucial, especially when studying dust attenuation of dusty star-forming galaxies. Shallower attenuation curves did not show a clear preference of compactness with attenuation, while the Calzetti attenuation curve preferred comparable spatial extent of unattenuated stellar light and dust emission. The evolution of the Re(UV)/Re(ALMA) ratio with redshift peeks around the cosmic noon in our sample of DSFGs, showing that this compactness is correlated with the cosmic SFR density of these dusty sources.

Most massive stars are believed to be born in close binary systems where they can exchange mass, which impacts the evolution of both binary components. Their evolution is of great interest in the search for the progenitors of gravitational waves. However, there are unknowns in the physics of mass transfer as observational examples are rare, especially at low metallicity. Nearby low-metallicity environments are particularly interesting hunting grounds for interacting systems as they act as the closest proxy for the early universe where we can resolve individual stars. Using multi-epoch spectroscopic data, we complete a consistent spectral and orbital analysis of the early-type massive binary SSN~7 hosting a ON3\,If$^\ast$+O5.5\,V((f)) star. Using these detailed results, we constrain an evolutionary scenario that can help us to understand binary evolution in low metallicity.} We were able to derive reliable radial velocities of the two components from the multi-epoch data, which were used to constrain the orbital parameters. The spectroscopic data covers the UV, optical, and near-IR, allowing a consistent analysis with the stellar atmosphere code, PoWR. Given the stellar and orbital parameters, we interpreted the results using binary evolutionary models. The two stars in the system have comparable luminosities of ${\log (L_1/L_{\odot}) = 5.75}$ and ${\log (L_2/L_{\odot}) = 5.78}$ for the primary and secondary, respectively, but have different temperatures (${T_1=43.6\,\mathrm{kK}}$ and ${T_2=38.7\,\mathrm{kK}}$). The primary ($32\,M_{\odot}$) is less massive than the secondary ($55\,M_{\odot}$), suggesting mass exchange. The mass estimates are confirmed by the orbital analysis. The revisited orbital period is $3\,\mathrm{d}$. Our evolutionary models also predict mass exchange. Currently, the system is a contact binary undergoing a slow Case A phase, making it the most massive [Abridged]

We present an analysis of the ultraviolet luminosity function (UV LF) and star formation rate density of distant galaxies ($7.5 < z < 13.5$) in the `blank' fields of the Prime Extragalactic Areas for Reionization Science (PEARLS) survey combined with Early Release Science (ERS) data from the CEERS, GLASS and NGDEEP surveys/fields. We use a combination of SED fitting tools and quality cuts to obtain a reliable selection and characterisation of high-redshift ($z>6.5$) galaxies from a consistently processed set of deep, near-infrared imaging. Within an area of 110 arcmin$^{2}$, we identify 214 candidate galaxies at redshifts $z>6.5$ and we use this sample to study the ultraviolet luminosity function (UV LF) in four redshift bins between $7.5<z<13.5$. The measured number density of galaxies at $z=8$ and $z=9$ match those of past observations undertaken by the em Hubble Space Telescope (HST). However, towards higher redshifts we find that the evolution of the UV LF is mild, resulting in higher measured number densities of UV luminous galaxies at $z=10.5$ and $z=12.5$ compared to predictions from simulations and past HST observations. When examining the star formation rate density of galaxies at this time period, our observations are still consistent with a constant star formation efficiency, are slightly lower than previous early estimations using JWST and support galaxy driven reionization at $z\sim8$.

We make use of classical lattice simulations in 3 + 1 dimensions to study the preheating stage of Higgs Inflation in Einstein-Cartan gravity. Focusing for concreteness on a simplified scenario involving the seminal Nieh-Yan term, we demonstrate the formation of dense and spatially localized oscillon configurations constituting up to 70% of the total energy density. The emergence of these meta-stable objects may lead to a prolonged period of matter domination, effectively modifying the post-inflationary history of the Universe as compared to the metric and Palatini counterparts. Notably, the creation of oscillons comes together with a significant gravitational wave signal, whose typical frequency lies, however, beyond the range accessible by existing and planned gravitational wave experiments. The impact of the Standard Model gauge bosons and fermions and the potential extension of our results to more general Einstein-Cartan settings is also discussed.

The quantum cosmology of the flat Friedmann-Lema{\^i}tre-Robertson-Walker Universe, filled with a scalar field, is considered in the de Broglie-Bohm (dBB) interpretation framework. A stiff-matter quantum bounce solution is obtained. The bouncing and subsequent pre-inflationary and inflationary dynamics are studied in details. We consider some representative primordial inflation models as examples, for which analytical expressions characterizing the dynamical quantities can be explicitly derived. The dependence of the inflationary dynamics on the quantum bounce parameters is then analyzed. The parameters emerging from our description are constrained by requiring the produced dynamics to be in accordance with some key cosmological quantities. The constraining conditions are also illustrated through regions of parameter space in terms of the bounce quantities.

Black-hole binary spin precession admits equilibrium solutions corresponding to systems with (anti-) aligned spins. Among these, binaries in the up-down configuration, where the spin of the heavier (lighter) black hole is co- (counter-) aligned with the orbital angular momentum, might be unstable to small perturbations of the spin directions. The occurrence of the up-down instability leads to gravitational-wave sources that formed with aligned spins but are detected with precessing spins. We present a Bayesian procedure based on the Savage-Dickey density ratio to test the up-down origin of gravitational-wave events. This is applied to both simulated signals, which indicate that achieving strong evidence is within the reach of current experiments, and the LIGO/Virgo events released to date, which indicate that current data are not informative enough.

If two particles move towards a black hole and collide in the vicinity of the horizon, under certain conditions their energy $E_{c.m.}$ in the center of mass frame can grow unbounded. This is the Ba\~{n}ados-Silk-West (BSW) effect. Usually, this effect is considered for extremal horizons and geodesic (or electrogedesic) trajectories. We study this effect in a more general context, when both geometric and dynamic factors are taken into account. We consider generic axially symmetric rotating black holes. The near-horizon behavior of metric coefficients is determined by three numbers $p,~q,$ $k$ that appear in the Taylor expansions for different types of a horizon$.$ This includes nonextremal, extremal and ultraextremal horizons. We also give general classification of possible trajectories that include so-called usual, subcritical, critical and ultracritical ones depending on the near-horizon behavior of the radial component of the four-velocity. We assume that particles move not freely but under the action of some unspecified force. We find when the finiteness of a force and the BSW effect are compatible with each other. The BSW effect implies that one of two particles has fine-tuned parameters. We show that such a particle always requires an infinite proper time for reaching the horizon. Otherwise, either a force becomes infinite or a horizon fails to be regular. This realizes the so-called principle of kinematic censorship that forbids literally infinite $E_{c.m.}$ in any act of collision. The obtained general results are illustrated for the Kerr-Newman-(anti-)de Sitter metric used as an example. The description of diversity of trajectories suggested in our work can be of use also in other contexts, beyond the BSW effect. In particular, we find the relation between a force and the type of a trajectory.

We demonstrate a radio-frequency parametric amplifier that exploits the gate-tunable quantum capacitance of an ultra high mobility two dimensional electron gas (2DEG) in a GaAs heterostructure at cryogenic temperatures. The prototype narrowband amplifier exhibits a gain greater than 20 dB up to an input power of - 66 dBm (1 dB compression), and a noise temperature TN of 1.3 K at 370 MHz. In contrast to superconducting amplifiers, the quantum capacitance parametric amplifier (QCPA) is operable at tesla-scale magnetic fields and temperatures ranging from milli kelvin to a few kelvin. These attributes, together with its low power (microwatt) operation when compared to conventional transistor amplifiers, suggest the QCPA may find utility in enabling on-chip integrated readout circuits for semiconductor qubits or in the context of space transceivers and radio astronomy instruments.

In direct dark matter (DM) detection via scattering off the electrons, the momentum transfer plays a crucial role. Previous work showed that for self-interacting DM, if the DM particle has a size (the so-called puffy DM), the radius effect could dominate the momentum transfer and become another source of velocity dependence for self-scattering cross section. In this work we investigate the direct detection of puffy DM particles with different radii through electron recoil. We find that comparing with the available experimental exclusion limits dominated by the mediator effect for XENON10, XENON100 and XENON1T, the constraints on the puffy DM-electron scattering cross-section become much weaker for large radius DM particles. For small-radius DM particles, the constraints remain similar to the point-like DM case.

We demonstrate how to constrain the degree of absolute alignment of the total angular momenta of LIGO-Virgo binary black holes, looking for a special direction in space that would break isotropy. We also allow for inhomogeneities in the distribution of black holes over the sky. Making use of dipolar models for the spatial distribution and orientation of the sources, we analyze 57 signals with false-alarm rates < 1/yr from the third LIGO-Virgo observing run. Accounting for selection biases, we find the population of LIGO-Virgo black holes to be fully consistent with both homogeneity and isotropy. We additionally find the data to constrain some directions of alignment more than others, and produce posteriors for the directions of total angular momentum of all binaries in our set. All code and data are made publicly available in https://github.com/maxisi/gwisotropy/.

Shibata and Sasaki (1999) introduced the so-called compaction function. Since then, it has been empirically established that the maximum value of this function (or its volume-averaged counterpart) in the long-wavelength solutions gives a very robust threshold of primordial black hole formation. In this paper, we show that in spite of initial intention, the Shibata-Sasaki compaction function cannot be interpreted as the ratio of the mass excess to the areal radius in the constant-mean-curvature slice of their choice but coincides with that in the {\it comoving} slice up to a constant factor depending on the equation of state. We also discuss the gauge-(in)dependence of the legitimate compaction function, i.e., the ratio of the mass excess to the areal radius, in the long-wavelength solutions.

Inflation correlators with one-loop massive exchange encode rich information about the dynamics of the massive loop particles. Their nonanalytic behavior in certain soft limits leads to characteristic oscillatory pattern, which is the leading signal of many particle models of cosmological collider physics. In this work, we investigate systematically such nonanalyticity for arbitrary one-particle-irreducible (1PI) one-loop correlators in various soft limits. With the partial Mellin-Barnes representation, we present and prove a factorization theorem and a cutting rule for arbitrary 1PI one-loop inflation correlators, which is reminiscent of the on-shell cutting rule for flat-space scattering amplitudes. We also show how to understand this factorization theorem from the viewpoint of operator product expansion on the future boundary. As an application of the one-loop factorization theorem, we derive new analytic and exact formulae for nonlocal cosmological collider signals for massive one-loop four-point inflation correlators of all possible 1PI topologies, including the bubble, the triangle, and the box graphs. Finally, we show how to push the computation of nonlocal signals to higher orders in the momentum ratio.

The study of charged particle dynamics in the combined gravitational and magnetic field can provide important theoretical insight into astrophysical processes around black holes. In this paper, we explore the charged particle dynamics in parabolic magnetic field configuration around Schwarzschild black hole, since the paraboloidal shapes of magnetic field lines around black holes are well motivated by the numerical simulations and supported by observations of relativistic jets. Analysing the stability of bounded orbits and using the effective potential approach, we show the possibility of existence of stable circular off-equatorial orbits around the symmetry axis. We also show the influence of radiation reaction force on the dynamics of charged particles, in particular on the chaoticity of the motion and Poincar\'{e} sections, oscillatory frequencies, and emitted electromagnetic spectrum. Applied to Keplerian accretion disks, we show that in parabolic magnetic field configuration, the thin accretion configurations can be either destroyed or transformed into a thick toroidal structure given the radiation reaction and electromagnetic-disk interactions included. Calculating the Fourier spectra for radiating charged particle trajectories, we find that the radiation reaction force does not affect the main frequency peaks, however, it lowers the higher harmonics making the spectrum more flat and diluted in high frequency range.