title: Simulation of Higher-Order Topological Phases and Related Topological Phase Transitions in a Superconducting Qubit

abstract: Higher-order topological phases give rise to new bulk and boundary physics, as well as new classes of topological phase transitions. While the realization of higher-order topological phases has been confirmed in many platforms by detecting the existence of gapless boundary modes, a direct determination of the higher-order topology and related topological phase transitions through the bulk in experiments has still been lacking. To bridge the gap, in this work we carry out the simulation of a two-dimensional second-order topological phase in a superconducting qubit. Owing to the great flexibility and controllability of the quantum simulator, we observe the realization of higher-order topology directly through the measurement of the pseudo-spin texture in momentum space of the bulk for the first time, in sharp contrast to previous experiments based on the detection of gapless boundary modes in real space. Also through the measurement of the evolution of pseudo-spin texture with parameters, we further observe novel topological phase transitions from the second-order topological phase to the trivial phase, as well as to the first-order topological phase with nonzero Chern number. Our work sheds new light on the study of higher-order topological phases and topological phase transitions.


title: Infinite-fold enhancement in communications capacity using pre-shared entanglement

abstract: Pre-shared entanglement can significantly boost communication rates in the regime of high thermal noise, and a low-brightness transmitter. In this regime, the ratio between the entanglement-assisted capacity and the Holevo capacity, the maximum reliable-communication rate permitted by quantum mechanics without any pre-shared entanglement as a resource, is known to scale as $\log(1/N_S)$, where $N_S \ll 1$ is the mean transmitted photon number per mode. This is especially promising in enabling a large boost to radio-frequency communications in the weak-transmit-power regime, by exploiting pre-shared optical-frequency entanglement, e.g., distributed by the quantum internet. In this paper, we propose a structured design of a quantum transmitter and receiver that leverages continuous-variable pre-shared entanglement from a downconversion source, which can harness this purported infinite-fold capacity enhancement---a problem open for over a decade. Finally, the implication of this result to the breaking of the well-known {\em square-root law} for covert communications, with pre-shared entanglement assistance, is discussed.


title: A multi-country dynamic factor model with stochastic volatility for euro area business cycle analysis

abstract: This paper develops a dynamic factor model that uses euro area (EA) country-specific information on output and inflation to estimate an area-wide measure of the output gap. Our model assumes that output and inflation can be decomposed into country-specific stochastic trends and a common cyclical component. Comovement in the trends is introduced by imposing a factor structure on the shocks to the latent states. We moreover introduce flexible stochastic volatility specifications to control for heteroscedasticity in the measurement errors and innovations to the latent states. Carefully specified shrinkage priors allow for pushing the model towards a homoscedastic specification, if supported by the data. Our measure of the output gap closely tracks other commonly adopted measures, with small differences in magnitudes and timing. To assess whether the model-based output gap helps in forecasting inflation, we perform an out-of-sample forecasting exercise. The findings indicate that our approach yields superior inflation forecasts, both in terms of point and density predictions.


title: Broadcasting Competitively against Adaptive Adversary in Multi-channel Radio Networks

abstract: Broadcasting in wireless networks is vulnerable to adversarial jamming. To thwart such behavior, \emph{resource competitive analysis} is proposed. In this framework, sending, listening, or jamming on one channel for one time slot costs one unit of energy. The adversary can employ arbitrary strategy to disrupt communication, but has a limited energy budget $T$. The honest nodes, on the other hand, aim to accomplish broadcast while spending only $o(T)$. Previous work has shown, in a $C$-channels network containing $n$ nodes, for large $T$ values, each node can receive the message in $\tilde{O}(T/C)$ time, while spending only $\tilde{O}(\sqrt{T/n})$ energy. However, these multi-channel algorithms only work for certain values of $n$ and $C$, and can only tolerate an oblivious adversary. In this work, we provide new upper and lower bounds for broadcasting in multi-channel radio networks, from the perspective of resource competitiveness. Our algorithms work for arbitrary $n,C$ values, require minimal prior knowledge, and can tolerate a powerful adaptive adversary. More specifically, in our algorithms, for large $T$ values, each node's runtime is $O(T/C)$, and each node's energy cost is $\tilde{O}(\sqrt{T/n})$. We also complement algorithmic results with lower bounds, proving both the time complexity and the energy complexity of our algorithms are optimal or near-optimal (within a poly-log factor). Our technical contributions lie in using "epidemic broadcast" to achieve time efficiency and resource competitiveness, and employing coupling techniques in the analysis to handle the adaptivity of the adversary. At the lower bound side, we first derive a new energy complexity lower bound for 1-to-1 communication in the multi-channel setting, and then apply simulation and reduction arguments to obtain the desired result.


title: Simulated Blockchains for Machine Learning Traceability and Transaction Values in the Monero Network

abstract: Monero is a popular crypto-currency which focuses on privacy. The blockchain uses cryptographic techniques to obscure transaction values as well as a `ring confidential transaction' which seeks to hide a real transaction among a variable number of spoofed transactions. We have developed training sets of simulated blockchains of 10 and 50 agents, for which we have control over the ground truth and keys, in order to test these claims. We featurize Monero transactions by characterizing the local structure of the public-facing blockchains and use labels obtained from the simulations to perform machine learning. Machine Learning of our features on the simulated blockchain shows that the technique can be used to aide in identifying individuals and groups, although it did not successfully reveal the hidden transaction values. We apply the technique on the real Monero blockchain to identify ShapeShift transactions, a cryptocurrency exchange that has leaked information through their API providing labels for themselves and their users.


title: Edge ideals with almost maximal finite index and their powers

abstract: A graded ideal $I$ in $\mathbb{K}[x_1,\ldots,x_n]$, where $\mathbb{K}$ is a field, is said to have almost maximal finite index if its minimal free resolution is linear up to the homological degree $\mathrm{pd}(I)-2$, while it is not linear at the homological degree $\mathrm{pd}(I)-1$, where $\mathrm{pd}(I)$ denotes the projective dimension of $I$. In this paper we classify the graphs whose edge ideals have this property. This in particular shows that for edge ideals the property of having almost maximal finite index does not depend on the characteristic of $\mathbb{K}$. We also compute the non-linear Betti numbers of these ideals. Finally, we show that for the edge ideal $I$ of a graph $G$ with almost maximal finite index, the ideal $I^s$ has a linear resolution for $s\geq 2$ if and only if the complementary graph $\bar{G}$ does not contain induced cycles of length $4$.


title: Enhanced Microgrid Power Flow Incorporating Hierarchical Control

abstract: An enhanced microgrid power flow (EMPF) is devised to incorporate hierarchical control effects. The new contributions are threefold: 1) an advanced-hierarchical-control-based Newton approach is established to accurately assess power sharing and voltage regulation effects; 2) a modified Jacobian matrix is derived to incorporate droop control and various secondary control modes; and 3) the secondary adjustment is calculated on top of the droop-control-based power flow results to ensure a robust Newton solution. Case studies validate that EMPF is efficacious and efficient and can serve as a powerful tool for microgrid operation and monitoring, especially for those highly meshed microgrids in urban areas.


title: Holomorphic Explosions

abstract: This paper concerns the problem of extending the parameter domain of holomorphic motions to include isolated boundary points, punctures of the domains. Supposing that the parameter domain has an isolated boundary point $\lambda^*$, we explore the extension properties of the holomorphic motion to $\lambda^*$.


title: Proof of Sun's supercongruence involving Catalan numbers

abstract: We confirm a conjectural supercongruence involving Catalan numbers, which is one of the 100 selected open conjectures on congruences of Sun. The proof makes use of hypergeometric series identities and symbolic summation method.


title: Securing Bring-Your-Own-Device (BYOD) Programming Exams

abstract: Traditional pen and paper exams are inadequate for modern university programming courses as they are misaligned with pedagogies and learning objectives that target practical coding ability. Unfortunately, many institutions lack the resources or space to be able to run assessments in dedicated computer labs. This has motivated the development of bring-your-own-device (BYOD) exam formats, allowing students to program in a similar environment to how they learnt, but presenting instructors with significant additional challenges in preventing plagiarism and cheating. In this paper, we describe a BYOD exam solution based on lockdown browsers, software which temporarily turns students' laptops into secure workstations with limited system or internet access. We combine the use of this technology with a learning management system and cloud-based programming tool to facilitate conceptual and practical programming questions that can be tackled in an interactive but controlled environment. We reflect on our experience of implementing this solution for a major undergraduate programming course, highlighting our principal lesson that policies and support mechanisms are as important to consider as the technology itself.


title: The case for a Casimir cosmology

abstract: The cosmological constant, also known as dark energy, was believed to be caused by vacuum fluctuations, but naive calculations give results in stark disagreement with fact. In the Casimir effect, vacuum fluctuations cause forces in dielectric media, which is very well described by Lifshitz theory. Recently, using the analogy between geometries and media, a cosmological constant of the correct order of magnitude was calculated with Lifshitz theory [U. Leonhardt, Ann. Phys. (New York) 411, 167973 (2019)]. This paper discusses the empirical evidence and the ideas behind the Lifshitz theory of the cosmological constant without requiring prior knowledge of cosmology and quantum field theory.


title: The Proximal Method of Multipliers for a Class of Nonsmooth Convex Optimization

abstract: This paper develops the proximal method of multipliers for a class of nonsmooth convex optimization. The method generates a sequence of minimization problems (subproblems). We show that the sequence of approximations to the solutions of the subproblems converges to a saddle point of the Lagrangian even if the original optimization problem may possess multiple solutions. The augmented Lagrangian due to Fortin appears in the subproblem. The remarkable property of the augmented Lagrangian over the standard Lagrangian is that it is always differentiable, and it is often semismoothly differentiable. This fact allows us to employ a nonsmooth Newton method for computing an approximation to the subproblem. The proximal term serves as the regularization of the objective function and guarantees the solvability of the Newton system without assuming strong convexity on the objective function. We exploit the theory of the nonsmooth Newton method to provide a rigorous proof for the global convergence of the proposed algorithm.


title: Commonly Knowing Whether

abstract: This paper introduces the notion of `commonly knowing whether', a non-standard version of standard common knowledge which is defined on the basis of `knowing whether', instead of standard `knowing that'. After giving five possible definitions of this notion, we explore the logical relations among them in the single-agent and multi-agent cases. We propose a sound and complete axiomatization. We investigate one of the five definitions in terms of expressivity via a strategy of modal comparison games.


title: Communications-Caching-Computing Tradeoff Analysis for Bidirectional Data Computation in Mobile Edge Networks

abstract: With the advent of the modern mobile traffic, e.g., online gaming, augmented reality delivery and etc., a novel bidirectional computation task model where the input data of each task consists of two parts, one generated at the mobile device in real-time and the other originated from the Internet proactively, is emerging as an important use case of 5G. In this paper, for ease of analytical analysis, we consider the homogeneous bidirectional computation task model in a mobile edge network which consists of one mobile edge computing (MEC) server and one mobile device, both enabled with computing and caching capabilities. Each task can be served via three mechanisms, i.e., local computing with local caching, local computing without local caching and computing at the MEC server. To minimize the average bandwidth, we formulate the joint caching and computing optimization problem under the latency, cache size and average power constraints. We derive the closed-form expressions for the optimal policy and the minimum bandwidth. The tradeoff among communications, computing and caching is illustrated both analytically and numerically, which provides insightful guideline for the network designers.


title: The Polyakov loop dependence of bulk viscosity of QCD matter

abstract: In this work, we show the dependence of bulk viscosity on Polyakov loop in 3+1 dimensional topologically massive model (TMM). This model contains equally massive non-Abelian gauge fields without spontaneous symmetry breaking. In earlier works, the bulk viscosity was found from the trace anomaly in massless $\phi^4$ model and Yang-Mills (YM) theory and its dependence on the quantum corrections was established. In TMM, the trace anomaly is absent due to the presence of kinetic term of a two-form field $B$ in the action. This model also provides the dependence of bulk viscosity on the mass of the gauge bosons. The mass of the gauge bosons in TMM acts as magnetic mass in the perturbative thermal field theory. This magnetic mass is gauge independent unlike what is found in massless YM theory. The ratio of bulk viscosity and entropy density vs temperature is plotted which shows the required behaviour near critical point. We also observe that the strong coupling constant has the same behaviour at high energy limit (i.e. asymptotic freedom) as that of massless YM theory at zero temperature.


title: Controlled transitions between phyllotactic states of repulsive particles confined on the surface of a cylinder

abstract: Phyllotactic states are regular lattice-like structures on cylinders and are a botanical classification scheme. In this communication, we report a sequence of transitions between phyllotactic states for particles with a repulsive particle-particle interaction on a cylindrical geometry at zero temperature. We can infer the transition points as a function of density via Monte Carlo simulations, as well as the mathematical descriptions of the ground states. The lattices we generate are described as phyllotactic states that fit onto the cylindrical surface as a set of helical chains. Our analysis shows how all state energies lie on the same parabola which we exploit to find the transitions.


title: Dispersal of protoplanetary disks by the combination of magnetically driven and photoevaporative winds

abstract: We investigate the roles of magnetically driven disk wind (MDW) and thermally driven photoevaporative wind (PEW) in the long-time evolution of protoplanetary disks. We start simulations from the early phase in which the disk mass is $0.118\,{\mathrm{M}_{\odot}}$ around a $1\,{\mathrm{M}_{\odot}}$ star and track the evolution until the disk is completely dispersed. We incorporate the mass loss by PEW and the mass loss and magnetic braking (wind torque) by MDW, in addition to the viscous accretion, viscous heating, and stellar irradiation. We find that MDW and PEW respectively have different roles: magnetically driven wind ejects materials from an inner disk in the early phase, whereas photoevaporation has a dominant role in the late phase in the outer ($\gtrsim1\,$au) disk. The disk lifetime, which depends on the combination of MDW, PEW, and viscous accretion, shows a large variation of $\sim1$-$20\,$Myr; the gas is dispersed mainly by the MDW and the PEW in the cases with a low viscosity and the lifetime is sensitive to the mass-loss rate and torque of the MDW, whereas the lifetime is insensitive to these parameters when the viscosity is high. Even in disks with very weak turbulence, the cooperation of MDW and PEW enables the disk dispersal within a few Myr.


title: The sharp higher order Lorentz--Poincar\'e and Lorentz--Sobolev inequalities in the hyperbolic spaces

abstract: In this paper, we study the sharp Poincar\'e inequality and the Sobolev inequalities in the higher order Lorentz--Sobolev spaces in the hyperbolic spaces. These results generalize the ones obtained in \cite{Nguyen2020a} to the higher order derivatives and seem to be new in the context of the Lorentz--Sobolev spaces defined in the hyperbolic spaces.


title: Towards Highly Efficient State Estimation with Nonlinear Measurements in Distribution Systems

abstract: This letter proposes a novel and highly efficient distribution system state estimation (DSSE) algorithm with nonlinear measurements from supervisory control and data acquisition (SCADA) systems. Conventional DSSE, i.e., a weighted least square (WLS)-based method, requires multiple Monte Carlo simulations for overall accuracy evaluation and has high calculation cost due to the nonlinear iterative process. The proposed method uses the Taylor series of voltages for constructing a linear DSSE model in the interval form and then solve this model by interval arithmetic. Compared with the nonlinear WLS-based method, the proposed method obtains accurate and robust estimates via a single random sampling of measurements and is computationally efficient. The comparative analysis of the IEEE 34-bus distribution system points to improved estimation results relative to the WLS-based method.


title: Channel Assignment in Uplink Wireless Communication using Machine Learning Approach

abstract: This letter investigates a channel assignment problem in uplink wireless communication systems. Our goal is to maximize the sum rate of all users subject to integer channel assignment constraints. A convex optimization based algorithm is provided to obtain the optimal channel assignment, where the closed-form solution is obtained in each step. Due to high computational complexity in the convex optimization based algorithm, machine learning approaches are employed to obtain computational efficient solutions. More specifically, the data are generated by using convex optimization based algorithm and the original problem is converted to a regression problem which is addressed by the integration of convolutional neural networks (CNNs), feed-forward neural networks (FNNs), random forest and gated recurrent unit networks (GRUs). The results demonstrate that the machine learning method largely reduces the computation time with slightly compromising of prediction accuracy.


title: Asymptotic inversion of the binomial and negative binomial cumulative distribution functions

abstract: The computation and inversion of the binomial and negative binomial cumulative distribution functions play a key role in many applications. In this paper, we explain how methods used for the central beta distribution function (described in [2]) can be used to obtain asymptotic representations of these functions, and also for their inversion. The performance of the asymptotic inversion methods is illustrated with numerical examples.


title: Pseudo-Nambu-Goldstone dark matter from gauged $U(1)_{B-L}$ symmetry

abstract: A pseudo-Nambu-Goldstone boson (pNGB) is an attractive candidate for dark matter since the current severe limits of dark matter direct detection experiments are naturally evaded by its nature. We construct a model with pNGB dark matter based on a {\it gauged} $U(1)_{B-L}$ symmetry, where no ad-hoc global symmetry is assumed. The model keeps natural suppression mechanism for the dark matter direct detection cross section. On the other hand, the pNGB can decay through the new high scale suppressed operators. We show that the pNGB has long enough lifetime to be a dark matter in the wide range of the parameter space of the model. The thermal relic abundance of pNGB dark matter can be fit with the observed value against the constraints on the dark matter decays from the cosmic-ray observations.


title: Aggregated Learning: A Vector-Quantization Approach to Learning Neural Network Classifiers

abstract: We consider the problem of learning a neural network classifier. Under the information bottleneck (IB) principle, we associate with this classification problem a representation learning problem, which we call "IB learning". We show that IB learning is, in fact, equivalent to a special class of the quantization problem. The classical results in rate-distortion theory then suggest that IB learning can benefit from a "vector quantization" approach, namely, simultaneously learning the representations of multiple input objects. Such an approach assisted with some variational techniques, result in a novel learning framework, "Aggregated Learning", for classification with neural network models. In this framework, several objects are jointly classified by a single neural network. The effectiveness of this framework is verified through extensive experiments on standard image recognition and text classification tasks.


title: Interpretable feature subset selection: A Shapley value based approach

abstract: For feature selection and related problems, we introduce the notion of classification game, a cooperative game, with features as players and hinge loss based characteristic function and relate a feature's contribution to Shapley value based error apportioning (SVEA) of total training error. Our major contribution is ($\star$) to show that for any dataset the threshold 0 on SVEA value identifies feature subset whose joint interactions for label prediction is significant or those features that span a subspace where the data is predominantly lying. In addition, our scheme ($\star$) identifies the features on which Bayes classifier doesn't depend but any surrogate loss function based finite sample classifier does; this contributes to the excess $0$-$1$ risk of such a classifier, ($\star$) estimates unknown true hinge risk of a feature, and ($\star$) relate the stability property of an allocation and negative valued SVEA by designing the analogue of core of classification game. Due to Shapley value's computationally expensive nature, we build on a known Monte Carlo based approximation algorithm that computes characteristic function (Linear Programs) only when needed. We address the potential sample bias problem in feature selection by providing interval estimates for SVEA values obtained from multiple sub-samples. We illustrate all the above aspects on various synthetic and real datasets and show that our scheme achieves better results than existing recursive feature elimination technique and ReliefF in most cases. Our theoretically grounded classification game in terms of well defined characteristic function offers interpretability (which we formalize in terms of final task) and explainability of our framework, including identification of important features.


title: Interpretation of the modulus spectra of organic field-effect transistors with electrode overlap and peripheral regions: determination of the electronic properties of the gate insulator and organic semiconductor

abstract: The modulus spectra of organic field-effect transistors (OFETs) with electrode overlap and peripheral regions have been experimentally and theoretically investigated. The complex impedance of regioregular poly(3-hexylthiophene-2,5-diyl) (P3HT) OFETs with electrode overlap and peripheral regions was measured with a frequency response analyzer. The complex modulus was derived from an equivalent circuit of OFETs with overlap and peripheral regions using a four-terminal matrix approach. The modulus spectra of the P3HT OFETs were successfully fitted by those calculated using the expression derived from the equivalent circuit. Three structures were found in the modulus spectra of the P3HT OFETs owing to the dielectric properties of the gate insulator, transport properties of the organic semiconductor, and contact resistance from the low to high frequency ranges. The resistivity of the gate insulators and the field-effect mobility of working OFETs were determined using the values of the circuit components of the equivalent circuit obtained by fitting.


title: Lyapunov spectrum properties and continuity of the lower joint spectral radius

abstract: We study ergodic optimization and multifractal behavior of Lyapunov exponents for matrix cocycles. We show the continuity of the entropy spectrum at the boundary of Lyapunov spectrum in the sense that $h_{top}(E(\alpha_{t}))\ \rightarrow h_{top}(E(\beta(\mathcal{A}))$ for generic cocycles, where $E(\alpha)=\{x\in X: \lim_{n\rightarrow \infty}\frac{1}{n}\log \|\mathcal{A}^{n}(x)\|=\alpha\}$. We also show that the Lyapunov spectrum is equal to the closure of the set where the entropy spectrum is positive for such cocycles over mixing subshifts of finite type. Moreover, we prove the restricted variational principle for such cocycles. We prove the continuity of the lower joint spectral radius for general cocycles under the assumption that linear cocycles satisfy a cone condition.


title: Average AoI in Multi-Source Systems with Source-Aware Packet Management

abstract: We study the information freshness under three different source aware packet management policies in a status update system consisting of two independent sources and one server. The packets of each source are generated according to the Poisson process and the packets are served according to an exponentially distributed service time. We derive the average age of information (AoI) of each source using the stochastic hybrid systems (SHS) technique for each packet management policy. In Policy 1, the queue can contain at most two waiting packets at the same time (in addition to the packet under service), one packet of source 1 and one packet of source 2. When the server is busy at an arrival of a packet, the possible packet of the same source waiting in the queue (hence, source-aware) is replaced by the arrived fresh packet. In Policy 2, the system (i.e., the waiting queue and the server) can contain at most two packets, one from each source. When the server is busy at an arrival of a packet, the possible packet of the same source in the system is replaced by the fresh packet. Policy 3 is similar to Policy 2 but it does not permit preemption in service, i.e., while a packet is under service all new arrivals from the same source are blocked and cleared. Numerical results are provided to assess the fairness between sources and the sum average AoI of the proposed policies.


title: Attention Flow: End-to-End Joint Attention Estimation

abstract: This paper addresses the problem of understanding joint attention in third-person social scene videos. Joint attention is the shared gaze behaviour of two or more individuals on an object or an area of interest and has a wide range of applications such as human-computer interaction, educational assessment, treatment of patients with attention disorders, and many more. Our method, Attention Flow, learns joint attention in an end-to-end fashion by using saliency-augmented attention maps and two novel convolutional attention mechanisms that determine to select relevant features and improve joint attention localization. We compare the effect of saliency maps and attention mechanisms and report quantitative and qualitative results on the detection and localization of joint attention in the VideoCoAtt dataset, which contains complex social scenes.


title: Local stationarity of exponential last passage percolation

abstract: We consider point to point last passage times to every vertex in a neighbourhood of size $\delta N^{\frac{2}{3}}$, distance $N$ away from the starting point. The increments of these last passage times in this neighbourhood are shown to be jointly equal to their stationary versions with high probability that depends on $\delta$ only. With the help of this result we show that 1) the $\text{Airy}_2$ process is locally close to a Brownian motion in total variation; 2) the tree of point to point geodesics starting from every vertex in a box of side length $\delta N^{\frac{2}{3}}$ going to a point at distance $N$ agree inside the box with the tree of infinite geodesics going in the same direction; 3) two geodesics starting from $N^{\frac{2}{3}}$ away from each other, to a point at distance $N$ will not coalesce too close to either endpoints on the macroscopic scale. Our main results rely on probabilistic methods only.


title: Computational Hardness of Multidimensional Subtraction Games

abstract: We study algorithmic complexity of solving subtraction games in a~fixed dimension with a finite difference set. We prove that there exists a game in this class such that any algorithm solving the game runs in exponential time. Also we prove an existence of a game in this class such that solving the game is PSPACE-hard. The results are based on the construction introduced by Larsson and W\"astlund. It relates subtraction games and cellular automata.


title: Properties and applications of some special integer number sequences

abstract: In this paper, we provide properties and applications of some special integer sequences. We generalize and give some properties of Pisano period. Moreover, we provide a new application in Cryptography and applications of some quaternion elements.


title: On the maintenance of an attached leading-edge vortex via model bird alula

abstract: Researchers have hypothesized that the post-stall lift benefit of bird's alular feathers, or alula, stems from the maintenance of an attached leading-edge vortex (LEV) over their thin-profiled, outer hand-wing. Here, we investigate the connection between the alula and LEV attachment via flow measurements in a wind tunnel. We show that a model alula, whose wetted area is 1% that of the wing, stabilizes a recirculatory aft-tilted LEV on a steadily-translating unswept wing inclined at post-stall incidences. The attached vortex is the result of the alula's ability to smoothly merge otherwise separate leading- and side-edge vortical flows. We identify two key processes that facilitate this merging: i) the steering of spanwise vorticity generated at the wing's leading edge back to the wing plane and ii) an aft-located wall-jet of high-magnitude root-to-tip spanwise flow (>80% that of the freestream velocity). The latter feature induces LEV roll-up while the former feature tilts LEV vorticity aft and evacuates this flow toward the wing tip via an outboard vorticity flux. We identify the alula's streamwise position (relative to the leading-edge of the thin wing) as important for vortex steering and the alula's cant angle as important for high-magnitude spanwise flow generation. These findings advance our understanding of the likely ways bird's leverage LEVs to augment slow flight.


title: Sensitivity analysis methods in the biomedical sciences

abstract: Sensitivity analysis is an important part of a mathematical modeller's toolbox for model analysis. In this review paper, we describe the most frequently used sensitivity techniques, discussing their advantages and limitations, before applying each method to a simple model. Also included is a summary of current software packages, as well as a modeller's guide for carrying out sensitivity analyses. Finally, we apply the popular Morris and Sobol methods to two models with biomedical applications, with the intention of providing a deeper understanding behind both the principles of these methods and the presentation of their results.


title: Quadratic soliton mode-locked degenerate optical parametric oscillator

abstract: By identifying the similarities between the coupled-wave equations and the parametrically driven nonlinear Schr\"odinger equation, we unveil the existence condition of the quadratic soliton mode-locked degenerate optical parametric oscillator in the previously unexplored parameter space of near-zero group velocity mismatch. We study the nature of the quadratic solitons and divide their dynamics into two distinctive branches depending on the system parameters. We find the nonlinear interaction between the resonant pump and signal results in phenomena that resemble the dispersive two-photon absorption and the dispersive Kerr effect. Origin of the quadratic soliton perturbation is identified and strategy to mitigate its detrimental effect is developed. Terahertz comb bandwidth and femtosecond pulse duration are attainable in an example periodically poled lithium niobate waveguide resonator in the short-wave infrared and an example orientation-patterned gallium arsenide free-space cavity in the long-wave infrared. The quadratic soliton mode-locking principle can be extended to other material platforms, making it a competitive ultrashort pulse and broadband comb source architecture at the mid-infrared.


title: Pricing Exchange Options under Stochastic Correlation

abstract: In this paper we study the pricing of exchange options when underlying assets have stochastic volatility and stochastic correlation. An approximation using a closed-form approximation based on a Taylor expansion of the conditional price is proposed. Numerical results are illustrated for exchanges between WTI and Brent type oil prices.


title: Distributed Fixed Point Method for Solving Systems of Linear Algebraic Equations

abstract: We present a class of iterative fully distributed fixed point methods to solve a system of linear equations, such that each agent in the network holds one of the equations of the system. Under a generic directed, strongly connected network, we prove a convergence result analogous to the one for fixed point methods in the classical, centralized, framework: the proposed method converges to the solution of the system of linear equations at a linear rate. We further explicitly quantify the rate in terms of the linear system and the network parameters. Next, we show that the algorithm provably works under time-varying directed networks provided that the underlying graph is connected over bounded iteration intervals, and we establish a linear convergence rate for this setting as well. A set of numerical results is presented, demonstrating practical benefits of the method over existing alternatives.


title: Stability of the standing waves of the concentrated NLSE in dimension two

abstract: In this paper we will continue the analysis of two dimensional Schr\"odinger equation with a fixed, pointwise, nonlinearity started in [2, 13]. In this model, the occurrence of a blow-up phenomenon has two peculiar features: the energy threshold under which all solutions blow up is strictly negative and coincides with the infimum of the energy of the standing waves; there is no critical power nonlinearity, i.e., for every power there exist blow-up solutions. Here we study the stability properties of stationary states to verify whether the anomalies mentioned before have any counterpart on the stability features.


title: How many ways a cell can move: the modes of self-propulsion of an active drop

abstract: Numerous physical models have been proposed to explain how cell motility emerges from internal activity, mostly focused on how crawling motion arises from internal processes. Here we offer a classification of self-propulsion mechanisms based on general physical principles, showing that crawling is not the only way for cells to move on a substrate. We consider a thin drop of active matter on a planar substrate and fully characterize its autonomous motion for all three possible sources of driving: (i) the stresses induced in the bulk by active components, which allow in particular tractionless motion, (ii) the self-propulsion of active components at the substrate, which gives rise to crawling motion, and (iii) a net capillary force, possibly self-generated, and coupled to internal activity. We determine travelling-wave solutions to the lubrication equations as a function of a dimensionless activity parameter for each mode of motion. Numerical simulations are used to characterize the drop motion over a wide range of activity magnitudes, and explicit analytical solutions in excellent agreement with the simulations are derived in the weak-activity regime.


title: Some applications of MV algebras

abstract: In this paper, some properties and applications of MV-algebras are provided. We define a Fibonacci sequence in an MV-algebra and we prove that such a stationary sequence gives us an idempotent element. Taking into account of the representation of a finite MV-agebra, by using Boolean elements of this algebra, we prove that such a sequence is always stationary and is not periodic as in the situation when such a sequence is studied on the group (Zn; +), the group of integers modulo n. Moreover, as an application in Coding Theory, to a Boolean algebra is attached a binary block code and is proved that, under some conditions, the converse is also true.


title: Multimode Single-Pass Spatio-temporal Squeezing

abstract: We present a single-pass source of broadband multimode squeezed light with potential application in quantum information and quantum metrology. The source is based on a type I parametric down-conversion (PDC) process inside a bulk nonlinear crystal in a non-collinear configuration. The generated squeezed light exhibits a spatiotemporal multimode behavior that is probed using a homodyne measurement with a local oscillator shaped both spatially and temporally. Finally we follow a covariance matrix based approach to reveal the distribution of the squeezing among several independent temporal and spatial modes. This unambiguously validates the multimode feature of our source.


title: Local Existence of Contact Discontinuities in Relativistic Magnetohydrodynamics

abstract: We study the free boundary problem for a contact discontinuity for the system of relativistic magnetohydrodynamics. A surface of contact discontinuity is a characteristic of this system with no flow across the discontinuity for which the pressure, the velocity and the magnetic field are continuous whereas the density, the entropy and the temperature may have a jump. For the two-dimensional case, we prove the local-in-time existence in Sobolev spaces of a unique solution of the free boundary problem provided that the Rayleigh--Taylor sign condition on the jump of the normal derivative of the pressure is satisfied at each point of the initial discontinuity.


title: High order semi-implicit multistep methods for time dependent partial differential equations

abstract: We consider the construction of semi-implicit linear multistep methods which can be applied to time dependent PDEs where the separation of scales in additive form, typically used in implicit-explicit (IMEX) methods, is not possible. As shown in Boscarino, Filbet and Russo (2016) for Runge-Kutta methods, these semi-implicit techniques give a great flexibility, and allows, in many cases, the construction of simple linearly implicit schemes with no need of iterative solvers. In this work we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype linear advection-diffusion equation and in the setting of strong stability preserving (SSP) methods. Our findings are demonstrated on several examples, including nonlinear reaction-diffusion and convection-diffusion problems.


title: Ultrafast two-photon emission in a doped semiconductor thin film

abstract: As a high-order quantum transition, two-photon emission has an extremely low occurrence rate compared to one-photon emission, thus having been considered a forbidden process. Here, we propose a scheme that allows ultrafast two-photon emission, leveraging highly confined surface plasmon polariton modes in a degenerately-doped, light-emitting semiconductor thin film. The surface plasmon polariton modes are tailored to have simultaneous spectral and spatial overlap with the two-photon emission in the semiconductor. Using degenerately-doped InSb as the prototype material, we show that the two-photon emission can be accelerated by 10 orders of magnitude: from tens of milliseconds to picoseconds, surpassing the one-photon emission rate. Our result provides a semiconductor platform for ultrafast single and entangled photon generation, with a tunable emission wavelength in the mid-infrared.


title: Linear programming bounds for quantum amplitude damping codes

abstract: Given that approximate quantum error-correcting (AQEC) codes have a potentially better performance than perfect quantum error correction codes, it is pertinent to quantify their performance. While quantum weight enumerators establish some of the best upper bounds on the minimum distance of quantum error-correcting codes, these bounds do not directly apply to AQEC codes. Herein, we introduce quantum weight enumerators for amplitude damping (AD) errors and work within the framework of approximate quantum error correction. In particular, we introduce an auxiliary exact weight enumerator that is intrinsic to a code space and moreover, we establish a linear relationship between the quantum weight enumerators for AD errors and this auxiliary exact weight enumerator. This allows us to establish a linear program that is infeasible only when AQEC AD codes with corresponding parameters do not exist. To illustrate our linear program, we numerically rule out the existence of three-qubit AD codes that are capable of correcting an arbitrary AD error.


title: Mobility-assisted Over-the-Air Computation for Backscatter Sensor Networks

abstract: Future intelligent systems will consist of a massive number of battery-less sensors, where quick and accurate aggregation of sensor data will be of paramount importance. Over-the-air computation (AirComp) is a promising technology wherein sensors concurrently transmit their measurements over the wireless channel, and a reader receives the noisy version of a function of measurements due to the superposition property. A key challenge in AirComp is the accurate power alignment of individual transmissions, addressed previously by using conventional precoding methods. In this paper, we investigate a UAVenabled backscatter communication framework, wherein UAV acts both as a power emitter and reader. The mobility of the reader is leveraged to replace the complicated precoding at sensors, where UAV first collects sum channel gains in the first flyover, and then, use these to estimate the actual aggregated sensor data in the second flyover. Our results demonstrate improvements of up to 10 dB in MSE compared to that of a benchmark case where UAV is incognizant of sum channel gains.


title: Nonequilibrium Layer in the Crust of Neutron Stars and Nonequilibrium Beta-Processes in Astrophysics

abstract: The formation of the chemical composition of neutron star envelopes, at high densities is considered. As hot matter is compressed in the process of collapse, which leads to the explosion of a core collapse supernova, the stage of nuclear equilibrium with free neutrino escape, kinetic equilibrium in \b{eta}-processes, and, as a result, the establishment of limited nuclear equilibrium with a fixed number of nuclei takes place. Cold matter is compressed at a fixed number of nuclei whose atomic weight initially does not change and subsequently decreases. A pycnonuclear reaction of the fusion of available nuclei and a decrease in their number begin at the end. The compression of cold matter is accompanied by an increase in the mass fraction of free neutrons. In this case, the chemical composition of the envelope differs significantly from the equilibrium one and contains a considerable store of nuclear energy. Nonequilibrium \b{eta}-reactions proceed at densities exceeding the upper bound for the non-equilibrium layer density, which lead to heating, nuclear energy release, and the possible attainment of a state of complete thermodynamic equilibrium. The thermodynamics of nonequilibrium \b{eta}-processes, which lead to the heating of matter as neutrinos escape freely, is considered.


title: A broadband X-ray view of the NLSy1 1E 0754.6+392.8

abstract: The soft X-ray band of many active galactic nuclei (AGNs) is affected by obscuration due to partially ionised matter crossing our line of sight. In this context, two past XMM-Newton observations (6 months apart) and a simultaneous NuSTAR-Swift ($\sim$8 years later) exposure of the Narrow Line Seyfert 1 galaxy 1E 0754.6+392.8 revealed an intense and variable WA and hints of additional absorbers in the Fe K$\alpha$ band. We present the first X-ray characterisation of this AGN discussing its broadband (0.3-79 keV) spectrum and temporal properties. We conduct a temporal and spectroscopic analysis on two $\sim$10 ks (net exposure) XMM-Newton snapshots performed in April and October 2006. We also study the high energy behaviour of 1E 0754.6+392.8 modelling its broadband spectrum using simultaneous Swift-NuSTAR data. Both phenomenological and physically motivated models are tested. We find the presence of flux variability ($\sim$150% and 30% for 0.3-2 and 2-10 keV bands, respectively) and spectral changes at months timescales ($\Delta\Gamma\sim$0.4). A reflection component that is consistent with being constant over years and arising from relatively cold material far from the central super massive black hole is detected. The main spectral feature shaping the 1E 0754.6+392.8 spectrum is a warm absorber. Such a component is persistent over the years and variability of its ionisation and column density is observed down on months in the ranges 3$\times10^{22} \rm cm^{-2}\lesssim$ N$_{\rm{H}}\lesssim7.2\times10^{22} \rm cm^{-2}$ and 1.5 $\lesssim\log(\xi/{\rm erg~s^{-1}~cm})\lesssim$2.1. Despite the short exposures, we find possible evidence of two additional highly ionised and high-velocity outflow components in absorption. Longer exposures are mandatory in order to characterise the complex outflow in this AGN.


title: AdS Robin solitons and their stability

abstract: We consider the four-dimensional Einstein-Klein-Gordon-AdS system with conformal mass subject to the Robin boundary conditions at infinity. Above a critical value of the Robin parameter, at which the AdS spacetime goes linearly unstable, we prove existence of a family of globally regular static solutions (that we call AdS Robin solitons) and discuss their properties.


title: Head and Tail Localization of C. elegans

abstract: C. elegans is commonly used in neuroscience for behaviour analysis because of it's compact nervous system with well-described connectivity. Localizing the animal and distinguishing between its head and tail are important tasks to track the worm during behavioural assays and to perform quantitative analyses. We demonstrate a neural network based approach to localize both the head and the tail of the worm in an image. To make empirical results in the paper reproducible and promote open source machine learning based solutions for C. elegans behavioural analysis, we also make our code publicly available.


title: Food Supply Chain and Business Model Innovation

abstract: This paper investigates the contribution of business model innovations in improvement of food supply chains. Through a systematic literature review, the notable business model innovations in the food industry are identified, surveyed, and evaluated. Findings reveal that the innovations in value proposition, value creation processes, and value delivery processes of business models are the successful strategies proposed in food industry. It is further disclosed that rural female entrepreneurs, social movements, and also urban conditions are the most important driving forces inducing the farmers to reconsider their business models. In addition, the new technologies and environmental factors are the secondary contributors in business model innovation for the food processors. It is concluded that digitalization has disruptively changed the food distributors models. E-commerce models and internet of things are reported as the essential factors imposing the retailers to innovate their business models. Furthermore, the consumption demand and the product quality are two main factors affecting the business models of all the firms operating in the food supply chain regardless of their positions in the chain. The findings of the current study provide an insight into the food industry to design a sustainable business model to bridge the gap between food supply and food demand.


title: One shot approach to lossy source coding under average distortion constraints

abstract: This paper presents a one shot analysis of the lossy compression problem under average distortion constraints. We calculate the exact expected distortion of a random code. The result is given as an integral formula using a newly defined functional $\tilde{D}(z,Q_Y)$ where $Q_Y$ is the random coding distribution and $z\in [0,1]$. When we plug in the code distribution as $Q_Y$, this functional produces the average distortion of the code, thus provide a converse result utilizing the same functional. Two alternative formulas are provided for $\tilde{D}(z,Q_Y)$, the first involves a supremum over some auxiliary distribution $Q_X$ which has resemblance to the channel coding meta-converse and the other involves an infimum over channels which resemble the well known Shannon distortion-rate function.


title: Estimating the Competitive Storage Model with Stochastic Trends in Commodity Prices

abstract: We propose a state-space model (SSM) for commodity prices that combines the competitive storage model with a stochastic trend. This approach fits into the economic rationality of storage decisions, and adds to previous deterministic trend specifications of the storage model. Parameters are estimated using a particle Markov chain Monte Carlo procedure. Empirical application to four commodity markets shows that the stochastic trend SSM is favored over deterministic trend specifications. The stochastic trend SSM identifies structural parameters that differ from those for deterministic trend specifications. In particular, the estimated price elasticities of demand are significantly larger under the stochastic trend SSM.


title: Unbiased and Efficient Log-Likelihood Estimation with Inverse Binomial Sampling

abstract: The fate of scientific hypotheses often relies on the ability of a computational model to explain the data, quantified in modern statistical approaches by the likelihood function. The log-likelihood is the key element for parameter estimation and model evaluation. However, the log-likelihood of complex models in fields such as computational biology and neuroscience is often intractable to compute analytically or numerically. In those cases, researchers can often only estimate the log-likelihood by comparing observed data with synthetic observations generated by model simulations. Standard techniques to approximate the likelihood via simulation either use summary statistics of the data or are at risk of producing severe biases in the estimate. Here, we explore another method, inverse binomial sampling (IBS), which can estimate the log-likelihood of an entire data set efficiently and without bias. For each observation, IBS draws samples from the simulator model until one matches the observation. The log-likelihood estimate is then a function of the number of samples drawn. The variance of this estimator is uniformly bounded, achieves the minimum variance for an unbiased estimator, and we can compute calibrated estimates of the variance. We provide theoretical arguments in favor of IBS and an empirical assessment of the method for maximum-likelihood estimation with simulation-based models. As case studies, we take three model-fitting problems of increasing complexity from computational and cognitive neuroscience. In all problems, IBS generally produces lower error in the estimated parameters and maximum log-likelihood values than alternative sampling methods with the same average number of samples. Our results demonstrate the potential of IBS as a practical, robust, and easy to implement method for log-likelihood evaluation when exact techniques are not available.

categories: cs.LG q-bio.NC q-bio.QM stat.CO stat.ME stat.ML


title: Non-perturbative signatures of non-linear Compton scattering

abstract: The probabilities of various elementary laser - photon - electron/positron interactions display in selected phase space and parameter regions typical non-perturbative dependencies such as $\propto {\cal P} \exp\{- a E_{crit} /E\}$, where ${\cal P}$ is a pre-exponential factor, $E_{crit}$ denotes the critical Sauter-Schwinger field strength, and $E$ characterizes the (laser) field strength. While the Schwinger process with $a = a_S \equiv \pi$ and the non-linear Breit-Wheeler process in the tunneling regime with $a = a_{n \ell BW} \equiv 4 m / 3 \omega'$ (with $\omega'$ the probe photon energy and $m$ the electron/positron mass) are famous results, we point out here that also the non-linear Compton scattering exhibits a similar behavior when focusing on high harmonics. Using a suitable cut-off $c > 0$, the factor $a$ becomes $a = a_{n \ell C} \equiv \frac23 c m /(p_0 + \sqrt{p_0^2 -m^2)}$. This opens the avenue towards a new signature of the boiling point of the vacuum even for field strengths $E$ below $E_{crit}$ by employing a high electron beam-energy $p_0$ to counter balance the large ratio $E_{crit} / E$ by a small factor $a$ to achieve $E / a \to E_{crit}$. In the weak-field regime, the cut-off facilitates a threshold leading to multi-photon signatures showing up in the total cross section at sub-threshold energies.


title: Calibration of ground based survey data using Gaia: Application to DES

abstract: The calibration of ground based optical imaging data to photometric accuracy of 10 mmag over the full survey area and to color uniformity better than 5 mmag on the scale of the VIS focal plane is a key science requirement for the Euclid mission. These accuracies enable stable photometric redshifts of galaxies and modeling of the color dependent VIS PSF for weak lensing studies. We use the Gaia photometry to calibrate the $g/r/i/z$ magnitudes of Dark Energy Survey (DES) data to meet the stringent Euclid requirements. The Gaia G band magnitude along with the BP-RP color information of stars observed in the DES single epoch (SE) exposures are used to derive the transformation from Gaia to DES photometry for individual CCDs and to characterize persistent photometric errors across the DECam focal plane. We use the color dependence of these persistent errors to characterize the $g/r/i/z$ bandpass variations across the DECam focal plane.


title: Domain Adaptive Bootstrap Aggregating

abstract: When there is a distributional shift between data used to train a predictive algorithm and current data, performance can suffer. This is known as the domain adaptation problem. Bootstrap aggregating, or bagging, is a popular method for improving stability of predictive algorithms, while reducing variance and protecting against over-fitting. This article proposes a domain adaptive bagging method coupled with a new iterative nearest neighbor sampler. The key idea is to draw bootstrap samples from the training data in such a manner that their distribution equals that of new testing data. The proposed approach provides a general ensemble framework that can be applied to arbitrary classifiers. We further modify the method to allow anomalous samples in the test data corresponding to outliers in the training data. Theoretical support is provided, and the approach is compared to alternatives in simulations and real data applications.


title: Machine-Learning Studies on Spin Models

abstract: With the recent developments in machine learning, Carrasquilla and Melko have proposed a paradigm that is complementary to the conventional approach for the study of spin models. As an alternative to investigating the thermal average of macroscopic physical quantities, they have used the spin configurations for the classification of the disordered and ordered phases of a phase transition through machine learning. We extend and generalize this method. We focus on the configuration of the long-range correlation function instead of the spin configuration itself, which enables us to provide the same treatment to multi-component systems and the systems with a vector order parameter. We analyze the Berezinskii-Kosterlitz-Thouless (BKT) transition with the same technique to classify three phases: the disordered, the BKT, and the ordered phases. We also present the classification of a model using the training data of a different model.


title: Thouless Energy Challenges Thermalization on the Ergodic Side of the Many-Body Localization Transition

abstract: We study the ergodic side of the many-body localization transition in its standard model, the disordered Heisenberg quantum spin chain. We show that the Thouless energy, extracted from long-range spectral statistics and the power-spectrum of the full momentum distribution fluctuations, is not large enough to guarantee thermalization. We find that both estimates coincide and behave non-monotonically, exhibiting a strong peak at an intermediate value of the disorder. Furthermore, we show that non-thermalizing initial conditions occur well within the ergodic phase with larger probability than expected. Finally, we propose a mechanism, driven by the Thouless energy and the presence of anomalous events, for the transition to the localized phase.


title: Comparison of time reversal symmetric and asymmetric nano-swimmers oriented with an electric field in soft matter

abstract: Using molecular dynamics simulations we compare the motion of a nano-swimmer based on Purcell suggested motor with a time asymmetrical cycle with the motion of the same molecular motor with a time symmetrical cycle. We show that Purcell theorem still holds at the nanoscale, despite the local structure and the medium fluctuations. Then, with the purpose of both orienting the swimmer displacement and increasing the breakdown of the theorem, we study the effect of an electric field on a polarized version of these swimmers. For small and large fields, the time asymmetrical swimmer is more efficient, as suggested by Purcell. However we find a field range for which Purcell theorem is broken for the time symmetric motor. Results suggest that the breakdown of the theorem is arising from the competition of the orientation field and Brownian forces, while for larger fields the field destroys the effect of fluctuations restoring the theorem.


title: Bag of Tricks for Retail Product Image Classification

abstract: Retail Product Image Classification is an important Computer Vision and Machine Learning problem for building real world systems like self-checkout stores and automated retail execution evaluation. In this work, we present various tricks to increase accuracy of Deep Learning models on different types of retail product image classification datasets. These tricks enable us to increase the accuracy of fine tuned convnets for retail product image classification by a large margin. As the most prominent trick, we introduce a new neural network layer called Local-Concepts-Accumulation (LCA) layer which gives consistent gains across multiple datasets. Two other tricks we find to increase accuracy on retail product identification are using an instagram-pretrained Convnet and using Maximum Entropy as an auxiliary loss for classification.


title: Derivation of the Landau-Pekar equations in a many-body mean-field limit

abstract: We consider the Fr\"ohlich Hamiltonian in a mean-field limit where many bosonic particles weakly couple to the quantized phonon field. For large particle number and suitably small coupling, we show that the dynamics of the system is approximately described by the Landau-Pekar equations. These describe a Bose-Einstein condensate interacting with a classical polarization field, whose dynamics is effected by the condensate, i.e., the back-reaction of the phonons that are created by the particles during the time evolution is of leading order.


title: Fast is better than free: Revisiting adversarial training

abstract: Adversarial training, a method for learning robust deep networks, is typically assumed to be more expensive than traditional training due to the necessity of constructing adversarial examples via a first-order method like projected gradient decent (PGD). In this paper, we make the surprising discovery that it is possible to train empirically robust models using a much weaker and cheaper adversary, an approach that was previously believed to be ineffective, rendering the method no more costly than standard training in practice. Specifically, we show that adversarial training with the fast gradient sign method (FGSM), when combined with random initialization, is as effective as PGD-based training but has significantly lower cost. Furthermore we show that FGSM adversarial training can be further accelerated by using standard techniques for efficient training of deep networks, allowing us to learn a robust CIFAR10 classifier with 45% robust accuracy to PGD attacks with $\epsilon=8/255$ in 6 minutes, and a robust ImageNet classifier with 43% robust accuracy at $\epsilon=2/255$ in 12 hours, in comparison to past work based on "free" adversarial training which took 10 and 50 hours to reach the same respective thresholds. Finally, we identify a failure mode referred to as "catastrophic overfitting" which may have caused previous attempts to use FGSM adversarial training to fail. All code for reproducing the experiments in this paper as well as pretrained model weights are at https://github.com/locuslab/fast_adversarial.


title: CMB Cold Spot in the Planck light

abstract: The Cold Spot, with an unusually cold region surrounded by a hot ring, is a statistically significant anomaly in the Cosmic Microwave Background (CMB) sky. In this work we assess whether different sets of multiple subvoids based on the 2dF-VST ATLAS Cold Spot galaxy redshift survey or a collapsing cosmic texture could have produced such an anomaly through a simultaneous search for their gravitational redshift and lensing signatures on the {\it Planck} CMB temperature anisotropies. We use patches with radii $R=10^\circ$ and $R=20^\circ$ to account for the inner cold region as well as the outer hot ring. As the void model, we explore two sets of $\Lambda$LTB templates characterized by different values of the model's free parameters, and a top-hat void template. We detect higher-than-expected gravitational redshift amplitudes for the first two sets, $\mathcal{A}_{\rm rs}=5.4\pm 1.4$ and $\mathcal{A}_{\rm rs}=14.4\pm 3.8$, and lower than expected for the top-hat model, $\mathcal{A}_{\rm rs}=0.3 \pm 0.1$. The amplitudes for the lensing imprint are consistent with zero for all these subvoid models. The estimated amplitude for the texture imprint from the gravitational redshift measurement implies the energy scale of the texture, parametrized by $\epsilon$, to be $\epsilon= (7.6\pm2.0)\times 10^{-5}$, with no detection of the lensing trace. We note that the deviation of the subvoid amplitudes from unity and the inability of the texture and some of the void profiles to reproduce the hot ring indicate theoretical insufficiencies, either in the construction of the model or in the assumed gravitational and cosmological framework leading to the imprints for the structures.


title: Specific properties of supercooled water in light of water anomalies

abstract: We review the effect of water anomalies on the properties of low temperature water. When supercooled, liquids dynamical properties change drastically. Supercooled liquids undergo an at least exponential decrease of their diffusion coefficient when temperature decreases while their structure merely does not change. We discuss how that still unexplained change of dynamical properties at low temperatures affect water differently from other liquids and what can be deduced from it.


title: Imaging entanglement correlations with a single-photon avalanche diode camera

abstract: Spatial correlations between two photons are the key resource in realising many quantum imaging schemes. Measurement of the bi-photon correlation map is typically performed using single-point scanning detectors or single-photon cameras based on CCD technology. However, both approaches are limited in speed due to the slow scanning and the low frame-rate of CCD-based cameras, resulting in data acquisition times on the order of many hours. Here we employ a high frame rate, single photon avalanche diode (SPAD) camera, to measure the spatial joint probability distribution of a bi-photon state produced by spontaneous parametric down-conversion, with statistics taken over $10^7$ frames acquired in just 140 seconds. We verified the presence of spatial entanglement between our photon pairs through the violation of an Einstein-Podolsky-Rosen criterion, with a confidence level of 227 sigmas. Our work demonstrates the potential of SPAD cameras in the rapid characterisation of photon correlations, leading the way towards quantum imaging in real-time.


title: Towards causality-aware predictions in static anticausal machine learning tasks: the linear structural causal model case

abstract: We propose a counterfactual approach to train ``causality-aware" predictive models that are able to leverage causal information in static anticausal machine learning tasks (i.e., prediction tasks where the outcome influences the features). In applications plagued by confounding, the approach can be used to generate predictions that are free from the influence of observed confounders. In applications involving observed mediators, the approach can be used to generate predictions that only capture the direct or the indirect causal influences. Mechanistically, we train supervised learners on (counterfactually) simulated features which retain only the associations generated by the causal relations of interest. We focus on linear models, where analytical results connecting covariances, causal effects, and prediction mean squared errors are readily available. Quite importantly, we show that our approach does not require knowledge of the full causal graph. It suffices to know which variables represent potential confounders and/or mediators. We discuss the stability of the method with respect to dataset shifts generated by selection biases and validate the approach using synthetic data experiments.


title: Spatial correlation of the thermally generated electromagnetic field in layered media

abstract: A general formulation for the cross-spectral density tensor enabling calculation of the spatial correlation of the thermally generated electromagnetic field in layered media is derived. The formulation is based on fluctuational electrodynamics, and is thus applicable in the near and far field of heat sources. The resulting cross-spectral density tensor is written in terms of a single integration over the parallel wavevector, as the angular integrations leading to numerical instability are evaluated analytically. Using this formulation, the spatial correlation length in the near field of a film made of silicon carbide (SiC) supporting surface phonon-polaritons (SPhPs) in the infrared is analyzed. It is shown that the spatial correlation length of a SiC heat source suspended in vacuum decreases substantially by decreasing its thickness owing to SPhP coupling. In the limit of a 10-nm-thick SiC film, the spatial correlation length is similar to that of a blackbody. The results also reveal that it is possible to control the spatial coherence of a thin SiC heat source via dielectric and metallic substrates, respectively allowing and preventing SPhP coupling. This suggests that active modulation of thermal emission via thin films supporting surface polaritons in the infrared is possible by using a phase change material substrate such as vanadium dioxide.


title: Evaluation of neutrinoless double beta decay: QCD running to sub-GeV scales

abstract: We evaluate QCD effects in the neutrinoless double beta ($0\nu\beta\beta$) decay, originating from new physics short-range mechanism in the form of five dimension-9 operators. For this, we employ the one-loop and two-loop renormalization group equations (RGEs) for the corresponding Wilson coefficients, performing the RGE-evolution from the new physics scales (estimated as $\Lambda ~ \sim 10^2$ GeV) to the typical spacelike $0\nu\beta\beta$-scale $Q \sim 0.1$ GeV. Since the latter scale is clearly nonperturbative, we apply various infrared-safe (IR-safe) variants of QCD where the running coupling has no Landau singularities at low spacelike $Q$. We point out that the correct treatment of the IR-safe analogs of the (noninteger) powers of the couplings is important. It turns out that in most cases of the considered operators the resulting QCD effects can be significant in this process, i.e., can be stronger than the effects of the present uncertainties in the nuclear matrix elements.


title: A comprehensive deep learning-based approach to reduced order modeling of nonlinear time-dependent parametrized PDEs

abstract: Traditional reduced order modeling techniques such as the reduced basis (RB) method (relying, e.g., on proper orthogonal decomposition (POD)) suffer from severe limitations when dealing with nonlinear time-dependent parametrized PDEs, because of the fundamental assumption of linear superimposition of modes they are based on. For this reason, in the case of problems featuring coherent structures that propagate over time such as transport, wave, or convection-dominated phenomena, the RB method usually yields inefficient reduced order models (ROMs) if one aims at obtaining reduced order approximations sufficiently accurate compared to the high-fidelity, full order model (FOM) solution. To overcome these limitations, in this work, we propose a new nonlinear approach to set reduced order models by exploiting deep learning (DL) algorithms. In the resulting nonlinear ROM, which we refer to as DL-ROM, both the nonlinear trial manifold (corresponding to the set of basis functions in a linear ROM) as well as the nonlinear reduced dynamics (corresponding to the projection stage in a linear ROM) are learned in a non-intrusive way by relying on DL algorithms; the latter are trained on a set of FOM solutions obtained for different parameter values. In this paper, we show how to construct a DL-ROM for both linear and nonlinear time-dependent parametrized PDEs; moreover, we assess its accuracy on test cases featuring different parametrized PDE problems. Numerical results indicate that DL-ROMs whose dimension is equal to the intrinsic dimensionality of the PDE solutions manifold are able to approximate the solution of parametrized PDEs in situations where a huge number of POD modes would be necessary to achieve the same degree of accuracy.


title: On the challenges ahead of spatial scientometrics focusing on the city level

abstract: Since the mid-1970s, it has become highly acknowledged to measure and evaluate changes in international research collaborations and the scientific performance of institutions and countries through the prism of bibliometric and scientometric data. Spatial bibliometrics and scientometrics (henceforward spatial scientometrics) have traditionally focused on examining both country and regional levels; however, in recent years, numerous spatial analyses on the city level have been carried out. While city-level scientometric analyses have gained popularity among policymakers and statistical/economic research organizations, researchers in the field of bibliometrics are divided regarding whether it is possible to observe the spatial unit 'city' through bibliometric and scientometric tools. After systematically scrutinizing relevant studies in the field, three major problems have been identified: 1) there is no standardized method of how cities should be defined and how metropolitan areas should be delineated, 2) there is no standardized method of how bibliometric and scientometric data on the city level should be collected and processed and 3) it is not clearly defined how cities can profit from the results of bibliometric and scientometric analysis focusing on them. This paper investigates major challenges ahead of spatial scientometrics, focusing on the city level and presents some possible solutions.


title: Towards the Small Quasi-Kernel Conjecture

abstract: Let $D=(V,A)$ be a digraph. A vertex set $K\subseteq V$ is a quasi-kernel of $D$ if $K$ is an independent set in $D$ and for every vertex $v\in V\setminus K$, $v$ is at most distance 2 from $K$. In 1974, Chv\'atal and Lov\'asz proved that every digraph has a quasi-kernel. P. L. Erd\H{o}s and L. A. Sz\'ekely in 1976 conjectured that if every vertex of $D$ has a positive indegree, then $D$ has a quasi-kernel of size at most $|V|/2$. This conjecture is only confirmed for narrow classes of digraphs, such as semicomplete multipartite, quasi-transitive, or locally demicomplete digraphs. In this note, we state a similar conjecture for all digraphs, show that the two conjectures are equivalent, and prove that both conjectures hold for a class of digraphs containing all orientations of 4-colorable graphs (in particular, of all planar graphs).


title: Tilting modules arising from knot invariants

abstract: We construct tilting modules over Jacobian algebras arising from knots. To a two-bridge knot $L[a_1,\ldots,a_n]$, we associate a quiver $Q$ with potential and its Jacobian algebra $A$. We construct a family of canonical indecomposable $A$-modules $M(i)$, each supported on a different specific subquiver $Q(i)$ of $Q$. Each of the $M(i)$ is expected to parametrize the Jones polynomial of the knot. We study the direct sum $M=\oplus_iM(i)$ of these indecomposables inside the module category of $A$ as well as in the cluster category. In this paper we consider the special case where the two-bridge knot is given by two parameters $a_1,a_2$. We show that the module $M$ is rigid and $\tau$-rigid, and we construct a completion of $M$ to a tilting (and $\tau$-tilting) $A$-module $T$. We show that the endomorphism algebra $\operatorname{End}_AT$ of $T$ is isomorphic to $A$, and that the mapping $T\mapsto A[1]$ induces a cluster automorphism of the cluster algebra $\mathcal{A}(Q)$. This automorphism is of order two. Moreover, we give a mutation sequence that realizes the cluster automorphism. In particular, we show that the quiver $Q$ is mutation equivalent to an acyclic quiver of type $T_{p,q,r}$ (a tree with three branches). This quiver is of finite type if $(a_1,a_2)=(a_1,2), (1,a_2),$ or $(2,3)$, it is tame for $(a_1,a_2)=(2,4)$ or $(3,3)$, and wild otherwise.


title: Heterogeneous Computation Assignments in Coded Elastic Computing

abstract: We study the optimal design of a heterogeneous coded elastic computing (CEC) network where machines have varying relative computation speeds. CEC introduced by Yang {\it et al.} is a framework which mitigates the impact of elastic events, where machines join and leave the network. A set of data is distributed among storage constrained machines using a Maximum Distance Separable (MDS) code such that any subset of machines of a specific size can perform the desired computations. This design eliminates the need to re-distribute the data after each elastic event. In this work, we develop a process for an arbitrary heterogeneous computing network to minimize the overall computation time by defining an optimal computation load, or number of computations assigned to each machine. We then present an algorithm to define a specific computation assignment among the machines that makes use of the MDS code and meets the optimal computation load.


title: An adjunction inequality obstruction to isotopy of embedded surfaces in 4-manifolds

abstract: Consider a smooth $4$-manifold $X$ and a diffeomorphism $f : X \to X$. We give an obstruction in the form of an adjunction inequality for an embedded surface in $X$ to be isotopic to its image under $f$. It follows that the minimal genus of a surface representing a given homology class and which is isotopic to its image under $f$ is generally larger than the minimal genus without the isotopy condition. We give examples where the inequality is strict. We use our obstruction to construct examples of infinitely many embedded surfaces which are all continuously isotopic but mutually non-isotopic smoothly.


title: Beam Tracking with Photon-Counting Detector Arrays in Free-Space Optical Communications

abstract: Optical beam center position on an array of detectors is an important (hidden) parameter that is essential not only from a tracking perspective, but is also important for optimal detection of Pulse Position Modulation symbols in free-space optical communications. In this paper, we have examined the beam position estimation problem for photon-counting detector arrays, and to this end, we have proposed and analyzed a number of non-Bayesian beam position estimators. These estimators are compared in terms of their mean-square error, bias and the probability of error performance. Additionally, the Cramer-Rao Lower Bounds (CRLB) of the tracking error is also derived, and the CRLB curves give us additional insights concerning the effect of number of detectors and the beam radius on mean-square error performance. Finally, the effect of beam position estimation on the probability of error performance is investigated, and our study concludes that the probability of error of the system is minimized when the beam position on the array is estimated as accurately as possible.


title: On stochastic equations with drift in $L_{d}$

abstract: For It\^o stochastic equations in $\mathbb{R}^{d}$ with drift in $L_{d}$ several results are discussed such as the existence of weak solutions, the existence of the corresponding Markov process, Aleksandrov type estimates of their Green's functions, which yield their summability to the power of $d/(d-1)$, the Fabes-Stroock type estimates which show that Green's functions are summable to a higher degree, the Fanghua Lin type estimates, which are one of the main tools in the $W^{2}_{p}$-theory of fully nonlinear elliptic equations, the fact that Green's functions are in the class $A_{\infty}$ of Muckenhoupt and a few other results.


title: Finite-Level Quantization Procedures for Construction and Decoding of Polar Codes

abstract: We consider finite-level, symmetric quantization procedures for construction and decoding of polar codes. Whether polarization occurs in the presence of quantization is not known in general. Hassani and Urbanke have shown that a simple three-level quantization procedure polarizes and they have proposed a calculation method to obtain a lower bound for achievable rates. We find an improved calculation method for achievable rates and also the exact asymptotic behavior of the block error probability under the aforementioned simple case. We then prove that certain D-level quantization schemes polarize and give a lower bound on achievable rates. Furthermore, we show that a broad class of quantization procedures result in a weaker form of the polarization phenomenon.


title: Adaptive Acquisition Schemes for Photon-Limited Free-Space Optical Communications

abstract: Acquisition and tracking systems form an important component of free-space optical communications due to directional nature of the optical signal. Acquisition subsystems are needed in order to search and locate the receiver terminal in an uncertainty/search region with very narrow laser beams. In this paper, we have proposed and analyzed two adaptive search schemes for acquisition systems that perform better---for the low probability of detection---than the spiral scanning approach. The first of these schemes, the adaptive spiral search, provides a better acquisition time performance by dividing the search region into a number of smaller subregions, and prioritizing search in regions of higher probability mass. The second technique---the shotgun approach---searches the region in a random manner by sampling the search region according to a Gaussian distribution. The adaptive spiral scheme outperforms the shotgun approach in terms of acquisition time, especially if the number of search subregions is large enough. However, a higher pointing accuracy is required by the adaptive spiral search in order to search the region precisely. On the other hand, the shotgun scanning approach does not require such stringent pointing accuracy.


title: Membership Inference Attacks Against Object Detection Models

abstract: Machine learning models can leak information regarding the dataset they have trained. In this paper, we present the first membership inference attack against black-boxed object detection models that determines whether the given data records are used in the training. To attack the object detection model, we devise a novel method named as called a canvas method, in which predicted bounding boxes are drawn on an empty image for the attack model input. Based on the experiments, we successfully reveal the membership status of privately sensitive data trained using one-stage and two-stage detection models. We then propose defense strategies and also conduct a transfer attack between the models and datasets. Our results show that object detection models are also vulnerable to inference attacks like other models.


title: A note on complex-hyperbolic Kleinian groups

abstract: Let $\Gamma$ be a discrete group of isometries acting on the complex hyperbolic $n$-space $\mathbb{H}^n_\mathbb{C}$. In this note, we prove that if $\Gamma$ is convex-cocompact, torsion-free, and the critical exponent $\delta(\Gamma)$ is strictly lesser than $2$, then the complex manifold $\mathbb{H}^n_\mathbb{C}/\Gamma$ is Stein. We also discuss several related conjectures.


title: Notes on Exact Power Calculations for t Tests and Analysis of Covariance

abstract: Tang derived the exact power formulae for t tests and analysis of covariance (ANCOVA) in superiority, noninferiority and equivalence trials. The power calculation in equivalence trials can be simplified by using Owen's Q function, which is available in standard statistical software. We extend the exact power determination method for ANCOVA to unstratified and stratified multi-arm randomized trials. The method is applied to the design of multi-arm trials and gold standard noninferiority trials.


title: Leveraging Quantum Annealing for Large MIMO Processing in Centralized Radio Access Networks

abstract: User demand for increasing amounts of wireless capacity continues to outpace supply, and so to meet this demand, significant progress has been made in new MIMO wireless physical layer techniques. Higher-performance systems now remain impractical largely only because their algorithms are extremely computationally demanding. For optimal performance, an amount of computation that increases at an exponential rate both with the number of users and with the data rate of each user is often required. The base station's computational capacity is thus becoming one of the key limiting factors on wireless capacity. QuAMax is the first large MIMO centralized radio access network design to address this issue by leveraging quantum annealing on the problem. We have implemented QuAMax on the 2,031 qubit D-Wave 2000Q quantum annealer, the state-of-the-art in the field. Our experimental results evaluate that implementation on real and synthetic MIMO channel traces, showing that 10~$\mu$s of compute time on the 2000Q can enable 48 user, 48 AP antenna BPSK communication at 20 dB SNR with a bit error rate of $10^{-6}$ and a 1,500 byte frame error rate of $10^{-4}$.


title: Active particles with fractional rotational Brownian motion

abstract: We study the two-dimensional overdamped motion of an active particle whose orientational dynamics is subject to fractional Brownian noise, whereas its position is affected by self-propulsion and Brownian fluctuations. From a Langevin-like model of active motion with constant swimming speed, we derive the corresponding Fokker-Planck equation, from which we find the angular probability density of the particle orientation for arbitrary values of the Hurst exponent that characterizes the fractional rotational noise. We provide analytical expressions for the velocity autocorrelation function and the translational mean-squared displacement, which show that active diffusion effectively emerges in the long-time limit for all values of the Hurst exponent. The corresponding expressions for the active diffusion coefficient and the effective rotational diffusion time are also derived. Our results are compared with numerical simulations of active particles with rotational motion driven by fractional Brownian noise, with which we find an excellent agreement.


title: TopoAna: A generic tool for the event type analysis of inclusive Monte-Carlo samples in high energy physics experiments

abstract: Inclusive Monte-Carlo samples are indispensable for signal selection and background suppression in many high energy physics experiments. A clear knowledge of the physics processes involved in the samples, including the types of processes and the number of processes in each type, is a great help to investigating signals and backgrounds. To help analysts obtain the physics process information from the truth information of the samples, we develop a physics process analysis program, TopoAna, with C++, ROOT, and LaTeX. The program implements the functionalities of component analysis and signal identification with many kinds of fine, customizable classification and matching algorithms. It tags physics processes in individual events accurately in the output root files, and exports the physics process information at the sample level clearly to the output plain text, tex source, and pdf files. Independent of specific software frameworks, the program is applicable to many experiments. At present, it has come into use in three $e^+e^-$ colliding experiments: the BESIII, Belle, and Belle II experiments. The use of the program in other similar experiments is also prospective.


title: The sharp Adams type inequalities in the hyperbolic spaces under the Lorentz-Sobolev norms

abstract: Let $2\leq m < n$ and $q \in (1,\infty)$, we denote by $W^mL^{\frac nm,q}(\mathbb H^n)$ the Lorentz-Sobolev space of order $m$ in the hyperbolic space $\mathbb H^n$. In this paper, we establish the following Adams inequality in the Lorentz-Sobolev space $W^m L^{\frac nm,q}(\mathbb H^n)$ \[ \sup_{u\in W^mL^{\frac nm,q}(\mathbb H^n),\, \|\nabla_g^m u\|_{\frac nm,q}\leq 1} \int_{\mathbb H^n} \Phi_{\frac nm,q}\big(\beta_{n,m}^{\frac q{q-1}} |u|^{\frac q{q-1}}\big) dV_g < \infty \] for $q \in (1,\infty)$ if $m$ is even, and $q \in (1,n/m)$ if $m$ is odd, where $\beta_{n,m}^{q/(q-1)}$ is the sharp exponent in the Adams inequality under Lorentz-Sobolev norm in the Euclidean space. To our knowledge, much less is known about the Adams inequality under the Lorentz-Sobolev norm in the hyperbolic spaces. We also prove an improved Adams inequality under the Lorentz-Sobolev norm provided that $q\geq 2n/(n-1)$ if $m$ is even and $2n/(n-1) \leq q \leq \frac nm$ if $m$ is odd, \[ \sup_{u\in W^mL^{\frac nm,q}(\mathbb H^n),\, \|\nabla_g^m u\|_{\frac nm,q}^q -\lambda \|u\|_{\frac nm,q}^q \leq 1} \int_{\mathbb H^n} \Phi_{\frac nm,q}\big(\beta_{n,m}^{\frac q{q-1}} |u|^{\frac q{q-1}}\big) dV_g < \infty \] for any $0< \lambda < C(n,m,n/m)^q$ where $C(n,m,n/m)^q$ is the sharp constant in the Lorentz-Poincar\'e inequality. Finally, we establish a Hardy-Adams inequality in the unit ball when $m\geq 3$, $n\geq 2m+1$ and $q \geq 2n/(n-1)$ if $m$ is even and $2n/(n-1) \leq q \leq n/m$ if $m$ is odd \[ \sup_{u\in W^mL^{\frac nm,q}(\mathbb H^n),\, \|\nabla_g^m u\|_{\frac nm,q}^q -C(n,m,\frac nm)^q \|u\|_{\frac nm,q}^q \leq 1} \int_{\mathbb B^n} \exp\big(\beta_{n,m}^{\frac q{q-1}} |u|^{\frac q{q-1}}\big) dx < \infty. \]


title: The sharp Sobolev type inequalities in the Lorentz--Sobolev spaces in the hyperbolic spaces

abstract: Let $W^1L^{p,q}(\mathbb H^n)$, $1\leq q,p < \infty$ denote the Lorentz-Sobolev spaces of order one in the hyperbolic spaces $\mathbb H^n$. Our aim in this paper is three-fold. First of all, we establish a sharp Poincar\'e inequality in $W^1L^{p,q}(\mathbb H^n)$ with $1\leq q \leq p$ which generalizes the result in \cite{NgoNguyenAMV} to the setting of Lorentz-Sobolev spaces. Second, we prove several sharp Poincar\'e-Sobolev type inequalities in $W^1L^{p,q}(\mathbb H^n)$ with $1\leq q \leq p < n$ which generalize the results in \cite{NguyenPS2018} to the setting of Lorentz-Sobolev spaces. Finally, we provide the improved Moser-Trudinger type inequalities in $W^1L^{n,q}(\mathbb{H}^n)$ in the critical case $p= n$ with $1\leq q \leq n$ which generalize the results in \cite{NguyenMT2018} and improve the results in \cite{YangLi2019}. In the proof of the main results, we shall prove a P\'olya--Szeg\"o type principle in $W^1 L^{p,q}(\mathbb H^n)$ with $1\leq q \leq p$ which maybe is of independent interest.


title: Tribimaximal Mixing in the $SU(5) \times \mathcal{T}_{13}$ Texture

abstract: We extend the recently proposed $SU(5) \times \mathcal{T}_{13}$ model for the asymmetric texture to the up-type quark and seesaw sectors. The hierarchical up-type quark masses are generated from higher-dimensional operators involving family-singlet Higgses, gauge-singlet familons, and vectorlike messengers. The complex-tribimaximal (TBM) seesaw mixing arises from the vacuum structure of a minimal number of familons, resulting in an alignment between the Yukawa and Majorana matrices of the seesaw formula. Introducing four right-handed neutrinos, normal ordering of the light neutrino masses is obtained, with $m_{\nu_1} = 27.6\ \mathrm{meV}$, $m_{\nu_2} = 28.9\ \mathrm{meV}$ and $m_{\nu_3} = 57.8\ \mathrm{meV}$. Their sum almost saturates Planck's cosmological upper bound ($120$ $\text{meV}$). The right-handed neutrino masses are expressed in terms of two parameters for a particular choice of familon vacuum alignment. We predict the $\require{cancel}\cancel{CP}$ Jarlskog-Greenberg invariant to be $|\mathcal{J}| = 0.028$, consistent with the current PDG estimate, and Majorana invariants $|\mathcal{I}_1| = 0.106$ and $|\mathcal{I}_2| = 0.011$. A sign ambiguity in the model parameters leads to two possibilities for the invariant mass parameter $|m_{\beta \beta}|$: $13.02$ or $25.21$ $\text{meV}$, both within an order of magnitude of the most rigorous experimental upper limit ($61$--$165$ $\text{meV}$).


title: LinearFold: linear-time approximate RNA folding by 5'-to-3' dynamic programming and beam search

abstract: Motivation: Predicting the secondary structure of an RNA sequence is useful in many applications. Existing algorithms (based on dynamic programming) suffer from a major limitation: their runtimes scale cubically with the RNA length, and this slowness limits their use in genome-wide applications. Results: We present a novel alternative $O(n^3)$-time dynamic programming algorithm for RNA folding that is amenable to heuristics that make it run in $O(n)$ time and $O(n)$ space, while producing a high-quality approximation to the optimal solution. Inspired by incremental parsing for context-free grammars in computational linguistics, our alternative dynamic programming algorithm scans the sequence in a left-to-right (5'-to-3') direction rather than in a bottom-up fashion, which allows us to employ the effective beam pruning heuristic. Our work, though inexact, is the first RNA folding algorithm to achieve linear runtime (and linear space) without imposing constraints on the output structure. Surprisingly, our approximate search results in even higher overall accuracy on a diverse database of sequences with known structures. More interestingly, it leads to significantly more accurate predictions on the longest sequence families in that database (16S and 23S Ribosomal RNAs), as well as improved accuracies for long-range base pairs (500+ nucleotides apart), both of which are well known to be challenging for the current models. Availability: Our source code is available at https://github.com/LinearFold/LinearFold, and our webserver is at http://linearfold.org (sequence limit: 100,000nt).

categories: q-bio.BM cs.DS math.CO physics.bio-ph q-bio.QM


title: Random iterations of maps on $\mathbb{R}^{k}$: asymptotic stability, synchronization and functional central limit theorem

abstract: We study independent and identically distributed random iterations of continuous maps defined on a connected closed subset $S$ of the Euclidean space $\mathbb{R}^{k}$. We assume the maps are monotone (with respect to a suitable partial order) and a "topological" condition on the maps. Then, we prove the existence of a pullback random attractor whose distribution is the unique stationary measure of the random iteration, and we obtain the synchronization of random orbits. As a consequence of the synchronization phenomenon, a functional central limit theorem is established.


title: A New Approach to a Theorem of Eng

abstract: The main aim of this work is to give a case-free algebraic proof for a theorem of Eng on the Poincar\'e polynomial of parabolic quotients of finite Coxeter groups evaluated at -1.


title: Modelling Orebody Structures: Block Merging Algorithms and Block Model Spatial Restructuring Strategies Given Mesh Surfaces of Geological Boundaries

abstract: This paper describes a framework for capturing geological structures in a 3D block model and improving its spatial fidelity given new mesh surfaces. Using surfaces that represent geological boundaries, the objectives are to identify areas where refinement is needed, increase spatial resolution to minimize surface approximation error, reduce redundancy to increase the compactness of the model and identify the geological domain on a block-by-block basis. These objectives are fulfilled by four system components which perform block-surface overlap detection, spatial structure decomposition, sub-blocks consolidation and block tagging, respectively. The main contributions are a coordinate-ascent merging algorithm and a flexible architecture for updating the spatial structure of a block model when given multiple surfaces, which emphasizes the ability to selectively retain or modify previously assigned block labels. The techniques employed include block-surface intersection analysis based on the separable axis theorem and ray-tracing for establishing the location of blocks relative to surfaces. To demonstrate the robustness and applicability of the proposed block merging strategy in a more narrow setting, it is used to reduce block fragmentation in an existing model where surfaces are not given and the minimum block size is fixed. To obtain further insight, a systematic comparison with octree subblocking subsequently illustrates the inherent constraints of dyadic hierarchical decomposition and the importance of inter-scale merging. The results show the proposed method produces merged blocks with less extreme aspect ratios and is highly amenable to parallel processing. The overall framework is applicable to orebody modelling given geological boundaries, and 3D segmentation more generally, where there is a need to delineate spatial regions using mesh surfaces within a block model.


title: A simpler winning strategy for Sim

abstract: We give a simple human-playable winning strategy for the second player in the game of Sim.


title: Universal Successor Features for Transfer Reinforcement Learning

abstract: Transfer in Reinforcement Learning (RL) refers to the idea of applying knowledge gained from previous tasks to solving related tasks. Learning a universal value function (Schaul et al., 2015), which generalizes over goals and states, has previously been shown to be useful for transfer. However, successor features are believed to be more suitable than values for transfer (Dayan, 1993; Barreto et al.,2017), even though they cannot directly generalize to new goals. In this paper, we propose (1) Universal Successor Features (USFs) to capture the underlying dynamics of the environment while allowing generalization to unseen goals and (2) a flexible end-to-end model of USFs that can be trained by interacting with the environment. We show that learning USFs is compatible with any RL algorithm that learns state values using a temporal difference method. Our experiments in a simple gridworld and with two MuJoCo environments show that USFs can greatly accelerate training when learning multiple tasks and can effectively transfer knowledge to new tasks.