title: Causal Inference for Comprehensive Cohort Studies

abstract: In a comprehensive cohort study of two competing treatments (say, A and B), clinically eligible individuals are first asked to enroll in a randomized trial and, if they refuse, are then asked to enroll in a parallel observational study in which they can choose treatment according to their own preference. We consider estimation of two estimands: (1) comprehensive cohort causal effect -- the difference in mean potential outcomes had all patients in the comprehensive cohort received treatment A vs. treatment B and (2) randomized trial causal effect -- the difference in mean potential outcomes had all patients enrolled in the randomized trial received treatment A vs. treatment B. For each estimand, we consider inference under various sets of unconfoundedness assumptions and construct semiparametric efficient and robust estimators. These estimators depend on nuisance functions, which we estimate, for illustrative purposes, using generalized additive models. Using the theory of sample splitting, we establish the asymptotic properties of our proposed estimators. We also illustrate our methodology using data from the Bypass Angioplasty Revascularization Investigation (BARI) randomized trial and observational registry to evaluate the effect of percutaneous transluminal coronary balloon angioplasty versus coronary artery bypass grafting on 5-year mortality. To evaluate the finite sample performance of our estimators, we use the BARI dataset as the basis of a realistic simulation study.


title: Dark matter search in a Beam-Dump eXperiment (BDX) at Jefferson Lab -- 2018 update to PR12-16-001

abstract: This document complements and completes what was submitted last year to PAC45 as an update to the proposal PR12-16-001 "Dark matter search in a Beam-Dump eXperiment (BDX)" at Jefferson Lab submitted to JLab-PAC44 in 2016. Following the suggestions contained in the PAC45 report, in coordination with the lab, we ran a test to assess the beam-related backgrounds and validate the simulation framework used to design the BDX experiment. Using a common Monte Carlo framework for the test and the proposed experiment, we optimized the selection cuts to maximize the reach considering simultaneously the signal, cosmic-ray background (assessed in Catania test with BDX-Proto) and beam-related backgrounds (irreducible NC and CC neutrino interactions as determined by simulation). Our results confirmed what was presented in the original proposal: with 285 days of a parasitic run at 65 $\mu$A (corresponding to $10^{22}$ EOT) the BDX experiment will lower the exclusion limits in the case of no signal by one to two orders of magnitude in the parameter space of dark-matter coupling versus mass.


title: Bethe phase including proton excitations

abstract: We evaluate the contribution of inelastic intermediate states (such as $p\to N^*$ excitations) to the phase between the one-photon-exchange and the `nuclear' high energy $pp$ scattering amplitudes as $t\to 0$, caused by multiphoton diagrams. It turns out to be rather small - much smaller than to have any influence on the experimental accuracy of the measurements of $\rho$, defined to be the ratio of the real to imaginary parts of the forward `nuclear' amplitude.


title: Accurate and Fast Retrieval for Complex Non-metric Data via Neighborhood Graphs

abstract: We demonstrate that a graph-based search algorithm-relying on the construction of an approximate neighborhood graph-can directly work with challenging non-metric and/or non-symmetric distances without resorting to metric-space mapping and/or distance symmetrization, which, in turn, lead to substantial performance degradation. Although the straightforward metrization and symmetrization is usually ineffective, we find that constructing an index using a modified, e.g., symmetrized, distance can improve performance. This observation paves a way to a new line of research of designing index-specific graph-construction distance functions.


title: Approximate frame representations via iterated operator systems

abstract: It is known that it is a very restrictive condition for a frame $\{f_k\}_{k=1}^\infty$ to have a representation $ \{T^n \varphi\}_{n=0}^\infty$ as the orbit of a bounded operator $T$ under a single generator $\varphi\in\mathcal{H}.$ In this paper we prove that, on the other hand, any frame can be approximated arbitrarily well by a suborbit $\{T^{\alpha(k)} \varphi\}_{k=1}^\infty$ of a bounded operator $T$. An important new aspect is that for certain important classes of frames, e.g., frames consisting of finitely supported vectors in $\ell^2(\mathbb{N}),$ we can be completely explicit about possible choices of the operator $T$ and the powers $\alpha(k),k\in \mathbb{N}.$ A similar approach carried out in $L^2(\mathbb{R})$ leads to an approximation of a frame using suborbits of two bounded operators. The results are illustrated with an application to Gabor frames generated by a compactly supported function. The paper is concluded with an appendix which collects general results about frame representations using multiple orbits of bounded operators.


title: Inverse Probability Weighted Estimators of Vaccine Effects Accommodating Partial Interference and Censoring

abstract: Estimating population-level effects of a vaccine is challenging because there may be interference, i.e., the outcome of one individual may depend on the vaccination status of another individual. Partial interference occurs when individuals can be partitioned into groups such that interference occurs only within groups. In the absence of interference, inverse probability weighted (IPW) estimators are commonly used to draw inference about causal effects of an exposure or treatment. Tchetgen Tchetgen and VanderWeele (2012) proposed a modified IPW estimator for causal effects in the presence of partial interference. Motivated by a cholera vaccine study in Bangladesh, this paper considers an extension of the Tchetgen Tchetgen and VanderWeele IPW estimator to the setting where the outcome is subject to right censoring using inverse probability of censoring weights (IPCW). Censoring weights are estimated using proportional hazards frailty models. The large sample properties of the IPCW estimators are derived, and simulation studies are presented demonstrating the estimators' performance in finite samples. The methods are then used to analyze data from the cholera vaccine study.


title: Sharp nonzero lower bounds for the Schur product theorem

abstract: By a result of Schur [J. Reine Angew. Math. 1911], the entrywise product $M \circ N$ of two positive semidefinite matrices $M,N$ is again positive. Vybiral [Adv. Math. 2020] improved on this by showing the uniform lower bound $M \circ \overline{M} \geq E_n / n$ for all $n \times n$ real or complex correlation matrices $M$, where $E_n$ is the all-ones matrix. This was applied to settle a conjecture of Novak [J. Complexity 1999] and to positive definite functions on groups. Vybiral (in his original preprint) asked if one can obtain similar uniform lower bounds for higher entrywise powers of $M$, or for $M \circ N$ when $N \neq M, \overline{M}$. A natural third question is to obtain a tighter lower bound that need not vanish as $n \to \infty$, i.e. over infinite-dimensional Hilbert spaces. In this note, we affirmatively answer all three questions by extending and refining Vybiral's result to lower-bound $M \circ N$, for arbitrary complex positive semidefinite matrices $M, N$. Specifically: we provide tight lower bounds, improving on Vybiral's bounds. Second, our proof is 'conceptual' (and self-contained), providing a natural interpretation of these improved bounds via tracial Cauchy-Schwarz inequalities. Third, we extend our tight lower bounds to Hilbert-Schmidt operators. As an application, we settle Open Problem 1 of Hinrichs-Krieg-Novak-Vybiral [J. Complexity, in press], which yields improvements in the error bounds in certain tensor product (integration) problems.


title: Overgroups of Levi subgroups I. The case of abelian unipotent radical

abstract: In the present paper we prove sandwich classification for the overgroups of the subsystem subgroup $E(\Delta,R)$ of the Chevalley group $G(\Phi,R)$ for the three types of pair $(\Phi,\Delta)$ (the root system and its subsystem) such that the group $G(\Delta,R)$ is (up to torus) a Levi subgroup of the parabolic subgroup with abelian unipotent radical. Namely we show that for any such an overgroup $H$ there exists a unique pair of ideals $\sigma$ of the ring $R$ such that $E(\Phi,\Delta,R,\sigma)\le H\le N_{G(\Phi,R)}(E(\Phi,\Delta,R,\sigma))$.


title: Pruning Algorithms for Low-Dimensional Non-metric k-NN Search: A Case Study

abstract: We focus on low-dimensional non-metric search, where tree-based approaches permit efficient and accurate retrieval while having short indexing time. These methods rely on space partitioning and require a pruning rule to avoid visiting unpromising parts. We consider two known data-driven approaches to extend these rules to non-metric spaces: TriGen and a piece-wise linear approximation of the pruning rule. We propose and evaluate two adaptations of TriGen to non-symmetric similarities (TriGen does not support non-symmetric distances). We also evaluate a hybrid of TriGen and the piece-wise linear approximation pruning. We find that this hybrid approach is often more effective than either of the pruning rules. We make our software publicly available.


title: Azimuthal variations of oxygen abundance profiles in star-forming regions of disc galaxies in the EAGLE simulations

abstract: The exploration of the spatial distribution of chemical abundances in star-forming regions in galactic discs provides clues to understand the complex interplay of physical processes that regulate the star formation activity and the chemical enrichment across a galaxy. We study the azimuthal variations of the normalized oxygen abundance profiles in the highest numerical resolution run of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) Project at $z=0$. We use young stellar populations to trace the abundances of star-forming regions. Oxygen profiles are estimated along different line of sights from a centrally located observer.The mean azimuthal variation in the EAGLE discs are $\sim 0.12 \pm 0.03$~dex~$R_{\rm eff}^{-1}$ for slopes and $\sim 0.12 \pm 0.03$~dex for the zero points, in agreement with previous works. Metallicity gradients measured along random directions correlate with those determine by averaging over the whole discs although with a large dispersion. We find a slight trend for higher azimuthal variations in the disc components of low star-forming and bulge-dominated galaxies. We also investigate the metallicity profiles of stellar populations with higher and lower levels of enrichment than the average metallicity profiles, and we find that high star-forming region with high metallicity tend to have slightly shallower metallicity slopes compared with the overall metallicity gradient. The simulated azimuthal variations in the EAGLE discs are in global agreement with observations, although the large variety of metallicity gradients would encourage further exploration of the metal mixing in numerical simulations.


title: Global $C^{2,\alpha}$ estimates for the Monge-Amp\`ere equation on polygonal domains in the plane

abstract: We classify global solutions of the Monge-Amp\`ere equation $\det D^2 u=1 $ on the first quadrant in the plane with quadratic boundary data. As an application, we obtain global $C^{2,\alpha}$ estimates for the non-degenerate Monge-Amp\`ere equation in convex polygonal domains in $\mathbb R^2$ provided a globally $C^2$, convex strict subsolution exists.


title: Canonical extensions of locally compact frames

abstract: Canonical extension of finitary ordered structures such as lattices, posets, proximity lattices, etc., is a certain completion which entirely describes the topological dual of the ordered structure and it does so in a purely algebraic and choice-free way. We adapt the general algebraic technique that constructs them to the theory of frames. As a result, we show that every locally compact frame embeds into a completely distributive lattice by a construction which generalises, among others, the canonical extensions for distributive lattices and proximity lattices. This construction also provides a new description of a construction by Marcel Ern\'e. Moreover, canonical extensions of frames enable us to frame-theoretically represent monotone maps with respect to the specialisation order.


title: An Analysis of the Milstein Scheme for SPDEs without a Commutative Noise Condition

abstract: In order to approximate solutions of stochastic partial differential equations (SPDEs) that do not possess commutative noise, one has to simulate the involved iterated stochastic integrals. Recently, two approximation methods for iterated stochastic integrals in infinite dimensions were introduced in C. Leonhard and A. R\"o{\ss}ler: Iterated stochastic integrals in infinite dimensions: approximation and error estimates, Stoch. Partial Differ. Equ. Anal. Comput., 7(2): 209-239 (2019). As a result of this, it is now possible to apply the Milstein scheme by Jentzen and R\"ockner: A Milstein scheme for SPDEs, Found. Comput. Math., 15(2): 313-362 (2015) to equations that need not fulfill the commutativity condition. We prove that the order of convergence of the Milstein scheme can be maintained when combined with one of the two approximation methods for iterated stochastic integrals. However, we also have to consider the computational cost and the corresponding effective order of convergence for a meaningful comparison with other schemes. An analysis of the computational cost shows that, in dependence on the equation, a combination of the Milstein scheme with both of the two methods may be the preferred choice. Further, the Milstein scheme is compared to the exponential Euler scheme and we show for different SPDEs depending on the parameters describing, e.g., the regularity of the equation, which one of the schemes achieves the highest effective order of convergence.


title: Find or Classify? Dual Strategy for Slot-Value Predictions on Multi-Domain Dialog State Tracking

abstract: Dialog state tracking (DST) is a core component in task-oriented dialog systems. Existing approaches for DST mainly fall into one of two categories, namely, ontology-based and ontology-free methods. An ontology-based method selects a value from a candidate-value list for each target slot, while an ontology-free method extracts spans from dialog contexts. Recent work introduced a BERT-based model to strike a balance between the two methods by pre-defining categorical and non-categorical slots. However, it is not clear enough which slots are better handled by either of the two slot types, and the way to use the pre-trained model has not been well investigated. In this paper, we propose a simple yet effective dual-strategy model for DST, by adapting a single BERT-style reading comprehension model to jointly handle both the categorical and non-categorical slots. Our experiments on the MultiWOZ datasets show that our method significantly outperforms the BERT-based counterpart, finding that the key is a deep interaction between the domain-slot and context information. When evaluated on noisy (MultiWOZ 2.0) and cleaner (MultiWOZ 2.1) settings, our method performs competitively and robustly across the two different settings. Our method sets the new state of the art in the noisy setting, while performing more robustly than the best model in the cleaner setting. We also conduct a comprehensive error analysis on the dataset, including the effects of the dual strategy for each slot, to facilitate future research.


title: AFeSe2 (A=Tl, K, Rb, or Cs): Iron-based superconducting analog of the cuprates

abstract: It has long been a challenging task to find compounds with similar crystal and electronic structures as cuprate superconductors with low dimensionality and strong antiferromagnetic fluctuations. The parent compounds of cuprate superconductors are Mott insulators with strong in-plane antiferromagnetic exchange interactions between Cu moments. Here we show, based on first-principles density functional calculations, that AFeSe2 (A=Tl, K, Rb, or Cs) exhibit many of the physical properties common to the cuprate parent compounds: (1) the FeSe2 layer in AFeSe2 is similar in crystalline and electronic structures to the CuO2 plane in cuprates, although Se atoms are not coplanar to the square Fe-lattice; (2) they are antiferromagnetic insulators, but with relatively small charge excitation gaps; (3) their ground states are Neel antiferromagnetic ordered, similar as in cuprates; and (4) the antiferromagnetic exchange interactions between Fe moments are larger than in other iron-based superconducting materials, but comparable to those in cuprates. Like cuprates, these compounds may become high-Tc superconductors upon doping of charge carriers either by chemical substitution or intercalation or by liquid or solid gating.


title: Non-singular Bounce Cosmology from Lagrange Multiplier $F(R)$ Gravity

abstract: In this work we study non-singular{ bounce cosmology} in the context of the Lagrange multiplier generalized $F(R)$ gravity theory of gravity. We specify our study by using a specific variant form of the well known matter{ bounce cosmology}, with scale factor $a(t) = \left(a_0t^2 + 1 \right)^n$, and we demonstrate that for $n < 1/2$, the primordial curvature perturbations are generated deeply in the contraction era. Particularly, we show explicitly that the perturbation modes exit the horizon at a large negative time during the contraction era, which in turn makes the ``low-curvature'' regime, the era for which the calculations of observational indices related to the primordial power spectrum can be considered reliable. Using the reconstruction techniques for the Lagrange multiplier $F(R)$ gravity, we construct the form of effective $F(R)$ gravity that can realize such a cosmological evolution, and we determine the power spectrum of the primordial curvature perturbations. Accordingly, we calculate the spectral index of the primordial curvature perturbations and the tensor-to-scalar ratio, and we confront these with the latest observational data. We also address the issue of stability of the primordial metric perturbations, and to this end, we determine the form of $F(R)$ which realizes the non-singular cosmology for the whole range of cosmic time $-\infty < t < \infty$, by solving the Friedmann equations without the ``low-curvature'' approximation. This study is performed numerically though, due to the high complexity of the resulting differential equations. We also investigate the energy conditions in the present context. The phenomenology of the non-singular bounce is also studied in the context of a standard $F(R)$ gravity.


title: Some results on higher eigenvalue optimization

abstract: In this paper we obtain several results concerning the optimization of higher Steklov eigenvalues both in two and higher dimensional cases. We first show that the normalized (by boundary length) $k$-th Steklov eigenvalue on the disk is not maximized for a smooth metric on the disk for $k\geq 3$. For $k=1$ the classical result of [W] shows that $\sigma_1$ is maximized by the standard metric on the round disk. For $k=2$ it was shown [GP1] that $\sigma_2$ is not maximized for a smooth metric. We also prove a local rigidity result for the critical catenoid and the critical M\"obius band as free boundary minimal surfaces in a ball under $C^2$ deformations. We next show that the first $k$ Steklov eigenvalues are continuous under certain degenerations of Riemannian manifolds in any dimension. Finally we show that for $k\geq 2$ the supremum of the $k$-th Steklov eigenvalue on the annulus over all metrics is strictly larger that that over $S^1$-invariant metrics. We prove this same result for metrics on the M\"obius band.


title: Multi-Source Domain Adaptation and Semi-Supervised Domain Adaptation with Focus on Visual Domain Adaptation Challenge 2019

abstract: This notebook paper presents an overview and comparative analysis of our systems designed for the following two tasks in Visual Domain Adaptation Challenge (VisDA-2019): multi-source domain adaptation and semi-supervised domain adaptation. Multi-Source Domain Adaptation: We investigate both pixel-level and feature-level adaptation for multi-source domain adaptation task, i.e., directly hallucinating labeled target sample via CycleGAN and learning domain-invariant feature representations through self-learning. Moreover, the mechanism of fusing features from different backbones is further studied to facilitate the learning of domain-invariant classifiers. Source code and pre-trained models are available at \url{https://github.com/Panda-Peter/visda2019-multisource}. Semi-Supervised Domain Adaptation: For this task, we adopt a standard self-learning framework to construct a classifier based on the labeled source and target data, and generate the pseudo labels for unlabeled target data. These target data with pseudo labels are then exploited to re-training the classifier in a following iteration. Furthermore, a prototype-based classification module is additionally utilized to strengthen the predictions. Source code and pre-trained models are available at \url{https://github.com/Panda-Peter/visda2019-semisupervised}.


title: Dilation theory in finite dimensions and matrix convexity

abstract: We establish a finite-dimensional version of the Arveson-Stinespring dilation theorem for unital completely positive maps on operator systems. This result can be seen as a general principle to deduce finite-dimensional dilation theorems from their classical infinite-dimensional counterparts. In addition to providing unified proofs of known finite-dimensional dilation theorems, we establish finite-dimensional versions of Agler's theorem on rational dilation on an annulus, of Berger's dilation theorem for operators of numerical radius at most $1$, and of the Putinar-Sandberg numerical range dilation theorem. As a key tool, we prove versions of Carath\'{e}odory's and of Minkowski's theorem for matrix convex sets.


title: Gravitational Lensing of the Cosmic Neutrino Background

abstract: We study gravitational lensing of the cosmic neutrino background. This signal is undetectable for the foreseeable future, but there is a rich trove of information available. At least some of the neutrinos from the early universe will be non-relativistic today, with a closer surface of last scattering (compared to the cosmic microwave background) and with larger angles of deflection. Lensing of massive neutrinos is strongly chromatic: both the amplitude of lensing and the cosmic time at which the potential is traversed depend on neutrino momentum, in principle giving access to our entire causal volume, not restricted to the light cone. As a concrete example, we focus on the case where the cosmic neutrino background would be strongly lensed when passing through halos of galaxy clusters and galaxies. We calculate the Einstein radius for cosmic neutrinos and investigate the impact of neutrino mass.


title: Quantum encryption with certified deletion

abstract: Given a ciphertext, is it possible to prove the deletion of the underlying plaintext? Since classical ciphertexts can be copied, clearly such a feat is impossible using classical information alone. In stark contrast to this, we show that quantum encodings enable certified deletion. More precisely, we show that it is possible to encrypt classical data into a quantum ciphertext such that the recipient of the ciphertext can produce a classical string which proves to the originator that the recipient has relinquished any chance of recovering the plaintext should the decryption key be revealed. Our scheme is feasible with current quantum technology: the honest parties only require quantum devices for single-qubit preparation and measurements; the scheme is also robust against noise in these devices. Furthermore, we provide an analysis that is suitable in the finite-key regime.


title: TorchBeast: A PyTorch Platform for Distributed RL

abstract: TorchBeast is a platform for reinforcement learning (RL) research in PyTorch. It implements a version of the popular IMPALA algorithm for fast, asynchronous, parallel training of RL agents. Additionally, TorchBeast has simplicity as an explicit design goal: We provide both a pure-Python implementation ("MonoBeast") as well as a multi-machine high-performance version ("PolyBeast"). In the latter, parts of the implementation are written in C++, but all parts pertaining to machine learning are kept in simple Python using PyTorch, with the environments provided using the OpenAI Gym interface. This enables researchers to conduct scalable RL research using TorchBeast without any programming knowledge beyond Python and PyTorch. In this paper, we describe the TorchBeast design principles and implementation and demonstrate that it performs on-par with IMPALA on Atari. TorchBeast is released as an open-source package under the Apache 2.0 license and is available at \url{https://github.com/facebookresearch/torchbeast}.


title: Differentially private anonymized histograms

abstract: For a dataset of label-count pairs, an anonymized histogram is the multiset of counts. Anonymized histograms appear in various potentially sensitive contexts such as password-frequency lists, degree distribution in social networks, and estimation of symmetric properties of discrete distributions. Motivated by these applications, we propose the first differentially private mechanism to release anonymized histograms that achieves near-optimal privacy utility trade-off both in terms of number of items and the privacy parameter. Further, if the underlying histogram is given in a compact format, the proposed algorithm runs in time sub-linear in the number of items. For anonymized histograms generated from unknown discrete distributions, we show that the released histogram can be directly used for estimating symmetric properties of the underlying distribution.


title: Nonlinear nonuniform $\mathcal{PT}$-symmetric Bragg grating structures

abstract: We explore the consequences of incorporating parity and time reversal ($\mathcal{PT}$) symmetries on the dynamics of nonreciprocal light propagation exhibited by a class of nonuniform periodic structures known as chirped $\mathcal{PT}$-symmetric fiber Bragg gratings (FBGs). The interplay among various grating parameters such as chirping, detuning, nonlinearities, and gain/loss gives rise to unique bi- and multi-stable states in the unbroken as well as broken $\mathcal{PT}$-symmetric regimes. The role of chirping on the steering dynamics of the hysteresis curve is influenced by the type of nonlinearities and the nature of detuning parameter. Also, incident directions of the input light robustly impact the steering dynamics of bistable and multistable states both in the unbroken and broken $\mathcal{PT}$-symmetric regimes. When the light launching direction is reversed, critical stable states are found to occur at very low intensities which opens up a new avenue for an additional way of controlling light with light. We also analyze the phenomenon of unidirectional wave transport and the reflective bi- and multi-stable characteristics at the so-called $\mathcal{PT}$-symmetry breaking point.


title: Failure analysis and reliability evaluation of modulation techniques for neutral point clamped inverters-A usage model approach

abstract: Up to now, many modulation techniques have been proposed for neutral point clamped (NPC) inverters. In this paper, for the first time, a general methodology is applied to calculate and compare the failure analysis and reliability of NPC inverter with most commonly used control strategies. Also, the mean time to failure (MTTF) of NPC inverter is derived for different control strategies. It is demonstrated that the key feature of control strategies in determining the reliability of inverter is their loss distribution among the switches. The failure rate of components that is relevant to this study and junction temperature calculation is developed, then conduction losses and switching losses of switches for different control strategies are calculated. Finally, the most reliable control strategy is identified. Experimental results obtained have promptly justified the theoretical analysis and outlined procedure.


title: Counterexamples on the monotonicity of delay optimal strategies for energy harvesting transmitters

abstract: We consider cross-layer design of delay optimal transmission strategies for energy harvesting transmitters where the data and energy arrival processes are stochastic. Using Markov decision theory, we show that the value function is weakly increasing in the queue state and weakly decreasing in the battery state. It is natural to expect that the delay optimal policy should be weakly increasing in the queue and battery states. We show via counterexamples that this is not the case. In fact, we show that for some sample scenarios the delay optimal policy may perform 5-13% better than the best monotone policy.


title: SUGAR-R: Robust Online Restoration Platform for SCADA-Absent Grid

abstract: Secure and fast grid restoration from a collapsed state is increasingly critical as blackouts are becoming more prevalent around the globe. Generally, the restoration of grid during a blackout is achieved with the help of Supervisory Control and Data Acquisition System (SCADA) based central control; however, with the threat of cyber-blackouts, this presumption of an available and secure SCADA system is not valid. This is also true for grids in developing countries as well as for many distribution and micro grids that lack SCADA. In this paper, we introduce an online framework for localized grid restoration that validates and updates a pre-defined crank path in real-time based on the vital grid states of voltages, currents and frequency. The proposed framework maintains an online network topology of the localized grid that can continuously sample local measurements and update the grid model, thereby circumventing SCADA based central control. In the results section we demonstrate the efficacy of this framework for blackstart by ensuring a feasible crank path with voltage and frequency states within bounds, while further assisting in synchronization of two disconnected sub-grids during the re-energization process using a distributed framework.


title: A Step by Step Mathematical Derivation and Tutorial on Kalman Filters

abstract: We present a step by step mathematical derivation of the Kalman filter using two different approaches. First, we consider the orthogonal projection method by means of vector-space optimization. Second, we derive the Kalman filter using Bayesian optimal filtering. We provide detailed proofs for both methods and each equation is expanded in detail.


title: Efficient implementation of adaptive order reconstructions

abstract: Including polynomials with small degree and stencil when designing very high order reconstructions is surely beneficial for their non oscillatory properties, but may bring loss of accuracy on smooth data unless special care is exerted. In this paper we address this issue with a new Central WENOZ (CWENOZ) approach, in which the reconstruction polynomial is computed from a single set of non linear weights, but the linear weights of the polynomials with very low degree (compared to the final desired accuracy) are infinitesimal with respect to the grid size. After proving general results that guide the choice of the CWENOZ parameters, we study a concrete example of a reconstruction that blends polynomials of degree six, four and two, mimicking already published Adaptive Order WENO reconstructions [Balsara, Garain and Shu, J.Comput.Phys., 2016], [Arbogast, Huang and Zhao, SIAM J.Numer.Anal., 2018]. The novel reconstruction yields similar accuracy and oscillations with respect to the previous ones, but saves up to 20% computational time since it does not rely on a hierarchic approach and thus does not compute multiple sets of nonlinear weights in each cell.


title: When Does Self-supervision Improve Few-shot Learning?

abstract: We investigate the role of self-supervised learning (SSL) in the context of few-shot learning. Although recent research has shown the benefits of SSL on large unlabeled datasets, its utility on small datasets is relatively unexplored. We find that SSL reduces the relative error rate of few-shot meta-learners by 4%-27%, even when the datasets are small and only utilizing images within the datasets. The improvements are greater when the training set is smaller or the task is more challenging. Although the benefits of SSL may increase with larger training sets, we observe that SSL can hurt the performance when the distributions of images used for meta-learning and SSL are different. We conduct a systematic study by varying the degree of domain shift and analyzing the performance of several meta-learners on a multitude of domains. Based on this analysis we present a technique that automatically selects images for SSL from a large, generic pool of unlabeled images for a given dataset that provides further improvements.


title: Deep Network Classification by Scattering and Homotopy Dictionary Learning

abstract: We introduce a sparse scattering deep convolutional neural network, which provides a simple model to analyze properties of deep representation learning for classification. Learning a single dictionary matrix with a classifier yields a higher classification accuracy than AlexNet over the ImageNet 2012 dataset. The network first applies a scattering transform that linearizes variabilities due to geometric transformations such as translations and small deformations. A sparse $\ell^1$ dictionary coding reduces intra-class variability while preserving class separation through projections over unions of linear spaces. It is implemented in a deep convolutional network with a homotopy algorithm having an exponential convergence. A convergence proof is given in a general framework that includes ALISTA. Classification results are analyzed on ImageNet.


title: Two-component Yang-Baxter maps and star-triangle relations

abstract: It is shown how Yang-Baxter maps may be directly obtained from classical counterparts of the star-triangle relations and quantum Yang-Baxter equations. This is based on reinterpreting the latter equation and its solutions which are given in terms of special functions, as a set-theoretical form of the Yang-Baxter equation whose solutions are given by quadrirational Yang-Baxter maps. The Yang-Baxter maps obtained through this approach are found to satisfy two different types of Yang-Baxter equations, one that is the usual equation involving a single map, and another equation that involves a pair of maps, which is a case of what is also known as an entwining Yang-Baxter equation. Apart from the elliptic case, each of these Yang-Baxter maps are quadrirational, but only maps that solve the former type of Yang-Baxter equation are reversible. The Yang-Baxter maps are expressed in terms of two-component variables, and two-component parameters, and have a natural QRT-like composition of separate maps for each component. Through this approach, sixteen different Yang-Baxter maps are derived from known solutions of the classical star-triangle relations.

categories: math-ph cond-mat.stat-mech math.MP math.QA nlin.SI


title: Recovering function fields from their integral $\ell$-adic cohomology with the Galois action

abstract: In this note, we consider function fields of higher-dimensional algebraic varieties defined over non-local fields, and show how the Galois action on the cohomology such function fields can be used to parameterize their divisorial valuations. By applying a recent theorem of Pop, we observe that, in dimension $\geq 3$, this information is enough to completely determine the function field and the base-field in question.


title: Timely Distributed Computation with Stragglers

abstract: We consider a status update system in which the update packets need to be processed to extract the embedded useful information. The source node sends the acquired information to a computation unit (CU) which consists of a master node and $n$ worker nodes. The master node distributes the received computation task to the worker nodes. Upon computation, the master node aggregates the results and sends them back to the source node to keep it \emph{updated}. We investigate the age performance of uncoded and coded (repetition coded, MDS coded, and multi-message MDS (MM-MDS) coded) schemes in the presence of stragglers under i.i.d.~exponential transmission delays and i.i.d~shifted exponential computation times. We show that asymptotically MM-MDS coded scheme outperforms the other schemes. Furthermore, we characterize the optimal codes such that the average age is minimized.


title: Constructing spherically symmetric Einstein-Dirac systems with multiple spinors: Ansatz, wormholes and other analytical solutions

abstract: In this paper we present a detailed calculation of an Ansatz that allows to obtain spherically symmetric Einstein-Dirac configurations in $d$-dimensions. We show that this is possible by combining $2^{\lfloor \frac{d-2}{2} \rfloor}$ Dirac fields, making use of the properties of the angular dependence of the spinors in a spherical background. By applying this Ansatz, we investigate some simple analytical solutions. One of them is a regular wormhole supported by the Dirac fields. Other solutions include a pathological black hole and a naked singularity. We analyze the domain of existence and properties of all these solutions.


title: Strong Lensing Time Delay Constraints on Dark Energy: a Forecast

abstract: Measurements of time delays between multiple quasar images produced by strong lensing are reaching a sensitivity that makes them a promising cosmological probe. Future surveys will provide significantly more measurements, reaching unprecedented depth in redshift, making strong lensing time delay (SLTD) observations competitive with other background probes. We forecast constraints on the nature of dark energy from upcoming SLTD surveys, simulating future catalogues with different numbers of lenses distributed up to redshift $z\sim 1$ and focusing on cosmological parameters such as the Hubble constant $H_0$ and parametrisations of the dark energy equation of state. We also explore the impact of our ability to precisely model the lens mass profile and its environment, on the forecasted constraints. We find that in the most optimistic cases, SLTD will constrain $H_0$ at the level of $\sim 0.1\%$, while the CPL equation of state parameters, $w_0$ and $w_a$, can be determined with errors $\sigma_{w_0}\sim 0.05$ and $\sigma_{w_a}\sim 0.3$, respectively. Furthermore, we investigate the bias introduced when a wrong cosmological model is assumed for the analysis. We find that the value of $H_0$ could be biased up to $10 \sigma$, assuming a perfect knowledge of the lens profile, when a $\Lambda$CDM model is used to analyse data that really belong to a $w$CDM cosmology with $w=-0.9$. Based on these findings, we identify a consistency check of the assumed cosmological model in future SLTD surveys, by splitting the dataset in several redshift bins. Depending on the characteristics of the survey, this could provide a smoking gun for dark energy.


title: Influence of exact Lorentz-violating mechanism on the critical exponents for massive $q$-deformed $\lambda\phi^{4}$ scalar field theory

abstract: We probe the influence of Lorentz-violating mechanism, treated exactly, on the radiative quantum corrections to critical exponents for massive $q$-deformed O($N$) $\lambda\phi^{4}$ scalar field theories. We attain that task by employing three distinct and independent field-theoretic renormalization group methods. Firstly we compute the critical exponents up to the finite next-to-leading order for then generalizing the results for any loop level. We show that the $q$-deformed critical exponents are insensible to the Lorentz symmetry breaking mechanism thus agreeing with the universality hypothesis.


title: Object-centric Forward Modeling for Model Predictive Control

abstract: We present an approach to learn an object-centric forward model, and show that this allows us to plan for sequences of actions to achieve distant desired goals. We propose to model a scene as a collection of objects, each with an explicit spatial location and implicit visual feature, and learn to model the effects of actions using random interaction data. Our model allows capturing the robot-object and object-object interactions, and leads to more sample-efficient and accurate predictions. We show that this learned model can be leveraged to search for action sequences that lead to desired goal configurations, and that in conjunction with a learned correction module, this allows for robust closed loop execution. We present experiments both in simulation and the real world, and show that our approach improves over alternate implicit or pixel-space forward models. Please see our project page (https://judyye.github.io/ocmpc/) for result videos.


title: Induced pairing of fermionic impurities in a one-dimensional strongly correlated Bose gas

abstract: We investigate numerically the problem of few (one, two) noninteracting spin$-1/2$ fermions in a shallow harmonic trap coupled via contact repulsive interactions to a uniform one-dimensional bath of lattice bosons, described by the Bose-Hubbard model. Through extensive density-matrix renormalization group calculations, we extract the binding energy and the effective mass of quasiparticles, including dressed impurities (polarons) and their two-body bound states (bipolarons), emerging from the effective non-local Casimir interaction between the impurities. We show that the mixture exhibits rather different pairing behaviors depending on the singlet \textit{vs}.~triplet spin state configurations of the two fermions. For opposite spin states, bipolarons are found for any finite value of the impurity-bath coupling. In particular, in the strong coupling regime their binding energy reduces to that of a single polaron, provided the boson-boson repulsion is not too weak. For equal spin states, we show that bipolarons emerge only beyond a critical strength of the Bose-Fermi interaction and their effective mass grows rapidly approaching the strong coupling regime.


title: Observing the tail of reionization: neutral islands in the $z=5.5$ Lyman-$\alpha$ Forest

abstract: Previous studies have noted difficulties in modeling the highest opacities of the $z>5.5$ Ly$\alpha$ forest, epitomized by the extreme Ly$\alpha$ trough observed towards quasar ULAS J0148+0600. One possibility is that the most opaque regions at these redshifts contain significant amounts of neutral hydrogen. This explanation, which abandons the common assumption that reionization ended before $z=6$, also reconciles evidence from independent observations of a significantly neutral Universe at $z=7.5$. Here we explore a model in which the neutral fraction is still $\approx 10\%$ at $z=5.5$. We confirm that this model can account for the observed scatter in Ly$\alpha$ forest opacities, as well as the observed Ly$\beta$ transmission in the J0148 trough. We contrast the model with a competing "earlier" reionization scenario characterized by a short mean free path and large fluctuations in the post-reionization ionizing background. We consider Ly$\alpha$ and Ly$\beta$ effective optical depths, their correlations, trough size distributions, dark pixel fractions, the IGM thermal history, and spatial distributions of Lyman-$\alpha$ emitters around forest sight lines. We find that the models are broadly similar in almost all of these statistics, suggesting that it may be difficult to distinguish between them definitively. We argue that improved constraints on the mean free path and the thermal history at $z>5$ could go a long way towards diagnosing the origin of the $z>5.5$ opacity fluctuations.


title: A versatile combinatorial approach of studying products of long cycles in symmetric groups

abstract: In symmetric groups, studies of permutation factorizations or triples of permutations satisfying certain conditions have a long history. One particular interesting case is when two of the involved permutations are long cycles, for which many surprisingly simple formulas have been obtained. Here we combinatorially enumerate the pairs of long cycles whose product has a given cycle-type and separates certain elements, extending several lines of studies, and we obtain general quantitative relations. As consequences, in a unified way, we recover a number of results expecting simple combinatorial proofs, including results of Boccara (1980), Zagier (1995), Stanley (2011), F\'{e}ray and Vassilieva (2012), as well as Hultman (2014). We obtain a number of new results as well. In particular, for the first time, given a partition of a set, we obtain an explicit formula for the number of pairs of long cycles on the set such that the product of the long cycles does not mix the elements from distinct blocks of the partition and has an independently prescribed number of cycles for each block of elements. As applications, we obtain new explicit formulas concerning factorizations of any even permutation into long cycles and the first nontrivial explicit formula for computing strong separation probabilities solving an open problem of Stanley (2010).


title: Bar Code and Janet-like division

abstract: Bar Codes are combinatorial objects encoding many properties of monomial ideals. In this paper we employ these objects to study Janet-like divisions. Given a finite set of terms U, from its Bar Code we can compute the Janet-like nonmultiplicative power of its elements and detect completeness of the set. Some observation on the computation of Janet-like bases conclude the work.


title: Fixed point results for a new metric space

abstract: In this paper, we introduce the neutrosophic contractive and neutrosophic mapping. We establish some results on fixed points of a neutrosophic mapping.


title: Data assimilation for a quasi-geostrophic model with circulation-preserving stochastic transport noise

abstract: This paper contains the latest installment of the authors' project on developing ensemble based data assimilation methodology for high dimensional fluid dynamics models. The algorithm presented here is a particle filter that combines model reduction, tempering, jittering, and nudging. The methodology is tested on a two-layer quasi-geostrophic model for a $\beta$-plane channel flow with $O(10^6)$ degrees of freedom out of which only a minute fraction are noisily observed. The model is reduced by following the stochastic variational approach for geophysical fluid dynamics introduced in Holm (Proc Roy Soc A, 2015) as a framework for deriving stochastic parametrisations for unresolved scales. The reduction is substantial: the computations are done only for $O(10^4)$ degrees of freedom. We introduce a stochastic time-stepping scheme for the two-layer model and prove its consistency in time. Then, we analyze the effect of the different procedures (tempering combined with jittering and nudging) on the performance of the data assimilation procedure using the reduced model, as well as how the dimension of the observational data (the number of "weather stations") and the data assimilation step affect the accuracy and uncertainty of the results.


title: Active-Code Replacement in the OODIDA Data Analytics Platform

abstract: OODIDA (On-board/Off-board Distributed Data Analytics) is a platform for distributing and executing concurrent data analytics tasks. It targets fleets of reference vehicles in the automotive industry and has a particular focus on rapid prototyping. Its underlying message-passing infrastructure has been implemented in Erlang/OTP. External Python applications perform data analytics tasks. Most work is performed by clients (on-board). A central cloud server performs supplementary tasks (off-board). OODIDA can be automatically packaged and deployed, which necessitates restarting parts of the system, or all of it. This is potentially disruptive. To address this issue, we added the ability to execute user-defined Python modules on clients as well as the server. These modules can be replaced without restarting any part of the system and they can even be replaced between iterations of an ongoing assignment. This facilitates use cases such as iterative A/B testing of machine learning algorithms or modifying experimental algorithms on-the-fly.


title: On the strong instability of the multi-layer Hele-Shaw flows

abstract: We study the effects of some injection policies used in oil recovery process. The Saffman-Taylor instability occurs when a less viscous fluid is displacing a more viscous one, in a rectangular Hele-Shaw cell. The injection of $N$ successive intermediate phases with constant viscosities (the multi-layer Hele-Shaw model) was studied in some recent papers, where a minimization of the Saffman-Taylor instability was obtained for large enough $N$. However, in this paper we get a particular eigenfunction of the linear stability system which leads to eigenvalues which become infinite for large wave numbers. We obtain a strong instability of the multi-layer Hele-Shaw displacement, even if $N$ is very large.


title: Dynamic Brain Functional Networks Guided By Anatomical Knowledge

abstract: Recently, the potential of dynamic brain networks as a neuroimaging biomarkers for mental illnesses is being increasingly recognized. However, there are several unmet challenges in developing such biomarkers, including the need for methods to model rapidly changing network states. In one of the first such efforts, we develop a novel approach for computing dynamic brain functional connectivity (FC), that is guided by brain structural connectivity (SC) computed from diffusion tensor imaging (DTI) data. The proposed approach involving dynamic Gaussian graphical models decomposes the time course into non-overlapping state phases determined by change points, each having a distinct network. We develop an optimization algorithm to implement the method such that the estimation of both the change points and the state-phase specific networks are fully data driven and unsupervised, and guided by SC information. The approach is scalable to large dimensions and extensive simulations illustrate its clear advantages over existing methods in terms of network estimation accuracy and detecting dynamic network changes. An application of the method to a posttraumatic stress disorder (PTSD) study reveals important dynamic resting state connections in regions of the brain previously implicated in PTSD. We also illustrate that the dynamic networks computed under the proposed method are able to better predict psychological resilience among trauma exposed individuals compared to existing dynamic and stationary connectivity approaches, which highlights its potential as a neuroimaging biomarker.


title: Stack Sorting with Increasing and Decreasing Stacks

abstract: We introduce a sorting machine consisting of $k+1$ stacks in series: the first $k$ stacks can only contain elements in decreasing order from top to bottom, while the last one has the opposite restriction. This device generalizes \cite{SM}, which studies the case $k=1$. Here we show that, for $k=2$, the set of sortable permutations is a class with infinite basis, by explicitly finding an antichain of minimal nonsortable permutations. This construction can easily be adapted to each $k \ge 3$. Next we describe an optimal sorting algorithm, again for the case $k=2$. We then analyze two types of left-greedy sorting procedures, obtaining complete results in one case and only some partial results in the other one. We close the paper by discussing a few open questions.


title: Graph-based Spatial-temporal Feature Learning for Neuromorphic Vision Sensing

abstract: Neuromorphic vision sensing (NVS)\ devices represent visual information as sequences of asynchronous discrete events (a.k.a., "spikes") in response to changes in scene reflectance. Unlike conventional active pixel sensing (APS), NVS allows for significantly higher event sampling rates at substantially increased energy efficiency and robustness to illumination changes. However, feature representation for NVS is far behind its APS-based counterparts, resulting in lower performance in high-level computer vision tasks. To fully utilize its sparse and asynchronous nature, we propose a compact graph representation for NVS, which allows for end-to-end learning with graph convolution neural networks. We couple this with a novel end-to-end feature learning framework that accommodates both appearance-based and motion-based tasks. The core of our framework comprises a spatial feature learning module, which utilizes residual-graph convolutional neural networks (RG-CNN), for end-to-end learning of appearance-based features directly from graphs. We extend this with our proposed Graph2Grid block and temporal feature learning module for efficiently modelling temporal dependencies over multiple graphs and a long temporal extent. We show how our framework can be configured for object classification, action recognition and action similarity labeling. Importantly, our approach preserves the spatial and temporal coherence of spike events, while requiring less computation and memory. The experimental validation shows that our proposed framework outperforms all recent methods on standard datasets. Finally, to address the absence of large real-world NVS datasets for complex recognition tasks, we introduce, evaluate and make available the American Sign Language letters (ASL-DVS), as well as human action dataset (UCF101-DVS, HMDB51-DVS and ASLAN-DVS).


title: The development of an optical design tool for atmospheric dispersion correction

abstract: In ground based astronomical observations, atmospheric dispersion shifts the image of the object at different wavelengths due to the wavelength-dependent index of refraction of the atmosphere. Thus, using an Atmospheric Dispersion Corrector (ADC) is mandatory in order to avoid any wavelength dependent losses. Typical ADC configurations, for high resolution astronomical instruments, are two counter-rotating prisms, a set of, at least, four prisms paired together. With the arrival of large telescopes with higher angular magnification, and spectrographs with higher resolution, the requirements on the dispersion correction are becoming more critical due to the impact on the produced science (e.g. radial velocity precision). We developed an ADC optical design tool in order to select the best set of glasses in terms of residuals, transmission, resulting image quality, Fresnel losses, taking into account the required spectral range and typical atmospheric conditions where the ADC will be working. A demonstration of the capabilities of the tool is presented with the analysis of the impact of different melt data, the effect of different glass Sellmeier coefficients between catalog and measured ones, that can create a difference in the residuals above few tens of milli-arcseconds (mas). The tool allows the investigation of critical steps on the ADC design phase and speeds up the glass selection process critical for the harder requirements of the future instruments/telescopes.


title: FedMD: Heterogenous Federated Learning via Model Distillation

abstract: Federated learning enables the creation of a powerful centralized model without compromising data privacy of multiple participants. While successful, it does not incorporate the case where each participant independently designs its own model. Due to intellectual property concerns and heterogeneous nature of tasks and data, this is a widespread requirement in applications of federated learning to areas such as health care and AI as a service. In this work, we use transfer learning and knowledge distillation to develop a universal framework that enables federated learning when each agent owns not only their private data, but also uniquely designed models. We test our framework on the MNIST/FEMNIST dataset and the CIFAR10/CIFAR100 dataset and observe fast improvement across all participating models. With 10 distinct participants, the final test accuracy of each model on average receives a 20% gain on top of what's possible without collaboration and is only a few percent lower than the performance each model would have obtained if all private datasets were pooled and made directly available for all participants.


title: Resonantly enhanced polariton wave mixing and Floquet parametric instability

abstract: We introduce a new theoretical approach for analyzing pump and probe experiments in non-linear acousto-optic systems. In our approach, the effect of coherently pumped polaritons is modeled as providing time-periodic modulation of the system parameters. Within this framework, propagation of the probe pulse is described by the Floquet version of Maxwell's equations and leads to such phenomena as frequency mixing and resonant parametric production of polariton pairs. We analyze light reflection from a slab of insulating material with a strongly excited phonon-polariton mode and obtain analytic expressions for the frequency-dependent reflection coefficient for the probe pulse. Our results are in agreement with recent experiments by Cartella et al. which demonstrated light amplification in resonantly excited SiC insulator. We show that, beyond a critical pumping strength, such systems should exhibit Floquet parametric instability, which corresponds to resonant scattering of the pump polaritons into pairs of finite momentum polaritons. We find that the parametric instability should be achievable in SiC using current experimental techniques and discuss its signatures, including the non-analytic frequency dependence of the reflection coefficient and the probe pulse afterglow. We discuss possible applications of the parametric instability phenomenon and suggest that similar types of instabilities can be present in other photoexcited non-linear systems.


title: Visualizing dissipative charge carrier dynamics at the nanoscale with superconducting charge qubit microscopy

abstract: The investigation of novel electronic phases in low-dimensional quantum materials demands for the concurrent development of new measurement techniques that combine surface sensitivity with high spatial resolution and high measurement accuracy. We propose a new quantum sensing imaging modality based on superconducting charge qubits to study dissipative charge carrier dynamics with nanometer spatial and high temporal resolution. Using analytical and numerical calculations we show that superconducting charge qubit microscopy (SCQM) has the potential to resolve temperature and resistivity changes in a sample as small as $\Delta T\leq0.1\;$mK and $\Delta\rho\leq1\cdot10^{4} \,\Omega\cdot$cm, respectively. Among other applications, SCQM will be especially suited to study the microscopic mechanisms underlying interaction driven quantum phase transitions, to investigate the boundary modes found in novel topological insulators and, in a broader context, to visualize the dissiaptive charge carrier dynamics occurring in mesoscopic and nanoscale devices.


title: Long-lived interacting phases of matter protected by multiple time-translation symmetries in quasiperiodically-driven systems

abstract: We show how a large family of interacting nonequilibrium phases of matter can arise from the presence of multiple time-translation symmetries, which occur by quasiperiodically driving an isolated quantum many-body system with two or more incommensurate frequencies. These phases are fundamentally different from those realizable in time-independent or periodically-driven (Floquet) settings. Focusing on high-frequency drives with smooth time-dependence, we rigorously establish general conditions for which these phases are stable in a parametrically long-lived `preheating' regime. We develop a formalism to analyze the effect of the multiple time-translation symmetries on the dynamics of the system, which we use to classify and construct explicit examples of the emergent phases. In particular, we discuss time quasi-crystals which spontaneously break the time-translation symmetries, as well as time-translation symmetry protected topological phases.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.stat-mech quant-ph


title: MHD Turbulence

abstract: We review the current status of research in MHD turbulence theory and numerical experiments and their applications to astrophysics and solar science. We introduce general tools for studying turbulence, basic turbulence models, MHD equations and their wave modes. Subsequently, we cover the theories and numerics of Alfvenic turbulence, imbalanced turbulence, small-scale dynamos and models and numerics for supersonic MHD turbulence.

categories: astro-ph.GA astro-ph.SR physics.flu-dyn physics.plasm-ph


title: An effective field theory approach to quarkonium at small transverse momentum

abstract: In this work we apply effective field theory (EFT) to observables in quarkonium production and decay that are sensitive to soft gluon radiation, in particular measurements that are sensitive to small transverse momentum. Within the EFT framework we study $\chi_Q$ decay to light quarks followed by the fragmentation of those quarks to light hadrons. We derive a factorization theorem that involves transverse momentum distribution (TMD) fragmentation functions and new quarkonium TMD shape functions. We derive renormalization group equations, both in rapidity and virtuality, which are used to evolve the different terms in the factorization theorem to resum large logarithms. This theoretical framework will provide a systematic treatment of quarkonium production and decay processes in TMD sensitive measurements.


title: Stellar Population Astrophysics (SPA) with the TNG: Identification of a Sulphur line at lambda(air) = 1063.600nm in GIANO-B stellar spectra

abstract: Context. In the advent of new infrared, high-resolution spectrometers, accurate and precise atomic data in the infrared is urgently needed. Identifications, wavelengths, strengths, broadening and hyper-fine splitting parameters of stellar lines in the near-IR are in many cases not accurate enough to model observed spectra, and in other cases even non existing. Some stellar features are unidentified. Aims. The aim with this work is to identify a spectral feature at lambda(vac) = 1063.891 nm or lambda(air) = 1063.600 nm seen in spectra of stars of different spectral types, observed with the GIANO-B spectrometer. Methods. Searching for spectral lines to match the unidentified feature in linelists from standard atomic databases was not successful. However, by investigating the original, published laboratory data we were able to identify the feature and solve the problem. To confirm its identification, we model the presumed stellar line in the solar intensity spectrum and find an excellent match. Results. We find that the observed spectral feature is a stellar line originating from the 4s'-4p' transition in S I, and that the reason for its absence in atomic line databases is a neglected air-to-vacuum correction in the original laboratory measurements from 1967 for this line only. From interpolation we determine the laboratory wavelength of the S I line to be lambda(vac) = 1063.8908 nm or lambda(air) = 1063.5993 nm, and the excitation energy of the upper level to be 9.74978 eV.


title: The Role of the Parker Instability in Structuring the Interstellar Medium

abstract: The Parker instability, a Rayleigh-Taylor like instability of thermal gas supported against gravity by magnetic fields and cosmic rays, is thought to be dynamically important for galaxy evolution, possibly promoting molecular cloud formation and the galactic dynamo. In previous work, we examined the effect of three different cosmic ray transport models on the Parker instability: decoupled ($\gamma_c = 0$), locked to the thermal gas ($\gamma_c = 4/3$) and coupled to the gas with streaming by self-confinement. We expand upon that work here by considering radiative cooling, a smooth gravitational potential, and simulations into the nonlinear regime. We determine that cosmic ray transport away from compression points, whether by diffusion or streaming, is the largest driver of the instability. Heating due to cosmic ray streaming is also destabilizing and especially affects the nonlinear regime. While cooling de-pressurizes the dense gas, streaming cosmic rays heat and inflate the diffuse extraplanar gas, greatly modifying the phase structure of the medium. In 3D, we find that the fastest growth favors short wavelength modes in the horizontal direction perpendicular to the background magnetic field; this is imprinted on Faraday rotation measure maps that may be used to detect the Parker instability. The modifications to the Parker instability that we observe in this work have large implications for the structure and evolution of galaxies, and they highlight the major role that cosmic rays play in shaping their environments.


title: Magnon Bose-Einstein Condensation and Superconductivity in a Frustrated Kondo Lattice

abstract: Motivated by recent experiments on magnetically frustrated heavy fermion metals, we theoretically study the phase diagram of the Kondo lattice model with a nonmagnetic valence bond solid ground state on a ladder. A similar physical setting may be naturally occurring in YbAl$_3$C$_3$, CeAgBi$_2$, and TmB$_4$ compounds. In the insulating limit, the application of a magnetic field drives a quantum phase transition to an easy-plane antiferromagnet, which is described by a Bose-Einstein condensation of magnons. Using a combination of field theoretical techniques and density matrix renormalization group calculations we demonstrate that in one dimension this transition is stable in the presence of a metallic Fermi sea and its universality class in the local magnetic response is unaffected by the itinerant gapless fermions. Moreover, we find that fluctuations about the valence bond solid ground state can mediate an attractive interaction that drives unconventional superconducting correlations. We discuss the extensions of our findings to higher dimensions and argue that, depending on the filling of conduction electrons, the magnon Bose-Einstein condensation transition can remain stable in a metal also in dimensions two and three.


title: Weighing the stellar constituents of the Galactic halo with APOGEE red giant stars

abstract: The stellar mass in the halo of the Milky Way is notoriously difficult to determine, owing to the paucity of its stars in the solar neighbourhood. With tentative evidence from $\it{Gaia}$ that the nearby stellar halo is dominated by a massive accretion event -- referred to as $\it{Gaia-Enceladus}$ or Sausage -- these constraints are now increasingly urgent. We measure the mass in kinematically selected mono-abundance populations (MAPs) of the stellar halo between $-3 <$ [Fe/H] $< -1$ and $0.0 <$ [Mg/Fe] $< 0.4$ using red giant star counts from APOGEE DR14. We find that MAPs are well fit by single power laws on triaxial ellipsoidal surfaces, and we show that that the power law slope $\alpha$ changes such that high $\mathrm{[Mg/Fe]}$ populations have $\alpha \sim 4$, whereas low [Mg/Fe] MAPs are more extended with shallow slopes, $\alpha \sim 2$. We estimate the total stellar mass to be $M_{*,\mathrm{tot}} = 1.3^{+0.3}_{-0.2}\times10^{9}\ \mathrm{M_{\odot}}$, of which we estimate $\sim 0.9^{+0.2}_{-0.1} \times 10^{9} \mathrm{M_{\odot}}$ to be accreted. We estimate that the mass of accreted stars with $e > 0.7$ is $M_{*,\mathrm{accreted},e>0.7} = 3 \pm 1\ \mathrm{(stat.)}\pm 1\ \mathrm{(syst.)}\times10^{8}\ \mathrm{M_\odot}$, or $\sim 30-50\%$ of the accreted halo mass. If the majority of these stars $\it{are}$ the progeny of a massive accreted dwarf, this places an upper limit on its stellar mass, and implies a halo mass for the progenitor of $\sim 10^{10.2\pm0.2}\ \mathrm{M_\odot}$. This constraint not only shows that the $\it{Gaia-Enceladus}$/Sausage progenitor may not be as massive as originally suggested, but that the majority of the Milky Way stellar halo was accreted. These measurements are an important step towards fully reconstructing the assembly history of the Milky Way.


title: Robust and efficient algorithms for high-dimensional black-box quantum optimization

abstract: Hybrid quantum-classical optimization using near-term quantum technology is an emerging direction for exploring quantum advantage in high-dimensional systems. However, precise characterization of all experimental parameters is often impractical and challenging. A viable approach is to use algorithms that rely only on black-box inference rather than analytical gradients. Here, we combine randomized perturbation gradient estimation with adaptive momentum gradient updates to create the AdamSPSA and AdamRSGF algorithms. We prove the asymptotic convergence of our algorithms in a convex setting, and we benchmark them against other gradient-based optimization algorithms on non-convex optimal control tasks. Our results show that these new algorithms accelerate the convergence rate, decrease the variance of loss trajectories, and efficiently tune up high-fidelity (above 99.9\%) Hann-window single-qubit gates from trivial initial conditions with twenty variables.


title: A PSF-based Approach to TESS High quality data Of Stellar clusters (PATHOS) -- I. Search for exoplanets and variable stars in the field of 47 Tuc

abstract: The TESS mission will survey ~85 % of the sky, giving us the opportunity of extracting high-precision light curves of millions of stars, including stellar cluster members. In this work, we present our project "A PSF-based Approach to TESS High quality data Of Stellar clusters" (PATHOS), aimed at searching and characterise candidate exoplanets and variable stars in stellar clusters using our innovative method for the extraction of high-precision light curves of stars located in crowded environments. Our technique of light-curve extraction involves the use of empirical Point Spread Functions (PSFs), an input catalogue and neighbour-subtraction. The PSF-based approach allows us to minimise the dilution effects in crowded environments and to extract high-precision photometry for stars in the faint regime (G>13). For this pilot project, we extracted, corrected, and analysed the light curves of 16641 stars located in a dense region centred on the globular cluster 47 Tuc. We were able to reach the TESS magnitude T~16.5 with a photometric precision of ~1 % on the 6.5-hour timescale; in the bright regime we were able to detect transits with depth of ~34 parts per million. We searched for variables and candidate transiting exoplanets. Our pipeline detected one planetary candidate orbiting a main sequence star in the Galactic field. We analysed the period-luminosity distribution for red-giant stars of 47 Tuc and the eclipsing binaries in the field. Light curves are uploaded on the Mikulski Archive for Space Telescopes under the project PATHOS.


title: MUSEQuBES: Calibrating the redshifts of Ly {\alpha} emitters using stacked circumgalactic medium absorption profiles

abstract: Lyman$-\alpha$ (Ly$\alpha$) emission lines are typically found to be redshifted with respect to the systemic redshifts of galaxies, likely due to resonant scattering of Ly$\alpha$ photons. Here we measure the average velocity offset for a sample of 96 $z\approx3.3$ Ly$\alpha$ emitters (LAEs) with a median Ly$\alpha$ flux (luminosity) of $\approx 10^{-17}~\rm erg~cm^{-2}~s^{-1}$ ($\approx10^{42}~\rm erg~s^{-1}$) and a median star formation rate (SFR) of $\approx1.3 \rm M_{\odot} yr^{-1}$ (not corrected for possible dust extinction), detected by the Multi-Unit Spectroscopic Explorer as part of our MUSEQuBES circumgalactic medium (CGM) survey. By postulating that the stacked CGM absorption profiles of these LAEs, probed by 8 background quasars, must be centered on the systemic redshift, we measure an average velocity offset, V$_{\rm offset} = 171 \pm 8$ $\rm km s^{-1}$, between the Ly$\alpha$ emission peak and the systemic redshift. The observed V$_{\rm offset}$ is lower by factors of $\approx1.4$ and $\approx2.6$ compared to the velocity offsets measured for narrow-band selected LAEs and Lyman break galaxies, respectively, which probe galaxies with higher masses and SFRs. Consistent with earlier studies based on direct measurements for individual objects, we find that the V$_{\rm offset}$ is correlated with the full width at half-maximum of the red peak of the Ly$\alpha$ line, and anti-correlated with the rest-frame equivalent width. Moreover, we find that $V_{\rm offset}$ is correlated with SFR with a sub-linear scaling relation, V$_{\rm offset}\propto \rm SFR^{0.16\pm0.03}$. Adopting the mass scaling for main sequence galaxies, such a relation suggests that V$_{\rm offset}$ scales with the circular velocity of the dark matter halos hosting the LAEs.


title: Distance Conjecture and De-Sitter Quantum Gravity

abstract: Motivated by the early discovery of the gigantic landscape of string theory vacua, in recent years people switched direction to try to find constraints on low energy effective field theories from UV-complete descriptions for example quantum gravity or string theory. Those constraints are called swampland criteria. Among the proposed swampland criteria the distance conjecture is the most well-motivated one from the universal behavior of string theory compactifications. It claims that for a type of scalar fields, called moduli, the effective field theory description will break down if it moves a large distance in field (moduli) space. In this work, we will show that in one proposal for a complete theory of quantum gravity including positive vacuum energy, the so called DS/dS correspondence, there is a strong signal of the validity of the distance conjecture.


title: The light CP-even MSSM Higgs mass including N$^3$LO+N$^3$LL QCD corrections

abstract: We present a calculation of the light neutral CP-even Higgs boson pole mass in the real MSSM which combines state-of-the-art EFT and fixed-order results, including the three-loop fixed-order QCD corrections as well as the resummation of logarithmic terms in the ratio of the weak to the SUSY scale up to fourth logarithmic order. This hybrid calculation should be valid for arbitrary SUSY scales above the weak scale. Comparison to the pure fixed-order and EFT results provides an estimate of their individual validity range.


title: Multi-wavelength properties of radio and machine-learning identified counterparts to submillimeter sources in S2COSMOS

abstract: We identify multi-wavelength counterparts to 1,147 submillimeter sources from the S2COSMOS SCUBA-2 survey of the COSMOS field by employing a recently developed radio$+$machine-learning method trained on a large sample of ALMA-identified submillimeter galaxies (SMGs), including 260 SMGs identified in the AS2COSMOS pilot survey. In total, we identify 1,222 optical/near-infrared(NIR)/radio counterparts to the 897 S2COSMOS submillimeter sources with S$_{850}$>1.6mJy, yielding an overall identification rate of ($78\pm9$)%. We find that ($22\pm5$)% of S2COSMOS sources have multiple identified counterparts. We estimate that roughly 27% of these multiple counterparts within the same SCUBA-2 error circles very likely arise from physically associated galaxies rather than line-of-sight projections by chance. The photometric redshift of our radio$+$machine-learning identified SMGs ranges from z=0.2 to 5.7 and peaks at $z=2.3\pm0.1$. The AGN fraction of our sample is ($19\pm4$)%, which is consistent with that of ALMA SMGs in the literature. Comparing with radio/NIR-detected field galaxy population in the COSMOS field, our radio+machine-learning identified counterparts of SMGs have the highest star-formation rates and stellar masses. These characteristics suggest that our identified counterparts of S2COSMOS sources are a representative sample of SMGs at z<3. We employ our machine-learning technique to the whole COSMOS field and identified 6,877 potential SMGs, most of which are expected to have submillimeter emission fainter than the confusion limit of our S2COSMOS surveys (S$_{850}$<1.5mJy). We study the clustering properties of SMGs based on this statistically large sample, finding that they reside in high-mass dark matter halos ($(1.2\pm0.3)\times10^{13}\,h^{-1}\,\rm M_{\odot}$), which suggests that SMGs may be the progenitors of massive ellipticals we see in the local Universe.


title: Baryonic effects on the matter bispectrum

abstract: The large-scale clustering of matter is impacted by baryonic physics, particularly AGN feedback. Modelling or mitigating this impact will be essential for making full use of upcoming measurements of cosmic shear and other large-scale structure probes. We study baryonic effects on the matter bispectrum, using measurements from a selection of state-of-the-art hydrodynamical simulations: IllustrisTNG, Illustris, EAGLE, and BAHAMAS. We identify a low-redshift enhancement of the bispectrum, peaking at $k\sim 3h\,{\rm Mpc}^{-1}$, that is present in several simulations, and discuss how it can be associated to the evolving nature of AGN feedback at late times. This enhancement does not appear in the matter power spectrum, and therefore represents a new source of degeneracy breaking between two- and three-point statistics. In addition, we provide physical interpretations for other aspects of these measurements, and make initial comparisons to predictions from perturbation theory, empirical fitting formulas, and the response function formalism. We publicly release our measurements (including estimates of their uncertainty due to sample variance) and bispectrum measurement code as resources for the community.


title: Forward Modeling of Double Neutron Stars: Insights from Highly-Offset Short Gamma-Ray Bursts

abstract: We present a detailed analysis of two well-localized, highly offset short gamma-ray bursts---GRB~070809 and GRB~090515---investigating the kinematic evolution of their progenitors from compact object formation until merger. Calibrating to observations of their most probable host galaxies, we construct semi-analytic galactic models that account for star formation history and galaxy growth over time. We pair detailed kinematic evolution with compact binary population modeling to infer viable post-supernova velocities and inspiral times. By populating binary tracers according to the star formation history of the host and kinematically evolving their post-supernova trajectories through the time-dependent galactic potential, we find that systems matching the observed offsets of the bursts require post-supernova systemic velocities of hundreds of kilometers per second. Marginalizing over uncertainties in the stellar mass--halo mass relation, we find that the second-born neutron star in the GRB~070809 and GRB~090515 progenitor systems received a natal kick of $\gtrsim 200~\mathrm{km\,s}^{-1}$ at the 78\% and 91\% credible levels, respectively. Applying our analysis to the full catalog of localized short gamma-ray bursts will provide unique constraints on their progenitors and help unravel the selection effects inherent to observing transients that are highly offset with respect to their hosts.


title: Emergent gravity from off-shell energy fixing

abstract: Off-shell processes do not preserve the Energy Momentum Tensor (EMT) in QFT. Fixing the EMT throughout off-shell processes, implies a graviton-like quantum field to emerge without dynamics. Its dynamics are generated through quantum corrections. This Fixed Off-Shell Energy Condition (FOSEC) implies the existence of a linear gravity-like theory, and in special cases the full Poincar\`e gauge theory. In this work it is shown that imposing the FOSEC in QFT implies the emergence of a viable quantum theory of gravity.


title: 3d N=4 Bootstrap and Mirror Symmetry

abstract: We investigate the non-BPS realm of 3d ${\cal N} = 4$ superconformal field theory by uniting the non-perturbative methods of the conformal bootstrap and supersymmetric localization, and utilizing special features of 3d ${\cal N} = 4$ theories such as mirror symmetry and a protected sector described by topological quantum mechanics (TQM). Supersymmetric localization allows for the exact determination of the conformal and flavor central charges, and the latter can be fed into the mini-bootstrap of the TQM to solve for a subset of the OPE data. We examine the implications of the $\mathbb{Z}_2$ mirror action for the SCFT single- and mixed-branch crossing equations for the moment map operators, and apply numerical bootstrap to obtain universal constraints on OPE data for given flavor symmetry groups. A key ingredient in applying the bootstrap analysis is the determination of the mixed-branch superconformal blocks. Among other results, we show that the simplest known self-mirror theory with $SU(2) \times SU(2)$ flavor symmetry saturates our bootstrap bounds, which allows us to extract the non-BPS data and examine the self-mirror $\mathbb{Z}_2$ symmetry thereof.


title: The Two LIGO/Virgo Binary Black Hole Mergers on 2019 August 28 Were Not Strongly Lensed

abstract: The LIGO/Virgo gravitational wave events S190828j and S190828l were detected only 21 minutes apart, from nearby regions of sky, and with the same source classifications (binary black hole mergers). It is therefore natural to speculate that the two signals are actually strongly lensed images of the same merger. However, an estimate of the separation of the (unknown) positions of the two events requires them to be >10 deg apart, much wider than the arcsecond-scale separations that usually arise in extragalactic lensing. The large separation is much more consistent with two independent, unrelated events that occurred close in time by chance. We quantify the overlap between simulated pairs of lensed events, and use frequentist hypothesis testing to reject S190828j/l as a lensed pair at 99.8% confidence.


title: One-Loop Yang-Mills Integrands from Scattering Equations

abstract: We investigate in the context of the scattering equations, how one-loop linear propagator integrands in gauge theories can be linked to integrands with quadratic propagators using a double forward limit. We illustrate our procedure through examples and demonstrate how the different parts of the derived quadratic integrand are consistent with cut-integrands derived from four-dimensional generalized unitarity. We also comment on applications and discuss possible further generalizations.


title: Sequences of $6d$ SCFTs on generic Riemann surfaces

abstract: We consider compactifications of $6d$ minimal $(D_{N+3},D_{N+3})$ type conformal matter SCFTs on a generic Riemann surface. We derive the theories corresponding to three punctured spheres (trinions) with three maximal punctures, from which one can construct models corresponding to generic surfaces. The trinion models are simple quiver theories with $\mathcal{N}=1$ $SU(2)$ gauge nodes. One of the three puncture non abelian symmetries is emergent in the IR. The derivation of the trinions proceeds by analyzing RG flows between conformal matter SCFTs with different values of $N$ and relations between their subsequent reductions to $4d$. In particular, using the flows we first derive trinions with two maximal and one minimal punctures, and then we argue that collections of $N$ minimal punctures can be interpreted as a maximal one. This suggestion is checked by matching the properties of the $4d$ models such as `t Hooft anomalies, symmetries, and the structure of the conformal manifold to the expectations from $6d$. We then use the understanding that collections of minimal punctures might be equivalent to maximal ones to construct trinions with three maximal punctures, and then $4d$ theories corresponding to arbitrary surfaces, for $6d$ models described by two $M5$ branes probing a $\mathbb{Z}_k$ singularity. This entails the introduction of a novel type of maximal puncture. Again, the suggestion is checked by matching anomalies, symmetries and the conformal manifold to expectations from six dimensions. These constructions thus give us a detailed understanding of compactifications of two sequences of six dimensional SCFTs to four dimensions.


title: Inhomogeneous chiral condensates in three-flavor quark matter

abstract: We investigate the effect of strange quark degrees of freedom on the formation of inhomogeneous chiral condensates in a three-flavor Nambu--Jona-Lasinio model in mean-field approximation. A Ginzburg-Landau study complemented by a stability analysis allow us to determine in a general way the location of the critical and Lifshitz points, together with the phase boundary where the (partially) chirally restored phase becomes unstable against developing inhomogeneities, without resorting to specific assumptions on the shape of the chiral condensate. We discuss the resulting phase structure and study the influence of the bare strange-quark mass $m_s$ and the axial anomaly on the size and location of the inhomogeneous phase compared to the first-order transition associated with homogeneous matter. We find that, as a consequence of the axial anomaly, critical and Lifshitz point split. For realistic strange-quark masses the effect is however very small and becomes sizeable only for small values of $m_s$.


title: Constraining the scatter in the galaxy-halo connection at Milky Way masses

abstract: We develop and implement two new methods for constraining the scatter in the relationship between galaxies and dark matter halos. These new techniques are sensitive to the scatter at low halo masses, making them complementary to previous constraints that are dependent on clustering amplitudes or rich galaxy groups, both of which are only sensitive to more massive halos. In both of our methods, we use a galaxy group finder to locate central galaxies in the SDSS main galaxy sample. Our first technique uses the small-scale cross-correlation of central galaxies with all lower mass galaxies. This quantity is sensitive to the satellite fraction of low-mass galaxies, which is in turn driven by the scatter between halos and galaxies. The second technique uses the kurtosis of the distribution of line-of-sight velocities between central galaxies and neighboring galaxies. This quantity is sensitive to the distribution of halo masses that contain the central galaxies at fixed stellar mass. Theoretical models are constructed using peak halo circular velocity, $V_{\rm peak}$, as our property to connect galaxies to halos. The cross-correlation technique yields a constraint of $\sigma[ M_\ast|V_{\rm peak}]=0.27\pm 0.05$ dex, corresponding to a scatter in $\log M_\ast$ at fixed $M_h$ of $\sigma[ M_\ast|M_h]=0.38\pm 0.06$ dex at $M_h=10^{11.8}$ Msun. The kurtosis technique yields $\sigma[ M_\ast|V_{\rm peak}]=0.30\pm0.03$, corresponding to $\sigma[ M_\ast|M_h]=0.34\pm 0.04$ at $M_h=10^{12.2}$ Msun. The values of $\sigma[ M_\ast|M_h]$ are significantly larger than the constraints at higher masses, in agreement with the results of hydrodynamic simulations. This increase is only partly due to the scatter between $V_{\rm peak}$ and $M_h$, and it represents an increase of nearly a factor of two relative to the values inferred from clustering and group studies at high masses.


title: O new physics, where art thou? A global search in the top sector

abstract: We provide a comprehensive global analysis of Run II top measurements at the LHC in terms of dimension-6 operators. A distinctive feature of the top sector as compared to the Higgs-electroweak sector is the large number of four-quark operators. We discuss in detail how they can be tested and how quadratic terms lead to a stable limit on each individual Wilson coefficient. Predictions for all observables are computed at NLO in QCD. Our SFitter analysis framework features a detailed error treatment and shows that theoretical uncertainties are a limiting factor.


title: 2D non-LTE modelling of a filament observed in the H_alpha line with the DST/IBIS spectropolarimeter

abstract: We study a fragment of a large quiescent filament observed on May 29, 2017 by the Interferometric BIdimensional Spectropolarimeter (IBIS) mounted at the Dunn Solar Telescope. We focus on its quiescent stage prior to its eruption. We analyse the spectral observations obtained in the H$\alpha$ line to derive the thermodynamic properties of the plasma of the observed fragment of the filament. We used a 2D filament model employing radiative transfer computations under conditions that depart from the local thermodynamic equilibrium. We employed a forward modelling technique in which we used the 2D model to producesynthetic H_alpha line profiles that we compared with the observations. We then found the set of model input parameters, which produces synthetic spectra with the best agreement with observations. Our analysis shows that one part of the observed fragment of the filament is cooler, denser, and more dynamic than its other part that is hotter, less dense, and more quiescent. The derived temperatures in the first part range from 6,000 K to 10,000$ K and in the latter part from 11,000 K to 14,000 K. The gas pressure is 0.2-0.4 dyn/cm}^{2} in the first part and around 0.15 dyn/cm}^{2} in the latter part. The more dynamic nature of the first part is characterised by the line-of-sight velocities with absolute values of 6-7 km/s and microturbulent velocities of 8-9 km/s. On the other hand, the latter part exhibits line-of-sight velocities with absolute values 0-2.5 km/s and microturbulent velocities of 4-6 km/s.


title: Chaotic cold accretion in giant elliptical galaxies heated by AGN cosmic rays

abstract: Black hole feedback plays a central role in shaping the circumgalactic medium (CGM) of elliptical galaxies. We systematically study the impact of plasma physics on the evolution of ellipticals by performing three-dimensional non-ideal magneto-hydrodynamic simulations of the interactions of active galactic nucleus (AGN) jets with the CGM including magnetic fields, and cosmic rays (CRs) and their transport processes. We find that the physics of feedback operating on large galactic scales depends very sensitively on plasma physics operating on small scales. Specifically, we demonstrate that: (i) in the purely hydrodynamical case, the AGN jets initially maintain the atmospheres in global thermal balance. However, local thermal instability generically leads to the formation of massive cold disks in the vicinity of the central black hole in disagreement with observations; (ii) including weak magnetic fields prevents the formation of the disks because local B-field amplification in the precipitating cold gas leads to strong magnetic breaking, which quickly extracts angular momentum from the accreting clouds. The magnetic fields transform the cold clouds into narrow filaments that do not fall ballistically; (iii) when plasma composition in the AGN jets is dominated by CRs, and CR transport is neglected, the atmospheres exhibit cooling catastrophes due to inefficient heat transfer from the AGN to CGM despite Coulomb/hadronic CR losses being present; (iv) including CR streaming and heating restores agreement with the observations, i.e., cooling catastrophes are prevented and massive cold central disks do not form. The AGN power is reduced as its energy is utilized efficiently.


title: Investigating the multiwavelength behaviour of the flat spectrum radio quasar CTA 102 during 2013-2017

abstract: We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.


title: Celestial fields on the string and the Schwarzian action

abstract: This paper describes the motion of a classical Nambu-Goto string in three-dimensional anti-de Sitter spacetime in terms of two `celestial' fields on the worldsheet. The fields correspond to retarded and advanced boundary times at which null rays emanating from the string reach the boundary. The formalism allows for a simple derivation of the Schwarzian action for near-AdS2 embeddings.


title: Reducing Service Deployment Cost Through VNF Sharing

abstract: Thanks to its computational and forwarding capabilities, the mobile network infrastructure can support several third-party ("vertical") services, each composed of a graph of virtual (network) functions (VNFs). Importantly, one or more VNFs are often common to multiple services, thus the services deployment cost could be reduced by letting the services share the same VNF instance instead of devoting a separate instance to each service. By doing that, however, it is critical that the target KPI (key performance indicators) of all services are met. To this end, we study the VNF sharing problem and make decisions on (i) when sharing VNFs among multiple services is possible, (ii) how to adapt the virtual machines running the shared VNFs to the combined load of the assigned services, and (iii) how to prioritize the services traffic within shared VNFs. All decisions aim to minimize the cost for the mobile operator, subject to requirements on end-to-end service performance, e.g., total delay. Notably, we show that the aforementioned priorities should be managed dynamically and vary across VNFs. We then propose the FlexShare algorithm to provide near-optimal VNF-sharing and priority assignment decisions in polynomial time. We prove that FlexShare is within a constant factor from the optimum and, using real-world VNF graphs, we show that it consistently outperforms baseline solutions.


title: Licci binomial edge ideals

abstract: We give a complete characterization of graphs whose binomial edge ideal is licci. An important tool is a new general upper bound for the regularity of binomial edge ideals.


title: Spatial coherence of spin-orbit-coupled Bose gases

abstract: Spin-orbit-coupled Bose-Einstein condensates (SOBECs) exhibit two new phases of matter, now known as the stripe and plane-wave phases. When two interacting spin components of a SOBEC spatially overlap, density modulations with periodicity given by the spin-orbit coupling strength appear. In equilibrium, these components fully overlap in the miscible stripe phase, and overlap only in a domain wall in the immiscible plane-wave phase. Here we probe the density modulation present in any overlapping region with optical Bragg scattering, and observe the sudden drop of Bragg scattering as the overlapping region shrinks. Using an atomic analogue of the Talbot effect, we demonstrate the existence of long-range coherence between the different spin components in the stripe phase and surprisingly even in the phase-separated plane-wave phase.


title: A year-long plateau in the late-time near-infrared light curves of Type Ia supernovae

abstract: The light curves of Type Ia supernovae are routinely used to constrain cosmology models. Driven by radioactive decay of 56Ni, the light curves steadily decline over time, but >150 days past explosion, the near-infrared portion is poorly characterized. We report a year-long plateau in the near-infrared light curve at 150-500 days, followed by a second decline phase accompanied by a possible appearance of [Fe I] emission lines. This near-infrared plateau contrasts sharply with Type IIP plateaus and requires a new physical mechanism. We suggest a such as masking of the "near-infrared catastrophe," a predicted yet unobserved sharp light-curve decline, by scattering of ultraviolet photons to longer wavelengths. The transition off the plateau could be due to a change in the dominant ionization state of the supernova ejecta. Our results shed new light on the complex radiative transfer processes that take place in Type Ia supernovae and enhance their use as "standard candles."


title: Non-Existence of Finite Order Solution of Non-homogeneous Second Order Linear Differential Equations

abstract: In this paper, we have considered second order non-homogeneous linear differential equations having entire coefficients. We have established conditions ensuring non-existence of finite order solution of such type of differential equations.