title: Excitonic emission in van-der-Waals nanotubes of transition metal dichalcogenides

abstract: Nanotubes (NTs) of transition metal dichalcogenides (TMDs), such as MoS2 and WS2, were first synthesized more than a quarter of a century ago; nevertheless, many of their properties have so far remained basically unknown. This review presents the state of the art in the knowledge of the optical properties of TMD NTs. We first evaluate general properties of multilayered TMD crystals, and analyze available data on electronic band structure and optical properties of related NTs. Then, the technology for the formation and the structural characteristics of TMD NTs are represented, focusing on the structures synthesized by chemical transport reaction. The core of this work is the presentation of the ability of TMD NTs to emit bright photoluminescence (PL), which has been discovered recently. By means of micro-PL spectroscopy of individual tubes we show that excitonic transitions relevant to both direct and indirect band gaps contribute to the emission spectra of the NTs despite the presence of dozens of monolayers in their walls. We highlight the performance of the tubes as efficient optical resonators, whose confined optical modes strongly affect the emission bands. Finally, a brief conclusion is presented, along with an outlook of the future studies of this novel member of the family of radiative NTs, which have unique potential for different nanophotonics applications.


title: On central complete Bell polynomials

abstract: In this paper, we consider central complete and incomplete Bell polynomials which are generalizations of the recently introduced central Bell polynomials and central analogues for the complete and incomplete Bell polynomials. We investigate some properties and identities for these polynomials. Especially, we give explicit formulas for the central complete and incomplete Bell polynomials related to central factorial numbers of the second kind.


title: Search for pair production of first-generation scalar leptoquarks at $\sqrt{s} =$ 13 TeV

abstract: A search for the pair production of first-generation scalar leptoquarks is performed using proton-proton collision data recorded at 13 TeV center-of-mass energy with the CMS detector at the LHC. The data correspond to an integrated luminosity of 35.9 fb$^{-1}$. The leptoquarks are assumed to decay to a quark, and either an electron or a neutrino with branching fractions $\beta$ and 1$-\beta$, respectively. The search targets the decay final states comprising two electrons, or one electron and large missing transverse momentum, along with two quarks that are detected as hadronic jets. First-generation scalar leptoquarks with masses below 1435 (1270) GeV are excluded for $\beta =$ 1.0 (0.5). These are the most stringent limits on the mass of first-generation scalar leptoquarks to date. The data are also interpreted to set exclusion limits in the context of an $R$-parity violating supersymmetric model, predicting promptly decaying top squarks with a similar dielectron final state.


title: A biconvex optimization for solving semidefinite programs via bilinear factorization

abstract: Many problems in machine learning can be reduced to learning a low-rank positive semidefinite matrix (denoted as $Z$), which encounters semidefinite program (SDP). Existing SDP solvers by classical convex optimization are expensive to solve large-scale problems. Employing the low rank of solution, Burer-Monteiro's method reformulated SDP as a nonconvex problem via the $quadratic$ factorization ($Z$ as $XX^\top$). However, this would lose the structure of problem in optimization. In this paper, we propose to convert SDP into a biconvex problem via the $bilinear$ factorization ($Z$ as $XY^\top$), and while adding the term $\frac{\gamma}{2}||X-Y||_F^2$ to penalize the difference of $X$ and $Y$. Thus, the biconvex structure (w.r.t. $X$ and $Y$) can be exploited naturally in optimization. As a theoretical result, we provide a bound to the penalty parameter $\gamma$ under the assumption of $L$-Lipschitz smoothness and $\sigma $-strongly biconvexity, such that, at stationary points, the proposed bilinear factorization is equivalent to Burer-Monteiro's factorization when the bound is arrived, that is $\gamma>\frac{1}{4}(L-\sigma)_+$. Our proposal opens up a new way to surrogate SDP by biconvex program. Experiments on two SDP-related applications demonstrate that the proposed method is effective as the state-of-the-art.


title: Concerning the Neural Code

abstract: The central problem with understanding brain and mind is the neural code issue: understanding the matter of our brain as basis for the phenomena of our mind. The richness with which our mind represents our environment, the parsimony of genetic data, the tremendous efficiency with which the brain learns from scant sensory input and the creativity with which our mind constructs mental worlds all speak in favor of mind as an emergent phenomenon. This raises the further issue of how the neural code supports these processes of organization. The central point of this communication is that the neural code has the form of structured net fragments that are formed by network self-organization, activate and de-activate on the functional time scale, and spontaneously combine to form larger nets with the same basic structure.


title: Proof of a rational Ramanujan-type series for $1/\pi$. The fastest one in level 3

abstract: Using a modular equation of level $3$ and degree $23$ due to Chan and Liaw, we prove the fastest convergent rational Ramanujan-type series for $1/\pi$ of level $3$.


title: Blocked All-Pairs Shortest Paths Algorithm on Intel Xeon Phi KNL Processor: A Case Study

abstract: Manycores are consolidating in HPC community as a way of improving performance while keeping power efficiency. Knights Landing is the recently released second generation of Intel Xeon Phi architecture. While optimizing applications on CPUs, GPUs and first Xeon Phi's has been largely studied in the last years, the new features in Knights Landing processors require the revision of programming and optimization techniques for these devices. In this work, we selected the Floyd-Warshall algorithm as a representative case study of graph and memory-bound applications. Starting from the default serial version, we show how data, thread and compiler level optimizations help the parallel implementation to reach 338 GFLOPS.


title: PT-symmetry a highly ordered quantum system

abstract: We argue by saying that due to conservation of energy ($\langle H\rangle_n \rightleftharpoons \langle K.E\rangle_n + \langle P.E\rangle_n$) PT-symmetry Hamiltonian $H = p^2 - (ix)^N$ is a highly ordered system. Further, it is found that $\langle K.E\rangle_n \gg \langle P.E\rangle_n$.


title: Logarithmic coefficients problems in families related to starlike and convex functions

abstract: Let $\es$ be the family of analytic and univalent functions $f$ in the unit disk $\D$ with the normalization $f(0)=f'(0)-1=0$, and let $\gamma_n(f)=\gamma_n$ denote the logarithmic coefficients of $f\in {\es}$. In this paper, we study bounds for the logarithmic coefficients for certain subfamilies of univalent functions. Also, we consider the families $\F(c)$ and $\G(\delta)$ of functions $f\in {\es}$ defined by $$ {\rm Re} \left ( 1+\frac{zf''(z)}{f'(z)}\right )>1-\frac{c}{2}\, \mbox{ and } \, {\rm Re} \left ( 1+\frac{zf''(z)}{f'(z)}\right )<1+\frac{\delta}{2},\quad z\in \D $$ for some $c\in(0,3]$ and $\delta\in (0,1]$, respectively. We obtain the sharp upper bound for $|\gamma_n|$ when $n=1,2,3$ and $f$ belongs to the classes $\F(c)$ and $\G(\delta)$, respectively. The paper concludes with the following two conjectures: \begin{itemize} \item If $f\in\F (-1/2)$, then $ \displaystyle |\gamma_n|\le \frac{1}{n}\left(1-\frac{1}{2^{n+1}}\right)$ for $n\ge 1$, and $$ \sum_{n=1}^{\infty}|\gamma_{n}|^{2} \leq \frac{\pi^2}{6}+\frac{1}{4} ~{\rm Li\,}_{2}\left(\frac{1}{4}\right) -{\rm Li\,}_{2}\left(\frac{1}{2}\right), $$ where ${\rm Li}_2(x)$ denotes the dilogarithm function. \item If $f\in \G(\delta)$, then $ \displaystyle |\gamma_n|\,\leq \,\frac{\delta}{2n(n+1)}$ for $n\ge 1$. \end{itemize}


title: Rapidity gaps and ancestry

abstract: The recently discovered correspondence between the distribution of rapidity gaps in electron-nucleus diffractive processes and the statistics of the height of genealogical trees in branching random walks is reviewed. In addition, a new comparison of numerical solutions of exact equations for diffraction on the one hand, and for ancestry on the other hand, both established in the framework of the color dipole model, is presented.


title: From Wigner-Yanase-Dyson conjecture to Carlen-Frank-Lieb conjecture (New title)

abstract: In this paper we study the joint convexity/concavity of the trace functions \[ \Psi_{p,q,s}(A,B)=\text{Tr}(B^{\frac{q}{2}}K^*A^{p}KB^{\frac{q}{2}})^s,~~p,q,s\in \mathbb{R}, \] where $A$ and $B$ are positive definite matrices and $K$ is any fixed invertible matrix. We will give full range of $(p,q,s)\in\mathbb{R}^3$ for $\Psi_{p,q,s}$ to be jointly convex/concave for all $K$. As a consequence, we confirm a conjecture of Carlen, Frank and Lieb. In particular, we confirm a weaker conjecture of Audenaert and Datta and obtain the full range of $(\alpha,z)$ for $\alpha$-$z$ R\'enyi relative entropies to be monotone under completely positive trace preserving maps. We also give simpler proofs of many known results, including the concavity of $\Psi_{p,0,1/p}$ for $0<p<1$ which was first proved by Epstein using complex analysis. The key is to reduce the problem to the joint convexity/concavity of the trace functions \[ \Psi_{p,1-p,1}(A,B)=\text{Tr} K^*A^{p}KB^{1-p},~~-1\le p\le 1, \] using a variational method.


title: DUNet: A deformable network for retinal vessel segmentation

abstract: Automatic segmentation of retinal vessels in fundus images plays an important role in the diagnosis of some diseases such as diabetes and hypertension. In this paper, we propose Deformable U-Net (DUNet), which exploits the retinal vessels' local features with a U-shape architecture, in an end to end manner for retinal vessel segmentation. Inspired by the recently introduced deformable convolutional networks, we integrate the deformable convolution into the proposed network. The DUNet, with upsampling operators to increase the output resolution, is designed to extract context information and enable precise localization by combining low-level feature maps with high-level ones. Furthermore, DUNet captures the retinal vessels at various shapes and scales by adaptively adjusting the receptive fields according to vessels' scales and shapes. Three public datasets DRIVE, STARE and CHASE_DB1 are used to train and test our model. Detailed comparisons between the proposed network and the deformable neural network, U-Net are provided in our study. Results show that more detailed vessels are extracted by DUNet and it exhibits state-of-the-art performance for retinal vessel segmentation with a global accuracy of 0.9697/0.9722/0.9724 and AUC of 0.9856/0.9868/0.9863 on DRIVE, STARE and CHASE_DB1 respectively. Moreover, to show the generalization ability of the DUNet, we used another two retinal vessel data sets, one is named WIDE and the other is a synthetic data set with diverse styles, named SYNTHE, to qualitatively and quantitatively analyzed and compared with other methods. Results indicates that DUNet outperforms other state-of-the-arts.


title: Hermite-Gaussian model for quantum states

abstract: In order to characterize quantum states within the context of information geometry, we propose a generalization of the Gaussian model, which we called the Hermite-Gaussian model. We obtain the Fisher-Rao metric and the scalar curvature for this model, and we show its relation with the one-dimensional quantum harmonic oscillator. Moreover, using this model we characterize some failies of states of the quantum harmonic oscillator. We find that for the eigenstates of the Hamiltonian, mixtures of eigenstates and even or odd superpositions of eienstates the associated Fisher-Rao metrics are diagonal.


title: Logarithmic coefficients of the inverse of univalent functions

abstract: Let $\es$ be the class of analytic and univalent functions in the unit disk $|z|<1$, that have a series of the form $f(z)=z+ \sum_{n=2}^{\infty}a_nz^n$. Let $F$ be the inverse of the function $f\in\es$ with the series expansion %in a disk of radius at least $1/4$ $F(w)=f^{-1}(w)=w+ \sum_{n=2}^{\infty}A_nw^n$ for $|w|<1/4$. The logarithmic inverse coefficients $\Gamma_n$ of $F$ are defined by the formula $\log\left(F(w)/w\right)\,=\,2\sum_{n=1}^{\infty}\Gamma_n(F)w^n$. % In this paper, we determine the logarithmic inverse coefficients bound of $F$ for the class In this paper, we first determine the sharp bound for the absolute value of $\Gamma_n(F)$ when $f$ belongs to $\es$ and for all $n \geq 1$. This result motivates us to carry forward similar problems for some of its important geometric subclasses. In some cases, we have managed to solve this question completely but in some other cases it is difficult to handle for $n\geq 4$. For example, in the case of convex functions $f$, we show that the logarithmic inverse coefficients $\Gamma_n(F)$ of $F$ satisfy the inequality \[ |\Gamma_n(F)|\,\le \, \frac{1}{2n} \mbox{ for } n\geq 1,2,3 \] and the estimates are sharp for the function $l(z)=z/(1-z)$. Although this cannot be true for $n\ge 10$, it is not clear whether this inequality could still be true for $4\leq n\leq 9$.


title: Optimal Rank and Select Queries on Dictionary-Compressed Text

abstract: We study the problem of supporting queries on a string $S$ of length $n$ within a space bounded by the size $\gamma$ of a string attractor for $S$. Recent works showed that random access on $S$ can be supported in optimal $O(\log(n/\gamma)/\log\log n)$ time within $O\left (\gamma\ \rm{polylog}\ n \right)$ space. In this paper, we extend this result to \emph{rank} and \emph{select} queries and provide lower bounds matching our upper bounds on alphabets of polylogarithmic size. Our solutions are given in the form of a space-time trade-off that is more general than the one previously known for grammars and that improves existing bounds on LZ77-compressed text by a $\log\log n$ time-factor in \emph{select} queries. We also provide matching lower and upper bounds for \emph{partial sum} and \emph{predecessor} queries within attractor-bounded space, and extend our lower bounds to encompass navigation of dictionary-compressed tree representations.


title: Tensions between the appearance data of T2K and NOvA

abstract: The long baseline neutrino experiments, T2K and NOvA, have taken significant amount of data in each of the four channels: (a) $\nu_\mu$ disappearance, (b) $\bar\nu_\mu$ disappearance (c) $\nu_e$ appearance and (d) $\bar\nu_e$ appearance. There is a mild tension between the disappearance and the appearance data sets of T2K. A more serious tension exists between the $\nu_e$ appearance data of T2K and the $\nu_e / \bar\nu_e$ appearance data of NOvA. This tension is significant enough that T2K rules out the best-fit point of NOvA at $95\%$ confidence level whereas NOvA rules out T2K best-fit point at $90\%$ confidence level. We explain the reason why these tensions arise. We also do a combined fit of T2K and NOvA data and comment on the results of this fit.


title: Reliable graph-based collaborative ranking

abstract: GRank is a recent graph-based recommendation approach the uses a novel heterogeneous information network to model users' priorities and analyze it to directly infer a recommendation list. Unfortunately, GRank neglects the semantics behind different types of paths in the network and during the process, it may use unreliable paths that are inconsistent with the general idea of similarity in neighborhood collaborative ranking. That negligence undermines the reliability of the recommendation list generated by GRank. This paper seeks to present a novel framework for reliable graph-based collaborative ranking, called ReGRank, that ranks items based on reliable recommendation paths that are in harmony with the semantics behind different approaches in neighborhood collaborative ranking. To our knowledge, ReGRank is the first unified framework for neighborhood collaborative ranking that in addition to traditional user-based collaborative ranking, can also be adapted for preference-based and representative-based collaborative ranking as well. Experimental results show that ReGRank significantly improves the state-of-the art neighborhood and graph-based collaborative ranking algorithms.


title: The distribution of the Lasso: Uniform control over sparse balls and adaptive parameter tuning

abstract: The Lasso is a popular regression method for high-dimensional problems in which the number of parameters $\theta_1,\dots,\theta_N$, is larger than the number $n$ of samples: $N>n$. A useful heuristics relates the statistical properties of the Lasso estimator to that of a simple soft-thresholding denoiser,in a denoising problem in which the parameters $(\theta_i)_{i\le N}$ are observed in Gaussian noise, with a carefully tuned variance. Earlier work confirmed this picture in the limit $n,N\to\infty$, pointwise in the parameters $\theta$, and in the value of the regularization parameter. Here, we consider a standard random design model and prove exponential concentration of its empirical distribution around the prediction provided by the Gaussian denoising model. Crucially, our results are uniform with respect to $\theta$ belonging to $\ell_q$ balls, $q\in [0,1]$, and with respect to the regularization parameter. This allows to derive sharp results for the performances of various data-driven procedures to tune the regularization. Our proofs make use of Gaussian comparison inequalities, and in particular of a version of Gordon's minimax theorem developed by Thrampoulidis, Oymak, and Hassibi, which controls the optimum value of the Lasso optimization problem. Crucially, we prove a stability property of the minimizer in Wasserstein distance, that allows to characterize properties of the minimizer itself.


title: Learning to Defend by Learning to Attack

abstract: Adversarial training provides a principled approach for training robust neural networks. From an optimization perspective, adversarial training is essentially solving a bilevel optimization problem. The leader problem is trying to learn a robust classifier, while the follower problem is trying to generate adversarial samples. Unfortunately, such a bilevel problem is difficult to solve due to its highly complicated structure. This work proposes a new adversarial training method based on a generic learning-to-learn (L2L) framework. Specifically, instead of applying existing hand-designed algorithms for the inner problem, we learn an optimizer, which is parametrized as a convolutional neural network. At the same time, a robust classifier is learned to defense the adversarial attack generated by the learned optimizer. Experiments over CIFAR-10 and CIFAR-100 datasets demonstrate that L2L outperforms existing adversarial training methods in both classification accuracy and computational efficiency. Moreover, our L2L framework can be extended to generative adversarial imitation learning and stabilize the training.


title: Bulk inversion asymmetry induced magnetogyrotropic reflection from quantum wells

abstract: Bulk inversion asymmetry (BIA) of III-V and II-VI semiconductor quantum wells is demonstrated by reflection experiments in magnetic field oriented in the structure plane. The linear in the magnetic field contribution to the reflection coefficients is measured at oblique incidence of $s$ and $p$ polarized light in vicinity of exciton resonances. We demonstrate that this contribution to the reflection is caused by magnetogyrotropy of quantum wells, i.e. by the terms in the optical response which are linear in both the magnetic field strength and light wavevector. Theory of magnetogyrotropic effects in light reflection is developed with account for linear in momentum BIA induced terms in the electron and hole effective Hamiltonians. Theoretical estimates agree with the experimental findings. We have found the electron BIA splitting constant in both GaAs and CdTe based quantum wells is about three times smaller than that for heavy holes.


title: Locality and renormalisation: universal properties and integrals on trees

abstract: The purpose of this paper is to build an algebraic framework suited to regularise branched structures emanating from rooted forests and which encodes the locality principle. This is achieved by means of the universal properties in the locality framework of properly decorated rooted forests. These universal properties are then applied to derive the multivariate regularisation of integrals indexed by rooted forests. We study their renormalisation, along the lines of Kreimer's toy model for Feynman integrals.


title: Learning sparse mixtures of rankings from noisy information

abstract: We study the problem of learning an unknown mixture of $k$ rankings over $n$ elements, given access to noisy samples drawn from the unknown mixture. We consider a range of different noise models, including natural variants of the "heat kernel" noise framework and the Mallows model. For each of these noise models we give an algorithm which, under mild assumptions, learns the unknown mixture to high accuracy and runs in $n^{O(\log k)}$ time. The best previous algorithms for closely related problems have running times which are exponential in $k$.


title: Squeezing of light field in a dissipative Jaynes-Cummings model

abstract: Based on the time-convolutionless master-equation approach, we investigate squeezing of light field in a dissipative Jaynes-Cummings model. The results show that squeezing light can be generated when the atom transits to a ground state from an excited state, and then a collapse-revival phenomenon will occur in the squeezing of light field due to atom-cavity coupling. Enhancing the atom-cavity coupling can increase the frequency of the collapse-revival of squeezing. The stronger the non-Markovian effect is, the more obvious the collapse-revival phenomenon is. The oscillatory frequency of the squeezing is dependents on the resonant frequency of the atom-cavity.


title: Harmonic oscillator chain in noncommutative phase space with rotational symmetry

abstract: We consider a quantum space with rotationally invariant noncommutative algebra of coordinates and momenta. The algebra contains tensors of noncommutativity constructed involving additional coordinates and momenta. In the rotationally invariant noncommutative phase space harmonic oscillator chain is studied. We obtain that noncommutativity affects on the frequencies of the system. In the case of a chain of particles with harmonic oscillator interaction we conclude that because of momentum noncommutativity the spectrum of the center-of-mass of the system is discrete and corresponds to the spectrum of harmonic oscillator.


title: Lattice sine-Gordon model

abstract: We obtain nonperturbative results on the sine-Gordon model using the lattice field technique. In particular, we employ the Fourier accelerated hybrid Monte Carlo algorithm for our studies. We find the critical temperature of the theory based on autocorrelation time, as well as the finite size scaling of the "thickness" observable used in an earlier lattice study by Hasenbusch et al. We study the entropy, which is smooth across all temperatures, supportive of an infinite order transition. This system has a well-known duality with the massive Thirring model, which can play the role of a toy models for Montonen-Olive duality in N=4 super-Yang-Mills theory, since it relates solitons to elementary field excitations. Our research lays a groundwork for such study on the lattice.


title: Sharp worst-case evaluation complexity bounds for arbitrary-order nonconvex optimization with inexpensive constraints

abstract: We provide sharp worst-case evaluation complexity bounds for nonconvex minimization problems with general inexpensive constraints, i.e.\ problems where the cost of evaluating/enforcing of the (possibly nonconvex or even disconnected) constraints, if any, is negligible compared to that of evaluating the objective function. These bounds unify, extend or improve all known upper and lower complexity bounds for unconstrained and convexly-constrained problems. It is shown that, given an accuracy level $\epsilon$, a degree of highest available Lipschitz continuous derivatives $p$ and a desired optimality order $q$ between one and $p$, a conceptual regularization algorithm requires no more than $O(\epsilon^{-\frac{p+1}{p-q+1}})$ evaluations of the objective function and its derivatives to compute a suitably approximate $q$-th order minimizer. With an appropriate choice of the regularization, a similar result also holds if the $p$-th derivative is merely H\"older rather than Lipschitz continuous. We provide an example that shows that the above complexity bound is sharp for unconstrained and a wide class of constrained problems, we also give reasons for the optimality of regularization methods from a worst-case complexity point of view, within a large class of algorithms that use the same derivative information.


title: Regulating spin reversal in dipolar systems by the quadratic Zeeman effect

abstract: A mechanism is advanced suggesting the resolution of the dichotomy of long-lived spin polarization storage versus fast spin reversal at the required time. A system of atoms or molecules is considered interacting through magnetic dipolar forces. The constituents are assumed to possess internal structure allowing for the generation of the alternating-current quadratic Zeeman effect, whose characteristics can be efficiently regulated by quasiresonant dressing. The sample is connected to an electric circuit producing a feedback field acting on spins. By switching on and off the alternating-current quadratic Zeeman effect it is possible to realize spin reversals with a required delay time. The suggested technique of regulated spin reversal can be used in quantum information processing and spintronics.


title: Time-Frequency Audio Features for Speech-Music Classification

abstract: Distinct striation patterns are observed in the spectrograms of speech and music. This motivated us to propose three novel time-frequency features for speech-music classification. These features are extracted in two stages. First, a preset number of prominent spectral peak locations are identified from the spectra of each frame. These important peak locations obtained from each frame are used to form Spectral peak sequences (SPS) for an audio interval. In second stage, these SPS are treated as time series data of frequency locations. The proposed features are extracted as periodicity, average frequency and statistical attributes of these spectral peak sequences. Speech-music categorization is performed by learning binary classifiers on these features. We have experimented with Gaussian mixture models, support vector machine and random forest classifiers. Our proposal is validated on four datasets and benchmarked against three baseline approaches. Experimental results establish the validity of our proposal.


title: Local stress analysis in Cu$_{50\%}$Zr$_{50\%}$ metallic glass under shear strain by means of first principle modeling

abstract: Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this paper using {\em ab initio} molecular dynamics simulations of the Cu$_{50}$Zr$_{50}$ metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the non-affine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.


title: Turing Degrees and Automorphism Groups of Substructure Lattices

abstract: The study of automorphisms of computable and other structures connects computability theory with classical group theory. Among the noncomputable countable structures, computably enumerable structures are one of the most important objects of investigation in computable model theory. In this paper, we focus on the lattice structure of computably enumerable substructures of a given canonical computable structure. In particular, for a Turing degree $\mathbf{d}$, we investigate the groups of $\mathbf{d}$ -computable automorphisms of the lattice of $\mathbf{d}$-computably enumerable vector spaces, of the interval Boolean algebra $\mathcal{B}_{\eta }$ of the ordered set of rationals, and of the lattice of $\mathbf{d}$ -computably enumerable subalgebras of $\mathcal{B}_{\eta }$. For these groups we show that Turing reducibility can be used to substitute the group-theoretic embedding. We also prove that the Turing degree of the isomorphism types for these groups is the second Turing jump of $\mathbf{d}$ , $\mathbf{d^{\prime \prime }}$.


title: CAAD 2018: Powerful None-Access Black-Box Attack Based on Adversarial Transformation Network

abstract: In this paper, we propose an improvement of Adversarial Transformation Networks(ATN) to generate adversarial examples, which can fool white-box models and black-box models with a state of the art performance and won the 2rd place in the non-target task in CAAD 2018.


title: Multidimensional segment trees can do range updates in poly-logarithmic time

abstract: Updating and querying on a range is a classical algorithmic problem with a multitude of applications. The Segment Tree data structure is particularly notable in handling the range query and update operations. A Segment Tree divides the range into disjoint segments and merges them together to perform range queries and range updates elegantly. Although this data structure is remarkably potent for 1-dimensional problems, it falls short in higher dimensions. Lazy Propagation enables the operations to be computed in $O(logn)$ time in a single dimension. However, the concept of lazy propagation could not be translated to higher-dimensional cases, which imposes a time complexity of $O(n^{k-1} \; logn)$ for operations on $k$-dimensional data. In this work, we have made an attempt to emulate the idea of lazy propagation differently so that it can be applied for 2-dimensional cases. Moreover, the proposed modification should be capable of performing most general aggregate functions similar to the original Segment Tree, and can also be extended to even higher dimensions. Our proposed algorithm manages to perform range sum queries and updates in $O(\log^2 n)$ time for a 2-dimensional problem, which becomes $O(\log^d n)$ for a $d$-dimensional situation.


title: Equitable voting rules

abstract: May's Theorem (1952), a celebrated result in social choice, provides the foundation for majority rule. May's crucial assumption of symmetry, often thought of as a procedural equity requirement, is violated by many choice procedures that grant voters identical roles. We show that a weakening of May's symmetry assumption allows for a far richer set of rules that still treat voters equally. We show that such rules can have minimal winning coalitions comprising a vanishing fraction of the population, but not less than the square root of the population size. Methodologically, we introduce techniques from group theory and illustrate their usefulness for the analysis of social choice questions.


title: Generalisation of the Magnetic Field Configuration of typical and atypical Confined Flares

abstract: Atypical flares cannot be naturally explained with standard models. To predict such flares, we need to define their physical characteristics, in particular, their magnetic environment, and identify pairs of reconnected loops. Here, we present in detail a case-study of a confined flare preceded by flux cancellation that leads to the formation of a filament. The slow rise of the non-eruptive filament favours the growth and reconnection of overlying loops. The flare is only of C5.0 class but it is a long duration event. The reason is that it is comprised of three successive stages of reconnection. A non-linear force-free field extrapolation and a magnetic topology analysis allow us to identify the loops involved in the reconnection process and build a reliable scenario for this atypical confined flare. The main result is that a curved magnetic polarity inversion line in active regions is a key ingredient for producing such atypical flares. A comparison with previous extrapolations for typical and atypical confined flares leads us to propose a cartoon for generalizing the concept


title: Farey boat I. Continued fractions and triangulations, modular group and polygon dissections

abstract: We reformulate several known results about continued fractions in combinatorial terms. Among them the theorem of Conway and Coxeter and that of Series, both relating continued fractions and triangulations. More general polygon dissections appear when extending these theorems for elements of the modular group $PSL(2,\mathbb{Z})$. These polygon dissections are interpreted as walks in the Farey tessellation. The combinatorial model of continued fractions can be further developed to obtain a canonical presentation of elements of $PSL(2,\mathbb{Z})$.


title: The Hermite ring conjecture and special linear groups for valuation rings

abstract: In this paper we prove that the Hermite ring conjecture holds for valuation rings $V$, and the special liner group $SL_n(V[x])$ coincides with the group generated by elementary matrices $E_n(V[x])$ for $n\geq3$. For any arithmetical ring $R$ and $n\geq3$, we show $SL_n(R[x]) = SL_n(R)\cdot E_n(R[x])$.


title: Inclusive breakup reaction of a two-cluster projectile on a two-fragment target: A genuine four-body problem

abstract: We develop a four-body model for the inclusive breakup of two-fragment halo projectiles colliding with two-fragment targets. In the case of a short lived projectiles, such as halo nuclei, on a deuteron target, the model allows the extraction of the neutron capture cross section of such projectiles. We supply examples.


title: Designing Dual-Band Absorbers by Graphene/Metallic Metasurfaces

abstract: This article presents a novel approach for designing dual-band absorbers based on graphene and metallic metasurfaces for terahertz and mid-infrared regimes, respectively. The absorbers are composed of a two-dimensional (2D) array of square patches deposited on a dielectric film terminated by a metal plate. Using an analytical circuit model, we obtain closed-form relations for different parameters of the structure to achieve the dual-band absorber. Two absorption bands with obtained absorptivity of 98% at 0.55 and 1.54 THz for the graphene-based structure and 7 and 25 THz for the metallic-based case are achieved. We demonstrate that the graphene-based absorber remains as dual-band for a wide range of the Fermi level. Furthermore, the proposed dual-band absorbers have a polarization-insensitive characteristic remaining over a wide range of incident angles. The most important advantage of this device is its simplicity compared to previously reported structures.


title: Deep Ad-hoc Beamforming

abstract: Far-field speech processing is an important and challenging problem. In this paper, we propose \textit{deep ad-hoc beamforming}, a deep-learning-based multichannel speech enhancement framework based on ad-hoc microphone arrays, to address the problem. It contains three novel components. First, it combines \textit{ad-hoc microphone arrays} with deep-learning-based multichannel speech enhancement, which reduces the probability of the occurrence of far-field acoustic environments significantly. Second, it groups the microphones around the speech source to a local microphone array by a supervised channel selection framework based on deep neural networks. Third, it develops a simple time synchronization framework to synchronize the channels that have different time delay. Besides the above novelties and advantages, the proposed model is also trained in a single-channel fashion, so that it can easily employ new development of speech processing techniques. Its test stage is also flexible in incorporating any number of microphones without retraining or modifying the framework. We have developed many implementations of the proposed framework and conducted an extensive experiment in scenarios where the locations of the speech sources are far-field, random, and blind to the microphones. Results on speech enhancement tasks show that our method outperforms its counterpart that works with linear microphone arrays by a considerable margin in both diffuse noise reverberant environments and point source noise reverberant environments. We have also tested the framework with different handcrafted features. Results show that although designing good features lead to high performance, they do not affect the conclusion on the effectiveness of the proposed framework.


title: On a topology and limits for inductive systems of $C^*$-algebras over partially ordered sets

abstract: Motivated by algebraic quantum field theory and our previous work we study properties of inductive systems of \ $C^*$-algebras over arbitrary partially ordered sets. A partially ordered set can be represented as the union of the family of its maximal upward directed subsets indexed by elements of a certain set. We consider a topology on the set of indices generated by a base of neighbourhoods. Examples of those topologies with different properties are given. An inductive system of $C^*$-algebras and its inductive limit arise naturally over each maximal upward directed subset. Using those inductive limits, we construct different types of $C^*$-algebras. In particular, for neighbourhoods of the topology on the set of indices we deal with the $C^*$-algebras which are the direct products of those inductive limits. The present paper is concerned with the above-mentioned topology and the algebras arising from an inductive system of $C^*$-algebras over a partially ordered set. We show that there exists a connection between properties of that topology and those $C^*$-algebras.


title: Hardness of computing and approximating predicates and functions with leaderless population protocols

abstract: Population protocols are a distributed computing model appropriate for describing massive numbers of agents with limited computational power. A population protocol "has an initial leader" if every valid initial configuration contains a single agent in a special "leader" state that helps to coordinate the computation. Although the class of predicates and functions computable with probability 1 is the same whether or not there is an initial leader (semilinear functions and predicates), it is not known whether a leader is necessary for fast computation. Efficient population protocols are generally defined as those computing in polylogarithmic in $n$ (parallel) time. We consider leaderless population protocols, regarding the computation finished when a configuration is reached from which a different output is no longer reachable. In this setting we show that a wide class of functions and predicates computable by population protocols are not efficiently computable (they require at least linear time to stabilize on a correct answer), nor are some linear functions even efficiently approximable. For example, the widely studied parity, majority, and equality predicates cannot be computed in sublinear time. Moreover, it requires at least linear time for a population protocol even to approximate any linear function with a coefficient outside of $\mathbb{N}$: for sufficiently small $\gamma > 0$, the output of a sublinear time protocol can stabilize outside the interval $f(m) (1 \pm \gamma)$ on infinitely many inputs $m$. We also show that it requires linear time to exactly compute a wide range of semilinear functions (e.g., $f(m)=m$ if $m$ is even and $2m$ if $m$ is odd). Finally, we show that with a sufficiently large value of $\gamma$, a population protocol can approximate any linear $f$ with nonnegative rational coefficients, within approximation factor $\gamma$, in $O(\log n)$ time.


title: Higher-order topological superconductivity: possible realization in Fermi gases and Sr$_2$RuO$_4$

abstract: We propose to realize second-order topological superconductivity in bilayer spin-polarized Fermi gas superfluids. We focus on systems with intralayer chiral $p$-wave pairing and with tunable interlayer hopping and interlayer interactions. Under appropriate circumstances, an interlayer even-parity $s$- or $d$-wave pairing may coexist with the intralayer $p$-wave. Our model supports localized Majorana zero modes not only at the corners of the system geometry, but also at the terminations of certain one-dimensional defects, such as lattice line defects and superfluid domain walls. We show how such topological phases and the Majorana zero modes therein can be manipulated in a multitude of ways by tuning the interlayer pairing and hopping. Generalized to spinful systems, we further propose that the putative $p$-wave superconductor Sr$_{2}$RuO$_{4}$, when subject to uniaxial strains, may also realize the desired topological phase.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.quant-gas cond-mat.str-el quant-ph


title: Relation Mention Extraction from Noisy Data with Hierarchical Reinforcement Learning

abstract: In this paper we address a task of relation mention extraction from noisy data: extracting representative phrases for a particular relation from noisy sentences that are collected via distant supervision. Despite its significance and value in many downstream applications, this task is less studied on noisy data. The major challenges exists in 1) the lack of annotation on mention phrases, and more severely, 2) handling noisy sentences which do not express a relation at all. To address the two challenges, we formulate the task as a semi-Markov decision process and propose a novel hierarchical reinforcement learning model. Our model consists of a top-level sentence selector to remove noisy sentences, a low-level mention extractor to extract relation mentions, and a reward estimator to provide signals to guide data denoising and mention extraction without explicit annotations. Experimental results show that our model is effective to extract relation mentions from noisy data.


title: Deep Learning based Computer-Aided Diagnosis Systems for Diabetic Retinopathy: A Survey

abstract: Diabetic retinopathy (DR) results in vision loss if not treated early. A computer-aided diagnosis (CAD) system based on retinal fundus images is an efficient and effective method for early DR diagnosis and assisting experts. A computer-aided diagnosis (CAD) system involves various stages like detection, segmentation and classification of lesions in fundus images. Many traditional machine-learning (ML) techniques based on hand-engineered features have been introduced. The recent emergence of deep learning (DL) and its decisive victory over traditional ML methods for various applications motivated the researchers to employ it for DR diagnosis, and many deep-learning-based methods have been introduced. In this paper, we review these methods, highlighting their pros and cons. In addition, we point out the challenges to be addressed in designing and learning about efficient, effective and robust deep-learning algorithms for various problems in DR diagnosis and draw attention to directions for future research.


title: Modern Data Analytics Approach to Predict Creep of High-Temperature Alloys

abstract: A breakthrough in alloy design often requires comprehensive understanding in complex multi-component/multi-phase systems to generate novel material hypotheses. We introduce a modern data analytics workflow that leverages high-quality experimental data augmented with advanced features obtained from high-fidelity models. Herein, we use an example of a consistently-measured creep dataset of developmental high-temperature alloy combined with scientific alloy features populated from a high-throughput computational thermodynamic approach. Extensive correlation analyses provide ranking insights for most impactful alloy features for creep resistance, evaluated from a large set of candidate features suggested by domain experts. We also show that we can accurately train machine learning models by integrating high-ranking features obtained from correlation analyses. The demonstrated approach can be extended beyond incorporating thermodynamic features, with input from domain experts used to compile lists of features from other alloy physics, such as diffusion kinetics and microstructure evolution.


title: Semiclassical sampling and discretization of certain linear inverse problems

abstract: We study sampling of Fourier Integral Operators $A$ at rates $sh$ with $s$ fixed and $h$ a small parameter. We show that the Nyquist sampling limit of $Af$ and $f$ are related by the canonical relation of $A$ using semiclassical analysis. We apply this analysis to the Radon transform in the parallel and the fan-beam coordinates. We explain and illustrate the optimal sampling rates for $Af$, the aliasing artifacts, and the effect of averaging (blurring) the data $Af$. We prove a Weyl type of estimate on the minimal number of sampling points to recover $f$ stably in terms of the volume of its semiclassical wave front set.


title: Wizard of Wikipedia: Knowledge-Powered Conversational agents

abstract: In open-domain dialogue intelligent agents should exhibit the use of knowledge, however there are few convincing demonstrations of this to date. The most popular sequence to sequence models typically "generate and hope" generic utterances that can be memorized in the weights of the model when mapping from input utterance(s) to output, rather than employing recalled knowledge as context. Use of knowledge has so far proved difficult, in part because of the lack of a supervised learning benchmark task which exhibits knowledgeable open dialogue with clear grounding. To that end we collect and release a large dataset with conversations directly grounded with knowledge retrieved from Wikipedia. We then design architectures capable of retrieving knowledge, reading and conditioning on it, and finally generating natural responses. Our best performing dialogue models are able to conduct knowledgeable discussions on open-domain topics as evaluated by automatic metrics and human evaluations, while our new benchmark allows for measuring further improvements in this important research direction.


title: Unfolding with Gaussian Processes

abstract: A method to perform unfolding with Gaussian processes (GPs) is presented. Using Bayesian regression, we define an estimator for the underlying truth distribution as the mode of the posterior. We show that in the case where the bin contents are distributed approximately according to a Gaussian, this estimator is equivalent to the mean function of a GP conditioned on the maximum likelihood estimator. Regularisation is introduced via the kernel function of the GP, which has a natural interpretation as the covariance of the underlying distribution. This novel approach allows for the regularisation to be informed by prior knowledge of the underlying distribution, and for it to be varied along the spectrum. In addition, the full statistical covariance matrix for the estimator is obtained as part of the result. The method is applied to two examples: a double-peaked bimodal distribution and a falling spectrum.


title: Sparse domination and the strong maximal function

abstract: We study the problem of dominating the dyadic strong maximal function by $(1, 1)$-type sparse forms based on rectangles with sides parallel to the axes, and show that such domination is impossible. Our proof relies on an explicit construction of a pair of maximally separated point sets with respect to an appropriately defined notion of distance


title: Regularity and stability analysis for a class of semilinear nonlocal differential equations in Hilbert spaces

abstract: We deal with a class of semilinear nonlocal differential equations in Hilbert spaces which is a general model for some anomalous diffusion equations. By using the theory of integral equations with completely positive kernel together with local estimates, some existence, regularity and stability results are established. An application to nonlocal partial differential equations is shown to demonstrate our abstract results.


title: Dipole of the luminosity distance as a test for dark energy models

abstract: The dependence of Hubble parameter on redshift can be determined directly from the dipole of luminosity distance to Supernovae Ia. We investigate the possibility of using the data on dipole of the luminosity distance obtained from the Supernovae Ia compilations SDSS, Union2.1, JLA and Pantheon to distinguish the dark energy models.


title: Maximal estimates for a generalized spherical mean Radon transform acting on radial functions

abstract: We study a generalized spherical means operator, viz. generalized spherical mean Radon transform, acting on radial functions. As the main results, we find conditions for the associated maximal operator and its local variant to be bounded on power weighted Lebesgue spaces. This translates, in particular, into almost everywhere convergence to radial initial data results for solutions to certain Cauchy problems for classical Euler-Poisson-Darboux and wave equations. Moreover, our results shed some new light to the interesting and important question of optimality of the yet known $L^p$ boundedness results for the maximal operator in the general non-radial case. It appears that these could still be notably improved, as indicated by our conjecture of the ultimate sharp result.


title: Stochastic Neighbor Embedding under f-divergences

abstract: The t-distributed Stochastic Neighbor Embedding (t-SNE) is a powerful and popular method for visualizing high-dimensional data. It minimizes the Kullback-Leibler (KL) divergence between the original and embedded data distributions. In this work, we propose extending this method to other f-divergences. We analytically and empirically evaluate the types of latent structure-manifold, cluster, and hierarchical-that are well-captured using both the original KL-divergence as well as the proposed f-divergence generalization, and find that different divergences perform better for different types of structure. A common concern with $t$-SNE criterion is that it is optimized using gradient descent, and can become stuck in poor local minima. We propose optimizing the f-divergence based loss criteria by minimizing a variational bound. This typically performs better than optimizing the primal form, and our experiments show that it can improve upon the embedding results obtained from the original $t$-SNE criterion as well.


title: Compressed Multiple Pattern Matching

abstract: Given $d$ strings over the alphabet $\{0,1,\ldots,\sigma{-}1\}$, the classical Aho--Corasick data structure allows us to find all $occ$ occurrences of the strings in any text $T$ in $O(|T| + occ)$ time using $O(m\log m)$ bits of space, where $m$ is the number of edges in the trie containing the strings. Fix any constant $\varepsilon \in (0, 2)$. We describe a compressed solution for the problem that, provided $\sigma \le m^\delta$ for a constant $\delta < 1$, works in $O(|T| \frac{1}{\varepsilon} \log\frac{1}{\varepsilon} + occ)$ time, which is $O(|T| + occ)$ since $\varepsilon$ is constant, and occupies $mH_k + 1.443 m + \varepsilon m + O(d\log\frac{m}{d})$ bits of space, for all $0 \le k \le \max\{0,\alpha\log_\sigma m - 2\}$ simultaneously, where $\alpha \in (0,1)$ is an arbitrary constant and $H_k$ is the $k$th-order empirical entropy of the trie. Hence, we reduce the $3.443m$ term in the space bounds of previously best succinct solutions to $(1.443 + \varepsilon)m$, thus solving an open problem posed by Belazzougui. Further, we notice that $L = \log\binom{\sigma (m+1)}{m} - O(\log(\sigma m))$ is a worst-case space lower bound for any solution of the problem and, for $d = o(m)$ and constant $\varepsilon$, our approach allows to achieve $L + \varepsilon m$ bits of space, which gives an evidence that, for $d = o(m)$, the space of our data structure is theoretically optimal up to the $\varepsilon m$ additive term and it is hardly possible to eliminate the term $1.443m$. In addition, we refine the space analysis of previous works by proposing a more appropriate definition for $H_k$. We also simplify the construction for practice adapting the fixed block compression boosting technique, then implement our data structure, and conduct a number of experiments showing that it is comparable to the state of the art in terms of time and is superior in space.


title: Dynamic Feature Acquisition Using Denoising Autoencoders

abstract: In real-world scenarios, different features have different acquisition costs at test-time which necessitates cost-aware methods to optimize the cost and performance trade-off. This paper introduces a novel and scalable approach for cost-aware feature acquisition at test-time. The method incrementally asks for features based on the available context that are known feature values. The proposed method is based on sensitivity analysis in neural networks and density estimation using denoising autoencoders with binary representation layers. In the proposed architecture, a denoising autoencoder is used to handle unknown features (i.e., features that are yet to be acquired), and the sensitivity of predictions with respect to each unknown feature is used as a context-dependent measure of informativeness. We evaluated the proposed method on eight different real-world datasets as well as one synthesized dataset and compared its performance with several other approaches in the literature. According to the results, the suggested method is capable of efficiently acquiring features at test-time in a cost- and context-aware fashion.


title: Global climate modeling of Saturn's atmosphere. Part II: multi-annual high-resolution dynamical simulations

abstract: The Cassini mission unveiled the intense and diverse activity in Saturn's atmosphere: banded jets, waves, vortices, equatorial oscillations. To set the path towards a better understanding of those phenomena, we performed high-resolution multi-annual numerical simulations of Saturn's atmospheric dynamics. We built a new Global Climate Model [GCM] for Saturn, named the Saturn DYNAMICO GCM, by combining a radiative-seasonal model tailored for Saturn to a hydrodynamical solver based on an icosahedral grid suitable for massively-parallel architectures. The impact of numerical dissipation, and the conservation of angular momentum, are examined in the model before a reference simulation employing the Saturn DYNAMICO GCM with a $1/2^{\circ}$ latitude-longitude resolution is considered for analysis. Mid-latitude banded jets showing similarity with observations are reproduced by our model. Those jets are accelerated and maintained by eddy momentum transfers to the mean flow, with the magnitude of momentum fluxes compliant with the observed values. The eddy activity is not regularly distributed with time, but appears as bursts; both barotropic and baroclinic instabilities could play a role in the eddy activity. The steady-state latitude of occurrence of jets is controlled by poleward migration during the spin-up of our model. At the equator, a weakly-superrotating tropospheric jet and vertically-stacked alternating stratospheric jets are obtained in our GCM simulations. The model produces Yanai (Rossby-gravity), Rossby and Kelvin waves at the equator, as well as extratropical Rossby waves, and large-scale vortices in polar regions. Challenges remain to reproduce Saturn's powerful superrotating jet and hexagon-shaped circumpolar jet in the troposphere, and downward-propagating equatorial oscillation in the stratosphere.


title: Multi-View Networks For Multi-Channel Audio Classification

abstract: In this paper we introduce the idea of multi-view networks for sound classification with multiple sensors. We show how one can build a multi-channel sound recognition model trained on a fixed number of channels, and deploy it to scenarios with arbitrary (and potentially dynamically changing) number of input channels and not observe degradation in performance. We demonstrate that at inference time you can safely provide this model all available channels as it can ignore noisy information and leverage new information better than standard baseline approaches. The model is evaluated in both an anechoic environment and in rooms generated by a room acoustics simulator. We demonstrate that this model can generalize to unseen numbers of channels as well as unseen room geometries.


title: Jordanable almost Abelian Lie algebras

abstract: We call a linear operator on a vector space over a field Jordanable if it has a Jordan canonical form. An almost Abelian Lie algebra has only one non-vanishing Lie bracket, which is given by a linear operator. If the latter is Jordanable then we call the almost Abelian Lie algebra Jordanable. The existence of canonical form allows to describe important structural properties of these Lie algebras in explicit terms. Lie subalgebras, ideals, automorphisms and derivations, as well as quadratic Casimir elements are given in explicitly in a suitable basis.


title: Systematic study of outflows in the Local Universe using CALIFA: I. Sample selection and main properties

abstract: We present a sample of 17 objects from the CALIFA survey where we find initial evidence of galactic winds based on their off-axis ionization properties. We identify the presence of outflows using various optical diagnostic diagrams (e.g., EW(H$\alpha$), [Nii]/H$\alpha$, [Sii]/H$\alpha$, [Oi]/H$\alpha$ line-ratio maps). We find that all 17 candidate outflow galaxies lie along the sequence of active star formation in the M$_\star$ vs. star-formation rate diagram, without a clear excess in the integrated SFR. The location of galaxies along the star-formation main sequence (SFMS) does not influence strongly the presence or not of outflows. The analysis of the star-formation rate density ($\Sigma_{\rm SFR}$) reveals that the CALIFA sources present higher values when compared with normal star-forming galaxies. The strength of this relation depends on the calibrator used to estimate the SFR. This excess in $\Sigma_{\rm SFR}$ is significant within the first effective radius supporting the idea that most outflows are driven by processes in the inner regions of a galaxy. We find that the molecular gas mass density ($\Sigma_\mathrm{gas}$) is a key parameter that plays an important role in the generation of outflows through its association with the local SFR. The canonical threshold reported for the generation of outflows -- $\Sigma_{\rm SFR}>0.1$ $\mathrm{M}_\odot \mathrm{yr}^{-1} \mathrm{kpc}^{-2}$ -- is only marginally exceeded in our sample. Within the Kennicutt-Schmidt diagram we propose a domain for galaxies hosting starburst-driven outflows defined by $\Sigma_{\rm SFR}>10^{-2} \,\mathrm{M}_\odot \mathrm{yr}^{-1} \mathrm{kpc}^{-2}$ and $\Sigma_\mathrm{gas}>10^{1.2} \, \mathrm{M}_\odot \mathrm{pc}^{-2}$ within a central kiloparcec region.


title: A Factor Graph Approach to Multi-Camera Extrinsic Calibration on Legged Robots

abstract: Legged robots are becoming popular not only in research, but also in industry, where they can demonstrate their superiority over wheeled machines in a variety of applications. Either when acting as mobile manipulators or just as all-terrain ground vehicles, these machines need to precisely track the desired base and end-effector trajectories, perform Simultaneous Localization and Mapping (SLAM), and move in challenging environments, all while keeping balance. A crucial aspect for these tasks is that all onboard sensors must be properly calibrated and synchronized to provide consistent signals for all the software modules they feed. In this paper, we focus on the problem of calibrating the relative pose between a set of cameras and the base link of a quadruped robot. This pose is fundamental to successfully perform sensor fusion, state estimation, mapping, and any other task requiring visual feedback. To solve this problem, we propose an approach based on factor graphs that jointly optimizes the mutual position of the cameras and the robot base using kinematics and fiducial markers. We also quantitatively compare its performance with other state-of-the-art methods on the hydraulic quadruped robot HyQ. The proposed approach is simple, modular, and independent from external devices other than the fiducial marker.


title: Adjoint method and inverse design for nonlinear nanophotonic devices

abstract: The development of inverse design, where computational optimization techniques are used to design devices based on certain specifications, has led to the discovery of many compact, non-intuitive structures with superior performance. Among various methods, large-scale, gradient-based optimization techniques have been one of the most important ways to design a structure containing a vast number of degrees of freedom. These techniques are made possible by the adjoint method, in which the gradient of an objective function with respect to all design degrees of freedom can be computed using only two full-field simulations. However, this approach has so far mostly been applied to linear photonic devices. Here, we present an extension of this method to modeling nonlinear devices in the frequency domain, with the nonlinear response directly included in the gradient computation. As illustrations, we use the method to devise compact photonic switches in a Kerr nonlinear material, in which low-power and high-power pulses are routed in different directions. Our technique may lead to the development of novel compact nonlinear photonic devices.


title: Automaticity and invariant measures of linear cellular automata

abstract: We show that spacetime diagrams of linear cellular automata $\Phi : {\mathbb F}_p^{\mathbb Z} \to {\mathbb F}_p^{\mathbb Z}$ with $(-p)$-automatic initial conditions are automatic. This extends existing results on initial conditions which are eventually constant. Each automatic spacetime diagram defines a $(\sigma, \Phi)$-invariant subset of ${\mathbb F}_p^{\mathbb Z}$, where $\sigma$ is the left shift map, and if the initial condition is not eventually periodic then this invariant set is nontrivial. For the Ledrappier cellular automaton we construct a family of nontrivial $(\sigma, \Phi)$-invariant measures on ${\mathbb F}_3^{\mathbb Z}$. Finally, given a linear cellular automaton $\Phi$, we construct a nontrivial $(\sigma, \Phi)$-invariant measure on ${\mathbb F}_p^{\mathbb Z}$ for all but finitely many $p$.


title: Recovery of compressively sensed ultrasound images with structured Sparse Bayesian Learning

abstract: In this paper, we consider the problem of recovering compressively sensed ultrasound images. We build on prior work, and consider a number of existing approaches that we consider to be the state-of-the-art. The methods we consider take advantage of a number of assumptions on the signals including those of temporal and spatial correlation, block structure, prior knowledge of the support, and non-Gaussianity. We conduct a series of intensive tests to quantify the performance of these methods. We find that by altering the parameters of the structured Sparse Bayesian Learning approaches considered, we can significantly improve the objective quality of the reconstructed images. The results we achieve are a significant improvement upon previously proposed reconstruction techniques. In addition, we further show that by careful choice of parameters, we can obtain near-optimal results whilst requiring only a small fraction of the computational time needed for the best reconstruction quality.


title: Aperiodic-Order-Induced Enhancement of Weak Nonlocality in Multilayered Dielectric Metamaterials

abstract: Recent studies on fully dielectric multilayered metamaterials have shown that the negligibly small nonlocal effects (spatial dispersion) typically observed in the limit of deeply subwavelength layers may be significantly enhanced by peculiar boundary effects occurring in certain critical parameter regimes. These phenomena, observed so far in periodic and randomly disordered geometries, are manifested as strong differences between the exact optical response of finite-size metamaterial samples and the prediction from conventional effective-theory-medium models based on mixing formulae. Here, with specific focus on the Thue-Morse geometry, we make a first step toward extending the studies above to the middle-ground of aperiodically ordered multilayers, lying in between perfect periodicity and disorder. We show that, also for these geometries, there exist critical parameter ranges that favor the buildup of boundary effects leading to strong enhancement of the (otherwise negligibly weak) nonlocality. However, the underlying mechanisms are fundamentally different from those observed in the periodic case, and exhibit typical footprints (e.g., fractal gaps, quasi-localized states) that are distinctive of aperiodic order. The outcomes of our study indicate that aperiodic order plays a key role in the buildup of the aforementioned boundary effects, and may also find potential applications to optical sensors, absorbers and lasers.


title: QuickXsort - A Fast Sorting Scheme in Theory and Practice

abstract: QuickXsort is a highly efficient in-place sequential sorting scheme that mixes Hoare's Quicksort algorithm with X, where X can be chosen from a wider range of other known sorting algorithms, like Heapsort, Insertionsort and Mergesort. Its major advantage is that QuickXsort can be in-place even if X is not. In this work we provide general transfer theorems expressing the number of comparisons of QuickXsort in terms of the number of comparisons of X. More specifically, if pivots are chosen as medians of (not too fast) growing size samples, the average number of comparisons of QuickXsort and X differ only by $o(n)$-terms. For median-of-$k$ pivot selection for some constant $k$, the difference is a linear term whose coefficient we compute precisely. For instance, median-of-three QuickMergesort uses at most $n \lg n - 0.8358n + O(\log n)$ comparisons. Furthermore, we examine the possibility of sorting base cases with some other algorithm using even less comparisons. By doing so the average-case number of comparisons can be reduced down to $n \lg n- 1.4106n + o(n)$ for a remaining gap of only $0.0321n$ comparisons to the known lower bound (while using only $O(\log n)$ additional space and $O(n \log n)$ time overall). Implementations of these sorting strategies show that the algorithms challenge well-established library implementations like Musser's Introsort.


title: Black Hole Parameter Estimation from Its Shadow

abstract: The Event Horizon Telescope (EHT), a global submillimeter wavelength very long baseline interferometry array, unveiled event-horizon-scale images of the supermassive black hole M87* as an asymmetric bright emission ring with a diameter of $ 42 \pm 3\; \mu$as, and it is consistent with the shadow of a Kerr black hole of general relativity. A Kerr black hole is also a solution of some alternative theories of gravity, while several modified theories of gravity admit non-Kerr black holes. While earlier estimates for the M87* black hole mass, depending on the method used, fall in the range $ \approx 3\times 10^9 M_\odot- 7 \times 10^9 $$ M_\odot $, the EHT data indicated a mass for the M87* black hole of $(6.5 \pm 0.7) \times 10^9 M_\odot $. This offers another promising tool to estimate black hole parameters and to probe theories of gravity in its most extreme region near the event horizon. The important question arises: Is it possible by a simple technique to estimate black hole parameters from its shadow, for arbitrary models? In this paper, we present observables, expressed in terms of ordinary integrals, characterizing a haphazard shadow shape to estimate the parameters associated with black holes, and then illustrate its relevance to four different models: Kerr, Kerr$-$Newman, and two rotating regular models. Our method is robust, accurate, and consistent with the results obtained from existing formalism, and it is applicable to more general shadow shapes that may not be circular due to noisy data.


title: Canonical Least Favorable Submodels:A New TMLE Procedure for Multidimensional Parameters

abstract: This paper is a fundamental addition to the world of targeted maximum likelihood estimation (TMLE) (or likewise, targeted minimum loss estimation) for simultaneous estimation of multi-dimensional parameters of interest. TMLE, as part of the targeted learning framework, offers a crucial step in constructing efficient plug-in estimators for nonparametric or semiparametric models. The so-called targeting step of targeted learning, involves fluctuating the initial fit of the model in a way that maximally adjusts the plug-in estimate per change in the log likelihood. Previously for multidimensional parameters of interest, iterative TMLE's were constructed using locally least favorable submodels as defined in van der Laan and Gruber, 2016, which are indexed by a multidimensional fluctuation parameter. In this paper we define a canonical least favorable submodel in terms of a single dimensional epsilon for a $d$-dimensional parameter of interest. One can view the clfm as the iterative analog to the one-step TMLE as constructed in van der Laan and Gruber, 2016. It is currently implemented in several software packages we provide in the last section. Using a single epsilon for the targeting step in TMLE could be useful for high dimensional parameters, where using a fluctuation parameter of the same dimension as the parameter of interest could suffer the consequences of curse of dimensionality. The clfm also enables placing the so-called clever covariate denominator as an inverse weight in an offset intercept model. It has been shown that such weighting mitigates the effect of large inverse weights sometimes caused by near positivity violations.


title: Active microrheology to determine viscoelastic parameters of Stokes-Oldroyd B fluids using optical tweezers

abstract: We use active microrheology to determine the frequency dependent moduli of a linear viscoelastic fluid in terms of the polymer time constant ($\lambda$), and the polymer ($\mu_p$) and solvent viscosity ($\mu_s$), respectively. We measure these parameters from the response function of an optically trapped Brownian probe in the fluid under an external perturbation, and at different dilutions of the viscoelastic component in the fluid. This is an improvement over bulk microrheology measurements in viscoelastic Stokes-Oldroyd B fluids which determine the complex elastic modulus $G(\omega)$ of the fluid, but do not, however, reveal the characteristics of the polymer chains and the Newtonian solvent of the complex fluid individually. In a recent work [Paul \textit{et al}., 2018 J. Phys. Condens. Matter \textbf{30} 345101], we linearized the Stokes-Oldroyd B fluid model and thereby explicitly formulated the frequency dependent moduli in terms of ($\mu_p$) and ($\mu_s$), which we now extend to account for an external sinusoidal force applied to the probe particle. We measure $\lambda$, $\mu_p$, and $\mu_s$ experimentally, and compare with existing the $\lambda$ values in the literature for the same fluid at some of the dilution levels, and obtain good agreement. Further, we use these parameters to calculate the complex elastic modulus of the fluid again at certain dilutions and verify successfully with existing data. This establishes our method as an alternate approach in the active microrheology of complex fluids which should reveal information about the composition of such fluids in significantly greater detail and high signal to noise.


title: Encoding-side-channel-free and measurement-device-independent quantum key distribution

abstract: We present a simple protocol where Alice and Bob only needs sending out a coherent state or not-sending out a coherent state to Charlie. There is no bases switching. We show that this protocol is both encoding-state-side-channel free to the source part and measurement-device-independent. We don't have to control exactly the whole space state of the light pulse, which is an impossible task in practice. The protocol is immune to all adverse due to encoding-state imperfections in side-channel space such as the photon frequency spectrum, emission time, propagation direction, spatial angular moment, and so on. Numerical simulation shows that our scheme can reach a side-channel-free result for quantum key distribution over a distance longer than 200 km given the single-photon-interference misalignment error rate of $20\%$.


title: Monolithic mixed-dimensional multigrid methods for single-phase flow in fractured porous media

abstract: This paper deals with the efficient numerical solution of single-phase flow problems in fractured porous media. A monolithic multigrid method is proposed for solving two-dimensional arbitrary fracture networks with vertical and/or horizontal possibly intersecting fractures. The key point is to combine two-dimensional multigrid components (smoother and inter-grid transfer operators) in the porous matrix with their one-dimensional counterparts within the fractures, giving rise to a mixed-dimensional multigrid method. This combination seems to be optimal since it provides an algorithm whose convergence matches the multigrid convergence factor for solving the Darcy problem. Several numerical experiments are presented to demonstrate the robustness of the monolithic mixed-dimensional multigrid method with respect to the permeability of the fractures, the grid size and the number of fractures in the network.


title: Lipschitz free $p$-spaces for $0<p<1$

abstract: This paper initiates the study of the structure of a new class of $p$-Banach spaces, $0<p<1$, namely the Lipschitz free $p$-spaces (alternatively called Arens-Eells $p$-spaces) $\mathcal{F}_{p}(\mathcal{M})$ over $p$-metric spaces. We systematically develop the theory and show that some results hold as in the case of $p=1$, while some new interesting phenomena appear in the case $0<p<1$ which have no analogue in the classical setting. For the former, we, e.g., show that the Lipschitz free $p$-space over a separable ultrametric space is isomorphic to $\ell_{p}$ for all $0<p\le 1$, or that $\ell_p$ isomorphically embeds into $\mathcal{F}_p(\mathcal{M})$ for any $p$-metric space $\mathcal{M}$. On the other hand, solving a problem by the first author and N. Kalton, there are metric spaces $\mathcal{N}\subset \mathcal{M}$ such that the natural embedding from $\mathcal{F}_p(\mathcal{N})$ to $\mathcal{F}_p(\mathcal{M})$ is not an isometry.


title: Lyman-$\alpha$ photons through rotating outflows

abstract: Outflows and rotation are two ubiquitous kinematic features in the gas kinematics of galaxies. Here we introduce a semi-analytic model to quantify how rotating outflows impact the morphology of the Lyman-$\alpha$ emission line. The model is contrasted against Monte Carlo radiative transfer simulations of outflowing gas with additional solid body rotation. We explore a range of neutral Hydrogen optical depth of $10^5 \leq \tau_{\mathrm{H}} \leq10^7$, rotational velocity $0 \leq v_{\mathrm{rot}}/\mathrm{km\ s}^{-1} \leq 100$ and outflow velocity $0\leq v_{\mathrm{out}}/\mathrm{km\ s}^{-1} \leq 50$. We find three consequences of rotation. First, it introduces a dependency with viewing angle; second it broadens the line and third it increases the flux at the line's center. For all viewing angles, the semi-analytic model reproduces the radiative transfer results for the line width and flux change at the line's center within a $7\%$ and $50\%$ precision for an optical depth of $\tau_{\mathrm{H}}=10^5$, respectively, and within $2\%$ and $1\%$ for an optical depth of $\tau_{\mathrm{H}}=10^7$. Using this model we also show that the peaks of integrated spectra taken from opposite sides of an edge-on rotating gas distribution should have a separation of $\frac{1}{2}v_{\mathrm{rot}}$. The semi-analytic model presented here is a convenient tool to introduce rotational kinematics as a post-processing step of idealized Monte Carlo simulations; it provides a framework to interpret \lya spectra in systems where rotation is expected or directly measured through kinematic maps.


title: Legible Normativity for AI Alignment: The Value of Silly Rules

abstract: It has become commonplace to assert that autonomous agents will have to be built to follow human rules of behavior--social norms and laws. But human laws and norms are complex and culturally varied systems, in many cases agents will have to learn the rules. This requires autonomous agents to have models of how human rule systems work so that they can make reliable predictions about rules. In this paper we contribute to the building of such models by analyzing an overlooked distinction between important rules and what we call silly rules--rules with no discernible direct impact on welfare. We show that silly rules render a normative system both more robust and more adaptable in response to shocks to perceived stability. They make normativity more legible for humans, and can increase legibility for AI systems as well. For AI systems to integrate into human normative systems, we suggest, it may be important for them to have models that include representations of silly rules.


title: ReXCam: Resource-Efficient, Cross-Camera Video Analytics at Scale

abstract: Enterprises are increasingly deploying large camera networks for video analytics. Many target applications entail a common problem template: searching for and tracking an object or activity of interest (e.g. a speeding vehicle, a break-in) through a large camera network in live video. Such cross-camera analytics is compute and data intensive, with cost growing with the number of cameras and time. To address this cost challenge, we present ReXCam, a new system for efficient cross-camera video analytics. ReXCam exploits spatial and temporal locality in the dynamics of real camera networks to guide its inference-time search for a query identity. In an offline profiling phase, ReXCam builds a cross-camera correlation model that encodes the locality observed in historical traffic patterns. At inference time, ReXCam applies this model to filter frames that are not spatially and temporally correlated with the query identity's current position. In the cases of occasional missed detections, ReXCam performs a fast-replay search on recently filtered video frames, enabling gracefully recovery. Together, these techniques allow ReXCam to reduce compute workload by 8.3x on an 8-camera dataset, and by 23x - 38x on a simulated 130-camera dataset. ReXCam has been implemented and deployed on a testbed of 5 AWS DeepLens cameras.


title: Hydrodynamic resistance matrices of colloidal particles with various shapes

abstract: The hydrodynamic resistance matrix is an important quantity for describing the dynamics of colloidal particles. This matrix encodes the shape- and size-dependent hydrodynamic properties of a particle suspended in a simple liquid at low Reynolds number and determines the particle's diffusion tensor. For this reason, the hydrodynamic resistance matrix is typically needed when modeling the motion of free purely Brownian, externally driven, or self-propelled colloidal particles or the behavior of dilute suspensions of such particles on the basis of Langevin equations, Smoluchowski equations, classical dynamical density functional theory, or other appropriate methods. So far, however, the hydrodynamic resistance matrix was available only for a few particle shapes. In this article, we therefore present the hydrodynamic resistance matrices for various particle shapes that are relevant for current research, including apolar and polar as well as convex and partially concave shapes. The elements of the hydrodynamic resistance matrices are given as functions of shape parameters like the aspect ratio of the corresponding particle so that the results apply not only to discrete but instead to continuous sets of particle shapes. This work shall stimulate and support future studies on colloidal particles with anisometric shapes.


title: Multi-Round Cooperative Search Games with Multiple Players

abstract: Assume that a treasure is placed in one of $M$ boxes according to a known distribution and that $k$ searchers are searching for it in parallel during $T$ rounds. We study the question of how to incentivize selfish players so that the success probability, namely, the probability that at least one player finds the treasure, would be maximized. We focus on congestion policies $C(s)$ that specify the reward that a player receives if it is one of $s$ players that (simultaneously) find the treasure for the first time. Our main technical contribution is proving that the exclusive policy, in which $C(1)=1$ and $C(s)=0$ for $s>1$, yields a price of anarchy of $(1-(1-{1}/{k})^{k})^{-1}$, and that this is the best possible price among all symmetric reward mechanisms. For this policy we also have an explicit description of a symmetric equilibrium, which is in some sense unique, and moreover enjoys the best success probability among all symmetric profiles. For general congestion policies, we show how to polynomially find, for any $e>0$, a symmetric multiplicative $(1+e)(1+C(k))$-equilibrium. Together with an appropriate reward policy, a central entity can suggest players to play a particular profile at equilibrium. As our main conceptual contribution, we advocate the use of symmetric equilibria for such purposes. Besides being fair, we argue that in many cases, despite the fact that some small fraction of players crash, symmetric equilibria remain efficient in terms of their group performances and, at the same time, serve as approximate equilibria. We show that this principle holds for a class of games, which we call monotonously scalable games. This applies in particular to our search game, assuming the sharing policy, in which $C(s)=1/s$. For the exclusive policy, this general result does not hold, but we show that the symmetric equilibrium is nevertheless robust under mild assumptions.


title: On a certain subclass of strongly starlike functions

abstract: Let $\mathcal{S}^*_t(\alpha_1,\alpha_2)$ denote the class of functions $f$ analytic in the open unit disc $\Delta$, normalized by the condition $f(0)=0=f'(0)-1$ and satisfying the following two--sided inequality: \begin{equation*} -\frac{\pi\alpha_1}{2}< \arg\left\{\frac{zf'(z)}{f(z)}\right\} <\frac{\pi\alpha_2}{2} \quad (z\in\Delta), \end{equation*} where $0<\alpha_1,\alpha_2\leq1$. The class $\mathcal{S}^*_t(\alpha_1,\alpha_2)$ is a subclass of strongly starlike functions of order $\beta$ where $\beta=\max\{\alpha_1,\alpha_2\}$. The object of the present paper is to derive some certain inequalities including (for example), upper and lower bounds for ${\rm Re}\{zf'(z)/f(z)\}$, growth theorem, logarithmic coefficient estimates and coefficient estimates for functions $f$ belonging to the class $\mathcal{S}^*_t(\alpha_1,\alpha_2)$.


title: Anomaly Detection in Paleoclimate Records using Permutation Entropy

abstract: Permutation entropy techniques can be useful in identifying anomalies in paleoclimate data records, including noise, outliers, and post-processing issues. We demonstrate this using weighted and unweighted permutation entropy of water-isotope records in a deep polar ice core. In one region of these isotope records, our previous calculations revealed an abrupt change in the complexity of the traces: specifically, in the amount of new information that appeared at every time step. We conjectured that this effect was due to noise introduced by an older laboratory instrument. In this paper, we validate that conjecture by re-analyzing a section of the ice core using a more-advanced version of the laboratory instrument. The anomalous noise levels are absent from the permutation entropy traces of the new data. In other sections of the core, we show that permutation entropy techniques can be used to identify anomalies in the raw data that are not associated with climatic or glaciological processes, but rather effects occurring during field work, laboratory analysis, or data post-processing. These examples make it clear that permutation entropy is a useful forensic tool for identifying sections of data that require targeted re-analysis---and can even be useful in guiding that analysis.

categories: physics.data-an cs.IT math.IT physics.ao-ph


title: Building a Winning Self-Driving Car in Six Months

abstract: The SAE AutoDrive Challenge is a three-year competition to develop a Level 4 autonomous vehicle by 2020. The first set of challenges were held in April of 2018 in Yuma, Arizona. Our team (aUToronto/Zeus) placed first. In this paper, we describe our complete system architecture and specialized algorithms that enabled us to win. We show that it is possible to develop a vehicle with basic autonomy features in just six months relying on simple, robust algorithms. We do not make use of a prior map. Instead, we have developed a multi-sensor visual localization solution. All of our algorithms run in real-time using CPUs only. We also highlight the closed-loop performance of our system in detail in several experiments.


title: NET map slope functions

abstract: This paper studies NET map slope functions. It establishes Lipschitz-type conditions for them. It relates Lipschitz-type conditions to the half-space theorem. It gives bounds on the number of slope function fixed points. It provides examples of rational NET maps with many matings.


title: Challenges in detecting evolutionary forces in language change using diachronic corpora

abstract: Newberry et al. (Detecting evolutionary forces in language change, Nature 551, 2017) tackle an important but difficult problem in linguistics, the testing of selective theories of language change against a null model of drift. Having applied a test from population genetics (the Frequency Increment Test) to a number of relevant examples, they suggest stochasticity has a previously under-appreciated role in language evolution. We replicate their results and find that while the overall observation holds, results produced by this approach on individual time series can be sensitive to how the corpus is organized into temporal segments (binning). Furthermore, we use a large set of simulations in conjunction with binning to systematically explore the range of applicability of the Frequency Increment Test. We conclude that care should be exercised with interpreting results of tests like the Frequency Increment Test on individual series, given the researcher degrees of freedom available when applying the test to corpus data, and fundamental differences between genetic and linguistic data. Our findings have implications for selection testing and temporal binning in general, as well as demonstrating the usefulness of simulations for evaluating methods newly introduced to the field.


title: Imaging Carrier Inhomogeneities in Ambipolar Tellurene Field Effect Transistors

abstract: Developing van der Waals (vdW) homojunction devices requires materials with narrow bandgaps and simultaneously high hole and electron mobilities for bipolar transport, as well as methods to image and study spatial variations in carrier type and associated conductivity with nanometer spatial resolution. Here we demonstrate the general capability of near-field scanning microwave microscopy (SMM) to image and study the local carrier type and associated conductivity in operando by studying ambiploar field effect transistors (FETs) of the 1D vdW material tellurium in 2D form. To quantitatively understand electronic variations across the device, we produce nanometer resolved maps of the local carrier equivalence backgate voltage. We show that the global device conductivity minimum determined from transport measurements does not arise from uniform carrier neutrality, but rather from the continued coexistence of p-type regions at the device edge and n-type regions in the interior of our micron-scale devices. This work both underscores and addresses the need to image and understand spatial variations in the electronic properties of nanoscale devices.


title: Optimizations of the Eigensolvers in the ELPA Library

abstract: The solution of (generalized) eigenvalue problems for symmetric or Hermitian matrices is a common subtask of many numerical calculations in electronic structure theory or materials science. Solving the eigenvalue problem can easily amount to a sizeable fraction of the whole numerical calculation. For researchers in the field of computational materials science, an efficient and scalable solution of the eigenvalue problem is thus of major importance. The ELPA-library is a well-established dense direct eigenvalue solver library, which has proven to be very efficient and scalable up to very large core counts. In this paper, we describe the latest optimizations of the ELPA-library for new HPC architectures of the Intel Skylake processor family with an AVX-512 SIMD instruction set, or for HPC systems accelerated with recent GPUs. We also describe a complete redesign of the API in a modern modular way, which, apart from a much simpler and more flexible usability, leads to a new path to access system-specific performance optimizations. In order to ensure optimal performance for a particular scientific setting or a specific HPC system, the new API allows the user to influence in straightforward way the internal details of the algorithms and of performance-critical parameters used in the ELPA-library. On top of that, we introduced an autotuning functionality, which allows for finding the best settings in a self-contained automated way. In situations where many eigenvalue problems with similar settings have to be solved consecutively, the autotuning process of the ELPA-library can be done "on-the-fly". Practical applications from materials science which rely on so-called self-consistency iterations can profit from the autotuning. On some examples of scientific interest, simulated with the FHI-aims application, the advantages of the latest optimizations of the ELPA-library are demonstrated.


title: Interaction of relativistically intense laser pulses with long-scale near critical plasmas for optimization of laser based sources of MeV elec-trons and gamma-rays

abstract: Experiments were performed to study electron acceleration by intense sub-picosecond laser pulses propagating in sub-mm long plasmas of near critical electron density (NCD). Low density foam layers of 300-500 um thickness were used as targets. The NCD-plasma was produced by a mechanism of a super-sonic ionization when a well-defined separate ns-pulse was sent onto the foam-target forerunning the relativistic main pulse. The effect of the relativistic laser pulse channeling and creation of quasi-static azimuthal magnetic and radial electric fields that keeps electrons in the channel ensured effective coupling of the laser energy into energetic electrons. Application of sub-mm thick low density foam layers provided substantial increase of the electron acceleration path in a NCD-plasma compared to the case of freely expanding plasmas created in the interaction of the ns-laser pulse with solid foils. Performed experiments on the electron heating by a 100J, 750 fs short laser pulse of (2-5)x1019 W/cm2 intensity demonstrated that the effective temperature of supra-thermal electrons increased from 1.5-2 MeV, in the case of the relativistic laser interaction with a metallic foil at high laser contrast, up to 13 MeV for the laser shots onto the pre-ionized foam. The observed tendency towards the strong increase of the mean electron energy and the number of ultra-relativistic laser-accelerated electrons is reinforced by the results of gamma-yield measurements that showed a 1000-fold increase of the measured doses. The experiment was supported by the 3D-PIC and FLUKA simulations made for used laser parameters and geometry of the experimental set-up. Both measurements and simulations show high directionality of the acceleration process close to the direction of the laser pulse propagation. The charge of super-ponderomotive electrons with E > 30 MeV reaches a high value of 78nC.


title: A Riemann-Hurwitz-Plucker formula

abstract: We prove a simultaneous generalization of the classical Riemann-Hurwitz and Plucker formulas, addressing the total inflection of a morphism from a (smooth, projective) curve to an arbitrary (smooth, projective) higher-dimensional variety. Our definition of ramification is relative to an algebraic family of divisors on the target variety, and our formula is obtained using the theory of refined top Chern classes. In assigning multiplicities to ramification points, we frequently have to consider excess degeneracy loci, but we are able to show nonetheless that the multiplicities are always nonnegative, and are positive under very mild hypotheses.


title: Nonparametric Spectral Methods for Multivariate Spatial and Spatial-Temporal Data

abstract: We propose computationally efficient methods for estimating stationary multivariate spatial and spatial-temporal spectra from incomplete gridded data. The methods are iterative and rely on successive imputation of data and updating of model estimates. Imputations are done according to a periodic model on an expanded domain. The periodicity of the imputations is a key feature that reduces edge effects in the periodogram and is facilitated by efficient circulant embedding techniques. In addition, we describe efficient methods for decomposing the estimated cross spectral density function into a linear model of coregionalization plus a residual process. The methods are applied to two storm datasets, one of which is from Hurricane Florence, which struck the souteastern United States in September 2018. The application demonstrates how fitted models from different datasets can be compared, and how the methods are computationally feasible on datasets with more than 200,000 total observations.


title: The Severi problem for abelian surfaces in the primitive case

abstract: We prove that the irreducible components of primitive class Severi varieties of general abelian surfaces are completely determined by the maximal factorization through an isogeny of the maps from the normalized curves.


title: Partitions of Matrix Spaces With an Application to $q$-Rook Polynomials

abstract: We study the row-space partition and the pivot partition on the matrix space $\mathbb{F}_q^{n \times m}$. We show that both these partitions are reflexive and that the row-space partition is self-dual. Moreover, using various combinatorial methods, we explicitly compute the Krawtchouk coefficients associated with these partitions. This establishes MacWilliams-type identities for the row-space and pivot enumerators of linear rank-metric codes. We then generalize the Singleton-like bound for rank-metric codes, and introduce two new concepts of code extremality. Both of them generalize the notion of MRD codes and are preserved by trace-duality. Moreover, codes that are extremal according to either notion satisfy strong rigidity properties analogous to those of MRD codes. As an application of our results to combinatorics, we give closed formulas for the $q$-rook polynomials associated with Ferrers diagram boards. Moreover, we exploit connections between matrices over finite fields and rook placements to prove that the number of matrices of rank $r$ over $\mathbb{F}_q$ supported on a Ferrers diagram is a polynomial in $q$, whose degree is strictly increasing in $r$. Finally, we investigate the natural analogues of the MacWilliams Extension Theorem for the rank, the row-space, and the pivot partitions.


title: The Arecibo Pisces-Perseus Supercluster Survey I: Harvesting ALFALFA

abstract: We report a multi-objective campaign of targeted 21 cm HI line observations of sources selected from the Arecibo Legacy Fast ALFA (Arecibo L-band Feed Array) survey (ALFALFA) and galaxies identified by their morphological and photometric properties in the Sloan Digital Sky Survey (SDSS). The aims of this program have been (1) to confirm the reality of some ALFALFA sources whose enigmatic nature suggest additional multiwavelength observations; (2) to probe the low signal-to-noise ratio regime, below the ALFALFA reliability limit; and (3) to explore the feasibility of using optical morphology, color and surface brightness to identify gas-rich objects in the region of the Pisces-Perseus Supercluster (PPS) whose HI fluxes are below the ALFALFA sensitivity limit at that distance. As expected, the reliability of ALFALFA detections depends strongly on the signal-to-noise ratio of the HI line signal and its coincidence with a probable stellar counterpart identified by its optical properties, suggestive of on-going star formation. The identification of low mass, star forming populations enables targeted HI line observations to detect galaxies with HI line fluxes below the ALFALFA sensitivity limits in fixed local volumes (D < 100 Mpc). The method explored here serves as the basis for extending the sample of gas-bearing objects as part of the on-going Arecibo Pisces-Perseus Supercluster Survey (APPSS).


title: New interaction potentials for alkali and alkaline-earth aluminosilicate glasses

abstract: We apply a recently developed optimization scheme to obtain effective potentials for alkali and alkaline-earth aluminosilicate glasses that contains lithium, sodium, potassium, or calcium as modifiers. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by means of ab initio molecular dynamics simulations and density and elastic modulus of glass at room temperature from experiments. The new interaction potentials are able to reproduce reliably the structure and various mechanical and vibrational properties over a wide range of compositions for binary silicates. We have tested these potentials for various ternary systems and find that they are transferable and can be mixed, thus allowing to reproduce and predict the structure and properties of multi-component glasses.