title: Spiral instabilities: Mechanism for recurrence

abstract: We argue that self-excited instabilities are the cause of spiral patterns in simulations of unperturbed stellar discs. In previous papers, we have found that spiral patterns were caused by a few concurrent waves, which we claimed were modes. The superposition of a few steadily rotating waves inevitably causes the appearance of the disc to change continuously, and creates the kind of shearing spiral patterns that have been widely reported. Although we have found that individual modes last for relatively few rotations, spiral activity persists because fresh instabilities appear, which we suspected were excited by the changes to the disc caused by previous disturbances. Here we confirm our suspicion by demonstrating that scattering at either of the Lindblad resonances seeds a new groove-type instability. With this logical gap closed, our understanding of the behaviour in the simulations is almost complete. We believe that our robust mechanism is a major cause of spiral patterns in the old stellar discs of galaxies, including the Milky Way where we have previously reported evidence for resonance scattering in the recently released Gaia data.


title: Cation Disorder and Lithium Insertion Mechanism of Wadsley--Roth Crystallographic Shear Phases

abstract: Wadsley--Roth crystallographic shear phases form a family of compounds that have attracted attention due to their excellent performance as lithium-ion battery electrodes. The complex crystallographic structure of these materials poses a challenge for first-principles computational modelling and hinders the understanding of their structural, electronic and dynamic properties. In this article, we study three different niobium-tungsten oxide crystallographic shear phases (Nb$_{12}$WO$_{33}$, Nb$_{14}$W$_{3}$O$_{44}$, Nb$_{16}$W$_5$O$_{55}$) using an enumeration-based approach and first-principles density-functional theory calculations. We report common principles governing the cation disorder, lithium insertion mechanism, and electronic structure of these materials. Tungsten preferentially occupies tetrahedral and block-central sites within the block-type crystal structures. The lithium insertion proceeds via a three-step mechanism, associated with an anisotropic evolution of the host lattice. Our calculations reveal an important connection between long-range and local structural changes: in the second step of the mechanism, the removal of local structural distortions leads to the contraction of the lattice along specific crystallographic directions, buffering the volume expansion of the material. Niobium-tungsten oxide shear structures host small amounts of localised electrons during initial lithium insertion due to the confining effect of the blocks, but quickly become metallic upon further lithiation. We argue that the combination of local, long-range, and electronic structural evolution over the course of lithiation is beneficial to the performance of these materials as battery electrodes. The mechanistic principles we establish arise from the compound-independent crystallographic shear structure, and are therefore likely to apply to Ti/Nb oxide or pure Nb oxide shear phases.


title: Adaptive optics benefit for quantum key distribution uplink from ground to a satellite

abstract: For quantum communications, the use of Earth-orbiting satellites to extend distances has gained significant attention in recent years, exemplified in particular by the launch of the Micius satellite in 2016. The performance of applied protocols such as quantum key distribution (QKD) depends significantly upon the transmission efficiency through the turbulent atmosphere, which is especially challenging for ground-to-satellite uplink scenarios. Adaptive optics (AO) techniques have been used in astronomical, communication, and other applications to reduce the detrimental effects of turbulence for many years, but their applicability to quantum protocols, and their requirements specifically in the uplink scenario, are not well established. Here, we model the effect of the atmosphere on link efficiency between an Earth station and a satellite using an optical uplink, and how AO can help recover from loss due to turbulence. Examining both low-Earth-orbit and geostationary uplink scenarios, we find that a modest link transmissivity improvement of about 3dB can be obtained in the case of a co-aligned downward beacon, while the link can be dramatically improved, up to 7dB, using an offset beacon, such as a laser guide star. AO coupled with a laser guide star would thus deliver a significant increase in the secret key generation rate of the QKD ground-to-space uplink system, especially as reductions of channel loss have favourably nonlinear key-rate response within this high-loss regime.


title: N\'{e}el and stripe ordering from spin-orbital entanglement in $\alpha$-Sr$_2$CrO$_4$

abstract: The rich phenomenology engendered by the coupling between the spin and orbital degrees of freedom has become appreciated as a key feature of many strongly-correlated electron systems. The resulting emergent physics is particularly prominent in a number of materials, from Fe-based unconventional superconductors to transition metal oxides, including manganites and vanadates. Here, we investigate the electronic ground states of $\alpha$-Sr$_2$CrO$_4$, a compound that is a rare embodiment of the spin-1 Kugel-Khomskii model on the square lattice -- a paradigmatic platform to capture the physics of coupled magnetic and orbital electronic orders. We have used resonant X-ray diffraction at the Cr-$K$ edge to reveal N\'{e}el magnetic order at the in-plane wavevector $\mathbf{Q}_N = (1/2, 1/2)$ below $T_N = 112$ K, as well as an additional electronic order at the 'stripe' wavevector $\mathbf{Q}_s = (1/2, 0)$ below T$_s$ $ \sim 50$ K. These findings are examined within the framework of the Kugel-Khomskii model by a combination of mean-field and Monte-Carlo approaches, which supports the stability of the spin N\'{e}el phase with subsequent lower-temperature stripe orbital ordering, revealing a candidate mechanism for the experimentally observed peak at $\mathbf{Q}_s$. On the basis of these findings, we propose that $\alpha$-Sr$_2$CrO$_4$ serves as a new platform in which to investigate multi-orbital physics and its role in the low-temperature phases of Mott insulators.


title: Revisiting the 16 Cygni planet host at unprecedented precision and exploring automated tools for precise abundances

abstract: The binary system 16 Cygni is key in studies of the planet-star chemical composition connection, as only one of the stars is known to host a planet. This allows us to better assess the possible influence of planet interactions on the chemical composition of stars that are born from the same cloud and thus, should have a similar abundance pattern. In our previous work, we found clear abundance differences for elements with Z$\leq30$ between both components of this system, and a trend of these abundances as a function of the condensation temperature (T$_{c}$), which suggests a spectral chemical signature related to planet formation. In this work we show that our previous findings are still consistent even if we include more species, like the volatile N and neutron capture elements (Z $>$ 30). We report a slope with T$_{c}$ of $1.56 \pm 0.24 \times 10^{-5}$ dex K$^{-1}$, that is good agreement with both our previous work and recent results by Nissen and collaborators. We also performed some tests using ARES and iSpec to automatic measure the equivalent width and found T$_c$ slopes in reasonable agreement with our results as well. In addition, we determine abundances for Li and Be by spectral synthesis, finding that 16 Cyg A is richer not only in Li but also in Be, when compared to its companion. This may be evidence of planet engulfment, indicating that the T$_{c}$ trend found in this binary system may be a chemical signature of planet accretion in the A component, rather than a imprint of the giant planet rocky core formation on 16 Cyg B.


title: Evolution of primordial black hole spin due to Hawking radiation

abstract: Near extremal Kerr black holes are subject to the Thorne limit $a<a^*_{\rm lim}=0.998$ in the case of thin disc accretion, or some generalized version of this in other disc geometries. However any limit that differs from the thermodynamics limit $a^* < 1$ can in principle be evaded in other astrophysical configurations, and in particular if the near extremal black holes are primordial and subject to evaporation by Hawking radiation only. We derive the lower mass limit above which Hawking radiation is slow enough so that a primordial black hole with a spin initially above some generalized Thorne limit can still be above this limit today. Thus, we point out that the observation of Kerr black holes with extremely high spin should be a hint of either exotic astrophysical mechanisms or primordial origin.


title: Gravitational-wave detection rates for compact binaries formed in isolation: LIGO/Virgo O3 and beyond

abstract: Using simulations performed with the population synthesis code MOBSE, we compute the merger rate densities and detection rates of compact binary mergers formed in isolation for second- and third-generation gravitational-wave detectors. We estimate how rates are affected by uncertainties on key stellar-physics parameters, namely common envelope evolution and natal kicks. We estimate how future upgrades will increase the size of the available catalog of merger events, and we discuss features of the merger rate density that will become accessible with third-generation detectors.


title: Classifying galaxies according to their HI content

abstract: We use machine learning to classify galaxies according to their HI content, based on both their optical photometry and environmental properties. The data used for our analyses are the outputs in the range $z = 0-1$ from MUFASA cosmological hydrodynamic simulation. In our previous paper, where we predicted the galaxy HI content using the same input features, HI rich galaxies were only selected for the training. In order for the predictions on real observation data to be more accurate, the classifiers built in this study will first establish if a galaxy is HI rich ($\rm{log(M_{HI}/M_{*})} > -2 $) before estimating its neutral hydrogen content using the regressors developed in the first paper. We resort to various machine learning algorithms and assess their performance with various metrics such as accuracy for instance. The performance of the classifiers gets better with increasing redshift and reaches their peak performance around $z = 1$. Random Forest method, the most robust among the classifiers when considering only the mock data for both training and test in this study, reaches an accuracy above $98.6 \%$ at $z = 0$ and above $99.0 \%$ at $z = 1$. We test our algorithms, trained with simulation data, on classification of the galaxies in RESOLVE, ALFALFA and GASS surveys. Interestingly, SVM algorithm, the best classifier for the tests, achieves a precision, the relevant metric for the tests, above $87.60\%$ and a specificity above $71.4\%$ with all the tests, indicating that the classifier is capable of learning from the simulated data to classify HI rich/HI poor galaxies from the real observation data. With the advent of large HI 21 cm surveys such as the SKA, this set of classifiers, together with the regressors developed in the first paper, will be part of a pipeline, a very useful tool, which is aimed at predicting HI content of galaxies.


title: Synthesis of Computable Regular Functions of Infinite Words

abstract: Regular functions from infinite words to infinite words can be equivalently specified by MSO-transducers, streaming $\omega$-string transducers as well as deterministic two-way transducers with look-ahead. In their one-way restriction, the latter transducers define the class of rational functions. Even though regular functions are robustly characterised by several finite-state devices, even the subclass of rational functions may contain functions which are not computable (by a Turing machine with infinite input). This paper proposes a decision procedure for the following synthesis problem: given a regular function $f$ (equivalently specified by one of the aforementioned transducer model), is $f$ computable and if it is, synthesize a Turing machine computing it. For regular functions, we show that computability is equivalent to continuity, and therefore the problem boils down to deciding continuity. We establish a generic characterisation of continuity for functions preserving regular languages under inverse image (such as regular functions). We exploit this characterisation to show the decidability of continuity (and hence computability) of rational and regular functions. For rational functions, we show that this can be done in $\mathsf{NLogSpace}$ (it was already known to be in $\mathsf{PTime}$ by Prieur). In a similar fashion, we also effectively characterise uniform continuity of regular functions, and relate it to the notion of uniform computability, which offers stronger efficiency guarantees.


title: D1-D5-P superstrata in 5 and 6 dimensions: separable wave equations and prepotentials

abstract: We construct the most general single-mode superstrata in 5 dimensions with ambipolar, two centered Gibbons Hawking bases, via dimensional reduction of superstrata in 6 dimensions. Previously, asymptotically $\text{AdS}_{3}\times \mathbb{S}^{2}$ 5-dimensional superstrata have been produced, giving microstate geometries of black strings in 5 dimensions. Our construction produces asymptotically $\text{AdS}_{2}\times \mathbb{S}^{3}$ geometries as well, the first instances of superstrata describing the microstate geometries of black holes in 5 dimensions. New examples of superstrata with separable massless wave equations in both 5 and 6 dimensions are uncovered. A $\mathbb{Z}_{2}$ symmetry which identifies distinct 6-dimensional superstrata when reduced to 5 dimensions is found. Finally we use the mathematical structure of the underlying hyper-K\"{a}hler bases to produce prepotentials for the superstrata fluxes in 5 dimensions and uplift them to apply in 6 dimensions as well.


title: End-to-End CAD Model Retrieval and 9DoF Alignment in 3D Scans

abstract: We present a novel, end-to-end approach to align CAD models to an 3D scan of a scene, enabling transformation of a noisy, incomplete 3D scan to a compact, CAD reconstruction with clean, complete object geometry. Our main contribution lies in formulating a differentiable Procrustes alignment that is paired with a symmetry-aware dense object correspondence prediction. To simultaneously align CAD models to all the objects of a scanned scene, our approach detects object locations, then predicts symmetry-aware dense object correspondences between scan and CAD geometry in a unified object space, as well as a nearest neighbor CAD model, both of which are then used to inform a differentiable Procrustes alignment. Our approach operates in a fully-convolutional fashion, enabling alignment of CAD models to the objects of a scan in a single forward pass. This enables our method to outperform state-of-the-art approaches by $19.04\%$ for CAD model alignment to scans, with $\approx 250\times$ faster runtime than previous data-driven approaches.


title: Altruism in populations at the extinction transition

abstract: We study the evolution of cooperation as a birth-death process in spatially extended populations. The benefit from the altruistic behavior of a cooperator is implemented by decreasing the death rate of its direct neighbors. The cost of cooperation is the increase of a cooperator's death rate proportional to the number of its neighbors. When cooperation has higher cost than benefit, cooperators disappear. Then the dynamics reduces to the contact process with defectors as the single particle type. Increasing the benefit-cost ratio above 1, the extinction transition of the contact process splits into a set of nonequilibrium transitions between four regimes when increasing the baseline death rate $p$ as a control parameter: (i) defection only, (ii) coexistence, (iii) cooperation only, (iv) extinction. We investigate the transitions between these regimes. As the main result, we find that full cooperation is established at the extinction transition as long as benefit is strictly larger than cost. Qualitatively identical phase diagrams are obtained for populations on square lattices and in pair approximation. Spatial correlations with nearest neighbors only are thus sufficient for sustained cooperation.


title: Antlia2's role in driving the ripples in the outer gas disk of the Galaxy

abstract: We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer HI disk of the Milky Way, using both test particle and Smoothed Particle Hydrodynamics simulations. We find that orbits with low pericenters, $\sim$ 10 kpc, produce disturbances that match the observed outer HI disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of $(\mu_{\alpha}cos\delta, \mu_{\delta}) = (-0.068,0.032) \pm (0.023,-0.031)~\rm mas/yr$, which agrees with results from Torrealba et al. (2019) within the errors, but gives lower mean pericenters, e.g., $\sim$ 15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis (Chakrabarti \& Blitz 2009) of the dwarf that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.


title: Dark matter capture in celestial objects: Improved treatment of multiple scattering and updated constraints from white dwarfs

abstract: We revisit dark matter (DM) capture in celestial objects, including the impact of multiple scattering, and obtain updated constraints on the DM-proton cross section using observations of white dwarfs. Considering a general form for the energy loss distribution in each scattering, we derive an exact formula for the capture probability through multiple scatterings. We estimate the maximum number of scatterings that $can$ take place, in contrast to the number $required$ to bring a dark matter particle to rest. We employ these results to compute a "dark" luminosity $L_{\rm DM}$, arising solely from the thermalized annihilation products of the captured dark matter. Demanding that $L_{\rm DM}$ not exceed the luminosity of the white dwarfs in the M4 globular cluster, we set a bound on the DM-proton cross section: $\sigma_{p} \lesssim 10^{-44} {\rm cm}^2$, almost independent of the dark matter mass between 100 GeV and 1 PeV and mildly weakening beyond. This is a stronger constraint than those obtained by direct detection experiments in both large mass $\left(M \gtrsim 5 \,\,\rm TeV\right)$ and small mass $\left(M \lesssim 10\,\, \rm GeV\right)$ regimes. For dark matter lighter than 350 MeV, which is beyond the sensitivity of present direct detection experiments, this is the strongest available constraint.


title: POSSIS: predicting spectra, light curves and polarization for multi-dimensional models of supernovae and kilonovae

abstract: We present POSSIS, a time-dependent three-dimensional Monte Carlo code for modelling radiation transport in supernovae and kilonovae. The code incorporates wavelength- and time-dependent opacities and predicts viewing-angle dependent spectra, light curves and polarization for both idealized and hydrodynamical explosion models. We apply the code to a kilonova model with two distinct ejecta components, one including lanthanide elements with relatively high opacities and the other devoid of lanthanides and characterized by lower opacities. We find that a model with total ejecta mass $M_\mathrm{ej}=0.04\,M_\odot$ and half-opening angle of the lanthanide-rich component $\Phi=30^\circ$ provides a good match to GW 170817 / AT 2017gfo for orientations near the polar axis (i.e. for a system viewed close to face-on). We then show how crucial is the use of self-consistent multi-dimensional models in place of combining one-dimensional models to infer important parameters as the ejecta masses. We finally explore the impact of $M_\mathrm{ej}$ and $\Phi$ on the synthetic observables and highlight how the relatively fast computation times of POSSIS make it well-suited to perform parameter-space studies and extract key properties of supernovae and kilonovae. Spectra calculated with POSSIS in this and future studies will be made publicly available.


title: Qualifying threshold of take off stage for successfully disseminated creative ideas

abstract: The creative process is essentially Darwinian and only a small proportion of creative ideas are selected for further development. However, the threshold that identifies this small fraction of successfully disseminated creative ideas at their early stage has not been thoroughly analyzed through the lens of Rogers innovation diffusion theory. Here, we take highly cited (top 1%) research papers as an example of the most successfully disseminated creative ideas and explore the time it takes and citations it receives at their take off stage, which play a crucial role in the dissemination of creativity. Results show the majority of highly cited papers will reach 10% and 25% of their total citations within two years and four years, respectively. Interestingly, our results also present a minimal number of articles that attract their first citation before publication. As for the discipline, number of references, and Price index, we find a significant difference exists: Clinical, Pre-Clinical & Health and Life Sciences are the first two disciplines to reach the C10% and C25% in a shorter amount of time. Highly cited papers with limited references usually take more time to reach 10% and 25% of their total citations. In addition, highly cited papers will attract citations rapidly when they cite more recent references. These results provide insights into the timespan and citations for a research paper to become highly cited at the take off stage in its diffusion process, as well as the factors that may influence it.


title: Dynamics of singular complex analytic vector fields with essential singularities II

abstract: The singular complex analytic vector fields $X$ on the Riemann sphere $\widehat{\mathbb C}_z$ belonging to the family ${\mathscr E}(r,d)=\left\{ X(z)=\frac{1}{P(z)} e^{E(z)}\frac{\partial }{\partial z}\ \Big\vert \ P, E\in\mathbb{C}[z]\right\}$, where $P$ is monic, $deg(P)=r$, $deg(E)=d$, $r+d\geq 1$, have a finite number of poles on the complex plane and an isolated essential singularity at infinity (for $d\geq 1$). Our aim is to describe geometrically $X$, particularly the singularity at infinity. We use the natural one to one correspondence between $X$, a global singular analytic distinguished parameter $\Psi_X(z)=\int^z P(\zeta) e^{-E(\zeta)}d\zeta$, and the Riemann surface ${\mathcal R}_X$ of this distinguished parameter. We introduce $(r,d)$-configuration trees which are weighted directed rooted trees. An $(r,d)$-configuration tree completely encodes the Riemann surface ${\mathcal R}_X$ and the singular flat metric associated on ${\mathcal R}_X$. The $(r,d)$-configuration trees provide "parameters" for the complex manifold ${\mathscr E}(r,d)$, which give explicit geometrical and dynamical information; a valuable tool for the analytic description of $X\in{\mathscr E}(r,d)$. Furthermore, given $X$, the phase portrait of the associated real vector field $Re(X)$ on the Riemann sphere is decomposed into $Re(X)$-invariant components: half planes and finite height strips. The germ of $X$ at infinity is described as a combinatorial word (consisting of hyperbolic, elliptic, parabolic and entire angular sectors having the point at infinity of $\widehat{\mathbb C}_z$ as center). The structural stability, under perturbation in ${\mathscr E}(r,d)$, of the phase portrait of $Re(X)$ is characterized by using the $(r,d)$-configuration trees. We provide explicit conditions, in terms of $r$ and $d$, as to when the number of topologically equivalent phase portraits of $Re(X)$ is unbounded.


title: Forecasted constraints on modified gravity from Sunyaev Zel'dovich tomography

abstract: Observational cosmology has become an important laboratory for testing General Relativity, with searches for modified gravity forming a significant portion of the science case for existing and future surveys. In this paper, we illustrate how future measurements of the Cosmic Microwave Background (CMB) temperature and polarization anisotropies can be combined with large galaxy surveys to improve constraints on modified gravity using the technique of Sunyaev Zel'dovich (SZ) tomography. SZ tomography uses the correlations between the kinetic/polarized SZ contributions to the small-angular scale CMB and the distribution of structure measured in a galaxy redshift survey to reconstruct the remote dipole and quadrupole fields, e.g. the CMB dipole and quadrupole observed throughout the Universe. We compute the effect of a class of modifications of gravity on the remote dipole and quadrupole fields, illustrating that these observables combine a number of the desirable features of existing probes. We then perform a fisher forecast of constraints on a two-parameter class of modifications of gravity for next-generation CMB experiments and galaxy surveys. By incorporating information from the reconstructed remote dipole and quadrupole fields, we find that it is possible to improve the constraints on this model by a factor of $\sim 2$ beyond what is possible with a galaxy survey alone. We conclude that SZ tomography is a promising method for testing gravity with future cosmological datasets.


title: A fast solver for the narrow capture and narrow escape problems in the sphere

abstract: We present an efficient method to solve the narrow capture and narrow escape problems for the sphere. The narrow capture problem models the equilibrium behavior of a Brownian particle in the exterior of a sphere whose surface is reflective, except for a collection of small absorbing patches. The narrow escape problem is the dual problem: it models the behavior of a Brownian particle confined to the interior of a sphere whose surface is reflective, except for a collection of small patches through which it can escape. Mathematically, these give rise to mixed Dirichlet/Neumann boundary value problems of the Poisson equation. They are numerically challenging for two main reasons: (1) the solutions are non-smooth at Dirichlet-Neumann interfaces, and (2) they involve adaptive mesh refinement and the solution of large, ill-conditioned linear systems when the number of small patches is large. By using the Neumann Green's functions for the sphere, we recast each boundary value problem as a system of first-kind integral equations on the collection of patches. A block-diagonal preconditioner together with a multiple scattering formalism leads to a well-conditioned system of second-kind integral equations and a very efficient approach to discretization. This system is solved iteratively using GMRES. We develop a hierarchical, fast multipole method-like algorithm to accelerate each matrix-vector product. Our method is insensitive to the patch size, and the total cost scales with the number N of patches as O(N log N), after a precomputation whose cost depends only on the patch size and not on the number or arrangement of patches. We demonstrate the method with several numerical examples, and are able to achieve highly accurate solutions with 100,000 patches in one hour on a 60-core workstation.


title: Network-based Fake News Detection: A Pattern-driven Approach

abstract: Fake news gains has gained significant momentum, strongly motivating the need for fake news research. Many fake news detection approaches have thus been proposed, where most of them heavily rely on news content. However, network-based clues revealed when analyzing news propagation on social networks is an information that has hardly been comprehensively explored or used for fake news detection. We bridge this gap by proposing a network-based pattern-driven fake news detection approach. We aim to study the patterns of fake news in social networks, which refer to the news being spread, spreaders of the news and relationships among the spreaders. Empirical evidence and interpretations on the existence of such patterns are provided based on social psychological theories. These patterns are then represented at various network levels (i.e., node-level, ego-level, triad-level, community-level and the overall network) for being further utilized to detect fake news. The proposed approach enhances the explainability in fake news feature engineering. Experiments conducted on real-world data demonstrate that the proposed approach can outperform the state of the arts.


title: Twisted and untwisted negativity spectrum of free fermions

abstract: A basic diagnostic of entanglement in mixed quantum states is known as the positive partial transpose (PT) criterion. Such criterion is based on the observation that the spectrum of the partially transposed density matrix of an entangled state contains negative eigenvalues, in turn, used to define an entanglement measure called the logarithmic negativity. Despite the great success of logarithmic negativity in characterizing bosonic many-body systems, generalizing the operation of PT to fermionic systems remained a technical challenge until recently when a more natural definition of PT for fermions that accounts for the Fermi statistics has been put forward. In this paper, we study the many-body spectrum of the reduced density matrix of two adjacent intervals for one-dimensional free fermions after applying the fermionic PT. We show that in general there is a freedom in the definition of such operation which leads to two different definitions of PT: the resulting density matrix is Hermitian in one case, while it becomes pseudo-Hermitian in the other case. Using the path-integral formalism, we analytically compute the leading order term of the moments in both cases and derive the distribution of the corresponding eigenvalues over the complex plane. We further verify our analytical findings by checking them against numerical lattice calculations.

categories: cond-mat.stat-mech cond-mat.str-el hep-th quant-ph


title: Ultra light Thomas-Fermi Dark Matter

abstract: We investigate the viability of a simple dark matter (DM) model consisting of a single fermion in the context of galactic dynamics. We use a consistent approach that does not presume a particular DM density profile but instead requires that the DM+baryon system is in hydrostatic equilibrium. Using a phenomenological baryon density profile, the model then predicts the DM distribution with a core like behavior close to the galactic center. The presence of supermassive black holes (SMBHs) in the center of large galaxies arises naturally in this framework. Using data from a set of large elliptical and spiral galaxies, and from a small set of dwarf galaxies, we find that the model can explain most of the bulk galactic properties, as well as some of the features observed in the rotation curves, provided the DM mass is in the $\mathcal{O}$(50 eV) range. More precise tests of the model require better modeling of the baryon profile and better control on the uncertainties in the data.


title: The Demand Query Model for Bipartite Matching

abstract: We introduce a `concrete complexity' model for studying algorithms for matching in bipartite graphs. The model is based on the "demand query" model used for combinatorial auctions. Most (but not all) known algorithms for bipartite matching seem to be translatable into this model including exact, approximate, sequential, parallel, and online ones. A perfect matching in a bipartite graph can be found in this model with O(n^{3/2}) demand queries (in a bipartite graph with n vertices on each side) and our main open problem is to either improve the upper bound or prove a lower bound. An improved upper bound could yield "normal" algorithms whose running time is better than the fastest ones known, while a lower bound would rule out a faster algorithm for bipartite matching from within a large class of algorithms. Our main result is a lower bound for finding an approximately maximum size matching in parallel: A deterministic algorithm that runs in n^{o(1)} rounds, where each round can make at most n^{1.99} demand queries cannot find a matching whose size is within n^{o(1)} factor of the maximum. This is in contrast to randomized algorithms that can find a matching whose size is $99\%$ of the maximum in O(\log n) rounds, each making n demand queries.


title: Topology Attack and Defense for Graph Neural Networks: An Optimization Perspective

abstract: Graph neural networks (GNNs) which apply the deep neural networks to graph data have achieved significant performance for the task of semi-supervised node classification. However, only few work has addressed the adversarial robustness of GNNs. In this paper, we first present a novel gradient-based attack method that facilitates the difficulty of tackling discrete graph data. When comparing to current adversarial attacks on GNNs, the results show that by only perturbing a small number of edge perturbations, including addition and deletion, our optimization-based attack can lead to a noticeable decrease in classification performance. Moreover, leveraging our gradient-based attack, we propose the first optimization-based adversarial training for GNNs. Our method yields higher robustness against both different gradient based and greedy attack methods without sacrificing classification accuracy on original graph.


title: Abelian tropical covers

abstract: Let $\mathfrak{A}$ be a finite abelian group. In this article, we classify harmonic $\mathfrak{A}$-covers of a tropical curve $\Gamma$ (which allow dilation along edges and at vertices) in terms of the cohomology group of a suitably defined sheaf on $\Gamma$. We give a realizability criterion for harmonic $\mathfrak{A}$-covers by patching local monodromy data in an extended homology group on $\Gamma$. As an explicit example, we work out the case $\mathfrak{A}=\mathbb{Z}/p\mathbb{Z}$ and explain how realizability for such covers is related to the nowhere-zero flow problem from graph theory.


title: Multiphase modeling of precipitation-induced membrane formation

abstract: We formulate a model for the dynamic growth of a membrane developing in a flow as the result of a precipitation reaction, a situation inspired by recent microfluidic experiments. The precipitating solid introduces additional forces on the fluid and eventually forms a membrane that is fixed in the flow due to adhesion with a substrate. A key challenge is that the location of the immobile membrane is unknown $\textit{a priori}$. To model this situation, we use a multiphase framework with fluid and membrane phases; the aqueous chemicals exist as scalar fields that react within the fluid to induce phase change. To verify that the model exhibits desired fluid-structure behaviors, we make a few simplifying assumptions to obtain a reduced form of the equations that is amenable to exact solution. This analysis demonstrates no-slip behavior on the developing membrane without $\textit{a priori}$ assumptions on its location. The model has applications towards precipitate reactions where the precipitate greatly affects the surrounding flow, a situation appearing in many laboratory and geophysical contexts including the hydrothermal vent theory for the origin of life. More generally, this model can be used to address fluid-structure interaction problems that feature the dynamic generation of structures.


title: Sequential Source Coding for Stochastic Systems Subject to Finite Rate Constraints

abstract: In this paper, we revisit the sequential source coding framework to analyze fundamental performance limitations of discrete-time stochastic control systems subject to feedback data-rate constraints in finite-time horizon. The basis of our results is a new characterization of the lower bound on the minimum total-rate achieved by sequential codes subject to a total (across time) distortion constraint and a computational algorithm that allocates optimally the rate-distortion for any fixed finite-time horizon. This characterization facilitates the derivation of analytical, non-asymptotic, and finite-dimensional lower and upper bounds in two control-related scenarios. (a) A parallel time-varying Gauss-Markov process with identically distributed spatial components that is quantized and transmitted through a noiseless channel to a minimum mean-squared error (MMSE) decoder. (b) A time-varying quantized LQG closed-loop control system, with identically distributed spatial components and with a random data-rate allocation. Our non-asymptotic lower bound on the quantized LQG control problem, reveals the absolute minimum data-rates for (mean square) stability of our time-varying plant for any fixed finite time horizon. We supplement our framework with illustrative simulation experiments.

categories: eess.SY cs.IT cs.SY math.DS math.IT math.OC


title: Magic entanglement renormalization for quantum fields

abstract: Continuous tensor networks are variational wavefunctions proposed in recent years to efficiently simulate quantum field theories (QFTs). Prominent examples include the continuous matrix product state (cMPS) and the continuous multi-scale entanglement renormalization ansatz (cMERA). While the cMPS can approximate ground states of a class of QFT Hamiltonians that are both local and interacting, cMERA is only well-understood for QFTs that are quasi-local and non-interacting. In this paper we propose the magic cMERA, a concrete realization of cMERA for a free boson QFT that simultaneously satisfies four remarkable properties: (i) it is the exact ground state of a strictly local Hamiltonian; (ii) in the massless case, its spectrum of scaling operators is exactly soluble in real space; (iii) it has the short-distance structure of a cMPS; (iv) it is generated by a quasi-local entangler that can be written as a continuous matrix product operator. None of these properties is fulfilled by previous cMERA proposals. Properties (iii)-(iv) establish a firm connection between cMERA and cMPS wavefunctionals, opening the path to applying powerful cMPS numerical techniques, valid for interacting QFTs, also to cMERA calculations.


title: GSTR: Secure Multi-hop Message Dissemination in Connected Vehicles using Social Trust Model

abstract: The emergence of connected vehicles paradigm has made secure communication a key concern amongst the connected vehicles. Communication between the vehicles and Road Side Units (RSUs) is critical to disseminate message among the vehicles. We focus on secure message transmission in connected vehicles using multi_hop social networks environment to deliver the message with varying trustworthiness. We proposed a Geographic Social Trust Routing (GSTR) approach; messages are propagated using multiple hops and by considering the various available users in the vehicular network. GSTR is proposed in an application perspective with an assumption that the users are socially connected. The users are selected based on trustworthiness as defined by social connectivity. The route to send a message is calculated based on the highest trust level of each node by using the nodes social network connections along the path in the network. GSTR determines the shortest route using the trusted nodes along the route for message dissemination. GSTR is made delay tolerant by introducing message storage in the cloud if a trustworthy node is unavailable to deliver the message. We compared the proposed approach with Geographic and Traffic Load based Routing (GTLR), Greedy Perimeter Stateless Routing (GPSR), Trust-based GPSR (T_GPSR). The performance results obtained show that GSTR ensures efficient resource utilization, lower packet losses at high vehicle densities.


title: k-Means Aperture Optimization Applied to Kepler K2 Time Series Photometry of Titan

abstract: Motivated by the Kepler K2 time series of Titan, we present an aperture optimization technique for extracting photometry of saturated moving targets with high temporally- and spatially-varying backgrounds. Our approach uses $k$-means clustering to identify interleaved families of images with similar Point-Spread Function and saturation properties, optimizes apertures for each family independently, then merges the time series through a normalization procedure. By applying $k$-means aperture optimization to the K2 Titan data, we achieve $\leq$0.33% photometric scatter in spite of background levels varying from 15% to 60% of the target's flux. We find no compelling evidence for signals attributable to atmospheric variation on the timescales sampled by these observations. We explore other potential applications of the $k$-means aperture optimization technique, including testing its performance on a saturated K2 eclipsing binary star. We conclude with a discussion of the potential for future continuous high-precision photometry campaigns for revealing the dynamical properties of Titan's atmosphere.


title: Twisted characters and holomorphic symmetries

abstract: We consider holomorphic twists of arbitrary supersymmetric theories in four dimensions. Working in the BV formalism, we rederive classical results characterizing the holomorphic twist of chiral and vector supermultiplets, computing the twist explicitly as a family over the space of nilpotent supercharges in minimal supersymmetry. The BV formalism allows one to work with or without auxiliary fields, according to preference; for chiral superfields, we show that the result of the twist is an identical BV theory, the holomorphic $\beta\gamma$ system with superpotential, independent of whether or not auxiliary fields are included. We compute the character of local operators in this holomorphic theory, demonstrating agreement of the free local operators with the usual index of free fields. The local operators with superpotential are computed via a spectral sequence, and are shown to agree with functions on a formal mapping space into the derived critical locus of the superpotential. We consider the holomorphic theory on various geometries, including Hopf manifolds and products of arbitrary pairs of Riemann surfaces, and offer some general remarks on dimensional reductions of holomorphic theories along the $(n-1)$-sphere to topological quantum mechanics. We also study an infinite-dimensional enhancement of the flavor symmetry in this example, to a recently-studied central extension of the derived holomorphic functions with values in the original Lie algebra that generalizes the familiar Kac--Moody enhancement in two-dimensional chiral theories.


title: Adaptative significance levels in linear regression models with known variance

abstract: The Full Bayesian Significance Test (FBST) for precise hypotheses was presented by Pereira and Stern [Entropy 1(4) (1999) 99-110] as a Bayesian alternative instead of the traditional significance test using p-value. The FBST is based on the evidence in favor of the null hypothesis (H). An important practical issue for the implementation of the FBST is the determination of how large the evidence must be in order to decide for its rejection. In the Classical significance tests, it is known that p-value decreases as sample size increases, so by setting a single significance level, it usually leads H rejection. In the FBST procedure, the evidence in favor of H exhibits the same behavior as the p-value when the sample size increases. This suggests that the cut-off point to define the rejection of H in the FBST should be a sample size function. In this work, the scenario of Linear Regression Models with known variance under the Bayesian approach is considered, and a method to find a cut-off value for the evidence in the FBST is presented by minimizing the linear combination of the averaged type I and type II error probabilities for a given sample size and also for a given dimension of the parametric space.


title: Accurate Low Rank Approximation at a Low Computational Cost

abstract: At sublinear computational cost we greatly decrease the spectral error norm of a crude Low Rank Approximation LRA) of a matrix that has fast decaying spectrum of singular values. Our formal estimates for this decrease is in good accordance with our numerical tests for synthetic and real-world inputs. Moreover, empirically, for a large subclass of inputs with fast decaying spectra, the entire LRA with our escalation was accurate and had sublinear cost.


title: Tuning morphology and thermal transport of asymmetric smart polymer blends by macromolecular engineering

abstract: A grand challenge in designing polymeric materials is to tune their properties by macromolecular engineering. In this context, one of the drawbacks that often limits broader applications under high temperature conditions is their poor thermal conductivity $\kappa$. Using molecular dynamics simulations, we establish a structure-property relationship in hydrogen bonded polymer blends for possible improvement of $\kappa$. For this purpose, we investigate two experimentally relevant hydrogen bonded systems -- one system consists of short poly({N}-acryloyl piperidine) (PAP) blended with longer chains of poly(acrylic acid) (PAA) and the second system is a mixture of PAA and short poly(acrylamide) (PAM) chains. Simulation results show that PAA-PAP blends are at the onset of phase separation over the full range of PAP monomer mole fraction $\phi_{PAP}$, which intensifies even more for $\phi_{PAP} > 0.3$. While PAA and PAP interact with preferential hydrogen bonding, phase separation is triggered by the dominant van der Waals attraction between the hydrophobic side groups of PAP. However, if PAP is replaced with PAM, which has a similar chemical structure as PAP without the hydrophobic side group, PAA-PAM blends show much improved solubility. Better solubility is due to the preferential hydrogen bonding between PAA and PAM. As a result, PAM oligomers act as cross-linking bridges between PAA chains resulting in a three dimensional highly cross-linked network. While $\kappa$ for PAA-PAP blends remain almost invariant with $\phi_{PAP}$, PAA-PAM systems show improved $\kappa$ with increasing PAM concentration and also with respect to PAA-PAP blends. Consistent with the theoretical prediction for the thermal transport of amorphous polymers, we show that $\kappa$ is proportional to the materials stiffness, i.e., the bulk modulus K and sound velocity v of PAA-PAM blends.


title: Label-Agnostic Sequence Labeling by Copying Nearest Neighbors

abstract: Retrieve-and-edit based approaches to structured prediction, where structures associated with retrieved neighbors are edited to form new structures, have recently attracted increased interest. However, much recent work merely conditions on retrieved structures (e.g., in a sequence-to-sequence framework), rather than explicitly manipulating them. We show we can perform accurate sequence labeling by explicitly (and only) copying labels from retrieved neighbors. Moreover, because this copying is label-agnostic, we can achieve impressive performance when transferring to new sequence-labeling tasks without retraining. We additionally consider a dynamic programming approach to sequence labeling in the presence of retrieved neighbors, which allows for controlling the number of distinct (copied) segments used to form a prediction, and leads to both more interpretable and accurate predictions.


title: FASTER Recurrent Networks for Efficient Video Classification

abstract: Typical video classification methods often divide a video into short clips, do inference on each clip independently, then aggregate the clip-level predictions to generate the video-level results. However, processing visually similar clips independently ignores the temporal structure of the video sequence, and increases the computational cost at inference time. In this paper, we propose a novel framework named FASTER, i.e., Feature Aggregation for Spatio-TEmporal Redundancy. FASTER aims to leverage the redundancy between neighboring clips and reduce the computational cost by learning to aggregate the predictions from models of different complexities. The FASTER framework can integrate high quality representations from expensive models to capture subtle motion information and lightweight representations from cheap models to cover scene changes in the video. A new recurrent network (i.e., FAST-GRU) is designed to aggregate the mixture of different representations. Compared with existing approaches, FASTER can reduce the FLOPs by over 10x? while maintaining the state-of-the-art accuracy across popular datasets, such as Kinetics, UCF-101 and HMDB-51.


title: Actions of solvable Baumslag-Solitar groups on hyperbolic metric spaces

abstract: We give a complete list of the cobounded actions of solvable Baumslag-Solitar groups on hyperbolic metric spaces up to a natural equivalence relation. The set of equivalence classes carries a natural partial order first introduced by Abbott-Balasubramanya-Osin, and we describe the resulting poset completely. There are finitely many equivalence classes of actions, and each equivalence class contains the action on a point, a tree, or the hyperbolic plane.


title: Photoproduction Reactions and Non-Strange Baryon Spectroscopy

abstract: We review the last two decades of using photon beams to measure the production of mesons, and in particular the information that can be obtained on the spectrum of light, non-strange baryons. This is a compendium of experimental results, which should be used as a complement to theoretical reviews of the subject. Lists of data sets are given, together with a comprehensive set of references. An indication of the impact of the data is presented with a summary of the results.


title: Psycholinguistics meets Continual Learning: Measuring Catastrophic Forgetting in Visual Question Answering

abstract: We study the issue of catastrophic forgetting in the context of neural multimodal approaches to Visual Question Answering (VQA). Motivated by evidence from psycholinguistics, we devise a set of linguistically-informed VQA tasks, which differ by the types of questions involved (Wh-questions and polar questions). We test what impact task difficulty has on continual learning, and whether the order in which a child acquires question types facilitates computational models. Our results show that dramatic forgetting is at play and that task difficulty and order matter. Two well-known current continual learning methods mitigate the problem only to a limiting degree.


title: Linear and Nonlinear Fractional Elliptic Problems

abstract: This paper surveys recent analytical and numerical research on linear problems for the integral fractional Laplacian, fractional obstacle problems, and fractional minimal graphs. The emphasis is on the interplay between regularity, including boundary behavior, and approximability by piecewise linear finite element methods. We discuss several error estimates on graded meshes, and computational challenges associated to implementing and solving efficiently the ensuing integral equations, along with numerical experiments.


title: Alzheimer's Disease Brain MRI Classification: Challenges and Insights

abstract: In recent years, many papers have reported state-of-the-art performance on Alzheimer's Disease classification with MRI scans from the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset using convolutional neural networks. However, we discover that when we split that data into training and testing sets at the subject level, we are not able to obtain similar performance, bringing the validity of many of the previous studies into question. Furthermore, we point out that previous works use different subsets of the ADNI data, making comparison across similar works tricky. In this study, we present the results of three splitting methods, discuss the motivations behind their validity, and report our results using all of the available subjects.


title: BowNet: Dilated Convolution Neural Network for Ultrasound Tongue Contour Extraction

abstract: Ultrasound imaging is safe, relatively affordable, and capable of real-time performance. One application of this technology is to visualize and to characterize human tongue shape and motion during a real-time speech to study healthy or impaired speech production. Due to the noisy nature of ultrasound images with low-contrast characteristic, it might require expertise for non-expert users to recognize organ shape such as tongue surface (dorsum). To alleviate this difficulty for quantitative analysis of tongue shape and motion, tongue surface can be extracted, tracked, and visualized instead of the whole tongue region. Delineating the tongue surface from each frame is a cumbersome, subjective, and error-prone task. Furthermore, the rapidity and complexity of tongue gestures have made it a challenging task, and manual segmentation is not a feasible solution for real-time applications. Employing the power of state-of-the-art deep neural network models and training techniques, it is feasible to implement new fully-automatic, accurate, and robust segmentation methods with the capability of real-time performance, applicable for tracking of the tongue contours during the speech. This paper presents two novel deep neural network models named BowNet and wBowNet benefits from the ability of global prediction of decoding-encoding models, with integrated multi-scale contextual information, and capability of full-resolution (local) extraction of dilated convolutions. Experimental results using several ultrasound tongue image datasets revealed that the combination of both localization and globalization searching could improve prediction result significantly. Assessment of BowNet models using both qualitatively and quantitatively studies showed them outstanding achievements in terms of accuracy and robustness in comparison with similar techniques.


title: Using generative modelling to produce varied intonation for speech synthesis

abstract: Unlike human speakers, typical text-to-speech (TTS) systems are unable to produce multiple distinct renditions of a given sentence. This has previously been addressed by adding explicit external control. In contrast, generative models are able to capture a distribution over multiple renditions and thus produce varied renditions using sampling. Typical neural TTS models learn the average of the data because they minimise mean squared error. In the context of prosody, taking the average produces flatter, more boring speech: an "average prosody". A generative model that can synthesise multiple prosodies will, by design, not model average prosody. We use variational autoencoders (VAEs) which explicitly place the most "average" data close to the mean of the Gaussian prior. We propose that by moving towards the tails of the prior distribution, the model will transition towards generating more idiosyncratic, varied renditions. Focusing here on intonation, we investigate the trade-off between naturalness and intonation variation and find that typical acoustic models can either be natural, or varied, but not both. However, sampling from the tails of the VAE prior produces much more varied intonation than the traditional approaches, whilst maintaining the same level of naturalness.


title: How much entanglement can be created in a closed system?

abstract: In a closed system, the total number of particles is fixed. We ask how much does this conservation law restrict the amount of entanglement that can be created. We derive a tight upper bound on the bipartite entanglement entropy in closed systems, and find what a maximally entangled state looks like in such a system. Finally, we illustrate numerically on an isolated system of one-dimensional fermionic gas, that the upper bound can be reached during its unitary evolution, when starting in a pure state that emulates a thermal state with high enough temperature. These results are in accordance with current experiments measuring R\'enyi-2 entanglement entropy, all of which employ a particle-conserving Hamiltonian, where our bound acts as a loose bound, and will become especially important for bounding the amount of entanglement that can be spontaneously created, once a direct measurement of entanglement entropy becomes feasible.


title: $f$-wave superfluidity from repulsive interaction in Rydberg-dressed Fermi gas

abstract: Interacting Fermi gas provides an ideal model system to understand unconventional pairing and intertwined orders relevant to a large class of quantum materials. Rydberg-dressed Fermi gas is a recent experimental system where the sign, strength, and range of the interaction can be controlled. The interaction in momentum space has a negative minimum at $q_c$ inversely proportional to the characteristic length-scale in real space, the soft-core radius $r_c$. We show theoretically that single-component (spinless) Rydberg-dressed Fermi gas in two dimensions has a rich phase diagram with novel superfluid and density wave orders due to the interplay of the Fermi momentum $p_F$, interaction range $r_c$, and interaction strength $u_0$. For repulsive bare interactions $u_0>0$, the dominant instability is $f$-wave superfluid for $p_Fr_c\lesssim 2$, and density wave for $p_Fr_c\gtrsim 4$. The $f$-wave pairing in this repulsive Fermi gas is reminiscent of the conventional Kohn-Luttinger mechanism, but has a much higher $T_c$. For attractive bare interactions $u_0<0$, the leading instability is $p$-wave pairing. The phase diagram is obtained from functional renormalization group that treats all competing many-body instabilities in the particle-particle and particle-hole channels on equal footing.


title: Identifying Visible Actions in Lifestyle Vlogs

abstract: We consider the task of identifying human actions visible in online videos. We focus on the widely spread genre of lifestyle vlogs, which consist of videos of people performing actions while verbally describing them. Our goal is to identify if actions mentioned in the speech description of a video are visually present. We construct a dataset with crowdsourced manual annotations of visible actions, and introduce a multimodal algorithm that leverages information derived from visual and linguistic clues to automatically infer which actions are visible in a video. We demonstrate that our multimodal algorithm outperforms algorithms based only on one modality at a time.


title: Multi-Level Taxonomy and Critical Review of Eco-Routing Methods

abstract: Routing decisions are initially based on minimizing travel time. Nevertheless, eco-routing considers the environmental aspect (e.g. emissions, fuel, and exposure) and was introduced to replace the initial routing concept to mitigate the undesirable impact of transportation systems on the environment. This review paper aims to provide a four level taxonomy and map eco-routing studies to the proposed taxonomy. Furthermore, the strengths and weaknesses of the presented models are summarized. In the literature reviewed, special emphasis was given to the role of vehicle connectivity in eco-routing. The main findings include: the microscopic level of aggregation of the flow and emission/fuel models was rarely employed for large case studies due to the associated complexity; one objective was optimized at a time for a majority of the studies; and all of the reviewed studies were applied in a centralized routing system environment.


title: Introducing the Hearthstone-AI Competition

abstract: The Hearthstone AI framework and competition motivates the development of artificial intelligence agents that can play collectible card games. A special feature of those games is the high variety of cards, which can be chosen by the players to create their own decks. In contrast to simpler card games, the value of many cards is determined by their possible synergies. The vast amount of possible decks, the randomness of the game, as well as the restricted information during the player's turn offer quite a hard challenge for the development of game-playing agents. This short paper introduces the competition framework and goes into more detail on the problems and challenges that need to be faced during the development process.


title: Pykg2vec: A Python Library for Knowledge Graph Embedding

abstract: Pykg2vec is an open-source Python library for learning the representations of the entities and relations in knowledge graphs. Pykg2vec's flexible and modular software architecture currently implements 16 state-of-the-art knowledge graph embedding algorithms, and is designed to easily incorporate new algorithms. The goal of pykg2vec is to provide a practical and educational platform to accelerate research in knowledge graph representation learning. Pykg2vec is built on top of TensorFlow and Python's multiprocessing framework and provides modules for batch generation, Bayesian hyperparameter optimization, mean rank evaluation, embedding, and result visualization. Pykg2vec is released under the MIT License and is also available in the Python Package Index (PyPI). The source code of pykg2vec is available at https://github.com/Sujit-O/pykg2vec.


title: Interpreting OWL Complex Classes in AutomationML based on Bidirectional Translation

abstract: The World Wide Web Consortium (W3C) has published several recommendations for building and storing ontologies, including the most recent OWL 2 Web Ontology Language (OWL). These initiatives have been followed by practical implementations that popularize OWL in various domains. For example, OWL has been used for conceptual modeling in industrial engineering, and its reasoning facilities are used to provide a wealth of services, e.g. model diagnosis, automated code generation, and semantic integration. More specifically, recent studies have shown that OWL is well suited for harmonizing information of engineering tools stored as AutomationML (AML) files. However, OWL and its tools can be cumbersome for direct use by engineers such that an ontology expert is often required in practice. Although much attention has been paid in the literature to overcome this issue by transforming OWL ontologies from/to AML models automatically, dealing with OWL complex classes remains an open research question. In this paper, we introduce the AML concept models for representing OWL complex classes in AutomationML, and present algorithms for the bidirectional translation between OWL complex classes and their corresponding AML concept models. We show that this approach provides an efficient and intuitive interface for nonexperts to visualize, modify, and create OWL complex classes.


title: Discovery of the first heavily obscured QSO candidate at $z>6$ in a close galaxy pair

abstract: While theoretical arguments predict that most of the early growth of supermassive black holes (SMBHs) happened during heavily obscured phases of accretion, current methods used for selecting $z>6$ quasars (QSOs) are strongly biased against obscured QSOs, thus considerably limiting our understanding of accreting SMBHs during the first Gyr of the Universe from an observational point of view. We report the $Chandra$ discovery of the first heavily obscured QSO candidate in the early universe, hosted by a close ($\approx5$ kpc) galaxy pair at $z=6.515$. One of the members is an optically classified type 1 QSO, PSO167-13. The companion galaxy was first detected as a [C II] emitter by ALMA. An X-ray source is significantly ($P=0.9996$) detected by $Chandra$ in the 2-5 keV band, with $<1.14$ net counts in the 0.5-2 keV band, although the current positional uncertainty does not allow a conclusive association with either PSO167-13 or its companion galaxy. From X-ray photometry and hardness-ratio arguments, we estimated an obscuring column density of $N_H>2\times10^{24}\,\mathrm{cm^{-2}}$ and $N_H>6\times10^{23}\,\mathrm{cm^{-2}}$ at $68\%$ and $90\%$ confidence levels, respectively. Thus, regardless of which of the two galaxies is associated with the X-ray emission, this source is the first heavily obscured QSO candidate at $z>6$.


title: The Regression Discontinuity Design

abstract: This handbook chapter gives an introduction to the sharp regression discontinuity design, covering identification, estimation, inference, and falsification methods.


title: Applying the tempered Lefschetz thimble method to the Hubbard model away from half-filling

abstract: The tempered Lefschetz thimble method is a parallel-tempering algorithm towards solving the numerical sign problem. It uses the flow time of the gradient flow as a tempering parameter and is expected to tame both the sign and multimodal problems simultaneously. In this paper, we further develop the algorithm so that the expectation values can be estimated precisely with a criterion ensuring global equilibrium and the sufficiency of the sample size. To demonstrate that this algorithm works well, we apply it to the quantum Monte Carlo simulation of the Hubbard model away from half-filling on a two-dimensional lattice of small size, and show that the numerical results agree nicely with exact values.


title: Feasible Surface Plasmon Routing Based on The Self-assembled InGaAs/GaAs Semiconductor Quantum Dot Located between Two Silver Metallic Waveguides

abstract: We proposed an experimentally feasible scheme of nano-plamonic switch and quantum router via the single self-assembled InGaAs/GaAs semiconductor quantum dot (SQD) with a V type three-level energy structure located between two silver metallic waveguides. We studied theoretically transmission and transfer rates of single plasmons in such a multi-ports system via the real-space approach, where our results showed that single plasmons from the input port could be switchable and redirected by controlling parameters, such as the intensity of the classical field, the detunings, and the interaction between the SQD and the waveguides. Our proposed scheme could be used not only in the design of quantum routers and quantum switches for the construction of quantum network, but also in quantum photonic integrated circuits.


title: Decay spectroscopy of $^{50}$Sc and $^{50m}$Sc to $^{50}$Ti

abstract: The $\beta$ decay of the isomeric and ground state of $^{50}$Sc to the semi-magic nucleus $^{50}_{22}$Ti$_{28}$ has been studied using a $^{50}$Ca beam delivered to the GRIFFIN $\gamma$-ray spectrometer at the TRIUMF-ISAC facility. $\beta$-decay branching ratios are reported to 16 excited states with a total of 38 $\gamma$-ray transitions linking them. These new data significantly expands the information available over previous studies. Relative intensities are measured to less than 0.001$\%$ that of the strongest transition with the majority of $\gamma$-ray transitions observed here in $\beta$ decay for the first time. The data are compared to shell-model calculations utilizing both phenomenologically-derived interactions employed in the ${\it pf}$ shell as well as a state-of-the-art, ${\it ab~initio}$ based interaction built in the valence-space in-medium similarity renormalization group framework.


title: Technical Preprint: Rationale and Design of a Planned Observational Study to Evaluate the Impact of Hydrocodone Rescheduling on Opioid Prescribing After Surgery

abstract: In October 2014, the US Drug Enforcement Agency (DEA) reclassified hydrocodone from Schedule III to Schedule II of the Controlled Substances Act, resulting in a prohibition on refills in the initial prescription. While this schedule change was associated with overall decreases in the rate of filled hydrocodone prescriptions and opioid dispensing, available studies conflict regarding its impact on acute opioid prescribing among surgical patients. Here, we present the rationale and design of a planned study to measure the effect of hydrocodone rescheduling using a difference-in-differences design that leverages anticipated variation in the relative impact of this policy on patients treated by surgeons that more or less frequently prescribed hydrocodone products versus other opioids prior to the schedule change. Additionally, we present findings from preliminary study conducted on a subset of our full planned sample to assess for potential differences in outcome trends over the 3 years prior to rescheduling among patients treated by surgeons who commonly prescribed hydrocodone versus those treated by surgeons who rarely prescribed hydrocodone.


title: Central charge and topological invariant of Calabi-Yau manifolds

abstract: F-theory, as a 12-dimensional theory that is a contender of the Theory of Everything, should be compactified into elliptically fibered threefolds or fourfolds of Calabi-Yau. Such manifolds have an elliptic curve as a fiber, and their bases may have singularities. We considered orbifold as simplest non-flat construction. Blow up modes of orbifold singularities can be considered as coordinates of complexified Kahler moduli space. Quiver diagrams are used for discribing D-branes near orbifold point. In this case it is possible to calculate Euler character defined through $\mbox{Ext}^i(A,B)$ groups and coherent sheaves $A, B$ over projective space, which are representations of orbifold space after blowing up procedure. These fractional sheaves are characterized by $Q_0$, $Q_2$ and $Q_4$ Ramon-Ramon charges, which have special type, calculated for $C^3/Z_3$ case. BPS central charge for $C^3/Z_3$ orbifold is calculated through Ramon-Ramon charges and Picard-Fuchs periods.


title: Compact inverse categories

abstract: The Ehresmann-Schein-Nambooripad theorem gives a structure theorem for inverse monoids: they are inductive groupoids. A particularly nice case due to Jarek is that commutative inverse monoids become semilattices of abelian groups. It has also been categorified by DeWolf-Pronk to a structure theorem for inverse categories as locally complete inductive groupoids. We show that in the case of compact inverse categories, this takes the particularly nice form of a semilattice of compact groupoids. Moreover, one-object compact inverse categories are exactly commutative inverse monoids. Compact groupoids, in turn, are determined in particularly simple terms of 3-cocycles by Baez-Lauda.


title: Transport in disordered systems: the single big jump approach

abstract: In a growing number of strongly disordered and dense systems, the dynamics of a particle pulled by an external force field exhibits super-diffusion. In the context of glass forming systems, super cooled glasses and contamination spreading in porous medium it was suggested to model this behavior with a biased continuous time random walk. Here we analyze the plume of particles far lagging behind the mean, with the single big jump principle. Revealing the mechanism of the anomaly, we show how a single trapping time, the largest one, is responsible for the rare fluctuations in the system. These non typical fluctuations still control the behavior of the mean square displacement, which is the most basic quantifier of the dynamics in many experimental setups. We show how the initial conditions, describing either stationary state or non-equilibrium case, persist for ever in the sense that the rare fluctuations are sensitive to the initial preparation. To describe the fluctuations of the largest trapping time, we modify Fr\'{e}chet's law from extreme value statistics, taking into consideration the fact that the large fluctuations are very different from those observed for independent and identically distributed random variables.


title: Schur ring and Codes for $S$-subgroups over $\Z_{2}^{n}$

abstract: In this paper the relationship between $S$-subgroups in $\Z_{2}^{n}$ and binary codes is shown. If the codes used are both $P(T)$-codes and $G$-codes, then the $S$-subgroup is free. The codes constructed are cyclic, decimated or symmetric and the $S$-subgroups obtained are free under the action the cyclic permutation subgroup, invariants under the action the decimated permutation subgroup and symmetric under the action of symmetric permutation subgroup, respectively. Also it is shows that there is no codes generating whole $\Z_{2}^{n}$ in any $\G_{n}(a)$-complete $S$-set of the $S$-ring $\mathfrak{S}(\Z_{2}^{n},S_{n})$.


title: Improving the communication for children with speech disorder using the smart toys

abstract: An attempt is made to develop a smart toy to help the children suffering with communication disorders. The children suffering with such disorders need additional attention and guidance to understand different types of social events and life activities. Various issues and features of the children with speech disorders are identified in this study and based on the inputs from the study, a working architecture is proposed with suitable policies. A prediction module with a checker component is designed in this work to produce alerts in at the time of abnormal behaviour of the child with communication disorder. The model is designed very sensitively to the behaviour of the child for a particular voice tone, based on which the smart toy will change to tones automatically. Such an arrangement proved to be helpful for the children to improve the communication with other due to the inclusion of continuous training for the smart toy from the prediction module.


title: SymNet: Symmetrical Filters in Convolutional Neural Networks

abstract: Symmetry is present in nature and science. In image processing, kernels for spatial filtering possess some symmetry (e.g. Sobel operators, Gaussian, Laplacian). Convolutional layers in artificial feed-forward neural networks have typically considered the kernel weights without any constraint. In this paper, we propose to investigate the impact of a symmetry constraint in convolutional layers for image classification tasks, taking our inspiration from the processes involved in the primary visual cortex and common image processing techniques. The goal is to assess the extent to which it is possible to enforce symmetrical constraints on the filters throughout the training process of a convolutional neural network (CNN) by modifying the weight update preformed during the backpropagation algorithm and to evaluate the change in performance. The main hypothesis of this paper is that the symmetrical constraint reduces the number of free parameters in the network, and it is able to achieve near identical performance to the modern methodology of training. In particular, we address the following cases: x/y-axis symmetry, point reflection, and anti-point reflection. The performance has been evaluated on four databases of images. The results support the conclusion that while random weights offer more freedom to the model, the symmetry constraint provides a similar level of performance while decreasing substantially the number of free parameters in the model. Such an approach can be valuable in phase-sensitive applications that require a linear phase property throughout the feature extraction process.


title: Growth Model Interpretation of Planet Size Distribution

abstract: The radii and orbital periods of 4000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. Yet an important question about the composition of planets ranging from 2 to 4 Earth radii still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multi-component, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 Earth masses, if half-ice and half-rock by mass, have radii of 2.5 Earth radii, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2-4 Earth radii range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are water worlds.


title: The $(\alpha, \beta)-$ramification invariants of a number field

abstract: Let $L$ be a number field. For a given prime $p$ we define integers $\alpha_{p}^{L}$ and $\beta_{p}^{L}$ with some interesting arithmetic properties. For instance, $\beta_{p}^{L}$ is equal to $1$ whenever $p$ does not ramify in $L$ and $\alpha_{p}^{L}$ is divisible by $p$ whenever $p$ is wildly ramified in $L$. The aforementioned properties, although interesting, follow easily from definitions; however a more interesting application of these invariants is the fact that they completely characterize the Dedekind zeta function of $L$. Moreover, if the residue class mod $p$ of $\alpha_{p}^{L}$ is not zero for all $p$ then such residues determine the genus of the integral trace.


title: Ultraviolet to Infrared Evolution and Nonperturbative Behavior of ${\rm SU}(N) \otimes {\rm SU}(N-4) \otimes {\rm U}(1)$ Chiral Gauge Theories

abstract: We analyze the ultraviolet to infrared evolution and nonperturbative properties of asymptotically free ${\rm SU}(N) \otimes {\rm SU}(N-4) \otimes {\rm U}(1)$ chiral gauge theories with $N_f$ copies of chiral fermions transforming according to $([2]_N,1)_{N-4} + ([\bar 1]_N,[\bar 1]_{N-4})_{-(N-2)} + (1,(2)_{N-4})_N$, where $[k]_N$ and $(k)_N$ denote the antisymmetric and symmetric rank-$k$ tensor representations of SU($N$) and the rightmost subscript is the U(1) charge. We give a detailed discussion for the lowest nondegenerate case, $N=6$. These theories can exhibit both self-breaking of a strongly coupled gauge symmetry and induced dynamical breaking of a weakly coupled gauge interaction symmetry due to fermion condensates produced by a strongly coupled gauge interaction. A connection with the dynamical breaking of ${\rm SU}(2)_L \otimes {\rm U}(1)_Y$ electroweak gauge symmetry by the quark condensates $\langle \bar q q\rangle$ due to color SU(3)$_c$ interactions is discussed. We also remark on direct-product chiral gauge theories with fermions in higher-rank tensor representations.


title: On the LoRa Modulation for IoT: Waveform Properties and Spectral Analysis

abstract: An important modulation technique for Internet of Things (IoT) is the one proposed by the LoRa allianceTM. In this paper we analyze the M-ary LoRa modulation in the time and frequency domains. First, we provide the signal description in the time domain, and show that LoRa is a memoryless continuous phase modulation. The cross-correlation between the transmitted waveforms is determined, proving that LoRa can be considered approximately an orthogonal modulation only for large M. Then, we investigate the spectral characteristics of the signal modulated by random data, obtaining a closed-form expression of the spectrum in terms of Fresnel functions. Quite surprisingly, we found that LoRa has both continuous and discrete spectra, with the discrete spectrum containing exactly a fraction 1/M of the total signal power.


title: Vector and axial-vector meson properties in a nonlocal SU(2) PNJL model

abstract: We study the features of a SU(2) Polyakov-Nambu-Jona-Lasinio model that includes wave function renormalization and nonlocal vector interactions. Within this framework we analyze, among other properties, the masses, width and decay constants of light vector and axial-vector mesons at finite temperature. Then we obtain the corresponding phase diagram in a finite density scenario, after characterizing the deconfinement and chiral restoration transitions.


title: Transport Triggered Array Processor for Vision Applications

abstract: Low-level sensory data processing in many Internet-of-Things (IoT) devices pursue energy efficiency by utilizing sleep modes or slowing the clocking to the minimum. To curb the share of stand-by power dissipation in those designs, near-threshold/sub-threshold operational points or ultra-low-leakage processes in fabrication are employed. Those limit the clocking rates significantly, reducing the computing throughputs of individual processing cores. In this contribution we explore compensating for the performance loss of operating in near-threshold region (Vdd =0.6V) through massive parallelization. Benefits of near-threshold operation and massive parallelism are optimum energy consumption per instruction operation and minimized memory roundtrips, respectively. The Processing Elements (PE) of the design are based on Transport Triggered Architecture. The fine grained programmable parallel solution allows for fast and efficient computation of learnable low-level features (e.g. local binary descriptors and convolutions). Other operations, including Max-pooling have also been implemented. The programmable design achieves excellent energy efficiency for Local Binary Patterns computations.


title: A physically-consistent, flexible and efficient strategy to convert local boundary conditions into nonlocal volume constraints

abstract: Nonlocal models provide exceptional simulation fidelity for a broad spectrum of scientific and engineering applications. However, wider deployment of nonlocal models is hindered by several modeling and numerical challenges. Among those, we focus on the nontrivial prescription of nonlocal boundary conditions, or volume constraints, that must be provided on a layer surrounding the domain where the nonlocal equations are posed. The challenge arises from the fact that, in general, data are provided on surfaces (as opposed to volumes) in the form of force or pressure data. In this paper we introduce an efficient, flexible and physically consistent technique for an automatic conversion of surface (local) data into volumetric data that does not have any constraints on the geometry of the domain and on the regularity of the nonlocal solution and that is not tied to any discretization. We show that our formulation is well-posed and that the limit of the nonlocal solution, as the nonlocality vanishes, is the local solution corresponding to the available surface data. Quadratic convergence rates are proved for the strong energy and L-2 convergence. We illustrate the theory with one dimensional numerical tests whose results provide the ground work for realistic simulations.


title: Polaron-transformed dissipative Lipkin-Meshkov-Glick Model

abstract: We investigate the Lipkin-Meshkov-Glick model coupled to a thermal bath. Since the isolated model itself exhibits a quantum phase transition, we explore the critical signatures of the open system. Starting from a system-reservoir interaction written in positive definite form, we find that the position of the critical point remains unchanged, in contrast to the popular mean-field prediction. Technically, we employ the polaron transform to be able to study the full crossover regime from the normal to the symmetry-broken phase, which allows us to investigate the fate of quantum-critical points subject to dissipative environments. The signatures of the phase transition are reflected in observables like stationary mode occupation or waiting-time distributions.


title: Examining Untempered Social Media: Analyzing Cascades of Polarized Conversations

abstract: Online social media, periodically serves as a platform for cascading polarizing topics of conversation. The inherent community structure present in online social networks (homophily) and the advent of fringe outlets like Gab have created online "echo chambers" that amplify the effects of polarization, which fuels detrimental behavior. Recently, in October 2018, Gab made headlines when it was revealed that Robert Bowers, the individual behind the Pittsburgh Synagogue massacre, was an active member of this social media site and used it to express his anti-Semitic views and discuss conspiracy theories. Thus to address the need of automated data-driven analyses of such fringe outlets, this research proposes novel methods to discover topics that are prevalent in Gab and how they cascade within the network. Specifically, using approximately 34 million posts, and 3.7 million cascading conversation threads with close to 300k users; we demonstrate that there are essentially five cascading patterns that manifest in Gab and the most "viral" ones begin with an echo-chamber pattern and grow out to the entire network. Also, we empirically show, through two models viz. Susceptible-Infected and Bass, how the cascades structurally evolve from one of the five patterns to the other based on the topic of the conversation with upto 84% accuracy.


title: Tests of Acoustic Scale Shifts in Halo-based Mock Galaxy Catalogues

abstract: We utilise mock catalogues from high-accuracy cosmological $N$-body simulations to quantify shifts in the recovery of the acoustic scale that could potentially result from galaxy clustering bias. The relationship between galaxies and dark matter halos presents a complicated source of systematic errors in modern redshift surveys, particularly when aiming to make cosmological measurements to sub-percent precision. Apart from a scalar, linear bias parameter accounting for the density contrast ratio between matter tracers and the true matter distribution, other types of galaxy bias, such as assembly and velocity biases, may also significantly alter clustering signals from small to large scales. We create mocks based on generalised halo occupation populations of 36 periodic boxes from the \abacuscosmos release with. In a total volume of $48 \, h^{-3} \mathrm{Gpc}^3$, we test various biased models along with an unbiased base case. Two reconstruction methods are applied to galaxy samples and the apparent acoustic scale is derived by fitting the two-point correlation function multipoles. With respect to the baseline, we find a $0.3\%$ shift in the line-of-sight acoustic scale for one variation in the satellite galaxy population, and we find an $0.7\%$ shift for an extreme level of velocity bias of the central galaxies. All other bias models are consistent with zero shift at the $0.2\%$ level after reconstruction. We note that the bias models explored are relatively large variations, producing sizeable and likely distinguishable changes in small-scale clustering, the modelling of which would further calibrate the BAO standard ruler.


title: Feedback Linearization for Quadrotors UAV

abstract: In the paper "Control Design for UAV Quadrotors via Embedded Model Control" [1], the authors designed a complete control unit for a UAV Quadrotor, based on the Embedded Model Control (EMC) methodology, in combination with the Feedback Linearization (FL); when applied to non-linear systems. Specifically, [1] proposes to use the FL as a novel way to design the internal model for the EMC state and disturbance predictor. To support the treatise in [1], in this report the feedback-linearized model of the UAV quadrotor leveraged in [1] is step-by-step derived.


title: Optimal In-field Routing for Full and Partial Field Coverage with Arbitrary Non-Convex Fields and Multiple Obstacle Areas

abstract: Within the context of optimising the logistics in agriculture this paper relates to optimal in-field routing for full and partial field coverage with arbitrary non-convex fields and multiple obstacle areas. It is distinguished between nine different in-field routing tasks: two for full-field coverage, seven for partial-field coverage and one for shortest path planning between any two vertices of the transition graph. It differentiates between equal or different start and end vertices for a task, coverage of only a subset of vertices, and a subset of edges or combinations. The proposed methods are developed primarily for applying sprays and fertilisers with larger operating widths and with fields where there is unique headland path. Partial field coverage where, e.g., only a specific subset of edges has to be covered is relevant for precision agriculture and also for optimised logistical operation of smaller-sized machinery with limited loading capacities. The result of this research is the proposition of two compatible algorithms for optimal full and partial field coverage path planning, respectively. These are evaluated on three real-world fields to demonstrate their characteristics and computational efficiency.


title: Radiation pressure clear-out of dusty photoevaporating discs

abstract: Theoretical models of protoplanetary disc dispersal predict a phase where photoevaporation has truncated the disc at several AU, creating a pressure trap which is dust-rich. Previous models predicted this phase could be long-lived (~Myr), contrary to the observational constraints. We show that dust in the pressure trap can be removed from the disc by radiation pressure exerting a significant acceleration, and hence radial velocity, on small dust particles that reside in the surface layers of the disc. The dust in the pressure trap is not subject to radial drift so it can grow to reach sizes large enough to fragment. Hence small particles removed from the surface layers are replaced by the fragments of larger particles. This link means radiation pressure can deplete the dust at all particle sizes. Through a combination of 1D and 2D models, along with secular models that follow the disc's long-term evolution, we show that radiation pressure can deplete dust from pressure traps created by photoevaporation in ~1e5 years, while the photoevaporation created cavity still resides at 10s of AU. After this phase of radiation pressure removal of dust, the disc is gas-rich and dust depleted and radially optically thin to stellar light, having observational signatures similar to a gas-rich, young debris disc. Indeed many of the young stars (~<10 Myr old) classified as hosting a debris disc may rather be discs that have undergone this process.


title: Phylogenetic correlations can suffice to infer protein partners from sequences

abstract: Determining which proteins interact together is crucial to a systems-level understanding of the cell. Recently, algorithms based on Direct Coupling Analysis (DCA) pairwise maximum-entropy models have allowed to identify interaction partners among paralogous proteins from sequence data. This success of DCA at predicting protein-protein interactions could be mainly based on its known ability to identify pairs of residues that are in contact in the three-dimensional structure of protein complexes and that coevolve to remain physicochemically complementary. However, interacting proteins possess similar evolutionary histories. What is the role of purely phylogenetic correlations in the performance of DCA-based methods to infer interaction partners? To address this question, we employ controlled synthetic data that only involve phylogeny and no interactions or contacts. We find that DCA accurately identifies the pairs of synthetic sequences that share evolutionary history. While phylogenetic correlations confound the identification of contacting residues by DCA, they are thus useful to predict interacting partners among paralogs. We find that DCA performs as well as phylogenetic methods to this end, and slightly better than them with large and accurate training sets. Employing DCA or phylogenetic methods within an Iterative Pairing Algorithm (IPA) allows to predict pairs of evolutionary partners without a training set. We demonstrate the ability of these various methods to correctly predict pairings among real paralogous proteins with genome proximity but no known physical interaction, illustrating the importance of phylogenetic correlations in natural data. However, for physically interacting and strongly coevolving proteins, DCA and mutual information outperform phylogenetic methods. We discuss how to distinguish physically interacting proteins from those only sharing evolutionary history.

categories: q-bio.BM cond-mat.stat-mech physics.bio-ph q-bio.QM


title: Beyond Folklore: A Scaling Calculus for the Design and Initialization of ReLU Networks

abstract: We propose a system for calculating a "scaling constant" for layers and weights of neural networks. We relate this scaling constant to two important quantities that relate to the optimizability of neural networks, and argue that a network that is "preconditioned" via scaling, in the sense that all weights have the same scaling constant, will be easier to train. This scaling calculus results in a number of consequences, among them the fact that the geometric mean of the fan-in and fan-out, rather than the fan-in, fan-out, or arithmetic mean, should be used for the initialization of the variance of weights in a neural network. Our system allows for the off-line design & engineering of ReLU neural networks, potentially replacing blind experimentation.


title: The study of unclassified B[e] stars and candidates in the Galaxy and Magellanic Clouds

abstract: We investigated 12 unclassified B[e] stars or candidates, 8 from the Galaxy, 2 from the Large Magellanic Cloud (LMC) and 2 from the Small Magellanic Cloud (SMC). Based on the analysis of high-resolution spectroscopic (FEROS) and photometric data, we confirmed the presence of the B[e] phenomenon for all objects of our sample, except for one (IRAS 07455-3143). We derived their effective temperature, spectral type, luminosity class, interstellar extinction and, using the distances from Gaia DR2, we obtained their bolometric magnitude, luminosity and radius. Modeling of the forbidden lines present in the FEROS spectra revealed information about the kinematics and geometry of the circumstellar medium of these objects. In addition, we analyzed the light curves of four stars, finding their most probable periods. The evolutionary stage of 11 stars of our sample is suggested from their position on the HR diagram, taking into account evolutionary tracks of stars with solar, LMC and SMC metallicities. As results, we identified B and B[e] supergiants, B[e] stars probably at the main sequence or close to its end, post-AGB and HAeB[e] candidates, and A[e] stars in the main sequence or in the pre-main sequence. However, our most remarkable results are the identification of the third A[e] supergiant (ARDB\,54, the first one in the LMC), and of an "LBV impostor" in the SMC (LHA 115-N82).


title: Analysis of the susceptible-infected-susceptible epidemic dynamics in networks via the non-backtracking matrix

abstract: We study the stochastic susceptible-infected-susceptible model of epidemic processes on finite directed and weighted networks with arbitrary structure. We present a new lower bound on the exponential rate at which the probabilities of nodes being infected decay over time. This bound is directly related to the leading eigenvalue of a matrix that depends on the non-backtracking and incidence matrices of the network. The dimension of this matrix is N+M, where N and M are the number of nodes and edges, respectively. We show that this new lower bound improves on an existing bound corresponding to the so-called quenched mean-field theory. Although the bound obtained from a recently developed second-order moment-closure technique requires the computation of the leading eigenvalue of an N^2 x N^2 matrix, we illustrate in our numerical simulations that the new bound is tighter, while being computationally less expensive for sparse networks. We also present the expression for the corresponding epidemic threshold in terms of the adjacency matrix of the line graph and the non-backtracking matrix of the given network.


title: Pure entropic regularization for metrical task systems

abstract: We show that on every $n$-point HST metric, there is a randomized online algorithm for metrical task systems (MTS) that is $1$-competitive for service costs and $O(\log n)$-competitive for movement costs. In general, these refined guarantees are optimal up to the implicit constant. While an $O(\log n)$-competitive algorithm for MTS on HST metrics was developed by Bubeck et al. (SODA 2019), that approach could only establish an $O((\log n)^2)$-competitive ratio when the service costs are required to be $O(1)$-competitive. Our algorithm can be viewed as an instantiation of online mirror descent with the regularizer derived from a multiscale conditional entropy. In fact, our algorithm satisfies a set of even more refined guarantees; we are able to exploit this property to combine it with known random embedding theorems and obtain, for any $n$-point metric space, a randomized algorithm that is $1$-competitive for service costs and $O((\log n)^2)$-competitive for movement costs.


title: Exact WKB and abelianization for the $T_3$ equation

abstract: We describe the exact WKB method from the point of view of abelianization, both for Schr\"odinger operators and for their higher-order analogues (opers). The main new example which we consider is the "$T_3$ equation," an order $3$ equation on the thrice-punctured sphere, with regular singularities at the punctures. In this case the exact WKB analysis leads to consideration of a new sort of Darboux coordinate system on a moduli space of flat $\mathrm{SL}(3)$-connections. We give the simplest example of such a coordinate system, and verify numerically that in these coordinates the monodromy of the $T_3$ equation has the expected asymptotic properties. We also briefly revisit the Schr\"odinger equation with cubic potential and the Mathieu equation from the point of view of abelianization.