title: Hyperspectral and multispectral image fusion under spectrally varying spatial blurs -- Application to high dimensional infrared astronomical imaging

abstract: Hyperspectral imaging has become a significant source of valuable data for astronomers over the past decades. Current instrumental and observing time constraints allow direct acquisition of multispectral images, with high spatial but low spectral resolution, and hyperspectral images, with low spatial but high spectral resolution. To enhance scientific interpretation of the data, we propose a data fusion method which combines the benefits of each image to recover a high spatio-spectral resolution datacube. The proposed inverse problem accounts for the specificities of astronomical instruments, such as spectrally variant blurs. We provide a fast implementation by solving the problem in the frequency domain and in a low-dimensional subspace to efficiently handle the convolution operators as well as the high dimensionality of the data. We conduct experiments on a realistic synthetic dataset of simulated observation of the upcoming James Webb Space Telescope, and we show that our fusion algorithm outperforms state-of-the-art methods commonly used in remote sensing for Earth observation.


title: A one-dimensional non-local singular SPDE

abstract: We examine in this article the one-dimensional, non-local, singular SPDE \begin{equation*} \partial_t u \;=\; -\, (-\Delta)^{1/2} u \,-\, \sinh(\gamma u) \,+\, \xi\;, \end{equation*} where $\gamma\in \mathbb{R}$, $(-\Delta)^{1/2}$ is the fractional Laplacian of order $1/2$, $\xi$ the space-time white noise in $\mathbb{R} \times \mathbb{T}$, and $\mathbb{T}$ the one-dimensional torus. We show that for $0<\gamma^2<\pi/7$ the Da Prato--Debussche method applies. One of the main difficulties lies in the derivation of a Schauder estimate for the semi-group associated to the fractional Laplacian due to the lack of smoothness resulting from the long range interaction.


title: Positively and negatively hydrated counterions in molecular dynamics simulations of DNA double helix

abstract: The DNA double helix is a polyanionic macromolecule that in water solutions is neutralized by metal ions (counterions). The property of the counterions to stabilize the water network (positive hydration) or to make it friable (negative hydration) is important in terms of the physical mechanisms of stabilization of the DNA double helix. In the present research, the effects of positive hydration of Na$^{+}$ counterions and negative hydration of K$^{+}$ and Cs$^{+}$ counterions, incorporated into the hydration shell of the DNA double helix have been studied using molecular dynamics simulations. The results have shown that the dynamics of the hydration shell of counterions depends on region of the double helix: minor groove, major groove, and outside the macromolecule. The longest average residence time has been observed for water molecules contacting with the counterions, localized in the minor groove of the double helix (about 50 ps for Na$^{+}$, and lower than 10 ps for K$^{+}$ and Cs$^{+}$). The estimated potentials of mean force for the hydration shells of the counterions show that the water molecules are constrained too strong, and, consequently, the effect of negative hydration for K$^{+}$ and Cs$^{+}$ counterions has not been observed in the simulations. The analysis has shown that the effects of counterion hydration can be described more accurately with water models having lower dipole moments.


title: Circular and Linear ${\rm e^+ e^-}$ Colliders: Another Story of Complementarity

abstract: The remarkable synergy and complementarity between the circular ${\rm e^+ e^-}$ and pp colliders has been extensively discussed. In this short document, we investigate the complementarity between the proposed circular and linear ${\rm e^+ e^-}$ colliders at the electroweak and TeV scale. This complementarity could be exploited on a world-wide scale, if both a large circular and a linear infrastructures were available. A possible implementation of such a complementary program is shown.


title: Vision and Language: from Visual Perception to Content Creation

abstract: Vision and language are two fundamental capabilities of human intelligence. Humans routinely perform tasks through the interactions between vision and language, supporting the uniquely human capacity to talk about what they see or hallucinate a picture on a natural-language description. The valid question of how language interacts with vision motivates us researchers to expand the horizons of computer vision area. In particular, "vision to language" is probably one of the most popular topics in the past five years, with a significant growth in both volume of publications and extensive applications, e.g., captioning, visual question answering, visual dialog, language navigation, etc. Such tasks boost visual perception with more comprehensive understanding and diverse linguistic representations. Going beyond the progresses made in "vision to language," language can also contribute to vision understanding and offer new possibilities of visual content creation, i.e., "language to vision." The process performs as a prism through which to create visual content conditioning on the language inputs. This paper reviews the recent advances along these two dimensions: "vision to language" and "language to vision." More concretely, the former mainly focuses on the development of image/video captioning, as well as typical encoder-decoder structures and benchmarks, while the latter summarizes the technologies of visual content creation. The real-world deployment or services of vision and language are elaborated as well.


title: Scaling relations and finite-size scaling in gravitationally correlated lattice percolation models

abstract: In some systems, the connecting probability (and thus the percolation process) between two sites depends on the geometric distance between them. To understand such process, we propose gravitationally correlated percolation models for link-adding networks on the two-dimensional lattice $G$ with two strategies $S_{\rm max}$ and $S_{\rm min}$, to add a link $l_{i,j}$ to connect site $i$ and site $j$ with mass $m_i$ and $m_j$, respectively; $m_i$ and $m_j$ are sizes of the clusters which contain site $i$ and site $j$, respectively. The probability to add the link $l_{i,j}$ is related to the generalized gravity $g_{ij} \equiv m_i m_j/r_{ij}^d$, where $r_{ij}$ is the geometric distance between $i$ and $j$, and $d$ is an adjustable decaying exponent. In the beginning of the simulation, all sites of $G$ are occupied and there is no link. In the simulation process, two inter-cluster links $l_{i,j}$ and $l_{k,n}$ are randomly chosen and the generalized gravities $g_{ij}$ and $g_{kn}$ are computed. In the strategy $S_{\rm max}$, the link with larger generalized gravity is added. In the strategy $S_{\rm min}$, the link with smaller generalized gravity is added, which include percolation on the Erd\H os-R\'enyi random graph and the Achlioptas process of explosive percolation as the limiting cases, $d \to \infty$ and $d \to 0$, respectively. Adjustable strategies facilitate or inhibit the network percolation in a generic view. We calculate percolation thresholds $T_c$ and critical exponents $\beta$ by numerical simulations. We also obtain various finite-size scaling functions for the node fractions in percolating clusters or arrival of saturation length with different intervening strategies.


title: Studying the localized $CP$ violation and the branching fraction of the $\bar{B}^0\rightarrow K^-\pi^+\pi^+\pi^-$ decay

abstract: In this work, we study the localized $CP$ violation and the branching fraction of the four-body decay $\bar{B}^0\rightarrow K^-\pi^+\pi^-\pi^+$ by employing a quasi-two-body QCD factorization approach. Considering the interference of $\bar{B}^0\rightarrow \bar{K}_0^*(700)\rho^0(770)\rightarrow K^-\pi^+\pi^-\pi^+$ and $\bar{B}^0\rightarrow \bar{K}^*(892)f_0(500)\rightarrow K^-\pi^+\pi^-\pi^+$ channels, we predict $\mathcal{A_{CP}}(\bar{B}^0\rightarrow K^-\pi^+\pi^-\pi^+)\in[0.15,0.28]$ and $\mathcal{B}(\bar{B}^0\rightarrow K^-\pi^+\pi^-\pi^+)\in[1.73,5.10]\times10^{-7}$, respectively, which shows that this two channels' interference mechanism can induce the localized $CP$ violation to this four-body decay. Meanwhile, within the two quark model framework for the scalar mesons $f_0(500)$ and $\bar{K}_0^*(700)$, we calculate the direct CP violations and branching fractions of the $\bar{B}^0\rightarrow \bar{K}_0^*(700)\rho^0(770)$ and $\bar{B}^0\rightarrow \bar{K}^*(892)f_0(500)$ decays, respectively. The corresponding results are $\mathcal{A_{CP}}(\bar{B}^0\rightarrow \bar{K}_0^*(700)\rho^0(770)) \in [0.20, 0.36]$, $\mathcal{A_{CP}}(\bar{B}^0\rightarrow \bar{K}^*(892)f_0(500))\in [0.08, 0.12]$, $\mathcal{B} (\bar{B}^0\rightarrow \bar{K}_0^*(700)\rho^0(770)\in [6.76, 18.93]\times10^{-8}$ and $\mathcal{B} (\bar{B}^0\rightarrow \bar{K}^*(892)f_0(500))\in [2.66, 4.80]\times10^{-6}$, respectively, indicating the $CP$ violations of these two two-body decays are both positive and the branching fractions are quite different. These studies provide a new way to investigate the aforementioned four-body decay and could be helpful in clarifying the configuration of the structure of light scalar meson.


title: Future Prospects for the Minimal Supersymmetric Standard Model

abstract: Experimental collaborations for the large hadron collider conducted many and various searches for supersymmetry. In the absence of signals, lower limits were put on sparticle masses but usually within frameworks with (over-) simplifications relative to the entire indications by supersymmetry models. For complementing current interpretations of experimental bounds, we introduce a 30-parameters version of the R-parity conserving minimal supersymmetric standard model (MSSM-30). Using a sample of the MSSM-30 which are in harmony with cold dark matter, flavour and precision electroweak constraints, we explicitly show the prospects for assessing neutralino candidate dark matter in contrast to future searches for supersymmetry. The MSSM-30 parameters regions that are beyond reach to dark matter direct detection experiments could be probed by future hadron-hadron colliders.


title: Dense Matter and Neutron Stars: Some Basic Notions

abstract: A number of properties of dense matter can be understood semiquantitatively in terms of simple physical arguments. We begin with the outer parts of neutron stars, and consider the density at which pressure ionization occurs, the density at which electrons become relativistic, the density at which neutrons drip out of nuclei, and the size of the equilibrium nucleus in dense matter. Subsequently, we treat the so-called "pasta" phases expected to occur at densities just below the density at which the transition from the crust to the liquid core of a neutron star occurs. We then consider aspects of superfluidity in dense matter. Estimates of pairing gaps in homogeneous nuclear matter are given, and the effect of the dense medium on the interaction between nucleons is described. Finally, we turn to superfluidity in the crust of neutron stars and especially the neutron superfluid density, an important quantity in the theory of sudden speedups of the rotation rate of some pulsars.


title: Double-ionization mechanisms of magnesium driven by electron impact

abstract: We study double ionization of Mg by electron impact through the vantage point of classical mechanics. We consider all electron-electron correlations in a Coulomb four-body problem, where three electrons belong to the atom and the fourth electron causes the impact ionization. From our model we compute the double-ionization probability of Mg for impact energies from 15, to 125 eV. Double ionization occurs through eight double-ionization mechanisms, which we classify into four categories: inner shell capture, direct, delay and ionized inner shell mechanisms. We show that delay and ionized inner shell mechanisms require electron-electron correlations among the four electrons, and are responsible for the second increase in the double-ionization probability. Furthermore, we show that our theoretical prediction about the relative prominence of certain double ionization mechanisms is in agreement with experimental results on the relative prominence of non-first- over first-order mechanisms.


title: Band gap opening at the Dirac point in Co/BiSbTeSe2(0001) system

abstract: Sub-angstrom Co coverage, being deposited on BiSbTeSe2(0001) surface at 200-330 C, opens a band gap at the Dirac point, with the shift of the Dirac point position caused by RT adsorbate pre-deposition. Temperature dependent measurements in 15-150 K range have shown no band gap width change. This fact indicates the nonmagnetic nature of the gap which may be attributed to the chemical hybridization of surface states upon the introduction of Co adatoms, which decrease crystallographic symmetry and eliminate topological protection of the surface states.


title: Exploring the Latest Pantheon SNIa Dataset by Using Three Kinds of Statistics Techniques

abstract: In this work, we explore the cosmological consequences of the latest Type Ia supernova (SN Ia) data-set, Pantheon, by adopting the $wCDM$ model. The Pantheon data-set is the largest SN Ia samples till now, which contains 1048 supernovae on the redshift range $0 < z < 2.3$. Here we take into account three kinds of SN Ia statistics techniques, including: 1. magnitude statistics (MS), which is the traditional SN Ia statistics technique; 2. flux statistics (FS), which bases on the flux-averaging (FA) method; 3. improved flux statistics (IFS), which combines the advantages of MS and FS. It should be mentioned that, The IFS technique need to scan the $(z_{cut},\Delta z)$ parameters plane, where $z_{cut}$ and $\Delta z$ are redshift cut-off and redshift interval of FA, respectively. The results are shown as follows. (1) Using SN data-set only, the best FA recipe for IFS is $(z_{cut},\Delta z)=(0.1,0.08)$; (2) Comparing to the old SN data-set, JLA, adopting Pantheon data-set can reduce the $2\sigma$ error bars of equation of state $w$ by 38\%, 47\% and 53\% for MS, FS and IFS, respectively; (3) FS gives closer results to other observations, such as Baryon acoustic oscillations and Cosmic microwave background; (4) Compared with FS and IFS, MS more favors a Universe that will end in a "big rip".


title: Necessary Conditions for Adverse Control Problems Expressed by Relaxed Derivatives

abstract: This paper provides a framework for deriving a new set of necessary conditions for adverse control problems among two players. The distinguish feature of such problems is that the first player has a priori knowledge on the second player strategy. A subclass of adverse control problems is the one of minimax control problems, which frequently arise in robust dynamic optimization. The conditions derived in this manuscript are expressed in terms of relaxed derivatives: the dual variables and the related functions are limits of computable sequences, obtained by considering a regularized version of the original problem and applying well known necessary condition. This topic was initially treated by J. Warga.


title: Arbitrary amplitude ion acoustic solitary structures in a collisionless magnetized plasma consisting of nonthermal and isothermal electrons

abstract: We have used the Sagdeev pseudo-potential technique to investigate the arbitrary amplitude ion acoustic solitons, double layers and supersolitons in a collisionless magnetized plasma consisting of adiabatic warm ions, isothermal cold electrons and nonthermal hot electrons immersed in an external uniform static magnetic field. We have used the phase portraits of the dynamical system describing the nonlinear behaviour of ion acoustic waves to confirm the existence of different solitary structures. We have also investigated the transition of different solitary structures: soliton (before the formation of double layer) $\rightarrow$ double layer $\rightarrow$ supersoliton $\rightarrow$ soliton (soliton after the formation of double layer) by considering the variation of $\theta$ only, where $\theta$ is the angle between the direction of the external uniform static magnetic field and the direction of propagation of the wave.


title: What can a heavy $U(1)_{\rm B-L}$ $Z^\prime$ boson do to the muon $(g-2)_\mu$ anomaly and to a new Higgs boson mass?

abstract: The minimal $U(1)_{\rm B-L}$ extension of the Standard Model (B-L-SM) offers an explanation for neutrino mass generation via a seesaw mechanism as well as contains two new physics states such as an extra Higgs boson and a new $Z^\prime$ gauge boson. The emergence of a second Higgs particle as well as a new $Z^\prime$ gauge boson, both linked to the breaking of a local $U(1)_{\rm B-L}$ symmetry, makes the B-L-SM rather constrained by direct searches at the Large Hadron Collider (LHC) experiments. We investigate the phenomenological status of the B-L-SM by confronting the new physics predictions with the LHC and electroweak precision data. Taking into account the current bounds from direct LHC searches, we demonstrate that the prediction for the muon $(g-2)_\mu$ anomaly in the B-L-SM yields at most a contribution of approximately $8.9 \times 10^{-12}$ which represents a tension of $3.28$ standard deviations, with the current $1\sigma$ uncertainty, by means of a $Z^\prime$ boson if its mass lies in a range of $6.3$ to $6.5$ TeV, within the reach of future LHC runs. This means that the B-L-SM, with heavy yet allowed $Z^\prime$ boson mass range, in practice does not resolve the tension between the observed anomaly in the muon $(g-2)_\mu$ and the theoretical prediction in the Standard Model. Such a heavy $Z^\prime$ boson also implies that the minimal value for a new Higgs mass is of the order of 400 GeV.


title: Mass constraints for 15 protoplanetary disks from HD 1-0

abstract: Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $\Delta_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time.


title: Heat flow and noncommutative quantum mechanics in phase-space

abstract: In this work, we investigate the heat flow of two interacting quantum systems on the perspective of noncommutativity phase-space effects and show that by controlling the new constants introduced in the quantum theory, due to a deformed Heisenberg-Weyl algebra, the heat flow from the hot to the cold system may be enhanced, thus decreasing the time required to reach thermal equilibrium. We also give a brief discussion on the robustness of the second law of thermodynamics in the context of noncommutative quantum mechanics


title: Giant Valley-Zeeman Splitting from Spin-Singlet and Spin-Triplet Interlayer Excitons in WSe2/MoSe2 Heterostructure

abstract: Transition metal dichalcogenides (TMDCs) heterostructure with a type II alignment hosts unique interlayer excitons with the possibility of spin-triplet and spin-singlet states. However, the associated spectroscopy signatures remain elusive, strongly hindering the understanding of the Moire potential modulation of the interlayer exciton. In this work, we unambiguously identify the spin-singlet and spin-triplet interlayer excitons in the WSe2/MoSe2 hetero-bilayer with a 60-degree twist angle through the gate- and magnetic field-dependent photoluminescence spectroscopy. Both the singlet and triplet interlayer excitons show giant valley-Zeeman splitting between the K and K' valleys, a result of the large Lande g-factor of the singlet interlayer exciton and triplet interlayer exciton, which are experimentally determined to be ~ 10.7 and ~ 15.2, respectively, in good agreement with theoretical expectation. The PL from the singlet and triplet interlayer excitons show opposite helicities, determined by the atomic registry. Helicity-resolved photoluminescence excitation (PLE) spectroscopy study shows that both singlet and triplet interlayer excitons are highly valley-polarized at the resonant excitation, with the valley polarization of the singlet interlayer exciton approaches unity at ~ 20 K. The highly valley-polarized singlet and triplet interlayer excitons with giant valley-Zeeman splitting inspire future applications in spintronics and valleytronics.


title: Microtoroidal resonators enhance long-distance dynamical entanglement generation in chiral quantum networks

abstract: Chiral quantum networks provide a promising route for realising quantum information processing and quantum communication. Here, we describe how two distant quantum nodes of chiral quantum network become dynamically entangled by a photon transfer through a common 1D chiral waveguide. We harness the directional asymmetry in chirally-coupled single-mode ring resonators to generate entangled state between two atoms. We report a concurrence of up to 0.969, a huge improvement over the 0.736 which was suggested and analyzed in great detail in Ref. [1]. This significant enhancement is achieved by introducing microtoroidal resonators which serve as efficient photonic interface between light and matter. Robustness of our protocol to experimental imperfections such as fluctuations in inter-nodal distance, imperfect chirality, various detunings and atomic spontaneous decay is demonstrated. Our proposal can be utilised for long-distance entanglement generation in quantum networks which is a key ingredient for many applications in quantum computing and quantum information processing.


title: Axion coupling in the hybrid Wannier representation

abstract: Many magnetic point-group symmetries induce a topological classification on crystalline insulators, dividing them into those that have a nonzero quantized Chern-Simons magnetoelectric coupling ("axion-odd" or "topological"), and those that do not ("axion-even" or "trivial"). For time-reversal or inversion symmetries, the resulting topological state is usually denoted as a "strong topological insulator" or an "axion insulator" respectively, but many other symmetries can also protect this "axion $Z_2$" index. Topological states are often insightfully characterized by choosing one crystallographic direction of interest, and inspecting the hybrid Wannier (or equivalently, the non-Abelian Wilson-loop) band structure, considered as a function of the two-dimensional Brillouin zone in the orthogonal directions. Here, we systematically classify the axion-quantizing symmetries, and explore the implications of such symmetries on the Wannier band structure. Conversely, we clarify the conditions under which the axion$Z_2$ index can be deduced from the Wannier band structure. In particular, we identify cases in which a counting of Dirac touchings between Wannier bands, or a calculation of the Chern number of certain Wannier bands, or both, allows for a direct determination of the axion $Z_2$ index. We also discuss when such symmetries impose a "flow" on the Wannier bands, such that they are always glued to higher and lower bands by degeneracies somewhere in the projected Brillouin zone, and the related question of when the corresponding surfaces can remain gapped, thus exhibiting a half-quantized surface anomalous Hall conductivity. Our formal arguments are confirmed and illustrated in the context of tight-binding models for several paradigmatic axion-odd symmetries including time reversal, inversion, simple mirror, and glide mirror symmetries.


title: W-PoseNet: Dense Correspondence Regularized Pixel Pair Pose Regression

abstract: Solving 6D pose estimation is non-trivial to cope with intrinsic appearance and shape variation and severe inter-object occlusion, and is made more challenging in light of extrinsic large illumination changes and low quality of the acquired data under an uncontrolled environment. This paper introduces a novel pose estimation algorithm W-PoseNet, which densely regresses from input data to 6D pose and also 3D coordinates in model space. In other words, local features learned for pose regression in our deep network are regularized by explicitly learning pixel-wise correspondence mapping onto 3D pose-sensitive coordinates as an auxiliary task. Moreover, a sparse pair combination of pixel-wise features and soft voting on pixel-pair pose predictions are designed to improve robustness to inconsistent and sparse local features. Experiment results on the popular YCB-Video and LineMOD benchmarks show that the proposed W-PoseNet consistently achieves superior performance to the state-of-the-art algorithms.


title: Skew braces of size $p^2q$

abstract: In this paper we enumerate the skew braces of size $p^2q$ for $p,q$ odd primes by the classification of regular subgroups of the holomorph of the groups of size $p^2q$. In particular, we provide explicit formulas for the skew braces of abelian type.


title: Cylindrically symmetric 2+1 gravity in terms of global variables:quantum dynamics

abstract: We perform quantization of a model in which gravity is coupled to a circular dust shell in 2+1 spacetime dimensions. Canonical analysis shows that momentum space of this model is ADS^2-space, and the global chart for it is provided by the Euler angles. In quantum kinematics, this results in non-commutativity in coordinate space and discreteness of the shell radius in timelike region, which includes the collapse point. At the level of quantum dynamics, we find transition amplitudes between zero and non-zero eigenvalues of the shell radius, which describe the rate of gravitational collapse (bounce). Their values are everywhere finite, which could be interpreted as resolution of the central singularity.


title: 3DFR: A Swift 3D Feature Reductionist Framework for Scene Independent Change Detection

abstract: In this paper we propose an end-to-end swift 3D feature reductionist framework (3DFR) for scene independent change detection. The 3DFR framework consists of three feature streams: a swift 3D feature reductionist stream (AvFeat), a contemporary feature stream (ConFeat) and a temporal median feature map. These multilateral foreground/background features are further refined through an encoder-decoder network. As a result, the proposed framework not only detects temporal changes but also learns high-level appearance features. Thus, it incorporates the object semantics for effective change detection. Furthermore, the proposed framework is validated through a scene independent evaluation scheme in order to demonstrate the robustness and generalization capability of the network. The performance of the proposed method is evaluated on the benchmark CDnet 2014 dataset. The experimental results show that the proposed 3DFR network outperforms the state-of-the-art approaches.


title: A new approach for the wobbling motion in the even-odd isotopes $^{161,163,165,167}$Lu

abstract: A new interpretation for the wobbling bands in the even-odd Lu isotopes is given within a particle-triaxial rotor semi-classical formalism. While in the previous papers the bands TSD1, TSD2, TSD3 and TSD4 are viewed as the ground, one, two and three phonon wobbling bands, here the corresponding experimental results are described as the ground band with spin equal to I=R+j, for R=0,2,4,...(TSD1), the ground band with I=R+j and R=1,3,5,...(TSD2), the one phonon excitations of TSD2 (TSD3), with the odd proton moving in the orbit $j=i_{13/2}$, and the ground band of I=R+j, with R=1,3,5,... and $j=h_{9/2}$ (TSD4). The moments of inertia (MoI) of the core for the first three bands are the same, and considered to be free parameters. Due to the core polarization effect caused by the particle-core coupling, the MoI's for TSD4 are different. The energies and the e.m. transitions are quantitatively well described. Also, the phase diagram of the odd system is drawn. In the parameter space, one indicates where the points associated with the fitted parameters are located, which is the region where the transversal wobbling mode might be possible, as well as where the wobbling motion is forbidden.


title: Mean Field Games with Branching

abstract: Mean field games are concerned with the limit of large-population stochastic differential games where the agents interact through their empirical distribution. In the classical setting, the number of players is large but fixed throughout the game. However, in various applications, such as population dynamics or economic growth, the number of players can vary across time which may lead to different Nash equilibria. For this reason, we introduce a branching mechanism in the population of agents and obtain a variation on the mean field game problem. As a first step, we study a simple model using a PDE approach to illustrate the main differences with the classical setting. We prove existence of a solution and show that it provides an approximate Nash-equilibrium for large population games. We also present a numerical example for a linear--quadratic model. Then we study the problem in a general setting by a probabilistic approach. It is based upon the relaxed formulation of stochastic control problems which allows us to obtain a general existence result.


title: Analysis of Reference and Citation Copying in Evolving Bibliographic Networks

abstract: Extensive literature demonstrates how the copying of references (links) can lead to the emergence of various structural properties (e.g., power-law degree distribution and bipartite cores) in bibliographic and other similar directed networks. However, it is also well known that the copying process is incapable of mimicking the number of directed triangles in such networks; neither does it have the power to explain the obsolescence of older papers. In this paper, we propose RefOrCite, a new model that allows for copying of both the references from (i.e., out-neighbors of) as well as the citations to (i.e., in-neighbors of) an existing node. In contrast, the standard copying model (CP) only copies references. While retaining its spirit, RefOrCite differs from the Forest Fire (FF) model in ways that makes RefOrCite amenable to mean-field analysis for degree distribution, triangle count, and densification. Empirically, RefOrCite gives the best overall agreement with observed degree distribution, triangle count, diameter, h-index, and the growth of citations to newer papers.


title: Positive Populations

abstract: A positive structure on the varieties of critical points of master functions for KZ equations is introduced. It comes as a combination of the ideas from classical works by G.Lusztig and a previous work by E.Mukhin and the second named author.


title: Efficient Training of Deep Classifiers for Wireless Source Identification using Test SNR Estimates

abstract: We study efficient deep learning training algorithms that process received wireless signals, if a test Signal to Noise Ratio (SNR) estimate is available. We focus on two tasks that facilitate source identification: 1- Identifying the modulation type, 2- Identifying the wireless technology and channel in the 2.4 GHz ISM band. For benchmarking, we rely on recent literature on testing deep learning algorithms against two well-known datasets. We first demonstrate that using training data corresponding only to the test SNR value leads to dramatic reductions in training time while incurring a small loss in average test accuracy, as it improves the accuracy for low SNR values. Further, we show that an erroneous test SNR estimate with a small positive offset is better for training than another having the same error magnitude with a negative offset. Secondly, we introduce a greedy training SNR Boosting algorithm that leads to uniform improvement in accuracy across all tested SNR values, while using a small subset of training SNR values at each test SNR. Finally, we demonstrate the potential of bootstrap aggregating (Bagging) based on training SNR values to improve generalization at low test SNR values with scarcity of training data.


title: Plenoptic x-ray microscopy

abstract: Plenoptic cameras use arrays of micro-lenses to capture multiple views of the same scene in a single compound image. They enable refocusing on different planes and depth estimation. However, until now, all types of plenoptic computational imaging have been limited to visible light. We demonstrate an x-ray plenoptic microscope that uses a concentrating micro-capillary array instead of a micro-lens array and can simultaneously acquire from one hundred to one thousand x-ray projections of imaged volumes that are located in the focal spot region of the micro-capillary array. Hence, tomographic slices at various depths near the focal plane can be reconstructed in a way similar to tomosynthesis, but from a single x-ray exposure. The microscope enables depth-resolved imaging of small subvolumes in large samples and can be used for imaging of weakly absorbing artificial and biological objects by means of propagation phase-contrast.


title: On a canonical lift of Artin's representation to loop braid groups

abstract: Each pointed topological space has an associated $\pi$-module, obtained from action of its first homotopy group on its second homotopy group. For the $3$-ball with a trivial link with $n$-components removed from its interior, its $\pi$-module $\mathcal{M}_n$ is of free type. In this paper we give an injection of the (extended) loop braid group into the group of automorphisms of $\mathcal{M}_n$. We give a topological interpretation of this injection, showing that it is both an extension of Artin's representation for braid groups and of Dahm's homomorphism for (extended) loop braid groups.


title: Convergence and sample complexity of gradient methods for the model-free linear quadratic regulator problem

abstract: Model-free reinforcement learning attempts to find an optimal control action for an unknown dynamical system by directly searching over the parameter space of controllers. The convergence behavior and statistical properties of these approaches are often poorly understood because of the nonconvex nature of the underlying optimization problems and the lack of exact gradient computation. In this paper, we take a step towards demystifying the performance and efficiency of such methods by focusing on the standard infinite-horizon linear quadratic regulator problem for continuous-time systems with unknown state-space parameters. We establish exponential stability for the ordinary differential equation (ODE) that governs the gradient-flow dynamics over the set of stabilizing feedback gains and show that a similar result holds for the gradient descent method that arises from the forward Euler discretization of the corresponding ODE. We also provide theoretical bounds on the convergence rate and sample complexity of the random search method with two-point gradient estimates. We prove that the required simulation time for achieving $\epsilon$-accuracy in the model-free setup and the total number of function evaluations both scale as $\log \, (1/\epsilon)$.

categories: math.OC cs.AI cs.LG cs.SY eess.SY physics.data-an


title: A Multilevel Stochastic Gradient method for PDE-constrained Optimal Control Problems with uncertain parameters

abstract: In this paper, we present a multilevel Monte Carlo (MLMC) version of the Stochastic Gradient (SG) method for optimization under uncertainty, in order to tackle Optimal Control Problems (OCP) where the constraints are described in the form of PDEs with random parameters. The deterministic control acts as a distributed forcing term in the random PDE and the objective function is an expected quadratic loss. We use a Stochastic Gradient approach to compute the optimal control, where the steepest descent direction of the expected loss, at each iteration, is replaced by independent MLMC estimators with increasing accuracy and computational cost. The refinement strategy is chosen a-priori such that the bias and, possibly, the variance of the MLMC estimator decays as a function of the iteration counter. Detailed convergence and complexity analyses of the proposed strategy are presented and asymptotically optimal decay rates are identified such that the total computational work is minimized. We also present and analyze an alternative version of the multilevel SG algorithm that uses a randomized MLMC estimator at each iteration. Our methodology is validated through a numerical example.


title: Random trigonometric polynomials: universality and non-universality of the variance for the number of real roots

abstract: In this paper, we study the number of real roots of random trigonometric polynomials with iid coefficients. When the coefficients have zero mean, unit variance and some finite high moments, we show that the variance of the number of real roots is asymptotically linear in terms of the expectation; furthermore, the multiplicative constant in this linear relationship depends only on the kurtosis of the common distribution of the polynomial's coefficients. This result is in sharp contrast to the classical Kac polynomials whose corresponding variance depends only on the first two moments. Our result is perhaps the first paper to establish the variance for general distribution of the coefficients including discrete ones, for a model of random polynomials outside the family of the Kac polynomials. Our method gives a fine comparison framework throughout Edgeworth expansion, asymptotic Kac-Rice formula and a detailed analysis of characteristic functions.


title: Minimization principle for shear alignment of liquid crystals

abstract: If a static perturbation is applied to a liquid crystal, the director configuration changes to minimize the free energy. If a shear flow is applied to a liquid crystal, one might ask: Does the director configuration change to minimize any effective potential? To address that question, we derive the Leslie-Ericksen equations for dissipative dynamics, and determine whether they can be expressed as relaxation toward a minimum. The answer may be yes or no, depending on the number of degrees of freedom. Using theory and simulations, we consider two specific examples, reverse tilt domains under simple shear flow and dowser configurations under plane Poiseuille flow, and demonstrate that each example shows relaxation toward the minimum of an effective potential.


title: A simple baseline for domain adaptation using rotation prediction

abstract: Recently, domain adaptation has become a hot research area with lots of applications. The goal is to adapt a model trained in one domain to another domain with scarce annotated data. We propose a simple yet effective method based on self-supervised learning that outperforms or is on par with most state-of-the-art algorithms, e.g. adversarial domain adaptation. Our method involves two phases: predicting random rotations (self-supervised) on the target domain along with correct labels for the source domain (supervised), and then using self-distillation on the target domain. Our simple method achieves state-of-the-art results on semi-supervised domain adaptation on DomainNet dataset. Further, we observe that the unlabeled target datasets of popular domain adaptation benchmarks do not contain any categories apart from testing categories. We believe this introduces a bias that does not exist in many real applications. We show that removing this bias from the unlabeled data results in a large drop in performance of state-of-the-art methods, while our simple method is relatively robust.


title: On Universal Scaling of Distributed Queues under Load Balancing

abstract: This paper considers the steady-state performance of load balancing algorithms in a many-server system with distributed queues. The system has $N$ servers, and each server maintains a local queue with buffer size $b-1,$ i.e. a server can hold at most one job in service and $b-1$ jobs in the queue. Jobs in the same queue are served according to the first-in-first-out (FIFO) order. The system is operated in a heavy-traffic regime such that the workload per server is $\lambda = 1 - N^{-\alpha}$ for $0.5\leq \alpha<1.$ We identify a set of algorithms such that the steady-state queues have the following universal scaling, where {\em universal} means that it holds for any $\alpha\in[0.5,1)$: (i) the number of of busy servers is $\lambda N-o(1);$ and (ii) the number of servers with two jobs (one in service and one in queue) is $O(N^{\alpha}\log N);$ and (iii) the number of servers with more than two jobs is $O\left(\frac{1}{N^{r(1-\alpha)-1}}\right),$ where $r$ can be any positive integer independent of $N.$ The set of load balancing algorithms that satisfy the sufficient condition includes join-the-shortest-queue (JSQ), idle-one-first (I1F), and power-of-$d$-choices (Po$d$) with $d\geq N^\alpha\log^2 N.$ We further argue that the waiting time of such an algorithm is near optimal order-wise.


title: Congruence pairs of principal MS-algebras and perfect extensions

abstract: The notion of a congruence pair for principal MS-algebras, simpler than the one given by Beazer for $K_2$-algebras \cite{6}, is introduced. It is proved that the congruences of the principal MS-algebras $L$ correspond to the MS-congruence pairs on simpler substructures $L^{\circ\circ}$ and $D(L)$ of $L$ that were associated to~$L$ in \cite{4}. An analogy of a well-known Gr\"atzer's problem \cite[Problem 57]{11} formulated for distributive p-algebras, which asks for a characterization of the congruence lattices in terms of the congruence pairs, is presented here for the principal MS-algebras (Problem 1). Unlike a recent solution to such a problem for the principal p-algebras in \cite{2}, it is demonstrated here on the class of principal MS-algebras, that a possible solution to the problem, though not very descriptive, can be simple and elegant. As a step to a more descriptive solution of Problem 1, a special case is then considered when a principal MS-algebra $L$ is a perfect extension of its greatest Stone subalgebra $L_{S}$. It is shown that this is exactly when de Morgan subalgebra $L^{\circ\circ}$ of $L$ is a perfect extension of the Boolean algebra $B(L)$. Two examples illustrating when this special case happens and when it does not are presented.


title: A Polynomial Time Algorithm for Computing the Strong Rainbow Connection Numbers of Odd Cacti

abstract: We consider the problem of computing the strong rainbow connection number $src(G)$ for cactus graphs $G$ in which all cycles have odd length. We present a formula to calculate $src(G)$ for such odd cacti which can be evaluated in linear time, as well as an algorithm for computing the corresponding optimal strong rainbow edge coloring, with polynomial worst case run time complexity. Although computing $src(G)$ is NP-hard in general, previous work has demonstrated that it may be computed in polynomial time for certain classes of graphs, including cycles, trees and block clique graphs. This work extends the class of graphs for which $src(G)$ may be computed in polynomial time.


title: A variational approach for many-body systems at finite temperature

abstract: We introduce a non-linear differential flow equation for density matrices that provides a monotonic decrease of the free energy and reaches a fixed point at the Gibbs thermal state. We use this equation to build a variational approach for analyzing equilibrium states of many-body systems and demonstrate that it can be applied to a broad class of states, including all bosonic and fermionic Gaussian states, as well as their generalizations obtained by unitary transformations, such as polaron transformations, in electron-phonon systems. We benchmark this method with a BCS lattice Hamiltonian and apply it to the Holstein model in two dimensions. For the latter, our approach reproduces the transition between the BCS pairing regime at weak interactions and the polaronic regime at stronger interactions, displaying phase separation between superconducting and charge-density wave phases.


title: A Novel Sharp Interface Immersed Boundary Framework for Viscous Flow Simulations at Arbitrary Mach Number Involving Complex and Moving Boundaries

abstract: This work presents a robust and efficient sharp interface immersed boundary (IBM) framework, which is applicable for all-speed flow regimes and is capable of handling arbitrarily complex bodies (stationary or moving). The work deploys an in-house, parallel, multi-block structured finite volume flow solver, which employs a 3D unsteady Favre averaged Navier Stokes equations in a generalized curvilinear coordinate system; while we employ a combination of HCIB (Hybrid Cartesian Immersed boundary) method and GC(Ghost-cell) for solution reconstruction near immersed boundary interface. A significant difficulty for these sharp interface approaches is of handling sharp features/edges of complex geometries. In this study, we observe that apart from the need for robust node classification strategy and higher order boundary formulations, the direction in which the reconstruction procedures are performed plays an important role in handling sharp edges. Taking this into account we present a versatile interface tracking procedure based on ray tracing algorithm and a novel three step solution reconstruction procedure that computes pseudo-normals in the regions where the normal is not well-defined and reconstructs the flow field along those directions. We demonstrate that this procedure enables solver to efficiently handle and accurately represent sharp-edged regions. A fifth-order weighted essentially non-oscillatory (WENO) scheme is used for capturing shock-induced discontinuities and complex fluid-solid interactions with high resolution. The developed IBM framework is applied to a wide range of flow phenomena encompassing all-speed regimes (M=0.001 to M = 2.0). A total of seven benchmark cases (three stationary and four moving bodies) are presented involving various geometries (cylinder, airfoil, wedge) and the predictions are found to be in excellent agreement with the published results.


title: Photon Orbits and Thermodynamic Phase Transition of Regular AdS Black Holes

abstract: We probe the phase structure of the regular AdS black holes using the null geodesics. The radius of photon orbit and minimum impact parameter shows a non-monotonous behaviour below the critical values of the temperature and the pressure, corresponding to the phase transition in extended phase space. The respective differences of the radius of unstable circular orbit and the minimum impact parameter can be seen as the order parameter for the small-large black hole phase transition, with a critical exponent $1/2$. Our study shows that there exists a close relationship between the gravity and thermodynamics for the regular AdS black holes.


title: Phase transitions for infinite products of large non-Hermitian random matrices

abstract: Products of $M$ i.i.d. non-Hermitian random matrices of size $N \times N$ relate Gaussian fluctuation of Lyapunov and stability exponents in dynamical systems (finite $N$ and large $M$) to local eigenvalue universality in random matrix theory (finite $M$ and large $N$). The remaining task is to study local eigenvalue statistics as $M$ and $N$ tend to infinity simultaneously, which lies at the heart of understanding two kinds of universal patterns. For products of i.i.d. complex Ginibre matrices, truncated unitary matrices and spherical ensembles, as $M+N\to \infty$ we prove that local statistics undergoes a transition when the relative ratio $M/N$ changes from $0$ to $\infty$: Ginibre statistics when $M/N \to 0$, normality when $M/N\to \infty$, and new critical phenomena when $M/N\to \gamma \in (0, \infty)$.


title: Graded-division algebras over arbitrary fields

abstract: A graded-division algebra is an algebra graded by a group such that all nonzero homogeneous elements are invertible. This includes division algebras equipped with an arbitrary group grading (including the trivial grading). We show that a classification of finite-dimensional graded-central graded-division algebras over an arbitrary field $\mathbb{F}$ can be reduced to the following three classifications, for each finite Galois extension $\mathbb{L}$ of $\mathbb{F}$: (1) finite-dimensional central division algebras over $\mathbb{L}$, up to isomorphism; (2) twisted group algebras of finite groups over $\mathbb{L}$, up to graded-isomorphism; (3) $\mathbb{F}$-forms of certain graded matrix algebras with coefficients in $\Delta\otimes_{\mathbb{L}}\mathcal{C}$ where $\Delta$ is as in (1) and $\mathcal{C}$ is as in (2). As an application, we classify, up to graded-isomorphism, the finite-dimensional graded-division algebras over the field of real numbers (or any real closed field) with an abelian grading group. We also discuss group gradings on fields.


abstract: We propose a trust region method for policy optimization that employs Quasi-Newton approximation for the Hessian, called Quasi-Newton Trust Region Policy Optimization QNTRPO. Gradient descent is the de facto algorithm for reinforcement learning tasks with continuous controls. The algorithm has achieved state-of-the-art performance when used in reinforcement learning across a wide range of tasks. However, the algorithm suffers from a number of drawbacks including: lack of stepsize selection criterion, and slow convergence. We investigate the use of a trust region method using dogleg step and a Quasi-Newton approximation for the Hessian for policy optimization. We demonstrate through numerical experiments over a wide range of challenging continuous control tasks that our particular choice is efficient in terms of number of samples and improves performance


title: Category-Level Articulated Object Pose Estimation

abstract: This project addresses the task of category-level pose estimation for articulated objects from a single depth image. We present a novel category-level approach that correctly accommodates object instances previously unseen during training. We introduce Articulation-aware Normalized Coordinate Space Hierarchy (ANCSH) - a canonical representation for different articulated objects in a given category. As the key to achieve intra-category generalization, the representation constructs a canonical object space as well as a set of canonical part spaces. The canonical object space normalizes the object orientation,scales and articulations (e.g. joint parameters and states) while each canonical part space further normalizes its part pose and scale. We develop a deep network based on PointNet++ that predicts ANCSH from a single depth point cloud, including part segmentation, normalized coordinates, and joint parameters in the canonical object space. By leveraging the canonicalized joints, we demonstrate: 1) improved performance in part pose and scale estimations using the induced kinematic constraints from joints; 2) high accuracy for joint parameter estimation in camera space.


title: Inverses of Matern Covariances on Grids

abstract: We conduct a study of the aliased spectral densities of Mat\'ern covariance functions on a regular grid of points, providing clarity on the properties of a popular approximation based on stochastic partial differential equations; while others have shown that it can approximate the covariance function well, we find that it assigns too much power at high frequencies and does not provide increasingly accurate approximations to the inverse as the grid spacing goes to zero, except in the one-dimensional exponential covariance case. We provide numerical results to support our theory, and in a simulation study, we investigate the implications for parameter estimation, finding that the SPDE approximation tends to overestimate spatial range parameters.


title: Minimization of atomic displacements as a guiding principle of the martensitic phase transformation

abstract: We present a unifying description for the martensitic transformation of steel that accounts for important experimentally observable features of the transformation namely, the Neumann bands, the interfacial (habit) plane between the transformed and untransformed phases and their orientation relationship (OR). It is obtained through a simple geometric minimization of the total distance traveled by all the atoms from the austenite (FCC or $\gamma$) phase to the martensite (BCC or $\alpha$) phase, without the need for any explicit energy minimization. Our description unites previously proposed mechanisms but it does not rely on assumptions and experimental knowledge regarding the shear planes and directions, or external adjustable parameters. We show how the Kurdjumov-Sach orientation relationship between the two phases and the $\{225\}_\gamma$ habit plane, which have both been extensively reported in experiments, naturally emerge from the distance minimization. We also propose an explanation for the occurrence of a different orientation relationship (Pitsch) in thin films.


title: Strong and weak lensing of Gravitational Waves: a semi-analytical approach

abstract: In this paper we study gravitational lensing of gravitational wave events. The main point of the present work is to introduce a semi-analytic approach so that each ingredient can be varied and tested individually. Our analytic model for the source population is motivated by a numerical study and we compare semi-analytical and numerical results. We determine the expected magnification for events seen at a given observed luminosity distance. We find that while the probability of significant magnification of the observed LIGO-Virgo events is very small, the probability distribution of the magnification has a significant tail to high $\mu$ such that e.g. the variance of the magnification is very large and even diverges in the geometric optics approximation. For the 10 binary black hole mergers observed by LIGO-Virgo in the O1+O2 observation campaigns, the probability that one of them has been magnified with magnification of 5 or bigger is $\mathcal{P}_{\text{obs}}(>5)\sim 0.01$ while the probability of magnification by 50 or bigger is $\mathcal{P}_{\text{obs}}(>50)\sim 0.005$.


title: Cancellative hypergraphs and Steiner triple systems

abstract: A triple system is cancellative if it does not contain three distinct sets $A,B,C$ such that the symmetric difference of $A$ and $B$ is contained in $C$. We show that every cancellative triple system $\mathcal{H}$ that satisfies certain inequality between the sizes of $\mathcal{H}$ and its shadow must be structurally close to the balanced blowup of some Steiner triple system. Our result contains a stability theorem for cancellative triple systems due to Keevash and Mubayi as a special case. It also implies that the boundary of the feasible region of cancellative triple systems has infinitely many local maxima, thus giving the first example showing this phenomenon.


title: Grassmanniennes affines tordues sur les entiers

abstract: We extend the work of Pappas-Rapoport-Zhu on twisted affine Grassmannians to wildly ramified, quasi-split, and residually split groups, assuming the maximal torus is induced. This relies on the construction, inspired by Tits, of certain smooth, affine, and connected $\mathbb{Z}[t]$-groups of parahoric type, which should be regarded as $\mathbb{Z}$-families of parahoric group schemes, and naturally extends a similar construction in the above articles after inverting $e$. The resulting $\mathbb{F}_p(t)$-groups are pseudo-reductive possibly non-standard in the sense of Conrad--Gabber--Prasad, and their $\mathbb{F}_p[[t]]$-models are parahoric in our generalized sense. We study their affine Grassmannians, establishing normality of Schubert varieties and Zhu's coherence theorem.


title: A new spectral method based on two classes of hat functions for solving systems of fractional differential equations and an application to respiratory syncytial virus infection

abstract: We propose a new spectral method, based on two classes of hat functions, for solving systems of fractional differential equations. The fractional derivative is considered in the Caputo sense. Properties of the basis functions, Caputo derivatives, and Riemann-Liouville fractional integrals, are used to reduce the main problem to a system of nonlinear algebraic equations. By analyzing in detail the resulting system, we show that the method needs few computational efforts. Two test problems are considered to illustrate the efficiency and accuracy of the proposed method. Finally, an application to a recent mathematical model in epidemiology is given, precisely to a system of fractional differential equations modeling the respiratory syncytial virus infection.


title: Diversity-Oriented Synthesis of Polymers of Intrinsic Microporosity with Explicit Solid Solvation Cages for Lithium Ions

abstract: Here, we describe a diversity-oriented synthetic strategy for microporous polymer membranes from which we identified those whose FVEs serve as solid solvation cages for lithium ions. Lead candidate membranes featuring such ion solvation cages exhibited both higher ionic conductivity and higher cation transference number than control membranes where FVEs were aspecific, which indicates conventional bounds for membrane permeability and selectivity for ion transport can be overcome. Such membranes show promise as dendrite-suppressing anode-electrolyte interlayers in high-voltage lithium-metal batteries for electric mobility.


title: The volumes of Miyauchi subgroups

abstract: Miyauchi described the $L$ and $\varepsilon$-factors attached to generic representations of the unramified unitary group of rank three in terms of local newforms defined by a sequence of subgroups. We calculate the volumes of these Miyauchi groups.


title: A proof of the Brown--Goodearl Conjecture for weak Hopf algebras

abstract: Let $H$ be a weak Hopf algebra that is a finitely generated module over its affine center. We show that $H$ has finite self-injective dimension and so the Brown--Goodearl Conjecture holds in this special weak Hopf setting.


title: Spectral Fluctuations in the Sachdev-Ye-Kitaev Model

abstract: We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being the scale fluctuations for which we give a simple estimate. Our numerical results for $N=32$ show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion is obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations using a nonlinear $\sigma$-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlator to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the $\sigma$-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for $q=4$ but also for the $q=2$ SYK model while the correct result in this case should have been Poisson statistics. We analyze the number variance and spectral form factor for $N=32$ and $q=4$ numerically. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings or exponentially short times, respectively.


title: Well-Posedness of Free Boundary Problem in Non-relativistic and Relativistic Ideal Compressible Magnetohydrodynamics

abstract: We consider the free boundary problem for non-relativistic and relativistic ideal compressible magnetohydrodynamics in two and three spatial dimensions with the total pressure vanishing on the plasma--vacuum interface. We establish the local-in-time existence and uniqueness of solutions to this nonlinear characteristic hyperbolic problem under the Rayleigh--Taylor sign condition on the total pressure. The proof is based on certain tame estimates in anisotropic Sobolev spaces for the linearized problem and a modification of the Nash--Moser iteration scheme. Our result is uniform in the light speed and appears to be the first well-posedness result for the free boundary problem in ideal compressible magnetohydrodynamics with zero total pressure on the moving boundary.


title: Scattering-based geometric shaping of photon-photon interactions

abstract: We construct an effective Hamiltonian of interacting bosons, based on scattered radiation off vibrational modes of designed molecular architectures. Making use of the infinite yet countable set of spatial modes representing the scattering of light, we obtain a variable photon-photon interaction in this basis. The effective Hamiltonian hermiticity is controlled by a geometric factor set by the overlaps of spatial modes. Using this mapping, we relate intensity measurements of the light to correlation functions of the interacting bosons evolving according to the effective Hamiltonian, rendering local as well as nonlocal observables accessible. This architecture may be used to simulate the dynamics of interacting bosons, as well as designing tool for multi-qubit photonic gates in quantum computing applications. Variable hopping, interaction and confinement of the active space of the bosons are demonstrated on a model system.


title: Parameter Free Clustering with Cluster Catch Digraphs (Technical Report)

abstract: We propose clustering algorithms based on a recently developed geometric digraph family called cluster catch digraphs (CCDs). These digraphs are used to devise clustering methods that are hybrids of density-based and graph-based clustering methods. CCDs are appealing digraphs for clustering, since they estimate the number of clusters; however, CCDs (and density-based methods in general) require some information on a parameter representing the \emph{intensity} of assumed clusters in the data set. We propose algorithms that are parameter free versions of the CCD algorithm and does not require a specification of the intensity parameter whose choice is often critical in finding an optimal partitioning of the data set. We estimate the number of convex clusters by borrowing a tool from spatial data analysis, namely Ripley's $K$ function. We call our new digraphs utilizing the $K$ function as RK-CCDs. We show that the minimum dominating sets of RK-CCDs estimate and distinguish the clusters from noise clusters in a data set, and hence allow the estimation of the correct number of clusters. Our robust clustering algorithms are comprised of methods that estimate both the number of clusters and the intensity parameter, making them completely parameter free. We conduct Monte Carlo simulations and use real life data sets to compare RK-CCDs with some commonly used density-based and prototype-based clustering methods.


title: Towards finding hay in a haystack: explicit tensors of border rank greater than $2.02m$ in $C^m\otimes C^m\otimes C^m$

abstract: We write down an explicit sequence of tensors in $C^m\otimes C^m\otimes C^m$, for all $m$ sufficiently large, having border rank at least $2.02m$, overcoming a longstanding barrier. We obtain our lower bounds via the border substitution method.


title: Meta-analysis of heterogeneous data: integrative sparse regression in high-dimensions

abstract: We consider the task of meta-analysis in high-dimensional settings in which the data sources are similar but non-identical. To borrow strength across such heterogeneous datasets, we introduce a global parameter that emphasizes interpretability and statistical efficiency in the presence of heterogeneity. We also propose a one-shot estimator of the global parameter that preserves the anonymity of the data sources and converges at a rate that depends on the size of the combined dataset. For high-dimensional linear model settings, we demonstrate the superiority of our identification restrictions in adapting to a previously seen data distribution as well as predicting for a new/unseen data distribution. Finally, we demonstrate the benefits of our approach on a large-scale drug treatment dataset involving several different cancer cell-lines.


title: GASOL: Gas Analysis and Optimization for Ethereum Smart Contracts

abstract: We present the main concepts, components, and usage of GASOL, a Gas AnalysiS and Optimization tooL for Ethereum smart contracts. GASOL offers a wide variety of cost models that allow inferring the gas consumption associated to selected types of EVM instructions and/or inferring the number of times that such types of bytecode instructions are executed. Among others, we have cost models to measure only storage opcodes, to measure a selected family of gas-consumption opcodes following the Ethereum's classification, to estimate the cost of a selected program line, etc. After choosing the desired cost model and the function of interest, GASOL returns to the user an upper bound of the cost for this function. As the gas consumption is often dominated by the instructions that access the storage, GASOL uses the gas analysis to detect under-optimized storage patterns, and includes an (optional) automatic optimization of the selected function. Our tool can be used within an Eclipse plugin for Solidity which displays the gas and instructions bounds and, when applicable, the gas-optimized Solidity function.


title: Strategies for the vectorized Block Conjugate Gradients method

abstract: Block Krylov methods have recently gained a lot of attraction. Due to their increased arithmetic intensity they offer a promising way to improve performance on modern hardware. Recently Frommer et al. presented a block Krylov framework that combines the advantages of block Krylov methods and data parallel methods. We review this framework and apply it on the Block Conjugate Gradients method,to solve linear systems with multiple right hand sides. In this course we consider challenges that occur on modern hardware, like a limited memory bandwidth, the use of SIMD instructions and the communication overhead. We present a performance model to predict the efficiency of different Block CG variants and compare these with experimental numerical results.


title: Sparse Optimization on General Atomic Sets: Greedy and Forward-Backward Algorithms

abstract: We consider the problem of sparse atomic optimization, where the notion of "sparsity" is generalized to meaning some linear combination of few atoms. The definition of atomic set is very broad; popular examples include the standard basis, low-rank matrices, overcomplete dictionaries, permutation matrices, orthogonal matrices, etc. The model of sparse atomic optimization therefore includes problems coming from many fields, including statistics, signal processing, machine learning, computer vision and so on. Specifically, we consider the problem of maximizing a restricted strongly convex (or concave), smooth function restricted to a sparse linear combination of atoms. We extend recent work that establish linear convergence rates of greedy algorithms on restricted strongly concave, smooth functions on sparse vectors to the realm of general atomic sets, where the convergence rate involves a novel quantity: the "sparse atomic condition number". This leads to the strongest known multiplicative approximation guarantees for various flavors of greedy algorithms for sparse atomic optimization; in particular, we show that in many settings of interest the greedy algorithm can attain strong approximation guarantees while maintaining sparsity. Furthermore, we introduce a scheme for forward-backward algorithms that achieves the same approximation guarantees. Secondly, we define an alternate notion of weak submodularity, which we show is tightly related to the more familiar version that has been used to prove earlier linear convergence rates. We prove analogous multiplicative approximation guarantees using this alternate weak submodularity, and establish its distinct identity and applications.


title: Skeleton Extraction from 3D Point Clouds by Decomposing the Object into Parts

abstract: Decomposing a point cloud into its components and extracting curve skeletons from point clouds are two related problems. Decomposition of a shape into its components is often obtained as a byproduct of skeleton extraction. In this work, we propose to extract curve skeletons, from unorganized point clouds, by decomposing the object into its parts, identifying part skeletons and then linking these part skeletons together to obtain the complete skeleton. We believe it is the most natural way to extract skeletons in the sense that this would be the way a human would approach the problem. Our parts are generalized cylinders (GCs). Since, the axis of a GC is an integral part of its definition, the parts have natural skeletal representations. We use translational symmetry, the fundamental property of GCs, to extract parts from point clouds. We demonstrate how this method can handle a large variety of shapes. We compare our method with state of the art methods and show how a part based approach can deal with some of the limitations of other methods. We present an improved version of an existing point set registration algorithm and demonstrate its utility in extracting parts from point clouds. We also show how this method can be used to extract skeletons from and identify parts of noisy point clouds. A part based approach also provides a natural and intuitive interface for user interaction. We demonstrate the ease with which mistakes, if any, can be fixed with minimal user interaction with the help of a graphical user interface.


title: Monotonicity considerations for stabilized DG cut cell schemes for the unsteady advection equation

abstract: For solving unsteady hyperbolic conservation laws on cut cell meshes, the so called small cell problem is a big issue: one would like to use a time step that is chosen with respect to the background mesh and use the same time step on the potentially arbitrarily small cut cells as well. For explicit time stepping schemes this leads to instabilities. In a recent preprint [arXiv:1906.05642], we propose penalty terms for stabilizing a DG space discretization to overcome this issue for the unsteady linear advection equation. The usage of the proposed stabilization terms results in stable schemes of first and second order in one and two space dimensions. In one dimension, for piecewise constant data in space and explicit Euler in time, the stabilized scheme can even be shown to be monotone. In this contribution, we will examine the conditions for monotonicity in more detail.


title: Bounded and Almost Periodic Solvability of Nonautonomous Quasilinear Hyperbolic Systems

abstract: The paper concerns boundary value problems for general nonautonomous first order quasilinear hyperbolic systems in a strip. We construct small global classical solutions, assuming that the right hand sides are small. In the case that all data of the quasilinear problem are almost periodic, we prove that the bounded solution is also almost periodic. For the nonhomogeneous version of a linearized problem, we provide stable dissipativity conditions ensuring a unique bounded continuous solution for any smooth right-hand sides. In the autonomous case, this solution is two times continuously differentiable. In the nonautonomous case, the continuous solution is differentiable under additional dissipativity conditions, which are essential. A crucial ingredient of our approach is a perturbation theorem for general linear hyperbolic systems. One of the technical complications we overcome is the "loss of smoothness" property of hyperbolic PDEs.


title: Secular Extragalactic Parallax: Measurement Methods and Predictions for Gaia

abstract: Secular extragalactic parallax caused by the solar system's velocity relative to the cosmic microwave background rest frame may be observable as a dipole proper motion field with amplitude $78~\mu$as yr$^{-1}$ Mpc. Nearby galaxies also exhibit proper motions caused by their transverse peculiar velocities that prevent detection of secular parallax for any single galaxy, although a statistical detection may be made instead. Such a detection could constrain the local Hubble parameter. We present methods to measure secular parallax using correlated extragalactic proper motions and find a first limit on the secular parallax amplitude using proper motions of 232 nearby galaxies from Gaia Data Release 2. The recovered dipole has insignificant upper limit of 3500 $\mu$as yr$^{-1}$ Mpc. This measurement will be improved by larger sample size and reduced proper motion uncertainties in future data releases. Using the local peculiar velocity field derived from Cosmicflows-3, we simulate galaxy proper motions and predict that a significant detection ($5-10\sigma$) of the secular parallax amplitude will be possible by Gaia's end of mission. The detection is contingent on proper motions of nearby ($<5$ Mpc), bright ($G<15$ mag) galaxies, and corresponds to an insignificant upper limit on the Hubble parameter. We further investigate the implications of our simulations for the study of transverse peculiar velocities, which we find to be consistent with large scale structure theory. The peculiar velocity field additionally results in low-multipole correlated proper motions on the order of $0.3~\mu$as yr$^{-1}$ that may be confounded with other cosmological proper motion measurements, such as limits on the gravitational wave background and the anisotropy of the Hubble expansion.


title: Smell Pittsburgh: Engaging Community Citizen Science for Air Quality

abstract: Urban air pollution has been linked to various human health concerns, including cardiopulmonary diseases. Communities who suffer from poor air quality often rely on experts to identify pollution sources due to the lack of accessible tools. Taking this into account, we developed Smell Pittsburgh, a system that enables community members to report odors and track where these odors are frequently concentrated. All smell report data are publicly accessible online. These reports are also sent to the local health department and visualized on a map along with air quality data from monitoring stations. This visualization provides a comprehensive overview of the local pollution landscape. Additionally, with these reports and air quality data, we developed a model to predict upcoming smell events and send push notifications to inform communities. We also applied regression analysis to identify statistically significant effects of push notifications on user engagement. Our evaluation of this system demonstrates that engaging residents in documenting their experiences with pollution odors can help identify local air pollution patterns, and can empower communities to advocate for better air quality. All citizen-contributed smell data are publicly accessible and can be downloaded from https://smellpgh.org.


title: Energy-Based Weak Measurement

abstract: When a well-localized photon is incident on a spatially superposed absorber but is not absorbed, the photon can still deliver energy to the absorber. It is shown that when the transferred energy is small relative to the energy uncertainty of the photon, this constitutes an unusual type of weak measurement of the absorber's energy, where the energy distribution of the unabsorbed photon acts as the measurement device, and the strongly disturbed state of the absorber becomes the effective pre-selection. Treating the final state of the absorber as the post-selection, it is shown that the absorber's energy increase is the weak value of its translational Hamiltonian, and the energy distribution of the photon shifts by the opposite amount. The basic case of non-scattering is examined, followed by the case of interaction-free energy transfer. Details and interpretations of the results are discussed.


title: On the projective description of spaces of ultradifferentiable functions of Roumieu type

abstract: We provide a projective description of the space $\mathcal{E}^{\{\mathfrak{M}\}}(\Omega)$ of ultradifferentiable functions of Roumieu type, where $\Omega$ is an arbitrary open set in $\mathbb{R}^d$ and $\mathfrak{M}$ is a weight matrix satisfying the analogue of Komatsu's condition $(M.2)'$. In particular, we obtain in a unified way projective descriptions of ultradifferentiable classes defined via a single weight sequence (Denjoy-Carleman approach) and via a weight function (Braun-Meise-Taylor approach) under considerably weaker assumptions than in earlier versions of these results.


title: On the Principle of Least Symmetry Breaking in Shallow ReLU Models

abstract: We consider the optimization problem associated with fitting two-layer ReLU networks with respect to the squared loss, where labels are assumed to be generated by a target network. Focusing first on standard Gaussian inputs, we show that the structure of spurious local minima detected by stochastic gradient descent (SGD) is, in a well-defined sense, the \emph{least loss of symmetry} with respect to the target weights. A closer look at the analysis indicates that this principle of least symmetry breaking may apply to a broader range of settings. Motivated by this, we conduct a series of experiments which corroborate this hypothesis for different classes of non-isotropic non-product distributions, smooth activation functions and networks with a few layers.


title: Towards Better Understanding of Adaptive Gradient Algorithms in Generative Adversarial Nets

abstract: Adaptive gradient algorithms perform gradient-based updates using the history of gradients and are ubiquitous in training deep neural networks. While adaptive gradient methods theory is well understood for minimization problems, the underlying factors driving their empirical success in min-max problems such as GANs remain unclear. In this paper, we aim at bridging this gap from both theoretical and empirical perspectives. First, we analyze a variant of Optimistic Stochastic Gradient (OSG) proposed in~\citep{daskalakis2017training} for solving a class of non-convex non-concave min-max problem and establish $O(\epsilon^{-4})$ complexity for finding $\epsilon$-first-order stationary point, in which the algorithm only requires invoking one stochastic first-order oracle while enjoying state-of-the-art iteration complexity achieved by stochastic extragradient method by~\citep{iusem2017extragradient}. Then we propose an adaptive variant of OSG named Optimistic Adagrad (OAdagrad) and reveal an \emph{improved} adaptive complexity $O\left(\epsilon^{-\frac{2}{1-\alpha}}\right)$, where $\alpha$ characterizes the growth rate of the cumulative stochastic gradient and $0\leq \alpha\leq 1/2$. To the best of our knowledge, this is the first work for establishing adaptive complexity in non-convex non-concave min-max optimization. Empirically, our experiments show that indeed adaptive gradient algorithms outperform their non-adaptive counterparts in GAN training. Moreover, this observation can be explained by the slow growth rate of the cumulative stochastic gradient, as observed empirically.


title: On the obstruction to the Hasse principle for multinorm equations

abstract: We investigate local-global principles for multinorm equations over a global field. To this extent, we generalize work of Drakokhrust and Platonov to provide explicit and computable formulae for the obstructions to the Hasse principle and weak approximation for multinorm equations. We illustrate the scope of this technique by extending results of Bayer-Fluckiger--Lee--Parimala, Demarche--Wei and Pollio.


title: On the Beilinson-Bloch-Kato conjecture for Rankin-Selberg motives

abstract: In this article, we study the Beilinson-Bloch-Kato conjecture for motives corresponding to the Rankin-Selberg product of conjugate self-dual automorphic representations, within the framework of the Gan-Gross-Prasad conjecture. We show that if the central critical value of the Rankin-Selberg $L$-function does not vanish, then the Bloch-Kato Selmer group with coefficients in a favorable field of the corresponding motive vanishes. We also show that if the class in the Bloch-Kato Selmer group constructed from certain diagonal cycle does not vanish, which is conjecturally equivalent to the nonvanishing of the central critical first derivative of the Rankin-Selberg $L$-function, then the Bloch-Kato Selmer group is of rank one.


title: De-biasing convex regularized estimators and interval estimation in linear models

abstract: New upper bounds are developed for the $L_2$ distance between $\xi/\text{Var}[\xi]^{1/2}$ and linear and quadratic functions of $z\sim N(0,I_n)$ for random variables of the form $\xi=bz^\top f(z) - \text{div} f(z)$. The linear approximation yields a central limit theorem when the squared norm of $f(z)$ dominates the squared Frobenius norm of $\nabla f(z)$ in expectation. Applications of this normal approximation are given for the asymptotic normality of de-biased estimators in linear regression with correlated design and convex penalty in the regime $p/n \le \gamma$ for constant $\gamma\in(0,{\infty})$. For the estimation of linear functions $\langle a_0,\beta\rangle$ of the unknown coefficient vector $\beta$, this analysis leads to asymptotic normality of the de-biased estimate for most normalized directions $a_0$, where ``most'' is quantified in a precise sense. This asymptotic normality holds for any convex penalty if $\gamma<1$ and for any strongly convex penalty if $\gamma\ge 1$. In particular the penalty needs not be separable or permutation invariant. By allowing arbitrary regularizers, the results vastly broaden the scope of applicability of de-biasing methodologies to obtain confidence intervals in high-dimensions. In the absence of strong convexity for $p>n$, asymptotic normality of the de-biased estimate is obtained for the Lasso and the group Lasso under additional conditions. For general convex penalties, our analysis also provides prediction and estimation error bounds of independent interest.


title: On the Reproducibility of Experiments of Indexing Repetitive Document Collections

abstract: This work introduces a companion reproducible paper with the aim of allowing the exact replication of the methods, experiments, and results discussed in a previous work [5]. In that parent paper, we proposed many and varied techniques for compressing indexes which exploit that highly repetitive collections are formed mostly of documents that are near-copies of others. More concretely, we describe a replication framework, called uiHRDC (universal indexes for Highly Repetitive Document Collections), that allows our original experimental setup to be easily replicated using various document collections. The corresponding experimentation is carefully explained, providing precise details about the parameters that can be tuned for each indexing solution. Finally, note that we also provide uiHRDC as reproducibility package.


title: On the Morality of Artificial Intelligence

abstract: Much of the existing research on the social and ethical impact of Artificial Intelligence has been focused on defining ethical principles and guidelines surrounding Machine Learning (ML) and other Artificial Intelligence (AI) algorithms [IEEE, 2017, Jobin et al., 2019]. While this is extremely useful for helping define the appropriate social norms of AI, we believe that it is equally important to discuss both the potential and risks of ML and to inspire the community to use ML for beneficial objectives. In the present article, which is specifically aimed at ML practitioners, we thus focus more on the latter, carrying out an overview of existing high-level ethical frameworks and guidelines, but above all proposing both conceptual and practical principles and guidelines for ML research and deployment, insisting on concrete actions that can be taken by practitioners to pursue a more ethical and moral practice of ML aimed at using AI for social good.


title: Stochastic Dynamic Cutting Plane for multistage stochastic convex programs

abstract: We introduce StoDCuP (Stochastic Dynamic Cutting Plane), an extension of the Stochastic Dual Dynamic Programming (SDDP) algorithm to solve multistage stochastic convex optimization problems. At each iteration, the algorithm builds lower affine functions not only for the cost-to-go functions, as SDDP does, but also for some or all nonlinear cost and constraint functions. We show the almost sure convergence of StoDCuP. We also introduce an inexact variant of StoDCuP where all subproblems are solved approximately (with bounded errors) and show the almost sure convergence of this variant for vanishing errors.


title: Colorectal Polyp Segmentation by U-Net with Dilation Convolution

abstract: Colorectal cancer (CRC) is one of the most commonly diagnosed cancers and a leading cause of cancer deaths in the United States. Colorectal polyps that grow on the intima of the colon or rectum is an important precursor for CRC. Currently, the most common way for colorectal polyp detection and precancerous pathology is the colonoscopy. Therefore, accurate colorectal polyp segmentation during the colonoscopy procedure has great clinical significance in CRC early detection and prevention. In this paper, we propose a novel end-to-end deep learning framework for the colorectal polyp segmentation. The model we design consists of an encoder to extract multi-scale semantic features and a decoder to expand the feature maps to a polyp segmentation map. We improve the feature representation ability of the encoder by introducing the dilated convolution to learn high-level semantic features without resolution reduction. We further design a simplified decoder which combines multi-scale semantic features with fewer parameters than the traditional architecture. Furthermore, we apply three post processing techniques on the output segmentation map to improve colorectal polyp detection performance. Our method achieves state-of-the-art results on CVC-ClinicDB and ETIS-Larib Polyp DB.


title: Decoherence Effects Break Reciprocity in Matter Transport

abstract: The decoherence of quantum states defines the transition between the quantum world and classical physics. Decoherence or, analogously, quantum mechanical collapse events pose fundamental questions regarding the interpretation of quantum mechanics and are technologically relevant because they limit the coherent information processing performed by quantum computers. We have discovered that the transition regime enables a novel type of matter transport. Applying this discovery, we present nanoscale devices in which decoherence, modeled by random quantum jumps, produces fundamentally novel phenomena by interrupting the unitary dynamics of electron wave packets. Noncentrosymmetric conductors with mesoscopic length scales act as two-terminal rectifiers with unique properties. In these devices, the inelastic interaction of itinerant electrons with impurities acting as electron trapping centers leads to a novel steady state characterized by partial charge separation between the two leads, or, in closed circuits to the generation of persistent currents. The interface between the quantum and the classical worlds therefore provides a novel transport regime of value for the realization of a new category of mesoscopic electronic devices.


title: Emergence of stochastic flocking for the discrete Cucker-Smale model with randomly switching topologies

abstract: We study emergent dynamics of the discrete Cucker-Smale (in short, DCS) model with randomly switching network topologies. For this, we provide a sufficient framework leading to the stochastic flocking with probability one. Our sufficient framework is formulated in terms of an admissible set of network topologies realized by digraphs and probability density function for random switching times. As examples for the law of switching times, we use the Poisson process and the geometric process and show that these two processes satisfy the required conditions in a given framework so that we have a stochastic flocking with probability one. As a corollary of our flocking analysis, we improve the earlier result [J.-G. Dong, S.-Y. Ha, J. Jung and D. Kim: On the stochastic flocking of the Cucker-Smale flock with randomly switching topologies. arXiv:1911.07390.] on the continuous C-S model.


title: Leptonic Source of Dark Matter and Radiative Majorana or Dirac Neutrino Mass

abstract: The notion of U(1) lepton number (which may only be softly broken) is applied to models of dark matter which interacts with leptons. Previous scotogenic models of Majorana or Dirac neutrino mass are shown to be derivable in this framework without additional symmetries. Only complete renormalizable theories are considered. An explicit class of models with $Z_n$ ($n \geq 5$) lepton and dark symmetry for Dirac neutrinos is derived, as well as an example of $Z_3$ dark symmetry.


title: EVA: An Encrypted Vector Arithmetic Language and Compiler for Efficient Homomorphic Computation

abstract: Fully-Homomorphic Encryption (FHE) offers powerful capabilities by enabling secure offloading of both storage and computation, and recent innovations in schemes and implementations have made it all the more attractive. At the same time, FHE is notoriously hard to use with a very constrained programming model, a very unusual performance profile, and many cryptographic constraints. Existing compilers for FHE either target simpler but less efficient FHE schemes or only support specific domains where they can rely on expert-provided high-level runtimes to hide complications. This paper presents a new FHE language called Encrypted Vector Arithmetic (EVA), which includes an optimizing compiler that generates correct and secure FHE programs, while hiding all the complexities of the target FHE scheme. Bolstered by our optimizing compiler, programmers can develop efficient general-purpose FHE applications directly in EVA. For example, we have developed image processing applications using EVA, with a very few lines of code. EVA is designed to also work as an intermediate representation that can be a target for compiling higher-level domain-specific languages. To demonstrate this, we have re-targeted CHET, an existing domain-specific compiler for neural network inference, onto EVA. Due to the novel optimizations in EVA, its programs are on average 5.3x faster than those generated by CHET. We believe that EVA would enable a wider adoption of FHE by making it easier to develop FHE applications and domain-specific FHE compilers.


title: Reduction of the energy-gap scaling by coherent catalysis in models of quantum annealing

abstract: Non-stoquastic drivers are known to improve the performance of quantum annealing by reducing first-order phase transitions into second-order ones in several mean-field-type model systems. Nevertheless, statistical-mechanical analysis shows that some target Hamiltonians still exhibit unavoidable first-order transitions even with non-stoquastic drivers, making them difficult for quantum annealing to solve. Recently, a mechanism called coherent catalysis was proposed by Durkin [Phys. Rev. A \textbf{99}, 032315 (2019)], in which he showed the existence of a particular point on the line of first-order phase transitions where the energy gap scales polynomially as expected for a second-order transition. We show by extensive numerical computations that this phenomenon is observed in a few additional mean-field-type optimization problems where non-stoquastic drivers fail to change the order of phase transition in the thermodynamic limit. This opens up the possibility of using coherent catalysis to search for exponential speedups in systems previously thought to be exponentially slow for quantum annealing to solve.


title: Apricot variety classification using image processing and machine learning approaches

abstract: Apricot which is a cultivated type of Zerdali (wild apricot) has an important place in human nutrition and its medical properties are essential for human health. The objective of this research was to obtain a model for apricot mass and separate apricot variety with image processing technology using external features of apricot fruit. In this study, five verities of apricot were used. In order to determine the size of the fruits, three mutually perpendicular axes were defined, length, width, and thickness. Measurements show that the effect of variety on all properties was statistically significant at the 1% probability level. Furthermore, there is no significant difference between the estimated dimensions by image processing approach and the actual dimensions. The developed system consists of a digital camera, a light diffusion chamber, a distance adjustment pedestal, and a personal computer. Images taken by the digital camera were stored as (RGB) for further analysis. The images were taken for a number of 49 samples of each cultivar in three directions. A linear equation is recommended to calculate the apricot mass based on the length and the width with R 2 = 0.97. In addition, ANFIS model with C-means was the best model for classifying the apricot varieties based on the physical features including length, width, thickness, mass, and projected area of three perpendicular surfaces. The accuracy of the model was 87.7.


title: Non-Cooperative Game Theory Based Rate Adaptation for Dynamic Video Streaming over HTTP

abstract: Dynamic Adaptive Streaming over HTTP (DASH) has demonstrated to be an emerging and promising multimedia streaming technique, owing to its capability of dealing with the variability of networks. Rate adaptation mechanism, a challenging and open issue, plays an important role in DASH based systems since it affects Quality of Experience (QoE) of users, network utilization, etc. In this paper, based on non-cooperative game theory, we propose a novel algorithm to optimally allocate the limited export bandwidth of the server to multi-users to maximize their QoE with fairness guaranteed. The proposed algorithm is proxy-free. Specifically, a novel user QoE model is derived by taking a variety of factors into account, like the received video quality, the reference buffer length, and user accumulated buffer lengths, etc. Then, the bandwidth competing problem is formulated as a non-cooperation game with the existence of Nash Equilibrium that is theoretically proven. Finally, a distributed iterative algorithm with stability analysis is proposed to find the Nash Equilibrium. Compared with state-of-the-art methods, extensive experimental results in terms of both simulated and realistic networking scenarios demonstrate that the proposed algorithm can produce higher QoE, and the actual buffer lengths of all users keep nearly optimal states, i.e., moving around the reference buffer all the time. Besides, the proposed algorithm produces no playback interruption.


title: Galois Symmetry Induced by Hecke Relations in Rational Conformal Field Theory and Associated Modular Tensor Categories

abstract: Hecke operators relate characters of rational conformal field theories (RCFTs) with different central charges, and extend the previously studied Galois symmetry of modular representations and fusion algebras. We show that the conductor $N$ of a RCFT and the quadratic residues mod $N$ play an important role in the computation and classification of Galois permutations. We establish a field correspondence in different theories through the picture of effective central charge, which combines Galois inner automorphisms and the structure of simple currents. We then make a first attempt to extend Hecke operators to the full data of modular tensor categories. The Galois symmetry encountered in the modular data appears in the fusion and the braiding matrices as well, and yields isomorphic structures in theories related by Hecke operators.


title: Study of Switched Max-Link Relay Selection for Cooperative Multiple-Antenna Systems

abstract: In this work, we present a switched relaying framework for multiple-input multiple-output (MIMO) relay systems where a source node may transmit directly to a destination node or aided by relays. We also investigate relay selection techniques for the proposed switched relaying framework, whose relays are equipped with buffers. In particular, we develop a novel relay selection protocol based on switching and the selection of the best link, denoted as Switched Max-Link. We then propose the Maximum Minimum Distance (MMD) relay selection criterion for MIMO systems, which is based on the optimal Maximum Likelihood (ML) principle and can provide significant performance gains over other criteria, along with algorithms that are incorporated into the proposed Switched Max-Link protocol. An analysis of the proposed Switched Max-Link protocol and the MMD relay selection criterion in terms of computational cost, pairwise error probability, sum-rate and average delay is carried out. Simulations show that Switched Max-Link using the MMD criterion outperforms previous works in terms of sum-rate, pairwise error probability, average delay and bit error rate.