title: 3D Turbulent Reconnection: Theory, Tests and Astrophysical Implications

abstract: Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws.


title: Broadband Ultrahigh-Resolution chip-scale Scanning Soliton Dual-Comb Spectroscopy

abstract: We present a chip-scale scanning dual-comb spectroscopy (SDCS) approach for broadband ultrahigh-resolution spectral acquisition. SDCS uses Si3N4 microring resonators that generate two single soliton micro-combs spanning 37 THz (300 nm) on the same chip from a single 1550-nm laser, forming a high-mutual-coherence dual-comb. We realize continuous tuning of the dual-comb system over the entire optical span of 37.5 THz with high precision using integrated microheater-based wavelength trackers. This continuous wavelength tuning is enabled by simultaneous tuning of the laser frequency and the two single soliton micro-combs over a full free spectral range of the microrings. We measure the SDCS resolution to be 319+-4.6 kHz. Using this SDCS system, we perform the molecular absorption spectroscopy of H13C14N over its 2.3 THz (18 nm)-wide overtone band, and show that the massively parallel heterodyning offered by the dual-comb expands the effective spectroscopic tuning speed of the laser by one order of magnitude. Our chip-scale SDCS opens the door to broadband spectrometry and massively parallel sensing with ultrahigh spectral resolution.


title: Distributed and Inexact Proximal Gradient Method for Online Convex Optimization

abstract: This paper develops and analyzes an online distributed proximal-gradient method (DPGM) for time-varying composite convex optimization problems. Each node of the network features a local cost that includes a smooth strongly convex function and a non-smooth convex function, both changing over time. By coordinating through a connected communication network, the nodes collaboratively track the trajectory of the minimizers without exchanging their local cost functions. The DPGM is implemented in an online fashion, that is, in a setting where only a limited number of steps are implemented before the function changes. Moreover, the algorithm is analyzed in an inexact scenario, that is, with a source of additive noise, that can represent e.g. communication noise or quantization. It is shown that the tracking error of the online inexact DPGM is upper-bounded by a convergent linear system, guaranteeing convergence within a neighborhood of the optimal solution.


title: Ergodic Sum Rate Analysis of UAV-Based Relay Networks with Mixed RF-FSO Channels

abstract: Unmanned aerial vehicle (UAV)-based communications is a promising new technology that can add a wide range of new capabilities to the current network infrastructure. Given the flexibility, cost-efficiency, and convenient use of UAVs, they can be deployed as temporary base stations (BSs) for on-demand situations like BS overloading or natural disasters. In this work, a UAV-based communication system with radio frequency (RF) access links to the mobile users (MUs) and a free-space optical (FSO) backhaul link to the ground station (GS) is considered. In particular, the RF and FSO channels in this network depend on the UAV's positioning and (in)stability. The relative position of the UAV with respect to the MUs impacts the likelihood of a line-of-sight (LOS) connection in the RF link and the instability of the hovering UAV affects the quality of the FSO channel. Thus, taking these effects into account, we analyze the end-to-end system performance of networks employing UAVs as buffer-aided (BA) and non-buffer-aided (non-BA) relays in terms of the ergodic sum rate. Simulation results validate the accuracy of the proposed analytical derivations and reveal the benefits of buffering for compensation of the random fluctuations caused by the UAV's instability. Our simulations also show that the ergodic sum rate of both BA and non-BA UAV-based relays can be enhanced considerably by optimizing the positioning of the UAV. We further study the impact of the MU density and the weather conditions on the end-to-end system performance.


title: Plasmonic sensors based on funneling light through nanophotonic structures

abstract: We present a refractometric sensor realized as a stack of metallic gratings with subwavelength features and embedded within a low-index dielectric medium. Light is strongly confined through funneling mechanisms and excites resonances that sense the analyte medium. Two terminations of the structure are compared. One of them has a dielectric medium in contact with the analyte and exploits the selective spectral transmission of the structure. The other design has a metallic continuous layer that generates surface plasmon resonances at the metal/analyte interface. Both designs respond with narrow spectral features that are sensible to the change in the refractive index of the analyte and can be used for sensing biomedical samples.


title: Onset of Cosmic Reionization: Evidence of An Ionized Bubble Merely 680 Myrs after the Big Bang

abstract: While most of the inter-galactic medium (IGM) today is permeated by ionized hydrogen, it was largely filled with neutral hydrogen for the first 700 million years after the Big Bang. The process that ionized the IGM (cosmic reionization) is expected to be spatially inhomogeneous, with fainter galaxies playing a significant role. However, we still have only a few direct constraints on the reionization process. Here we report the first spectroscopic confirmation of two galaxies and very likely a third galaxy in a group (hereafter EGS77) at redshift z = 7.7, merely 680 Myrs after the Big Bang. The physical separation among the three members is < 0.7 Mpc. We estimate the radius of ionized bubble of the brightest galaxy to be about 1.02 Mpc, and show that the individual ionized bubbles formed by all three galaxies likely overlap significantly, forming a large yet localized ionized region, which leads to the spatial inhomogeneity in the reionization process. It is striking that two of three galaxies in EGS77 are quite faint in the continuum, thanks to our selection of reionizing sources using their Lyman-alpha line emission. Indeed, one is the faintest spectroscopically confirmed galaxy yet discovered at such high redshifts. Our observations provide direct constraints in the process of cosmic reionization, and allow us to investigate the properties of sources responsible for reionizing the universe.


title: Drug delivery in tumors is enhanced by bacterial proteolytic activity in a size and binding affinity dependent manner -A mechanistic understanding

abstract: The use of bacteria has been attractive to cancer researchers as drug delivery vehicle because motile bacteria are able to penetrate in tumors. In particular, the combination of therapeutic bacteria and conventional chemotherapy leads to dramatically high anti-tumor efficay. However, the mechanisms of the synergy, in part, remain unclear. To aim for understanding the mechanisms of the synergy of the combination therapy, simultaneous delivery of C. novyi-NT and chemotherapeutic agents in tumors is mathematically modeled from porous media approach. Simulated doxorubicin concentration in tumors after Doxil administration with or without bacteria agreed reasonably well with experimental literature. The simulated doxorubicin concentration in tumors by the combination of Doxil and C. novyi-NT is over twice higher than that of Doxil alone, as observed in previous experimental literature. This enhanced concentration is because of the degradation of extracellular matrix of collagen by bacterial proteolytic activity, which reduced interstitial fluid pressure in tumors by increasing hydraulic conductivity of interstitium, and thus increases convection through vessel walls. Additionally, solid stress alleviation caused by collagen degradation increases vessel density by decompressing blood vessels. On the other hand, the simulated doxorubicin concentration in tumors for non-liposomal doxorubicin is not enhanced by C. novyi-NT because vascular permeability of free-doxorubicin is larger than Doxil, and thus, increased but relatively small convection across vessel walls is outweighed by the efflux due to increased interstitial flow. A strategy to further enhance this combination therapy is discussed with sensitivity analysis.


title: Stahl--Totik regularity for continuum Schr\"odinger operators

abstract: We develop a theory of regularity for continuum Schr\"odinger operators based on the Martin compactification of the complement of the essential spectrum. This theory is inspired by Stahl--Totik regularity for orthogonal polynomials, but requires a different approach, since Stahl--Totik regularity is formulated in terms of the potential theoretic Green function with a pole at $\infty$, logarithmic capacity, and the equilibrium measure for the support of the measure, notions which do not extend to the case of unbounded spectra. For any half-line Schr\"odinger operator with a bounded potential (in a locally $L^1$ sense), we prove that its essential spectrum obeys the Akhiezer--Levin condition, and moreover, that the Martin function at $\infty$ obeys the two-term asymptotic expansion $\sqrt{-z} + \frac{a}{2\sqrt{-z}} + o(\frac 1{\sqrt{-z}})$ as $z \to -\infty$. The constant $a$ in that expansion plays the role of a renormalized Robin constant suited for Schr\"odinger operators and enters a universal inequality $a \le \liminf_{x\to\infty} \frac 1x \int_0^x V(t)dt$. This leads to a notion of regularity, with connections to the root asymptotics of Dirichlet solutions and zero counting measures. We also present applications to decaying and ergodic potentials.


title: Estimation of the shear strength of heterogeneous granular soils with matrix

abstract: The determination of the mechanical properties of soils containing particles larger than the allowable size by standard laboratory equipment is problematic. It is indeed necessary to eliminate the coarsest fraction to carry out the tests, which poses a problem of representativeness of the results at the scale of the structure. Parallel gradation and scalping methods are presented as a means of estimating the shear strength of heterogeneous granular soils from tests on their finer fraction. This study focuses on the effect of fines content on the estimation by these methods, which is not well known. The results showed that the parallel gradation method can predict the angle of friction of the initial material despite a high percentage of fines. However, the cohesion of the initial material was overestimated. As for the scalping method proposed by Fragaszy et al., the higher fines content in the modelled material resulted in an overestimation of the cohesion and an underestimation of the friction angle of the initial material.


title: BAT AGN Spectroscopic Survey -- XV: The High Frequency Radio Cores of Ultra-hard X-ray Selected AGN

abstract: We have conducted 22 GHz radio imaging at 1" resolution of 100 low-redshift AGN selected at 14-195 keV by the Swift-BAT. We find a radio core detection fraction of 96%, much higher than lower-frequency radio surveys. Of the 96 radio-detected AGN, 55 have compact morphologies, 30 have morphologies consistent with nuclear star formation, and 11 have sub-kpc to kpc-scale jets. We find that the total radio power does not distinguish between nuclear star formation and jets as the origin of the radio emission. For 87 objects, we use optical spectroscopy to test whether AGN physical parameters are distinct between radio morphological types. We find that X-ray luminosities tend to be higher if the 22 GHz morphology is jet-like, but find no significant difference in other physical parameters. We find that the relationship between the X-ray and core radio luminosities is consistent with the $L_R/L_X \sim 10^{-5}$ of coronally active stars. We further find that the canonical fundamental planes of black hole activity systematically over-predict our radio luminosities, particularly for objects with star formation morphologies.


title: Predicting competitions by combining conditional logistic regression and subjective Bayes: An Academy Awards case study

abstract: Predicting the outcome of elections, sporting events, entertainment awards, and other competitions has long captured the human imagination. Such prediction is growing in sophistication in these areas, especially in the rapidly growing field of data-driven journalism intended for a general audience as the availability of historical information rapidly balloons. Providing statistical methodology to probabilistically predict competition outcomes faces two main challenges. First, a suitably general modeling approach is necessary to assign probabilities to competitors. Second, the modeling framework must be able to accommodate expert opinion, which is usually available but difficult to fully encapsulate in typical data sets. We overcome these challenges with a combined conditional logistic regression/subjective Bayes approach. To illustrate the method, we re-analyze data from a recent Time.com piece in which the authors attempted to predict the 2019 Best Picture Academy Award winner using standard logistic regression. Towards engaging and educating a broad readership, we discuss strategies to deploy the proposed method via an online application.


title: The Robotilter: An Automated Lens / CCD Alignment System for the Evryscope

abstract: Camera lenses are increasingly used in wide-field astronomical surveys due to their high performance, wide field-of-view (FOV) unreachable from traditional telescope optics, and modest cost. The machining and assembly tolerances for commercially available optical systems cause a slight misalignment (tilt) between the lens and CCD, resulting in PSF degradation. We have built an automated alignment system (Robotilters) to solve this challenge, optimizing 4 degrees of freedom - 2 tilt axes, a separation axis (the distance between the CCD and lens), and the lens focus (the built-in focus of the lens by turning the lens focusing ring which moves the optical elements relative to one another) in a compact and low-cost package. The Robotilters remove tilt and optimize focus at the sub 10 micron level, are completely automated, take 2 hours to run, and remain stable for multiple years once aligned. The Robotilters were built for the Evryscope telescope (a 780 MPix 22-camera array with an 8150 sq.deg. field of view and continuous 2-minute cadence) designed to detect short timescale events across extremely large sky areas simultaneously. Variance in quality across the image field, especially the corners and edges compared to the center, is a significant challenge in wide-field astronomical surveys like the Evryscope. The individual star PSFs (which typically extend only a few pixels) are highly susceptible to slight increases in optical aberrations in this situation. The Robotilter solution resulted in a limiting magnitude improvement of .5 mag in the center of the image and 1.0 mag in the corners for typical Evryscope cameras, with less distorted and smaller PSFs (half the extent in the corners and edges in many cases). In this paper we describe the Robotilter mechanical and software design, camera alignment results, long term stability, and image improvement.


title: Moser-Tardos resampling algorithm, entropy compression method and the subset gas

abstract: We establish a connection between the entropy compression method and the Moser-Tardos algorithmic version of the Lov\'asz local lemma through the cluster expansion of the subset gas. We also show that the Moser-Tardos resampling algorithm and the entropy compression bactracking algorithm produce identical bounds.


title: Standing waves of the quintic NLS equation on the tadpole graph

abstract: The tadpole graph consists of a circle and a half-line attached at a vertex. We analyze standing waves of the nonlinear Schr\"{o}dinger equation with quintic power nonlinearity equipped with the Neumann-Kirchhoff boundary conditions at the vertex. The profile of the standing wave with the frequency $\omega\in (-\infty,0)$ is characterized as a global minimizer of the quadratic part of energy constrained to the unit sphere in $L^6$. The set of minimizers includes the set of ground states of the system, which are the global minimizers of the energy at constant mass ($L^2$-norm), but it is actually wider. While ground states exist only for a certain interval of masses, the standing waves exist for every $\omega \in (-\infty,0)$ and correspond to a bigger interval of masses. It is shown that there exist critical frequencies $\omega_0$ and $\omega_1$ such that the standing waves are the ground states for $\omega \in [\omega_0,0)$, local minimizers of the energy at constant mass for $\omega \in (\omega_1,\omega_0)$, and saddle points of the energy at constant mass for $\omega \in (-\infty,\omega_1)$. Proofs make use of both the variational methods and the analytical theory for differential equations.


title: Exponential bounds for inhomogeneous random graphs in a Gaussian case

abstract: Rank 1 inhomogeneous random graphs are a natural generalization of Erd\H{o}s R\'enyi random graphs. In this generalization each node is given a weight. Then the probability that an edge is present depends on the product of the weights of the nodes it is connecting. In this article, we give precise and uniform exponential bounds on the size, weight and surplus of rank 1 inhomogeneous random graphs where the weights of the nodes behave like a random variable with finite fourth moment. We focus on the case where the mean degree of a random node is slightly larger than 1, we call that case the barely supercritical regime. These bounds will be used in follow up articles to study a general class of random minimum spanning trees. They are also of independent interest since they show that these inhomogeneous random graphs behave like Erd\H{o}s R\'enyi random graphs even in a barely supercritical regime. The proof relies on novel concentration bounds for sampling without replacement and a careful study of the exploration process.


title: Double parton distributions of the pion in the NJL model

abstract: We evaluate the valence double parton distribution (dPDF) of the pion in the Nambu--Jona-Lasinio model. At the low-energy quark-model scale and in the chiral limit a particularly simple factorized form $D(x_1,x_2, \vec{q}) = \delta(1-x_1-x_2) F(\vec{q})$ follows, where $x_{1,2}$ denote the longitudinal momentum fractions of the valence quark and antiquark, and $\vec{q}$ is their relative transverse momentum. For $\vec{q}=\vec{0}$ our result complies to the Gaunt-Sterling sum rules. We carry out the necessary dDGLAP evolution to higher scales via the Mellin moments and explore its impact on the correlation defined as the ratio of dPDF to the product of single parton distributions, $D(x_1,x_2, \vec{q}=\vec{0})/D(x_1)D(x_2)$. Since the ratios of the valence Mellin moments $\langle x_1^n x_2^m \rangle / \langle x_1^n \rangle \langle x_2^m \rangle $ are invariants of the dDGLAP evolution, they may serve as robust measures of these correlations. Model predictions, which can be tested in the upcoming lattice simulations, are provided. We also discuss the transverse form factor related to the dPDF of the pion.


title: Improving Grid Computing Performance by Optimally Reducing Checkpointing Effect

abstract: Grid computing is a collection of computer resources that are gathered together from various areas to give computational resources such as storage, data or application services. This is to permit clients to access this huge measure of processing resources without the need to know where these might be found and what technology such as, hardware equipment and operating system was used. Dependability and performance are among the key difficulties faced in a grid computing environment. Various systems have been proposed in the literature to handle recouping from resource failure in Grid computing environment. One case of such system is checkpointing. Checkpointing is a system that endures faults when resources failed. Checkpointing method has the upside of lessening the work lost because of resource faults. However, checkpointing presents a huge runtime overhead. In this paper, we propose an improved checkpointing system to bring down runtime overhead. A replica is added to ensure the availability of resources. This replicates all checkpointing files to other machines as opposed to having dedicated checkpointing machine. The results of simulation employing GridSim noted that retaining the number of resources fixed and changing the number of gridlets, gains of up to 11%, 9%, and 11% on makespan, throughput, and turnaround time respectively were realized while changing the number of resources and preserving the number of gridlets fixed, increases of up to 11%, 8%, and 9% on makespan, throughput, and turnaround time respectively, were realized


title: Off-shell Ward identities for N-gluon amplitudes

abstract: Off-shell Ward identities in non-abelian gauge theory continue to be a subject of active research, since they are, in general, inhomogeneous and their form depends on the chosen gauge-fixing procedure. For the three-gluon and four-gluon vertices, it is known that a relatively simple form of the Ward identity can be achieved using the pinch technique or, equivalently, the background-field method with quantum Feynman gauge. The latter is also the gauge-fixing underlying the string-inspired formalism, and here we use this formalism to derive the corresponding form of the Ward identity for the one-loop N - gluon amplitudes.


title: Constraining the gauge and scalar sectors of the doublet left-right symmetric model

abstract: We consider a left-right symmetric extension of the Standard Model where the spontaneous breakdown of the left-right symmetry is triggered by doublets. The electroweak $\rho$ parameter is protected from large corrections in this Doublet Left-Right Model (DLRM), contrary to the triplet case. This allows in principle for more diverse patterns of symmetry breaking. We consider several constraints on the gauge and scalar sectors of DLRM: the unitarity of scattering processes involving gauge bosons with longitudinal polarisations, the radiative corrections to the muon $\Delta r$ parameter and the electroweak precision observables measured at the $Z$ pole and at low energies. Combining these constraints within the frequentist CKMfitter approach, we see that the fit pushes the scale of left-right symmetry breaking up to a few TeV, while favouring an electroweak symmetry breaking triggered not only by the $ SU(2)_L \times SU(2)_R $ bi-doublet, which is the case most commonly considered in the literature, but also by the $ SU(2)_L $ doublet.


title: Tuning Multigrid Methods with Robust Optimization

abstract: Local Fourier analysis is a useful tool for predicting and analyzing the performance of many efficient algorithms for the solution of discretized PDEs, such as multigrid and domain decomposition methods. The crucial aspect of local Fourier analysis is that it can be used to minimize an estimate of the spectral radius of a stationary iteration, or the condition number of a preconditioned system, in terms of a symbol representation of the algorithm. In practice, this is a "minimax" problem, minimizing with respect to solver parameters the appropriate measure of work, which involves maximizing over the Fourier frequency. Often, several algorithmic parameters may be determined by local Fourier analysis in order to obtain efficient algorithms. Analytical solutions to minimax problems are rarely possible beyond simple problems; the status quo in local Fourier analysis involves grid sampling, which is prohibitively expensive in high dimensions. In this paper, we propose and explore optimization algorithms to solve these problems efficiently. Several examples, with known and unknown analytical solutions, are presented to show the effectiveness of these approaches.


title: Towards Scalable Dataframe Systems

abstract: Dataframes are a popular abstraction to represent, prepare, and analyze data. Despite the remarkable success of dataframe libraries in Rand Python, dataframes face performance issues even on moderately large datasets. Moreover, there is significant ambiguity regarding dataframe semantics. In this paper we lay out a vision and roadmap for scalable dataframe systems. To demonstrate the potential in this area, we report on our experience building MODIN, a scaled-up implementation of the most widely-used and complex dataframe API today, Python's pandas. With pandas as a reference, we propose a simple data model and algebra for dataframes to ground discussion in the field. Given this foundation, we lay out an agenda of open research opportunities where the distinct features of dataframes will require extending the state of the art in many dimensions of data management. We discuss the implications of signature data-frame features including flexible schemas, ordering, row/column equivalence, and data/metadata fluidity, as well as the piecemeal, trial-and-error-based approach to interacting with dataframes.


title: Judicial Favoritism of Politicians: Evidence from Small Claims Court

abstract: Multiple studies have documented racial, gender, political ideology, or ethnical biases in comparative judicial systems. Supplementing this literature, we investigate whether judges rule cases differently when one of the litigants is a politician. We suggest a theory of power collusion, according to which judges might use rulings to buy cooperation or threaten members of the other branches of government. We test this theory using a sample of small claims cases in the state of S\~ao Paulo, Brazil, where no collusion should exist. The results show a negative bias of 3.7 percentage points against litigant politicians, indicating that judges punish, rather than favor, politicians in court. This punishment in low-salience cases serves as a warning sign for politicians not to cross the judiciary when exercising checks and balances, suggesting yet another barrier to judicial independence in development settings.


title: The calculation of the neutrinoless double-beta decay matrix element within the realistic shell model

abstract: We approach the calculation of the nuclear matrix element of the neutrinoless double-beta decay process, considering the light-neutrino-exchange channel, by way of the realistic shell model. To this end, we start from a realistic nucleon-nucleon potential and then derive the effective shell-model Hamiltonian and neutrinoless double-beta decay operator within the many-body perturbation theory. We focus on investigating the perturbative properties of the effective shell-model operator of such a decay process, aiming to establish the degree of reliability of our predictions. The contributions of the so-called short-range correlations and of the correction of Pauli-principle violations to the effective shell-model operator, the latter introduced in many-valence nucleon systems, are also taken into account. The subjects of our study are a few candidates to the neutrinoless double-beta decay detection, in a mass interval ranging from A=48 up to A=136, whose spin- and spin-isospin-dependent decay properties we have studied in previous works. Our results will be finally compared with shell-model calculations for the same set of nuclei.


title: Two-Level Transformer and Auxiliary Coherence Modeling for Improved Text Segmentation

abstract: Breaking down the structure of long texts into semantically coherent segments makes the texts more readable and supports downstream applications like summarization and retrieval. Starting from an apparent link between text coherence and segmentation, we introduce a novel supervised model for text segmentation with simple but explicit coherence modeling. Our model -- a neural architecture consisting of two hierarchically connected Transformer networks -- is a multi-task learning model that couples the sentence-level segmentation objective with the coherence objective that differentiates correct sequences of sentences from corrupt ones. The proposed model, dubbed Coherence-Aware Text Segmentation (CATS), yields state-of-the-art segmentation performance on a collection of benchmark datasets. Furthermore, by coupling CATS with cross-lingual word embeddings, we demonstrate its effectiveness in zero-shot language transfer: it can successfully segment texts in languages unseen in training.


title: Exploration of Interaction Techniques for Graph-based Modelling in Virtual Reality

abstract: Editing and manipulating graph-based models within immersive environments is largely unexplored and certain design activities could benefit from using those technologies. For example, in the case study of architectural modelling, the 3D context of Virtual Reality naturally matches the intended output product, i.e. a 3D architectural geometry. Since both the state of the art and the state of the practice are lacking, we explore the field of VR-based interactive modelling, and provide insights as to how to implement proper interactions in that context, with broadly available devices. We consequently produce several open-source software prototypes for manipulating graph-based models in VR.


title: Aleatoric and Epistemic Uncertainty with Random Forests

abstract: Due to the steadily increasing relevance of machine learning for practical applications, many of which are coming with safety requirements, the notion of uncertainty has received increasing attention in machine learning research in the last couple of years. In particular, the idea of distinguishing between two important types of uncertainty, often refereed to as aleatoric and epistemic, has recently been studied in the setting of supervised learning. In this paper, we propose to quantify these uncertainties with random forests. More specifically, we show how two general approaches for measuring the learner's aleatoric and epistemic uncertainty in a prediction can be instantiated with decision trees and random forests as learning algorithms in a classification setting. In this regard, we also compare random forests with deep neural networks, which have been used for a similar purpose.


title: Submodular Matroid Secretary Problem with Shortlists

abstract: In the matroid secretary problem, the elements of a matroid $\mathcal{M}$ arrive in random order. Once we observe an item we need to irrevocably decide whether or not to accept it. The set of selected elements should form an independent set of the matroid. The goal is to maximize the total sum of the values assigned to these elements. We introduce a version of this problem motivated by the shortlist model in [Agrawal et al.]. In this setting, the algorithm is allowed to choose a subset of items as part of a shortlist, possibly more than $k=rk(\mathcal{M})$ items. Then, after seeing the entire input, the algorithm can choose an independent subset from the shortlist. Furthermore we generalize the objective function to any monotone submodular function. Is there an online algorithm achieve a constant competitive ratio using a shortlist of size $O(k)$? We design an algorithm that achieves a $\frac{1}{2}(1-1/e^2-\epsilon-O(1/k))$ competitive ratio for any constant $\epsilon>0$, using a shortlist of size $O(k)$. This is especially surprising considering that the best known competitive ratio for the matroid secretary problem is $O(\log \log k)$. An important application of our algorithm is for the random order streaming of submodular functions. We show that our algorithm can be implemented in the streaming setting using $O(k)$ memory. It achieves a $\frac{1}{2}(1-1/e^2-\epsilon-O(1/k))$ approximation. The previously best known approximation ratio for streaming submodular maximization under matroid constraint is 0.25 (adversarial order) due to [Feldman et al.], [Chekuri et al.], and [Chakrabarti et al.]. Moreover, we generalize our results to the case of p-matchoid constraints and give a $\frac{1}{p+1}(1-1/e^{p+1}-\epsilon-O(1/k))$ approximation using $O(k)$ memory, which asymptotically approaches the best known offline guarantee $\frac{1}{p+1}$ [Nemhauser et al.]


title: A method to deal with the critical case in stochastic population dynamics

abstract: In numerous papers, the behaviour of stochastic population models is investigated through the sign of a real quantity which is the growth rate of the population near the extinction set. In many cases, it is proven that when this growth rate is positive, the process is persistent in the long run, while if it is negative, the process converges to extinction. However, the critical case when the growth rate is null is rarely treated. The aim of this paper is to provide a method that can be applied in many situations to prove that in the critical case, the process congerves in temporal average to extinction. A number of applications are given, for Stochastic Differential Equations and Piecewise Deterministic Markov Processes modelling prey-predator, epidemilogical or structured population dynamics.


title: Non Gaussian information of heterogeneity in Soft Matter

abstract: Heterogeneity in dynamics in the form of non-Gaussian molecular displacement distributions appears ubiquitously in soft matter. We address the quantification of such heterogeneity using an information-theoretic measure of the distance between the actual displacement distribution and its nearest Gaussian estimation. We explore the usefulness of this measure in two generic scenarios of random walkers in heterogeneous media. We show that our proposed measure leads to a better quantification of non-Gaussianity than the conventional ones based on moment ratios.


title: Electronic, optical and thermoelectric properties of sodium pnictogen chalcogenides: A first principles study

abstract: Ternary chalcogenides have been of recent investigation for applications in photovoltaic and thermoelectric devices. We study the structural, electronic, optical, and thermoelectric properties of nine ternary chalcogenides, NaAX$_{2}$, where A represents pnictogens (As, Sb, and Te) and X represents chalcogens (S, Se, and Te). Calculations based on density functional theory yield the following results: (i) phonon dispersion curves predict three of the compounds, NaAsS$_{2}$, NaSbS$_{2}$, and NaSbSe$_{2}$, to be dynamically stable in the monoclinic, C2/c, structure, (ii) the layered atomistic configuration causes the corresponding electronic and optical properties to display a high degree of anisotropy, (iii) A-X electronic bonding features vary significantly with structural distortions arising from atomic size mismatch, therefore directly influencing stability, (iv) strong absorption is observed in the stable compounds, with coefficients ranging from 10$^{4}$ to 10$^{5}$ cm$^{-1}$ in the visible-UV range, and (v) remarkably high Seebeck coefficients exceeding 500 ${\mu}$V/K at carrier concentrations commonly achieved in such materials. From these results, we conclude that NaAsS$_{2}$, NaSbS$_{2}$, and NaSbSe$_{2}$ are suitable candidates for both photovoltaic, particularly in tandem solar cells, and thermoelectric applications. Experimental synthesis and verification are suggested.


title: Second-Order Cone Relaxations of the Optimal Power Flow for Active Distribution Grids

abstract: Convex relaxations of the AC Optimal Power Flow (OPF) problem are essential not only for identifying the globally optimal solution but also for enabling the use of OPF formulations in Bilevel Programming and Mathematical Programs with Equilibrium Constraints (MPEC), which are required for solving problems such as the coordination between transmission and distribution system operator (TSO/DSO) or optimal network investment. Focusing on active distribution grids and radial networks, this paper introduces a framework that collects and compares, for the first time to our knowledge, the performance of the most promising convex OPF formulations for practical applications. Our goal is to establish a solid basis that will inform the selection of the most appropriate algorithm for different applications. This paper (i) introduces a unified mathematical and simulation framework, (ii) extends existing methods to retain exactness in a wider number of cases and (iii) consider reactive power injections. We conduct simulations on the IEEE 34 and 123 radial test feeders with distributed energy resources (DERs), using yearly solar irradiation and load data.


title: The Orion Protostellar Explosion and Runaway Stars Revisited: Stellar Masses, Disk Retention, and an Outflow from BN

abstract: The proper motions of the three stars ejected from Orion's OMC1 cloud core are combined with the requirement that their center of mass is gravitationally bound to OMC1 to show that radio source I (Src I) is likely to have a mass around 15 Solar masses consistent with recent measurements. Src I, the star with the smallest proper motion, is suspected to be either an AU-scale binary or a protostellar merger remnant produced by a dynamic interaction ~550 years ago. Near-infrared 2.2 um images spanning ~21 years confirm the ~55 km/s motion of `source x' (Src x) away from the site of stellar ejection and point of origin of the explosive OMC1 protostellar outflow. The radial velocities and masses of the Becklin-Neugebauer (BN) object and Src I constrain the radial velocity of Src x to be V_{LSR} = -28 +/-10 km/s . Several high proper-motion radio sources near BN, including Zapata 11 ([ZRK2004] 11) and a diffuse source near IRc 23, may trace a slow bipolar outflow from BN. The massive disk around Src I is likely the surviving portion of a disk that existed prior to the stellar ejection. Though highly perturbed, shocked, and re-oriented by the N-body interaction, enough time has elapsed to allow the disk to relax with its spin axis roughly orthogonal to the proper motion.


title: An eco-epidemiological model with general functional response of predator to prey

abstract: We consider a nonautonomous eco-epidemiological model with general functions for predation on infected and uninfected preys as well as general functions associated to the vital dynamics of the susceptible prey and predator populations. We obtain persistence and extinction results for the infected prey based on assumptions on systems related to the dynamics in the absence of infected preys. We apply our results to eco-epidemiological models constructed from predator-prey models existent in the literature. Some illustrative simulation is undertaken.


title: Low-energy theorem for $\gamma\to 3\pi$: surface terms against $\pi a_1$-mixing

abstract: We reconsider the contribution due to $\pi a_1$-mixing to the anomalous $\gamma\to\pi^+\pi^0\pi^-$ amplitude from the standpoint of the low-energy theorem $F^{\pi}=e f_\pi^2 F^{3\pi}$, which relates the electromagnetic form factor $F_{\pi^0\to\gamma\gamma}=F^\pi$ with the form factor $F_{\gamma\to\pi^+\pi^0\pi^-}=F^{3\pi}$ both taken at vanishing momenta of mesons. Our approach is based on a recently proposed covariant diagonalization of $\pi a_1$-mixing within a standard effective QCD-inspired meson Lagrangian obtained in the framework of the Nambu-Jona-Lasinio model. We show that the two surface terms appearing in the calculation of the anomalous triangle quark diagrams or AVV- and AAA-type amplitudes are uniquely fixed by this theorem. As a result, both form factors $F^\pi$ and $F^{3\pi}$ are not affected by the $\pi a_1$-mixing, but the concept of vector meson dominance (VMD) fails for $\gamma\to\pi^+\pi^0\pi^-$.


title: On phase asymmetries in oscillatory pipe flow

abstract: We present results from direct numerical simulations (DNS) of oscillatory pipe flow at several dimensionless frequencies Wo = {6.5, 13, 26} and one fixed shear Reynolds number Re_{\tau} = 1440. Starting from a fully-developed turbulent velocity field at that Re_{\tau}, the oscillatory flow either relaminarises or reaches a conditionally turbulent or strongly asymmetric state depending on Wo. The numerical method is validated by demonstrating excellent agreement of our DNS results with experimental data and analytical predictions from literature for the limiting cases of non-oscillating but turbulent and oscillating but laminar pipe flow. For an oscillating turbulent pipe flow we further found a very good agreement between qualitative descriptions of the characteristic flow features observed in experiments and our DNS. Here, we focus on the observation of a strongly asymmetric behaviour between the positive and the negative half-cycles of the oscillatory pipe flow at Wo=6.5.


title: Characterizing the Uncertainty in Cluster Magnetic Fields derived from Rotation Measures

abstract: Magnetic fields play vital roles in intracluster media (ICMs), but estimating their strengths and distributions from observations is a major challenge. Faraday rotation measures (RMs) are widely applied to this task, so it is critical to understand inherent uncertainties in RM analysis. In this paper, we seek to characterize those uncertainties given the types of information available today, independent of the specific technique used. We conduct synthetic RM observations through the ICM of a galaxy cluster drawn from an MHD cosmological simulation in which the magnetic field is known. We analyze the synthetic RM observations using an analytical formalism based on commonly used model assumptions allowing us to relate model physical variables to outcome uncertainties. Despite the simplicity of some assumptions, and unknown physical parameters, we are able to extract an approximate magnitude of the central magnetic field within an apparently irreducible uncertain factor approximately 3. Principal, largely irreducible, uncertainties come from the unknown depth along the line of sight of embedded polarized sources, the lack of robust coherence lengths from area-constrained polarization sampling, and the unknown scaling between ICM electron density and magnetic field strength. The RM-estimated central magnetic field strengths span more than an order of magnitude including the full range of synthetic experiments.


title: Optimization of Mean-field Spin Glasses

abstract: Mean-field spin glasses are families of random energy functions (Hamiltonians) on high-dimensional product spaces. In this paper we consider the case of Ising mixed $p$-spin models, namely Hamiltonians $H_N:\Sigma_N\to {\mathbb R}$ on the Hamming hypercube $\Sigma_N = \{\pm 1\}^N$, which are defined by the property that $\{H_N({\boldsymbol \sigma})\}_{{\boldsymbol \sigma}\in \Sigma_N}$ is a centered Gaussian process with covariance ${\mathbb E}\{H_N({\boldsymbol \sigma}_1)H_N({\boldsymbol \sigma}_2)\}$ depending only on the scalar product $\langle {\boldsymbol \sigma}_1,{\boldsymbol \sigma}_2\rangle$. The asymptotic value of the optimum $\max_{{\boldsymbol \sigma}\in \Sigma_N}H_N({\boldsymbol \sigma})$ was characterized in terms of a variational principle known as the Parisi formula, first proved by Talagrand and, in a more general setting, by Panchenko. The structure of superlevel sets is extremely rich and has been studied by a number of authors. Here we ask whether a near optimal configuration ${\boldsymbol \sigma}$ can be computed in polynomial time. We develop a message passing algorithm whose complexity per-iteration is of the same order as the complexity of evaluating the gradient of $H_N$, and characterize the typical energy value it achieves. When the $p$-spin model $H_N$ satisfies a certain no-overlap gap assumption, for any $\varepsilon>0$, the algorithm outputs ${\boldsymbol \sigma}\in\Sigma_N$ such that $H_N({\boldsymbol \sigma})\ge (1-\varepsilon)\max_{{\boldsymbol \sigma}'} H_N({\boldsymbol \sigma}')$, with high probability. The number of iterations is bounded in $N$ and depends uniquely on $\varepsilon$. More generally, regardless of whether the no-overlap gap assumption holds, the energy achieved is given by an extended variational principle, which generalizes the Parisi formula.


title: Convergence and limits of finite trees

abstract: Motivated by the work of Lov\'asz and Szegedy on the convergence and limits of dense graph sequences, we investigate the convergence and limits of finite trees with respect to sampling in normalized distance. Based on separable real trees, we introduce the notion of a dendron and show that the limits of finite trees are exactly the dendrons. We also prove that the limit dendron is unique.


title: Coulomb Blockade Effects in a Topological Insulator Grown on a High-Tc Cuprate Superconductor

abstract: The evidence for proximity-induced superconductivity in heterostructures of topological insulators and high-Tc cuprates has been intensely debated. We use molecular beam epitaxy to grow thin films of topological insulator Bi2Te3 on a cuprate Bi2Sr2CaCu2O8+x, and study the surface of Bi2Te3 using low-temperature scanning tunneling microscopy and spectroscopy. In few unit-cell thick Bi2Te3 films, we find a V-shaped gap-like feature at the Fermi energy in dI/dV spectra. By reducing the coverage of Bi2Te3 films to create nanoscale islands, we discover that this spectral feature dramatically evolves into a much larger hard gap, which can be understood as a Coulomb blockade gap. This conclusion is supported by the evolution of dI/dV spectra with the lateral size of Bi2Te3 islands, as well as by topographic measurements that show an additional barrier separating Bi2Te3 and Bi2Sr2CaCu2O8+x. We conclude that the prominent gap-like feature in dI/dV spectra in Bi2Te3 films is not a proximity-induced superconducting gap. Instead, it can be explained by Coulomb blockade effects, which take into account additional resistive and capacitive coupling at the interface. Our experiments provide a fresh insight into the tunneling measurements of complex heterostructures with buried interfaces.


title: Selfish Algorithm and Emergence of Collective Intelligence

abstract: We propose a model for demonstrating spontaneous emergence of collective intelligent behavior from selfish individual agents. Agents' behavior is modeled using our proposed selfish algorithm ($SA$) with three learning mechanisms: reinforced learning ($SAL$), trust ($SAT$) and connection ($SAC$). Each of these mechanisms provides a distinctly different way an agent can increase the individual benefit accrued through playing the prisoner's dilemma game ($PDG$) with other agents. The $SA$ provides a generalization of the self-organized temporal criticality ($SOTC$) model and shows that self-interested individuals can simultaneously produce maximum social benefit from their decisions. The mechanisms in the $SA$ are self-tuned by the internal dynamics and without having a pre-established network structure. Our results demonstrate emergence of mutual cooperation, emergence of dynamic networks, and adaptation and resilience of social systems after perturbations. The implications and applications of the $SA$ are discussed.


title: The generic combinatorial simplex

abstract: We employ projective Fra\"iss\'e theory to define the "generic combinatorial $n$-simplex" as the pro-finite, simplicial complex that is canonically associated with a family of simply defined selection maps between finite triangulations of the simplex. The generic combinatorial $n$-simplex is a combinatorial object that can be used to define the geometric realization of a simplicial complex without any reference to the Euclidean space. It also reflects dynamical properties of its homeomorphism group down to finite combinatorics. As part of our study of the generic combinatorial simplex, we define and prove results on domination closure for Fra\"iss\'e classes, and we develop further the theories of stellar moves and cellular maps. We prove that the domination closure of selection maps contains the class of face-preserving simplicial maps that are cellular on each face of the $n$-simplex and is contained in the class of simplicial, face-preserving near-homeomorphisms. Under the PL-Poincar\'e conjecture, this gives a characterization of the domination closure of selections.


title: Inner privacy of conscious experiences and quantum information

abstract: The human mind is constituted by inner, subjective, private, first-person conscious experiences that cannot be measured with physical devices or observed from an external, objective, public, third-person perspective. The qualitative, phenomenal nature of conscious experiences also cannot be communicated to others in the form of a message composed of classical bits of information. Because in a classical world everything physical is observable and communicable, it is a daunting task to explain how an empirically unobservable, incommunicable consciousness could have any physical substrates such as neurons composed of biochemical molecules, water, and electrolytes. The challenges encountered by classical physics are exemplified by a number of thought experiments including the inverted qualia argument, the private language argument, the beetle in the box argument and the knowledge argument. These thought experiments, however, do not imply that our consciousness is nonphysical and our introspective conscious testimonies are untrustworthy. The principles of classical physics have been superseded by modern quantum physics, which contains two fundamentally different kinds of physical objects: unobservable quantum state vectors, which define what physically exists, and quantum operators (observables), which define what can physically be observed. Identifying consciousness with the unobservable quantum information contained by quantum physical brain states allows for application of quantum information theorems to resolve possible paradoxes created by the inner privacy of conscious experiences, and explains how the observable brain is constructed by accessible bits of classical information that are bound by Holevo's theorem and extracted from the physically existing quantum brain upon measurement with physical devices.


title: Colloidal interactions and unusual crystallization versus de-mixing of elastic multipoles formed by gold mesoflowers

abstract: Colloidal interactions in nematic liquid crystals can be described as interactions between elastic multipoles that depend on particle shape, topology, chirality, boundary conditions and induced topological defects. Here, we describe a nematic colloidal system consisting of mesostructures of gold capable of inducing elastic multipoles of different order. Elastic monopoles are formed by relatively large asymmetric mesoflower particles, for which gravity and elastic torque balancing yields monopole-type interactions. High-order multipoles are instead formed by smaller mesoflowers with a myriad of shapes corresponding to multipoles of different orders, consistent with our computer simulations based on free energy minimization. We reveal unexpected many-body interactions in this colloidal system, ranging from de-mixing of elastic monopoles to a zoo of unusual colloidal crystals formed by high-order multipoles like hexadecapoles. Our findings show that gold mesoflowers may serve as a designer toolkit for engineering colloidal interaction and self-assembly, potentially exceeding that in atomic and molecular systems.


title: Comment on "Dirac fermions in Som-Raychaudhuri space-time with scalar and vector potential and the energy momentum distributions"

abstract: We point out a misleading treatment and incorrect expressions in a recent paper published in this Journal [Eur. Phys. J. C (2019) 79: 541] regarding solutions for the Dirac equation in presence of scalar and vector potentials in a class of flat G\"odel-type space-time called Som-Raychaudhuri space-time. Following the appropriate procedure we obtain the solution for this system.


title: Room temperature skyrmions at zero field in exchange-biased ultrathin films

abstract: We demonstrate that magnetic skyrmions with a mean diameter around 60 nm can be stabilized at room temperature and zero external magnetic field in an exchange-biased Pt/Co/NiFe/IrMn multilayer stack. This is achieved through an advanced optimization of the multilayer stack composition in order to balance the different magnetic energies controlling the skyrmion size and stability. Magnetic imaging is performed both with magnetic force microscopy and scanning Nitrogen-Vacancy magnetometry, the latter providing unambiguous measurements at zero external magnetic field. In such samples, we show that exchange bias provides an immunity of the skyrmion spin texture to moderate external magnetic field, in the tens of mT range, which is an important feature for applications as memory devices. These results establish exchange-biased multilayer stacks as a promising platform towards the effective realization of memory and logic devices based on magnetic skyrmions.


title: Equal cost of computation for truth and falsity of experimental quantum propositions necessitates quantum parallel computing

abstract: Notwithstanding interest and excitement building around quantum computing in the last decades, a concise statement saying where this computing can truly help is still missing. As it is shown in the present paper, equal cost of computation for truth and falsity of experimental quantum propositions (required in order to infer a conclusion from a premise) cannot be achieved with classical computing. On the other hand, this equality might be realized with quantum parallel computing provided that the efficiency of such computing can be greater than 1.


title: Novel Wake-up Scheme for Energy-Efficient Low-Latency Mobile Devices in 5G Networks

abstract: Improved mobile device battery lifetime and latency mini-mization are critical requirements for enhancing the mobile broadband services and user experience. Long-term evolution (LTE) networks have adopted discontinuous reception (DRX) as the baseline solution for prolonged battery lifetime. However, in every DRX cycle, the mobile device baseband processing unit monitors and decodes the control signaling, and thus all instances without any actual data allocation leads to unnecessary energy consumption. This fact together with the long start-up and power-down times can prevent adopting frequent wake-up instants, which in turn leads to considerable latency. In this work,a novel wake-up scheme is described and studied, to tackle the trade-off between latency and battery lifetime in future 5G networks, seeking thus to facilitate an always-available experience, rather than always-on. Analytical and simulation-based results show that the proposed scheme is a promising approach to control the user plane latency and energy consumption, when the device is operating in the power saving mode. The aim of this article is to describe the overall wake-up system operating principle and the associated signaling methods,receiver processing solutions and essential implementation aspects. Additionally, the advantages compared to DRX-based systems are shown and demonstrated, through the analysis of the system energy-efficiency and latency characteristics, with special emphasis on future 5G-grade mobile devices


title: Local polynomial regression for pooled response data

abstract: We propose local polynomial estimators for the conditional mean of a continuous response when only pooled response data are collected under different pooling designs. Asymptotic properties of these estimators are investigated and compared. Extensive simulation studies are carried out to compare finite sample performance of the proposed estimators under various model settings and pooling strategies. We apply the proposed local polynomial regression methods to two real-life applications to illustrate practical implementation and performance of the estimators for the mean function.


title: Deep Learning-Based Intrusion Detection System for Advanced Metering Infrastructure

abstract: Smart grid is an alternative solution of the conventional power grid which harnesses the power of the information technology to save the energy and meet today's environment requirements. Due to the inherent vulnerabilities in the information technology, the smart grid is exposed to a wide variety of threats that could be translated into cyber-attacks. In this paper, we develop a deep learning-based intrusion detection system to defend against cyber-attacks in the advanced metering infrastructure network. The proposed machine learning approach is trained and tested extensively on an empirical industrial dataset which is composed of several attack categories including the scanning, buffer overflow, and denial of service attacks. Then, an experimental comparison in terms of detection accuracy is conducted to evaluate the performance of the proposed approach with Naive Bayes, Support Vector Machine, and Random Forest. The obtained results suggest that the proposed approaches produce optimal results comparing to the other algorithms. Finally, we propose a network architecture to deploy the proposed anomaly-based intrusion detection system across the Advanced Metering Infrastructure network. In addition, we propose a network security architecture composed of two types of Intrusion detection system types, Host and Network-based, deployed across the Advanced Metering Infrastructure network to inspect the traffic and detect the malicious one at all the levels.


title: A Performance Comparison of Data Mining Algorithms Based Intrusion Detection System for Smart Grid

abstract: Smart grid is an emerging and promising technology. It uses the power of information technologies to deliver intelligently the electrical power to customers, and it allows the integration of the green technology to meet the environmental requirements. Unfortunately, information technologies have its inherent vulnerabilities and weaknesses that expose the smart grid to a wide variety of security risks. The Intrusion detection system (IDS) plays an important role in securing smart grid networks and detecting malicious activity, yet it suffers from several limitations. Many research papers have been published to address these issues using several algorithms and techniques. Therefore, a detailed comparison between these algorithms is needed. This paper presents an overview of four data mining algorithms used by IDS in Smart Grid. An evaluation of performance of these algorithms is conducted based on several metrics including the probability of detection, probability of false alarm, probability of miss detection, efficiency, and processing time. Results show that Random Forest outperforms the other three algorithms in detecting attacks with higher probability of detection, lower probability of false alarm, lower probability of miss detection, and higher accuracy.


title: Fairness in Multi-agent Reinforcement Learning for Stock Trading

abstract: Unfair stock trading strategies have been shown to be one of the most negative perceptions that customers can have concerning trading and may result in long-term losses for a company. Investment banks usually place trading orders for multiple clients with the same target assets but different order sizes and diverse requirements such as time frame and risk aversion level, thereby total earning and individual earning cannot be optimized at the same time. Orders executed earlier would affect the market price level, so late execution usually means additional implementation cost. In this paper, we propose a novel scheme that utilizes multi-agent reinforcement learning systems to derive stock trading strategies for all clients which keep a balance between revenue and fairness. First, we demonstrate that Reinforcement learning (RL) is able to learn from experience and adapt the trading strategies to the complex market environment. Secondly, we show that the Multi-agent RL system allows developing trading strategies for all clients individually, thus optimizing individual revenue. Thirdly, we use the Generalized Gini Index (GGI) aggregation function to control the fairness level of the revenue across all clients. Lastly, we empirically demonstrate the superiority of the novel scheme in improving fairness meanwhile maintaining optimization of revenue.


title: Competitive equilibria between staking and on-chain lending

abstract: Proof of Stake (PoS) is a burgeoning Sybil resistance mechanism that aims to have a digital asset ("token") serve as security collateral in crypto networks. However, PoS has so far eluded a comprehensive threat model that encompasses both Byzantine attacks from distributed systems and financial attacks that arise from the dual usage of the token as a means of payment and a Sybil resistance mechanism. In particular, the existence of derivatives markets makes malicious coordination among validators easier to execute than in Proof of Work systems. We demonstrate that it is also possible for on-chain lending smart contracts to cannibalize network security in PoS systems. When the yield provided by these contracts is more attractive than the inflation rate provided from staking, stakers will tend to remove their staked tokens and lend them out, thus reducing network security. In this paper, we provide a simple stochastic model that describes how rational validators with varying risk preferences react to changes in staking and lending returns. For a particular configuration of this model, we provide a formal proof of a phase transition between equilibria in which tokens are predominantly staked and those in which they are predominantly lent. We further validate this emergent adversarial behavior (e.g. reduced staked token supply) with agent-based simulations that sample transitions under more realistic conditions. Our results illustrate that rational, non-adversarial actors can dramatically reduce PoS network security if block rewards are not calibrated appropriately above the expected yields of on-chain lending.


title: Estimation of the yield curve for Costa Rica using combinatorial optimization metaheuristics applied to nonlinear regression

abstract: The term structure of interest rates or yield curve is a function relating the interest rate with its own term. Nonlinear regression models of Nelson-Siegel and Svensson were used to estimate the yield curve using a sample of historical data supplied by the National Stock Exchange of Costa Rica. The optimization problem involved in the estimation process of model parameters is addressed by the use of four well known combinatorial optimization metaheuristics: Ant colony optimization, Genetic algorithm, Particle swarm optimization and Simulated annealing. The aim of the study is to improve the local minima obtained by a classical quasi-Newton optimization method using a descent direction. Good results with at least two metaheuristics are achieved, Particle swarm optimization and Simulated annealing. Keywords: Yield curve, nonlinear regression, Nelson-


title: Wide Neural Networks with Bottlenecks are Deep Gaussian Processes

abstract: There has recently been much work on the "wide limit" of neural networks, where Bayesian neural networks (BNNs) are shown to converge to a Gaussian process (GP) as all hidden layers are sent to infinite width. However, these results do not apply to architectures that require one or more of the hidden layers to remain narrow. In this paper, we consider the wide limit of BNNs where some hidden layers, called "bottlenecks", are held at finite width. The result is a composition of GPs that we term a "bottleneck neural network Gaussian process" (bottleneck NNGP). Although intuitive, the subtlety of the proof is in showing that the wide limit of a composition of networks is in fact the composition of the limiting GPs. We also analyze theoretically a single-bottleneck NNGP, finding that the bottleneck induces dependence between the outputs of a multi-output network that persists through extreme post-bottleneck depths, and prevents the kernel of the network from losing discriminative power at extreme post-bottleneck depths.


title: Mapping large-scale diffuse gamma-ray emission in 10-100 TeV band with Cherenkov telescopes

abstract: Measurement of diffuse gamma-ray emission from the Milky Way with Imaging Atmospheric Cherenkov Telescopes (IACT) is difficult because of high level of charged cosmic ray background and small field-of-view. We show that such a measurement is nevertheless possible in the energy band 10-100 TeV. The minimal charged particle background for IACTs is achieved in event selections used for the analyses of the cosmic ray electrons. Tight cuts on the event quality in these event selections allow to achieve the background level sufficiently low to allow the measurements of the diffuse Galactic gamma-ray flux above 10 TeV. We calculate the sensitivities of different types of IACT arrays for the Galactic diffuse emission measurements and compare them with the diffuse gamma-ray flux from different parts of the sky measured by the Fermi Large Area Telescope below 3 TeV and with the astrophysical neutrino signal measured by IceCube telescope. We show that deep exposure of existing IACT systems is sufficient for detection of the diffuse flux from all the Galactic Plane up to Galactic latitude |b|~ degrees. Medium Size Telescope array of CTA will be able to detect the diffuse flux up 30 degree Galactic latitude. Its sensitivity will be sufficient for detection of the gamma-ray counterpart of the Galactic component of the IceCube astrophysical neutrino signal above 10 TeV. We also propose that a dedicated IACT system composed of small but wide field-of-view telescopes could be used for mapping the 10-100 TeV the diffuse gamma-ray emission from across all the sky.


title: A Joint Fermi-GBM and LIGO/Virgo Analysis of Compact Binary Mergers From the First and Second Gravitational-wave Observing Runs

abstract: We present results from offline searches of Fermi Gamma-ray Burst Monitor (GBM) data for gamma-ray transients coincident with the compact binary coalescences observed by the gravitational-wave (GW) detectors Advanced LIGO and Advanced Virgo during their first and second observing runs. In particular, we perform follow-up for both confirmed events and low significance candidates reported in the LIGO/Virgo catalog GWTC-1. We search for temporal coincidences between these GW signals and GBM triggered gamma-ray bursts (GRBs). We also use the GBM Untargeted and Targeted subthreshold searches to find coincident gamma-rays below the on-board triggering threshold. This work implements a refined statistical approach by incorporating GW astrophysical source probabilities and GBM visibilities of LIGO/Virgo sky localizations to search for cumulative signatures of coincident subthreshold gamma-rays. All search methods recover the short gamma-ray burst GRB 170817A occurring ~1.7 s after the binary neutron star merger GW170817. We also present results from a new search seeking GBM counterparts to LIGO single-interferometer triggers. This search finds a candidate joint event, but given the nature of the GBM signal and localization, as well as the high joint false alarm rate of $1.1 \times 10^{-6}$ Hz, we do not consider it an astrophysical association. We find no additional coincidences.


title: Nucleosynthesis of "Light" Heavy Nuclei in Neutrino-driven Winds. Role of ($\alpha,n$) reactions

abstract: Neutrino-driven winds following core collapse supernovae have been proposed as a suitable site where the so-called light heavy elements (between Sr to Ag) can be synthetized. For moderately neutron-rich winds, ($\alpha,n$) reactions play a critical role in the weak r process, becoming the main mechanism to drive nuclear matter towards heavier elements. In this paper we summarize the sensitivity of network-calculated abundances to the astrophysical conditions, and to uncertainties in the ($\alpha,n$) reaction rates. A list of few ($\alpha,n$) reactions were identified to dominate the uncertainty in the calculated elemental abundances. Measurements of these reactions will allow to identify the astrophysical conditions of the weak r process by comparing calculated/observed abundances in r-limited stars.


title: McKean-Vlasov optimal control: limit theory and equivalence between different formulations

abstract: We study a McKean-Vlasov optimal control problem with common noise, in order to establish the corresponding limit theory, as well as the equivalence between different formulations, including the strong, weak and relaxed formulation. In contrast to the strong formulation, where the problem is formulated on a fixed probability space equipped with two Brownian filtrations, the weak formulation is obtained by considering a more general probability space with two filtrations satisfying an $(H)$-hypothesis type condition from the theory of enlargement of filtrations. When the common noise is uncontrolled, our relaxed formulation is obtained by considering a suitable controlled martingale problem. As for classical optimal control problems, we prove that the set of all relaxed controls is the closure of the set of all strong controls, when considered as probability measures on the canonical space. Consequently, we obtain the equivalence of the different formulations of the control problem, under additional mild regularity conditions on the reward functions. This is also a crucial technical step to prove the limit theory of the McKean-Vlasov control problem, that is to say proving that it consists in the limit of a large population control problem with common noise.


title: Learning Accurate Integer Transformer Machine-Translation Models

abstract: We describe a method for training accurate Transformer machine-translation models to run inference using 8-bit integer (INT8) hardware matrix multipliers, as opposed to the more costly single-precision floating-point (FP32) hardware. Unlike previous work, which converted only 85 Transformer matrix multiplications to INT8, leaving 48 out of 133 of them in FP32 because of unacceptable accuracy loss, we convert them all to INT8 without compromising accuracy. Tested on the newstest2014 English-to-German translation task, our INT8 Transformer Base and Transformer Big models yield BLEU scores that are 99.3% to 100% relative to those of the corresponding FP32 models. Our approach converts all matrix-multiplication tensors from an existing FP32 model into INT8 tensors by automatically making range-precision trade-offs during training. To demonstrate the robustness of this approach, we also include results from INT6 Transformer models.


title: A Probability Density Theory for Spin-Glass Systems

abstract: Spin-glass systems are universal models for representing many-body phenomena in statistical physics and computer science. High quality solutions of NP-hard combinatorial optimization problems can be encoded into low energy states of spin-glass systems. In general, evaluating the relevant physical and computational properties of such models is difficult due to critical slowing down near a phase transition. Ideally, one could use recent advances in deep learning for characterizing the low-energy properties of these complex systems. Unfortunately, many of the most promising machine learning approaches are only valid for distributions over continuous variables and thus cannot be directly applied to discrete spin-glass models. To this end, we develop a continuous probability density theory for spin-glass systems with arbitrary dimensions, interactions, and local fields. We show how our formulation geometrically encodes key physical and computational properties of the spin-glass in an instance-wise fashion without the need for quenched disorder averaging. We show that our approach is beyond the mean-field theory and identify a transition from a convex to non-convex energy landscape as the temperature is lowered past a critical temperature. We apply our formalism to a number of spin-glass models including the Sherrington-Kirkpatrick (SK) model, spins on random Erd\H{o}s-R\'enyi graphs, and random restricted Boltzmann machines.


title: Tracking the evolution of magmas from heterogeneous mantle sources to eruption

abstract: This contribution reviews the effects of source heterogeneities, melt-rock reactions and intracrustal differentiation on magma chemistry across mid-ocean ridges, intraplate settings and subduction zones using experimental studies and natural data. We compare melting behaviors of pyroxenites and peridotites and their relative contributions to magmas as functions of composition, mantle potential temperatures and lithospheric thickness. We also discuss the fate of chemically distinct melts derived from heterogeneities as they travel through a peridotitic mantle. Using nearly 60,000 natural major element compositions of volcanic rocks, melt inclusions, and crystalline cumulates, we assess broad petrogenetic trends in as large of a global dataset as possible. Consistent with previous studies, major element chemistry of mid-ocean ridge basalts (MORBs) and their cumulates favor a first-order control of intracrustal crystal-liquid segregation, while trace element studies emphasize the role of melt-rock reactions, highlighting the decoupling between the two. Ocean island basalts (OIB) show a larger compositional variability than MORB, partly attributed to large variations of pyroxenite proportions in the mantle source. However, the estimated proportions vary considerably with heterogeneity composition, melting model and thermal structure of the mantle. For arcs, we highlight current views on the role of the downgoing slab into the source of primary arc magmas, and the role of the overriding lithosphere as a magmatic chemical filter and as the repository of voluminous arc cumulates. Our approach of simultaneously looking at a large database of volcanic + deep crustal rocks across diverse tectonic settings underscores the challenge of deciphering the source signal versus intracrustal/lithospheric processes.


title: Excepional Points from Hamiltonians of hybrid physical systems: Squeezing and anti-Squeezing

abstract: We study the appearance of Exceptional Points in a hybrid system composed of a superconducting flux-qubit and an ensemble of nitrogen-vacancy colour centres in diamond. We discuss the possibility of controlling the generation of Exceptional Points, by the analysis of the model space parameters. One of the characteristic features of the presence of Exceptional Points, it is the departure from the exponential decay behaviour of the observables as a function of time. We study the time evolution of different initial states, in the presence of the hybrid system, by computing the reduced density matrix of each subsystem. We present the results we have obtained for the steady behaviour of different observables. We analyse the appearance of Squeezed Spin States and of anti-Squeezed Spin States.


title: Observational evidence for the origin of high-energy neutrinos in parsec-scale nuclei of radio-bright active galaxies

abstract: Observational information on high-energy astrophysical neutrinos is being continuously collected by the IceCube observatory. However, the sources of neutrinos are still unknown. In this study, we use radio very-long-baseline interferometry (VLBI) data for a complete VLBI-flux-density limited sample of active galactic nuclei (AGN). We address the problem of the origin of astrophysical neutrinos with energies above 200 TeV in a statistical manner. It is found that AGN positionally associated with IceCube events have typically stronger parsec-scale cores than the rest of the sample. The post-trial probability of a chance coincidence is 0.2%. We select the four strongest AGN as highly probable associations: 3C 279, NRAO 530, PKS 1741-038, and PKS 2145+067. Moreover, we find an increase of radio emission at frequencies above 10 GHz around neutrino arrival times for several other VLBI-selected AGN on the basis of RATAN-600 monitoring. The most pronounced example of such behavior is PKS 1502+106. We conclude that AGN with bright Doppler-boosted jets constitute an important population of neutrino sources. High-energy neutrinos are produced in their central parsec-scale regions, probably in proton-photon interactions at or around the accretion disk. Radio-bright AGN that are likely associated with neutrinos have very diverse gamma-ray properties suggesting that gamma-rays and neutrinos may be produced in different regions of AGN and not directly related. A small viewing angle of the jet-disk axis is, however, required to detect either of them.


title: Asymmetric Drift of Andromeda Analogs in the IllustrisTNG Simulation

abstract: We analyze the kinematics as a function of stellar age for Andromeda (M31) mass analogs from the IllustrisTNG cosmological simulation. We divide the star particles into four age groups: less than 1 Gyr, 1 to 5 Gyr, 5 to 10 Gyr, and greater 10 Gyr, and compare the kinematics of these groups to that of the neutral gas cells. We calculate rotation curves for the stellar and gaseous components of each analog from 2 kpc to 20 kpc from the center of mass. We find that the lag, or asymmetric drift (AD), between the gas rotation curve and the stellar rotation curve on average increases with stellar age. This finding is consistent with observational measurements of AD in the disk of the Andromeda galaxy. When the M31 analogs are separated into groups based on merger history, we find that there is a difference in the AD of the analogs that have had a 4:1 merger the last 4 Gyr, 8 Gyr, or 12 Gyr compared to analogs that have not experienced a 4:1 merger in the same time frame. The subset of analogs that have had a 4:1 merger within the last 4 Gyr are also similar to AD measurements of stars in the disk of M31, providing evidence that M31 may in fact have recently merged with a galaxy nearly one fourth of its mass. Further work using high resolution zoom in simulations is required to explore the contribution of internal heating to AD.


title: The analogy of K-correction in the topic of gamma-ray bursts

abstract: It is well-known that there are two types of gamma-ray bursts (GRBs): short/hard and long/soft ones, respectively. The long GRBs are coupled to supernovae, but the short ones are associated with the so called macronovae (also known as kilonovae), which can serve as the sources of gravitational waves as well. The kilonovae can arise from the merging of two neutron-stars. The neutron stars can be substituded by more massive black holes as well. Hence, the topic of gamma-ray bursts (mainly the topic of short ones) and the topic of massive binaries, are strongly connected. In this contribution, the redshifts of GRBs are studied. The surprising result - namely that the apparently fainter GRBs can be in average at smaller distances - is discussed again. In essence, the results of M\'esz\'aros et al. (2011) are studied again using newer samples of GRBs. The former result is confirmed by the newer data.


title: Gamma-ray bursts: A brief survey of the diversity

abstract: The separation of the gamma-ray bursts (GRBs) into short/hard and long/soft subclasses, respectively, is well supported both theoretically and observationally. The long ones are coupled to supernovae type Ib/Ic - the short ones are connected to the merging of two neutron stars, where one or even both neutron stars can be substituted by black holes. These short GRBs - as merging binaries - can also serve as the sources of gravitation waves, and are observable as the recently detected macronovae. Since 1998 there are several statistical studies suggesting the existence of more than two subgroups. There can be a subgroup having an intermediate durations; there can be a subgroup with ultra-long durations; the long/soft subgroup itself can be divided into two subclasses with respect to the luminosity of GRBs. The authors with other collaborators provided several statistical studies in this topic. This field of the GRB-diversity is briefly surveyed in this contribution.


title: Pedagogical comments about nonperturbative Ward-constrained melonic renormalization group flow

abstract: This paper, in addition to our recent works, intends to explore the behavior of the Wetterich flow equations in the portion of the theory space spanned by non-branching melons constrained with Ward-identities. We focus on a rank-5 just-renormalizable tensorial group field theory and consider a non-trivial extension of the local potential approximation namely effective vertex expansion for just-renormalizable quartic melonic interactions, disregarding effects coming from disconnected interactions. Investigating the dynamically constrained flow, we show explicitly that results weakly rely on the number of quartic interactions involved in the classical action. In particular, the predictions for the fully connected model are essentially the same as for the single colored model. Finally, closing the flow equations using Ward identities without additional assumptions to compute integrals involved in the effective vertex expansion, we do not find reliable fixed point in the unconstrained theory space connected with the Gaussian region.


title: Gravitational Theory of Cosmology, Galaxies and Galaxy Clusters

abstract: A modified gravitational theory explains early universe and late time cosmology, galaxy and galaxy cluster dynamics. The modified gravity (MOG) theory extends general relativity (GR) by three extra degrees of freedom: a scalar field $G$, enhancing the strength of the Newtonian gravitational constant $G_N$, a gravitational, spin 1 vector graviton field $\phi_\mu$, and the effective mass $\mu$ of the ultralight spin 1 graviton. For $t < t_{\rm rec}$, where $t_{\rm rec}$ denotes the time of recombination and re-ionization, the density of the vector graviton $\rho_\phi > \rho_b$, where $\rho_b$ is the density of baryons, while for $t > t_{\rm rec}$ we have $\rho_b > \rho_\phi$. The matter density is parameterized by $\Omega_M=\Omega_b+\Omega_\phi+\Omega_r$ where $\Omega_r=\Omega_\gamma+\Omega_\nu$. For the cosmological parameter values obtained by the Planck Collaboration, the CMB acoustical oscillation power spectrum, polarization and lensing data can be fitted as in the $\Lambda$CDM model. When the baryon density $\rho_b$ dominates the late time universe, MOG explains galaxy rotation curves, the dynamics of galaxy clusters, galaxy lensing and the galaxy clusters matter power spectrum without dominant dark matter.


title: Binary Kripke Semantics for a Strong Logic for Naive Truth

abstract: I show that the logic $\textsf{TJK}^{d+}$, one of the strongest logics currently known to support the naive theory of truth, is obtained from the Kripke semantics for constant domain intuitionistic logic by (i) dropping the requirement that the accessibility relation is reflexive and (ii) only allowing reflexive worlds to serve as counterexamples to logical consequence. In addition, I provide a simplified natural deduction system for $\textsf{TJK}^{d+}$, in which a restricted form of conditional proof is used to establish conditionals.


title: Resilient Multi-Dimensional Consensus in Adversarial Environment

abstract: This paper considers the multi-dimensional consensus in networked systems, where some of the agents might be misbehaving (or faulty). Despite the influence of these misbehaviors, the benign agents aim to reach an agreement while avoiding being seriously influenced by the faulty ones. To this end, this paper first considers a general class of consensus algorithms, where each benign agent computes an "auxiliary point" based on the received values and moves its state toward this point. Concerning this generic form, we present conditions for achieving resilient consensus and obtain a lower bound on the exponential convergence rate. Assuming that the number of malicious agents is upper bounded, two specific resilient consensus algorithms are further developed based on the obtained conditions. Particularly, the first solution, based on Helly's Theorem, achieves the consensus within the convex hull formed by the benign agents' initial states, where the auxiliary point can be efficiently computed through linear programming. On the other hand, the second algorithm serves as a "built-in" security guarantee for standard average consensus algorithms, in the sense that its performance coincides exactly with that of the standard ones in the absence of faulty nodes while also resisting the serious influence of the misbehaving ones in adversarial environment. Some numerical examples are provided in the end to verify the theoretical results.


title: Torsion Discriminance for Stability of Linear Time-Invariant Systems

abstract: This paper proposes a new approach to describe the stability of linear time-invariant systems via the torsion $\tau(t)$ of the state trajectory. For a system $\dot{r}(t)=Ar(t)$ where $A$ is invertible, we show that (1) if there exists a measurable set $E_1$ with positive Lebesgue measure, such that $r(0)\in E_1$ implies that $\lim\limits_{t\to+\infty}\tau(t)\neq0$ or $\lim\limits_{t\to+\infty}\tau(t)$ does not exist, then the zero solution of the system is stable; (2) if there exists a measurable set $E_2$ with positive Lebesgue measure, such that $r(0)\in E_2$ implies that $\lim\limits_{t\to+\infty}\tau(t)=+\infty$, then the zero solution of the system is asymptotically stable. Furthermore, we establish a relationship between the $i$th curvature $(i=1,2,\cdots)$ of the trajectory and the stability of the zero solution when $A$ is similar to a real diagonal matrix.


title: Relative Flatness and Generalization

abstract: Flatness of the loss curve is conjectured to be connected to the generalization ability of machine learning models, in particular neural networks. While it has been empirically observed that flatness measures consistently correlate strongly with generalization, it is still an open theoretical problem why and under which circumstances flatness is connected to generalization, in particular in light of reparameterizations that change certain flatness measures but leave generalization unchanged. We investigate the connection between flatness and generalization by relating it to the interpolation from representative data, deriving notions of representativeness, and feature robustness. The notions allow us to rigorously connect flatness and generalization and to identify conditions under which the connection holds. Moreover, they give rise to a novel, but natural relative flatness measure that correlates strongly with generalization, simplifies to ridge regression for ordinary least squares, and solves the reparameterization issue.


title: Numerical modeling of thin anisotropic composite membrane under dynamic load

abstract: This work aims to describe a mathematical model and a numerical method to simulate a thin anisotropic composite membrane moving and deforming in 3D space under a dynamic load of an arbitrary time and space profile. The model and the method allow to consider problems when quasi static approximation is not valid, and elastic waves caused by the impact should be calculated. The model and the method can be used for numerical study of different processes in thin composite layers, such as shock load, ultrasound propagation, non-destructive testing procedures, vibrations. The thin membrane is considered as 2D object in 3D space, this approach allows to reduce computational time still having an arbitrary material rheology and load profile.


title: The LSST DESC Data Challenge 1: Generation and Analysis of Synthetic Images for Next Generation Surveys

abstract: Data Challenge 1 (DC1) is the first synthetic dataset produced by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). DC1 is designed to develop and validate data reduction and analysis and to study the impact of systematic effects that will affect the LSST dataset. DC1 is comprised of $r$-band observations of 40 deg$^{2}$ to 10-year LSST depth. We present each stage of the simulation and analysis process: a) generation, by synthesizing sources from cosmological N-body simulations in individual sensor-visit images with different observing conditions; b) reduction using a development version of the LSST Science Pipelines; and c) matching to the input cosmological catalog for validation and testing. We verify that testable LSST requirements pass within the fidelity of DC1. We establish a selection procedure that produces a sufficiently clean extragalactic sample for clustering analyses and we discuss residual sample contamination, including contributions from inefficiency in star-galaxy separation and imperfect deblending. We compute the galaxy power spectrum on the simulated field and conclude that: i) survey properties have an impact of 50\% of the statistical uncertainty for the scales and models used in DC1 ii) a selection to eliminate artifacts in the catalogs is necessary to avoid biases in the measured clustering; iii) the presence of bright objects has a significant impact (2- to 6-$\sigma$) in the estimated power spectra at small scales ($\ell > 1200$), highlighting the impact of blending in studies at small angular scales in LSST;


title: Simple explanation of Landauer's bound and its ineffectiveness for multivalued logic

abstract: We discuss, using recent results on the Landauer's bound in multivalued logic, the difficulties and pitfalls of how to apply this principle. The presentation is based on Szilard's version of Maxwell's demon experiment and use of equilibrium Thermodynamics. Different versions of thermodynamical/mechanical memory are presented - one-hot encoding version and the implementation based on reversed Szilard's experiment. Relation of the Landauer's principle to Galois connection is explained in detail.

categories: cond-mat.stat-mech cond-mat.mtrl-sci cs.IT math-ph math.IT math.MP


title: On the Finite Optimal Convergence of Logic-Based Benders' Decomposition in Solving 0-1 Min-max Regret Optimization Problems with Interval Costs

abstract: This paper addresses a class of problems under interval data uncertainty composed of min-max regret versions of classical 0-1 optimization problems with interval costs. We refer to them as interval 0-1 min-max regret problems. The state-of-the-art exact algorithms for this class of problems work by solving a corresponding mixed integer linear programming formulation in a Benders' decomposition fashion. Each of the possibly exponentially many Benders' cuts is separated on the fly through the resolution of an instance of the classical 0-1 optimization problem counterpart. Since these separation subproblems may be NP-hard, not all of them can be modeled by means of linear programming, unless P = NP. In these cases, the convergence of the aforementioned algorithms are not guaranteed in a straightforward manner. In fact, to the best of our knowledge, their finite convergence has not been explicitly proved for any interval 0-1 min-max regret problem. In this work, we formally describe these algorithms through the definition of a logic-based Benders' decomposition framework and prove their convergence to an optimal solution in a finite number of iterations. As this framework is applicable to any interval 0-1 min-max regret problem, its finite optimal convergence also holds in the cases where the separation subproblems are NP-hard.


title: On the definition of a concentration function relevant to the ROC curve

abstract: This is a reader's reaction to a recent paper by E. Schechtman and G. Schechtman (Metron, 2019) about the correct definition of a concentration function for the diagnostic, i.e. supervised classification, problem. We propose and motivate a different definition and refer to the relevant literature.


title: Equivalence of the Fermat potential and the lensing potential approaches to computing the integrated Sachs-Wolfe effect

abstract: We show in detail that the recently derived expression for evaluating the integrated Sachs-Wolfe (ISW) temperature shift in the cosmic microwave background (CMB) caused by individual embedded (compensated) lenses is equivalent to the conventional approach for flat background cosmologies. The conventional approach requires evaluating an integral of the time derivative of the lensing potential, whereas the new Fermat potential approach is simpler and only requires taking a derivative of the potential part of the time delay.


title: Double-End Queues with Non-Poisson Inputs and Their Effective Algorithms

abstract: It is interesting and challenging to study double-ended queues with First-Come-First-Match discipline under customers' impatient behavior and non-Poisson inputs. The system stability can be guaranteed by the customers' impatient behavior, while the existence of impatient customers makes analysis of such double-ended queues more difficult or even impossible to find an explicitly analytic solution, thus it becomes more and more important to develop effective numerical methods in a variety of practical matching problems. This paper studies a block-structured double-ended queue, whose block structure comes from two independent Markovian arrival processes (MAPs), which are non-Poisson inputs. We show that such a queue can be expressed as a new bilateral quasi birth-and-death (QBD) process which has its own interest. Based on this, we provide a detailed analysis for both the bilateral QBD process and the double-ended queue, including the system stability, the queue size distributions, the average stationary queue lengths, and the sojourn time of any arriving customers. Furthermore, we develop three effective algorithms for computing the performance measures (i.e., the probabilities of stationary queue lengths, the average stationary queue lengths, and the average sojourn times) of the double-ended queue with non-Poisson inputs. Finally, we use some numerical examples in tabular and graphical to illustrate how the performance measures are influenced by some key system parameters. We believe that the methodology and results described in this paper can be applicable to deal with more general double-ended queues in practice, and develop some effective algorithms for the purpose of many actual uses.


title: Scanning Transmission Electron Microscopy of Oxide Interfaces and Heterostructures

abstract: Thin film oxides are a source of endless fascination for the materials scientist. These materials are highly flexible, can be integrated into almost limitless combinations, and exhibit many useful functionalities for device applications. While precision synthesis techniques, such as molecular beam epitaxy (MBE) and pulsed laser deposition (PLD), provide a high degree of control over these systems, there remains a disconnect between ideal and realized materials. Because thin films adopt structures and chemistries distinct from their bulk counterparts, it is often difficult to predict what properties will emerge. The complex energy landscape of the synthesis process is also strongly influenced by non-equilibrium growth conditions imposed by the substrate, as well as the kinetics of thin film crystallization and fluctuations in process variables, all of which can lead to significant deviations from targeted outcomes. High-resolution structural and chemical characterization techniques, as described in this volume, are needed to verify growth models, bound theoretical calculations, and guide materials design. While many characterization options exist, most are spatially-averaged or indirect, providing only partial insight into the complex behavior of these systems. Over the past several decades, scanning transmission electron microscopy (STEM) has become a cornerstone of oxide heterostructure characterization owing to its ability to simultaneously resolve structure, chemistry, and defects at the highest spatial resolution. STEM methods are an essential complement to averaged scattering techniques, offering a direct picture of resulting materials that can inform and refine the growth process to achieve targeted properties. There is arguably no other technique that can provide such a broad array of information at the atomic-scale, all within a single experimental session.


title: Speeding Up Dark Matter With Solar Neutrinos

abstract: We present a novel mechanism of using solar neutrinos to speed up dark matter, inspired by the fact that neutrinos are the most energetic particles from the Sun with a well-understood spectrum. In a neutrino portal dark sector model, we show that dark matter with sub-GeV mass could be accelerated by the $pp$ neutrinos to velocities well above $10^{-3}c$ and capable of depositing large enough energy at direct detection experiments. A crucial ingredient of this mechanism is the dissociation of stable dark matter bound states that exist in Nature. The resulting dark matter velocity distribution bears a strong resemblance in shape to the solar neutrino spectrum. As an application, we derive a leading limit on light dark matter interaction by reinterpreting a recent PICO experiment result.