title: Optimizing Precision for Open-World Website Fingerprinting

abstract: Traffic analysis attacks to identify which web page a client is browsing, using only her packet metadata --- known as website fingerprinting --- has been proven effective in closed-world experiments against privacy technologies like Tor. However, due to the base rate fallacy, these attacks have failed in large open-world settings against clients that visit sensitive pages with a low base rate. We find that this is because they have poor precision as they were designed to maximize recall. In this work, we argue that precision is more important than recall for open-world website fingerprinting. For this reason, we develop three classes of {\em precision optimizers}, based on confidence, distance, and ensemble learning, that can be applied to any classifier to increase precision. We test them on known website fingerprinting attacks and show significant improvements in precision. Against a difficult scenario, where the attacker wants to monitor and distinguish 100 sensitive pages each with a low mean base rate of 0.00001, our best optimized classifier can achieve a precision of 0.78; the highest precision of any known attack before optimization was 0.014. We use precise classifiers to tackle realistic objectives in website fingerprinting, including selection, identification, and defeating website fingerprinting defenses.


title: Collision of eigenvalues for matrix-valued processes

abstract: We examine the probability that at least two eigenvalues of an Hermitian matrix-valued Gaussian process, collide. In particular, we determine sharp conditions under which such probability is zero. As an application, we show that the eigenvalues of a real symmetric matrix-valued fractional Brownian motion of Hurst parameter $H$, collide when $H<1/2$ and don't collide when $H>\frac{1}{2}$, while those of a complex Hermitian fractional Brownian motion collide when $H<\frac{1}{3}$ and don't collide when $H>\frac{1}{3}$. Our approach is based on the relation between hitting probabilities for Gaussian processes with the capacity and Hausdorff dimension of measurable sets.


title: Selecting the Best in GANs Family: a Post Selection Inference Framework

abstract: "Which Generative Adversarial Networks (GANs) generates the most plausible images?" has been a frequently asked question among researchers. To address this problem, we first propose an \emph{incomplete} U-statistics estimate of maximum mean discrepancy $\mathrm{MMD}_{inc}$ to measure the distribution discrepancy between generated and real images. $\mathrm{MMD}_{inc}$ enjoys the advantages of asymptotic normality, computation efficiency, and model agnosticity. We then propose a GANs analysis framework to select and test the "best" member in GANs family using the Post Selection Inference (PSI) with $\mathrm{MMD}_{inc}$. In the experiments, we adopt the proposed framework on 7 GANs variants and compare their $\mathrm{MMD}_{inc}$ scores.


title: NtMalDetect: A Machine Learning Approach to Malware Detection Using Native API System Calls

abstract: As computing systems become increasingly advanced and as users increasingly engage themselves in technology, security has never been a greater concern. In malware detection, static analysis, the method of analyzing potentially malicious files, has been the prominent approach. This approach, however, quickly falls short as malicious programs become more advanced and adopt the capabilities of obfuscating its binaries to execute the same malicious functions, making static analysis extremely difficult for newer variants. The approach assessed in this paper is a novel dynamic malware analysis method, which may generalize better than static analysis to newer variants. Inspired by recent successes in Natural Language Processing (NLP), widely used document classification techniques were assessed in detecting malware by doing such analysis on system calls, which contain useful information about the operation of a program as requests that the program makes of the kernel. Features considered are extracted from system call traces of benign and malicious programs, and the task to classify these traces is treated as a binary document classification task of system call traces. The system call traces were processed to remove the parameters to only leave the system call function names. The features were grouped into various n-grams and weighted with Term Frequency-Inverse Document Frequency. This paper shows that Linear Support Vector Machines (SVM) optimized by Stochastic Gradient Descent and the traditional Coordinate Descent on the Wolfe Dual form of the SVM are effective in this approach, achieving a highest of 96% accuracy with 95% recall score. Additional contributions include the identification of significant system call sequences that could be avenues for further research.


title: Inverse Gauss curvature flows with free boundaries in a cone

abstract: We consider strictly convex hypersurfaces with the boundary which meets a strictly convex cone perpendicularly. We prove that if these hypersurfaces expand inside this cone, driven by the power of the Gauss curvature, then the evolution exists for all the time and the evolving hypersurfaces converge smoothly to a piece of round sphere after rescaling.


title: Efficiency of Centrifugal Mechanism in Producing PeV Neutrinos From Active Galactic Nuclei

abstract: A several-step theoretical model is constructed to trace the origin of ultra high energy (UHE) $[1-2]$PeV neutrinos detected, recently, by the IceCube collaboration. Protons in the AGN magnetosphere, experiencing different gravitational centrifugal force, provide free energy for the parametric excitation of Langmuir waves via a generalized two-stream instability. Landau damping of these waves, outside the AGN magnetosphere, can accelerate protons to ultra high energies. The ultimate source for this mechanism, the Langmuir-Landau-Centrifugal-Drive (LLCD), is the gravitational energy of the compact object. The LLCD generated UHE protons provide the essential ingredient in the creation of UHE neutrinos via appropriate hadronic reactions; protons of energy~ $10^{17}$eV can be generated in the plasmas surrounding AGN with bolometric luminosities of the order of $10^{43}$ergs s$^{-1}$. By estimating the diffusive energy flux of extragalactic neutrinos in the energy interval $[1-2]$PeV, we find that an acceptably small fraction $0.003\%$ of the total bolometric luminosity will suffice to create the observed fluxes of extragalactic ultra-high energy neutrinos.


title: Teaching Machines to Code: Neural Markup Generation with Visual Attention

abstract: We present a neural transducer model with visual attention that learns to generate LaTeX markup of a real-world math formula given its image. Applying sequence modeling and transduction techniques that have been very successful across modalities such as natural language, image, handwriting, speech and audio; we construct an image-to-markup model that learns to produce syntactically and semantically correct LaTeX markup code over 150 words long and achieves a BLEU score of 89%; improving upon the previous state-of-art for the Im2Latex problem. We also demonstrate with heat-map visualization how attention helps in interpreting the model and can pinpoint (detect and localize) symbols on the image accurately despite having been trained without any bounding box data.


title: Magnetic moments of the lowest-lying singly heavy baryons

abstract: A light baryon is viewed as $N_c$ valence quarks bound by meson mean fields in the large $N_c$ limit. In much the same way a singly heavy baryon is regarded as $N_c-1$ valence quarks bound by the same mean fields, which makes it possible to use the properties of light baryons to investigate those of the heavy baryons. A heavy quark being regarded as a static color source in the limit of the infinitely heavy quark mass, the magnetic moments of the heavy baryon are determined entirely by the chiral soliton consisting of a light-quark pair. The magnetic moments of the baryon sextet are obtained by using the parameters fixed in the light-baryon sector. In this mean-field approach, the numerical results of the magnetic moments of the baryon sextet with spin $3/2$ are just 3/2 larger than those with spin $1/2$. The magnetic moments of the bottom baryons are the same as those of the corresponding charmed baryons.


title: Extraction of $|V_{cb}|$ from two-body hadronic B decays

abstract: We propose a method of extracting the Cabibbo-Kobayashi-Maskawa matrix element $|V_{cb}|$ from two-body hadronic decay processes of $B\to DK$ with precisely determined form factors of $B$ meson semi-leptonic decays. The amplitude $\mathcal{M}(\bar{B}^0 \to D^+ K^-)$ which does not include the effect of hadronic final state interactions can be theoretically evaluated by using factorization and form factors of semi-leptonic B decays. We can obtain all the amplitudes in an isospin relation $\mathcal{A}(B^-\to D^0K^-) =\mathcal{A}(\bar{B}^0\to D^+K^{-})+\mathcal{A}(\bar{B}^0\to D^{0}\bar{K}^0)$ including the effect of hadronic final state interactions as well as $|V_{cb}|$ using the experimental data of branching fractions of these three processes with a truncation of the states which contribute to the hadronic final state interactions. The extracted value of $|V_{cb}|$ is $(37\pm 6)\times 10^{-3}$. The decay processes of $B\to DK^{*}$ and $B\to D^{*}K$ can also be used in the same way and the extracted values of $|V_{cb}|$ are $(41\pm 7)\times 10^{-3}$ and $(42\pm 9)\times 10^{-3}$, respectively. This method becomes possible by virtue of recent precise determinations of the form factors of semi-leptonic B decays. The uncertainties of $|V_{cb}|$ by this method are expected to be reduced by the results of future B-factory experiments and lattice calculations.


title: Dual Garside structures and Coxeter sortable elements

abstract: In Artin-Tits groups attached to Coxeter groups of spherical type, we give a combinatorial formula to express the simple elements of the dual braid monoids in the classical Artin generators. Every simple dual braid is obtained by lifting an $S$-reduced expression of its image in the Coxeter group, in a way which involves Reading's $c$-sortable elements. It has as an immediate consequence that simple dual braids are Mikado braids (the known proofs of this result either require topological realizations of the Artin groups or categorification techniques), and hence that their images in the Iwahori-Hecke algebras have positivity properties. In the classical types, this requires to give an explicit description of the inverse of Reading's bijection from $c$-sortable elements to noncrossing partitions of a Coxeter element $c$, which might be of independent interest. The bijections are described in terms of the noncrossing partition models in these types. While the proof of the formula is case-by-case, it is entirely combinatorial and we develop an approach which reduces a uniform proof to uniformly proving a lemma about inversion sets of $c$-sortable elements.


title: Seismic response in modern cities

abstract: The proposed homogeneous flat-faced layer-like model of a city (termed overlayer), covering what is generally considered to be a dangerous site (from the point of view of seismic hazard) lends itself to an explicit theoretical analysis of its response to a seismic body wave radiated by distant sources. This study is carried out for: ground response of the complete site/overlayer configuration which is compared to the response of the configuration in which the overlayer is absent, response at the top of the layer for various layer thicknesses, and determination, as a function of frequency, of the fraction of incident flux that is dissipated in the overlayer, the underlying layer and, by radiation damping, in the hard half space. It is shown that all of these entities are highly frequency-dependent and even large in certain frequency intervals, without any resonant (in the sense of mode excitation) phenomena coming into play. The results of this study also show that transfer functions do not necessarily reflect the global response in the built component of a city and that more-appropriate energy-related functions, termed spectral absorptance (in the blocks of the city or their layer-like surrogate at the characteristic frequency of the seismic pulse) and absorptance (integral over frequency of the spectral absorptance), can increase with increasing city density or increasing city height. In fact, it is shown that more than a third of the incident seismic energy can be sent into, and therefore cause serious damage to, the built component. On the basis of these findings, it appears that the probable evolution of the morphology and constitutive properties of cities with time will make the latter more vulnerable to damage and destruction when submitted to seismic waves.


title: On the Power-of-d-choices with Least Loaded Server Selection

abstract: Motivated by distributed schedulers that combine the power-of-d-choices with late binding and systems that use replication with cancellation-on-start, we study the performance of the LL(d) policy which assigns a job to a server that currently has the least workload among d randomly selected servers in large-scale homogeneous clusters. We consider general service time distributions and propose a partial integro-differential equation to describe the evolution of the system. This equation relies on the earlier proven ansatz for LL(d) which asserts that the workload distribution of any finite set of queues becomes independent of one another as the number of servers tends to infinity. Based on this equation we propose a fixed point iteration for the limiting workload distribution and study its convergence. For exponential job sizes we present a simple closed form expression for the limiting workload distribution that is valid for any work-conserving service discipline as well as for the limiting response time distribution in case of first-come-first-served scheduling. We further show that for phase-type distributed job sizes the limiting workload and response time distribution can be expressed via the unique solution of a simple set of ordinary differential equations. Numerical and analytical results that compare response time of the classic power-of-d-choices algorithm and the LL(d) policy are also presented and the accuracy of the limiting response time distribution for finite systems is illustrated using simulation.


title: CADDeLaG: Framework for distributed anomaly detection in large dense graph sequences

abstract: Random walk based distance measures for graphs such as commute-time distance are useful in a variety of graph algorithms, such as clustering, anomaly detection, and creating low dimensional embeddings. Since such measures hinge on the spectral decomposition of the graph, the computation becomes a bottleneck for large graphs and do not scale easily to graphs that cannot be loaded in memory. Most existing graph mining libraries for large graphs either resort to sampling or exploit the sparsity structure of such graphs for spectral analysis. However, such methods do not work for dense graphs constructed for studying pairwise relationships among entities in a data set. Examples of such studies include analyzing pairwise locations in gridded climate data for discovering long distance climate phenomena. These graphs representations are fully connected by construction and cannot be sparsified without loss of meaningful information. In this paper we describe CADDeLaG, a framework for scalable computation of commute-time distance based anomaly detection in large dense graphs without the need to load the entire graph in memory. The framework relies on Apache Spark's memory-centric cluster-computing infrastructure and consists of two building blocks: a decomposable algorithm for commute time distance computation and a distributed linear system solver. We illustrate the scalability of CADDeLaG and its dependency on various factors using both synthetic and real world data sets. We demonstrate the usefulness of CADDeLaG in identifying anomalies in a climate graph sequence, that have been historically missed due to ad hoc graph sparsification and on an election donation data set.


title: Photoionization of the bound systems at high energies

abstract: We consider photoionization of a system bound by the central potential $V(r)$. We demonstrate that the high energy nonrelativistic asymptotics of the photoionization cross section can be obtained without solving the wave equation. The asymptotics can be expressed in terms of the Fourier transform of the potential by employing the Lippmann--Schwinger equation. We find the asymptotics for the screened Coulomb field. We demonstrate that the leading corrections to this asymptotics are described by the universal factor. The high energy nonrelativistic asymptotics is found to be determined by the analytic properties of the potential $V(r)$. We show that the energy dependence of the asymptotics of photoionization cross sections of fullerenes is to large extent model dependent. We demonstrate that if the fullerene field $V(r)$ is approximated by the function with singularities in the complex plane, the power drop of the asymptotics is reached at the energies which as so high that the cross section becomes unobservably small. The preasymptotic behavior with a faster drop of the cross sections becomes important in these cases.


title: The Meissner Effect for axially symmetric charged black holes

abstract: In our previous work [N. G\"urlebeck, M. Scholtz, Phys. Rev. D 95 064010 (2017)], we have shown that electric and magnetic fields are expelled from the horizons of extremal, stationary and axially symmetric uncharged black holes; this is called the Meissner effect for black holes. Here, we generalize this result in several directions. First, we allow that the black hole carries charge, which requires a generalization of the definition of the Meissner effect. Next, we introduce the notion of almost isolated horizons, which is weaker than the usual notion of isolated horizons, since the geometry of the former is not necessarily completely time-independent. Moreover, we allow the horizon to be pierced by strings, thereby violating the usual assumption on the spherical topology made in the definition of the weakly isolated horizon. Finally, we spell out in detail all assumptions entering the proof and show that the Meissner effect is an inherent property of black holes even in full non-linear theory.


title: Threshold response and bistability in gene regulation by small noncoding RNA

abstract: In this paper, we study through mathematical modelling the combined effect of transcriptional and translational regulation by proteins and small noncoding RNAs (sRNA) in a genetic feedback motif that has an important role in the survival of E.coli under stress associated with oxygen and energy availability. We show that subtle changes in this motif can bring in drastically different effects on the gene expression. In particular, we show that a threshold response in the gene expression changes to a bistable response as the regulation on sRNA synthesis or degradation is altered. These results are obtained under deterministic conditions. Next, we study how the gene expression is altered by additive and multiplicative noise which might arise due to probabilistic occurrences of different biochemical events. Using the Fokker-Planck formulation, we obtain steady state probability distributions for sRNA concentration for the network motifs displaying bistability. The probability distributions are found to be bimodal with two peaks at low and high concentrations of sRNAs. We further study the variations in the probability distributions under different values of noise strength and correlations. The results presented here might be of interest for designing synthetic network for artificial control.


abstract: In this paper, we prove the Miyaoka-Yau inequality for compact K\"ahler manifolds with semi-positive canonical bundle. The key point of the proof is the estimate for the $L^2$-norm of the scalar curvature along the K\"ahler-Ricci flow.


title: Accelerating Adaptive Cubic Regularization of Newton's Method via Random Sampling

abstract: In this paper, we consider an unconstrained optimization model where the objective is a sum of a large number of possibly nonconvex functions, though overall the objective is assumed to be smooth and convex. Our bid to solving such model uses the framework of cubic regularization of Newton's method. As well known, the crux in cubic regularization is its utilization of the Hessian information, which may be computationally expensive for large-scale problems. To tackle this, we resort to approximating the Hessian matrix via sub-sampling. In particular, we propose to compute an approximated Hessian matrix by either \textit{uniformly}\/ or \textit{non-uniformly}\/ sub-sampling the components of the objective. Based upon such sampling strategy, we develop accelerated adaptive cubic regularization approaches and provide theoretical guarantees on global iteration complexity of $O(\epsilon^{-1/3})$ with high probability, which matches that of the original accelerated cubic regularization methods \cite{Jiang-2017-Unified} using the \textit{full}\/ Hessian information. Interestingly, we show that in the worst case scenario our algorithm still achieves an $O\left(\log(\epsilon^{-1})\epsilon^{-5/6}\right)$ iteration complexity bound. The performances of the proposed methods on the regularized logistic regression problems show a clear effect of acceleration in terms of the epoch counts on several real data sets.


title: Implementation of Massive MIMO Uplink Receiver on RaPro Prototyping Platform

abstract: The updated physical layer standard of the fifth generation wireless communication suggests the necessity of a rapid prototyping platform. To this end, we develop RaPro, a multi-core general purpose processor-based massive multiple-input-multiple-output (MIMO) prototyping platform. To enhance RaPro, high performance detection and beamforming are needed, whereas both of them request for accurate channel state information (CSI). In this paper, linear minimum mean square error (LMMSE)-based channel estimator is adopted and encapsulated inside RaPro to gain more accurate CSI. Considering the high comlexity and unknown of channel statistics, we design low-complexity LMMSE channel estimator to alleviate the rising complexity along with increasing antenna number and set more computational resource aside for massive MIMO uplink detection and downlink beamforming. Simulation results indicate the high mean square error performance and robustness of designed low-complexity method. Indoor and corridor scenario tests show prominent improvement in bit error rate performance. Time cost analysis proves the practical use and real-time transmission ability of the implemented uplink receiver on RaPro.


title: Quantum Circuit Design for Training Perceptron Models

abstract: Perceptron model is a fundamental linear classifier in machine learning and also the building block of artificial neural networks. Recently, Wiebe et al. (arXiv:1602.04799) proposed that the training of a perceptron can be quadratically speeded using Grover search with a quantum computer, which has potentially important big-data applications. In this paper, we design a quantum circuit for implementing this algorithm. The Grover oracle, the central part of the circuit, is realized by Quantum-Fourier-Transform based arithmetics that specifies whether an input weight vector can correctly classify all training data samples. We also analyze the required number of qubits and universal gates for the algorithm, as well as the success probability using uniform sampling, showing that it has higher possibility than spherical Gaussian distribution $N(0,1)$. The feasibility of the circuit is demonstrated by a testing example using the IBM-Q cloud quantum computer, where 16 qubits are used to classify four data samples.


title: Blind Source Separation with Optimal Transport Non-negative Matrix Factorization

abstract: Optimal transport as a loss for machine learning optimization problems has recently gained a lot of attention. Building upon recent advances in computational optimal transport, we develop an optimal transport non-negative matrix factorization (NMF) algorithm for supervised speech blind source separation (BSS). Optimal transport allows us to design and leverage a cost between short-time Fourier transform (STFT) spectrogram frequencies, which takes into account how humans perceive sound. We give empirical evidence that using our proposed optimal transport NMF leads to perceptually better results than Euclidean NMF, for both isolated voice reconstruction and BSS tasks. Finally, we demonstrate how to use optimal transport for cross domain sound processing tasks, where frequencies represented in the input spectrograms may be different from one spectrogram to another.


title: Virtual walks in spin space: a study in a family of two-parameter models

abstract: We investigate the dynamics of classical spins mapped as walkers in a virtual "spin" space using a generalised two-parameter family of spin models characterized by parameters $y$ and $z$ [M. J. de Oliveira, J. F. F. Mendes and M. A. Santos, J. Phys. A Math. Gen. \textbf{26}, 2317 (1993)]. The behavior of $S(x,t)$, the probability that the walker is at position $x$ at time $t$ is studied in detail. In general $S(x,t) \sim t^{-\alpha}f(x/t^{\alpha})$ with $\alpha \simeq 1$ or $0.5$ at large times depending on the parameters. In particular, $S(x,t)$ for the point $y=1, z=0.5$ corresponding to the voter model shows a crossover in time; associated with this crossover, two timescales can be defined which vary with the system size $L$ as $L^2\log L$. We also show that as the voter model point is approached from the disordered regions along different directions, the width of the Gaussian distribution $S(x,t)$ diverges in a power law manner with different exponents. For the majority voter case, the results indicate that the the virtual walk can detect the phase transition perhaps more efficiently compared to other non-equilibrium methods.


title: On the Theory of Variance Reduction for Stochastic Gradient Monte Carlo

abstract: We provide convergence guarantees in Wasserstein distance for a variety of variance-reduction methods: SAGA Langevin diffusion, SVRG Langevin diffusion and control-variate underdamped Langevin diffusion. We analyze these methods under a uniform set of assumptions on the log-posterior distribution, assuming it to be smooth, strongly convex and Hessian Lipschitz. This is achieved by a new proof technique combining ideas from finite-sum optimization and the analysis of sampling methods. Our sharp theoretical bounds allow us to identify regimes of interest where each method performs better than the others. Our theory is verified with experiments on real-world and synthetic datasets.


title: On the timing behavior of PSR B1259-63 under the propeller torque from a transient accretion disc

abstract: The $\gamma$-ray pulsar binary system PSR B1259-63 flares in GeV after each periastron. The origin of those flares is still under debate. Recently, [2017ApJ...844..114Y] proposed a mechanism that might explain the GeV flares. In that model, a transient accretion disc is expected to be formed from the matter which was gravity-captured by the neutron star from the main sequence companion's circumstellar disc. The transient accretion disc exerts a spin down torque on the neutron star (propeller effect), which might be traceable via pulsar timing observation of PSR B1259-63. In this paper, we phenomenologically consider the propeller effect with a parameter $\chi$, which describes the coupling between the disc matter and the neutron star. Comparing the expected timing residuals against the recent observation in [2014MNRAS.437.3255S], we conclude that the angular momentum transfer is very weak (with the coupling parameter $\chi\le10^{-4}$).


title: Characterising the observational properties of {\delta} Sct stars in the era of space photometry from the Kepler mission

abstract: The {\delta} Sct stars are a diverse group of intermediate-mass pulsating stars located on and near the main sequence within the classical instability strip in the Hertzsprung-Russell diagram. Many of these stars are hybrid stars pulsating simultaneously with pressure and gravity modes that probe the physics at different depths within a star's interior. Using two large ensembles of {\delta} Sct stars observed by the Kepler Space Telescope, the instrumental biases inherent to Kepler mission data and the statistical properties of these stars are investigated. An important focus of this work is an analysis of the relationships between the pulsational and stellar parameters, and their distribution within the classical instability strip. It is found that a non-negligible fraction of main sequence {\delta} Sct stars exist outside theoretical predictions of the classical instability boundaries, which indicates the necessity of a mass-dependent mixing length parameter to simultaneously explain low- and high-radial order pressure modes in {\delta} Sct stars within the Hertzsprung-Russell diagram. Furthermore, a search for regularities in the amplitude spectra of these stars is also presented, specifically the frequency difference between pressure modes of consecutive radial order. In this work, it is demonstrated that an ensemble-based approach using space photometry from the Kepler mission is not only plausible for {\delta} Sct stars, but that it is a valuable method for identifying the most promising stars for mode identification and asteroseismic modelling. The full scientific potential of studying {\delta} Sct stars is as yet unrealised. The ensembles discussed in this paper represent a high-quality data set for future studies of rotation and angular momentum transport inside A and F stars using asteroseismology.


title: Modeling Venus-Like Worlds Through Time

abstract: We explore the atmospheric and surface history of a hypothetical paleo-Venus climate using a 3-D General Circulation Model. We constrain our model with the in-situ and remote sensing Venus data available today. Given that Venus and Earth are believed to be similar geochemically some aspects of Earth's history are also utilized. We demonstrate that it is possible for ancient Venus and Venus-like exoplanetary worlds to exist within the liquid water habitable zone with insolations up to nearly 2 times that of modern Earth.


title: Analysis of the Web Graph Aggregated by Host and Pay-Level Domain

abstract: In this paper the web is analyzed as a graph aggregated by host and pay-level domain (PLD). The web graph datasets, publicly available, have been released by the Common Crawl Foundation and are based on a web crawl performed during the period May-June-July 2017. The host graph has $\sim$1.3 billion nodes and $\sim$5.3 billion arcs. The PLD graph has $\sim$91 million nodes and $\sim$1.1 billion arcs. We study the distributions of degree and sizes of strongly/weakly connected components (SCC/WCC) focusing on power laws detection using statistical methods. The statistical plausibility of the power law model is compared with that of several alternative distributions. While there is no evidence of power law tails on host level, they emerge on PLD aggregation for indegree, SCC and WCC size distributions. Finally, we analyze distance-related features by studying the cumulative distributions of the shortest path lengths, and give an estimation of the diameters of the graphs.


title: Detection of fractional solitons in quantum spin Hall systems

abstract: We propose two experimental setups that allow for the implementation and the detection of fractional solitons of the Goldstone-Wilczek type. The first setup is based on two magnetic barriers at the edge of a quantum spin Hall system for generating the fractional soliton. If then a quantum point contact is created with the other edge, the linear conductance shows evidence of the fractional soliton. The second setup consists of a single magnetic barrier covering both edges and implementing a long quantum point contact. In this case, the fractional soliton can unambiguously be detected as a dip in the conductance without the need to control the magnetization of the barrier.


title: Optimization of the light intensity for Photodetectors calibration

abstract: In this article we present an evaluation of the uncertainty in the average number of photoelectrons, which is important for the calibration of photodetectors. We show that the statistical uncertainty depends on light intensity, and on the method of evaluation. For some cases there is optimal light intensity where the accuracy reaches its optimal value with fixed statistics. A method of photoelectron evaluation based on the extraction of pedestal's (zero) probability gives the best accuracy at approximately 1.6 photoelectrons for a noiseless photodetector and shifts out to higher values with the presence of noise. In the general case, estimation of the average number of photoelectrons is biased and might need special consideration.


title: Mean Field Multi-Agent Reinforcement Learning

abstract: Existing multi-agent reinforcement learning methods are limited typically to a small number of agents. When the agent number increases largely, the learning becomes intractable due to the curse of the dimensionality and the exponential growth of agent interactions. In this paper, we present \emph{Mean Field Reinforcement Learning} where the interactions within the population of agents are approximated by those between a single agent and the average effect from the overall population or neighboring agents; the interplay between the two entities is mutually reinforced: the learning of the individual agent's optimal policy depends on the dynamics of the population, while the dynamics of the population change according to the collective patterns of the individual policies. We develop practical mean field Q-learning and mean field Actor-Critic algorithms and analyze the convergence of the solution to Nash equilibrium. Experiments on Gaussian squeeze, Ising model, and battle games justify the learning effectiveness of our mean field approaches. In addition, we report the first result to solve the Ising model via model-free reinforcement learning methods.


title: Mitigation of dynamical instabilities in laser arrays via non-Hermitian coupling

abstract: Arrays of coupled semiconductor lasers are systems possessing complex dynamical behavior that are of major interest in photonics and laser science. Dynamical instabilities, arising from supermode competition and slow carrier dynamics, are known to prevent stable phase locking in a wide range of parameter space, requiring special methods to realize stable laser operation. Inspired by recent concepts of parity-time ($\mathcal{PT}$) and non-Hermitian photonics, in this work we consider non-Hermitian coupling engineering in laser arrays in a ring geometry and show, both analytically and numerically, that non-Hermitian coupling can help to mitigate the onset of dynamical laser instabilities. In particular, we consider in details two kinds of nearest-neighbor non-Hermitian couplings: symmetric but complex mode coupling (type-I non-Hermitian coupling) and asymmetric mode coupling (type-II non-Hermitian coupling). Suppression of dynamical instabilities can be realized in both coupling schemes, resulting in stable phase-locking laser emission with the lasers emitting in phase (for type-I coupling) or with $\pi/2$ phase gradient (for type-II coupling), resulting in a vortex far-field beam. In type-II non-Hermitian coupling, chirality induced by asymmetric mode coupling enables laser phase locking even in presence of moderate disorder in the resonance frequencies of the lasers.


title: Improving the Decoding Threshold of Tailbiting Spatially Coupled LDPC Codes by Energy Shaping

abstract: We show how the iterative decoding threshold of tailbiting spatially coupled (SC) low-density parity-check (LDPC) code ensembles can be improved over the binary input additive white Gaussian noise channel by allowing the use of different transmission energies for the codeword bits. We refer to the proposed approach as energy shaping. We focus on the special case where the transmission energy of a bit is selected among two values, and where a contiguous portion of the codeword is transmitted with the largest one. Given these constraints, an optimal energy boosting policy is derived by means of protograph extrinsic information transfer analysis. We show that the threshold of tailbiting SC-LDPC code ensembles can be made close to that of terminated code ensembles while avoiding the rate loss (due to termination). The analysis is complemented by Monte Carlo simulations, which confirm the viability of the approach.


abstract: This is the paper "Niels Henrik Abel and the birth of fractional calculus", Podlubny, I., Magin, R. L., Trymorush I., Fractional Calculus and Applied Analysis, vol.20, no.5, pp.1068-1075, 2017 (https://doi.org/10.1515/fca-2017-0057) with the supplements to it. Hereby translation to English of the first part of the important Abel's paper (Abel, N. H., Oplosning af et par opgaver ved hjelp af bestemte integraler, Magazin for Naturvidenskaberne, Aargang I, Bind 2, Christiania, 1823) and another important Abel's paper (Abel, N. H., Auflosung einer mechanischen ausgabe. Journal fur die Reine und Angewandte Mathematik, Bind I, 153-157, Berlin, 1826) are provided. To preserve the original flavor and the historical development of mathematical notation, the formulas are typeset in the original manner.


title: Strain-assisted optomechanical coupling of polariton condensate spin to a micromechanical resonator

abstract: We report spin and intensity coupling of an exciton-polariton condensate to the mechanical vibrations of a circular membrane microcavity. We optically drive the microcavity resonator at the lowest mechanical resonance frequency while creating an optically-trapped spin-polarized polariton condensate in different locations on the microcavity, and observe spin and intensity oscillations of the condensate at the vibration frequency of the resonator. Spin oscillations are induced by vibrational strain driving, whilst the modulation of the optical trap due to the displacement of the membrane causes intensity oscillations in the condensate emission. Our results demonstrate spin-phonon coupling in a macroscopically coherent condensate.


title: Genomics as a Service: a Joint Computing and Networking Perspective

abstract: This paper provides a global picture about the deployment of networked processing services for genomic data sets. Many current research make an extensive use genomic data, which are massive and rapidly increasing over time. They are typically stored in remote databases, accessible by using Internet. For this reason, a significant issue for effectively handling genomic data through data networks consists of the available network services. A first contribution of this paper consists of identifying the still unexploited features of genomic data that could allow optimizing their networked management. The second and main contribution of this survey consists of a methodological classification of computing and networking alternatives which can be used to offer what we call the Genomic-as-a-Service (GaaS) paradigm. In more detail, we analyze the main genomic processing applications, and classify not only the main computing alternatives to run genomics workflows in either a local machine or a distributed cloud environment, but also the main software technologies available to develop genomic processing services. Since an analysis encompassing only the computing aspects would provide only a partial view of the issues for deploying GaaS system, we present also the main networking technologies that are available to efficiently support a GaaS solution. We first focus on existing service platforms, and analyze them in terms of service features, such as scalability, flexibility, and efficiency. Then, we present a taxonomy for both wide area and datacenter network technologies that may fit the GaaS requirements. It emerges that virtualization, both in computing and networking, is the key for a successful large-scale exploitation of genomic data, by pushing ahead the adoption of the GaaS paradigm. Finally, the paper illustrates a short and long-term vision on future research challenges in the field.


title: A Weighted Likelihood Approach Based on Statistical Data Depths

abstract: We propose a general approach to construct weighted likelihood estimating equations with the aim of obtain robust estimates. The weight, attached to each score contribution, is evaluated by comparing the statistical data depth at the model with that of the sample in a given point. Observations are considered regular when the ratio of these two depths is close to one, whereas, when the ratio is large the corresponding score contribution may be downweigthed. Details and examples are provided for the robust estimation of the parameters in the multivariate normal model. Because of the form of the weights, we expect that, there will be no downweighting under the true model leading to highly efficient estimators. Robustness is illustrated using two real data sets.


title: Flat-space Holography and Correlators of Robinson-Trautman Stress tensor

abstract: We propose a quasi-local stress tensor for the four-dimensional asymptotically flat Robinson-Trautman geometries by taking the flat-space limit from the corresponding asymptotically AdS solutions. This stress tensor results in the correct charges of the generators of BMS symmetry if we define conformal infinity by an anisotropic scaling of the metric components. Using flat-space holography this stress tensor is related to expectation values of the stress tensor in a dual field theory called BMS-invariant field theory (BMSFT). We also calculate the two and three point functions of the proposed stress tensor.


title: A 2x2 quantum dot array with controllable inter-dot tunnel couplings

abstract: The interaction between electrons in arrays of electrostatically defined quantum dots is naturally described by a Fermi-Hubbard Hamiltonian. Moreover, the high degree of tunability of these systems make them a powerful platform to simulate different regimes of the Hubbard model. However, most quantum dot array implementations have been limited to one-dimensional linear arrays. In this letter, we present a square lattice unit cell of 2$\times$2 quantum dots defined electrostatically in a AlGaAs/GaAs heterostructure using a double-layer gate technique. We probe the properties of the array using nearby quantum dots operated as charge sensors. We show that we can deterministically and dynamically control the charge occupation in each quantum dot in the single- to few-electron regime. Additionally, we achieve simultaneous individual control of the nearest-neighbor tunnel couplings over a range 0-40~$\mu$eV. Finally, we demonstrate fast ($\sim 1$~$\mu$s) single-shot readout of the spin state of electrons in the dots, through spin-to-charge conversion via Pauli spin blockade. These advances pave the way to analog quantum simulations in two dimensions, not previously accessible in quantum dot systems.


title: History PCA: A New Algorithm for Streaming PCA

abstract: In this paper we propose a new algorithm for streaming principal component analysis. With limited memory, small devices cannot store all the samples in the high-dimensional regime. Streaming principal component analysis aims to find the $k$-dimensional subspace which can explain the most variation of the $d$-dimensional data points that come into memory sequentially. In order to deal with large $d$ and large $N$ (number of samples), most streaming PCA algorithms update the current model using only the incoming sample and then dump the information right away to save memory. However the information contained in previously streamed data could be useful. Motivated by this idea, we develop a new streaming PCA algorithm called History PCA that achieves this goal. By using $O(Bd)$ memory with $B\approx 10$ being the block size, our algorithm converges much faster than existing streaming PCA algorithms. By changing the number of inner iterations, the memory usage can be further reduced to $O(d)$ while maintaining a comparable convergence speed. We provide theoretical guarantees for the convergence of our algorithm along with the rate of convergence. We also demonstrate on synthetic and real world data sets that our algorithm compares favorably with other state-of-the-art streaming PCA methods in terms of the convergence speed and performance.


title: Discrepancy-based Evolutionary Diversity Optimization

abstract: Diversity plays a crucial role in evolutionary computation. While diversity has been mainly used to prevent the population of an evolutionary algorithm from premature convergence, the use of evolutionary algorithms to obtain a diverse set of solutions has gained increasing attention in recent years. Diversity optimization in terms of features on the underlying problem allows to obtain a better understanding of possible solutions to the problem at hand and can be used for algorithm selection when dealing with combinatorial optimization problems such as the Traveling Salesperson Problem. We explore the use of the star-discrepancy measure to guide the diversity optimization process of an evolutionary algorithm. In our experimental investigations, we consider our discrepancy-based diversity optimization approaches for evolving diverse sets of images as well as instances of the Traveling Salesperson problem where a local search is not able to find near optimal solutions. Our experimental investigations comparing three diversity optimization approaches show that a discrepancy-based diversity optimization approach using a tie-breaking rule based on weighted differences to surrounding feature points provides the best results in terms of the star discrepancy measure.


title: A range condition for polyconvex variational regularization

abstract: In the context of convex variational regularization it is a known result that, under suitable differentiability assumptions, source conditions in the form of variational inequalities imply range conditions, while the converse implication only holds under an additional restriction on the operator. In this article we prove the analogous result for polyconvex regularization. More precisely, we show that the variational inequality derived by Kirisits, Scherzer (2017) implies that the derivative of the regularization functional must lie in the range of the dual-adjoint of the derivative of the operator. In addition, we show how to adapt the restriction on the operator in order to obtain the converse implication.


title: Hierarchical multi-stage MCMC follow-up of continuous gravitational wave candidates

abstract: Leveraging Markov chain Monte Carlo (MCMC) optimization of the F-statistic, we introduce a method for the hierarchical follow-up of continuous gravitational wave candidates identified by wide-parameter space semi-coherent searches. We demonstrate parameter estimation for continuous wave sources and develop a framework and tools to understand and control the effective size of the parameter space, critical to the success of the method. Monte Carlo tests of simulated signals in noise demonstrate that this method is close to the theoretical optimal performance.


title: Mapping Images to Scene Graphs with Permutation-Invariant Structured Prediction

abstract: Machine understanding of complex images is a key goal of artificial intelligence. One challenge underlying this task is that visual scenes contain multiple inter-related objects, and that global context plays an important role in interpreting the scene. A natural modeling framework for capturing such effects is structured prediction, which optimizes over complex labels, while modeling within-label interactions. However, it is unclear what principles should guide the design of a structured prediction model that utilizes the power of deep learning components. Here we propose a design principle for such architectures that follows from a natural requirement of permutation invariance. We prove a necessary and sufficient characterization for architectures that follow this invariance, and discuss its implication on model design. Finally, we show that the resulting model achieves new state of the art results on the Visual Genome scene graph labeling benchmark, outperforming all recent approaches.


title: A dynamical approach in exploring the unknown mass in the Solar system using pulsar timing arrays

abstract: The error in the Solar system ephemeris will lead to dipolar correlations in the residuals of pulsar timing array for widely separated pulsars. In this paper, we utilize such correlated signals, and construct a Bayesian data-analysis framework to detect the unknown mass in the Solar system and to measure the orbital parameters. The algorithm is designed to calculate the waveform of the induced pulsar-timing residuals due to the unmodelled objects following the Keplerian orbits in the Solar system. The algorithm incorporates a Bayesian-analysis suit used to simultaneously analyse the pulsar-timing data of multiple pulsars to search for coherent waveforms, evaluate the detection significance of unknown objects, and to measure their parameters. When the object is not detectable, our algorithm can be used to place upper limits on the mass. The algorithm is verified using simulated data sets, and cross-checked with analytical calculations. We also investigate the capability of future pulsar-timing-array experiments in detecting the unknown objects. We expect that the future pulsar timing data can limit the unknown massive objects in the Solar system to be lighter than $10^{-11}$ to $10^{-12}$ $M_{\odot}$, or measure the mass of Jovian system to fractional precision of $10^{-8}$ to $10^{-9}$.


title: Black Hole Metric: Overcoming the PageRank Normalization Problem

abstract: In network science, there is often the need to sort the graph nodes. While the sorting strategy may be different, in general sorting is performed by exploiting the network structure. In particular, the metric PageRank has been used in the past decade in different ways to produce a ranking based on how many neighbors point to a specific node. PageRank is simple, easy to compute and effective in many applications, however it comes with a price: as PageRank is an application of the random walker, the arc weights need to be normalized. This normalization, while necessary, introduces a series of unwanted side-effects. In this paper, we propose a generalization of PageRank named Black Hole Metric which mitigates the problem. We devise a scenario in which the side-effects are particularily impactful on the ranking, test the new metric in both real and synthetic networks, and show the results.


title: Bifurcations, robustness and shape theory of attractors of discrete dynamical systems

abstract: We study in this paper global properties, mainly of topological nature, of attractors of discrete dynamical systems. We consider the Andronov-Hopf bifurcation for homeomorphisms of the plane and establish some robustness properties for attractors of such homeomorphisms. We also give relations between attractors of flows and quasi-attractors of homeomorphisms in Rn. Finally, we give a result on the shape (in the sense of Borsuk) of invariant sets of IFSs on the plane, and make some remarks about the recent theory of Conley attractors for IFS.


title: Remarks on hypergeometric Cauchy numbers

abstract: Hypergeometric numbers can be recognized as one of the most natural extensions of the classical Cauchy numbers in terms of determinants, though many kinds of generalizations of the Cauchy numbers have been considered by many authors. In addition, there are some relations between the hypergeometric Cauchy numbers and the classical Cauchy numbers. In this paper, we give the determinant expressions of hypergeometric Cauchy numbers and their generalizations, and show some interesting expressions of hypergeometric Cauchy numbers.


title: Transport of polymer particles in a oil-water flow in porous media: enhancing oil recovery

abstract: We study a heuristic, core-scale model for the transport of polymer particles in a two phase (oil and water) porous medium. We are motivated by recent experimental observations which report increased oil recovery when polymers are injected after the initial waterflood. The recovery mechanism is believed to be microscopic diversion of the flow, where injected particles can accumulate in narrow pore throats and clog it, in a process known as a log-jamming effect. The blockage of the narrow pore channels lead to a microscopic diversion of the water flow, causing a redistribution of the local pressure, which again can lead to the mobilization of trapped oil, enhancing its recovery. Our objective herein is to develop a core-scale model that is consistent with the observed production profiles. We show that previously obtained experimental results can be qualitatively explained by a simple two-phase flow model with an additional transport equation for the polymer particles. A key aspect of the formulation is that the microscopic heterogeneity of the rock and a dynamic altering of the permeability must be taken into account in the rate equations.


title: Computational Image Enhancement for Frequency Modulated Continuous Wave (FMCW) THz Image

abstract: In this paper, a novel method to enhance Frequency Modulated Continuous Wave (FMCW) THz imaging resolution beyond its diffraction limit is proposed. Our method comprises two stages. Firstly, we reconstruct the signal in depth-direction using a sinc-envelope, yielding a significant improvement in depth estimation and signal parameter extraction. The resulting high precision depth estimate is used to deduce an accurate reflection intensity THz image. This image is fed in the second stage of our method to a 2D blind deconvolution procedure, adopted to enhance the lateral THz image resolution beyond the diffraction limit. Experimental data acquired with a FMCW system operating at 577 GHz with a bandwidth of 126 GHz shows that the proposed method enhances the lateral resolution by a factor of 2.29 to 346.2um with respect to the diffraction limit. The depth accuracy is 91um. Interestingly, the lateral resolution enhancement achieved with this blind deconvolution concept leads to better results in comparison to conventional gaussian deconvolution. Experimental data on a PCB resolution target is presented, in order to quantify the resolution enhancement and to compare the performance with established image enhancement approaches. The presented technique allows exposure of the interwoven fibre reinforced embedded structures of the PCB test sample.


title: Point systems in Games for Health: A bibliometric scoping study

abstract: Very few details about point systems used in games for health are reported in scientific literature. To shed some light on this topic a bibliometric study, analyzing the papers containing terms related to games for health and point systems was performed and a mini taxonomy was derived. The search string game* AND health AND (point* OR score) AND system* in a Scopus bibliographic database was used to produce the corpus. We limited the search to articles, reviews and conference papers written in English and to topics related to medical, health and social subjects. The corpus papers abstracts and titles were analysed by VOSviewer and a scientific landscape was generated. The search resulted in a corpus consisting of 354 papers. The derived taxonomy contains three objects; video games, serious games and educational games. The biblimetric mapping and taxonomy revealed some interesting conclusions: (1) the video games have mostly negative effects on health, (2) the serious games might have both a direct positive health effects on users and also indirect effects by improved competencies of health professionals, and (3) the research is concerned not only to computer based educational games, but also to traditional table games and sporting games. Based on the derived taxonomy we can conclude that point systems should reward physical activity and healthy living style and punish sedentary activities.


title: Documentation of TauSpinner algorithms -- program for simulating spin effects in tau-lepton production at LHC

abstract: The tau-lepton plays an important role in the physics program at the Large Hadron Collider (LHC). It offers a probe in searches for New Physics and to measure parameters of the Standard Model. Spin of tau lepton is an interesting phenomenological quantity for separation of signal from background or in measuring properties of particles. A proper treatment of tau spin effects in the Monte Carlo simulations is important for understanding the detector acceptance as well as for the measurements/use of tau polarization and tau spin correlations. The TauSpinner is a tool to modify tau$spin effects in any sample. Also production matrix elements can be modified with its help. The general principle of TauSpinner algorithm rely on the kinematic of outgoing particles (or jets) only, and to average over all possible initial states accordingly to assumed parton distributions and cross-sections. Incoming on mass-shell parton four-momenta are reconstructed in approximation, but energy-momentum conservation is obeyed. Samples of events featuring tau leptons with their decay products, can be used as an input. The information on the polarization and spin correlations, and production/decay matrix elements is reconstructed from the kinematics of the tau lepton(s) (also nu_tau in case of W-mediated processes) and tau decay products (also jets). No other information is required. By calculating spin (or production/decay matrix element) weights, attributed to each event, it enables evaluation of modifications without regenerating events. Elements of program functionalities together with explanation of the use are given in Appendices. Present publication collect program documentation, previously scattered over several publiications.


title: Shape invariance and equivalence relations for pseudowronskians of Laguerre and Jacobi polynomials

abstract: In a previous paper we derived equivalence relations for pseudo-Wronskian determinants of Hermite polynomials. In this paper we obtain the analogous result for Laguerre and Jacobi polynomials. The equivalence formulas are richer in this case since rational Darboux transformations can be defined for four families of seed functions, as opposed to only two families in the Hermite case. The pseudo-Wronskian determinants of Laguerre and Jacobi type will thus depend on two Maya diagrams, while Hermite pseudo-Wronskians depend on just one Maya diagram. We show that these equivalence relations can be interpreted as the general transcription of shape invariance and specific discrete symmetries acting on the parameters of the isotonic oscillator and Darboux-Poschl-Teller potential.


title: Supermixed labyrinth fractals

abstract: Labyrinth fractals are dendrites in the unit square. They were introduced and studied in the last decade first in the self-similar case [Cristea & Steinsky (2009,2011)], then in the mixed case [Cristea & Steinsky (2017), Cristea & Leobacher (2017)]. Supermixed fractals constitute a significant generalisation of mixed labyrinth fractals: each step of the iterative construction is done according to not just one labyrinth pattern, but possibly to several different patterns. In this paper we introduce and study supermixed labyrinth fractals and the corresponding prefractals, called supermixed labyrinth sets, with focus on the aspects that were previously studied for the self-similar and mixed case: topological properties and properties of the arcs between points in the fractal. The facts and formulae found here extend results proven in the above mentioned cases. One of the main results is a sufficient condition for infinite length of arcs in mixed labyrinth fractals.


title: Geometric Properties of Bessel function derivatives

abstract: In this paper our aim is to find the radii of starlikeness and convexity of Bessel function derivatives for three different kind of normalization. The key tools in the proof of our main results are the Mittag-Leffler expansion for nth derivative of Bessel function and properties of real zeros of it. In addition, by using the Euler-Rayleigh inequalities we obtain some tight lower and upper bounds for the radii of starlikeness and convexity of order zero for the normalized nth derivative of Bessel function. The main results of the paper are natural extensions of some known results on classical Bessel functions of the first kind.


title: Spontaneous symmetry breaking of domain walls in phase-competing regions

abstract: We study the nature of domain walls in an ordered phase in the phase-competing region of two Ising-type order parameters. Considering a two-component $\phi^4$ theory, we show that the domain wall of the ground-state (primary) order parameter shows a second-order phase transition associated with the secondary order parameter of the competing phase; the effective theory of the phase transition is given by the Landau theory of Ising-type phase transition. We find that the phase boundary of this phase transition is different from the spinodal line of the competing order. Experimentally, the phase transition is detected by the divergence of the susceptibility corresponding to the secondary order when the temperature is quenched to introduce the domain walls.


title: Estimates for a general fractional relaxation equation and application to an inverse source problem

abstract: A general fractional relaxation equation is considered with a convolutional derivative in time introduced by A. Kochubei (Integr. Equ. Oper. Theory 71 (2011), 583-600). This equation generalizes the single-term, multi-term and distributed-order fractional relaxation equations. The fundamental and the impulse-response solutions are studied in detail. Properties such as analyticity and subordination identities are established and employed in the proof of an upper and a lower bound. The obtained results extend some known properties of the Mittag-Leffler functions. As an application of the estimates, uniqueness and conditional stability are established for an inverse source problem for the general time-fractional diffusion equation on a bounded domain.


title: Massivizing Computer Systems: a Vision to Understand, Design, and Engineer Computer Ecosystems through and beyond Modern Distributed Systems

abstract: Our society is digital: industry, science, governance, and individuals depend, often transparently, on the inter-operation of large numbers of distributed computer systems. Although the society takes them almost for granted, these computer ecosystems are not available for all, may not be affordable for long, and raise numerous other research challenges. Inspired by these challenges and by our experience with distributed computer systems, we envision Massivizing Computer Systems, a domain of computer science focusing on understanding, controlling, and evolving successfully such ecosystems. Beyond establishing and growing a body of knowledge about computer ecosystems and their constituent systems, the community in this domain should also aim to educate many about design and engineering for this domain, and all people about its principles. This is a call to the entire community: there is much to discover and achieve.


title: 2D materials coupled to hybrid metal-dielectric waveguides for THz technology

abstract: In this letter, we propose hybrid metal-dielectric waveguides coupled to 2D materials that provide strong light-matter interaction at THz frequencies. We investigate the properties of the fundamental propagating modes and show that the strength of in-plane electric field components is maximized at the top of the dielectric strip on which the 2D material is deposited. Our simulation predicts 100 % modulation of THz light by tuning the Fermi level of a graphene sheet deposited onto a 1mm-long waveguide. We also show the potential of graphene multilayers coupled to these waveguides for achieving lasing at THz frequency. Our approach is compatible with CMOS or THz quantum cascade laser technologies.


title: Nonlinear characterisation of a silicon integrated Bragg waveguide filter

abstract: Bragg waveguides are promising optical filters for pump suppression in spontaneous Four-Wave Mixing (FWM) photon sources. In this work, we investigate the generation of unwanted photon pairs in the filter itself. We do this by taking advantage of the relation between spontaneous and classical FWM, which allows for the precise characterisation of the nonlinear response of the device. The pair generation rate estimated from the classical measurement is compared with the theoretical value calculated by means of a full quantum model of the filter, which also allows to investigate the spectral properties of the generated pairs. We find a good agreement between theory and experiment, confirming that stimulated FWM is a valuable approach to characterise the nonlinear response of an integrated filter, and that the pairs generated in a Bragg waveguide are not a serious issue for the operation of a fully integrated nonclassical source.


title: Digital Cultural Heritage imaging via osmosis filtering

abstract: In Cultural Heritage (CH) imaging, data acquired within different spectral regions are often used to inspect surface and sub-surface features. Due to the experimental setup, these images may suffer from intensity inhomogeneities, which may prevent conservators from distinguishing the physical properties of the object under restoration. Furthermore, in multi-modal imaging, the transfer of information between one modality to another is often used to integrate image contents. In this paper, we apply the image osmosis model proposed in (Weickert et al. 2013) to solve similar problems arising when using diagnostic CH imaging techniques based on reflectance, emission and fluorescence mode in the optical and thermal range. For an efficient computation, we use stable operator splitting techniques. We test our methods on real artwork datasets: the thermal measurements of the mural painting "Monocromo" by Leonardo Da Vinci, the UV-VIS-IR imaging of an ancient Russian icon and the Archimedes Palimpsest dataset.


title: Second order optimality on orthogonal Stiefel manifolds

abstract: The main tool to study a second order optimality problem is the Hessian operator associated to the cost function that defines the optimization problem. By regarding an orthogonal Stiefel manifold as a constraint manifold embedded in an Euclidean space we obtain a concise matrix formula for the Hessian of a cost function defined on such a manifold. We introduce an explicit local frame on an orthogonal Stiefel manifold in order to compute the components of the Hessian matrix of a cost function. We present some important properties of this frame. As applications we rediscover second order conditions of optimality for the Procrustes and the Penrose regression problems (previously found in the literature). For the Brockett problem we find necessary and sufficient conditions for a critical point to be a local minimum. Since many optimization problems are approached using numerical algorithms, we give an explicit description of the Newton algorithm on orthogonal Stiefel manifolds.


title: Chromatic symmetric functions via the group algebra of $S_n$

abstract: We prove some Schur positivity results for the chromatic symmetric function $X_G$ of a (hyper)graph $G$, using connections to the group algebra of the symmetric group. The first such connection works for (hyper)forests $F$: we describe the Schur coefficients of $X_F$ in terms of eigenvalues of a product of Hermitian idempotents in the group algebra, one factor for each edge (a more general formula of similar shape holds for all chordal graphs). Our main application of this technique is to prove a conjecture of Taylor on the Schur positivity of certain $X_F$, which implies Schur positivity of the formal group laws associated to various combinatorial generating functions. We also introduce the pointed chromatic symmetric function $X_{G,v}$ associated to a rooted graph $(G,v)$. We prove that if $X_{G,v}$ and $X_{H,w}$ are positive in the generalized Schur basis of Strahov, then the chromatic symmetric function of the wedge sum of $(G,v)$ and $(H,w)$ is Schur positive.


title: Smooth heaps and a dual view of self-adjusting data structures

abstract: We present a new connection between self-adjusting binary search trees (BSTs) and heaps, two fundamental, extensively studied, and practically relevant families of data structures. Roughly speaking, we map an arbitrary heap algorithm within a natural model, to a corresponding BST algorithm with the same cost on a dual sequence of operations (i.e. the same sequence with the roles of time and key-space switched). This is the first general transformation between the two families of data structures. There is a rich theory of dynamic optimality for BSTs (i.e. the theory of competitiveness between BST algorithms). The lack of an analogous theory for heaps has been noted in the literature. Through our connection, we transfer all instance-specific lower bounds known for BSTs to a general model of heaps, initiating a theory of dynamic optimality for heaps. On the algorithmic side, we obtain a new, simple and efficient heap algorithm, which we call the smooth heap. We show the smooth heap to be the heap-counterpart of Greedy, the BST algorithm with the strongest proven and conjectured properties from the literature, widely believed to be instance-optimal. Assuming the optimality of Greedy, the smooth heap is also optimal within our model of heap algorithms. As corollaries of results known for Greedy, we obtain instance-specific upper bounds for the smooth heap, with applications in adaptive sorting. Intriguingly, the smooth heap, although derived from a non-practical BST algorithm, is simple and easy to implement (e.g. it stores no auxiliary data besides the keys and tree pointers). It can be seen as a variation on the popular pairing heap data structure, extending it with a "power-of-two-choices" type of heuristic.


title: Admissible Time Series Motif Discovery with Missing Data

abstract: The discovery of time series motifs has emerged as one of the most useful primitives in time series data mining. Researchers have shown its utility for exploratory data mining, summarization, visualization, segmentation, classification, clustering, and rule discovery. Although there has been more than a decade of extensive research, there is still no technique to allow the discovery of time series motifs in the presence of missing data, despite the well-documented ubiquity of missing data in scientific, industrial, and medical datasets. In this work, we introduce a technique for motif discovery in the presence of missing data. We formally prove that our method is admissible, producing no false negatives. We also show that our method can piggy-back off the fastest known motif discovery method with a small constant factor time/space overhead. We will demonstrate our approach on diverse datasets with varying amounts of missing data


title: Beta - and gamma-decay transition rates in the two-group shell model

abstract: In the framework of the two-group configuration model we obtain formulas for the reduced transition rates for beta- and gamma-transitions in even-even, odd-odd, even-odd, and odd-even nuclei. We explored dependencies of the transition rates on the occupancies of the involved subshells, as well as on the spin values of the initial and final states. The obtained formulas are useful for the qualitative spectroscopic analysis of experimental data, particulary in the regions of magicity, including the regions of the "remote" nuclei.


title: Gauge invariant determination of charged hadron masses

abstract: In this paper we show, for the first time, that charged-hadron masses can be calculated on the lattice without relying on gauge fixing at any stage of the calculations. In our simulations we follow a recent proposal and formulate full QCD+QED on a finite volume, without spoiling locality, by imposing C-periodic boundary conditions in the spatial directions. Electrically charged states are interpolated with a class of operators, originally suggested by Dirac and built as functionals of the photon field, that are invariant under local gauge transformations. We show that the quality of the numerical signal of charged-hadron masses is the same as in the neutral sector and that charged-neutral mass splittings can be calculated with satisfactory accuracy in this setup. We also discuss how to describe states of charged hadrons with real photons in a fully gauge-invariant way by providing a first evidence that the proposed strategy can be numerically viable.


title: Robust and sparse Gaussian graphical modeling under cell-wise contamination

abstract: Graphical modeling explores dependences among a collection of variables by inferring a graph that encodes pairwise conditional independences. For jointly Gaussian variables, this translates into detecting the support of the precision matrix. Many modern applications feature high-dimensional and contaminated data that complicate this task. In particular, traditional robust methods that down-weight entire observation vectors are often inappropriate as high-dimensional data may feature partial contamination in many observations. We tackle this problem by giving a robust method for sparse precision matrix estimation based on the $\gamma$-divergence under a cell-wise contamination model. Simulation studies demonstrate that our procedure outperforms existing methods especially for highly contaminated data.


title: Quantum Walks of kicked Bose-Einstein condensates

abstract: We analytically investigate the recently proposed and implemented discrete-time quantum walk based on kicked ultra-cold atoms. We show how the internal level structure of the kicked atoms leads to the emergence of a relative light-shift phase immediately relevant for the experimental realization. Analytical solutions are provided for the momentum distribution for both the case of quantum resonance and the near-resonant quasimomenta.


title: Approximate quantum Markov chains

abstract: This book is an introduction to quantum Markov chains and explains how this concept is connected to the question of how well a lost quantum mechanical system can be recovered from a correlated subsystem. To achieve this goal, we strengthen the data-processing inequality such that it reveals a statement about the reconstruction of lost information. The main difficulty in order to understand the behavior of quantum Markov chains arises from the fact that quantum mechanical operators do not commute in general. As a result we start by explaining two techniques of how to deal with non-commuting matrices: the spectral pinching method and complex interpolation theory. Once the reader is familiar with these techniques a novel inequality is presented that extends the celebrated Golden-Thompson inequality to arbitrarily many matrices. This inequality is the key ingredient in understanding approximate quantum Markov chains and it answers a question from matrix analysis that was open since 1973, i.e., if Lieb's triple matrix inequality can be extended to more than three matrices. Finally, we carefully discuss the properties of approximate quantum Markov chains and their implications.


title: Enhanced non-Markovian behavior in quantum walks with Markovian disorder

abstract: Non-Markovian quantum effects are typically observed in systems interacting with structured reservoirs. Discrete-time quantum walks are prime example of such systems in which, quantum memory arises due to the controlled interaction between the coin and position degrees of freedom. Here we show that the information backflow that quantifies memory effects can be enhanced when the particle is subjected to uncorrelated static or dynamic disorder. The presence of disorder in the system leads to localization effects in 1-dimensional quantum walks. We shown that it is possible to infer about the nature of localization in position space by monitoring the information backflow in the reduced system. Further, we study other useful properties of quantum walk such as entanglement, interference and its connection to quantum non-Markovianity.


title: Voros symbols as cluster coordinates

abstract: We show that the Borel sums of the Voros symbols considered in the theory of exact WKB analysis arise naturally as Fock-Goncharov coordinates of framed $PGL_2(\mathbb{C})$-local systems on a marked bordered surface. Using this result, we show that these Borel sums can be meromorphically continued to any point of $\mathbb{C}^*$, and we prove an asymptotic property of the monodromy map introduced in collaboration with Tom Bridgeland.


title: Evolution of Images with Diversity and Constraints Using a Generator Network

abstract: Evolutionary search has been extensively used to generate artistic images. Raw images have high dimensionality which makes a direct search for an image challenging. In previous work this problem has been addressed by using compact symbolic encodings or by constraining images with priors. Recent developments in deep learning have enabled a generation of compelling artistic images using generative networks that encode images with lower-dimensional latent spaces. To date this work has focused on the generation of images concordant with one or more classes and transfer of artistic styles. There is currently no work which uses search in this latent space to generate images scoring high or low aesthetic measures. In this paper we use evolutionary methods to search for images in two datasets, faces and butterflies, and demonstrate the effect of optimising aesthetic feature scores in one or two dimensions. The work gives a preliminary indication of which feature measures promote the most interesting images and how some of these measures interact.


title: Chaos in Kuramoto oscillator networks

abstract: Kuramoto oscillators are widely used to explain collective phenomena in networks of coupled oscillatory units. We show that simple networks of two populations with a generic coupling scheme, where both coupling strengths and phase lags between and within populations are distinct, can exhibit chaotic dynamics as conjectured by Ott and Antonsen [Chaos, 18, 037113 (2008)]. These chaotic mean-field dynamics arise universally across network size, from the continuum limit of infinitely many oscillators down to very small networks with just two oscillators per population. Hence, complicated dynamics are expected even in the simplest description of oscillator networks.


title: DMRG simulations of SU(N) Heisenberg chains using standard Young tableaux: fundamental representation and comparison with finite-size Bethe ansatz

abstract: We develop an efficient method to perform density matrix renormalization group simulations of the SU(N) Heisenberg chain with open boundary conditions taking full advantage of the SU(N) symmetry of the problem. This method is an extension of the method previously developed for exact diagonalizations and relies on a systematic use of the basis of standard Young tableaux. Concentrating on the model with the fundamental representation at each site (i.e. one particle per site in the fermionic formulation), we have benchmarked our results for the ground state energy up to N = 8 and up to 420 sites by comparing them with Bethe ansatz results on open chains, for which we have derived and solved the Bethe ansatz equations. The agreement for the ground state energy is excellent for SU(3) (12 digits). It decreases with N , but it is still satisfactory for N = 8 (6 digits). Central charges c are also extracted from the entanglement entropy using the Calabrese-Cardy formula, and agree with the theoretical values expected from the SU(N)1 Wess-Zumino-Witten CFTs.


title: Extensive study of electron acceleration by relativistic surface plasmons

abstract: The excitation of surface plasmons with ultra-intense ($I\sim 5\times 10^{19}$ W/cm$^2$), high contrast ($\sim 10^{12}$) laser pulses on periodically-modulated solid targets has been recently demonstrated to produce collimated bunches of energetic electrons along the target surface [Fedeli et al., Phys. Rev. Lett. 116, 5001 (2016)]. Here we report an extensive experimental and numerical study aimed to a complete characterization of the acceleration mechanism, demonstrating its robustness and promising characteristics for an electron source. By comparing different grating structures, we identify the relevant parameters to optimize the acceleration and obtain bunches of $\sim 650$ pC of charge at several MeV of energy with blazed gratings.


title: On the P vs NP question: a proof of inequality

abstract: The analysis discussed in this paper is based on a well-known NP-complete problem which is called satisfiability problem or SAT. From SAT a new NP-complete problem is derived, which is described by a Boolean function called core function. In this paper it is proved that the cost of the minimal implementation of core function increases with n exponentially. Since the synthesis of core function is an NP-complete problem, this result is equivalent to proving that P and NP do not coincide.


title: Electron drag in ferromagnetic structures separated by an insulating interface

abstract: We consider electron drag in a system of two ferromagnetic layers separated by an insulating interface. The source of it is expected to be magnon-electron interactions. Namely, we assume that the external voltage is applied to the "active" layer stimulating electric current through this layer. In its turn, the scattering of the current-carrying electrons by magnons leads to a magnon drag current within this layer. The 3-magnons interactions between magnons in the two layers (being of non-local nature) lead to magnon drag within the "passive" layer which, correspondingly, produce electron drag current via processes of magnon-electron scattering. We estimate the drag current and compare it to the phonon-induced one.


title: An autonomous single-piston engine with a quantum rotor

abstract: Pistons are elementary components of a wide variety of thermal engines, allowing to convert input fuel into rotational motion. Here, we propose a single-piston engine where the rotational degree of freedom is effectively realized by the flux of a Josephson loop -- a quantum rotor -- while the working volume corresponds to the effective length of a superconducting resonator. Our autonomous design implements a Carnot cycle, relies solely on standard thermal baths and can be implemented with circuit quantum electrodynamics. We demonstrate how the engine is able to extract a net positive work via its built-in synchronicity using a filter cavity as an effective valve, eliminating the need for external control.


title: Shannon entropy and particle decays

abstract: We deploy Shannon's information entropy to the distribution of branching fractions in a particle decay. This serves to quantify how important a given new reported decay channel is, from the point of view of the information that it adds to the already known ones. Because the entropy is additive, one can subdivide the set of channels and discuss, for example, how much information the discovery of a new decay branching would add; or subdivide the decay distribution down to the level of individual quantum states (which can be quickly counted by the phase space). We illustrate the concept with some examples of experimentally known particle decay distributions.


title: Phase and micromotion of Bose-Einstein condensates in a time-averaged ring trap

abstract: Rapidly scanning magnetic and optical dipole traps have been widely utilised to form time-averaged potentials for ultracold quantum gas experiments. Here we theoretically and experimentally characterise the dynamic properties of Bose-Einstein condensates in ring-shaped potentials that are formed by scanning an optical dipole beam in a circular trajectory. We find that unidirectional scanning leads to a non-trivial phase profile of the condensate that can be approximated analytically using the concept of phase imprinting. While the phase profile is not accessible through in-trap imaging, time-of-flight expansion manifests clear density signatures of an in-trap phase step in the condensate, coincident with the instantaneous position of the scanning beam. The phase step remains significant even when scanning the beam at frequencies two orders of magnitude larger than the characteristic frequency of the trap. We map out the phase and density properties of the condensate in the scanning trap, both experimentally and using numerical simulations, and find excellent agreement. Furthermore, we demonstrate that bidirectional scanning eliminated the phase gradient, rendering the system more suitable for coherent matter wave interferometry.


title: Optimal control of the customer dynamics based on marketing policy

abstract: We consider an optimal control problem for a non-autonomous model of ODEs that describes the evolution of the number of customers in some firm. Namely we study the best marketing strategy. Considering a $L^2$ cost functional, we establish the existence and uniqueness of optimal solutions, using an inductive argument to obtain uniqueness on the whole interval from local uniqueness. We also present some simulation results, based on our model, and compare them with results we obtain for an $L^1$ cost functional. For the $L^1$ cost functional the optimal solutions are of bang-bang type and thus easier to implement, because at every moment possible actions are chosen from a finite set of possibilities. For the autonomous case of $L^2$ problem, we show the effectiveness of the optimal control strategy against other formulations of the problem with simpler controls.


title: Grammar-based Compression of Unranked Trees

abstract: We introduce forest straight-line programs (FSLPs) as a compressed representation of unranked ordered node-labelled trees. FSLPs are based on the operations of forest algebra and generalize tree straight-line programs. We compare the succinctness of FSLPs with two other compression schemes for unranked trees: top dags and tree straight-line programs of first-child/next sibling encodings. Efficient translations between these formalisms are provided. Finally, we show that equality of unranked trees in the setting where certain symbols are associative or commutative can be tested in polynomial time. This generalizes previous results for testing isomorphism of compressed unordered ranked trees.


title: Biorthogonal systems on unit interval and zeta dilation operators

abstract: An elementary 'quantum-mechanical' derivation of the conditions for a system of functions to form a Reisz basis of a Hilbert space on a finite interval is presented.


title: Measurement of the Branching Fraction For the Semi-Leptonic Decays $D^{0(+)}\to \pi^{-(0)}\mu^+\nu_\mu$ and Test of Lepton Flavor Universality

abstract: Using a data sample corresponding to an integrated luminosity of $2.93\,\rm fb^{-1}$ taken at a center-of-mass energy of 3.773\,GeV with the BESIII detector operated at the BEPCII collider, we perform an analysis of the semi-leptonic decays $D^{0(+)}\to \pi^{-(0)}\mu^+\nu_\mu$. The branching fractions of $D^0\to \pi^-\mu^+\nu_\mu$ and $D^+\to \pi^0\mu^+\nu_\mu$ are measured to be $(0.272 \pm 0.008_{\rm stat.} \pm 0.006_{\rm syst.})\%$ and $(0.350 \pm 0.011_{\rm stat.} \pm 0.010_{\rm syst.})\%$, respectively, where the former is of much improved precision compared to previous results and the latter is determined for the first time. Using these results along with previous BESIII measurements of $D^{0(+)}\to \pi^{-(0)}e^+\nu_e$, we calculate the branching fraction ratios to be ${\mathcal R}^0\equiv {\mathcal B}_{D^{0}\to \pi^{-}\mu^+\nu_\mu}/{\mathcal B}_{D^{0}\to \pi^{-}e^+\nu_e}=0.922\pm 0.030_{\rm stat.}\pm0.022_{\rm syst.}$ and ${\mathcal R}^+\equiv {\mathcal B}_{D^{+}\to \pi^{0}\mu^+\nu_\mu}/{\mathcal B}_{D^{+}\to \pi^{0}e^+\nu_e}=0.964\pm 0.037_{\rm stat.}\pm0.026_{\rm syst.}$, which are compatible with the theoretical expectation of lepton flavor universality within $1.7\sigma$ and $0.5\sigma$, respectively. We also examine the branching fraction ratios in different four-momentum transfer square regions, and find no significant deviations from the standard model predictions.


title: A (2+1)-dimensional Anisotropic KPZ growth model with a smooth phase

abstract: Stochastic growth processes in dimension $(2+1)$ were conjectured by D. Wolf, on the basis of renormalization-group arguments, to fall into two distinct universality classes, according to whether the Hessian $H_\rho$ of the speed of growth $v(\rho)$ as a function of the average slope $\rho$ satisfies $\det H_\rho>0$ ("isotropic KPZ class") or $\det H_\rho\le 0$ ("anisotropic KPZ (AKPZ)" class). The former is characterized by strictly positive growth and roughness exponents, while in the AKPZ class fluctuations are logarithmic in time and space. It is natural to ask (a) if one can exhibit interesting growth models with "smooth" stationary states, i.e., with $O(1)$ fluctuations (instead of logarithmically or power-like growing, as in Wolf's picture) and (b) what new phenomena arise when $v(\cdot)$ is not smooth, so that $H_\rho$ is not defined. The two questions are actually related and here we provide an answer to both, in a specific framework. We define a $(2+1)$-dimensional interface growth process, based on the so-called shuffling algorithm for domino tilings. The stationary, non-reversible measures are translation-invariant Gibbs measures on perfect matchings of $\mathbb Z^2$, with $2$-periodic weights. If $\rho\ne0$, fluctuations are known to grow logarithmically in space and to behave like a two-dimensional GFF. We prove that fluctuations grow at most logarithmically in time and that $\det H_\rho<0$: the model belongs to the AKPZ class. When $\rho=0$, instead, the stationary state is "smooth", with correlations uniformly bounded in space and time; correspondingly, $v(\cdot)$ is not differentiable at $\rho=0$ and we extract the singularity of the eigenvalues of $H_\rho$ for $\rho\sim 0$.


title: Non-idempotent types for classical calculi in natural deduction style

abstract: In the first part of this paper, we define two resource aware typing systems for the {\lambda}{\mu}-calculus based on non-idempotent intersection and union types. The non-idempotent approach provides very simple combinatorial arguments-based on decreasing measures of type derivations-to characterize head and strongly normalizing terms. Moreover, typability provides upper bounds for the lengths of the head reduction and the maximal reduction sequences to normal-form. In the second part of this paper, the {\lambda}{\mu}-calculus is refined to a small-step calculus called {\lambda}{\mu}s, which is inspired by the substitution at a distance paradigm. The {\lambda}{\mu}s-calculus turns out to be compatible with a natural extensionof the non-idempotent interpretations of {\lambda}{\mu}, i.e., {\lambda}{\mu}s-reduction preserves and decreases typing derivations in an extended appropriate typing system. We thus derive a simple arithmetical characterization of strongly {\lambda}{\mu}s-normalizing terms by means of typing.


title: How Much Data Do You Need? An Operational, Pre-Asymptotic Metric for Fat-tailedness

abstract: This note presents an operational measure of fat-tailedness for univariate probability distributions, in $[0,1]$ where 0 is maximally thin-tailed (Gaussian) and 1 is maximally fat-tailed. Among others,1) it helps assess the sample size needed to establish a comparative $n$ needed for statistical significance, 2) allows practical comparisons across classes of fat-tailed distributions, 3) helps understand some inconsistent attributes of the lognormal, pending on the parametrization of its scale parameter. The literature is rich for what concerns asymptotic behavior, but there is a large void for finite values of $n$, those needed for operational purposes. Conventional measures of fat-tailedness, namely 1) the tail index for the power law class, and 2) Kurtosis for finite moment distributions fail to apply to some distributions, and do not allow comparisons across classes and parametrization, that is between power laws outside the Levy-Stable basin, or power laws to distributions in other classes, or power laws for different number of summands. How can one compare a sum of 100 Student T distributed random variables with 3 degrees of freedom to one in a Levy-Stable or a Lognormal class? How can one compare a sum of 100 Student T with 3 degrees of freedom to a single Student T with 2 degrees of freedom? We propose an operational and heuristic measure that allow us to compare $n$-summed independent variables under all distributions with finite first moment. The method is based on the rate of convergence of the Law of Large numbers for finite sums, $n$-summands specifically. We get either explicit expressions or simulation results and bounds for the lognormal, exponential, Pareto, and the Student T distributions in their various calibrations --in addition to the general Pearson classes.


title: Extremals For Fractional Order Hardy-Sobolev-Maz'ya Inequality

abstract: In this article, we derive the existence of positive solutions of a semi-linear, non-local elliptic PDE, involving a singular perturbation of the fractional laplacian, coming from the fractional Hardy-Sobolev-Maz'ya inequality, derived in this paper. We also derive symmetry properties and a precise asymptotic behaviour of solutions.


title: Nano-Opto-Electro-Mechanical Systems

abstract: A new class of hybrid systems that couple optical, electrical and mechanical degrees of freedom in nanoscale devices is under development in laboratories worldwide. These nano-opto-electro-mechanical systems (NOEMS) offer unprecedented opportunities to dynamically control the flow of light in nanophotonic structures, at high speed and low power consumption. Drawing on conceptual and technological advances from cavity optomechanics, they also bear the potential for highly efficient, low-noise transducers between microwave and optical signals, both in the classical and quantum domains. This Progress Article discusses the fundamental physical limits of NOEMS, reviews the recent progress in their implementation, and suggests potential avenues for further developments in this field.


title: Dynamics of high-order solitons in the nonlocal nonlinear Schr\"{o}dinger equations

abstract: A study of high-order solitons in three nonlocal nonlinear Schr\"{o}dinger equations is presented, which includes the \PT-symmetric, reverse-time, and reverse-space-time nonlocal nonlinear Schr\"{o}dinger equations. General high-order solitons in three different equations are derived from the same Riemann-Hilbert solutions of the AKNS hierarchy, except for the difference in the corresponding symmetry relations on the "perturbed" scattering data. Dynamics of general high-order solitons in these equations is further analyzed. It is shown that the high-order fundamental-soliton is always moving on several different trajectories in nearly equal velocities, and they can be nonsingular or repeatedly collapsing, depending on the choices of the parameters. It is also shown that high-order multi-solitons could have more complicated wave structures and behave very differently from high-order fundamental solitons. More interesting is the high-order hybrid-pattern solitons, which are derived from combination of different size of block matrix in the Riemann-Hilbert solutions and thus they can describe a nonlinear interaction between several types of solitons.


title: On the $L^p$ norm of the torsion unction

abstract: Bounds are obtained for the $L^p$ norm of the torsion function $v_{\Omega}$, i.e. the solution of $-\Delta v=1,\, v\in H_0^1(\Omega),$ in terms of the Lebesgue measure of $\Omega$ and the principal eigenvalue $\lambda_1(\Omega)$ of the Dirichlet Laplacian acting in $L^2(\Omega)$. We show that these bounds are sharp for $1\le p\le 2$.


title: 2PI effective action for the SYK model and tensor field theories

abstract: We discuss the two-particle irreducible (2PI) effective action for the SYK model and for tensor field theories. For the SYK model the 2PI effective action reproduces the bilocal reformulation of the model without using replicas. In general tensor field theories the 2PI formalism is the only way to obtain a bilocal reformulation of the theory, and as such is a precious instrument for the identification of soft modes and for possible holographic interpretations. We compute the 2PI action for several models, and push it up to fourth order in the $1/N$ expansion for the model proposed by Witten in arXiv:1610.09758, uncovering a one-loop structure in terms of an auxiliary bilocal action.


title: Finding small-width connected path decompositions in polynomial time

abstract: A connected path decomposition of a simple graph $G$ is a path decomposition $(X_1,\ldots,X_l)$ such that the subgraph of $G$ induced by $X_1\cup\cdots\cup X_i$ is connected for each $i\in\{1,\ldots,l\}$. The connected pathwidth of $G$ is then the minimum width over all connected path decompositions of $G$. We prove that for each fixed $k$, the connected pathwidth of any input graph can be computed in polynomial-time. This answers an open question raised by Fedor V. Fomin during the GRASTA 2017 workshop, since connected pathwidth is equivalent to the connected (monotone) node search game.


title: Sharp Lower Bounds for the First Eigenvalues of the Bi-Laplace Operator

abstract: We obtain sharp lower bounds for the first eigenvalue of four types of eigenvalue problem defined by the bi-Laplace operator on compact manifolds with boundary and determine all the eigenvalues and the corresponding eigenfunctions of a Wentzell-type bi-Laplace problem on Euclidean balls.