title: Targeted Learning Ensembles for Optimal Individualized Treatment Rules with Time-to-Event Outcomes

abstract: We consider estimation of an optimal individualized treatment rule from observational and randomized studies when a high-dimensional vector of baseline variables is available. Our optimality criterion is with respect to delaying expected time to occurrence of an event of interest (e.g., death or relapse of cancer). We leverage semiparametric efficiency theory to construct estimators with desirable properties such as double robustness. We propose two estimators of the optimal rule, which arise from considering two loss functions aimed at (i) directly estimating the conditional treatment effect (also know as the blip function), and (ii) recasting the problem as a weighted classification problem that uses the 0-1 loss function. Our estimated rules are super learning ensembles that minimize the cross-validated risk of a linear combination in a user-supplied library of candidate estimators. We prove oracle inequalities bounding the finite sample excess risk of the estimator. The bounds depend on the excess risk of the oracle selector and a doubly robust term related to estimation of the nuisance parameters. We discuss some important implications of these oracle inequalities such as the convergence rates of the value of our estimator to that of the oracle selector. We illustrate our methods in the analysis of a phase III randomized study testing the efficacy of a new therapy for the treatment of breast cancer.


title: Distributed deep learning on edge-devices: feasibility via adaptive compression

abstract: A large portion of data mining and analytic services use modern machine learning techniques, such as deep learning. The state-of-the-art results by deep learning come at the price of an intensive use of computing resources. The leading frameworks (e.g., TensorFlow) are executed on GPUs or on high-end servers in datacenters. On the other end, there is a proliferation of personal devices with possibly free CPU cycles; this can enable services to run in users' homes, embedding machine learning operations. In this paper, we ask the following question: Is distributed deep learning computation on WAN connected devices feasible, in spite of the traffic caused by learning tasks? We show that such a setup rises some important challenges, most notably the ingress traffic that the servers hosting the up-to-date model have to sustain. In order to reduce this stress, we propose adaComp, a novel algorithm for compressing worker updates to the model on the server. Applicable to stochastic gradient descent based approaches, it combines efficient gradient selection and learning rate modulation. We then experiment and measure the impact of compression, device heterogeneity and reliability on the accuracy of learned models, with an emulator platform that embeds TensorFlow into Linux containers. We report a reduction of the total amount of data sent by workers to the server by two order of magnitude (e.g., 191-fold reduction for a convolutional network on the MNIST dataset), when compared to a standard asynchronous stochastic gradient descent, while preserving model accuracy.


title: Nearest Labelset Using Double Distances for Multi-label Classification

abstract: Multi-label classification is a type of supervised learning where an instance may belong to multiple labels simultaneously. Predicting each label independently has been criticized for not exploiting any correlation between labels. In this paper we propose a novel approach, Nearest Labelset using Double Distances (NLDD), that predicts the labelset observed in the training data that minimizes a weighted sum of the distances in both the feature space and the label space to the new instance. The weights specify the relative tradeoff between the two distances. The weights are estimated from a binomial regression of the number of misclassified labels as a function of the two distances. Model parameters are estimated by maximum likelihood. NLDD only considers labelsets observed in the training data, thus implicitly taking into account label dependencies. Experiments on benchmark multi-label data sets show that the proposed method on average outperforms other well-known approaches in terms of Hamming loss, 0/1 loss, and multi-label accuracy and ranks second after ECC on the F-measure.


title: Mini-AOD: A New Analysis Data Format for CMS

abstract: The CMS experiment has developed a new analysis object format ("Mini-AOD") targeting approximately 10% of the size of the Run 1 AOD format. The motivation for the Mini-AOD format is to have a small and quickly derived data format from which the majority of CMS analysis users can start their analysis work. This format is targeted at having sufficient information to serve about 80% of CMS analysis, while dramatically simplifying the disk and I/O resources needed for analysis. Such large reductions were achieved using a number of techniques, including defining light-weight physics-object candidate representations, increasing transverse momentum thresholds for storing physics-object candidates, and reduced numerical precision when it is not required at the analysis level. In this contribution we discuss the critical components of the Mini-AOD format, our experience with its deployment and the planned physics analysis flow for Run 2 based on the Mini-AOD.


title: Support Vector Machines and generalisation in HEP

abstract: We review the concept of Support Vector Machines (SVMs) and discuss examples of their use in a number of scenarios. Several SVM implementations have been used in HEP and we exemplify this algorithm using the Toolkit for Multivariate Analysis (TMVA) implementation. We discuss examples relevant to HEP including background suppression for $H\to\tau^+\tau^-$ at the LHC with several different kernel functions. Performance benchmarking leads to the issue of generalisation of hyper-parameter selection. The avoidance of fine tuning (over training or over fitting) in MVA hyper-parameter optimisation, i.e. the ability to ensure generalised performance of an MVA that is independent of the training, validation and test samples, is of utmost importance. We discuss this issue and compare and contrast performance of hold-out and k-fold cross-validation. We have extended the SVM functionality and introduced tools to facilitate cross validation in TMVA and present results based on these improvements.


title: A geometrical assessment of spinorial energy conditions

abstract: We consider the problem of energy for spinor fields coupled to their surrounding curved-twisted space-time, and we show that when treated geometrically we cannot even be certain that there is a problem for the energy in the first place.


title: A density version for H\"aggstr\"om's theorem

abstract: Given invariant percolation on a regular tree, where the probability of an edge to be open equals $p$, is it always possible to find an infinite self-avoiding path along which the density of open edges is bigger then $p$?


title: Subsonic Flow over a Thin Airfoil in Ground Effect

abstract: In this paper, the problem of compressible flow over a thin airfoil located near the ground is studied. A singular integral equation, also known as Possio equation, that relates the pressure jump along the airfoil to its downwash is derived. The derivation of the equation utilizes Laplace transform, Fourier transform, method of images, and theory of Mikhlin multipliers. The existence and uniqueness of solution to the Possio equation is verified for the steady state case and an approximate solution is obtained. The aerodynamic loads are then calculated based on the approximate solution. Moreover, the divergence speed of a continuum wing structure located near the ground is obtained based on the derived expressions for the aerodynamic loads.


title: Simple rules for complex decisions

abstract: From doctors diagnosing patients to judges setting bail, experts often base their decisions on experience and intuition rather than on statistical models. While understandable, relying on intuition over models has often been found to result in inferior outcomes. Here we present a new method, select-regress-and-round, for constructing simple rules that perform well for complex decisions. These rules take the form of a weighted checklist, can be applied mentally, and nonetheless rival the performance of modern machine learning algorithms. Our method for creating these rules is itself simple, and can be carried out by practitioners with basic statistics knowledge. We demonstrate this technique with a detailed case study of judicial decisions to release or detain defendants while they await trial. In this application, as in many policy settings, the effects of proposed decision rules cannot be directly observed from historical data: if a rule recommends releasing a defendant that the judge in reality detained, we do not observe what would have happened under the proposed action. We address this key counterfactual estimation problem by drawing on tools from causal inference. We find that simple rules significantly outperform judges and are on par with decisions derived from random forests trained on all available features. Generalizing to 22 varied decision-making domains, we find this basic result replicates. We conclude with an analytical framework that helps explain why these simple decision rules perform as well as they do.


title: Resonant properties of dipole skyrmions in amorphous Fe/Gd multilayers

abstract: The dynamic response of dipole skyrmions in Fe/Gd multilayer films is investigated by ferromagnetic resonance measurements and compared to micromagnetic simulations. We detail thickness and temperature dependent studies of the observed modes as well as the effects of magnetic field history on the resonant spectra. Correlation between the modes and the magnetic phase maps constructed from real-space imaging and scattering patterns allows us to conclude the resonant modes arise from local topological features such as dipole skyrmions but does not depend on the collective response of a closed packed lattice of these chiral textures. Using, micromagnetic modeling, we are able to quantitatively reproduce our experimental observations which suggests the existence of localized spin-wave modes that are dependent on the helicity of the dipole skyrmion. We identify four localized spin wave excitations for the skyrmions that are excited under either in-plane or out-of-plane r.f. fields. Lastly we show that dipole skyrmions and non-chiral bubble domains exhibit qualitatively different localized spin wave modes.


title: Bright-field microscopy of transparent objects: a ray tracing approach

abstract: Formation of a bright-field microscopic image of a transparent phase object is described in terms of elementary geometrical optics. Our approach is based on the premise that image replicates the intensity distribution (real or virtual) at the front focal plane of the objective. The task is therefore reduced to finding the change in intensity at the focal plane caused by the object. This can be done by ray tracing complemented with the requirement of conservation of the number of rays. Despite major simplifications involved in such an analysis, it reproduces some results from the paraxial wave theory. Additionally, our analysis suggests two ways of extracting quantitative phase information from bright-field images: by vertically shifting the focal plane (the approach used in the transport-of-intensity analysis) or by varying the angle of illumination. In principle, information thus obtained should allow reconstruction of the object morphology.


title: Is your software on dope? Formal analysis of surreptitiously "enhanced" programs

abstract: Usually, it is the software manufacturer who employs verification or testing to ensure that the software embedded in a device meets its main objectives. However, these days we are confronted with the situation that economical or technological reasons might make a manufacturer become interested in the software slightly deviating from its main objective for dubious reasons. Examples include lock-in strategies and the $\mathrm{NO}_x$ emission scandals in automotive industry. This phenomenon is what we call software doping. It is turning more widespread as software is embedded in ever more devices of daily use. The primary contribution of this article is to provide a hierarchy of simple but solid formal definitions that enable to distinguish whether a program is clean or doped. Moreover, we show that these characterisations provide an immediate framework for analysis by using already existing verification techniques. We exemplify this by applying self-composition on sequential programs and model checking of HyperLTL formulas on reactive models.


title: A note on complete classification of $(\delta+\alpha u^2)$-constacyclic codes of length $p^k$ over $\F_{p^m}+u\F_{p^m}+u^2\F_{p^m}$

abstract: For units $\delta$ and $\alpha$ in $\F_{p^m}$, the structure of $(\delta+\alpha u^2)$-constacyclic codes of length $p^k$ over $\mathbb{F}_{p^m}+u\mathbb{F}_{p^m}+u^2\mathbb{F}_{p^m}$ is studied and self-dual $(\delta+\alpha u^2)$-constacyclic codes are analyzed.


title: Motility of active nematic films driven by "active anchoring"

abstract: We provide a minimal model for an active nematic film in contact with both a solid substrate and a passive isotropic fluid, and explore its dynamics in one and two dimensions using a combination of hybrid Lattice Boltzmann simulations and analytical calculations. By imposing nematic anchoring at the substrate while active flows induce a preferred alignment at the interface ("active anchoring"), we demonstrate that directed fluid flow spontaneously emerges in cases where the two anchoring types are opposing. In one dimension, our model reduces to an analogue of a loaded elastic column. Here, the transition from a stationary to a motile state is akin to the buckling bifurcation, but offers the possilibity to reverse the flow direction for a given set of parameters and boundary conditions solely by changing initial conditions. The two-dimensional variant of our model allows for additional tangential instabilities, leading to self-assembled propagating surface waves for intermediate activity and for a continously deforming irregular surface at high activity. Our results might be relevant for designing active microfluidic geometries, but also for curvature-guided self-assembly or switchable diffraction gratings.


title: On the nature of the generating series of walks in the quarter plane

abstract: In the present paper, we introduce a new approach, relying on the Galois theory of difference equations, to study the nature of the generating series of walks in the quarter plane. Using this approach, we are not only able to recover many of the recent results about these series, but also to go beyond them. For instance, we give for the first time hypertranscendency results, {\it i.e.}, we prove that certain of these generating series do not satisfy any nontrivial nonlinear algebraic differential equation with rational coefficients.


title: Measurements of entanglement over a kilometric distance to test superluminal models of Quantum Mechanics: preliminary results

abstract: As shown in the \emph{EPR} paper (Einstein, Podolsky e Rosen, 1935), Quantum Mechanics is a non-local Theory. The Bell theorem and the successive experiments ruled out the possibility of explaining quantum correlations using only local hidden variables models. Some authors suggested that quantum correlations could be due to superluminal communications that propagate isotropically with velocity \emph{$v_{t}>c$} in a preferred reference frame. For finite values of \emph{$v_{t}$} and in some special cases, Quantum Mechanics and superluminal models lead to different predictions. So far, no deviations from the predictions of Quantum Mechanics have been detected and only lower bounds for the superluminal velocities \emph{$v_{t}$} have been established. Here we describe a new experiment that increases the maximum detectable superluminal velocities and we give some preliminary results.


title: A characterization of a class of convex log-Sobolev inequalities on the real line

abstract: We give a sufficient and necessary condition for a probability measure $\mu$ on the real line to satisfy the logarithmic Sobolev inequality for convex functions. The condition is expressed in terms of the unique left-continuous and non-decreasing map transporting the symmetric exponential measure onto $\mu$. The main tool in the proof is the theory of weak transport costs. As a consequence, we obtain dimension-free concentration bounds for the lower and upper tails of convex functions of independent random variables which satisfy the convex log-Sobolev inequality.


title: Model Predictive Control for Distributed Microgrid Battery Energy Storage Systems

abstract: This paper proposes a new convex model predictive control strategy for dynamic optimal power flow between battery energy storage systems distributed in an AC microgrid. The proposed control strategy uses a new problem formulation, based on a linear d-q reference frame voltage-current model and linearised power flow approximations. This allows the optimal power flows to be solved as a convex optimisation problem, for which fast and robust solvers exist. The proposed method does not assume real and reactive power flows are decoupled, allowing line losses, voltage constraints and converter current constraints to be addressed. In addition, non-linear variations in the charge and discharge efficiencies of lithium ion batteries are analysed and included in the control strategy. Real-time digital simulations were carried out for an islanded microgrid based on the IEEE 13 bus prototypical feeder, with distributed battery energy storage systems and intermittent photovoltaic generation. It is shown that the proposed control strategy approaches the performance of a strategy based on non-convex optimisation, while reducing the required computation time by a factor of 1000, making it suitable for a real-time model predictive control implementation.


title: Target assignment for robots constrained by limited communication range

abstract: This paper investigates the task assignment problem for multiple dispersed robots constrained by limited communication range. The robots are initially randomly distributed and need to visit several target locations while minimizing the total travel time. A centralized rendezvous-based algorithm is proposed, under which all the robots first move towards a rendezvous position until communication paths are established between every pair of robots either directly or through intermediate peers, and then one robot is chosen as the leader to make a centralized task assignment for the other robots. Furthermore, we propose a decentralized algorithm based on a single-traveling-salesman tour, which does not require all the robots to be connected through communication. We investigate the variation of the quality of the assignment solutions as the level of information sharing increases and as the communication range grows, respectively. The proposed algorithms are compared with a centralized algorithm with shared global information and a decentralized greedy algorithm respectively. Monte Carlo simulation results show the satisfying performance of the proposed algorithms.


title: Effects of transients in LIGO suspensions on searches for gravitational waves

abstract: This paper presents an analysis of the transient behavior of the Advanced LIGO suspensions used to seismically isolate the optics. We have characterized the transients in the longitudinal motion of the quadruple suspensions during Advanced LIGO's first observing run. Propagation of transients between stages is consistent with modelled transfer functions, such that transient motion originating at the top of the suspension chain is significantly reduced in amplitude at the test mass. We find that there are transients seen by the longitudinal motion monitors of quadruple suspensions, but they are not significantly correlated with transient motion above the noise floor in the gravitational wave strain data, and therefore do not present a dominant source of background noise in the searches for transient gravitational wave signals.


title: Accretion disk coronae of Intermediate Polar Cataclysmic Variables - 3D MagnetoHydro-Dynamic modeling and thermal X-ray emission

abstract: IPCVs contain a magnetic, rotating white dwarf surrounded by a magnetically truncated accretion disk. To explain their strong flickering X-ray emission, accretion has been successfully taken into account. Nevertheless, observations suggest that accretion phenomena could not be the only process behind it. An intense flaring activity occurring on the surface of the disk may generate a corona, contribute to the thermal X-ray emission and influence the system stability. Our purposes are: investigating the formation of an extended corona above the accretion disk, due to an intense flaring activity occurring on the disk surface; studying its effects on the disk and stellar magnetosphere; assessing its contribution to the observed X-ray flux. We have developed a 3D MHD model of a IPCV. The model takes into account gravity, disk viscosity, thermal conduction, radiative losses and coronal flare heating. To perform a parameter space exploration, several system conditions have been considered, with different magnetic field intensity and disk density values. From the results of the evolution of the model, we have synthesized the thermal X-ray emission. The simulations show the formation of an extended corona, linking disk and star. The flaring activity is capable of strongly influencing the disk configuration and its stability, effectively deforming the magnetic field lines. Hot plasma evaporation phenomena occur in the layer immediately above the disk. The flaring activity gives rise to a thermal X-ray emission in both the [0.1-2.0] keV and the [2.0-10] keV bands. An intense coronal activity occurring on the disk surface of an IPCV can affect the structure of the disk depending noticeably on the density of the disk and the magnetic field of the central object. Moreover, the synthesis of the thermal X-ray fluxes shows that this flaring activity may contribute to the observed thermal X-ray emission.


title: Jamming Resistant Receivers for Massive MIMO

abstract: We design jamming resistant receivers to enhance the robustness of a massive MIMO uplink channel against jamming. In the pilot phase, we estimate not only the desired channel, but also the jamming channel by exploiting purposely unused pilot sequences. The jamming channel estimate is used to construct the linear receive filter to reduce impact that jamming has on the achievable rates. The performance of the proposed scheme is analytically and numerically evaluated. These results show that the proposed scheme greatly improves the rates, as compared to conventional receivers. Moreover, the proposed schemes still work well with stronger jamming power.


title: Quantum groups and noncommutative spacetimes with cosmological constant

abstract: Noncommutative spacetimes are widely believed to model some properties of the quantum structure of spacetime at the Planck regime. In this contribution the construction of (anti-)de Sitter noncommutative spacetimes obtained through quantum groups is reviewed. In this approach the quantum deformation parameter $z$ is related to a Planck scale, and the cosmological constant $\Lambda$ plays the role of a second deformation parameter of geometric nature, whose limit $\Lambda\to 0$ provides the corresponding noncommutative Minkowski spacetimes.


title: Periods of meromorphic quadratic differentials and Goldman bracket

abstract: We study symplectic properties of monodromy map for second order linear equation with meromorphic potential having only simple poles on a Riemann surface. We show that the canonical symplectic structure on the cotangent bundle $T^*M_{g,n}$ implies the Goldman bracket on the corresponding character variety under the monodromy map, thereby extending the recent results of the paper of M.Bertola, C.Norton and the author from the case of holomorphic to meromorphic potentials with simple poles.


title: Parity and Time Reversal Symmetry in Hanbury Brown-Twiss Effect

abstract: The current manuscript employs the parity and time reversal symmetry in the Hanbury Brown-Twiss experiment. For this purpose, we develop a general scattering matrix framework founded on the concatenation of many individual compounded scattering processes on the setup. In this way, we derive the general scattering matrix of a parity and time reversal symmetric Hanbury Brown-Twiss experiment (HBT-PT). Within such scattering formulation, we propose a theoretical framework which provides how to measure the symmetry of the system through the correlation function of a pair of particles transmitted through the leads. The correlation function naturally reveal the quantum statistics of both bosons and fermions and demonstrate a very preponderant role of PT symmetry on the HBT experiment. We indicate the formation of both quantum and classical universal Turing machine depending on controllable parameters of the apparatus.


title: Comparison of Polar Decoders with Existing Low-Density Parity-Check and Turbo Decoders

abstract: Polar codes are a recently proposed family of provably capacity-achieving error-correction codes that received a lot of attention. While their theoretical properties render them interesting, their practicality compared to other types of codes has not been thoroughly studied. Towards this end, in this paper, we perform a comparison of polar decoders against LDPC and Turbo decoders that are used in existing communications standards. More specifically, we compare both the error-correction performance and the hardware efficiency of the corresponding hardware implementations. This comparison enables us to identify applications where polar codes are superior to existing error-correction coding solutions as well as to determine the most promising research direction in terms of the hardware implementation of polar decoders.


title: On correlations between class numbers of imaginary quadratic fields

abstract: Let $h(-n)$ be the class number of the imaginary quadratic field with discriminant $-n$. We establish an asymtotic formula for correlations involving $h(-n)$ and $h(-n-l)$, over fundamental discriminants that avoid the congruence class $1\pmod{8}$. Our result is uniform in the shift $l$, and the proof uses an identity of Gauss relating $h(-n)$ to representations of integers as sums of three squares. We also prove analogous results on correlations involving $r_Q(n)$, the number of representations of an integer $n$ by an integral positive definite quadratic form $Q$.


title: Low-resolution simulations of vesicle suspensions in 2D

abstract: Vesicle suspensions appear in many biological and industrial applications. These suspensions are characterized by rich and complex dynamics of vesicles due to their interaction with the bulk fluid, and their large deformations and nonlinear elastic properties. Many existing state-of-the-art numerical schemes can resolve such complex vesicle flows. However, even when using provably optimal algorithms, these simulations can be computationally expensive, especially for suspensions with a large number of vesicles. In this paper, we investigate the effect of a number of algorithmic empirical fixes in an attempt to make low-resolution simulations more stable and more predictive. Based on our empirical studies for a number of flow configurations, we propose a scheme that attempts to integrate these fixes in a systematic way. This low-resolution scheme is an extension of our previous work Quaife and Biros (2014) and Quaife and Biros (2016). Our low-resolution correction algorithms (LRCA) include anti-aliasing and membrane reparametrization for avoiding spurious oscillations in vesicles' membranes, adaptive time stepping and a repulsion force for handling vesicle collisions and, correction of vesicles' area and arc-length for maintaining physical vesicle shapes. We perform a systematic error analysis by comparing the low-resolution simulations of dilute and dense suspensions with their high-fidelity, fully resolved, counterparts. We observe that the LRCA enables both efficient and statistically accurate low-resolution simulations of vesicle suspensions, while it can be 10x to 100x faster.


title: Multi-Task Convolutional Neural Network for Pose-Invariant Face Recognition

abstract: This paper explores multi-task learning (MTL) for face recognition. We answer the questions of how and why MTL can improve the face recognition performance. First, we propose a multi-task Convolutional Neural Network (CNN) for face recognition where identity classification is the main task and pose, illumination, and expression estimations are the side tasks. Second, we develop a dynamic-weighting scheme to automatically assign the loss weight to each side task, which is a crucial problem in MTL. Third, we propose a pose-directed multi-task CNN by grouping different poses to learn pose-specific identity features, simultaneously across all poses. Last but not least, we propose an energy-based weight analysis method to explore how CNN-based MTL works. We observe that the side tasks serve as regularizations to disentangle the variations from the learnt identity features. Extensive experiments on the entire Multi-PIE dataset demonstrate the effectiveness of the proposed approach. To the best of our knowledge, this is the first work using all data in Multi-PIE for face recognition. Our approach is also applicable to in-the-wild datasets for pose-invariant face recognition and achieves comparable or better performance than state of the art on LFW, CFP, and IJB-A datasets.


title: Quantized Compressed Sensing for Partial Random Circulant Matrices

abstract: We provide the first analysis of a non-trivial quantization scheme for compressed sensing measurements arising from structured measurements. Specifically, our analysis studies compressed sensing matrices consisting of rows selected at random, without replacement, from a circulant matrix generated by a random subgaussian vector. We quantize the measurements using stable, possibly one-bit, Sigma-Delta schemes, and use a reconstruction method based on convex optimization. We show that the part of the reconstruction error due to quantization decays polynomially in the number of measurements. This is in line with analogous results on Sigma-Delta quantization associated with random Gaussian or subgaussian matrices, and significantly better than results associated with the widely assumed memoryless scalar quantization. Moreover, we prove that our approach is stable and robust; i.e., the reconstruction error degrades gracefully in the presence of non-quantization noise and when the underlying signal is not strictly sparse. The analysis relies on results concerning subgaussian chaos processes as well as a variation of McDiarmid's inequality.


title: Invariants of noncommutative projective schemes

abstract: In this note we compute several invariants (e.g. algebraic K-theory, cyclic homology and topological Hochschild homology) of the noncommutative projective schemes associated to Koszul algebras of finite global dimension.


title: Enhanced Quantization: The Right way to Quantize Everything

abstract: Canonical quantization relies on Cartesian, canonical, phase-space coordinates to promote to Hermitian operators, which also become the principal ingredients in the quantum Hamiltonian. While generally appropriate, this procedure can also fail, e.g., for covariant, quartic, scalar fields in five-and-more spacetime dimensions (and possibly four spacetime dimensions as well), which become trivial; such failures are normally blamed on the `problem' rather than on the 'quantization procedure'. In Enhanced Quantization the association of $c$-numbers to $q$-numbers is chosen very differently such that: (i) there is no need to seek classical, Cartesian, phase-space coordinates; (ii) every classical, contact transformation is applicable and no change of the quantum operators arises; (iii) a new understanding of the importance of 'Cartesian coordinates' is established; and (iv) although discussed elsewhere in detail, the procedures of enhanced quantization offer fully acceptable solutions yielding non-trivial results for quartic scalar fields in four-and-more spacetime dimensions. In early sections, this paper offers a wide-audience approach to the basic principles of Enhanced Quantization using simple examples; later, several significant examples are cited for a deeper understanding. An historical note concludes the paper.

categories: quant-ph hep-th math-ph math.MP physics.hist-ph


title: Bosons with incommensurate potential and spin-orbit coupling

abstract: We chart out the phase diagram of ultracold `spin-half' bosons in a one-dimensional optical lattice in the presence of Aubry-Andr\'e (AA) potential and with spin-orbit (SO) and Raman couplings investigating the transition from superfluid (SF) to localized phases and the existence of density wave phase for nearest-neighbor interaction (NNI). We show that the presence of SO coupling and AA potential leads to a novel spin-split momentum distribution of the bosons in the localized phase near the boundary with the SF phase, which can act as a signature of such a transition. We also obtain the level statistics of the bosons in the superfluid phase with finite NNI and demonstrate its change from Gaussian Unitary Ensemble (GUE) to Gaussian Orthogonal Ensemble (GOE) as a function of the Raman coupling. We discuss experiments which can test our theory.


title: Simflowny 2: An upgraded platform for scientific modeling and simulation

abstract: Simflowny is an open platform which automatically generates parallel code of scientific dynamical models for different simulation frameworks. Here we present major upgrades on this software to support an extended set of families of models, in particular: i) a new generic family for partial differential equations, which can include spatial derivatives of any order, ii) a new family for agent based models to study complex phenomena --either on a spatial domain or on a graph--. Additionally we introduce a flexible graphical user interface (GUI) to accommodate these and future families of equations. This paper describes the new GUI architecture and summarizes the formal representation and implementation of these new families, providing several validation results.


title: Stability of Accelerated Expansion in Nonlinear Electrodynamics

abstract: This paper is devoted to study the phase space analysis of isotropic and homogenous universe model by taking a noninteracting mixture of electromagnetic and viscous radiating fluids whose viscous pressure satisfies a nonlinear version of the Israel-Stewart transport equation. We establish an autonomous system of equations by introducing normalized dimensionless variables. In order to analyze stability of the system, we find corresponding critical points for different values of the parameters. We also evaluate power-law scale factor whose behavior indicates different phases of the universe model. It is concluded that bulk viscosity as well as electromagnetic field enhances the stability of accelerated expansion of the isotropic and homogeneous universe model.


title: A quaternary diophantine inequality by prime numbers of a special type

abstract: Let $1<c<832/825$. For large real numbers $N>0$ and a small constant $\vartheta>0$, the inequality \begin{equation*} |p_1^c+p_2^c+p_3^c+p_4^c-N|<\vartheta \end{equation*} has a solution in prime numbers $p_1,\,p_2,\,p_3,\,p_4$ such that, for each $i\in\{1,2,3,4\}$, $p_i+2$ has at most $32$ prime factors.


title: Spectral methods for Langevin dynamics and associated error estimates

abstract: We prove the consistency of Galerkin methods to solve Poisson equations where the differential operator under consideration is the generator of the Langevin dynamics. We show in particular how the hypocoercive nature of this operator can be used at the discrete level to first prove the invertibility of the rigidity matrix, and next provide error bounds on the approximation of the solution of the Poisson equation. We present general convergence results in an abstract setting, as well as explicit convergence rates for a simple one-dimensional example discretized using a tensor basis. Our theoretical findings are illustrated by numerical simulations.


title: Evaluation of Trace Alignment Quality and its Application in Medical Process Mining

abstract: Trace alignment algorithms have been used in process mining for discovering the consensus treatment procedures and process deviations. Different alignment algorithms, however, may produce very different results. No widely-adopted method exists for evaluating the results of trace alignment. Existing reference-free evaluation methods cannot adequately and comprehensively assess the alignment quality. We analyzed and compared the existing evaluation methods, identifying their limitations, and introduced improvements in two reference-free evaluation methods. Our approach assesses the alignment result globally instead of locally, and therefore helps the algorithm to optimize overall alignment quality. We also introduced a novel metric to measure the alignment complexity, which can be used as a constraint on alignment algorithm optimization. We tested our evaluation methods on a trauma resuscitation dataset and provided the medical explanation of the activities and patterns identified as deviations using our proposed evaluation methods.


title: Dynamics of coupled phantom and tachyon fields

abstract: In this paper, we apply the dynamical analysis to a coupled phantom field with scaling potential taking particular forms of the coupling (linear and combination of linear), and present phase space analysis. We investigate if there exist late time accelerated scaling attractor that has the ratio of dark energy and dark matter densities of the order one. We observe that the scrutinized couplings cannot alleviate the coincidence problem, however acquire stable late time accelerated solutions. We also discuss coupled tachyon field with inverse square potential assuming linear coupling.


title: Bo\"otes-HiZELS: an optical to near-infrared survey of emission-line galaxies at $\bf z=0.4-4.7$

abstract: We present a sample of $\sim 1000$ emission line galaxies at $z=0.4-4.7$ from the $\sim0.7$deg$^2$ High-$z$ Emission Line Survey (HiZELS) in the Bo\"otes field identified with a suite of six narrow-band filters at $\approx 0.4-2.1$ $\mu$m. These galaxies have been selected on their Ly$\alpha$ (73), {\sc [Oii]} (285), H$\beta$/{\sc [Oiii]} (387) or H$\alpha$ (362) emission-line, and have been classified with optical to near-infrared colours. A subsample of 98 sources have reliable redshifts from multiple narrow-band (e.g. [O{\sc ii}]-H$\alpha$) detections and/or spectroscopy. In this survey paper, we present the observations, selection and catalogs of emitters. We measure number densities of Ly$\alpha$, [O{\sc ii}], H$\beta$/{\sc [Oiii]} and H$\alpha$ and confirm strong luminosity evolution in star-forming galaxies from $z\sim0.4$ to $\sim 5$, in agreement with previous results. To demonstrate the usefulness of dual-line emitters, we use the sample of dual [O{\sc ii}]-H$\alpha$ emitters to measure the observed [O{\sc ii}]/H$\alpha$ ratio at $z=1.47$. The observed [O{\sc ii}]/H$\alpha$ ratio increases significantly from 0.40$\pm0.01$ at $z=0.1$ to 0.52$\pm0.05$ at $z=1.47$, which we attribute to either decreasing dust attenuation with redshift, or due to a bias in the (typically) fiber-measurements in the local Universe which only measure the central kpc regions. At the bright end, we find that both the H$\alpha$ and Ly$\alpha$ number densities at $z\approx2.2$ deviate significantly from a Schechter form, following a power-law. We show that this is driven entirely by an increasing X-ray/AGN fraction with line-luminosity, which reaches $\approx 100$ \% at line-luminosities $L\gtrsim3\times10^{44}$ erg s$^{-1}$.


title: The Evolution and Properties of Rotating Massive Star Populations

abstract: We investigate the integrated properties of massive (>10 Msun), rotating, single-star stellar populations for a variety of initial rotation rates (v/vcrit=0.0, 0.2, 0.4, 0.5, and 0.6). We couple the new MESA Isochrone and Stellar Tracks (MIST) models to the Flexible Stellar Population Synthesis (FSPS) package, extending the stellar population synthesis models to include the contributions from very massive stars (>100 Msun), which can be significant in the first ~4 Myr after a starburst. These models predict ionizing luminosities that are consistent with recent observations of young nuclear star clusters. We also construct composite stellar populations assuming a distribution of initial rotation rates. Even in low-metallicity environments where rotation has a significant effect on the evolution of massive stars, we find that stellar population models require a significant contribution from fast-rotating (v/vcrit>0.4) stars in order to sustain the production of ionizing photons beyond a few Myr following a starburst. These results have potentially important implications for cosmic reionization by massive stars and the interpretation of nebular emission lines in high-redshift star-forming galaxies.


title: QED$_3$ with quenched disorder: quantum critical states with interactions and disorder

abstract: Quantum electrodynamics in 2+1-dimensions (QED$_3$) is a strongly coupled conformal field theory (CFT) of a U(1) gauge field coupled to $2N$ two-component massless fermions. The $N=2$ CFT has been proposed as a ground state of the spin-1/2 kagome Heisenberg antiferromagnet. We study QED$_3$ in the presence of weak quenched disorder in its two spatial directions. When the disorder explicitly breaks the fermion flavor symmetry from SU($2N$)$\rightarrow$U(1)$\times$SU($N$) but preserves time-reversal symmetry, we find that the theory flows to a non-trivial fixed line at non-zero disorder with a continuously varying dynamical critical exponent $z>1$. We determine the zero-temperature flavor (spin) conductivity along the critical line. Our calculations are performed in the large-$N$ limit, and the disorder is handled using the replica method.


title: On the waveforms of gravitationally lensed gravitational waves

abstract: Strong lensing by intervening galaxies can produce multiple images of gravitational waves from sources at cosmological distances. These images acquire additional phase-shifts as the over-focused wavefront passes through itself along the line of sight. Time domain waveforms of Type-II images (associated with saddle points of the time delay) exhibit a non-trivial distortion from the unlensed waveforms. This phenomenon is in addition to the usual frequency-independent magnification, and happens even in the geometric limit where the wavelength is much shorter than the deflector's gravitational length scale. Similarly, Type-III images preserve the original waveform's shape but exhibit a sign flip. We show that for non-precessing binaries undergoing circular inspiral and merger, these distortions are equivalent to rotating the line of sight about the normal to the orbital plane by $45^\circ$ (Type II) and $90^\circ$ (Type III). This effect will enable us to distinguish between the different topological types among a set of multiple images, and give us valuable insight into the lens model. Furthermore, we show that for eccentric binaries, the waveform of a Type-II image is distorted in a manner that is inequivalent to a change of the source's orbital parameters.


title: CMB-induced radio quenching of high-redshift jetted AGNs with highly magnetic hotspots

abstract: In an effort to understand the cause of the apparent depletion in the number density of radio-loud AGNs at $z>3$, this work investigates the viability of the so-called Cosmic Microwave Background (CMB) quenching mechanism of intrinsically jetted, high-z AGNs, whereby Inverse Compton scattering of CMB photons off electrons within the extended lobes results in a substantial dimming of the lobe synchrotron emission at GHz frequencies, while simultaneously boosting their diffuse X-ray signal. We focus on five $z>3.5$ radio galaxies that have sufficiently deep Chandra exposure (> 50 ks) to warrant a meaningful investigation of any extended X-ray emission. For those objects with evidence for statistically significant extended X-ray lobes (4C 41.17 and 4C 03.24), we combine the Chandra measurements with literature data at lower frequencies to assemble the systems' Spectral Energy Distributions (SEDs), and utilize state-of-the-art SED modelling (Ghisellini et al. 2015) -- including emission from the disk, torus, jet, hotspots, and lobes -- to infer their physical parameters. For both radio galaxies, the magnetic energy density in the hotspots is found to exceed the energy density in CMB photons, wheres the opposite is true for the lobes. This implies that any extended synchrotron emission likely originates from the hotspots themselves, rather than the lobes. Conversely, Inverse Compton scattering of CMB photons dominates the extended X-ray emission from the lobes, which are effectively "radio-quenched". As a result, CMB quenching is effective in these systems in spite of the fact that the observed X-ray to radio luminosity ratio does not bear the signature $(1+z)^4$ dependence of the CMB energy density.


title: A measurement of the z = 0 UV background from H$\alpha$ fluorescence

abstract: We report the detection of extended Halpha emission from the tip of the HI disk of the nearby edge-on galaxy UGC 7321, observed with the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope. The Halpha surface brightness fades rapidly where the HI column density drops below N(HI) = 10^19 cm^-2 , consistent with fluorescence arising at the ionisation front from gas that is photoionized by the extragalactic ultraviolet background (UVB). The surface brightness measured at this location is (1.2 +/- 0.5)x10^-19 erg/s/cm^2/arcsec^2, where the error is mostly systematic and results from the proximity of the signal to the edge of the MUSE field of view, and from the presence of a sky line next to the redshifted Halpha wavelength. By combining the Halpha and the HI 21 cm maps with a radiative transfer calculation of an exponential disk illuminated by the UVB, we derive a value for the HI photoionization rate of Gamma ~ (6-8)x10^-14 1/s . This value is consistent with transmission statistics of the Lyalpha forest and with recent models of a UVB which is dominated by quasars.


title: The Dragonfly Nearby Galaxies Survey. III. The Luminosity Function of the M101 Group

abstract: We obtained follow-up HST observations of the seven low surface brightness galaxies discovered with the Dragonfly Telephoto Array in the field of the massive spiral galaxy M101. Out of the seven galaxies, only three were resolved into stars and are potentially associated with the M101 group at $D=7\text{ Mpc}$. Based on HST ACS photometry in the broad F606W and F814W filters, we use a maximum likelihood algorithm to locate the Tip of the Red Giant Branch (TRGB) in galaxy color-magnitude diagrams. Distances are $6.38^{+0.35}_{-0.35}, 6.87^{+0.21}_{-0.30}$ and $6.52^{+0.25}_{-0.27} \text{ Mpc}$ and we confirm that they are members of the M101 group. Combining the three confirmed low luminosity satellites with previous results for brighter group members, we find the M101 galaxy group to be a sparsely populated galaxy group consisting of seven group members, down to $M_V = -9.2 \text{ mag}$. We compare the M101 cumulative luminosity function to that of the Milky Way and M31. We find that they are remarkably similar; In fact, the cumulative luminosity function of the M101 group gets even flatter for fainter magnitudes, and we show that the M101 group might exhibit the two known small-scale flaws in the $\Lambda\textrm{CDM}$ model, namely `the missing satellite' problem and the `too big to fail' problem. Kinematic measurements of M101$'$s satellite galaxies are required to determine whether the `too big to fail' problem does in fact exist in the M101 group.


title: CHIME FRB: An application of FFT beamforming for a radio telescope

abstract: We have developed FFT beamforming techniques for the CHIME radio telescope, to search for and localize the astrophysical signals from Fast Radio Bursts (FRBs) over a large instantaneous field-of-view (FOV) while maintaining the full angular resolution of CHIME. We implement a hybrid beamforming pipeline in a GPU correlator, synthesizing 256 FFT-formed beams in the North-South direction by four formed beams along East-West via exact phasing, tiling a sky area of ~250 square degrees. A zero-padding approximation is employed to improve chromatic beam alignment across the wide bandwidth of 400 to 800 MHz. We up-channelize the data in order to achieve fine spectral resolution of $\Delta\nu$=24 kHz and time cadence of 0.983 ms, desirable for detecting transient and dispersed signals such as those from FRBs.

categories: astro-ph.IM astro-ph.CO astro-ph.HE astro-ph.SR


title: Evolution of Interstellar Medium, Star Formation, and Accretion at High Redshift

abstract: ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.


title: New constraints and discovery potential of sub-GeV dark matter with xenon detectors

abstract: Existing xenon dark matter (DM) direct detection experiments can probe the DM-nucleon interaction of DM with a sub-GeV mass through a search for photon emission from the recoiling xenon atom. We show that LUX's constraints on sub-GeV DM, which utilise the scintillation (S1) and ionisation (S2) signals, are approximately three orders of magnitude more stringent than previous xenon constraints in this mass range, derived from the XENON10 and XENON100 S2-only searches. The new LUX constraints provide the most stringent direct detection constraints for DM particles with a mass below 0.5 GeV. In addition, the photon emission signal in LUX and its successor LZ maintain the discrimination between background and signal events so that an unambiguous discovery of sub-GeV DM is possible. We show that LZ has the potential to reconstruct the DM mass with 20% accuracy for particles lighter than 0.5 GeV.


title: The Grism Lens-Amplified Survey from Space (GLASS). XI. Detection of CIV in Multiple Images of $z=6.11$ Ly$\alpha$ Emitter Behind RXCJ2248.7-4431

abstract: The CIII] and CIV rest-frame UV emission lines are powerful probes of the ionizations states of galaxies. They have furthermore been suggested as alternatives for spectroscopic redshift confirmation of objects at the epoch of reionization ($z>6$), where the most frequently used redshift indicator, Ly$\alpha$, is attenuated by the high fraction of neutral hydrogen in the inter-galactic medium. However, currently only very few confirmations of carbon UV lines at these high redshifts exist, making it challenging to quantify these claims. Here, we present the detection of CIV$\lambda\lambda$1548,1551\AA\ in \HST\ slitless grism spectroscopy obtained by GLASS of a Ly$\alpha$ emitter at $z=6.11$ multiply imaged by the massive foreground galaxy cluster RXJ2248. The CIV emission is detected at the 3--5$\sigma$ level in two images of the source, with marginal detection in two other images. We do not detect significant CIII]$\lambda\lambda$1907,1909\AA\ emission implying an equivalent width EW$_\textrm{CIII]}<20$\AA\ (1$\sigma$) and $\textrm{CIV/CIII}>0.7$ (2$\sigma$). Combined with limits on the rest-frame UV flux from the HeII$\lambda$1640\AA\ emission line and the OIII]$\lambda\lambda$1661,1666\AA\ doublet, we put constraints on the metallicity and the ionization state of the galaxy. The estimated line ratios and equivalent widths do not support a scenario where an AGN is responsible for ionizing the carbon atoms. SED fits including nebular emission lines imply a source with a mass of log(M/M$_\odot)\sim9$, SFR of around 10M$_\odot$/yr, and a young stellar population $<50$Myr old. The source shows a stronger ionizing radiation field than objects with detected CIV emission at $z<2$ and adds to the growing sample of low-mass (log(M/M$_\odot)\lesssim9$) galaxies at the epoch of reionization with strong radiation fields from star formation.


title: SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code

abstract: Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.


title: Majorana Quasi-Particles Protected by $\mathbb{Z}_2$ Angular Momentum Conservation

abstract: We show how angular momentum conservation can stabilise a symmetry-protected quasi-topological phase of matter supporting Majorana quasi-particles as edge modes in one-dimensional cold atom gases. We investigate a number-conserving four-species Hubbard model in the presence of spin-orbit coupling. The latter reduces the global spin symmetry to an angular momentum parity symmetry, which provides an extremely robust protection mechanism that does not rely on any coupling to additional reservoirs. The emergence of Majorana edge modes is elucidated using field theory techniques, and corroborated by density-matrix-renormalization-group simulations. Our results pave the way toward the observation of Majorana edge modes with alkaline-earth-like fermions in optical lattices, where all basic ingredients for our recipe - spin-orbit coupling and strong inter-orbital interactions - have been experimentally realized over the last two years.


title: The Oblique Orbit of WASP-107b from K2 Photometry

abstract: Observations of nine transits of WASP-107 during the {\it K2} mission reveal three separate occasions when the planet crossed in front of a starspot. The data confirm the stellar rotation period to be 17 days --- approximately three times the planet's orbital period --- and suggest that large spots persist for at least one full rotation. If the star had a low obliquity, at least two additional spot crossings should have been observed. They were not observed, giving evidence for a high obliquity. We use a simple geometric model to show that the obliquity is likely in the range 40-140$^\circ$, i.e., both spin-orbit alignment and anti-alignment can be ruled out. WASP-107 thereby joins the small collection of relatively low-mass stars hosting a giant planet with a high obliquity. Most such stars have been observed to have low obliquities; all the exceptions, including WASP-107, involve planets with relatively wide orbits ("warm Jupiters", with $a_{\rm min}/R_\star \gtrsim 8$). This demonstrates a connection between stellar obliquity and planet properties, in contradiction to some theories for obliquity excitation.


title: Fulde-Ferrell-Larkin-Ovchinnikov state to topological superfluidity transition in bilayer spin-orbit coupled degenerate Fermi gas

abstract: Recently a scheme has been proposed for generating the 2D Rashba-type spin-orbit coupling (SOC) for ultracold atomic bosons in a bilayer geometry [S.-W. Su et al, Phys. Rev. A \textbf{93}, 053630 (2016)]. Here we investigate the superfluidity properties of a degenerate Fermi gas affected by the SOC in such a bilayer system. We demonstrate that a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state appears in the regime of small to moderate atom-light coupling. In contrast to the ordinary SOC, the FFLO state emerges in the bilayer system without adding any external fields or spin polarization. As the atom-light coupling increases, the system can transit from the FFLO state to a topological superfluid state. These findings are also confirmed by the BdG simulations with a weak harmonic trap added.


title: Constraints on the location of the gamma-ray emission region for the gamma-ray-loud radio source GB 1310+487

abstract: We employ a single-zone leptonic jet model, with synchrotron, synchrotron self-Compton (SSC) and external Compton (EC) process, to reproduce the quasi-simultaneous multi-wavelength spectral energy distributions in active and quiescent states of the narrow-line gamma-ray-loud radio source GB 1310+487. In the case of the EC process, the external seed photons from both broad line region (BLR) and dust torus are considered by assuming that the gamma-ray emission region is located at the outside boundary of the BLR and inside the dust torus. Comparing the energy density of external photon fields $U_{\rm BLR}$ obtained by model fitting with that constrained from the BLR observations. We find that the location of the gamma-ray emitting region of GB 1310+487 can be tightly constrained at the outer edge of the BLR (the dissipation distance of the $\gamma$-ray emission region from central black hole $r_{\rm diss} \sim {\rm~a~few~times~of}~R_{\rm BLR}$). The ratio of magnetic energy and emitting-electron energy in the radiation blob ($\epsilon_{B}=L_{B}/L_{\rm e}$) is gradually increased from Flare 1, Flare 2 to Post-Flare, where the magnetic energy increase and matter energy decrease. These results suggest that the conversion of the magnetic field and the matter (radiation electrons) energy and the location of the $\gamma$-ray emission region (or ambient photon field) may play an important role in different radiation states of GB 1310+487.


title: Approximate reversal of quantum Gaussian dynamics

abstract: Recently, there has been focus on determining the conditions under which the data processing inequality for quantum relative entropy is satisfied with approximate equality. The solution of the exact equality case is due to Petz, who showed that the quantum relative entropy between two quantum states stays the same after the action of a quantum channel if and only if there is a \textit{reversal channel} that recovers the original states after the channel acts. Furthermore, this reversal channel can be constructed explicitly and is now called the \textit{Petz recovery map}. Recent developments have shown that a variation of the Petz recovery map works well for recovery in the case of approximate equality of the data processing inequality. Our main contribution here is a proof that bosonic Gaussian states and channels possess a particular closure property, namely, that the Petz recovery map associated to a bosonic Gaussian state $\sigma$ and a bosonic Gaussian channel $\mathcal{N}$ is itself a bosonic Gaussian channel. We furthermore give an explicit construction of the Petz recovery map in this case, in terms of the mean vector and covariance matrix of the state $\sigma$ and the Gaussian specification of the channel $\mathcal{N}$.


title: Nonmaximal $\theta_{23}$ mixing at NOvA from neutrino decoherence

abstract: In study of muon neutrino disappearance at 810 km, the NOvA experiment finds flavor mixing of the atmospheric sector to deviate from maximal ($\sin^2\theta_{23} = 0.5$) by 2.6 $\sigma$. The result is in tension with the 295-km baseline measurements of T2K which are consistent with maximal mixing. We propose that $\theta_{23}$ is in fact maximal, and that the disagreement is harbinger of environmentally-induced decoherence. The departure from maximal mixing can be accounted for by an energy-independent decoherence of strength $\Gamma = (2.3 \pm 1.1) \times 10^{-23}$ GeV.


title: An Efficient Parallel Data Clustering Algorithm Using Isoperimetric Number of Trees

abstract: We propose a parallel graph-based data clustering algorithm using CUDA GPU, based on exact clustering of the minimum spanning tree in terms of a minimum isoperimetric criteria. We also provide a comparative performance analysis of our algorithm with other related ones which demonstrates the general superiority of this parallel algorithm over other competing algorithms in terms of accuracy and speed.


title: 4d N=1 from 6d N=(1,0) on a torus with fluxes

abstract: Compactifying N=(1,0) theories on a torus, with additional fluxes for global symmetries, we obtain N=1 supersymmetric theories in four dimensions. It is shown that for many choices of flux these models are toric quiver gauge theories with singlet fields. In particular we compare the anomalies deduced from the description of the six-dimensional theory and the anomalies of the quiver gauge theories. We also give predictions for anomalies of four-dimensional theories corresponding to general compactifications of M5-branes probing C_2/Z_k singularities.


title: Variational formulation of the earth's elastic-gravitational deformations under low regularity conditions

abstract: We present a construction of the action, in the framework of the calculus of variations and Sobolev spaces, describing deformations and the oscillations of a uniformly rotating, elastic and self-gravitating earth. We establish the Fr\'{e}chet differentiability of the action under minimal regularity assumptions, which constrain the possible composition of an earth model. Thus we obtain well-defined Euler-Lagrange equations, weakly and strongly, that is, the system of elastic-gravitational equations.


title: Quantitative uniqueness of solutions to second order elliptic equations with singular lower order terms

abstract: In this article, we study the quantitative uniqueness of solutions to second order elliptic equations with singular lower order terms. We quantify the strong unique continuation property by estimating the maximal vanishing order of solutions. That is, when $u$ is a non-trivial solution to $\triangle u + W \cdot \nabla u + V u = 0$ in some open, connected subset of $\mathbb R^n$, where $n \geq 3$, we characterize the vanishing order of solutions in terms of the norms of $V$ and $W$ in their respective Lebesgue spaces. Using these maximal order of vanishing estimates, we also establish quantitative unique continuation at infinity results for solutions to $\triangle u + W \cdot \nabla u + V u = 0$ in $\mathbb R^n$. The main tools in our work are new versions of $L^p\to L^q$ Carleman estimates for a range of $p$- and $q$-values.


title: Theory of charge density wave depinning by electromechanical effect

abstract: We discuss the first theory for the depinning of low dimensional, incommensurate, charge density waves (CDWs) in the strong electron-phonon (e-p) regime. Arguing that most real CDWs systems invariably develop a gigantic dielectric constant (GDC) at very low frequencies, we propose an electromechanical mechanism which is based on a local field effect. At zero electric field and large enough e-p coupling the structures are naturally pinned by the lattice due to its discreteness, and develop modulation functions which are characterized by discontinuities. When the electric field is turned on, we show that it exists a finite threshold value for the electric field above which the discontinuities of the modulation functions vanish due to CDW deformation. The CDW is then free to move. The signature of this pinning/depinning transition as a function of the increasing electric field can be directly observed in the phonon spectrum by using inelastic neutrons or X-rays experiments.


title: Thermochemical modelling of brown dwarf discs

abstract: The physical properties of brown dwarf discs, in terms of their shapes and sizes, are still largely unexplored by observations. To what extent brown dwarf discs are similar to scaled-down T Tauri discs is currently unknown, and this work is a step towards establishing a relationship through the eventual modelling of future observations. We use observations of the brown dwarf disc $\rho$ Oph 102 to infer a fiducial model around which we build a small grid of brown dwarf disc models, in order to model the CO, HCN, and HCO+ line fluxes and the chemistry which drives their abundances. These are the first brown dwarf models to be published which relate detailed, 2D radiation thermochemical disc models to observational data. We predict that moderately extended ALMA antenna configurations will spatially resolve CO line emission around brown dwarf discs, and that HCN and HCO+ will be detectable in integrated flux, following our conclusion that the flux ratios of these molecules to CO emission are comparable to that of T Tauri discs. These molecules have not yet been observed in sub-mm wavelengths in a brown dwarf disc, yet they are crucial tracers of the warm surface-layer gas and of ionization in the outer parts of the disc. We present the prediction that if the physical and chemical processes in brown dwarf discs are similar to those that occur in T Tauri discs -- as our models suggest -- then the same diagnostics that are used for T Tauri discs can be used for brown dwarf discs (such as HCN and HCO+ lines that have not yet been observed in the sub-mm), and that these lines should be observable with ALMA. Through future observations, either confirmation (or refutation) of these ideas about brown dwarf disc chemistry will have strong implications for our understanding of disc chemistry, structure, and subsequent planet formation in brown dwarf discs.


title: On The Limiting Distributions of the Total Height On Families of Trees

abstract: A symbolic-computational algorithm, fully implemented in Maple, is described, that computes explicit expressions for generating functions that enable the efficient computations of the expectation, variance, and higher moments, of the random variable `sum of distances to the root', defined on any given family of rooted ordered trees (defined by degree restrictions). Taking limits, we confirm, via elementary methods, the fact, due to David Aldous, and expanded by Svante Janson and others, that the limiting (scaled) distributions are all the same, and coincide with the limiting distribution of the same random variable, when it is defined on labeled rooted trees.


title: Spectral Algorithms for Temporal Graph Cuts

abstract: The sparsest cut problem consists of identifying a small set of edges that breaks the graph into balanced sets of vertices. The normalized cut problem balances the total degree, instead of the size, of the resulting sets. Applications of graph cuts include community detection and computer vision. However, cut problems were originally proposed for static graphs, an assumption that does not hold in many modern applications where graphs are highly dynamic. In this paper, we introduce the sparsest and normalized cut problems in temporal graphs, which generalize their standard definitions by enforcing the smoothness of cuts over time. We propose novel formulations and algorithms for computing temporal cuts using spectral graph theory, multiplex graphs, divide-and-conquer and low-rank matrix approximation. Furthermore, we extend our formulation to dynamic graph signals, where cuts also capture node values, as graph wavelets. Experiments show that our solutions are accurate and scalable, enabling the discovery of dynamic communities and the analysis of dynamic graph processes.


title: Variational principle for the parallel-symplectic representation of electromagnetic gyrokinetic theory

abstract: The nonlinear (full-$f$) electromagnetic gyrokinetic Vlasov-Maxwell equations are derived in the parallel-symplectic representation from an Eulerian gyrokinetic variational principle. The gyrokinetic Vlasov-Maxwell equations are shown to possess an exact energy conservation law, which is derived by Noether method from the gyrokinetic variational principle. Here, the gyrocenter Poisson bracket and the gyrocenter Jacobian contain contributions from the perturbed magnetic field. In the full-$f$ formulation of the gyrokinetic Vlasov-Maxwell theory presented here, the gyrocenter parallel-Amp\`{e}re equation contains a second-order contribution to the gyrocenter current density that is derived from the second-order gyrocenter ponderomotive Hamiltonian.


title: An Improved Dictatorship Test with Perfect Completeness

abstract: A Boolean function $f:\{0,1\}^n\rightarrow \{0,1\}$ is called a dictator if it depends on exactly one variable i.e $f(x_1, x_2, \ldots, x_n) = x_i$ for some $i\in [n]$. In this work, we study a $k$-query dictatorship test. Dictatorship tests are central in proving many hardness results for constraint satisfaction problems. The dictatorship test is said to have {\em perfect completeness} if it accepts any dictator function. The {\em soundness} of a test is the maximum probability with which it accepts any function far from a dictator. Our main result is a $k$-query dictatorship test with perfect completeness and soundness $ \frac{2k + 1}{2^k}$, where $k$ is of the form $2^t -1$ for any integer $t > 2$. This improves upon the result of \cite{TY15} which gave a dictatorship test with soundness $ \frac{2k + 3}{2^k}$.


title: Joint Routing, Scheduling and Power Control Providing Hard Deadline in Wireless Multihop Networks

abstract: We consider optimal/efficient power allocation policies in a single/multihop wireless network in the presence of hard end-to-end deadline delay constraints on the transmitted packets. Such constraints can be useful for real time voice and video. Power is consumed in only transmission of the data. We consider the case when the power used in transmission is a convex function of the data transmitted. We develop a computationally efficient online algorithm, which minimizes the average power for the single hop. We model this problem as dynamic program (DP) and obtain the optimal solution. Next, we generalize it to the multiuser, multihop scenario when there are multiple real time streams with different hard deadline constraints.


title: How do packing defects modify the cooperative motions in supercooled liquids ?

abstract: We use molecular dynamic simulations to investigate the relation between the presence of packing defects in a glass-former and the spontaneous cooperative motions called dynamic heterogeneity. For that purpose we use a simple diatomic glass-former and add a small number of larger or smaller diatomic probes. The diluted probes modify locally the packing, inducing structural defects in the liquid, while we find that the number of defects is small enough not to disturb the average structure. We find that a small packing modification around a few molecules can deeply influence the dynamics of the whole liquid, when supercooled. When we use small probe molecules, the dynamics accelerates and the dynamic heterogeneity decreases. In contrast, for large probes the dynamics slows down and the dynamic heterogeneity increases. The induced heterogeneities and transport coefficient modification increase when the temperature decreases and disappear around the onset temperature of the cage dynamics.


title: Codimension one Fano distributions on Fano manifolds

abstract: In this paper we investigate codimension one Fano distributions on Fano manifolds with Picard number one. We classify Fano distributions of maximal index on complete intersections in weighted projective spaces, Fano contact manifolds, Grassmannians of lines and their linear sections, and describe their moduli spaces. As a consequence, we obtain a classification of codimension one del Pezzo distributions on Fano manifolds with Picard number one.


title: Measurement-device-independent entanglement and randomness estimation in quantum networks

abstract: Detection of entanglement in quantum networks consisting of many parties is one of the important steps towards building quantum communication and computation networks. We consider a scenario where the measurement devices used for this certification are uncharacterised. In this case, it is well known that by using quantum states as inputs for the measurement devices it is possible to detect any entangled state (a situation known as measurement device-independent entanglement witnessing). Here we go beyond entanglement detection and provide methods to estimate the amount of entanglement in a quantum network. We also consider the task of randomness certification and show that randomness can be certified in a variety of cases, including single-partite experiments or setups using only separable states.


title: Chiral vortical effect with finite rotation, temperature, and curvature

abstract: We perform an explicit calculation of the axial current at finite rotation and temperature in curved space. We find that finite curvature and mass corrections to the chiral vortical effect satisfy a relation of the chiral gap effect, that is, a fermion mass-shift by a scalar curvature. We also point out that a product term of the angular velocity and the scalar curvature shares the same coefficient as the mixed gravitational chiral anomaly. We discuss possible applications of the curvature induced chiral vortical effect to rotating astrophysical compact objects described by the Kerr metric. Instead of direct calculation we assume that the Chern-Simons current can approximate the physical axial current. We make a proposal that the chiral vortical current from rotating compact objects could provide a novel microscopic mechanism behind the generation of collimated jets.


title: Gauge Freedom in complex holomorphic systems

abstract: The aim of this paper is to introduce and analyze a new gauge symmetry that appears in complex holomorphic systems. This symmetry allow us to project the system, using different gauge conditions, to several real systems which are connect by gauge transformations in the complex space. We prove that the space of solutions of one system is related to the other by the gauge transformation. The gauge transformations are in some cases canonical transformations. However, in other cases are more general transformations that change the symplectic structure, but there is still a map between the systems. In this way our construction extend the group of canonical transformations in classical mechanics. Also, we show how to extend the analysis to the quantum case using path integrals by means of the Batalin-Fradkin-Vilkovisky theorem and within the canonical formalism, where we show explicitly that solutions of the Schr\"odinger equation are gauge related.


title: Estimating Individualized Treatment Rules for Ordinal Treatments

abstract: Precision medicine is an emerging scientific topic for disease treatment and prevention that takes into account individual patient characteristics. It is an important direction for clinical research, and many statistical methods have been recently proposed. One of the primary goals of precision medicine is to obtain an optimal individual treatment rule (ITR), which can help make decisions on treatment selection according to each patient's specific characteristics. Recently, outcome weighted learning (OWL) has been proposed to estimate such an optimal ITR in a binary treatment setting by maximizing the expected clinical outcome. However, for ordinal treatment settings, such as individualized dose finding, it is unclear how to use OWL. In this paper, we propose a new technique for estimating ITR with ordinal treatments. In particular, we propose a data duplication technique with a piecewise convex loss function. We establish Fisher consistency for the resulting estimated ITR under certain conditions, and obtain the convergence and risk bound properties. Simulated examples and two applications to datasets from an irritable bowel problem and a type 2 diabetes mellitus observational study demonstrate the highly competitive performance of the proposed method compared to existing alternatives.


title: Demonstrating an in-situ topological band transition in cylindrical granular chains

abstract: We numerically investigate and experimentally demonstrate an in-situ topological band transition in a highly tunable mechanical system made of cylindrical granular particles. This system allows us to tune its inter-particle stiffness in a controllable way, simply by changing the contact angles between the cylinders. The spatial variation of particles' stiffness results in an in-situ transition of the system's topology. This manifests as the emergence of a boundary mode in the finite system, which we observe experimentally via laser Doppler vibrometry. When two topologically different systems are placed adjacently, we analytically predict and computationally and experimentally demonstrate the existence of a finite-frequency topologically protected mode at their interface.


title: Determining the finite subgraphs of curve graphs

abstract: We prove that there is an algorithm to determine if a given finite graph is an induced subgraph of a given curve graph.


title: From nothing to something II: nonlinear systems via consistent correlated bang

abstract: Chinese ancient sage Laozi said everything comes from \emph{\bf \em "nothing"}. \rm In the first letter (Chin. Phys. Lett. 30 (2013) 080202), infinitely many discrete integrable systems have been obtained from "nothing" via simple principles (Dao). In this second letter, a new idea, the consistent correlated bang, is introduced to obtain nonlinear dynamic systems including some integrable ones such as the continuous nonlinear Schr\"odinger equation (NLS), the (potential) Korteweg de Vries (KdV) equation, the (potential) Kadomtsev-Petviashvili (KP) equation and the sine-Gordon (sG) equation. These nonlinear systems are derived from nothing via suitable "Dao", the shifted parity, the charge conjugate, the delayed time reversal, the shifted exchange, the shifted-parity-rotation and so on.


title: Radiative Two-Loop Neutrino Masses with Dark Matter

abstract: Using the Weinberg operator, we present a full collection of genuine two-loop models for neutrino mass generation, which contain a dark matter particle as one of the internal messengers. These models can be constructed simply by adding new fields that are singlets or doublets of SU(2)L. We ensure the stability of the dark matter candidate by the addition of a Z2 symmetry that will also be used to forbid tree level or one-loop diagrams. Thus we only present models where the main contribution for neutrinos masses is generated from the corresponding two-loop diagram. We also discuss a short outline corresponding to some phenomenological characteristics of these models.


title: Digit Preserving Multiplication in Continued Fraction Representations

abstract: A permutiple is a number which is an integer multiple of some permutation of its digits. A well-known example is 9801 since it is an integer multiple of its reversal, 1089. In this paper, we consider the permutiple problem in an entirely different setting: continued fractions. We pose the question of when the simple continued fraction representation of a rational number is an integer multiple of a permutation of its partial quotients (or digits, as we shall call them). We develop some general results and apply them to finding new examples. In doing so, we attempt to classify all 2, 3, and 4-digit continued fraction permutiples in terms of basic permutiple types which we discover along the way. We also generate new examples from old by finding conditions which guarantee that digit-string concatenation yields other permutiples.


title: Spaser quenching by off-resonant plasmon modes

abstract: We study the effect of off-resonant plasmon modes on spaser threshold in nanoparticle-based spasers. We develop an analytical semiclassical model and derive spaser threshold condition accounting for gain coupling to higher-order plasmons. We show that such a coupling originates from inhomogeneity of gain distribution near the metal surface and leads to an upward shift of spaser frequency and population inversion threshold. This effect is similar, albeit significantly weaker, to quenching of plasmon-enhanced fluorescence near metal nanostructures due to excitation of off-resonant modes with wide spectral band. We also show that spaser quenching is suppressed for high gain concentrations and establish a simple criterion for quenching onset, which we support by numerical calculations for spherical geometry.


title: Photometric Observations of Supernova 2013cq Associated with GRB 130427A

abstract: We observed the afterglow of GRB 130427A with the RATIR instrument on the 1.5-m Harold L. Johnson telescope of the Observatorio Astron\'omico Nacional on Sierra San Pedro M\'artir. Our homogenous $griZYJH$ photometry extends from the night of burst to three years later. We fit a model for the afterglow. There is a significant positive residual which matches the behavior of SN 1998bw in the $griZ$ filters; we suggest that this is a photometric signature of the supernova SN 2013cq associated with the GRB. The peak absolute magnitude of the supernova is $M_r=-18.43\pm0.11$.


title: Semi-hyperbolic rational maps and size of Fatou components

abstract: Recently Merenkov and Sabitova introduced the notion of a homogeneous planar set. Using this notion they proved a result for Sierpi${\'n}$ski carpet Julia sets of hyperbolic rational maps that relates the diameters of the peripheral circles to the Hausdorff dimension of the Julia set. We extend this theorem to Julia sets (not necessarily Sierpi${\'n}$ski carpets) of semi-hyperbolic rational maps, and prove a stronger version of the theorem that was conjectured by Merenkov and Sabitova.


title: Optimal order Jackson type inequality for scaled Shepard approximation

abstract: We study a variation of the Shepard approximation scheme by introducing a dilation factor into the base function, which synchronizes with the Hausdorff distance between the data set and the domain. The novelty enables us to establish an optimal order Jackson \cite{jackson} type error estimate (with an explicit constant) for bounded continuous functions on any given convex domain. We also improve en route an upper bound estimate due to Narcowich and Ward for the numbers of well-separated points in thin annuli, which is of independent interest.


title: Hot magnetized nuclear matter: Thermodynamic and Saturation Properties

abstract: We have used a realistic nuclear potential, AV18, and a many body technique, the lowest order constraint variational (LOCV) approach, to calculate the properties of hot magnetized nuclear matter. By investigating the free energy, spin polarization parameter, and symmetry energy, we have studied the temperature and magnetic field dependence of the saturation properties of magnetized nuclear matter. In addition, we have calculated the equation of state of magnetized nuclear matter at different temperatures and magnetic fields. It was found that the flashing temperature of nuclear matter decreases by increasing the magnetic field. In addition, we have studied the effect of the magnetic field on liquid gas phase transition of nuclear matter. The liquid gas coexistence curves, the order parameter of the liquid gas phase transition, and the properties of critical point at different magnetic fields have been calculated.


title: Quantum dilogarithm identities for the square product of A-type Dynkin quivers

abstract: The famous pentagon identity for quantum dilogarithms has a generalization for every Dynkin quiver, due to Reineke. A more advanced generalization is associated with a pair of alternating Dynkin quivers, due to Keller. The description and proof of Keller's identities involves cluster algebras and cluster categories, and the statement of the identity is implicit. In this paper we describe Keller's identities explicitly, and prove them by a dimension counting argument. Namely, we consider quiver representations $\boldsymbol{\mathrm{Rep}}_\gamma$ together with a superpotential function $W_\gamma$, and calculate the Betti numbers of the equivariant $W_\gamma$ rapid decay cohomology algebra of $\boldsymbol{\mathrm{Rep}}_\gamma$ in two different ways corresponding to two natural stratifications of $\boldsymbol{\mathrm{Rep}}_\gamma$. This approach is suggested by Kontsevich and Soibelman in relation with the Cohomological Hall Algebra of quivers, and the associated Donaldson-Thomas invariants.


title: Linear Time Computation of Moments in Sum-Product Networks

abstract: Bayesian online algorithms for Sum-Product Networks (SPNs) need to update their posterior distribution after seeing one single additional instance. To do so, they must compute moments of the model parameters under this distribution. The best existing method for computing such moments scales quadratically in the size of the SPN, although it scales linearly for trees. This unfortunate scaling makes Bayesian online algorithms prohibitively expensive, except for small or tree-structured SPNs. We propose an optimal linear-time algorithm that works even when the SPN is a general directed acyclic graph (DAG), which significantly broadens the applicability of Bayesian online algorithms for SPNs. There are three key ingredients in the design and analysis of our algorithm: 1). For each edge in the graph, we construct a linear time reduction from the moment computation problem to a joint inference problem in SPNs. 2). Using the property that each SPN computes a multilinear polynomial, we give an efficient procedure for polynomial evaluation by differentiation without expanding the network that may contain exponentially many monomials. 3). We propose a dynamic programming method to further reduce the computation of the moments of all the edges in the graph from quadratic to linear. We demonstrate the usefulness of our linear time algorithm by applying it to develop a linear time assume density filter (ADF) for SPNs.


title: Symplectic integrators for second-order linear non-autonomous equations

abstract: Two families of symplectic methods specially designed for second-order time-dependent linear systems are presented. Both are obtained from the Magnus expansion of the corresponding first-order equation, but otherwise they differ in significant aspects. The first family is addressed to problems with low to moderate dimension, whereas the second is more appropriate when the dimension is large, in particular when the system corresponds to a linear wave equation previously discretised in space. Several numerical experiments illustrate the main features of the new schemes.


title: Monadic Second Order Logic with Measure and Category Quantifiers

abstract: We investigate the extension of Monadic Second Order logic, interpreted over infinite words and trees, with generalized "for almost all" quantifiers interpreted using the notions of Baire category and Lebesgue measure.


title: Training Language Models Using Target-Propagation

abstract: While Truncated Back-Propagation through Time (BPTT) is the most popular approach to training Recurrent Neural Networks (RNNs), it suffers from being inherently sequential (making parallelization difficult) and from truncating gradient flow between distant time-steps. We investigate whether Target Propagation (TPROP) style approaches can address these shortcomings. Unfortunately, extensive experiments suggest that TPROP generally underperforms BPTT, and we end with an analysis of this phenomenon, and suggestions for future work.


title: Remarks on the construction of worm algorithms for lattice field theories in worldline representation

abstract: We introduce a generalized worldline model where the partition function is a sum over configurations of a conserved flux on a d-dimensional lattice. The weights for the configurations of the corresponding worldlines have factors living on the links of the lattice, as well as terms which live on the sites x and depend on all fluxes attached to x. The model represents a general class of worldline systems, among them the dual representation of the relativistic Bose gas at finite density. We construct a suitable worm algorithm and show how to correctly distribute the site weights in the various Metropolis probabilities that determine the worm. We analyze the algorithm in detail and give a proof of detailed balance. Our algorithm admits the introduction of an amplitude parameter A that can be chosen freely. Using a numerical simulation of the relativistic Bose gas we demonstrate that A allows one to influence the starting and terminating probabilities and thus the average length and the efficiency of the worm.


title: The population of SNe/SNRs in the starburst galaxy Arp 220. A self-consistent analysis of 20 years of VLBI monitoring

abstract: The nearby ultra-luminous infrared galaxy (ULIRG) Arp 220 is an excellent laboratory for studies of extreme astrophysical environments. For 20 years, Very Long Baseline Interferometry (VLBI) has been used to monitor a population of compact sources thought to be supernovae (SNe), supernova remnants (SNRs) and possibly active galactic nuclei (AGNs). Using new and archival VLBI data spanning 20 years, we obtain 23 high-resolution radio images of Arp 220 at wavelengths from 18 cm to 2 cm. From model-fitting to the images we obtain estimates of flux densities and sizes of all detected sources. We detect radio continuum emission from 97 compact sources and present flux densities and sizes for all analysed observation epochs. We find evidence for a LD-relation within Arp 220, with larger sources being less luminous. We find a compact source LF $n(L)\propto L^\beta$ with $\beta=-2.19\pm0.15$, similar to SNRs in normal galaxies. Based on simulations we argue that there are many relatively large and weak sources below our detection threshold. The observations can be explained by a mixed population of SNe and SNRs, where the former expand in a dense circumstellar medium (CSM) and the latter interact with the surrounding interstellar medium (ISM). Nine sources are likely luminous, type IIn SNe. This number of luminous SNe correspond to few percent of the total number of SNe in Arp 220 which is consistent with a total SN-rate of 4 yr$^{-1}$ as inferred from the total radio emission given a normal stellar initial mass function (IMF). Based on the fitted luminosity function, we argue that emission from all compact sources, also below our detection threshold, make up at most 20\% of the total radio emission at GHz frequencies.


title: The VLT-FLAMES Tarantula Survey XXVI: Properties of the O-dwarf population in 30 Doradus

abstract: The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical techniques, we applied the Bayesian BONNSAI tool to estimate evolutionary masses. We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (vsini higher than 300 km/s ) in our sample have masses below 25 MSun and intermediate rotation-corrected gravities (log gc between 3.9 and 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below 20 MSun. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with luminosity lower than logL/LSun about 5.1. Ultraviolet spectroscopy is needed to undertake a detailed investigation of the wind properties of these dwarfs.


title: An Ando-Choi-Effros lifting theorem respecting subspaces

abstract: We prove a version of the Ando-Choi-Effros lifting theorem respecting subspaces, which in turn relies on Oja's principle of local reflexivity respecting subspaces. To achieve this, we first develop a theory of pairs of $M$-ideals. As a first consequence we get a version respecting subspaces of the Michael-Pe{\l}czy\'nski extension theorem. Other applications are related to linear and Lipschitz bounded approximation properties for a pair consisting of a Banach space and a subspace. We show that in the separable case, the BAP for such a pair is equivalent to the simultaneous splitting of an associated pair of short exact sequences given by a construction of Lusky. We define a Lipschitz version of the BAP for pairs, and study its relationship to the (linear) BAP for pairs. The two properties are not equivalent in general, but they are when the pair has an additional Lipschitz-lifting property in the style of Godefroy and Kalton. We also characterize, in the separable case, those pairs of a metric space and a subset whose corresponding pair of Lipschitz-free spaces has the BAP.


title: Bayesian Additive Adaptive Basis Tensor Product Models for Modeling High Dimensional Surfaces: An application to high-throughput toxicity testing

abstract: Many modern data sets are sampled with error from complex high-dimensional surfaces. Methods such as tensor product splines or Gaussian processes are effective/well suited for characterizing a surface in two or three dimensions but may suffer from difficulties when representing higher dimensional surfaces. Motivated by high throughput toxicity testing where observed dose-response curves are cross sections of a surface defined by a chemical's structural properties, a model is developed to characterize this surface to predict untested chemicals' dose-responses. This manuscript proposes a novel approach that models the multidimensional surface as a sum of learned basis functions formed as the tensor product of lower dimensional functions, which are themselves representable by a basis expansion learned from the data. The model is described, a Gibbs sampling algorithm proposed, and is investigated in a simulation study as well as data taken from the US EPA's ToxCast high throughput toxicity testing platform.