title: Influence of plasma background on 3D scrape-off layer filaments

abstract: This paper presents the effect of self-consistent plasma backgrounds including plasma-neutral interactions, on the dynamics of filament propagation. The principle focus is on the influence of the neutrals on the filament through both direct interactions and through their influence on the plasma background. Both direct and indirect interactions influence the motion of filaments. A monotonic increase of filament peak velocity with upstream electron temperature is observed, while a decrease with increasing electron density is observed. If ordered by the target temperature, the density dependence disappears and the filament velocity is only a function of the target temperature. Smaller filaments keep a density dependence, as a result of the density dependence of the plasma viscosity. The critical size $\delta^*$, where filaments are fastest, is shifted to larger sizes for higher densities, due to the plasma viscosity. If the density dependence of the plasma viscosity is removed, $\delta^*$ has no temperature dependence, but rather a density dependence.


title: Simultaneous Signal Subspace Rank and Model Selection with an Application to Single-snapshot Source Localization

abstract: This paper proposes a novel method for model selection in linear regression by utilizing the solution path of $\ell_1$ regularized least-squares (LS) approach (i.e., Lasso). This method applies the complex-valued least angle regression and shrinkage (c-LARS) algorithm coupled with a generalized information criterion (GIC) and referred to as the c-LARS-GIC method. c-LARS-GIC is a two-stage procedure, where firstly precise values of the regularization parameter, called knots, at which a new predictor variable enters (or leaves) the active set are computed in the Lasso solution path. Active sets provide a nested sequence of regression models and GIC then selects the best model. The sparsity order of the chosen model serves as an estimate of the model order and the LS fit based only on the active set of the model provides an estimate of the regression parameter vector. We then consider a source localization problem, where the aim is to detect the number of impinging source waveforms at a sensor array as well to estimate their direction-of-arrivals (DoA-s) using only a single-snapshot measurement. We illustrate via simulations that, after formulating the problem as a grid-based sparse signal reconstruction problem, the proposed c-LARS-GIC method detects the number of sources with high probability while at the same time it provides accurate estimates of source locations.


title: Antiferromagnetic-ferromagnetic phase domain development in nanopatterned FeRh islands

abstract: The antiferromagnetic to ferromagnetic phase transition in B2-ordered FeRh is imaged in laterally confined nanopatterned islands using photoemission electron microscopy with x-ray magnetic circular dichroism contrast. The resulting magnetic images directly detail the progression in the shape and size of the FM phase domains during heating and cooling through the transition. In 5 um square islands this domain development during heating is shown to proceed in three distinct modes: nucleation, growth, and merging, each with subsequently greater energy costs. In 0.5 um islands, which are smaller than the typical final domain size, the growth mode is stunted and the transition temperature was found to be reduced by 20 K. The modification to the transition temperature is found by high resolution scanning transmission electron microscopy to be due to a 100 nm chemically disordered edge grain present as a result of ion implantation damage during the patterning. FeRh has unique possibilities for magnetic memory applications; the inevitable changes to its magnetic properties due to subtractive nanofabrication will need to be addressed in future work in order to progress from sheet films to suitable patterned devices.


title: Magnetothermal disk winds in X-ray binaries: poloidal magnetic fields suppress thermal winds

abstract: Magnetic, radiation pressure, and thermal driving are the three mechanisms capable of launching accretion disk winds. In X-ray binaries, radiation pressure is often not significant, as in many systems the luminosity is too low for driving due to continuum transitions yet too high for driving due to line transitions. This leaves thermal and magnetic driving as the contender launching mechanisms in these systems. Using Athena++, we perform axisymmetric ideal MHD simulations that include radiative heating and cooling processes appropriate for Compton heated winds to show that the inclusion of magnetic fields into a thermally driven wind has the opposite effect of what one might expect: rather than provide a velocity boost, the thermal wind is suppressed in low plasma beta regions where the field lines are strong enough to reshape the direction of the flow. Our analysis reveals that magneto-centrifugal launching is present but weak, while the reduction in wind velocity is not due to the change in gravitational potential through the magnetically imposed streamline geometry, but rather due to the increased flow tube area just above the surface of the disk, which is less conducive to acceleration. Our results suggest that for magnetothermal wind models to be successful at producing fast dense outflows in low mass X-ray binaries, the winds must be magnetically launched well within the Compton radius.


title: Asymptotic behavior of solutions to the generalized KdV-Burgers equation with slowly decaying data

abstract: We consider the asymptotic behavior of the global solutions to the initial value problem for the generalized KdV-Burgers equation. It is known that the solution to this problem converges to a self-similar solution to the Burgers equation called a nonlinear diffusion wave. In this paper, we derive the optimal asymptotic rate to the nonlinear diffusion wave when the initial data decays slowly at spatial infinity. In particular, we investigate that how the change of the decay rate of the initial value affects the asymptotic rate to the nonlinear diffusion wave.


title: Multi-agent Gaussian Process Motion Planning via Probabilistic Inference

abstract: This paper deals with motion planning for multiple agents by representing the problem as a simultaneous optimization of every agent's trajectory. Each trajectory is considered as a sample from a one-dimensional continuous-time Gaussian process (GP) generated by a linear time-varying stochastic differential equation driven by white noise. By formulating the planning problem as probabilistic inference on a factor graph, the structure of the pertaining GP can be exploited to find the solution efficiently using numerical optimization. In contrast to planning each agent's trajectory individually, where only the current poses of other agents are taken into account, we propose simultaneous planning of multiple trajectories that works in a predictive manner. It takes into account the information about each agent's whereabouts at every future time instant, since full trajectories of each agent are found jointly during a single optimization procedure. We compare the proposed method to an individual trajectory planning approach, demonstrating significant improvement in both success rate and computational efficiency.


title: Partial regularity for manifold constrained p(x)-harmonic maps

abstract: We prove that manifold constrained $p(x)$-harmonic maps are $C^{1,\beta}$-regular outside a set of zero $n$-dimensional Lebesgue's measure, for some $\beta \in (0,1)$. We also provide an estimate from above of the Hausdorff dimension of the singular set.


title: Simulations of photo-nuclear dijets with Pythia 8 and their sensitivity to nuclear PDFs

abstract: In ultra-peripheral heavy-ion collisions the nuclei pass with such a large impact parameter that no strong interaction can take place. However, a photon emitted by one nucleus can still interact either with the other nucleus or with a photon from the another nucleus. The former event type corresponds to photoproduction, earlier studied e.g. in electron-proton collisions at HERA, but now with a nuclear target. This provides an unique opportunity to study the nuclear modifications of the PDFs with the photo-nuclear processes measured in lead-lead collisions at the LHC. Here we present our recent implementation of photoproduction processes in the Pythia 8 general-purpose Monte-Carlo event generator. The resulting simulations are compared to charged-hadron and dijet photoproduction data from the HERA experiments, and the theoretical uncertainties of the framework are studied. We discuss how the relevant part of the photon flux from heavy ions is obtained using the equivalent photon approximation and present predictions for photo-nuclear dijet cross section at the LHC. To estimate the potential of this data to further constrain the nuclear PDFs, the expected statistical uncertainty based on LHC luminosity is compared to uncertainties in the current nuclear PDF analyses.


title: Active attenuation of a trailing vortex inspired by a parabolized stability analysis

abstract: Designing effective control for complex three-dimensional flow fields proves to be non-trivial. Oftentimes, intuitive control strategies lead to suboptimal control. To navigate the control space, we utilize a linear parabolized stability analysis to guide the design of a control scheme for a trailing vortex flow field aft of a NACA0012 half-wing at an angle of attack $\alpha = 5^\circ$ and a chord-based Reynolds number $Re = 1000$. The stability results show that the unstable mode with the smallest growth rate (fifth wake mode) provides a pathway to excite a vortex instability, whereas the principal unstable mode remains in the wake of the wing. Inspired by this finding, we perform direct numerical simulations that excite each mode with body forces matching the shape function from the stability analysis. Relative to the baseline uncontrolled case, the principal wake mode reduces the vortex length, while the fifth wake mode further shortens the tip vortex. Analogously, the streamwise circulation of the trailing vortex is found to be significantly reduced. From these results, we conclude that a rudimentary linear stability analysis can provide key insight into the underlying physics and help engineers design more effective control.


title: Hybrid Planning and Control for Multiple Fixed-Wing Aircraft under Input Constraints

abstract: This paper presents a novel hybrid control protocol for de-conflicting multiple vehicles with constraints on control inputs. We consider turning rate and linear speed constraints to represent fixed-wing or car-like vehicles. A set of state-feedback controllers along with a state-dependent switching logic are synthesized in a hybrid system to generate collision-free trajectories that converge to the desired destinations of the vehicles. The switching law is designed so that the safety can be guaranteed while no Zeno behavior can occur. A novel temporary goal assignment technique is also designed to guarantee convergence. We analyze the individual modes for safety and the closed-loop hybrid system for convergence. The theoretical developments are demonstrated via simulation results.


title: Total-variation methods for gravitational-wave denoising: performance tests on Advanced LIGO data

abstract: We assess total-variation methods to denoise gravitational-wave signals in real noise conditions, by injecting numerical-relativity waveforms from core-collapse supernovae and binary black hole mergers in data from the first observing run of Advanced LIGO. This work is an extension of our previous investigation where only Gaussian noise was used. Since the quality of the results depends on the regularization parameter of the model, we perform an heuristic search for the value that produces the best results. We discuss various approaches for the selection of this parameter, either based on the optimal, mean, or multiple values, and compare the results of the denoising upon these choices. Moreover, we also present a machine-learning-informed approach to obtain the Lagrange multiplier of the method through an automatic search. Our results provide further evidence that total-variation methods can be useful in the field of Gravitational-Wave Astronomy as a tool to remove noise.


title: Rolling Skyrmions and the Nuclear Spin-Orbit Force

abstract: We compute the nuclear spin-orbit coupling from the Skyrme model. Previous attempts to do this were based on the product ansatz, and as such were limited to a system of two well-separated nuclei. Our calculation utilises a new method, and is applicable to the phenomenologically important situation of a single nucleon orbiting a large nucleus. We find that, to second order in perturbation theory, the coefficient of the spin-orbit coupling induced by pion field interactions has the wrong sign, but as the strength of the pion-nucleon interactions increases the correct sign is recovered non-perturbatively.


title: Motion in time-periodic backgrounds with applications to ultralight dark matter haloes at galactic centers

abstract: We consider motion in spherically symmetric but time-dependent backgrounds. This problem is of interest, for example, in the context of ultralight dark matter, where galactic haloes produce a time-dependent and periodic gravitational potential. We study the properties of motion of stars in such spacetimes, for different field strengths and frequency, and including dissipative effects. We show that orbital resonances may occur and that spectroscopic emission lines from stars in these geometries exhibit characteristic, periodic modulation patterns. In addition, we work out a fully relativistic and weak-field description of a special class of time-periodic geometries, that of scalar oscillatons. When applied to the galactic center, our results indicate that the motion of S2-like stars may carry distinguishable observational imprints of ultra-light dark matter.

categories: gr-qc astro-ph.CO astro-ph.GA hep-ph physics.class-ph


title: Quantifying Operations with an Application to Coherence

abstract: To describe certain facets of non-classicality, it is necessary to quantify properties of operations instead of states. This is the case if one wants to quantify how well an operation detects non-classicality, which is a necessary prerequisite for its use in quantum technologies. To do so rigorously, we build resource theories on the level of operations, exploiting the concept of resource destroying maps. We discuss the two basic ingredients of these resource theories, the free operations and the free super-operations, which are sequential and parallel concatenations with free operations. This leads to defining properties of functionals that are well suited to quantify the resources of operations. We introduce these concepts at the example of coherence. In particular, we present two measures quantifying the ability of an operation to detect, i.e. to use, coherence, one of them with an operational interpretation, and provide methods to evaluate them.


title: Rate-Memory Trade-Off for Caching and Delivery of Correlated Sources

abstract: This paper studies the fundamental limits of content delivery in a cache-aided broadcast network for correlated content generated by a discrete memoryless source with arbitrary joint distribution. Each receiver is equipped with a cache of equal capacity, and the requested files are delivered over a shared error-free broadcast link. A class of achievable correlation-aware schemes based on a two-step source coding approach is proposed. Library files are first compressed, and then cached and delivered using a combination of correlation-unaware multiple-request cache-aided coded multicast schemes. The first step uses Gray-Wyner source coding to represent the library via private descriptions and descriptions that are common to more than one file. The second step then becomes a multiple-request caching problem, where the demand structure is dictated by the configuration of the compressed library, and it is interesting in its own right. The performance of the proposed two-step scheme is evaluated by comparing its achievable rate with a lower bound on the optimal peak and average rate-memory tradeoffs in a two-file multiple-receiver network, and in a three-file two-receiver network. Specifically, in a network with two files and two receivers, the achievable rate matches the lower bound for a significant memory regime and it is within half of the conditional entropy of files for all other memory values. In the three-file two-receiver network, the two-step strategy achieves the lower bound for large cache capacities, and it is within half of the joint entropy of two of the sources conditioned on the third one for all other cache sizes.


title: Movement-efficient Sensor Deployment in Wireless Sensor Networks with Limited Communication Range

abstract: We study a mobile wireless sensor network (MWSN) consisting of multiple mobile sensors or robots. Three key factors in MWSNs, sensing quality, energy consumption, and connectivity, have attracted plenty of attention, but the interaction of these factors is not well studied. To take all the three factors into consideration, we model the sensor deployment problem as a constrained source coding problem. %, which can be applied to different coverage tasks, such as area coverage, target coverage, and barrier coverage. Our goal is to find an optimal sensor deployment (or relocation) to optimize the sensing quality with a limited communication range and a specific network lifetime constraint. We derive necessary conditions for the optimal sensor deployment in both homogeneous and heterogeneous MWSNs. According to our derivation, some sensors are idle in the optimal deployment of heterogeneous MWSNs. Using these necessary conditions, we design both centralized and distributed algorithms to provide a flexible and explicit trade-off between sensing uncertainty and network lifetime. The proposed algorithms are successfully extended to more applications, such as area coverage and target coverage, via properly selected density functions. Simulation results show that our algorithms outperform the existing relocation algorithms.


title: $C^{1,\alpha}$ isometric extensions

abstract: In this paper we consider the Cauchy problem for isometric immersions. More precisely, given a smooth isometric immersion of a codimension one submanifold we construct $C^{1,\alpha}$ isometric extensions for any $\alpha<\frac{1}{n(n+1)+1}$ via the method of convex integration.


title: Neural Code Comprehension: A Learnable Representation of Code Semantics

abstract: With the recent success of embeddings in natural language processing, research has been conducted into applying similar methods to code analysis. Most works attempt to process the code directly or use a syntactic tree representation, treating it like sentences written in a natural language. However, none of the existing methods are sufficient to comprehend program semantics robustly, due to structural features such as function calls, branching, and interchangeable order of statements. In this paper, we propose a novel processing technique to learn code semantics, and apply it to a variety of program analysis tasks. In particular, we stipulate that a robust distributional hypothesis of code applies to both human- and machine-generated programs. Following this hypothesis, we define an embedding space, inst2vec, based on an Intermediate Representation (IR) of the code that is independent of the source programming language. We provide a novel definition of contextual flow for this IR, leveraging both the underlying data- and control-flow of the program. We then analyze the embeddings qualitatively using analogies and clustering, and evaluate the learned representation on three different high-level tasks. We show that even without fine-tuning, a single RNN architecture and fixed inst2vec embeddings outperform specialized approaches for performance prediction (compute device mapping, optimal thread coarsening); and algorithm classification from raw code (104 classes), where we set a new state-of-the-art.


title: Optimal measurements for quantum multiparameter estimation with general states

abstract: We generalize the approach by Braunstein and Caves [Phys. Rev. Lett. 72, 3439 (1994)] to quantum multi-parameter estimation with general states. We derive a matrix bound of the classical Fisher information matrix due to each measurement operator. The saturation of all these bounds results in the saturation of the matrix Helstrom Cram\'er-Rao bound. Remarkably, the saturation of the matrix bound is equivalent to the saturation of the scalar bound with respect to any given positive definite weight matrix. Necessary and sufficient conditions are obtained for the optimal measurements that give rise to the Helstrom Cram\'er-Rao bound associated with a general quantum state. To saturate the Helstrom bound with separable measurements or collective measurement entangling only a small number of identical states, we find it is necessary for the symmetric logarithmic derivatives to commute on the support of the state. As an important application of our results, we construct several local optimal measurements for the problem of estimating the three-dimensional separation of two incoherent optical point sources.


title: It Takes a Village: Documenting the Contributions of Non-Scientific Staff to Scientific Research

abstract: Science, especially large-scale basic research, is a collaborative endeavor, often drawing on the skills of people from a wide variety of disciplines. These people include not just scientists, but also administrators, engineers, and many others. Fermilab, a Department of Energy National Laboratory and the United States' premier particle physics laboratory, exemplifies this kind of research; many of its high-energy physics experiments involve hundreds of collaborators from all over the world. The Fermilab Archives seeks to document the history of the lab and the unique scientific research its staff and visitors perform. Adequately documenting the lab's work often requires us to go far beyond things like the writings and correspondence of scientists to also capture the administrative and social histories of the experiments and the context in which they were performed. At Fermilab, we have sought to capture these elements of the lab's activities through an oral history program that focuses on support staff as well as physicists and collection development choices that recognize the importance of records documenting the cultural life of the lab. These materials are not merely supplementary, but rather essential documentation of the many types of labor that go into the planning and execution of an experiment or the construction of an accelerator and the context in which this work is performed.


title: Topologically robust zero-sum games and Pfaffian orientation -- How network topology determines the long-time dynamics of the antisymmetric Lotka-Volterra equation

abstract: To explore how the topology of interaction networks determines the robustness of dynamical systems, we study the antisymmetric Lotka-Volterra equation (ALVE). The ALVE is the replicator equation of zero-sum games in evolutionary game theory, in which the strengths of pairwise interactions between strategies are defined by an antisymmetric matrix such that typically some strategies go extinct over time. Here we show that there also exist topologically robust zero-sum games, such as the rock-paper-scissors game, for which all strategies coexist for all choices of interaction strengths. We refer to such zero-sum games as coexistence networks and construct coexistence networks with an arbitrary number of strategies. By mapping the long-time dynamics of the ALVE to the algebra of antisymmetric matrices, we identify simple graph-theoretical rules by which coexistence networks are constructed. Examples are triangulations of cycles characterized by the golden ratio $\varphi = 1.6180...$, cycles with complete subnetworks, and non-Hamiltonian networks. In graph-theoretical terms, we extend the concept of a Pfaffian orientation from even-sized to odd-sized networks. Our results show that the topology of interaction networks alone can determine the long-time behavior of nonlinear dynamical systems, and may help to identify robust network motifs arising, for example, in ecology.


title: Radiative instability of quantum electrodynamics in chiral matter

abstract: Modification of the photon dispersion relation in chiral matter enables $1\to 2$ scattering. As a result, the single fermion and photon states are unstable to photon radiation and pair production respectively. In particular, a fast fermion moving through chiral matter can spontaneously radiate a photon, while a photon can spontaneously radiate a fast fermion and anti-fermion pair. The corresponding spectra are derived in the ultra-relativistic approximation. It is shown that the polarization of the produced and decayed photons is determined by the sign of the chiral conductivity. Impact of a flat thin domain wall on the spectra is computed.


title: Improving Protein Gamma-Turn Prediction Using Inception Capsule Networks

abstract: Protein gamma-turn prediction is useful in protein function studies and experimental design. Several methods for gamma-turn prediction have been developed, but the results were unsatisfactory with Matthew correlation coefficients (MCC) around 0.2-0.4. One reason for the low prediction accuracy is the limited capacity of the methods; in particular, the traditional machine-learning methods like SVM may not extract high-level features well to distinguish between turn or non-turn. Hence, it is worthwhile exploring new machine-learning methods for the prediction. A cutting-edge deep neural network, named Capsule Network (CapsuleNet), provides a new opportunity for gamma-turn prediction. Even when the number of input samples is relatively small, the capsules from CapsuleNet are very effective to extract high-level features for classification tasks. Here, we propose a deep inception capsule network for gamma-turn prediction. Its performance on the gamma-turn benchmark GT320 achieved an MCC of 0.45, which significantly outperformed the previous best method with an MCC of 0.38. This is the first gamma-turn prediction method utilizing deep neural networks. Also, to our knowledge, it is the first published bioinformatics application utilizing capsule network, which will provides a useful example for the community.


title: A Reputation System for Artificial Societies

abstract: One approach to achieving artificial general intelligence (AGI) is through the emergence of complex structures and dynamic properties arising from decentralized networks of interacting artificial intelligence (AI) agents. Understanding the principles of consensus in societies and finding ways to make consensus more reliable becomes critically important as connectivity and interaction speed increase in modern distributed systems of hybrid collective intelligences, which include both humans and computer systems. We propose a new form of reputation-based consensus with greater resistance to reputation gaming than current systems have. We discuss options for its implementation, and provide initial practical results.


title: Quantum Nash equilibrium in the thermodynamic limit

abstract: The quantum Nash equilibrium in the thermodynamic limit is studied for games like quantum Prisoner's dilemma and the quantum game of chicken. A phase transition is seen in both games as a function of the entanglement in the game. We observe that for maximal entanglement irrespective of the classical payoffs, a majority of players choose Quantum strategy over Defect in the thermodynamic limit.

categories: quant-ph cond-mat.stat-mech cs.GT econ.EM physics.soc-ph


title: Reducing Property Graph Queries to Relational Algebra for Incremental View Maintenance

abstract: The property graph data model of modern graph database systems is increasingly adapted for storing and processing heterogeneous datasets like networks. Many challenging applications with near real-time requirements -- e.g. financial fraud detection, recommendation systems, and on-the-fly validation -- can be captured with graph queries, which are evaluated repeatedly. To ensure quick response time for a changing data set, these applications would benefit from applying incremental view maintenance (IVM) techniques, which can perform continuous evaluation of queries and calculate the changes in the result set upon updates. However, currently, no graph databases provide support for incremental views. While IVM problems have been studied extensively over relational databases, views on property graph queries require operators outside the scope of standard relational algebra. Hence, tackling this problem requires the integration of numerous existing IVM techniques and possibly further extensions. In this paper, we present an approach to perform IVM on property graphs, using a nested relational algebraic representation for property graphs and graph operations. Then we define a chain of transformations to reduce most property graph queries to flat relational algebra and use techniques from discrimination networks (used in rule-based expert systems) to evaluate them. We demonstrate the approach using our prototype tool, ingraph, which uses openCypher, an open graph query language specified as part of an industry initiative. However, several aspects of our approach can be generalised to other graph query languages such as G-CORE and PGQL.


title: Hydrodynamic Phonon Transport Perpendicular to Diffuse-Gray Boundaries

abstract: In this paper, we examine the application of an ideal phonon-hydrodynamic material as the heat transfer medium between two non-hydrodynamic contacts with a finite temperature difference. We use the integral-equation approach to solve a modified phonon Boltzmann transport equation with the displaced Bose-Einstein distribution as the equilibrium distribution between two boundaries perpendicular to the heat transfer direction. When the distance between the boundaries is smaller than the phonon normal scattering mean free path, our solution converges to the ballistic limit as expected. In the other limit, we find that, although the local thermal conductivity in the bulk of the hydrodynamic material approaches infinity, the thermal boundary resistance at the hydrodynamic/non-hydrodynamic interfaces becomes dominant. Our study provides insights to both the steady-state thermal characterization of phonon-hydrodynamic materials and the practical application of phonon-hydrodynamic materials for thermal management.


title: A Scalable Machine Learning Approach for Inferring Probabilistic US-LI-RADS Categorization

abstract: We propose a scalable computerized approach for large-scale inference of Liver Imaging Reporting and Data System (LI-RADS) final assessment categories in narrative ultrasound (US) reports. Although our model was trained on reports created using a LI-RADS template, it was also able to infer LI-RADS scoring for unstructured reports that were created before the LI-RADS guidelines were established. No human-labelled data was required in any step of this study; for training, LI-RADS scores were automatically extracted from those reports that contained structured LI-RADS scores, and it translated the derived knowledge to reasoning on unstructured radiology reports. By providing automated LI-RADS categorization, our approach may enable standardizing screening recommendations and treatment planning of patients at risk for hepatocellular carcinoma, and it may facilitate AI-based healthcare research with US images by offering large scale text mining and data gathering opportunities from standard hospital clinical data repositories.


title: Couplings for determinantal point processes and their reduced Palm distributions with a view to quantifying repulsiveness

abstract: For a determinantal point process $X$ with a kernel $K$ whose spectrum is strictly less than one, Andr{\'e} Goldman has established a coupling to its reduced Palm process $X^u$ at a point $u$ with $K(u,u)>0$ so that almost surely $X^u$ is obtained by removing a finite number of points from $X$. We sharpen this result, assuming weaker conditions and establishing that $X^u$ can be obtained by removing at most one point from $X$, where we specify the distribution of the difference $\xi_u:=X\setminus X^u$. This is used for discussing the degree of repulsiveness in DPPs in terms of $\xi_u$, including Ginibre point processes and other specific parametric models for DPPs.


title: Statistical Optimal Transport via Factored Couplings

abstract: We propose a new method to estimate Wasserstein distances and optimal transport plans between two probability distributions from samples in high dimension. Unlike plug-in rules that simply replace the true distributions by their empirical counterparts, our method promotes couplings with low transport rank, a new structural assumption that is similar to the nonnegative rank of a matrix. Regularizing based on this assumption leads to drastic improvements on high-dimensional data for various tasks, including domain adaptation in single-cell RNA sequencing data. These findings are supported by a theoretical analysis that indicates that the transport rank is key in overcoming the curse of dimensionality inherent to data-driven optimal transport.


title: Motion Planning for a Humanoid Mobile Manipulator System

abstract: A high redundant non-holonomic humanoid mobile dual-arm manipulator system is presented in this paper where the motion planning to realize "human-like" autonomous navigation and manipulation tasks is studied. Firstly, an improved MaxiMin NSGA-II algorithm, which optimizes five objective functions to solve the problems of singularity, redundancy, and coupling between mobile base and manipulator simultaneously, is proposed to design the optimal pose to manipulate the target object. Then, in order to link the initial pose and that optimal pose, an off-line motion planning algorithm is designed. In detail, an efficient direct-connect bidirectional RRT and gradient descent algorithm is proposed to reduce the sampled nodes largely, and a geometric optimization method is proposed for path pruning. Besides, head forward behaviors are realized by calculating the reasonable orientations and assigning them to the mobile base to improve the quality of human-robot interaction. Thirdly, the extension to on-line planning is done by introducing real-time sensing, collision-test and control cycles to update robotic motion in dynamic environments. Fourthly, an EEs' via-point-based multi-objective genetic algorithm is proposed to design the "human-like" via-poses by optimizing four objective functions. Finally, numerous simulations are presented to validate the effectiveness of proposed algorithms.


title: Dark-exciton based strain sensing in transition metal dichalcogenides

abstract: The trend towards ever smaller high-performance devices in modern technology requires novel materials with new functionalities. The recent emergence of atomically thin two-dimensional (2D) materials has opened up possibilities for the design of ultra-thin and flexible nanoelectronic devices. As truly 2D materials, they exhibit an optimal surface-to-volume ratio, which results in an extremely high sensitivity to external changes. This makes these materials optimal candidates for sensing applications. Here, we exploit the remarkably diverse exciton landscape in monolayer transition metal dichalcogenides to propose a novel dark-exciton-based concept for ultra sensitive strain sensors. We demonstrate that the dark-bright-exciton separation can be controlled by strain, which has a crucial impact on the activation of dark excitonic states. This results in a pronounced intensity change of dark excitons in photoluminescence spectra, when only 0.05 $\%$ strain is applied. The predicted extremely high optical gauge factors of up to 8000 are promising for the design of optical strain sensors.


title: Opportunistic Scheduling in Underlay Cognitive Radio based Systems: User Selection Probability Analysis

abstract: In this paper, an underlay cognitive radio (CR) system is considered with multiple cognitive or secondary users contending to transmit their information to the cognitive destination (e.g., eNodeB) using the spectral resource of a primary user. The novel closed-form expressions are derived for the selection probabilities of cognitive users with opportunistic scheduling wherein an optimal metric is employed for opportunistic transmission. The analytical results corroborated by the Monte Carlo simulations, can be used to demonstrate the fairness achieved in opportunistic scheduling. It is shown that the fairness in terms of equal chance for transmission amongst all cognitive users can only be seen for the scenarios when the fraction of distances between the cognitive transmitter and cognitive receiver, and cognitive transmitter and primary receiver is identical for each of the cognitive transmitters.


title: Sub-Weyl bounds for $GL(2)$ $L$-functions

abstract: In this paper we obtain a sub-Weyl bound for $L(1/2+it,f)$ for $f$ a Hecke modular form.


title: Faster SGD training by minibatch persistency

abstract: It is well known that, for most datasets, the use of large-size minibatches for Stochastic Gradient Descent (SGD) typically leads to slow convergence and poor generalization. On the other hand, large minibatches are of great practical interest as they allow for a better exploitation of modern GPUs. Previous literature on the subject concentrated on how to adjust the main SGD parameters (in particular, the learning rate) when using large minibatches. In this work we introduce an additional feature, that we call minibatch persistency, that consists in reusing the same minibatch for K consecutive SGD iterations. The computational conjecture here is that a large minibatch contains a significant sample of the training set, so one can afford to slightly overfitting it without worsening generalization too much. The approach is intended to speedup SGD convergence, and also has the advantage of reducing the overhead related to data loading on the internal GPU memory. We present computational results on CIFAR-10 with an AlexNet architecture, showing that even small persistency values (K=2 or 5) already lead to a significantly faster convergence and to a comparable (or even better) generalization than the standard "disposable minibatch" approach (K=1), in particular when large minibatches are used. The lesson learned is that minibatch persistency can be a simple yet effective way to deal with large minibatches.


title: From Planetesimal to Planet in Turbulent Disks. II. Formation of Gas Giant Planets

abstract: In the core accretion scenario, gas giant planets are formed form solid cores with several Earth masses via gas accretion. We investigate the formation of such cores via collisional growth from kilometer-sized planetesimals in turbulent disks. The stirring by forming cores induces collisional fragmentation and surrounding planetesimals are ground down until radial drift. The core growth is therefore stalled by the depletion of surrounding planetesiamls due to collisional fragmentation and radial drift. The collisional strength of planetesimals determines the planetesimal-depletion timescale, which is prolonged for large planetesiamls. The size of planetesiamls around growing cores is determined by the planetesimal size distribution at the onset of runaway growth. Strong turbulence delays the onset of runaway growth, resulting in large planetesimals. Therefore, the core mass evolution depends on turbulent parameter $\alpha$; the formation of cores massive enough without significant depletion of surrounding planetesimals needs strong turbulence of $\alpha \gtrsim 10^{-3}$. However, the strong turbulence with $\alpha \gtrsim 10^{-3}$ leads to a significant delay of the onset of runaway growth and prevents the formation of massive cores within the disk lifetime. The formation of cores massive enough within several millions years therefore requires the several times enhancement of the solid surface densities, which is achieved in the inner disk $\lesssim 10$AU due to pile-up of drifting dust aggregates. In addition, the collisional strength $Q_{\rm D}^*$ even for kilometer-sized or smaller bodies affects the growth of cores; $Q_{\rm D}^* \gtrsim 10^7 \,{\rm erg/g}$ for bodies $\lesssim 1$ km is likely for this gas giant formation.


title: Search for the Higgs boson produced in association with a vector boson and decaying into two spin-zero particles in the $H \rightarrow aa \rightarrow 4b$ channel in $pp$ collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

abstract: A search for exotic decays of the Higgs boson into a pair of spin-zero particles, $H \rightarrow aa$, where the $a$-boson decays into $b$-quarks promptly or with a mean proper lifetime $c\tau_a$ up to 6 mm and has a mass in the range of $20-60$ GeV, is presented. The search is performed in events where the Higgs boson is produced in association with a $W$ or $Z$ boson, giving rise to a signature of one or two charged leptons (electrons or muons) and multiple jets from $b$-quark decays. The analysis is based on the dataset of proton-proton collisions at $\sqrt{s} = 13$ TeV recorded in 2015 and 2016 by the ATLAS detector at the CERN Large Hadron Collider, corresponding to an integrated luminosity of 36.1 fb$^{-1}$. No significant excess of events above the Standard Model background prediction is observed, and 95% confidence-level upper limits are derived for the production cross-sections for $pp \rightarrow WH$, $ZH$ and their combination, times the branching ratio of the decay chain $H \rightarrow aa \rightarrow 4b$. For $a$-bosons which decay promptly, the upper limit on the combination of cross-sections for $WH$ and $ZH$ times the branching ratio of $H \rightarrow aa \rightarrow 4b$ ranges from 3.0 pb for $m_a= 20$ GeV to 1.3 pb for $m_a = 60$ GeV, assuming that the ratio of $WH$ to $ZH$ cross-sections follows the Standard Model prediction. For $a$-bosons with longer proper lifetimes, the most stringent limits are 1.8 pb and 0.6 pb, respectively, at $c\tau_a \sim 0.4$ mm.


title: Deterministic $O(1)$-Approximation Algorithms to 1-Center Clustering with Outliers

abstract: The 1-center clustering with outliers problem asks about identifying a prototypical robust statistic that approximates the location of a cluster of points. Given some constant $0 < \alpha < 1$ and $n$ points such that $\alpha n$ of them are in some (unknown) ball of radius $r,$ the goal is to compute a ball of radius $O(r)$ that also contains $\alpha n$ points. This problem can be formulated with the points in a normed vector space such as $\mathbb{R}^d$ or in a general metric space. The problem has a simple randomized solution: a randomly selected point is a correct solution with constant probability, and its correctness can be verified in linear time. However, the deterministic complexity of this problem was not known. In this paper, for any $\ell_p$ vector space, we show an $O(nd)$-time solution with a ball of radius $O(r)$ for a fixed $\alpha > \frac{1}{2},$ and for any normed vector space, we show an $O(nd)$-time solution with a ball of radius $O(r)$ when $\alpha > \frac{1}{2}$ as well as an $O (nd \log^{(k)}(n))$-time solution with a ball of radius $O(r)$ for all $\alpha > 0, k \in \mathbb{N},$ where $\log^{(k)}(n)$ represents the $k$th iterated logarithm, assuming distance computation and vector space operations take $O(d)$ time. For an arbitrary metric space, we show for any $C \in \mathbb{N}$ an $O(n^{1+1/C})$-time solution that finds a ball of radius $2Cr,$ assuming distance computation between any pair of points takes $O(1)$-time. Moreover, this algorithm is optimal for general metric spaces, as we show that for any fixed $\alpha, C,$ there is no $o(n^{1+1/C})$-query and thus no $o(n^{1+1/C})$-time solution that deterministically finds a ball of radius $2Cr$.


title: Records from partial comparisons and discrete approximations

abstract: In this paper we study records obtained from partial comparisons within a sequence of independent and identically distributed (i.i.d.) random variables, indexed by positive integers, with a common density~\(f.\) Our main result is that if the comparison sets along a subsequence of the indices satisfy a certain compatibility property, then the corresponding record events are independent. Moreover, the record event probabilities do not depend on the density~\(f\) and we obtain closed form expressions for the distribution of~\(r^{th}\) record value for any integer~\(r \geq 1.\) Our proof techniques extend to the discrete case as well and we estimate the difference in record event probabilities associated with a continuous random variable~\(X\) and its discrete approximations.


title: No Threshold graphs are cospectral

abstract: A threshold graph G on n vertices is defined by binary sequence of length n. In this paper we present an explicit formula for computing the characteristic polynomial of a threshold graph from its binary sequence. Applications include obtaining a formula for the determinant of adjacency matrix of a threshold graph and showing that no two nonisomorphic threshold graphs are cospectral.


title: Numerical Simulation of magnetized jet creation using a hollow ring of laser beams

abstract: Three dimensional FLASH magneto-hydrodynamics(MHD) modeling is carried out to interpret the OMEGA laser experiments of strongly magnetized, highly collimated jets driven by a ring of 20 OMEGA beams. The predicted optical Thomson scattering spectra and proton images are in good agreement with a subset of the experimental data. Magnetic fields generated via the Biermann battery term are amplified at the boundary between the core and the surrounding of the jet. The simulation predicts multiple axially aligned magnetic flux ropes with alternating poloidal component. Future applications of the hollow ring configuration in laboratory astrophysics are discussed.


title: A simple test for ideal memristors

abstract: An ideal memristor is defined as a resistor with memory that, when subject to a time-dependent current, $I(t)$, its resistance $R_M(q)$ depends {\it only} on the charge $q$ that has flowed through it, so that its voltage response is $V(t)=R_M(q)I(t)$. It has been argued that a clear fingerprint of these ideal memristors is a pinched hysteresis loop in their I-V curves. However, a pinched I-V hysteresis loop is not a definitive test of whether a resistor with memory is truly an ideal memristor because such a property is shared also by other resistors whose memory depends on additional internal state variables, other than the charge. Here, we introduce a very simple and {\it unambiguous} test that can be utilized to check experimentally if a resistor with memory is indeed an ideal memristor. Our test is based on the duality property of a capacitor-memristor circuit whereby, for any initial resistance states of the memristor and any form of the applied voltage, the final state of an ideal memristor must be identical to its initial state, if the capacitor charge finally returns to its initial value. In actual experiments, a sufficiently wide range of voltage amplitudes and initial states are enough to perform the test. The proposed test can help resolve some long-standing controversies still existing in the literature about whether an ideal memristor does actually exist or it is a purely mathematical concept.


title: Comparing the emission spectra of U and Th hollow cathode lamps and a new U line-list

abstract: Thorium hollow cathode lamps (HCLs) are used as frequency calibrators for many high resolution astronomical spectrographs, some of which aim for Doppler precision at the 1 m/s level. We aim to determine the most suitable combination of elements (Th or U, Ar or Ne) for wavelength calibration of astronomical spectrographs, to characterize differences between similar HCLs, and to provide a new U line-list. We record high resolution spectra of different HCLs using a Fourier transform spectrograph: (i) U-Ne, U-Ar, Th-Ne, and Th-Ar lamps in the spectral range from 500 to 1000 nm and U-Ne and U-Ar from 1000 to 1700 nm; (ii) we systematically compare the number of emission lines and the line intensity ratio for a set of 12 U-Ne HCLs; and (iii) we record a master spectrum of U-Ne to create a new U line-list. Uranium lamps show more lines suitable for calibration than Th lamps from 500 to 1000 nm. The filling gas of the lamps significantly affects their performance because Ar and Ne lines contaminate different spectral regions. We find differences (up to 88 %) in the line intensity of U lines in different lamps from the same batch. We find 8239 isolated lines between 500 and 1700 nm that we attribute to U, 3379 of which were not contained in earlier line-lists. The U line-list is available at the http://www.astro.physik.uni-goettingen.de/research/U_atlas . We suggest using a combination of U-Ne and U-Ar lamps to wavelength-calibrate astronomical spectrographs up to 1000 nm. From 1000 to 1700 nm, U-Ne shows better properties. The differences in line strength between different HCLs underline the importance of characterizing HCLs in the laboratory. The new 3379 U lines can significantly improve the radial velocity precision of astronomical spectrographs.


title: Quadratic Approximation of Generalized Tribonacci Sequences

abstract: In this paper, we give quadratic approximation of generalized Tribonacci sequence $\{V_{n}\}_{n\geq0}$ defined by Eq. (\ref{eq:7}) and use this result to give the matrix form of the $n$-th power of a companion matrix of $\{V_{n}\}_{n\geq0}$. Then we re-prove the cubic identity or Cassini-type formula for $\{V_{n}\}_{n\geq0}$ and the Binet's formula of the generalized Tribonacci quaternions.


title: GOE Statistics for Levy Matrices

abstract: In this paper we establish eigenvector delocalization and bulk universality for L\'{e}vy matrices, which are real, symmetric, $N \times N$ random matrices $\textbf{H}$ whose upper triangular entries are independent, identically distributed $\alpha$-stable laws. First, if $\alpha \in (1, 2)$ and $E \in \mathbb{R}$ is any energy bounded away from $0$, we show that every eigenvector of $\textbf{H}$ corresponding to an eigenvalue near $E$ is completely delocalized and that the local spectral statistics of $\textbf{H}$ around $E$ converge to those of the Gaussian Orthogonal Ensemble (GOE) as $N$ tends to $\infty$. Second, we show for almost all $\alpha \in (0, 2)$, there exists a constant $c(\alpha) > 0$ such that the same statements hold if $|E| < c (\alpha)$.


title: Majority vote model with ancillary noise in complex networks

abstract: We analyze the properties of the majority-vote (MV) model with an additional noise in which a local spin can be changed independently of its neighborhood. In the standard MV, one of the simplest nonequilibrium systems exhibiting an order-disorder phase transition, spins are aligned with their local majority with probability $1-f$, and with complementary probability $f$, the majority rule is not followed. In the noisy MV (NMV), a random spin flip is succeeded with probability $p$ (with complementary $1-p$ the usual MV rule is accomplished). Such extra ingredient was considered by Vieira and Crokidakis [Physica A {\bf 450}, 30 (2016)] for the square lattice. Here, we generalize the NMV for arbitrary networks, including homogeneous [random regular (RR) and Erd\"os Renyi (ER)] and heterogeneous [Barabasi-Albert (BA)] structures, through mean-field calculations and numerical simulations. Results coming from both approaches are in excellent agreement with each other, revealing that the presence of additional noise does not affect the classification of phase transition, which remains continuous irrespective of the network degree and its distribution. The critical point and the threshold probability $p_t$ marking the disappearance of the ordered phase depend on the node distribution and increase with the connectivity $k$. The critical behavior, investigated numerically, exhibits a common set of critical exponents for RR and ER topologies, but different from BA and regular lattices. Finally, our results indicate that (in contrary to a previous proposition) there is no first-order transition in the NMV for large $k$.


title: State equation from the spectral structure of human brain activity

abstract: Neural electromagnetic (EM) signals recorded non-invasively from individual human subjects vary in complexity and magnitude. Nonetheless, variation in neural activity has been difficult to quantify and interpret, due to complex, broad-band features in the frequency domain. Studying signals recorded with magnetoencephalography (MEG) from healthy young adult subjects while in resting and active states, a systematic framework inspired by thermodynamics is applied to neural EM signals. Despite considerable inter-subject variation in terms of spectral entropy and energy across time epochs, data support the existence of a robust and linear relationship defining an effective state equation, with higher energy and lower entropy in the resting state compared to active, consistently across subjects. Mechanisms underlying the emergence of relationships between empirically measured effective state functions are further investigated using a model network of coupled oscillators, suggesting an interplay between noise and coupling strength can account for coherent variation of empirically observed quantities. Taken together, the results show macroscopic neural observables follow a robust, non-trivial conservation rule for energy modulation and information generation.


title: Neural Ordinary Differential Equations

abstract: We introduce a new family of deep neural network models. Instead of specifying a discrete sequence of hidden layers, we parameterize the derivative of the hidden state using a neural network. The output of the network is computed using a black-box differential equation solver. These continuous-depth models have constant memory cost, adapt their evaluation strategy to each input, and can explicitly trade numerical precision for speed. We demonstrate these properties in continuous-depth residual networks and continuous-time latent variable models. We also construct continuous normalizing flows, a generative model that can train by maximum likelihood, without partitioning or ordering the data dimensions. For training, we show how to scalably backpropagate through any ODE solver, without access to its internal operations. This allows end-to-end training of ODEs within larger models.


title: Observation of Quadrupole Transitions and Edge Mode Topology in an LC network

abstract: High-order topological insulators are a recent development extending the topological theory of charge polarization to higher multipole moments. Since their theoretical proposal, several experimental realizations of high-order topological insulators have been reported. However, high order topological transitions have not been observed. In this letter, we report on the observation of a high-order topological transition in a quadrupole topological insulator implemented in an LC circuit with nonlinear couplings. This system presents the ability to confine electromagnetic energy in its corner states, with a localization length that can be tuned over a broad range through the use of an external bias voltage. Additionally, we provide an experimental characterization and an effective theory for the boundary states, further corroborating their topological nature by direct measurement of the winding number.


title: Relating the cut distance and the weak* topology for graphons

abstract: The theory of graphons is ultimately connected with the so-called cut norm. In this paper, we approach the cut norm topology via the weak* topology (when considering a predual of $L^{1}$-functions). We prove that a sequence $W_1,W_2,W_3,\ldots$ of graphons converges in the cut distance if and only if we have equality of the sets of weak* accumulation points and of weak* limit points of all sequences of graphons $W_1',W_2',W_3',\ldots$ that are weakly isomorphic to $W_1,W_2,W_3,\ldots$. We further give a short descriptive set theoretic argument that each sequence of graphons contains a subsequence with the property above. This in particular provides an alternative proof of the theorem of Lov\'asz and Szegedy about compactness of the space of graphons. We connect these results to "multiway cut" characterization of cut distance convergence from [Ann. of Math. (2) 176 (2012), no. 1, 151-219]. These results are more naturally phrased in the Vietoris hyperspace $K$ over graphons with the weak* topology. We show that graphons with the cut distance topology are homeomorphic to a closed subset of $K$, and deduce several consequences of this fact. From these concepts a new order on the space of graphons emerges. This order allows to compare how structured two graphons are. We establish basic properties of this "structurdness order".


title: A Graph-Theoretic Analysis of Distributed Replicator Dynamic

abstract: This paper attempts to develop a graph-theoretic multi-agent perspective of population games to study the \quotes{truncation} behavior. The proposed method considers fitness of the population as a dynamical system to address the issue of restrictive description of this behavior which pertains to the underlying population dynamic. The fitness dynamic resembles an agreement protocol that enables comments on the steady-state characteristics of the graph that represents the population structure. The structural attributes of the underlying graph and the truncation behavior are emphasized by exploiting the spectral properties of the associated Laplacian matrix. The asymptotic stability of the fitness agreement protocol has been shown to be sufficient in concluding the stability of population dynamics. Simulation results validating the proposed hypothesis have been discussed.


title: Constructing Fast Network through Deconstruction of Convolution

abstract: Convolutional neural networks have achieved great success in various vision tasks; however, they incur heavy resource costs. By using deeper and wider networks, network accuracy can be improved rapidly. However, in an environment with limited resources (e.g., mobile applications), heavy networks may not be usable. This study shows that naive convolution can be deconstructed into a shift operation and pointwise convolution. To cope with various convolutions, we propose a new shift operation called active shift layer (ASL) that formulates the amount of shift as a learnable function with shift parameters. This new layer can be optimized end-to-end through backpropagation and it can provide optimal shift values. Finally, we apply this layer to a light and fast network that surpasses existing state-of-the-art networks.


title: Object-Oriented Dynamics Predictor

abstract: Generalization has been one of the major challenges for learning dynamics models in model-based reinforcement learning. However, previous work on action-conditioned dynamics prediction focuses on learning the pixel-level motion and thus does not generalize well to novel environments with different object layouts. In this paper, we present a novel object-oriented framework, called object-oriented dynamics predictor (OODP), which decomposes the environment into objects and predicts the dynamics of objects conditioned on both actions and object-to-object relations. It is an end-to-end neural network and can be trained in an unsupervised manner. To enable the generalization ability of dynamics learning, we design a novel CNN-based relation mechanism that is class-specific (rather than object-specific) and exploits the locality principle. Empirical results show that OODP significantly outperforms previous methods in terms of generalization over novel environments with various object layouts. OODP is able to learn from very few environments and accurately predict dynamics in a large number of unseen environments. In addition, OODP learns semantically and visually interpretable dynamics models.


title: Learning Unit State Recognition Based on Multi-channel Data Fusion

abstract: Despite recent advances in MOOC, the current e-learning systems have advantages of alleviating barriers by time differences, and geographically spatial separation between teachers and students. However, there has been a 'lack of supervision' problem that e-learner's learning unit state(LUS) can't be supervised automatically. In this paper, we present a fusion framework considering three channel data sources: 1) videos/images from a camera, 2) eye movement information tracked by a low solution eye tracker and 3) mouse movement. Based on these data modalities, we propose a novel approach of multi-channel data fusion to explore the learning unit state recognition. We also propose a method to build a learning state recognition model to avoid manually labeling image data. The experiments were carried on our designed online learning prototype system, and we choose CART, Random Forest and GBDT regression model to predict e-learner's learning state. The results show that multi-channel data fusion model have a better recognition performance in comparison with single channel model. In addition, a best recognition performance can be reached when image, eye movement and mouse movement features are fused.


title: Few-Shot Segmentation Propagation with Guided Networks

abstract: Learning-based methods for visual segmentation have made progress on particular types of segmentation tasks, but are limited by the necessary supervision, the narrow definitions of fixed tasks, and the lack of control during inference for correcting errors. To remedy the rigidity and annotation burden of standard approaches, we address the problem of few-shot segmentation: given few image and few pixel supervision, segment any images accordingly. We propose guided networks, which extract a latent task representation from any amount of supervision, and optimize our architecture end-to-end for fast, accurate few-shot segmentation. Our method can switch tasks without further optimization and quickly update when given more guidance. We report the first results for segmentation from one pixel per concept and show real-time interactive video segmentation. Our unified approach propagates pixel annotations across space for interactive segmentation, across time for video segmentation, and across scenes for semantic segmentation. Our guided segmentor is state-of-the-art in accuracy for the amount of annotation and time. See http://github.com/shelhamer/revolver for code, models, and more details.


title: Multilingual Scene Character Recognition System using Sparse Auto-Encoder for Efficient Local Features Representation in Bag of Features

abstract: The recognition of texts existing in camera-captured images has become an important issue for a great deal of research during the past few decades. This give birth to Scene Character Recognition (SCR) which is an important step in scene text recognition pipeline. In this paper, we extended the Bag of Features (BoF)-based model using deep learning for representing features for accurate SCR of different languages. In the features coding step, a deep Sparse Auto-encoder (SAE)-based strategy was applied to enhance the representative and discriminative abilities of image features. This deep learning architecture provides more efficient features representation and therefore a better recognition accuracy. Our system was evaluated extensively on all the scene character datasets of five different languages. The experimental results proved the efficiency of our system for a multilingual SCR.


title: Distinguishing Refracted Features using Light Field Cameras with Application to Structure from Motion

abstract: Robots must reliably interact with refractive objects in many applications; however, refractive objects can cause many robotic vision algorithms to become unreliable or even fail, particularly feature-based matching applications, such as structure-from-motion. We propose a method to distinguish between refracted and Lambertian image features using a light field camera. Specifically, we propose to use textural cross-correlation to characterise apparent feature motion in a single light field, and compare this motion to its Lambertian equivalent based on 4D light field geometry. Our refracted feature distinguisher has a 34.3% higher rate of detection compared to state-of-the-art for light fields captured with large baselines relative to the refractive object. Our method also applies to light field cameras with much smaller baselines than previously considered, yielding up to 2 times better detection for 2D-refractive objects, such as a sphere, and up to 8 times better for 1D-refractive objects, such as a cylinder. For structure from motion, we demonstrate that rejecting refracted features using our distinguisher yields up to 42.4% lower reprojection error, and lower failure rate when the robot is approaching refractive objects. Our method lead to more robust robot vision in the presence of refractive objects.


title: Semantic Analysis of (Reflectional) Visual Symmetry: A Human-Centred Computational Model for Declarative Explainability

abstract: We present a computational model for the semantic interpretation of symmetry in naturalistic scenes. Key features include a human-centred representation, and a declarative, explainable interpretation model supporting deep semantic question-answering founded on an integration of methods in knowledge representation and deep learning based computer vision. In the backdrop of the visual arts, we showcase the framework's capability to generate human-centred, queryable, relational structures, also evaluating the framework with an empirical study on the human perception of visual symmetry. Our framework represents and is driven by the application of foundational, integrated Vision and Knowledge Representation and Reasoning methods for applications in the arts, and the psychological and social sciences.


title: Transfer Learning for Related Reinforcement Learning Tasks via Image-to-Image Translation

abstract: Despite the remarkable success of Deep RL in learning control policies from raw pixels, the resulting models do not generalize. We demonstrate that a trained agent fails completely when facing small visual changes, and that fine-tuning---the common transfer learning paradigm---fails to adapt to these changes, to the extent that it is faster to re-train the model from scratch. We show that by separating the visual transfer task from the control policy we achieve substantially better sample efficiency and transfer behavior, allowing an agent trained on the source task to transfer well to the target tasks. The visual mapping from the target to the source domain is performed using unaligned GANs, resulting in a control policy that can be further improved using imitation learning from imperfect demonstrations. We demonstrate the approach on synthetic visual variants of the Breakout game, as well as on transfer between subsequent levels of Road Fighter, a Nintendo car-driving game. A visualization of our approach can be seen in https://youtu.be/4mnkzYyXMn4 and https://youtu.be/KCGTrQi6Ogo .


title: Localizing and Quantifying Damage in Social Media Images

abstract: Traditional post-disaster assessment of damage heavily relies on expensive GIS data, especially remote sensing image data. In recent years, social media has become a rich source of disaster information that may be useful in assessing damage at a lower cost. Such information includes text (e.g., tweets) or images posted by eyewitnesses of a disaster. Most of the existing research explores the use of text in identifying situational awareness information useful for disaster response teams. The use of social media images to assess disaster damage is limited. In this paper, we propose a novel approach, based on convolutional neural networks and class activation maps, to locate damage in a disaster image and to quantify the degree of the damage. Our proposed approach enables the use of social network images for post-disaster damage assessment and provides an inexpensive and feasible alternative to the more expensive GIS approach.


title: DeepTerramechanics: Terrain Classification and Slip Estimation for Ground Robots via Deep Learning

abstract: Terramechanics plays a critical role in the areas of ground vehicles and ground mobile robots since understanding and estimating the variables influencing the vehicle-terrain interaction may mean the success or the failure of an entire mission. This research applies state-of-the-art algorithms in deep learning to two key problems: estimating wheel slip and classifying the terrain being traversed by a ground robot. Three data sets collected by ground robotic platforms (MIT single-wheel testbed, MSL Curiosity rover, and tracked robot Fitorobot) are employed in order to compare the performance of traditional machine learning methods (i.e. Support Vector Machine (SVM) and Multi-layer Perceptron (MLP)) against Deep Neural Networks (DNNs) and Convolutional Neural Networks (CNNs). This work also shows the impact that certain tuning parameters and the network architecture (MLP, DNN and CNN) play on the performance of those methods. This paper also contributes a deep discussion with the lessons learned in the implementation of DNNs and CNNs and how these methods can be extended to solve other problems.


title: Deep Sequence Learning with Auxiliary Information for Traffic Prediction

abstract: Predicting traffic conditions from online route queries is a challenging task as there are many complicated interactions over the roads and crowds involved. In this paper, we intend to improve traffic prediction by appropriate integration of three kinds of implicit but essential factors encoded in auxiliary information. We do this within an encoder-decoder sequence learning framework that integrates the following data: 1) offline geographical and social attributes. For example, the geographical structure of roads or public social events such as national celebrations; 2) road intersection information. In general, traffic congestion occurs at major junctions; 3) online crowd queries. For example, when many online queries issued for the same destination due to a public performance, the traffic around the destination will potentially become heavier at this location after a while. Qualitative and quantitative experiments on a real-world dataset from Baidu have demonstrated the effectiveness of our framework.


title: G2D: from GTA to Data

abstract: This document describes G2D, a software that enables capturing videos from Grand Theft Auto V (GTA V), a popular role playing game set in an expansive virtual city. The target users of our software are computer vision researchers who wish to collect hyper-realistic computer-generated imagery of a city from the street level, under controlled 6DOF camera poses and varying environmental conditions (weather, season, time of day, traffic density, etc.). G2D accesses/calls the native functions of the game; hence users can directly interact with G2D while playing the game. Specifically, G2D enables users to manipulate conditions of the virtual environment on the fly, while the gameplay camera is set to automatically retrace a predetermined 6DOF camera pose trajectory within the game coordinate system. Concurrently, automatic screen capture is executed while the virtual environment is being explored. G2D and its source code are publicly available at https://goo.gl/SS7fS6 In addition, we demonstrate an application of G2D to generate a large-scale dataset with groundtruth camera poses for testing structure-from-motion (SfM) algorithms. The dataset and generated 3D point clouds are also made available at https://goo.gl/DNzxHx


title: In situ TensorView: In situ Visualization of Convolutional Neural Networks

abstract: Convolutional Neural Networks(CNNs) are complex systems. They are trained so they can adapt their internal connections to recognize images, texts and more. It is both interesting and helpful to visualize the dynamics within such deep artificial neural networks so that people can understand how these artificial networks are learning and making predictions. In the field of scientific simulations, visualization tools like Paraview have long been utilized to provide insights and understandings. We present in situ TensorView to visualize the training and functioning of CNNs as if they are systems of scientific simulations. In situ TensorView is a loosely coupled in situ visualization open framework that provides multiple viewers to help users to visualize and understand their networks. It leverages the capability of co-processing from Paraview to provide real-time visualization during training and predicting phases. This avoid heavy I/O overhead for visualizing large dynamic systems. Only a small number of lines of codes are injected in TensorFlow framework. The visualization can provide guidance to adjust the architecture of networks, or compress the pre-trained networks. We showcase visualizing the training of LeNet-5 and VGG16 using in situ TensorView.


title: Two Stream Self-Supervised Learning for Action Recognition

abstract: We present a self-supervised approach using spatio-temporal signals between video frames for action recognition. A two-stream architecture is leveraged to tangle spatial and temporal representation learning. Our task is formulated as both a sequence verification and spatio-temporal alignment tasks. The former task requires motion temporal structure understanding while the latter couples the learned motion with the spatial representation. The self-supervised pre-trained weights effectiveness is validated on the action recognition task. Quantitative evaluation shows the self-supervised approach competence on three datasets: HMDB51, UCF101, and Honda driving dataset (HDD). Further investigations to boost performance and generalize validity are still required.


title: The Tully-Fisher relation for flat galaxies

abstract: We construct a multiparametric Tully-Fisher (TF) relation for a large sample of edge-on galaxies from the Revised Flat Galaxy Catalog using HI data from the EDD database and parameters from the EGIS catalog. We incorporate a variety of additional parameters including structural parameters of edge-on galaxies in different bandpasses. Besides the rotation curve maximum, only the HI-to-optical luminosity ratio and optical colours play a statistically significant role in the multiparametric TF relation. We are able to decrease the standard deviation of the multiparametric TF relation down to 0.32 mag, which is at the level of best modern samples of galaxies used for studies of the matter motion in the Universe via the TF-relation.


title: Detecting and interpreting myocardial infarction using fully convolutional neural networks

abstract: Objective: We aim to provide an algorithm for the detection of myocardial infarction that operates directly on ECG data without any preprocessing and to investigate its decision criteria. Approach: We train an ensemble of fully convolutional neural networks on the PTB ECG dataset and apply state-of-the-art attribution methods. Main results: Our classifier reaches 93.3% sensitivity and 89.7% specificity evaluated using 10-fold cross-validation with sampling based on patients. The presented method outperforms state-of-the-art approaches and reaches the performance level of human cardiologists for detection of myocardial infarction. We are able to discriminate channel-specific regions that contribute most significantly to the neural network's decision. Interestingly, the network's decision is influenced by signs also recognized by human cardiologists as indicative of myocardial infarction. Significance: Our results demonstrate the high prospects of algorithmic ECG analysis for future clinical applications considering both its quantitative performance as well as the possibility of assessing decision criteria on a per-example basis, which enhances the comprehensibility of the approach.


title: Optimizing optical tweezing with directional scattering in composite microspheres

abstract: Trapping of microspheres with a single focused laser beam is usually limited to materials with relative refractive indexes slightly larger than one. We show that directional light scattering can be employed to optically trap high-index materials. For this purpose, we propose a material platform to achieve zero backward scattering (ZBS), also known as the first Kerker condition, in a composite media containing spherical inclusions of silica embedded in a SiC microsphere. By tuning the volume filling fraction of inclusions and the microsphere radius, stable trapping can be achieved, provided that ZBS is combined with the condition for destructive interference between the fields reflected at the external and internal interfaces of the microsphere when located at the focal point. We show that our proposal also holds even in the presence of a significant amount of spherical aberration, which is a common condition in most optical tweezers setups. In this case, achieving ZBS is essential for trapping high-index materials.


title: On parametric Gevrey asymptotics for some initial value problems in two asymmetric complex time variables

abstract: We study a family of nonlinear initial value partial differential equations in the complex domain under the action of two asymmetric time variables. Different Gevrey bounds and multisummability results are obtain depending on each element of the family, providing a more complete picture on the asymptotic behavior of the solutions of PDEs in the complex domain in several complex variables. The main results lean on a fixed point argument in certain Banach space in the Borel plane, together with a Borel summability procedure and the action of different Ramis-Sibuya type theorems.


title: Gravitational Radiation from Color-Kinematics Duality

abstract: We perturbatively calculate classical radiation in Yang-Mills theory and dilaton gravity, to next-to-leading order in couplings. The radiation is sourced by the scattering of two relativistic massive scalar sources with the dynamical effect taken into account, corresponding to the post-Minkowskian regime in gravity. We show how to arrange the Yang-Mills radiation such that the duality between colors and kinematics is manifest, including the three-term Jacobi identity. The search for duality-satisfying expressions exploits an auxiliary bi-adjoint scalar theory as a guide for locality. The double copy is obtained by replacing the color factors in Yang-Mills with kinematic counterparts, following Bern-Carrasco-Johansson construction in S- matrix. On the gravity side, the radiation is directly computed at the third post-Minkowskian order with massive sources. We find perfect agreement between observables in dilaton gravity and the Yang-Mills double copy. This non-trivially generalizes the leading-order rules by Goldberger and Ridgway. For the first time, the kinematic Jacobi identity appears beyond field-theory S-matrix, suggesting that the color-kinematics duality holds more generally. Our results offer a path for simplifying analytical calculations in post-Minkowskian regime.


title: Density profiles of ultracompact minihalos: Implications for constraining the primordial power spectrum

abstract: Enhanced density fluctuations on small scales would lead to the formation of numerous dark matter minihalos, so limits on the minihalo abundance can place upper bounds on the small-scale primordial power spectrum. In particular, the ultracompact minihalo (UCMH), a dark matter structure hypothesized to possess a $\rho\propto r^{-9/4}$ density profile due to its formation at $z\geq 1000$, has been used to establish an upper bound on the primordial power spectrum at scales smaller than 2 Mpc. The extreme slope of this density profile amplifies the observational signals of UCMHs. However, we recently showed via N-body simulations that the $\rho\propto r^{-9/4}$ density profile does not develop in realistic formation scenarios, throwing UCMH-derived power spectrum constraints into question. Instead, minihalos develop shallower inner profiles with power-law indices between $-3/2$ and $-1$. In this paper, we expand on that result and discuss its implications. Using a model that is calibrated to simulation results and predicts halo structures in spiked power spectra based on their formation times, we calculate new upper bounds on the primordial power spectrum based on limits on the dark matter annihilation rate within the Galaxy. We find that despite assuming shallower profiles, this minihalo model actually yields stronger constraints than the previous UCMH picture owing to its inclusion of all minihalos instead of only the earliest-forming ones.


title: Fornax: a Flexible Code for Multiphysics Astrophysical Simulations

abstract: This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes. Instead of focusing exclusively on tests of neutrino radiation hydrodynamics relevant to the core-collapse supernova problem for which Fornax is primarily intended, we present here classical and rigorous demonstrations of code performance relevant to a broad range of multi-dimensional hydrodynamic and multi-group radiation hydrodynamic problems. Our code solves the comoving-frame radiation moment equations using the M1 closure, utilizes conservative high-order reconstruction, employs semi-explicit matter and radiation transport via a high-order time stepping scheme, and is suitable for application to a wide range of astrophysical problems. To this end, we first describe the philosophy, algorithms, and methodologies of Fornax and then perform numerous stringent code tests, that collectively and vigorously exercise the code, demonstrate the excellent numerical fidelity with which it captures the many physical effects of radiation hydrodynamics, and show excellent strong scaling well above 100k MPI tasks.

categories: astro-ph.IM astro-ph.HE astro-ph.SR physics.comp-ph


title: The Measurement of Polarization in Radio Astronomy

abstract: Modern dual-polarization receivers allow a radio telescope to characterize the full polarization state of incoming insterstellar radio waves. Many astronomers incorrectly consider a polarimeter to be the "backend" of the telescope. We go to lengths to dissuade the reader of this concept: the backend is the least complicated component of the radio telescope when it comes to measuring polarization. The feed, telescope structure, dish surface, coaxial cables, optical fibers, and electronics can each alter the polarization state of the received astronomical signal. We begin with an overview of polarized radiation, introducing Jones and Stokes vectors, and then discuss construction of digitized pseudo-Stokes vectors from the outputs of modern correlators. We describe the measurement and calibration process for polarization observations and illustrate how instrumental polarization can affect a measurement. Finally, we draw attention to the confusion generated by various polarization conventions and highlight the need for observers to state all adopted conventions when reporting polarization results.


title: Chasing Lyman alpha-emitting galaxies at z = 8.8

abstract: With a total integration time of 168 hours and a narrowband (NB) filter tuned to Lyman alpha at z = 8.8, the UltraVISTA survey has set out to find some of the most distant galaxies, on the verge of the Epoch of Reionization. Previous calculations of the expected number of detected Lya-emitting galaxies (LAEs) at this redshift did not explicitly take into account the radiative transfer (RT) of Lya. In this work we combine a theoretical model for the halo mass function with numerical results from high-res cosmological hydrosimulations with LyC+Lya RT, assessing the visibility of LAEs residing in these halos. Uncertainties such as cosmic variance and the anisotropic escape of Lya are taken into account, and it is predicted that once the survey has finished, the probabilities of detecting none, one, or more than one are ~90%, ~10%, and ~1%; a significantly smaller success rate compared to earlier predictions, due to the combined effect of a highly neutral IGM scattering Lya to such large distances from the galaxy that they fall outside the observational aperture, and to the actual depth of the survey being less than predicted. Because the IGM affects NB and broadband (BB) magnitudes differently, we argue for a relaxed color selection criterion of NB - BB ~ +0.85. But since the flux is continuum-dominated, even if a galaxy is detectable in the NB its probability of being selected as a NB excess object is <~35%. Various properties of galaxies at this redshift are predicted, e.g. UV and Lya LFs, M*-Mh relation, spectral shape, optimal aperture, and the anisotropic escape of Lya through both a dusty ISM and a partly neutral IGM. Finally, we describe and publish a fast numerical code for adding numbers with asymmetric uncertainties ("x_{-sigma_1}^{+sigma_2}") proving to be significantly better than the standard, but wrong, way of adding upper and lower uncertainties in quadrature separately.


title: X-shooter and ALMA spectroscopy of GRB 161023A - A study of metals and molecules in the line of sight towards a luminous GRB

abstract: Long gamma-ray bursts are produced during the dramatic deaths of massive stars with very short lifetimes, meaning that they explode close to the birth place of their progenitors. During a short period they become the most luminous objects observable in the Universe, being perfect beacons to study high-redshift star-forming regions. To use the afterglow of GRB 161023A at a redshift $z=2.710$ as a background source to study the environment of the explosion and the intervening systems along its line-of-sight. r the first time, we complement UV/Optical/NIR spectroscopy with millimetre spectroscopy using ALMA, which allows us to probe the molecular content of the host galaxy. The X-shooter spectrum shows a plethora of absorption features including fine-structure and metastable transitions of Fe, Ni, Si, C and O. We present photometry ranging from 43 s to over 500 days after the burst. We infer a host-galaxy metallicity of [Zn/H] $=-1.11\pm0.07$, which corrected for dust depletion results in [X/H] $=-0.94\pm0.08$. We do not detect molecular features in the ALMA data, but we derive limits on the molecular content of $log(N_{CO}/cm^{-2})<15.7$ and $log(N_{HCO+}/cm^{-2})<13.2$, which are consistent with those that we obtain from the optical spectra, $log(N_{H_2}/cm^{-2})<15.2$ and $log(N_{CO}/cm^{-2})<14.5$. Within the host galaxy we detect three velocity systems through UV/Optical/NIR absorption spectroscopy, all with levels that were excited by the GRB afterglow. We determine the distance from these systems to the GRB to be in the range between 0.7 and 1.0 kpc. The sight-line to GRB 161023A shows 9 independent intervening systems, most of them with multiple components. (Abridged)


title: The formation of high-mass binary star systems

abstract: We develop a semi-analytic model to investigate how accretion onto wide low-mass binary stars can result in a close high-mass binary system. The key ingredient is to allow mass accretion while limiting the gain in angular momentum. We envision this process as being regulated by an external magnetic field during infall. Molecular clouds are made to collapse spherically with material either accreting onto the stars or settling in a disk. Our aim is to determine what initial conditions are needed for the resulting binary to be both massive and close. Whether material accretes, and what happens to the binary separation as a result, depends on the relative size of its specific angular momentum, compared to the specific angular momentum of the binary. When we add a magnetic field we are introducing a torque to the system which is capable of stripping the molecular cloud of some of its angular momentum, and consequently easing the formation of high-mass binaries. Our results suggest that clouds in excess of 1000 M$_\odot$ and radii of 0.5 pc or larger, can easily form binary systems with masses in excess of 25 M$_\odot$ and separations of order 10 R$_\odot$ with magnetic fields of order 100 {\mu}G (mass-to-flux ratios of order 5).


title: A Detection of Sgr A* in the far infrared

abstract: We report the first detection of the Galactic Centre massive black hole, Sgr~A*, in the far infrared. Our measurements were obtained with PACS on board the \emph{Herschel} satellite at $100~\mathrm{\mu m}$ and $160~\mathrm{\mu m}$. While the warm dust in the Galactic Centre is too bright to allow for a direct detection of Sgr~A*, we measure a significant and simultaneous variation of its flux of $\Delta F_{\nu\widehat{=}160 ~\mathrm{\mu m}} = (0.27\pm0.06)~\mathrm{Jy}$ and $\Delta F_{\nu\widehat{=}100 ~\mathrm{\mu m}}= (0.16\pm0.10)~\mathrm{Jy}$ during one observation. The significance level of the $160 ~\mathrm{\mu m}$ band variability is $4.5\sigma$ and the corresponding $100 ~\mathrm{\mu m}$ band variability is significant at $1.6\sigma$. We find no example of an equally significant false positive detection. Conservatively assuming a variability of $25\%$ in the FIR, we can provide upper limits to the flux. Comparing the latter with theoretical models we find that 1D RIAF models have difficulties explaining the observed faintness. However, the upper limits are consistent with modern ALMA and VLA observations. Our upper limits provide further evidence for a spectral peak at $\sim 10^{12} ~ \mathrm{Hz}$ and constrain the number density of $\gamma \sim 100$ electrons in the accretion disk and or outflow.


title: Long-Lived Particles at the Energy Frontier: The MATHUSLA Physics Case

abstract: We examine the theoretical motivations for long-lived particle (LLP) signals at the LHC in a comprehensive survey of Standard Model (SM) extensions. LLPs are a common prediction of a wide range of theories that address unsolved fundamental mysteries such as naturalness, dark matter, baryogenesis and neutrino masses, and represent a natural and generic possibility for physics beyond the SM (BSM). In most cases the LLP lifetime can be treated as a free parameter from the $\mu$m scale up to the Big Bang Nucleosynthesis limit of $\sim 10^7$m. Neutral LLPs with lifetimes above $\sim$ 100m are particularly difficult to probe, as the sensitivity of the LHC main detectors is limited by challenging backgrounds, triggers, and small acceptances. MATHUSLA is a proposal for a minimally instrumented, large-volume surface detector near ATLAS or CMS. It would search for neutral LLPs produced in HL-LHC collisions by reconstructing displaced vertices (DVs) in a low-background environment, extending the sensitivity of the main detectors by orders of magnitude in the long-lifetime regime. In this white paper we study the LLP physics opportunities afforded by a MATHUSLA-like detector at the HL-LHC. We develop a model-independent approach to describe the sensitivity of MATHUSLA to BSM LLP signals, and compare it to DV and missing energy searches at ATLAS or CMS. We then explore the BSM motivations for LLPs in considerable detail, presenting a large number of new sensitivity studies. While our discussion is especially oriented towards the long-lifetime regime at MATHUSLA, this survey underlines the importance of a varied LLP search program at the LHC in general. By synthesizing these results into a general discussion of the top-down and bottom-up motivations for LLP searches, it is our aim to demonstrate the exceptional strength and breadth of the physics case for the construction of the MATHUSLA detector.


title: Mass Assembly of Stellar Systems and their Evolution with the SMA - 1.3 mm Subcompact Data Release

abstract: We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .


title: Manifestly Soft Gauge Invariant Formulation of vNRQCD

abstract: Homogeneous power counting in Non-Relativistic QCD (NRQCD) implies the simultaneous existence of both soft and ultrasoft gluons. In the velocity renormalization group (vNRQCD) formalism we show that operators involving soft fields interacting with heavy potential quarks can be put in a manifestly gauge invariant form by utilizing gluon and quark building blocks which contain Wilson lines, and are analogous to those used in the soft collinear effective theory. This leads to several simplifications, in particular significantly reducing the size of the operator basis, which stream-lines matching and anomalous dimension calculations at subleading order in the velocity expansion. Also, soft ghosts no longer couple via potential like interactions to the heavy quark fields, and hence do not appear until two loops. Furthermore the the color structures that appear at each order in $\alpha_s$ in the static potential are only those which should arise according to non-Abelian exponentiation. We also discuss how zero-bin subtractions clarify the role of $i\varepsilon$ terms in the heavy quark propagator poles carrying soft momenta. Even though the choice of the direction of the soft Wilson lines in our building blocks does not matter, there is still a limited set of consistent choices, and the predicted form of the eikonal poles does influence results at two-loops and beyond.


title: A critical re-evaluation of the Thorne-\.Zytkow object candidate HV 2112

abstract: It has been argued in the literature that the star HV~2112 in the Small Magellanic Cloud is the first known example of a T\.ZO, a Red Supergiant with a degenerate neutron core. This claim is based on the star having a high luminosity ($\log (L/L_\odot)$~> 5), an extremely cool effective temperature, and a surface enriched in in lithium, calcium and various $irp$-process elements. In this paper we re-examine this evidence, and present new measurements of the stellar properties. By compiling archival photometry from blue to mid-IR for HV~2112 and integrating under it's spectral energy distribution we find a bolometric luminosity in the range of $\log (L/L_\odot)$=4.70-4.91, lower than that found in previous work and comparable to bright asymptotic giant branch (AGB) stars. We compare a VLT+XSHOOTER spectrum of HV~2112 to other late type, luminous SMC stars, finding no evidence for enhancements in Rb, Ca or K, though there does seem to be an enrichment in Li. We therefore conclude that a much more likely explanation for HV~2112 is that it is an intermediate mass($\sim$5M$_\odot$) AGB star. However, from our sample of comparison stars we identify a new T\.ZO candidate, HV~11417, which seems to be enriched in Rb but for which we cannot determine a Li abundance.


title: Joint halo mass function for modified gravity and massive neutrinos I: simulations and cosmological forecasts

abstract: We present a halo mass function accurate over the full relevant Hu-Sawicki $f(R)$ parameter space based on spherical collapse calculations and calibrated to a suite of modified gravity $N$-body simulations that include massive neutrinos. We investigate the ability of current and forthcoming galaxy cluster observations to detect deviations from general relativity while constraining the total neutrino mass and including systematic uncertainties. Our results indicate that the degeneracy between massive neutrino and modify gravity effects is a limiting factor for the current searches for new gravitational physics with clusters of galaxies, but future surveys will be able to break the degeneracy.


title: Entangling Bosonic Modes via an Engineered Exchange Interaction

abstract: The realization of robust universal quantum computation with any platform ultimately requires both the coherent storage of quantum information and (at least) one entangling operation between individual elements. The use of continuous-variable bosonic modes as the quantum element is a promising route to preserve the coherence of quantum information against naturally-occurring errors. However, operations between bosonic modes can be challenging. In analogy to the exchange interaction between discrete-variable spin systems, the exponential-SWAP unitary [$\mathbf{U}_{\mathrm{E}}\left(\theta_c\right)$] can coherently transfer the states between two bosonic modes, regardless of the chosen encoding, realizing a deterministic entangling operation for certain $\theta_c$. Here, we develop an efficient circuit to implement $\mathbf{U}_{\mathrm{E}}\left(\theta_c\right)$ and realize the operation in a three-dimensional circuit QED architecture. We demonstrate high-quality deterministic entanglement between two cavity modes with several different encodings. Our results provide a crucial primitive necessary for universal quantum computation using bosonic modes.