Papers for Thursday, Sep 30 2021

Interpolating unstructured data using barycentric coordinates becomes infeasible at high dimensions due to the prohibitive memory requirements of building a Delaunay triangulation. We present a new algorithm to construct ad-hoc simplices that are empirically guaranteed to contain the target coordinates, based on a nearest neighbor heuristic and an iterative dimensionality reduction through projection. We use these simplices to interpolate the astrophysical cooling function $\Lambda$ and show that this new approach clearly outperforms our previous implementation at high dimensions.

Massive stellar origin black hole binaries (SBHBs), originating from stars above the pair-instability mass gap, are primary candidate for multi-band gravitational wave (GW) observations. Here we study the possibility to use them as effective dark standard sirens to constrain cosmological parameters. The long lasting inspiral signal emitted by these systems is accessible by the future $Laser \; Interferometer \; Space \; Antenna$ (LISA), while the late inspiral and merger are eventually detected by third generation ground-based telescopes such as the $Einstein \; Telescope$ (ET). The direct measurement of the luminosity distance and the sky position to the source, together with the inhomogeneous redshift distribution of possible host galaxies, allow to infer cosmological parameters by probabilistic means. The efficiency of this statistical method relies in high parameter estimation performances, and we show that this multi-band approach allows a precise determination of the Hubble constant H$_0$ with just ${\cal O}(10)$ detected sources. For selected SBHB population models, assuming $4$ ($10$) years of LISA observations, we find that H$_0$ is tipically determined at $\sim 2\%$ ($\sim 1.5\%$), whereas $\Omega_m$ is only mildly constrained with a typical precision of $30\%$ ($20\%$). The inference procedure sometimes leads to inconsistent results due to the appearance of spurious peaks in the posterior of the parameters. We discuss the origin of this issue and ways to mitigate it.

We measure the near-UV (rest-frame $\sim 2400$\AA) to optical color for early-type galaxies in 12 clusters at $0.3 < z < 1.0$. We show that this is a suitable proxy for the more common far-ultraviolet bandpass used to measure the ultraviolet upturn and find that the upturn is detected to $z=0.6$ in these data, in agreement with previous work. We find evidence that the strength of the upturn starts to wane beyond this redshift and largely disappears at $z=1$. Our data is most consistent with models where early-type galaxies contain minority stellar populations with non-cosmological helium abundances, up to around 46\%, formed at $z \geq 3$, resembling multiple stellar population globular clusters in our Galaxy. This suggests that elliptical galaxies and globular clusters share similar chemical evolution and star formation histories. The vast majority of the stellar mass in these galaxies also must have been in place at $z > 3$.

The cosmological constraints from the kinetic Sunyaev-Zeldovich experiments are degenerate with the optical depth measurement, which is commonly known as the optical-depth degeneracy. In this work, we introduce a new statistic based on the first moment of relative velocity between pairs in a triplet, which is capable of constraining cosmological parameters independent of the optical depth, and $\sigma_8$. Using 22,000 $N$-body simulations from the Quijote suite, we quantify the information content in the new statistic using Fisher matrix forecast. We find that it is able to obtain strong constraints on the cosmological parameters, particularly on the summed neutrino mass. The constraints have a factor of 6.2-12.9, and 2.3-5.7 improvement on all cosmological model parameters when compared to those obtained from the mean pairwise velocity, and the redshift-space halo power spectrum, respectively. Thus the new statistic paves a way forward to constrain cosmological parameters independent of the optical depth and $\sigma_8$ using data from future kinetic Sunyaev-Zeldovich experiments alone.

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $\lambda_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects.

Tidal Disruption Events (TDEs) have been found to show a preference for post-starburst (PS) and quiescent Balmer-strong (QBS) galaxies. This preference can be used to help find TDEs in transient surveys. But what other transients might "contaminate" such a search, and by how much? We examine all reported transients coincident with the centers of galaxies in the French & Zabludoff (2018) catalog of spectroscopically-confirmed PS and QBS galaxies and photometrically-identified PS and QBS galaxy candidates. We find that TDEs and Type Ia supernovae (SNe) are the only types of transients classified in the centers of these galaxies (aside from one AGN flare), with Type Ia SNe being 8.3+-0.2 times more prevalent than TDEs (1-sigma confidence bounds). This factor is ~2.7 times lower than in a control sample of quiescent galaxies. Narrowing the sample to spectroscopically-confirmed QBS galaxies doesn't change these statistics much. In spectroscopically-confirmed PS galaxies, on the other hand, TDEs are the ones that outnumber Type Ia SNe 2+-0.6 to 1. However, there are few such galaxies in the catalog. By classifying transients from the entire catalog, three times more TDEs are expected to be found, but with a ~16-times larger Type Ia SN contamination. We use the public ZTF photometric archive to search for possibly missed TDEs in the French & Zabludoff (2018) galaxies. We find three unclassified clear transients - neither of which are likely missed TDEs based on their light-curve colors.

We present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at field level using the mass-weighted halo field in the AbacusSummit simulations at $z=0.5$. Unlike the bias expansion method widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity $\delta_1$ in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of $\delta_1$ and other quantities, and shows a plateau at high $\delta_1$. Especially for massive halos with $6\times10^{12}\ h^{-1}\ M_\odot$, the halo-to-mass ratio starts soaring up at $\delta_1>0$, substantially different from the predictions of the bias expansion. We show that for $M>3\times10^{11}\ h^{-1}\ M_\odot$ halos, a non-parametric halo-to-mass ratio as a function of $\delta_1$ and its local derivative $\nabla^2\delta_1$ can recover the halo power spectra to sub-percent level accuracy for wavenumbers $k=0.01-0.1\ h\ {\rm Mpc}^{-1}$ given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At $k<0.01\ h\ {\rm Mpc}^{-1}$ and for $M>6\times10^{12}\ h^{-1}\ M_\odot$ halos, our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.

Metal-rich near-Earth asteroids (NEAs) represent a small fraction of the NEA population that is mostly dominated by S- and C-type asteroids. Because of this, their identification and study provide us with a unique opportunity to learn more about the formation and evolution of this particular type of bodies, as well as their relationship with meteorites found on Earth. We present near-infrared (NIR) spectroscopic data of NEAs 6178 (1986 DA) and 2016 ED85. We found that the spectral characteristics of these objects are consistent with those of metal-rich asteroids, showing red slopes, convex shapes, and a weak pyroxene absorption band at $\sim$0.93 $\mu$m. The compositional analysis showed that they have a pyroxene chemistry of Fs$_{40.6\pm3.3}$Wo$_{8.9\pm1.1}$ and a mineral abundance of $\sim$15% pyroxene and 85% metal. We determined that these objects were likely transported to the near-Earth space via the 5:2 mean motion resonance with Jupiter. Asteroid spectra were compared with the spectra of mesosiderites and bencubbinites. Differences in the NIR spectra and pyroxene chemistry suggest that bencubbinites are not good meteorite analogs. Mesosiderites were found to have a similar pyroxene chemistry and produced a good spectral match when metal was added to the silicate component. We estimated that the amounts of Fe, Ni, Co, and the platinum group metals present in 1986 DA could exceed the reserves worldwide.

Lithium is one of the few elements produced during the Big Bang Nucleosynthesis in the early universe. Moreover, its fragility makes it useful as a proxy for stellar environmental conditions. As such, the lithium abundance in old systems is at the core of different astrophysical problems. Stars in the lower red giant branch allow studying globular clusters where main sequence stars are too faint to be observed. We use these stars to analyze the initial Li content of the clusters and compare it to cosmological predictions, to measure spreads in Li between different stellar populations, and to study signs of extra depletion in these giants. We use GIRAFFE spectra to measure the lithium and sodium abundances of lower red giant branch stars in 5 globular clusters. These cover an extensive range in metallicity, from [Fe/H]$\sim-0.7$ to [Fe/H]$\sim-2.3$ dex. We find that the lithium abundance in these lower red giant branch stars forms a plateau, with values from $\mathrm{A(Li)_{NLTE}}=0.84$ to $1.03$ dex, showing no clear correlation with metallicity. When using stellar evolutionary models to calculate the primordial abundance of these clusters, we recover values $\mathrm{A(Li)_{NLTE}}=2.1-2.3$ dex, consistent with the constant value observed in warm metal-poor halo stars, the Spite plateau. Additionally, we find no difference in the lithium abundance of first and second population stars in each cluster. We also report the discovery of a Li-rich giant in the cluster NGC3201, with $\mathrm{A(Li)_{NLTE}}=1.63\pm0.18$ dex, where the enrichment mechanism is probably pollution from external sources.

The ionising feedback of young massive stars is well known to influence the dynamics of the birth environment and hence plays an important role in regulating the star formation process in molecular clouds. For this reason, modern hydrodynamics codes adopt a variety of techniques accounting for these radiative effects. A key problem hampering these efforts is that the hydrodynamics are often solved using smoothed particle hydrodynamics (SPH), whereas radiative transfer is typically solved on a grid. Here we present a radiation-hydrodynamics (RHD) scheme combining the SPH code Phantom and the Monte Carlo Radiative Transfer (MCRT) code CMacIonize, using the particle distribution to construct a Voronoi grid on which the MCRT is performed. We demonstrate that the scheme successfully reproduces the well-studied problem of D-type H II region expansion in a uniform density medium. Furthermore, we use this simulation setup to study the robustness of the RHD code with varying choice of grid structure, density mapping method, and mass and temporal resolution. To test the scheme under more realistic conditions, we apply it to a simulated star-forming cloud reminiscing those in the Central Molecular Zone of our galaxy, in order to estimate the amount of ionised material that a single source could create. We find that a stellar population of several $10^3~\rm{M_{\odot}}$ is needed to noticeably ionise the cloud. Based on our results, we formulate a set of recommendations to guide the numerical setup of future and more complex simulations of star forming clouds.

Spectroscopic surveys are providing large samples of lithium (Li) measurements in evolved stars. A seemingly unexpected result from this work has been the apparent detection of Li at a higher rate in core helium-burning stars than in luminous shell hydrogen-burning stars, which has been interpreted as evidence for ubiquitous Li production on the upper red giant branch or at helium ignition. This is distinct from the "Li-rich giant" problem and reflects bulk red clump star properties. We provide an analysis of the GALAH Li data that accounts for the distribution of progenitor masses of field red clump stars observed today. Using standard models of the post-main sequence evolution of low-mass stars, we show that the observed distribution of Li among the bulk of field clump giants is natural, and that observations of clump and red giants in the field should not be compared without correcting for population effects. For typical stellar population distributions, moderate Li abundances among the bulk of field clump giants are expected without the need for a new Li production mechanism. Our model predicts a large fraction of very low Li abundances from low mass progenitors, with higher Li abundances from higher mass ones. Moreover, there should be a large number of upper limits for red clump stars, and higher abundances should correspond to higher masses. The most recent GALAH data indeed confirm the presence of large numbers of upper limits, and a much lower mean Li abundance in clump stars, which is concordant with our interpretation.

High contrast imaging (HCI) systems rely on active wavefront control (WFC) to deliver deep raw contrast in the focal plane, and on calibration techniques to further enhance contrast by identifying planet light within the residual speckle halo. Both functions can be combined in an HCI system and we discuss a path toward designing HCI systems capable of calibrating residual starlight at the fundamental contrast limit imposed by photon noise. We highlight the value of deploying multiple high-efficiency wavefront sensors (WFSs) covering a wide spectral range and spanning multiple optical locations. We show how their combined information can be leveraged to simultaneously improve WFS sensitivity and residual starlight calibration, ideally making it impossible for an image plane speckle to hide from WFS telemetry. We demonstrate residual starlight calibration in the laboratory and on-sky, using both a coronagraphic setup, and a nulling spectro-interferometer. In both case, we show that bright starlight can calibrate residual starlight.

Stock 2 is a little-studied open cluster that shows an extended main-sequence turnoff (eMSTO). In order to investigate this phenomenon and characterise the cluster itself we performed high-resolution spectroscopy in the framework of the Stellar Population Astrophysics (SPA) project. We employed the High Accuracy Radial velocity Planet Searcher in North hemisphere spectrograph (HARPS-N) at the Telescopio Nazionale Galileo (TNG). We completed our observations with additional spectra taken with the Catania Astrophysical Observatory Spectrograph (CAOS). In total we observed 46 stars (dwarfs and giants), which represent, by far, the largest sample collected for this cluster to date. We provide the stellar parameters, extinction, radial and projected rotational velocities for most of the stars. Chemical abundances for 21 species with atomic numbers up to 56 have also been derived. We notice a differential reddening in the cluster field whose average value is 0.27 mag. It seems to be the main responsible for the observed eMSTO, since it cannot be explained as the result of different rotational velocities, as found in other clusters. We estimate an age for Stock 2 of 450$\pm$150 Ma which corresponds to a MSTO stellar mass of $\approx$2.8 M$_{\odot}$. The cluster mean radial velocity is around 8.0 km s$^{-1}$. We find a solar-like metallicity for the cluster, [Fe/H]=$-$0.07$\pm$0.06, compatible with its Galactocentric distance. MS stars and giants show chemical abundances compatible within the errors, with the exceptions of Barium and Strontium, which are clearly overabundant in giants, and Cobalt, which is only marginally overabundant. Finally, Stock 2 presents a chemical composition fully compatible with that observed in other open clusters of the Galactic thin disc.

The standard model of stellar evolution in Close Binary Systems assumes that during mass transfer episodes the system is in a synchronised and circularised state. Remarkably, the redback system PSR J1723-2837 has an orbital period derivative $\dot{P}_{orb}$ too large to be explained by this model. Motivated by this fact, we investigate the action of tidal forces in between two consecutive mass transfer episodes for a system under irradiation feedback, which is a plausible progenitor for PSR J1723-2837. We base our analysis on Hut's treatment of equilibrium tidal evolution, generalised by considering the donor as a two layers object that may not rotate as a rigid body. We also analyse three different relations for the viscosity with the tidal forcing frequency. We found that the large value measured for $\dot{P}_{orb}$ can be reached by systems where the donor star rotates slower (by few percent) than the orbit just after mass transfer episodes. Van Staden & Antoniadis have observed this object and reported a lack of synchronism, opposite to that required by the Hut's theory to account for the observed $\dot{P}_{orb}$. Motivated by this discrepancy, we analyse photometric data obtained by the spacecraft Kepler second mission K2, with the purpose of identifying the periods present in PSR J1723-2837. We notice several periods close to those of the orbit and the rotation. The obtained periods pattern reveals the presence of a more complex phenomenology, which would not be well described in the frame of the weak friction model of equilibrium tides.

The phasing parameter F determines the relative phasing between satellites in different orbital planes and thereby affects the relative position of the satellites in a constellation. The collisions between satellites within the constellation can be avoided if the minimum distance among them is large. From among the possible values of F in a constellation, a value of F is desired that leads to the maximum value of the minimum distance between satellites. We investigate F for two biggest upcoming satellite constellations including Starlink Phase 1 Version 3 and Kuiper Shell 2. No existing work or FCC filing provides a value of F that is suitable for these two constellations. We look for the best value of F in these constellations that provides the maximum value of the minimum distance to ensure intra-constellation avoidance of collisions between satellites. To this end, we simulate each constellation for each value of F to find its best value based on ranking. Out of the 22 and 36 possible values of F for Starlink Phase 1 Version 3 and Kuiper Shell 2, respectively, it is observed that the best value of F with highest ranking is 17 and 11 that leads to the largest minimum distance between satellites of 61.83 km and 55.89 km in these constellations, respectively.

We present updated radial-velocity (RV) analyses of the AU Mic system. AU Mic is a young (22 Myr) early M dwarf known to host two transiting planets - $P_{b}\sim8.46$ days, $R_{b}=4.38_{-0.18}^{+0.18}\ R_{\oplus}$, $P_{c}\sim18.86$ days, $R_{c}=3.51_{-0.16}^{+0.16}\ R_{\oplus}$. With visible RVs from CARMENES-VIS, CHIRON, HARPS, HIRES, {\sc {\textsc{Minerva}}}-Australis, and TRES, as well as near-infrared (NIR) RVs from CARMENES-NIR, CSHELL, IRD, iSHELL, NIRSPEC, and SPIRou, we provide a $5\sigma$ upper limit to the mass of AU Mic c of $M_{c}\leq20.13\ M_{\oplus}$ and present a refined mass of AU Mic b of $M_{b}=20.12_{-1.57}^{+1.72}\ M_{\oplus}$. Used in our analyses is a new RV modeling toolkit to exploit the wavelength dependence of stellar activity present in our RVs via wavelength-dependent Gaussian processes. By obtaining near-simultaneous visible and near-infrared RVs, we also compute the temporal evolution of RV-``color'' and introduce a regressional method to aid in isolating Keplerian from stellar activity signals when modeling RVs in future works. Using a multi-wavelength Gaussian process model, we demonstrate the ability to recover injected planets at $5\sigma$ significance with semi-amplitudes down to $\approx$ 10\,m\,s$^{-1}$ with a known ephemeris, more than an order of magnitude below the stellar activity amplitude. However, we find that the accuracy of the recovered semi-amplitudes is $\sim$50\% for such signals with our model.

We present a novel method to produce empirical generative models of all kinds of astronomical transients from datasets of unlabeled light curves. Our hybrid model, that we call ParSNIP, uses a neural network to model the unknown intrinsic diversity of different transients and an explicit physics-based model of how light from the transient propagates through the universe and is observed. The ParSNIP model predicts the time-varying spectra of transients despite only being trained on photometric observations. With a three-dimensional intrinsic model, we are able to fit out-of-sample multiband light curves of many different kinds of transients with model uncertainties of 0.04-0.06 mag. The representation learned by the ParSNIP model is invariant to redshift, so it can be used to perform photometric classification of transients even with heavily biased training sets. Our classification techniques significantly outperform state-of-the-art methods on both simulated (PLAsTiCC) and real (PS1) datasets with 2.3$\times$ and 2$\times$ less contamination respectively for classification of Type~Ia supernovae. We demonstrate how our model can identify previously-unobserved kinds of transients and produce a sample that is 90% pure. The ParSNIP model can also estimate distances to Type Ia supernovae in the PS1 dataset with an RMS of 0.150 $\pm$ 0.007 mag compared to 0.155 $\pm$ 0.008 mag for the SALT2 model on the same sample. We discuss how our model could be used to produce distance estimates for supernova cosmology without the need for explicit classification.

We present results from timing and spectral analysis of the HMXB X-ray pulsar IGR J19294+1816 observed using \asr during its recent Type-I outburst in October, 2019. AstroSat observations sampled the outburst at the decline phase right after the outburst peak. We carried out timing analysis on the light curves obtained using the Large Area X-ray Proportional Counter (LAXPC) instrument on board AstroSat and measured a spin period of 12.485065$\pm$0.000015 s. The peak in the power density spectrum (PDS) corresponding to the spin period of 12.48 s also shows a broadened base. We also detected a Quasi Periodic Oscillation (QPO) feature at 0.032$\pm$0.002 Hz with an RMS fractional amplitude of $\sim$18% in the PDS. We further carried out a joint spectral analysis using both the Soft X-ray Telescope (SXT) and the LAXPC instruments and detected a Cyclotron Resonant Scattering Feature (CRSF) at 42.7$\pm$0.9 keV and an Fe emission line at 6.4$\pm$0.1 keV. IGR J19294+1816, being an intermediate spin pulsar, has exhibited a plethora of spectral and timing features during its most recent 2019 outburst, adding it to the list of transients that exhibit both a QPO as well as a CRSF.

We present the radio properties of optically selected quasars with $z\geq3$. The complete sample consists of 102 quasars with a flux density level $S_{1.4}\geq100$ mJy in a declination range -35$^{\circ}$ $\leq$ Dec $\leq$ +49$^{\circ}$. The observations were obtained in 2017-2020 using the radio telescope RATAN-600. We measured flux densities at six frequencies 1.2, 2.3, 4.7, 8.2, 11.2, and 22 GHz quasi-simultaneously with uncertainties of 9-31 %. The detection rate is 100, 89, and 46 % at 4.7, 11.2, and 22 GHz, respectively. We have analysed the averaged radio spectra of the quasars based on the RATAN and literature data. We classify 46 % of radio spectra as peaked-spectrum, 24 % as flat, and none as ultra-steep spectra ($\alpha\leq-1.1$). The multifrequency data reveal that a peaked spectral shape (PS) is a common feature for bright high-redshift quasars. This indicates the dominance of bright compact core emission and the insignificant contribution of extended optically thin kpc-scale components in observed radio spectra. Using these new radio data, the radio loudness $\log~R$ was estimated for 71 objects with a median value of 3.5, showing that the majority of the quasars are highly radio-loud with $\log~R>2.5$. We have not found any significant correlation between $z$ and $\alpha$. Several new megahertz-peaked spectrum (MPS) and gigahertz-peaked spectrum (GPS) candidates are suggested. Further studies of their variability and additional low-frequency observations are needed to classify them precisely.

The geometric characteristics of dust clouds provide important information on the physical processes that structure such clouds. One of such characteristics is the $2D$ fractal dimension $D$ of a cloud projected onto the sky plane. In previous studies, which were mostly based on infrared (IR) data, the fractal dimension of individual clouds was found to be in a range from 1.1 to 1.7 with a preferred value of 1.2--1.4. In the present work, we use data from Stripe82 of the Sloan Digital Sky Survey to measure the fractal dimension of the cirrus clouds. This is done here for the first time for optical data with significantly better resolution as compared to IR data. To determine the fractal dimension, the perimeter-area method is employed. We also consider IR (IRAS and Herschel) counterparts of the corresponding optical fields to compare the results between the optical and IR. We find that the averaged fractal dimension across all clouds in the optical is $\langle D \rangle =1.69^{+0.05}_{-0.05}$ which is significantly larger then the fractal dimension of its IR counterparts $\langle D\rangle=1.38^{+0.07}_{-0.06}$. We examine several reasons for this discrepancy (choice of masking and minimal contour level, image and angular resolution, etc.) and find that for approximately half of our fields the different angular resolution (point spread function) of the optical and IR data can explain the difference between the corresponding fractal dimensions. For the other half of the fields, the fractal dimensions of the IR and visual data remain inconsistent, which can be associated with physical properties of the clouds, but further physical simulations are required to prove it.

Faraday Rotation Measures (RM) should be interpreted with caution because there could be multiple magneto-ionized medium components that contribute to the net Faraday rotation along sight-lines. We introduce a simple test using Galactic diffuse polarised emission that evaluates whether structures evident in RM observations are associated with distant circumgalactic medium (CGM) or foreground interstellar medium (ISM). We focus on the Magellanic Leading Arm region where a clear excess of RM was previously reported. There are two gaseous objects standing out in this direction: the distant Magellanic Leading Arm and the nearby Antlia supernova remnant (SNR). We recognized narrow depolarised filaments in the $2.3\,\rm GHz$ S-band Polarization All Sky Survey (S-PASS) image that overlaps with the reported RM excess. We suggest that there is a steep gradient in Faraday rotation in a foreground screen arising from the Antlia SNR. The estimated strength of the line-of-sight component of the magnetic field is $B_{\parallel}\sim 5\,\rm\mu G$, assuming that the excess of RM is entirely an outcome of the magnetized supernova shell. Our analysis indicates that the overlap between the RM excess and the Magellanic Leading Arm is only a remarkable coincidence. We suggest for future RM grid studies that checking Galactic diffuse polarisation maps is a convenient way to identify local Faraday screens.

Kinematic weak lensing describes the distortion of a galaxy's projected velocity field due to lensing shear, an effect recently reported for the first time by Gurri et al. based on a sample of 18 galaxies at $z \sim 0.1$. In this paper, we develop a new formalism that combines the shape information from imaging surveys with the kinematic information from resolved spectroscopy to better constrain the lensing distortion of source galaxies and to potentially address systematic errors that affect conventional weak-lensing analyses. Using a Bayesian forward model applied to mock galaxy observations, we model distortions in the source galaxy's velocity field simultaneously with the apparent shear-induced offset between the kinematic and photometric major axes. We show that this combination dramatically reduces the statistical uncertainty on the inferred shear, yielding statistical error gains of a factor of 2--6 compared to kinematics alone. While we have not accounted for errors from intrinsic kinematic irregularities, our approach opens kinematic lensing studies to higher redshifts where resolved spectroscopy is more challenging. For example, we show that ground-based integral-field spectroscopy of background galaxies at $z \sim 0.7$ can deliver gravitational shear measurements with S/N $\sim 1$ per source galaxy at 1 arcminute separations from a galaxy cluster at $z \sim 0.3$. This suggests that even modest samples observed with existing instruments could deliver improved galaxy cluster mass measurements and well-sampled probes of their halo mass profiles to large radii.

The Widefield Arecibo Virgo Extragalactic Survey (WAVES) was an ongoing HI survey of the Virgo Cluster with the Arecibo Observatory's 305m William E. Gordon Telescope at the time of its structural failure. The full 20 square degrees of the southern field and 10 of the planned 35 square degrees of the northern field had been observed to full depth, adding to 25 square degrees observed to the same depth in the cluster by the Arecibo Galaxy Environment Survey. We here review what WAVES reveals about four optically dark HI structures that were previously discovered in the survey area, including two that are not seen despite being well above our detection limit.

During the survey phase of the Kepler mission, several thousands of stars were observed in short cadence, allowing the detection of solar-like oscillations in more than 500 main-sequence and sub-giant stars. Later, the Kepler Science Office discovered an issue in the calibration that affected half of the short-cadence data, leading to a new data release (DR25) with improved corrections. We re-analyze the one-month time series of the Kepler survey phase to search for new solar-like oscillations. We study the seismic parameters of 99 stars (46 targets with new reported solar-like oscillations) increasing by around 8% the known sample of solar-like stars with asteroseismic analysis of the short-cadence data from Kepler. We compute the masses and radii using seismic scaling relations and find that this new sample populates the massive stars (above 1.2Ms and up to 2Ms) and subgiant phase. We determine the granulation parameters and amplitude of the modes, which agree with previously derived scaling relations. The stars studied here are slightly fainter than the previously known sample of main-sequence and subgiants with asteroseismic detections. We also study the surface rotation and magnetic activity levels of those stars. Our sample of has similar levels of activity compared to the previously known sample and in the same range as the Sun between the minimum and maximum of its activity cycle. We find that for 7 stars, a possible blend could be the reason for the previous non detection. We compare the radii obtained from the scaling relations with the Gaia ones and find that the Gaia radii are overestimated by 4.4% on average compared to the seismic radii and a decreasing trend with evolutionary stage. We re-analyze the DR25 of the main-sequence and sub-giant stars with solar-like oscillations previously detected and provide their global seismic parameters for a total of 526 stars.

We present the results of the COS Intragroup Medium (COS-IGrM) Survey that used the Cosmic Origins Spectrograph on the Hubble Space Telescope to observe a sample of 18 UV bright quasars, each probing the intragroup medium (IGrM) of a galaxy group. We detect Ly$\alpha$, C II, N V, Si II, Si III, and O VI in multiple sightlines. The highest ionization species detected in our data is O VI, which was detected in 8 out of 18 quasar sightlines. The wide range of ionization states observed provide evidence that the IGrM is patchy and multiphase. We find that the O VI detections generally align with radiatively cooling gas between $10^{5.8}$ and $10^6$ K. The lack of O VI detections in 10 of the 18 groups illustrates that O VI may not be the ideal tracer of the volume filling component of the IGrM. Instead, it either exists at trace levels in a hot IGrM or is generated in the boundary between the hotter IGrM and cooler gas.

The two most successful methods for exoplanet detection rely on the detection of planetary signals in photometric and radial velocity time-series. This depends on numerical techniques that exploit the synergy between data and theory to estimate planetary, orbital, and/or stellar parameters. In this work, we present a new version of the exoplanet modelling code pyaneti. This new release has a special emphasis on the modelling of stellar signals in radial velocity time-series. The code has a built-in multi-dimensional Gaussian process approach to modelling radial velocity and activity indicator time-series with different underlying covariance functions. This new version of the code also allows multi-band and single transit modelling; it runs on Python 3, and features overall improvements in performance. We describe the new implementation and provide tests to validate the new routines that have direct application to validate the new routines that have direct application to exoplanet detection and characterisation. We have made the code public and freely available at https://github.com/oscaribv/pyaneti. We also present the codes citlalicue and citlalatonac that allow one to create synthetic photometric and spectroscopic time-series, respectively, with planetary and stellar-like signals.

We perform the first study of the principal core $g$-mode oscillation in hybrid stars containing quark matter, utilizing a crossover model for the hadron-to-quark transition inspired by lattice QCD. The ensuing results are compared with our recent findings of $g$-mode frequencies in hybrid stars with a first-order phase transition using Gibbs constructions. We find that models using Gibbs construction yield $g$-mode amplitudes and the associated gravitational energy radiated that dominate over those of the chosen crossover model owing to the distinct behaviors of the equilibrium and adiabatic sound speeds in the various models. Based on our results, we conclude that were $g$-modes to be detected in upgraded LIGO and Virgo detectors it would indicate a first-order phase transition akin to a Gibbs construction.

We present XMM-Newton X-ray observations of nine confirmed lensed quasars at $1 \lesssim z \lesssim 3$ identified by the Gaia Gravitational Lens program. Eight systems are strongly detected, with 0.3--8.0 keV fluxes $F_{0.3-8.0} \gtrsim 5 \times 10^{-14}\ {\rm erg}\ {\rm cm}^{-2}\ {\rm s}^{-1}$. Modeling the X-ray spectra with an absorbed power law, we derive power law photon indices and 2--10 keV luminosities for the eight detected quasars. In addition to presenting sample properties for larger quasar population studies and for use in planning for future caustic crossing events, we also identify three quasars of interest: a quasar that shows evidence of flux variability from previous ROSAT observations, the most closely-separated individual lensed sources resolved by XMM-Newton, and one of the X-ray brightest quasars known at $z>3$. These sources represent the tip of discovery that will be enabled by SRG/eROSITA.

We explore the possibility that radio loud gamma-ray bursts (GRBs) result from the collapse of a massive star in an interacting binary system, while radio quiet GRBs are produced by the collapse of a single massive star. A binary collapsar system can have the necessary angular momentum and energy budget to explain the longer prompt gamma-ray durations and higher isotropic energies seen in the the radio loud sub-sample of long GRBs. Additionally, tidal interactions between the stars in binary systems can lead to rich and extended circumstellar environments that allow for the presence of the long-lived radio afterglows seen in the radio loud systems. Finally, the relative fraction of stars in binary systems versus single star systems appears consistent with the fraction of radio loud versus radio quiet GRBs.

Recently, two photons with energy of about 1 PeV have been detected by LHAASO from the Crab nebula, opening an ultra-high energy window for studying the pulsar wind nebulae (PWNe). Remarkably, the LHAASO spectrum at the highest-energy end shows a possible hardening, which could indicate the presence of a new component. A two-component scenario with a main electron component and a secondary proton component has been proposed to explain the whole spectrum of the Crab Nebula, requiring a proton energy of $10^{46}-10^{47}{\rm ergs}$ remaining in the present Crab Nebula. In this paper, we study the energy content of relativistic protons in pulsar winds with the LHAASO data of the Crab Nebula, considering the effect of diffusive escape of relativistic protons. Depending on the extent of the escape of relativistic protons, the total energy of protons lost in the pulsar wind could be 10-100 times larger than that remaining in the nebula presently. We find that the current LHAASO data allows up to $(10-50)\%$ of the spindown energy of pulsars being converted into relativistic protons. The escaping protons from PWNe could make a considerable contribution to the cosmic-ray flux of 10-100 PeV. We also discuss the leptonic scenario for the possible spectral hardening at PeV energies.

We measure the star formation rate (SFR) per unit gas mass and the star formation efficiency (SFE$_{\rm gas}$ for total gas, SFE$_{\rm mol}$ for the molecular gas) in 81 nearby galaxies selected from the EDGE-CALIFA survey, using $^{12}$CO(J=1-0) and optical IFU data. For this analysis we stack CO spectra coherently by using the velocities of H$\alpha$ detections to detect fainter CO emission out to galactocentric radii $r_{\rm gal} \sim 1.2 r_{25}$ ($\sim 3 R_{\rm e}$), and include the effects of metallicity and high surface densities in the CO-to-H$_2$ conversion. We determine the scale lengths for the molecular and stellar components, finding a close to 1:1 relation between them. This result indicates that CO emission and star formation activity are closely related. We examine the radial dependence of SFE$_{\rm gas}$ on physical parameters such as galactocentric radius, stellar surface density $\Sigma_{\star}$, dynamical equilibrium pressure $P_{\rm DE}$, orbital timescale $\tau_{\rm orb}$, and the Toomre $Q$ stability parameter (including star and gas $Q_{\rm star+gas}$). We observe a generally smooth, continuous exponential decline in the SFE$_{\rm gas}$ with $r_{\rm gal}$. The SFE$_{\rm gas}$ dependence on most of the physical quantities appears to be well described by a power-law. Our results also show a flattening in the SFE$_{\rm gas}$-$\tau_{\rm orb}$ relation at $\log[\tau_{\rm orb}]\sim 7.9-8.1$ and a morphological dependence of the SFE$_{\rm gas}$ per orbital time, which may reflect star formation quenching due to the presence of a bulge component. We do not find a clear correlation between SFE$_{\rm gas}$ and $Q_{\rm star+gas}$.

The James Webb Space Telescope (JWST) is estimated to observe objects as far as z = 15 with 100 times the sensitivity of the Hubble Space Telescope (HST). Several previous simulations have predicted the characteristics of the JWST deep field images using assumptions of universe geometry, the galaxy number densities, and the evolutions of high-z galaxies. While the assumptions made by these previous simulations are based on initial conditions defined by recent observations, the corresponding ranges of uncertainty in their measurements can lead to substantially different results. This study presents a novel geometric-focused deep field simulation using an ensemble approach to demonstrate the high variability in results due to the uncertainty ranges of the measurements used as initial conditions. A parameter sensitivity ensemble was run perturbing each initial condition individually through its range of uncertainty to determine which initial condition range of uncertainty has the largest effect on the resulting simulation. A 1000-member ensemble was run perturbing the initial conditions through their ranges of uncertainty obtained from recent estimates. A galaxy coverage percentage was calculated for each ensemble member and averaged together. The Apparent Galaxy Wall (AGW) effect is introduced and defined as >= 50% of a deep field image occupied by galaxies. A one-way one-sample t-test was conducted to conclude the JWST is likely to observe the AGW effect with an estimated galaxy coverage percentage of 55.12 +/- 30.30%. This study finds the most sensitive parameter to changes within its range of uncertainty is the estimated number of unseen galaxies in the HUDF. A discussion is included on the potential impacts of the AGW effect being observed and its potential to form a pseudo-cosmological horizon that may inhibit the effectiveness of future observatories.

We investigate the dark matter halos of 256 star-forming disc-like galaxies at $z\sim 1$ using the KMOS redshift one spectroscopic survey (KROSS). This sample covers the redshifts $0.6 \leq z \leq 1.04$, effective radii $0.69 \leq R_e [\mathrm{kpc}] \leq 7.76$, and total stellar masses $8.7 \leq log(M_{\mathrm{star}} \ [\mathrm{M_\odot}]) \leq 11.32$. We present a mass modelling approach to study the rotation curves of these galaxies, which allow us to dynamically calculate the physical properties associated with the baryons and the dark matter halo. For the former we assume a Freeman disc, while for the latter we employ the NFW and the Burkert halo profiles, separately. At the end, we compare the results of both cases with state-of-the-art cosmological galaxy simulations (EAGLE, TNG100 and TNG50). We find that the {\em cored} dark matter halo emerged as the dominant quantity from a radius 1-3 times the effective radius. Its fraction to the total mass is in good agreement with the outcome of hydrodynamical galaxy simulations. Remarkably, we found that the dark matter core of $z\sim 1$ star-forming galaxies are smaller and denser than their local counterparts. We conclude that dark matter halos have gradually expanded over the past 6.5 Gyrs. That is, observations are capable of capturing the dark matter response to the baryonic processes (e.g. feedbacks), and thus giving us the first empirical evidence of {\em gravitational potential fluctuations} in the inner region of galaxies, which can be verified with deep surveys and future missions.

We analyze the $^{134}_{55}$Cs $\to {}^{134}_{56}$Ba and $^{135}_{55}$Cs$\to {}^{135}_{56}$Ba $\beta^-$ decays, which are crucial production channels for Ba isotopes in Asymptotic Giant Branch (AGB) stars. We reckon, from relativistic quantum mechanis, the effects of multichannel scattering onto weak decays, including nuclear and electronic excited states populated above 10 keV, for both parent and daughter nuclei. We find a significant increase (by more than a factor 3 for 134Cs) of the half-lives with respect to previous recommendations [1, 2], and we discuss our method in view of these last calculations based on general systematics. The major impact on half-lives comes from nuclear excited state decays, while including electronic temperatures yields a 20% increase, at energies typical of low- and intermediate-mass AGB stars (M<=8Msun). Our predictions strongly modify branching ratios along the s-process path, and allow nucleosynthesis models to account well for the isotopic admixtures of Ba in presolar SiC grains.

We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the survey speed of the telescope. With this upgraded instrument, an imaging survey covering an area of 2300 deg2 is being performed which will deliver both continuum and spectral line data sets, of which the first data has been publicly released. In addition, a time domain transient and pulsar survey covering 15,000 deg2 is in progress. An overview of the Apertif science drivers, hardware and software of the upgraded telescope is presented, along with its key performance characteristics.

We suggest a new method to reconstruct, within canonical single-field inflation, the inflaton potential directly from the primordial power spectrum which may deviate significantly from near scale-invariance. Our approach relies on a more generalized slow-roll approximation than the standard one, and can probe the properties of the inflaton potential reliably. We give a few examples for reconstructing potential and discuss the validity of our method.

Magnetar Giant Flares (MGFs) have been long proposed to contribute at least a sub-sample of the observed short Gamma-ray Bursts (GRBs). The recent discovery of the short GRB 200415A in the nearby galaxy NGC 253 established a textbook-version connection between these two phenomena. Unlike previous observations of the Galactic MGFs, the unsaturated instrument spectra of GRB 200415A provide for the first time an opportunity to test the theoretical models with the observed $\gamma$-ray photons. This paper proposed a new readily fit-able model for the MGFs, which invokes an expanding fireball Comptonized by the relativistic magnetar wind at photosphere radius. In this model, a large amount of energy is released from the magnetar crust due to the magnetic reconnection or the starquakes of the star surface and is injected into confined field lines, forming a trapped fireball bubble. After breaking through the shackles and expanding to the photospheric radius, the thermal photons of the fireball are eventually Comptonized by the relativistic $e^{\pm}$ pairs in the magnetar wind region, which produces additional higher-energy gamma-ray emission. The model predicts a modified thermal-like spectrum characterized by a low-energy component in the Rayleigh-Jeans regime, a smooth component affected by coherent Compton scattering in the intermediate energy range, and a high-energy tail due to the inverse Compton process. By performing a Monte-Carlo fit to the observational spectra of GRB 200415A, we found that the observation of the burst is entirely consistent with our model predictions.

We study the radiatively driven fluid jets around a non-rotating black hole. The radiation arising from the inner compact corona and outer sub-Keplerian part of the disc accelerates the jets. We obtain the steady state, semi-analytical, radiatively driven outflow solutions. The thermodynamics of the outflow is described by a variable adiabatic index equation of state. We develop a TVD routine to investigate the time dependent behaviour of the radiatively driven bipolar outflow. We inject with flow variables from the steady state outflow solutions in the TVD code and allow the code to settle into steady state and match the numerical results with the steady state solution. The radiation arising out of the accretion disc can provide a wide range of jet solutions, depending upon parameters like the intensity of disc, location of the inner corona etc. We induce the time dependence of the radiation field by inducing oscillation of the inner corona of the accretion disc. The radiation field then makes the bipolar outflow time dependent. We show that a non-steady radiation field arising out of disc oscillations can generate the internal shocks closer to the jet base. Depending on the disc geometry, there might be transient shocks in the jet and there might be multiple non-stationary shocks in the jet, which are of much interest in jet physics.

When very-high-energy gamma rays interact high in the Earth's atmosphere, they produce cascades of particles that induce flashes of Cherenkov light. Imaging Atmospheric Cherenkov Telescopes (IACTs) detect these flashes and convert them into shower images that can be analyzed to extract the properties of the primary gamma ray. The dominant background for IACTs is comprised of air shower images produced by cosmic hadrons, with typical noise-to-signal ratios of several orders of magnitude. The standard technique adopted to differentiate between images initiated by gamma rays and those initiated by hadrons is based on classical machine learning algorithms, such as Random Forests, that operate on a set of handcrafted parameters extracted from the images. Likewise, the inference of the energy and the arrival direction of the primary gamma ray is performed using those parameters. State-of-the-art deep learning techniques based on convolutional neural networks (CNNs) have the potential to enhance the event reconstruction performance, since they are able to autonomously extract features from raw images, exploiting the pixel-wise information washed out during the parametrization process. Here we present the results obtained by applying deep learning techniques to the reconstruction of Monte Carlo simulated events from a single, next-generation IACT, the Large-Sized Telescope (LST) of the Cherenkov Telescope Array (CTA). We use CNNs to separate the gamma-ray-induced events from hadronic events and to reconstruct the properties of the former, comparing their performance to the standard reconstruction technique. Three independent implementations of CNN-based event reconstruction models have been utilized in this work, producing consistent results.

In this paper, we have studied the convectively coupled equatorially trapped waves in rotating stars, with and without magnetic field. The equatorial trapped HD and MHD Poincar\'e, Rossby, mixed Rossby-Poincar\'e, and Kelvin waves were identified. The effects of stratification and non-traditional Coriolis force terms have been investigated. When the flow is strongly stratified, the wave frequencies of the convectively coupled model are almost the same as those of shallow water model. However, when the flow is weakly stratified, the wave frequencies are constrained by the buoyancy frequency. The non-traditional Coriolis terms affect the widths and phases of the equatorial waves. The width increases with the increasing non-traditional Coriolis parameter. Phase shift occurs when the non-traditional Coriolis parameter is included. Magnetic effect is significant when the magnetic field is strong. We have applied the model in the solar atmosphere and solar tachocline to explain the Rieger type periodicities. For the solar atmosphere, when magnetic effect is taken into account, we find that the magnetic field should be smaller than $5G$ in the solar photosphere. Otherwise, the Rieger type periodicities can be only attributed to long Rossby waves. For the solar tachocline, we find that magnetic field of the solar tachocline should be smaller than $50kG$ to observe the 160 days Rieger period. In addition, we find that the effect of the non-traditional Coriolis terms is not obvious in the solar photosphere, but its effect on the tachocline is significant.

Extragalactic black hole X-ray binaries LMC X-1 and LMC X-3 are the persistent sources, usually found in soft spectral state. We present the results from spectral and timing study of these sources using MAXI, NICER, NuSTAR and AstroSat observations carried out during 2014-2020. Study of long-term MAXI lightcurve shows that the fractional variability of flux in 2-10keV is moderate ($\sim$20%) in LMCX-1 and high ($\sim$50%) in LMCX-3. The energy spectra of both sources are characterized by a disc blackbody and a Comptonization component with LMCX-1 having an additional Fe-line emission feature. NICER (0.3-10keV), NuSTAR (3.0-40keV) and AstroSat (0.5-20keV) spectral analysis collectively show that LMC X-1 remained in the soft state (disc flux contribution $f_{disc}>80$%, photon index$\Gamma\sim2.06-4.08$) throughout 2014-2020. Mass accretion rate, $\dot{M}$ of LMC X-1 calculated from bolometric luminosity (0.1-50keV) is found to be within $0.07-0.24\dot{M}_{Edd}$ (Eddington mass accretion rate). Although LMC X-3 remained in the soft state ($f_{disc}>95\%,\Gamma\sim2.3$) during most of the time, it exhibits transition into intermediate ($f_{disc}=47-73\%,\Gamma\sim2.02-2.36$) and hard state ($f_{disc}\sim26\%,\Gamma\sim1.6$). $\dot{M}$ of LMC X-3 through different spectral states varies within $0.01-0.42\dot{M}_{Edd}$. Temporal study show that the Power Density Spectra (PDS) in 0.3-10keV follow a red-noise with rms of 2% for LMC X-1 and in case of LMC X-3, it is 0.08-2.35% during the soft state, but relatively high in the intermediate(3.05-4.91%) and hard state($\sim$17.06%). From continuum-fitting method we constrain spin of LMC X-1 and LMC X-3 to be within 0.85-0.94 and 0.16-0.33 respectively and from Fe-line fitting method, spin of LMC X-1 is found to be 0.93-0.94. Finally, we discuss the implication of our findings in the context of accretion disc dynamics around the vicinity of the BHs.

We report the first evidence of intrinsic alignment (IA) of red galaxies at $z>1$. We measure the gravitational shear-intrinsic ellipticity (GI) cross-correlation function at $z\sim1.3$ using galaxy positions from the FastSound spectroscopic survey and galaxy shapes from Canada-Hawaii-France telescope lensing survey data. Adopting the non-linear alignment model, we obtain a $2.4\sigma$ level detection of the IA amplitude, $A^{\rm LA}=27.48_{-11.54}^{+11.53}$, larger than the value extrapolated from the constraints obtained at lower redshifts. Our measured IA is translated into a $\sim 20\%$ contamination to the weak lensing power spectrum for the red galaxies. This marginal detection of IA for red galaxies at $z>1$ motivates the continuing investigation of the nature of IA for weak lensing studies. Furthermore, our result provides the first step to utilize IA measurements in future high-$z$ surveys as a cosmological probe, complementary to galaxy clustering and lensing.

In this work we present the results of a direct imaging survey for brown dwarf companions around the nearest stars at the mid-infrared 10 micron range ($\lambda_{c}$=8.7$\mu$m, $\Delta\lambda$=1.1$\mu$m) using the CanariCam instrument at the 10.4 m Gran Telescopio Canarias (GTC). We imaged the 25 nearest stellar systems within 5 pc of the Sun at declinations $\delta > -25^{\circ}$ (at least half have planets from radial velocity), reaching a mean detection limit of 11.3$\pm$0.2 mag (1.5 mJy) in the Si-2 8.7$\mu$m band over a range of angular separations from 1 to 10 arcsec. This would have allowed us to uncover substellar companions at projected orbital separations between $\sim$2 and 50 au, with effective temperatures down to 600 K and masses greater than 30 $M_{Jup}$ assuming an average age of 5 Gyr and down to the deuterium-burning mass limit for objects with ages $<$1 Gyr. From the non-detection of such companions, we determined upper limits on their occurrence rate at depths and orbital separations yet unexplored by deep imaging programs. For the M dwarfs, main components of our sample, we found with a 90% confidence level that less than 20% of these low-mass stars have L and T-type brown dwarf companions with $m \gtrsim 30 M_{Jup}$ and $T_{eff} \gtrsim$ 600 K at $\sim$3.5--35 au projected orbital separations.

The relationship between bipolar magnetic regions (BMRs) and their sunspots is an important property of the solar magnetic field, but it is not well constrained. One consequence is that it is a challenge for surface flux transport models (SFTMs) based on sunspot observations to determine the details of BMR emergence, which they require as input, from such data. We aimed to establish the relationship between the amount of magnetic flux in newly emerged BMRs and the area of the enclosed sunspots. Earlier attempts to constrain BMR magnetic flux were hindered by the fact that there is no proper database of the magnetic and physical properties of newly emerged BMRs currently available. We made use of the empirical model of the relationship between the disc-integrated facular and network magnetic flux and the total surface coverage by sunspots reported in a recent study. The structure of the model is such that it enabled us to establish, from these disc-integrated quantities, an empirical relationship between the magnetic flux and sunspot area of individual newly emerged BMRs, circumventing the lack of any proper BMR database. Applying the constraint on BMR magnetic flux derived here to an established SFTM retained its ability to replicate various independent datasets and the correlation between the model output polar field at the end of each cycle and the observed strength of the following cycle. The SFTM output indicates that facular and network magnetic flux rises with increasing sunspot magnetic flux at a slowing rate such that it appears to gradually saturate, analogous to earlier studies. The activity dependence of the ratio of facular and network flux to sunspot flux is consistent with the findings of recent studies: although the Sun is faculae-dominated, it is only marginally so as facular and network brightening and sunspot darkening appear to be closely balanced.

Most of the high-contrast imaging (HCI) data-processing techniques used over the last 15 years have relied on the angular differential imaging (ADI) observing strategy, along with subtraction of a reference point spread function (PSF) to generate exoplanet detection maps. Recently, a new algorithm called regime switching model (RSM) map has been proposed to take advantage of these numerous PSF-subtraction techniques; RSM uses several of these techniques to generate a single probability map. Selection of the optimal parameters for these PSF-subtraction techniques as well as for the RSM map is not straightforward, is time consuming, and can be biased by assumptions made as to the underlying data set. We propose a novel optimisation procedure that can be applied to each of the PSF-subtraction techniques alone, or to the entire RSM framework. The optimisation procedure consists of three main steps: (i) definition of the optimal set of parameters for the PSF-subtraction techniques using the contrast as performance metric, (ii) optimisation of the RSM algorithm, and (iii) selection of the optimal set of PSF-subtraction techniques and ADI sequences used to generate the final RSM probability map. The optimisation procedure is applied to the data sets of the exoplanet imaging data challenge (EIDC), which provides tools to compare the performance of HCI data-processing techniques. The data sets consist of ADI sequences obtained with three state-of-the-art HCI instruments: SPHERE, NIRC2, and LMIRCam. The results of our analysis demonstrate the interest of the proposed optimisation procedure, with better performance metrics compared to the earlier version of RSM, as well as to other HCI data-processing techniques.

Stripped-envelope supernovae (SE SNe) of Type Ib and Type Ic are thought to result from explosions of massive stars having lost their outer envelopes. The favoured explosion mechanism is by core-collapse, with the shock later revived by neutrino heating. However, there is an upper limit to the amount of radioactive Nickel-56 that such models can accomplish. Recent literature point to a tension between the maximum luminosity from such simulations and observations. We use a well characterized sample of SE SNe from the Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS) to collect a sample of spectroscopically classified normal Type Ibc SNe for which we use the ZTF light curves to determine the maximum luminosity. We cull the sample further based on data quality, light-curve shape, distance and colors. The methodology of the sample construction from this BTS sample can be used for many other future investigations. In total we use 129 Type Ib or Type Ic BTS SNe with an initial rough luminosity distribution peaked at Mr = -17.61 +- 0.72, and where 36% are apparently brighter than the theoretically predicted maximum brightness of Mr = -17.8. When we further cull this sample to ensure that the SNe are normal Type Ibc with good LC data within the Hubble flow, the sample of 94 objects has Mr = -17.64 +- 0.54. A main uncertainty in absolute magnitude determinations for SNe is the host galaxy extinction correction, but the reddened objects only get more luminous after corrections. If we simply exclude objects with red, unusual or uncertain colors, we are left with 14 objects at M = -17.90 +- 0.73, whereof a handful are most certainly brighter than the suggested theoretical limit. The main result of this study is thus that normal SNe Ibc do indeed reach luminosities above $10^{42.6}$ erg/s, apparently in conflict with existing explosion models.

The progress in the construction and operation of the Baikal Gigaton Volume Detector in Lake Baikal is reported. The detector is designed for search for high energy neutrinos whose sources are not yet reliably identified. It currently includes 2304 optical modules arranged on 64 strings, providing an effective volume of 0.4 km3 for cascades with energy above 100 TeV. We review the scientific case for Baikal-GVD, the construction plan, and first results from the partially built experiment, which is currently the largest neutrino telescope in the Northern Hemisphere and still growing up.

MAXI J1820+070 (ASSASN-18ey) is a Black hole X-ray binary discovered in 2018. The brightness of the source triggered multi-wavelength campaigns of this source from different observatories. We analyse the Power Density Spectra obtained from NICER high cadence observations of the source in the hard state. We obtain the evolution of the characteristic frequencies by modelling the PDS. We interpret the characteristic frequencies of various PDS components (both QPOs and broad band noise components) as variability occurring at a particular radius, and explain them in the context of the Relativistic Precession Model. We estimate the dimensionless spin of the black hole at $0.799^{+0.016}_{-0.015}$ by fitting the Relativistic Precession Model.

We discuss the prospects for identifying nearest isolated black holes (IBHs) in our Galaxy. IBHs accreting gas from the interstellar medium (ISM) likely form magnetically arrested disks (MADs). We show that thermal electrons in the MADs emit optical signals through the cyclo-synchrotron process while non-thermal electrons accelerated via magnetic reconnections emit a flat-spectrum synchrotron radiation in the X-ray to MeV gamma-ray ranges. The Gaia catalog will include at most a thousand of IBHs within $\lesssim 1$ kpc that are distributed on and around the cooling sequence of white dwarfs (WDs) in the Hertzsprung-Russell diagram. These IBH candidates should also be detected by eROSITA, with which they can be distinguished from isolated WDs and neutron stars. Followup observations with hard X-ray and MeV gamma-ray satellites will be useful to unambiguously identify IBHs.

We investigate a sample of seven edge-on spiral galaxies using Chandra observations. Edge-on spiral galaxies allow us to clearly separate source associated with their star-forming regions versus the outer edges of the system; offering a clear advantage over other systems. We uncover a number of X-ray point sources across these galaxies, and after eliminating contaminating foreground and background sources, we identify 12 candidate ultraluminous X-ray sources. All of these sources are projected onto the central regions, implying that the majority of ULXs in this sample of spiral galaxies are disk/bulge, and thus not halo sources. This also includes two transient ULXs, which may be long-duration transients and low mass X-ray binaries. This finding illustrates the need for further studies of transient ULXs.

Throughout an activity cycle, magnetic structures rise to the stellar surface, evolve and decay. Tracing their evolution on a stellar corona allows us to characterize the X-ray cycles. However, directly mapping magnetic structures is feasible only for the Sun, while such structures are spatially unresolved with present-day X-ray instruments on stellar coronae. We present here a new method, implemented by us, that indirectly reproduces the stellar X-ray spectrum and its variability with solar magnetic structures. The technique converts solar corona observations into a format identical to that of stellar X-ray observations and, specifically, XMM-Newton spectra. From matching these synthetic spectra with those observed for a star of interest, a fractional surface coverage with solar magnetic structures can be associated to each X-ray observation. We apply this method to two young solar-like stars: $\epsilon$ Eri ($\sim 400$ Myr), the youngest star to display a coronal cycle ($\sim 3$ yr), and Kepler 63 ($\sim 200$ Myr), for which the X-ray monitoring did not reveal a cyclic variability. We found that even during the cycle minimum a large portion of $\epsilon$ Eri's corona is covered with active structures. Therefore, there is little space for additional magnetic regions during the maximum, explaining the small observed cycle amplitude ($\Delta f \sim 0.12$) in terms of the X-ray luminosity. Kepler 63 displays an even higher coverage with magnetic structure than the corona of $\epsilon$ Eri. This supports the hypothesis that for stars younger than $<400$ Myr the X-ray cycles are inhibited by a massive presence of coronal regions.

To a first approximation, the microlensing phenomenon is achromatic, and great advancement has been achieved in the interpretation of the achromatic signals, which among other achievements has led to the discovery and characterization of well above $100$ new exoplanets. At higher order accuracy in the observations, microlensing has a chromatic component (a color term) which has so far been much less explored. Here, we analyze the chromatic microlensing effect of $4$ different physical phenomena, which have the potential to add important new knowledge about the stellar properties not easily reachable with other methods of observations. Our simulation is limited to the case of main-sequence source stars. Microlensing is particularly sensitive to giant and sub-giant stars near the Galactic center. While this population can be studied in short snapshots by use of the largest telescopes in the world, a general monitoring and characterization of the population can be achieved by use of more accessible medium-sized telescopes with specialized equipments through dual-color monitoring from observatories at sites with excellent seeing. We quantify our results to what will be achievable from the Danish $1.54$m telescope at La Silla observatory by use of the existing dual-color lucky imaging camera. Such potential monitoring programs of the bulge population from medium-sized telescopes include the characterization of starspots, limb-darkening, the frequency of close-in giant planet companions, and gravity darkening for blended source stars. We conclude our simulations with quantifying the likelihood of detecting these different phenomena per object where they are present to be $\sim 60\%$ and $\sim 30\%$ for the mentioned phenomena, when monitored during high magnification and caustic crossings, respectively.

Asteroids and other small bodies in the solar system tend to have irregular shapes, owing to their low gravity. This irregularity does not only apply to the topology, but also to the underlying geology, potentially containing regions of different densities and materials. The topology can be derived from optical observations, while the mass density distribution of an object is only observable, to some extent, in its gravitational field. In a companion paper, we presented geodesyNets, a neural network approach to infer the mass density distribution of an object from measurements of its gravitational field. In the present work, we apply this approach to the asteroid Bennu using real data from the Osiris Rex mission. The mission measured the trajectories of not only the Osiris Rex spacecraft itself, but also of numerous pebble-sized rock particles which temporarily orbited Bennu. From these trajectory data, we obtain a representation of Bennu's mass density and validate it by propagating, in the resulting gravity field, multiple pebbles not used in the training process. The performance is comparable to that of a polyhedral gravity model of uniform density, but does not require a shape model. As little additional information is needed, we see this as a step towards autonomous on-board inversion of gravitational fields.

Meteor physics can provide new clues about the size, structure, and density of cometary disintegration products, establishing a bridge between different research fields. From meteor magnitude data we have estimated the mass distribution of meteoroids from different cometary streams by using the relation between the luminosity and the mass obtained by Verniani (1973). These mass distributions are in the range observed for dust particles released from comets 1P/Halley and 81P/Wild 2 as measured from spacecraft. From the derived mass distributions, we have integrated the incoming mass for the most significant meteor showers. By comparing the mass of the collected Interplanetary Dust Particles (IDPs) with that derived for cometary meteoroids a gap of several orders of magnitude is encountered. The largest examples of fluffy particles are clusters of IDPs no larger than 100 micrometers in size (or 5x10^-7 g in mass) while the largest cometary meteoroids are centimeter-sized objects. Such gaps can be explained by the fragmentation in the interstellar medium or in the atmosphere of the original cometary particles. As an application of the mass distribution computations we describe the significance of the disruption of fragile comets in close approaches to Earth as a more efficient (and probably more frequent) way to deliver volatiles than direct impacts. We finally apply our model to quantify the flux of meteoroids from different meteoroid streams, and to describe the main physical processes contributing to the progressive decay of cometary meteoroids in the interplanetary medium.

The inner-most regions of circumbinary discs are unstable to a parametric instability whose non-linear evolution is hydrodynamical turbulence. This results in significant particle stirring, impacting on planetary growth processes such as the streaming instability or pebble accretion. In this paper, we present the results of three-dimensional, inviscid global hydrodynamical simulations of circumbinary discs with embedded particles of 1 cm size. Hydrodynamical turbulence develops in the disc, and we examine the effect of the particle back-reaction on vertical dust. We find that higher solid-to-gas ratios lead to smaller gas vertical velocity fluctuations, and therefore to smaller dust scale heights. For a metallicity $Z=0.1$, the dust scale height near the edge of the tidally-truncated cavity is $\sim 80\%$ of the gas scale height, such that growing a Ceres-mass object to a 10 $M_\oplus$ core via pebble accretion would take longer than the disc lifetime. Collision velocities for small particles are also higher than the critical velocity for fragmentation, which precludes grain growth and the possibility of forming a massive planetesimal seed for pebble accretion. At larger distances from the binary, turbulence is weak enough to enable not only efficient pebble accretion but also grain growth to sizes required to trigger the streaming instability. In these regions, an in-situ formation scenario of circumbinary planets involving the streaming instability to form a massive planetesimal followed by pebble accretion onto this core is viable. In that case, planetary migration has to be invoked to explain the presence of circumbinary planets at their observed locations.

Interplanetary Coronal Mass Ejections (ICMEs) are known to modify the structure of the solar wind as well as interact with the space environment of planetary systems. Their large magnetic structures have been shown to interact with galactic cosmic rays, leading to the Forbush decrease (FD) phenomenon. We revisit in the present article the 17 years of Advanced Composition Explorer spacecraft ICME detection along with two neutron monitors (McMurdo and Oulu) with a superposed epoch analysis to further analyze the role of the magnetic ejecta in driving FDs. We investigate in the following the role of the sheath and the magnetic ejecta in driving FDs, and we further show that for ICMEs without a sheath, a magnetic ejecta only is able to drive significant FDs of comparable intensities. Furthermore, a comparison of samples with and without a sheath with similar speed profiles enable us to show that the magnetic field intensity, rather than its fluctuations, is the main driver for the FD. Finally, the recovery phase of the FD for isolated magnetic ejecta shows an anisotropy in the level of the GCRs. We relate this finding at 1 au to the gradient of the GCR flux found at different heliospheric distances from several interplanetary missions.

Elemental abundances of the most metal-poor stars reflect the conditions in the early Galaxy and the properties of the first stars. We present a spectroscopic follow-up of two ultra metal-poor stars ([Fe/H]<-4.0) identified by the survey {\em Pristine}: Pristine 221.8781+9.7844 and Pristine 237.8588+12.5660 (hereafter Pr 221 and Pr 237, respectively). Combining data with earlier observations, we find a radial velocity of -149.25 $\pm$ 0.27 and -3.18 $\pm$ 0.19 km/s for Pr 221 and Pr 237, respectively, with no evidence of variability between 2018 and 2020. From a one-dimensional (1D) local thermodynamic equilibrium (LTE) analysis, we measure [Fe/H]$_{\rm LTE}$=-4.79 $\pm$ 0.14 for Pr 221 and [Fe/H]$_{\rm LTE}$=-4.22 $\pm$ 0.12 for Pr 237, in good agreement with previous studies. Abundances of Li, Na, Mg, Al, Si, Ca, Ti, Fe, and Sr were derived based on the non-LTE (NLTE) line formation calculations. When NLTE effects are included, we measure slightly higher metallicities: [Fe/H]$_{\rm NLTE}$=-4.40 $\pm$ 0.13 and [Fe/H]$_{\rm NLTE}$=-3.93 $\pm$ 0.12, for Pr 221 and Pr 237, respectively. Analysis of the G-band yields [C/Fe]$_{\rm 1D-LTE} \leq$ +2.3 and [C/Fe]$_{\rm 1D-LTE} \leq$ +2.0 for Pr 221 and Pr 237. Both stars belong to the low-carbon band. Upper limits on nitrogen abundances are also derived. Abundances for other elements exhibit good agreement with those of stars with similar parameters. Finally, to get insight into the properties of their progenitors, we compare NLTE abundances to theoretical yields of zero-metallicity supernovae. This suggests that the supernovae progenitors had masses ranging from 10.6 to 14.4 M$_{\odot}$ and low-energy explosions with 0.3-1.2 $\times$ 10$^{51}$ erg.

In this paper, we implement a new model-independent method to test the invariance of the speed of light $c$ as a function of redshift, by combining the measurements of galaxy cluster gas mass fraction, $H(z)$ from cosmic chronometers, and type-Ia supernovae (SNe Ia). In our analyses, we consider both a constant depletion factor (which corresponds to the ratio by which $f_{gas}$ is depleted with respect to the universal baryonic mean) and one varying with redshift. We also consider the influence of different $H_0$ estimates on our results. We look for a variation of $c$, given by $c(z)=c_0(1+c_1z)$. It is found a degeneracy between our final results on $c$ variation and the assumptions on gas mass fraction depletion factor. Most of our analyses indicate negligible variation of the speed of light.

Despite being bright ($V=12$) and nearby ($d=212$ pc) ASAS J071404+7004.3 has only recently been identified as a nova-like cataclysmic variable. We present time-resolved optical spectroscopy obtained at the Isaac Newton Telescope together with $\textit{Swift}$ X-ray and ultraviolet observations. We combined these with $\textit{TESS}$ photometry and find a period of 3.27h and a mass transfer rate of $4-9 \times 10^{-9} M_{sun}/yr$. Historical photometry shows at least one low state establishing the system as a VY Scl star. Our high-cadence spectroscopy also revealed rapidly changing winds emanating from the accretion disc. We have modelled these using the Monte Carlo PYTHON code and shown that all the emission lines could emanate from the wind - which would explain the lack of double-peaked lines in such systems. In passing,we discuss the effect of variability on the position of cataclysmic variables in the $\textit{Gaia}$ Hertzsprung-Russell diagram.

Many future small satellite missions are aimed to provide low-cost remote sensing data at unprecedented revisit rates, with a ground resolution of less than one meter. This requires high resolution, fast and sensitive line-scan imagers operating at low power consumption and ideally featuring spectral sensitivity. In this paper we present comprehensive characterization results of our 7 band Back-Side Illuminated (BSI) CCD-in-CMOS sensor with a pixel pitch of 5.4 um. We have extensively characterized the key performance parameters of our CCD-in-CMOS sensor, such as quantum efficiency (QE), full well capacity (FWC), read noise, conversion gain, non-linearity, dark current etc. Novelty of this device is the combination of 7 TDI bands on the same imager allowing simultaneous multispectral TDI capture. Glass-based broadband filters with a typical band-pass width of about 100 nm have been developed and glued together to form a filter assembly of 6 band-pass filters and one panchromatic channel. Multispectral capability of this sensor is particularly interesting for Low Earth Observation (LEO) applications such as environmental monitoring, precision agriculture, disaster detection and monitoring. To highlight its ad-vantages for use in vegetation observation, we demonstrated a fake leaf and a real leaf imaging using a 7 band BSI sensor with integrated filters operating in 7-band mode at 15 kHz.

The SPHERE project studies primary cosmic rays by detection of the Cherenkov light of extensive air showers reflected from the snowy surface of the earth. Measurements with the aerial-based detector SPHERE-2 were performed in 2011-2013. The detector was lifted by the balloon at altitudes up to 900 m above snowed surface of Lake Baikal, Russia. The results of the experiment are summarized now in a series of papers that opens with this article. An overview of the SPHERE-2 detector telemetry monitoring systems is presented along with the analysis of the measurements conditions including atmosphere profile. The analysis of the detector state and environment atmosphere conditions monitoring provided various cross-checks of detector calibration, positioning and performance.

Gravitational microlensing is one of the methods to detect exoplanets; planets outside our solar system. Here we focus on theoretical modeling of three lens systems and in particular circumbinary systems. Circumbinary systems include two stars and a planet and are estimated to make up a sizable portion of all exoplanets. Extending a method developed for binary lenses to the three lens case, we explore the parameter space of circumbinary systems producing exact magnification maps and light curves.

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, for example, a cutoff power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. Previous works have presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this solution that is both fast and accurate enough to be fitted to observed spectra is available as a model in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during the $\chi^2$-minimization which simplifies the application of the model considerably. We demonstrate this approach as well a study of possible parameter degeneracies with a comprehensive Markov-chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3. The derived accretion-flow structure features a small column radius of $\sim$63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

The division between stripped-envelope supernovae (SE-SNe) and superluminous supernovae (SLSNe) is not well defined in either photometric or spectroscopic space. While a sharp luminosity threshold has been suggested, there remains an increasing number of transitional objects that reach this threshold without the spectroscopic signatures common to SLSNe. In this work we present data and analysis on four SNe transitional between SE-SNe and SLSNe; the He-poor SNe 2019dwa and 2019cri, and the He-rich SNe 2019hge and 2019unb. Each object displays long-lived and variable photometric evolution with luminosities around the SLSN threshold of $M_r < -19.8$ mag. Spectroscopically however, these objects are similar to SE-SNe, with line velocities lower than either SE-SNe and SLSNe, and thus represent an interesting case of rare transitional events.

Although the sample of exoplanets in binaries has been greatly expanded, the sample heterogeneity and observational bias are obstacles toward a clear figure of exoplanet demographics in the binary environment. To overcome the obstacles, we conduct a statistical study that focuses on S-type planetary systems detected by the Radial Velocity method. We try to account for observational biases by estimating, from available RV data, planet detection efficiencies for each individual system. Our main results are as follows: (1) Single (multiple) planetary systems are mostly found in close (wide) binaries with separation aB<(>)100-300 AU. (2) In binaries, single and multiple-planet systems are similar in 1-D distributions of mass and period as well as eccentricity (in contrast to the "eccentricity dichotomy" found in single-star systems) but different in the 2-D period-mass diagram. Specifically, there is a rectangular-shaped gap in the period-mass diagram of ingle-planet systems but not for multiples. This gap also depends on binary separation and is more prominent in close binaries. (3) There is a rising upper envelope in the period-mass diagram for planets in wide binaries as well as in single stars but not in close binaries. More specifically, there is a population of massive short period planets in close binaries but almost absent in wide binaries or single stars. We suggest that enhanced planetary migration, collision, and /or ejection in close binaries could be the potential underlying explanation for these three features.

This paper studies the magnetic topology of successively erupting active regions (ARs) 11429 and 12371. Employing vector magnetic field observations from Helioseismic and Magnetic Imager, the pre-eruptive magnetic structure is reconstructed by a model of non-linear force-free field (NLFFF). For all the five CMEs from these ARs, the pre-eruptive magnetic structure identifies an inverse-S sigmoid consistent with the coronal plasma tracers in EUV observations. In all the eruption cases, the quasi-separatrix layers (QSLs) of Large Q values are continuously enclosing core field bipolar regions in which inverse-S shaped flare ribbons are observed. These QSLs essentially represent the large connectivity gradients between the domains of twisted core flux within the inner bipolar region and the surrounding potential like arcade. It is consistent with the observed field structure largely with the sheared arcade. The QSL maps in the chromosphere are compared with the flare-ribbons observed at the peak time of the flares. The flare ribbons are largely inverse-S shape morphology with their continuity of visibility is missing in the observations. For the CMEs in the AR 12371, the QSLs outline the flare ribbons as a combination of two inverse J-shape sections with their straight parts being separated. These QSLs are typical with the weakly twisted flux rope. Similarly, for the CMEs in the AR 11429, the QSLs are co-spatial with the flare ribbons both in the middle of the PIL and in the hook sections. In the frame work of standard model of eruptions, the observed flare ribbons are the characteristic of the pre-eruptive magnetic structure being sigmoid which is reproduced by {\bf the} NLFFF model with a weakly twisted flux rope at the core.

The Orion Nebula Cluster (ONC) is the nearest dense star-forming region at $\sim$400 pc away, making it an ideal target to study the impact of high stellar density and proximity to massive stars (the Trapezium) on protoplanetary disk evolution. The OMC1 molecular cloud is a region of high extinction situated behind the Trapezium in which actively forming stars are shielded from the Trapezium's strong radiation. In this work, we survey disks at high resolution with ALMA at three wavelengths with resolutions of 0.095\arcsec (3 mm; Band 3), 0.048\arcsec (1.3 mm; Band 6), and 0.030\arcsec (0.85 mm; Band 7) centered on radio Source I. We detect 127 sources, including 15 new sources that have not previously been detected at any wavelength. 72 sources are spatially resolved at 3 mm, with sizes from $\sim$8 - 100 AU. We classify 76 infrared-detected sources as foreground ONC disks and the remainder as embedded OMC1 disks. The two samples have similar disk sizes, but the OMC1 sources have a dense and centrally concentrated spatial distribution, indicating they may constitute a spatially distinct subcluster. We find smaller disk sizes and a lack of large (>75 AU) disks in both our samples compared to other nearby star-forming regions, indicating that environmental disk truncation processes are significant. While photoevaporation from nearby massive Trapezium stars may account for the smaller disks in the ONC, the embedded sources in OMC1 are hidden from this radiation and thus must truncated by some other mechanism, possibly dynamical truncation or accretion-driven contraction.

Transmission spectroscopy of transiting exoplanets is a proven technique that can yield information on the composition and structure of a planet's atmosphere. However, transmission spectra may be compromised by inhomogeneities in the stellar photosphere. The sub-Neptune-sized habitable zone planet K2-18 b has water absorption detected in its atmosphere using data from the Hubble Space Telescope (HST). Herein, we examine whether the reported planetary atmospheric signal seen from HST transmission spectroscopy of K2-18 b could instead be induced by time-varying star spots. We built a time-variable spectral model of K2-18 that is designed to match the variability amplitude seen in K2 photometric data, and used this model to simulate 1000 HST data-sets that follow the K2-18 b observation strategy. More than 1% of these provide a better fit to the data than the best-fitting exoplanet atmosphere model. After resampling our simulations to generate synthetic HST observations, we find that 40% of random draws would produce an atmospheric detection at a level at least as significant as that seen in the actual HST data of K2-18 b. This work illustrates that the inferred detection of an atmosphere on K2-18 b may alternatively be explained by stellar spectral contamination due to the inhomogeneous photosphere of K2-18. We do not rule out a detection of water in the planet's atmosphere, but provide a plausible alternative that should be considered, and conclude that more observations are needed to fully rule out stellar contamination.

We calculate the quasi-normal mode complex frequencies of the Kerr-Newman black hole with arbitrary values of spin and charge, for the modes typically dominant during a binary black hole coalescence, $(\ell,m,n) = \{(2,2,0), (2,2,1), (3,3,0) \}$. Building analytical fits of the black hole spectrum, we construct a template to model the post-merger phase of a binary black hole coalescence in the presence of a remnant $U(1)$ charge. Aside from astrophysical electric charge, our template can accommodate extensions of the Standard Model, such as a dark photon. Applying the model to LIGO-Virgo detections, we find that we are unable to distinguish between the charged and uncharged hypotheses from a purely post-merger analysis of the current events. However, restricting the mass and spin to values compatible with the analysis of the full signal, we obtain a 90th percentile bound $\bar{q} < 0.33$ on the black hole charge-to-mass ratio, for the most favorable case of GW150914. Under similar assumptions, by simulating a typical loud signal observed by the LIGO-Virgo network at its design sensitivity, we assess that this model can provide a robust measurement of the charge-to-mass ratio only for values $\bar{q} \gtrsim 0.5$; here we also assume that the mode amplitudes are similar to the uncharged case in creating our simulated signal. Lower values, down to $\bar{q} \sim 0.3$, could instead be detected when evaluating the consistency of the pre-merger and post-merger emission.

Valuable theoretical predictions of nuclear dipole excitations in the whole nuclear chart are of great interest for different applications, including in particular nuclear astrophysics. We present here the systematic study of the electric dipole (E1) photon strength functions (PSFs) combining the microscopic Hartree-Fock-Bogoliubov plus Quasiparticle Random Phase Approximation (HFB+QRPA) model and the parametrizations constrained by the available experimental giant dipole resonance (GDR) data. For about 10000 nuclei with 8<Z<124 lying between the proton and the neutron drip-lines on nuclear chart, the particle-hole strength distributions are computed using the HFB+QRPA model under the assumption of spherical symmetry and making use of the BSk27 Skyrme effective interaction derived from the most accurate HFB mass model (HFB-27) so far achieved. Large-scale calculations of the BSk27+QRPA E1 PSFs are performed in the framework of a specific folding procedure, in which three phenomenological improvements are considered. First, two interference factors are introduced and adjusted to reproduce at best the available experimental GDR data. Second, an empirical expression accounting for the deformation effect is applied to describe the peak splitting of the strength function. Third, the width of the strength function is corrected by a temperature-dependent term, which effectively increases the de-excitation photon strength function at low-energy. The E1 PSFs as well as the extracted GDR peaks and widths are compared with available experimental data. A relatively good agreement with data indicates the reliability of the calculations. Eventually, the astrophysical (n,g) rates for all the 10000 nuclei with 8<Z<124 are estimated using the present E1 PSFs. The resulting reaction rates are compared with previous BSk7+QRPA results and Gogny-HFB+QRPA predictions based on the D1M interaction.

We calculate the energy-dependent cross section of the $np\leftrightarrow d\gamma$ process in chiral effective field theory and apply state-of-the-art tools for quantification of theory uncertainty. We focus on the low-energy regime, where the magnetic dipole and the electric dipole transitions cross over, including the range relevant for big-bang nucleosynthesis. Working with the leading one- and two-body electromagnetic currents, we study the order-by-order convergence of this observable in the chiral expansion of the nuclear potential. We find that the Gaussian process error model describes the observed convergence very well, allowing us to present Bayesian credible intervals for the truncation error with correlations between the cross sections at different energies taken into account. We obtain a 1$\sigma$ estimate of about 0.2\% for the uncertainty from the truncation of the nuclear potential. This is an important step towards calculations with statistically interpretable uncertainties for astrophysical reactions involving light nuclei.

The search of violation of Lorentz symmetry, or Lorentz violation (LV), is an active research field. The effects of LV are expected to be very small and special systems are often used to search it. High-energy astrophysical neutrinos offer a unique system to search signatures of LV due to the three factors: high neutrino energy, long propagation distance, and the presence of quantum mechanical interference. In this brief review, we introduce tests of LV and summarize existing searches of LV using atmospheric and astrophysical neutrinos.

We calculate fast conversions of two flavor neutrinos by considering Boltzmann collisions of neutrino scatterings. In an idealized angular distribution of neutrinos with electron-lepton number crossing, we find that the collision terms of the neutrino scattering enhance the transition probability of fast flavor conversions as in the previous study. We analyze the dynamics of fast flavor conversions with collisions in detail based on the motion of polarization vectors in cylindrical coordinate analogous to a pendulum motion. The phase of the polarization vector is partially synchronized, and the phase deviation from the Hamiltonian governs the dynamics. The collision terms break a closed orbit and gradually make the phase space smaller. The flavor conversions are enhanced during this limit cycle. After the significant flavor conversion, all of the neutrino polarization vectors start to align with the z-axis owing to the collision effect within the time scale of the collision term irrespective of neutrino scattering angles. Though our analysis does not fully understand the dynamics of fast flavor conversion, the framework gives a new insight into this complicated phenomenon in further study.

Extreme-mass-ratio inspirals will be prized sources for the upcoming space-based gravitational-wave observatory LISA. The hunt for these is beset by many open theoretical and computational problems in both source modeling and data analysis. We draw attention here to one of the most poorly understood: the phenomenon of non-local correlations in the space of extreme-mass-ratio-inspiral signals. Such correlations are ubiquitous in the continuum of possible signals (degeneracy), and severely hinder the search for actual signals in LISA data. However, they are unlikely to manifest in a realistic set of putative signals (confusion). We develop an inventory of new analysis tools in order to conduct an extensive qualitative study of degeneracy - its nature, causes, and implications. Previously proposed search strategies for extreme-mass-ratio inspirals are reviewed in the light of our results, and additional guidelines are suggested for the scientific analysis of such sources.

Owing to it's non-minimal nature of coupling with gravity, the scalar field in a scalar-tensor theory is covariantly conserved with more than one choice of its energy-momentum tensor (EMT). Apparently this ambiguity in the choice of the symmetric EMT results from the way the conserved EMT is identified by different algebraic manipulations of the gravity field equation in such non-minimally coupled (NMC) theories. In this paper, we demonstrate that %using , one can arrive at these different choices of EMT, not by mere algebraic manipulation, but by demanding the equivalence of a NMC theory in Jordan frame with a minimally coupled theory in Einstein frame at the level of physical laws.

One prominent feature in the atmospheres of Jupiter and Saturn is the appearance of large-scale vortices. However, the sustaining mechanism of these large-scale vortices remains unclear. One possible mechanism is that these large-scale vortices are driven by rotating convection. Here we present numerical simulation results on rapidly rotating Rayleigh-B\'enard convection at a small Prandtl number $Pr=0.1$ (close to the turbulent Prandtl numbers of Jupiter and Saturn). We have identified four flow regimes in our simulation: multiple small vortices, coexisted large-scale cyclone and anticyclone, large-scale cyclone, and turbulence. The formation of large-scale vortices requires two conditions to be satisfied: the vertical Reynolds number is large ($Re_{z}\ge 400$), and the Rossby number is small ($Ro\leq 0.4$). Large-scale cyclone first appears when $Ro$ decreases to be smaller than 0.4. When $Ro$ further decreases to be smaller than 0.1, coexisted large-scale anticyclone emerges. We have studied the heat transfer in rapidly rotating convection. The result reveals that the heat transfer is more efficient in the anticyclonic region than in the cyclonic region. Besides, we find that 2D effect increases and 3D effect decreases in transporting convective flux as rotation rate increases. We find that aspect ratio has an effect on the critical Rossby number for the emergence of large-scale vortices. Our results provide helpful insights on understanding the dynamics of large-scale vortices in gas giants.

The study of stellar burning began just over 100 years ago. Nonetheless, we do not yet have a detailed picture of the nucleosynthesis within stars and how nucleosynthesis impacts stellar structure and the remnants of stellar evolution. Achieving this understanding will require precise direct measurements of the nuclear reactions involved. This report summarizes the status of direct measurements for stellar burning, focusing on developments of the last couple of decades, and offering a prospectus of near-future developments.

A search for magnetised quark nuggets (MQN) is reported using acoustic signals from hydrophones placed in the Great Salt Lake (GSL) in the USA. No events satisfying the expected signature were seen. This observation allows limits to be set on the flux of MQNs penetrating the atmosphere and depositing energy in the GSL. The expected signature of the events was derived from pressure pulses caused by high-explosive cords between the lake surface and bottom at various locations in the GSL. The limits obtained from this search are compared with those obtained from previous searches and are compared to models for the formation of MQNs.

The mass and spin properties of binary black holes (BBHs) inferred from their gravitational-wave signatures reveal important clues about how these systems form. BBHs originating from isolated binary evolution are expected to have spins preferentially aligned with their orbital angular momentum, whereas there is no such preference in binaries formed via dynamical assembly. The fidelity with which near-future gravitational-wave detectors can measure off-axis spins will have implications for the study of BBH formation channels. In this work, we examine the degree to which the Advanced LIGO Plus (A+) and Advanced Virgo Plus (AdV+) interferometric detectors can measure both aligned and misaligned spins. We compare spin resolution between the LIGO-Virgo network operating at either A+/AdV+ (``Plus'') sensitivity or Advanced-era design (``Design'') sensitivity using simulated BBH gravitational-wave signals injected into synthetic detector noise. The signals are distributed over the mass-spin parameter space of likely BBH systems, accounting for the effects of precession and higher-order modes. We find that the Plus upgrades yield significant improvements in spin estimation for systems with unequal masses and moderate or large spins. Using simulated signals modelled after different types of hierarchical BBH mergers, we also conclude that the Plus detector network will yield substantially improved spin estimates for 1G+2G binaries compared to the Design network.

The tau neutrino is the least well measured particle in the Standard Model. Most notably, the tau neutrino row of the lepton mixing matrix is quite poorly constrained when unitarity is not assumed. In this paper we identify several new and overlooked data sets involving tau neutrinos that improve our understanding of the tau neutrino part of the mixing matrix. We present new results on the unitarity of the tau row leveraging existing constraints on the electron and muon rows for the cases of unitarity violation, with and without kinematically accessible steriles. We also show the expected sensitivity due to upcoming experiments and demonstrate that the tau neutrino row precision is expected to be comparable to the muon neutrino row in a careful combined fit.