Search

Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the polarized sky patches of the size of around ten-by-ten arcminutes, at a distance of up to 20 degrees from the Moon, on bright-Moon nights as a wide-field linear polarimetric flat source. We observed 19 patches of the sky adjacent to the bright-Moon with the RoboPol instrument in the SDSS-r broadband filter. These were observed on five nights within two days of the full-Moon across two RoboPol observing seasons. We find that for 18 of the 19 patches, the uniformity in the measured normalized Stokes parameters $q$ and $u$ is within 0.2 %, with 12 patches exhibiting uniformity within 0.07 % or better for both $q$ and $u$ simultaneously, making them reliable and stable wide-field linear polarization flats. We demonstrate that the sky on bright-Moon nights is an excellent wide-field linear polarization flat source. Various combinations of the normalized Stokes parameters $q$ and $u$ can be obtained by choosing suitable locations of the sky patch with respect to the Moon, Comment: 8 pages including appendix, 6 figures and 3 tables. Submitted to Astronomy and Astrophysics for review. Comments are welcome

We present an analysis of the theoretical and observed light curve parameters of the fundamental mode (FU) classical Cepheids in the Magellanic Clouds in $V$- and $I$- photometric bands. The state-of-the-art 1D non-linear radial stellar pulsation (RSP) code in MESA (\textsc{mesa-rsp}) has been utilized to generate the theoretical light curves using four sets of convection parameters. Theoretical light curves with two chemical compositions: $Z=0.008$ and $Z=0.004$ appropriate for the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC), respectively, covered a wide range of periods ($3

1$) and all periods. The multiphase relations obtained from theoretical and observed light curves in the PL/PC/AC plane are found to be dynamic in nature, with the effect more pronounced at $\Phi \sim 0.75-0.85$. Furthermore, a contrasting behaviour of the theoretical/observed multiphase PL and PC relations between the short and long periods has been found for both LMC and SMC. The analysis shows that multiphase PL relations are more stringent to test the models with observations over the FPs. Distances to the LMC/SMC determined using long period Cepheids are found to be in good agreement with the literature values when the term $R_{21}$ is added to the PL relation., Comment: 17 pages, 11 Figures and 9 Tables. Accepted for publication by MNRAS

Wildfires pose a significant threat to the safe and reliable operations of the electric grid. To mitigate wildfire risk, system operators resort to public safety power shutoffs, or PSPS, that shed load for a subset of customers. As wildfire risk forecasts are stochastic, such decision-making may often be sub-optimal. This paper proposes a two-stage topology control problem that jointly minimizes generation and load-shedding costs in the face of uncertain fire risk. Compared to existing work, we include preand post-event topology control actions and consider scenarios where the wildfire risk is known with low and high confidence. The effectiveness of the proposed approach is demonstrated using a benchmark test system, artificially geo-located in Southern California, and using stochastic wildfire risk data that exists in the literature. Our work provides a crucial study of the comparative benefits of pre-event versus post-event control and the effects of wildfire risk accuracy on each control strategy., Comment: 6 pages, 6 figures

We present new results on community recovery based on the PageRank Nibble algorithm on a sparse directed stochastic block model (dSBM). Our results are based on a characterization of the local weak limit of the dSBM and the limiting PageRank distribution. This characterization allows us to estimate the probability of misclassification for any given connection kernel and any given number of seeds (vertices whose community label is known). The fact that PageRank is a local algorithm that can be efficiently computed in both a distributed and asynchronous fashion, makes it an appealing method for identifying members of a given community in very large networks where the identity of some vertices is known., Comment: PageRank Nibble, directed stochastic block model, local weak convergence, community detection

The number counts of homogeneous samples of radio sources are a tried and true method of probing the large scale structure of the Universe, as most radio sources outside the galactic plane are at cosmological distances. As such they are expected to trace the cosmic radio dipole, an anisotropy analogous to the dipole seen in the cosmic microwave background (CMB). Results have shown that although the cosmic radio dipole matches the direction of the CMB dipole, it has a significantly larger amplitude. This result challenges our assumption of the Universe being isotropic, which can have large repercussions for the current cosmological paradigm. Though significant measurements have been made, sensitivity to the radio dipole is generally hampered by systematic effects that can cause large biases in the measurement. Here we assess these systematics with data from the MeerKAT Absorption Line Survey (MALS). We present the analysis of ten MALS pointings, focusing on systematic effects that could lead to an inhomogeneous catalogue. We describe the calibration and creation of full band continuum images and catalogues, producing a combined catalogue containing 16,313 sources and covering 37.5 square degrees of sky down to a sensitivity of 10 $\mu$Jy/beam. We measure the completeness, purity, and flux recovery statistics for these catalogues using simulated data. We investigate different source populations in the catalogues by looking at flux densities and spectral indices, and how they might influence source counts. Using the noise characteristics of the pointings, we find global measures that can be used to correct for the incompleteness of the catalogue, producing corrected number counts down to 100 - 200 $\mu$Jy. We show that we can homogenise the catalogues and properly account for systematic effects. We determine that we can measure the dipole to $3\sigma$ significance with 100 MALS pointings., Comment: 30 pages, 25 figures. Accepted for publication in Astronomy & Astrophysics

The Gaussian Process (GP) based Chance-Constrained Optimal Power Flow (CC-OPF) is an open-source Python code developed for solving economic dispatch (ED) problem in modern power grids. In recent years, integrating a significant amount of renewables into a power grid causes high fluctuations and thus brings a lot of uncertainty to power grid operations. This fact makes the conventional model-based CC-OPF problem non-convex and computationally complex to solve. The developed tool presents a novel data-driven approach based on the GP regression model for solving the CC-OPF problem with a trade-off between complexity and accuracy. The proposed approach and developed software can help system operators to effectively perform ED optimization in the presence of large uncertainties in the power grid., Comment: 6 pages, 2 figures

We obtain local weak limits in probability for Collapsed Branching Processes (CBP), which are directed random networks obtained by collapsing random-sized families of individuals in a general continuous-time branching process. The local weak limit of a given CBP, as the network grows, is shown to be a related continuous-time branching process stopped at an independent exponential time. This is done through an explicit coupling of the in-components of vertices with the limiting object. We also show that the in-components of a finite collection of uniformly chosen vertices locally weakly converge (in probability) to i.i.d. copies of the above limit, reminiscent of propagation of chaos in interacting particle systems. We obtain as special cases novel descriptions of the local weak limits of directed preferential and uniform attachment models. We also outline some applications of our results for analyzing the limiting in-degree and PageRank distributions., Comment: 23 pages, 3 figures

Humans naturally exploit haptic feedback during contact-rich tasks like loading a dishwasher or stocking a bookshelf. Current robotic systems focus on avoiding unexpected contact, often relying on strategically placed environment sensors. Recently, contact-exploiting manipulation policies have been trained in simulation and deployed on real robots. However, they require some form of real-world adaptation to bridge the sim-to-real gap, which might not be feasible in all scenarios. In this paper we train a contact-exploiting manipulation policy in simulation for the contact-rich household task of loading plates into a slotted holder, which transfers without any fine-tuning to the real robot. We investigate various factors necessary for this zero-shot transfer, like time delay modeling, memory representation, and domain randomization. Our policy transfers with minimal sim-to-real gap and significantly outperforms heuristic and learnt baselines. It also generalizes to plates of different sizes and weights. Demonstration videos and code are available at https://sites.google.com/view/compliant-object-insertion., Comment: Accepted for publication at 2023 IEEE International Conference on Robotics and Automation (ICRA)

No direct experimental constraints exist on Leptoquark (LQ) couplings with quarks and right-handed neutrinos (RHNs). If a LQ dominantly couples to RHNs, it can leave unique signatures at the LHC. The RHNs can be produced copiously from LQ decays as long as they are lighter than the LQs. LQ-induced RHN production has never been searched for in experiments. This channel can act as a simultaneous probe for RHNs and LQs that dominantly couple to RHNs. In this paper, we consider all possible charge-$2/3$ and $1/3$ scalar and vector LQs that dominantly couple to second-generation quarks and RHN. We study the pair and single productions of TeV-scale LQs and their subsequent decay to sub-TeV RHNs, realised in the inverse seesaw framework. We also consider RHN pair production through a $t$-channel LQ exchange. The single LQ production and $t$-channel contributions can be significant for large LQ-RHN-quark couplings. We systematically combine events from these processes leading to a pair of RHNs plus jets to study the prospects of LQ-assisted RHN pair production. We analyse the monolepton and opposite-sign dilepton final states and estimate the discovery reach at the high-luminosity LHC., Comment: 11 pages, 4 figures, 3 tables. Matches with the journal version

Partitions wherein the even parts appear in two different colours are known as cubic partitions. Recently, Merca introduced and studied the function $A(n)$, which is defined as the difference between the number of cubic partitions of $n$ into an even number of parts and the number of cubic partitions of $n$ into an odd number of parts. In particular, using Smoot's \textsf{RaduRK} Mathematica package, Merca proved the following congruences by finding the exact generating functions of the respective functions. For all $n\ge0$, \begin{align*}A(9n+5)\equiv 0\pmod 3,\\ A(27n+26)\equiv 0\pmod 3. \end{align*} By using generating function manipulations and dissections, da Silva and Sellers proved these congruences and two infinite families of congruences modulo 3 arising from these congruences. In this paper, by employing Ramanujan's theta function identities, we present simplified formulas of the generating functions from which proofs of the congruences of Merca as well as those of da Silva and Sellers follow quite naturally. We also study analogous partition functions wherein multiples of $k$ appear in two different colours, where $k\in\{3,5,7,23\}$., Comment: 20 pages

Precision observables are well known for constraining most of the Beyond Standard Model (BSM) scenarios tightly. We present here a simple and comprehensive fitting framework for various BSM scenarios to these observables. We start with the fit of $S$, $T$ and $V$ parameter and their correlations using the Electroweak Precision Observables (EWPO) including the recent $m_W$ measurement from CDF-II. Utilizing these observables, we also fit various New Physics (NP) scenarios consisting of different subsets of dimension-6 Standard Model Effective Field Theory (SMEFT) operators in the Warsaw basis out of a total of 10 appearing at tree level in EWPO. To further constrain these scenarios, we augment these observables with $\Delta_{CKM}$ measurement using 1-loop matching of the Low Energy Effective Field Theory (LEFT) to SMEFT operators at the Z-pole. We show that the inclusion of $\Delta_{CKM}$ constraint indeed results in stronger bounds on the SMEFT Wilson Coefficients. We also constrain the UV parameters of BSM extensions like Vectorlike leptons (VLL) and find out that such a minimal extension is in tension with the forward-backward asymmetry in $b$-sector ($A_b^{FB}$) and the recent measurement of $M_W$. In order to lift the two blind directions, which one encounters while fitting all the 10 SMEFT WCs at tree-level, we also include the LEP-II observables pertaining to the $WW$ production and present the results for the fits with and without $\Delta_{CKM}$ constraint., Comment: 31 pages, 6 figures, 19 tables

We report the detection of stimulated hydrogen radio recombination line (RRL) emission from ionized gas in a $z=0.89$ galaxy using 580--1670 MHz observations from the MeerKAT Absorption Line Survey (MALS). The RRL emission originates in a galaxy that intercepts and strongly lenses the radio blazar PKS 1830-211 ($z=2.5$). This is the second detection of RRLs outside of the local universe and the first clearly associated with hydrogen. We detect effective H144$\alpha$ (and H163$\alpha$) transitions at observed frequencies of 1156 (798) MHz by stacking 17 (27) RRLs with 21$\sigma$ (14$\sigma$) significance. The RRL emission contains two main velocity components and is coincident in velocity with HI 21 cm and OH 18 cm absorption. We use the RRL spectral line energy distribution and a Bayesian analysis to constrain the density ($n_e$) and the volume-averaged pathlength ($\ell$) of the ionized gas. We determine $\log( n_e ) = 2.0_{-0.7}^{+1.0}$ cm$^{-3}$ and $\log( \ell ) = -0.7_{-1.1}^{+1.1}$ pc towards the north east (NE) lensed image, likely tracing the diffuse thermal phase of the ionized ISM in a thin disk. Towards the south west (SW) lensed image, we determine $\log( n_e ) = 3.2_{-1.0}^{+0.4}$ cm$^{-3}$ and $\log( \ell ) = -2.7_{-0.2}^{+1.8}$ pc, tracing gas that is more reminiscent of H II regions. We estimate a star formation (surface density) rate of $\Sigma_{\mathrm{SFR}} \sim 0.6$ M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$ or SFR $\sim 50$ M$_{\odot}$ yr$^{-1}$, consistent with a star-forming main sequence galaxy of $M_{\star} \sim 10^{11}$ M$_{\odot}$. The discovery presented here opens up the possibility of studying ionized gas at high redshifts using RRL observations from current and future (e.g., SKA and ngVLA) radio facilities., Comment: 17 pages, 10 figures, 2 tables, accepted in ApJ

We learn optimal instance-specific heuristics for the global minimization of nonconvex quadratically-constrained quadratic programs (QCQPs). Specifically, we consider partitioning-based convex mixed-integer programming relaxations for nonconvex QCQPs and propose the novel problem of strong partitioning to optimally partition variable domains without sacrificing global optimality. Since solving this max-min strong partitioning problem exactly can be very challenging, we design a local optimization method that leverages generalized gradients of the value function of its inner-minimization problem. However, even solving the strong partitioning problem to local optimality can be time-consuming. To address this, we propose a simple and practical machine learning (ML) approximation for homogeneous families of QCQPs. We conduct a detailed computational study on randomly generated QCQP families, including instances of the pooling problem, using the open-source global solver Alpine. Numerical experiments demonstrate that our ML approximation of strong partitioning reduces Alpine's solution time by a factor of 2 to 4.5 on average, with a maximum reduction factor of 10 to 200 across the different QCQP families.

We introduce the Conditional Independence Regression CovariancE (CIRCE), a measure of conditional independence for multivariate continuous-valued variables. CIRCE applies as a regularizer in settings where we wish to learn neural features $\varphi(X)$ of data $X$ to estimate a target $Y$, while being conditionally independent of a distractor $Z$ given $Y$. Both $Z$ and $Y$ are assumed to be continuous-valued but relatively low dimensional, whereas $X$ and its features may be complex and high dimensional. Relevant settings include domain-invariant learning, fairness, and causal learning. The procedure requires just a single ridge regression from $Y$ to kernelized features of $Z$, which can be done in advance. It is then only necessary to enforce independence of $\varphi(X)$ from residuals of this regression, which is possible with attractive estimation properties and consistency guarantees. By contrast, earlier measures of conditional feature dependence require multiple regressions for each step of feature learning, resulting in more severe bias and variance, and greater computational cost. When sufficiently rich features are used, we establish that CIRCE is zero if and only if $\varphi(X) \perp \!\!\! \perp Z \mid Y$. In experiments, we show superior performance to previous methods on challenging benchmarks, including learning conditionally invariant image features., Comment: ICLR 2023

The BL Lac PKS 1413+135 was observed by the Large Survey Project "MeerKAT Absorption Line Survey" (MALS) in the L-band, at 1139 MHz and 1293-1379 MHz, targeting the HI and OH lines in absorption at z = 0.24671. The radio continuum is thought to come from a background object at redshift lower than 0.5, as suggested by the absence of gravitational images. The HI absorption line is detected at high signal-to-noise, with a narrow central component, and a red wing, confirming previous results. The OH 1720 MHz line is clearly detected in (maser) emission, peaking at a velocity shifted by -10 to -15 km/s with respect to the HI peak. The 1612 MHz line is lost due to radio interferences. The OH 1667 MHz main line is tentatively detected in absorption, but not the 1665 MHz one. Over 30 years, a high variability is observed in optical depths, due to the rapid changes of the line of sight, caused by the superluminal motions of the radio knots. The HI line has varied by 20 per cent in depth, while the OH-1720 MHz depth has varied by a factor 4. The position of the central velocity and the widths also varied. The absorbing galaxy is an early-type spiral (maybe S0) seen edge-on, with a prominent dust lane, covering the whole disk. Given the measured mass concentration, and the radio continuum size at centimeter wavelengths (100 mas corresponding to 400 pc at z = 0.25), the width of absorption lines from the nuclear regions are expected up to 250 km/S. The narrowness of the observed lines (< 15 km/s) suggest that the absorption comes from an outer gas ring, as frequently observed in S0 galaxies. The millimetric lines are even narrower (< 1 km/s), which corresponds to the continuum size restricted to the core. The core source is covered by individual 1 pc molecular clouds, of column density a few 10^22 cm-2, which is compatible with the gas screen detected in X-rays., Comment: 7 pages, 6 figures, accepted in A&A

The anisotropic diffusion equation is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and its interplay with the ambient magnetic field. This diffusion term contributes to the highly stiff nature of the CR transport equation. In order to conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods that directly compute the exponential of the matrix to solve the linear anisotropic diffusion equation. These methods allow us to take even larger step sizes; in certain cases, we are able to choose a step size as large as the simulation time, i.e., only one time step. This can substantially speed-up the simulations whilst generating highly accurate solutions (l2 error $\leq 10^{-10}$). Additionally, we test an approach based on extracting a constant diffusion coefficient from the anisotropic diffusion equation, where the constant coefficient term is solved implicitly or exponentially and the remainder is treated using some explicit method. We find that this approach, for homogeneous linear problems, is unable to improve on the exponential-based methods that directly evaluate the matrix exponential., Comment: in submission

Krylov-based algorithms have long been preferred to compute the matrix exponential and exponential-like functions appearing in exponential integrators. Of late, direct polynomial interpolation of the action of these exponential-like functions have been shown to be competitive with the Krylov methods. We analyse the performance of the state-of-the-art Krylov algorithm, KIOPS, and the method of polynomial interpolation at Leja points for a number of exponential integrators for various test problems and with varying amounts of stiffness. Additionally, we investigate the performance of an iterative scheme that combines both the KIOPS and Leja approach, named LeKry, that shows substantial improvements over both the Leja- and Krylov-based methods for certain exponential integrators. Whilst we do manage to single out a favoured iterative scheme for each of the exponential integrators that we consider in this study, we do not find any conclusive evidence for preferring either KIOPS or Leja for different classes of exponential integrators. We are unable to identify a superior exponential integrator, one that performs better than all others, for most, if not all of the problems under consideration. We, however, do find that the performance significantly depends on the interplay between the iterative scheme and the specific exponential integrator under consideration., Comment: in submission; attached supplemental material

We introduce a method, MMD-B-Fair, to learn fair representations of data via kernel two-sample testing. We find neural features of our data where a maximum mean discrepancy (MMD) test cannot distinguish between representations of different sensitive groups, while preserving information about the target attributes. Minimizing the power of an MMD test is more difficult than maximizing it (as done in previous work), because the test threshold's complex behavior cannot be simply ignored. Our method exploits the simple asymptotics of block testing schemes to efficiently find fair representations without requiring complex adversarial optimization or generative modelling schemes widely used by existing work on fair representation learning. We evaluate our approach on various datasets, showing its ability to ``hide'' information about sensitive attributes, and its effectiveness in downstream transfer tasks.

Interfacial two-dimensional electron gas (2DEG), especially the SrTiO3-based ones at the unexpected interface of insulators, have emerged to be a promising candidate for efficient charge-spin current interconversion. In this article, to gain insight into the mechanism of the charge-spin current interconversion at the oxide-based 2DEG, we focused on conducting interfaces between insulating SrTiO3 and two types of aluminium-based amorphous insulators, namely SrTiO3/AlN and SrTiO3/Al2O3, and estimated their charge-spin conversion efficiency, {\theta}_cs. The two types of amorphous insulators were selected to explicitly probe the overlooked contribution of oxygen vacancy to the {\theta}_cs. We proposed a mechanism to explain results of spin-torque ferromagnetic resonance (ST-FMR) measurements and developed an analysis protocol to reliably estimate the {\theta}_cs of the oxide based 2DEG. The resultant {\theta}_cs/t, where t is the thickness of the 2DEG, were estimated to be 0.244 nm-1 and 0.101 nm-1 for the SrTiO3/AlN and SrTiO3/Al2O3, respectively, and they are strikingly comparable to their crystalline counterparts. Furthermore, we also observe a large direct current modulation of resonance linewidth in SrTiO3/AlN samples, confirming its high {\theta}_cs and attesting an oxygen-vacancy-enabled charge-spin conversion. Our findings emphasize the defects' contribution to the charge-spin interconversion, especially in the oxide-based low dimensional systems, and provide a way to create and enhance charge-spin interconversion via defect engineering.