title: A Novel GPU-based Parallel Implementation Scheme and Performance Analysis of Robot Forward Dynamics Algorithms

abstract: We propose a novel unifying scheme for parallel implementation of articulated robot dynamics algorithms. It is based on a unified Lie group notation for deriving the equations of motion of articulated robots, where various well-known forward algorithms differ only by their joint inertia matrix inversion strategies. This new scheme leads to a unified abstraction of state-of-the-art forward dynamics algorithms into combinations of block bi-diagonal and/or block tri-diagonal systems, which may be efficiently solved by parallel all-prefix-sum operations (scan) and parallel odd-even elimination (OEE) respectively. We implement the proposed scheme on a Nvidia CUDA GPU platform for the comparative study of three algorithms, namely the hybrid articulated-body inertia algorithm (ABIA), the parallel joint space inertia inversion algorithm (JSIIA) and the constrained force algorithm (CFA), and the performances are analyzed.


title: A zero-one Law for improvements to Dirichlet's Theorem

abstract: We give an integrability condition on a function $\psi$ guaranteeing that for almost all (or almost no) $x\in\mathbb{R}$, the system $|qx-p|\leq \psi(t)$, $|q|<t$ is solvable in $p\in \mathbb{Z}$, $q\in \mathbb{Z}\smallsetminus \{0\}$ for sufficiently large $t$. Along the way, we characterize such $x$ in terms of the growth of their continued fraction entries, and we establish that Dirichlet's Approximation Theorem is sharp in a very strong sense. Higher-dimensional generalizations are discussed at the end of the paper.


title: The Validity of the Weizsacker-Williams Approximation and the Analysis of Beam Dump Experiments: Production of a New Scalar Boson

abstract: Beam dump experiments have been used to search for new particles with null results interpreted in terms of limits on masses $m_\phi$ and coupling constants $\epsilon$. However these limits have been obtained by using approximations [including the Weizs\"{a}cker-Williams (WW) approximation] or Monte-Carlo simulations. We display methods, using a new scalar boson as an example, to obtain the cross section and the resulting particle production numbers without using approximations or Monte-Carlo simulations. We show that the approximations cannot be used to obtain accurate values of cross sections. The corresponding exclusion plots differ by substantial amounts when seen on a linear scale. In the event of a discovery, we generate pseudodata (assuming given values of $m_\phi$ and $\epsilon$) in the currently allowed regions of parameter space. The use of approximations to analyze the pseudodata for the future experiments is shown to lead to considerable errors in determining the parameters. Furthermore, a new region of parameter space can be explored without using one of the common approximations, $m_\phi\gg m_e$. Our method can be used as a consistency check for Monte-Carlo simulations.


title: Deep Learning for Video Classification and Captioning

abstract: Accelerated by the tremendous increase in Internet bandwidth and storage space, video data has been generated, published and spread explosively, becoming an indispensable part of today's big data. In this paper, we focus on reviewing two lines of research aiming to stimulate the comprehension of videos with deep learning: video classification and video captioning. While video classification concentrates on automatically labeling video clips based on their semantic contents like human actions or complex events, video captioning attempts to generate a complete and natural sentence, enriching the single label as in video classification, to capture the most informative dynamics in videos. In addition, we also provide a review of popular benchmarks and competitions, which are critical for evaluating the technical progress of this vibrant field.


title: Nonparametric Bayesian Topic Modelling with the Hierarchical Pitman-Yor Processes

abstract: The Dirichlet process and its extension, the Pitman-Yor process, are stochastic processes that take probability distributions as a parameter. These processes can be stacked up to form a hierarchical nonparametric Bayesian model. In this article, we present efficient methods for the use of these processes in this hierarchical context, and apply them to latent variable models for text analytics. In particular, we propose a general framework for designing these Bayesian models, which are called topic models in the computer science community. We then propose a specific nonparametric Bayesian topic model for modelling text from social media. We focus on tweets (posts on Twitter) in this article due to their ease of access. We find that our nonparametric model performs better than existing parametric models in both goodness of fit and real world applications.


title: Multisource thermal model to the transverse momentum spectra in pp collisions at RHIC and LHC energies

abstract: In this paper, an improved multi-source thermal model is used to analyze the transverse momentum spectra in $pp$ collisions at high energies ranging from $\sqrt{\mathrm{\it s_{NN}}}$ = 62.4 GeV to 7 TeV. We give a detailed comparison between the theoretical results and experimental data at RHIC and LHC energies. It is shown that the excitation factors of emission sources depend linearly on ln$\sqrt{\mathrm{\it s_{NN}}}$ in the framework. Based on the variation regularity of the source-excitation factors, transverse momentum spectra are predicted in $pp$ collisions at higher energies, potential future $pp$ colliders operating at $\sqrt{\mathrm{\it s_{NN}}}$ = 33 and 100 TeV.


title: Homotopy theory of monoid actions via group actions and an Elmendorf style theorem

abstract: Let $M$ be a monoid and $G:\mathbf{Mon} \to \mathbf{Grp}$ be the group completion functor from monoids to groups. Given a collection $\mathcal{X}$ of submonoids of $M$ and for each $N\in \mathcal{X}$ a collection $\mathcal{Y}_N$ of subgroups of $G(N)$, we construct a model structure on the category of $M$-spaces and $M$-equivariant maps, called the $(\mathcal X,\mathcal Y)$-model structure, in which weak equivalences and fibrations are induced from the standard $\mathcal{Y}_N$-model structures on $G(N)$-spaces for all $N\in \mathcal{X}$. We also show that for a pair of collections $(\mathcal{X},\mathcal{Y})$ there is a small category ${\mathbf O}_{(\mathcal{X},\mathcal{Y})}$ whose objects are $M$-spaces $M\times_NG(N)/H$ for each $N\in \mathcal X$ and $H\in \mathcal Y_N$ and morphisms are $M$-equivariant maps, such that the $(\mathcal X,\mathcal Y)$-model structure on the category of $M$-spaces is Quillen equivalent to the projective model structure on the category of contravariant ${\mathbf O}_{(\mathcal{X},\mathcal{Y})}$-diagrams of spaces.


title: Five-loop massive tadpoles

abstract: We provide an update on a long-term project that aims at evaluating massive vacuum integrals at the five-loop frontier, with high precision and in various space-time dimensions. A number of applications are sketched, mainly concerning the determination of anomalous dimensions, for quantum field theories in four, three and two dimensions.


title: Observation of Magnetic reconnection at a 3D null point associated with a solar eruption

abstract: Magnetic null has long been recognized as a special structure serving as a preferential site for magnetic reconnection (MR). However, the direct observational study of MR at null-points is largely lacking. Here, we show the observations of MR around a magnetic null associated with an eruption that resulted in an M1.7 flare and a coronal mass ejection. The GOES X- ray profile of the flare exhibited two peaks at 02:23 UT and 02:40 UT on 2012 November 8, respectively. Based on the imaging observations, we find that the first and also primary X- ray peak was originated from MR in the current sheet underneath the erupting magnetic flux rope (MFR). On the other hand, the second and also weaker X-ray peak was caused by MR around a null-point located above the pre-eruption MFR. The interaction of the null-point and the erupting MFR can be described as a two-step process. During the first step, the erupting and fast expanding MFR passed through the null-point, resulting in a significant displacement of the magnetic field surrounding the null. During the second step, the displaced magnetic field started to move back, resulting in a converging inflow and subsequently the MR around the null. The null-point reconnection is a different process from the current sheet reconnection in this flare; the latter is the cause of the main peak of the flare, while the former is the cause of the secondary peak of the flare and the conspicuous high-lying cusp structure.


title: Modules close to SSP- and SIP-modules

abstract: In this paper, we investigate some properties of SIP, SSP and CS-Rickart modules. We give equivalent conditions for SIP and SSP modules; establish connections between the class of semisimple artinian rings and the class of SIP rings. It shows that $R$ is a semisimple artinian ring if and only if $R_R$ is SIP and every right $R$-module has a SIP-cover. We also prove that $R$ is a semiregular ring and $J(R) = Z(R_R)$ if only if every finitely generated projective module is a SIP-CS module which is also a $C2$ module.


title: Krigings Over Space and Time Based on Latent Low-Dimensional Structures

abstract: We propose a new approach to represent nonparametrically the linear dependence structure of a spatio-temporal process in terms of latent common factors. Though it is formally similar to the existing reduced rank approximation methods (Section 7.1.3 of Cressie and Wikle, 2011), the fundamental difference is that the low-dimensional structure is completely unknown in our setting, which is learned from the data collected irregularly over space but regularly over time. Furthermore a graph Laplacian is incorporated in the learning in order to take the advantage of the continuity over space, and a new aggregation method via randomly partitioning space is introduced to improve the efficiency. We do not impose any stationarity conditions over space either, as the learning is facilitated by the stationarity in time. Krigings over space and time are carried out based on the learned low-dimensional structure, which is scalable to the cases when the data are taken over a large number of locations and/or over a long time period. Asymptotic properties of the proposed methods are established. Illustration with both simulated and real data sets is also reported.


title: Gallium Bismuth Halides GaBi-X2 (X= I, Br, Cl) Monolayers with Distorted Hexagonal Framework: Novel Room-Temperature Quantum Spin Hall Insulators

abstract: Quantum Spin Hall (QSH) insulators with a large topologically nontrivial bulk gap are crucial for future applications of the QSH effect. Among these, group III-V monolayers and their halides with chair structure (regular hexagonal framework, RHF) were widely studied. Using first-principles calculations, we propose a new structure model for the functionalized group III-V monolayers, which consist of rectangular GaBi-X2 (X=I, Br, Cl) monolayers with a distorted hexagonal framework (DHF). These structures have a much lower energy than the GaBi-X2 monolayers with chair structure. Remarkably, the DHF GaBi-X2 monolayers are all QSH insulators, which exhibit sizeable nontrivial band gaps ranging from 0.17 eV to 0.39 eV. Those band gaps can be widely tuned by applying different spin-orbit coupling (SOC) strengths, resulting in a distorted Dirac cone.


title: Twitter-Network Topic Model: A Full Bayesian Treatment for Social Network and Text Modeling

abstract: Twitter data is extremely noisy -- each tweet is short, unstructured and with informal language, a challenge for current topic modeling. On the other hand, tweets are accompanied by extra information such as authorship, hashtags and the user-follower network. Exploiting this additional information, we propose the Twitter-Network (TN) topic model to jointly model the text and the social network in a full Bayesian nonparametric way. The TN topic model employs the hierarchical Poisson-Dirichlet processes (PDP) for text modeling and a Gaussian process random function model for social network modeling. We show that the TN topic model significantly outperforms several existing nonparametric models due to its flexibility. Moreover, the TN topic model enables additional informative inference such as authors' interests, hashtag analysis, as well as leading to further applications such as author recommendation, automatic topic labeling and hashtag suggestion. Note our general inference framework can readily be applied to other topic models with embedded PDP nodes.


title: Central Runge-Kutta discontinuous Galerkin methods for the special relativistic hydrodynamics

abstract: This paper developes Runge-Kutta $P^K$-based central discontinuous Galerkin (CDG) methods with WENO limiter to the one- and two-dimensional special relativistic hydrodynamical (RHD) equations, $K=1,2,3$. Different from the non-central DG methods, the \CDG{} have to find two approximate solutions defined on mutually dual meshes. For each mesh, the CDG approximate solutions on its dual mesh are used to calculate the flux values in the cell and on the cell boundary so that the approximate solutions on mutually dual meshes are coupled with each other, and the use of numerical flux may be avoided. The WENO limiter is adaptively implemented via two steps: the "troubled" cells are first identified by using a modified TVB minmod function, and then the WENO technique is used to locally reconstruct new polynomials of degree $(2K+1)$ replacing the CDG solutions inside the "troubled' cells by the cell average values of the CDG solutions in the neighboring cells as well as the original cell averages of the "troubled" cells. Because the WENO limiter is only employed for finite "troubled" cells, the computational cost can be as little as possible. The accuracy of the CDG without the numerical dissipation is analyzed and calculation of the flux integrals over the cells is also addressed. Several test problems in one and two dimensions are solved by using our \CDG{} with WENO limiter. The computations demonstrate that our methods are stable, accurate, and robust in solving complex RHD problems.


title: Injectivity of generalized Wronski maps

abstract: We study linear projections on Pluecker space whose restriction to the Grassmannian is a non-trivial branched cover. When an automorphism of the Grassmannian preserves the fibers, we show that the Grassmannian is necessarily of m-dimensional linear subspaces in a symplectic vector space of dimension 2m, and the linear map is the Lagrangian involution. The Wronski map for a self-adjoint linear differential operator and pole placement map for symmetric linear systems are natural examples.


title: Black Hole Quantum Vacuum Polarization in Higher Dimensions

abstract: The goal of this paper is to extend to higher dimensionality the methods and computations of vacuum polarization effects in black hole spacetimes. We focus our attention on the case of five dimensional Schwarzschild-Tangherlini black holes, for which we adapt the general method initially developed by Candelas and later refined by Anderson and others. We make use of point splitting regularization and of the WKB approximation to extract the divergences occuring in the coincidence limit of the Green function and, after calculating the counter-terms using the Schwinger - De Witt expansion, we explicitly prove the cancellation of the divergences and the regularity of the vacuum polarization once counter-terms are added up. We finally handle numerically the renormalized expression of the vacuum polarization. As a check on the method we also prove the regularity of the vacuum polarization in the six dimensional case in the large mass limit.


title: Particle-Filter-Enabled Real-Time Sensor Fault Detection Without a Model of Faults

abstract: We are experiencing an explosion in the amount of sensors measuring our activities and the world around us. These sensors are spread throughout the built environment and can help us perform state estimation and control of related systems, but they are often built and/or maintained by third parties or system users. As a result, by outsourcing system measurement to third parties, the controller must accept their measurements without being able to directly verify the sensors' correct operation. Instead, detection and rejection of measurements from faulty sensors must be done with the raw data only. Towards this goal, we present a method of detecting possibly faulty behavior of sensors. The method does not require that the control designer have any model of faulty sensor behavior. As we discuss, it turns out that the widely-used particle filter state estimation algorithm provides the ingredients necessary for a hypothesis test against all ranges of correct operating behavior, obviating the need for a fault model to compare measurements. We demonstrate the applicability of our method by demonstrating its ability to reject faulty measurements and improve state estimation accuracy in a nonlinear vehicle traffic model without information of generated faulty measurements' characteristics. In our test, we correctly identify nearly 90% of measurements as faulty or non-faulty without having any fault model. This leads to only a 3% increase in state estimation error over a theoretical 100%-accurate fault detector.


title: Dispersive Jaynes-Cummings Hamiltonian describing a two-level atom interacting with a two-level single mode field

abstract: We investigate the time evolution of statistical properties of a single mode radiation field after its interaction with a two-level atom. The entire system is described by a dispersive Jaynes-Cummings Hamiltonian assuming the atomic state evolving from an initial superposition of its excited and ground states, $\vert e\rangle +\vert g\rangle ,$ and the field evolving from an initial superposition of two excited levels, $\vert n_{1}\rangle+ \vert n_{2}\rangle$. It is found that the field evolution is periodic, the period depending on the ratio $n_{2}/n_{1}.$ The energy excitation oscillates between these two states and the statististics can be either sub- or super-Poissonian, depending on the values $n_{1},$ $n_{2}$.


title: Manipulation of magnetic Skyrmions with a Scanning Tunneling Microscope

abstract: The dynamics of a single magnetic Skyrmion in an atomic spin system under the influence of Scanning Tunneling Microscope is investigated by computer simulations solving the Landau-Lifshitz-Gilbert equation. Two possible scenarios are described: manipulation with aid of a spin-polarized tunneling current and by an electric field created by the scanning tunneling microscope. The dynamics during the creation and annihilation process is studied and the possibility to move single Skyrmions is showed.


title: Environment-Protected Solid State Based Distributed Charge Qubit

abstract: A novel solid state based charge qubit is presented. The system consists of a one-dimensional wire with a pair of qubits embedded at its center. It is shown that the system supports collective states localized in the left and right sides of the wire and therefore, as a whole, performs as a single qubit. The couplings between the ground and excited states of the two central qubits are inversely proportional making them fully asynchronized and allowing for coherent manipulation and gate operations. Initialization and measurement devices, such as leads and charge detectors, connected to the edges of the wire are modeled by a continuum of energy states. The coupling to the continuum is discussed using the effective non-Hermitian Hamiltonian. At weak continuum coupling, all internal states uniformly acquire small decay widths. This changes dramatically as the coupling strength increases: the width distribution undergoes a sharp restructuring and is no longer uniformly divided among the eigenstates. Two broad resonances localized at the ends of the wire are formed. These superradiant states (analogous to Dicke states in quantum optics), effectively protect the remaining internal states from decaying into the continuum and hence increase the lifetime of the qubit. Environmental noise is introduced by considering random Gaussian fluctuations of electronic energies. The interplay between decoherence and superradiance is studied by solving the stochastic Liouville equation. In addition to increasing the lifetime, the emergence of the superradiant states increases the qubit coherence.


title: Simulating single photons with realistic photon sources

abstract: Quantum information processing provides remarkable advantages over its classical counterpart. Quantum optical systems are proved to be sufficient for realizing general quantum tasks, which however often rely on single photon sources. In practice, imperfect single photon sources, such as weak coherent state source, are used instead, which will inevitably limit the power in demonstrating quantum effects. For instance, with imperfect photon sources, the key rate of the BB84 quantum key distribution protocol will be very low, which fortunately can be resolved by utilizing the decoy state method. As a generalization, we investigate an efficient way to simulate single photons with imperfect ones to an arbitrary desired accuracy when the number of photonic inputs is small. Based on this simulator, we can thus replace the tasks that involve only a few single photon inputs with the ones that only make use of imperfect photon sources. In addition, our method also provides a quantum simulator to quantum computation based on quantum optics. In the main context, we take phase randomized coherent state as an example for analysis. A general photon source applies similarly and may provide some further advantages for certain tasks.


title: Non-formality of the odd dimensional framed little balls operads

abstract: We prove that the chain operad of the framed little balls (or disks) operad is not formal as a non-symmetric operad over the rationals if the dimension of their balls is odd and greater than 4.


title: Numerical search for a potential planet sculpting the young disc of HD 115600

abstract: Radial and azimuthal features (such as disc offsets and eccentric rings) seen in high resolution images of debris discs, provide us with the unique opportunity of finding potential planetary companions which betray their presence by gravitationally sculpting such asymmetric features. The young debris disc around HD 115600, imaged recently by the Gemini Planet Imager, is such a disc with an eccentricity 0.1<e<0.2 and a projected offset from the star of 4 AU. Using our modified N-body code which incorporates radiation forces, we firstly aim to determine the orbit of a hidden planetary companion potentially responsible for shaping the disc. We run a suite of simulations covering a broad range of planetary parameters using a Monte Carlo Markov Chain sampling method and create synthetic images from which we extract the geometric disc parameters to be compared with the observed and model-derived quantities. We then repeat the study using a traditional grid to explore the planetary parameter space and aim secondly to compare the efficiency of both sampling methods. We find a planet of 7.8 Jupiter mass orbiting at 30 AU with an eccentricity of e = 0.2 to be the best fit to the observations of HD 115600. Technically, such planet has a contrast detectable by direct imaging, however the system's orientation does not favour such detection. In this study, at equal number of explored planetary configurations, the Monte Carlo Markov Chain not only converges faster but provides a better fit than a traditional grid.


title: Measurement of the B^0 -> D^*- pi^+ pi^- pi^+ branching fraction

abstract: Using a sample of (470.9 +- 2.8) x 10^6 BB-bar pairs, we measure the decay branching fraction B(B^0 -> D^*- pi^+ pi^- pi^+) = (7.26 +- 0.11 +- 0.31) x 10^-3, where the first uncertainty is statistical and the second is systematic. Our measurement will be helpful in studies of lepton universality by measuring B(B^0 -> D^*- tau^+ nu_tau) using tau^+ -> pi^+ pi^- pi^+ nu-bar_tau decays, normalized to B(B^0 -> D^*- pi^+ pi^- pi^+.


title: Fully pseudospectral solution of the conformally invariant wave equation near the cylinder at spacelike infinity. II: Schwarzschild background

abstract: It has recently been demonstrated (Class. Quantum Grav. 31, 085010, 2014) that the conformally invariant wave equation on a Minkowski background can be solved with a fully pseudospectral numerical method. In particular, it is possible to include spacelike infinity into the numerical domain, which is appropriately represented as a cylinder, and highly accurate numerical solutions can be obtained with a moderate number of gridpoints. In this paper, we generalise these considerations to the spherically-symmetric wave equation on a Schwarzschild background. In the Minkowski case, a logarithmic singularity at the future boundary is present at leading order, which can easily be removed to obtain completely regular solutions. An important new feature of the Schwarzschild background is that the corresponding solutions develop logarithmic singularities at infinitely many orders. This behaviour seems to be characteristic for massive space-times. In this sense this work is indicative of properties of the solutions of the Einstein equations near spatial infinity. The use of fully pseudospectral methods allows us to still obtain very accurate numerical solutions, and the convergence properties of the spectral approximations reveal details about the singular nature of the solutions on spacelike and null infinity. These results seem to be impossible to achieve with other current numerical methods. Moreover, we describe how to impose conditions on the asymptotic behaviour of initial data so that the leading-order logarithmic terms are avoided, which further improves the numerical accuracy.


title: Decoupled Asynchronous Proximal Stochastic Gradient Descent with Variance Reduction

abstract: In the era of big data, optimizing large scale machine learning problems becomes a challenging task and draws significant attention. Asynchronous optimization algorithms come out as a promising solution. Recently, decoupled asynchronous proximal stochastic gradient descent (DAP-SGD) is proposed to minimize a composite function. It is claimed to be able to off-loads the computation bottleneck from server to workers by allowing workers to evaluate the proximal operators, therefore, server just need to do element-wise operations. However, it still suffers from slow convergence rate because of the variance of stochastic gradient is nonzero. In this paper, we propose a faster method, decoupled asynchronous proximal stochastic variance reduced gradient descent method (DAP-SVRG). We prove that our method has linear convergence for strongly convex problem. Large-scale experiments are also conducted in this paper, and results demonstrate our theoretical analysis.


title: Most recent changepoint detection in Panel data

abstract: Detecting recent changepoints in time-series can be important for short-term prediction, as we can then base predictions just on the data since the changepoint. In many applications we have panel data, consisting of many related univariate time-series. We present a novel approach to detect sets of most recent changepoints in such panel data which aims to pool information across time-series, so that we preferentially infer a most recent change at the same time-point in multiple series. Our approach is computationally efficient as it involves analysing each time-series independently to obtain a profile likelihood like quantity that summarises the evidence for the series having either no change or a specific value for its most recent changepoint. We then post-process this output from each time-series to obtain a potentially small set of times for the most recent changepoints, and, for each time, the set of series which has their most recent changepoint at that time. We demonstrate the usefulness of this method on two data sets: forecasting events in a telecommunications network and inference about changes in the net asset ratio for a panel of US firms.


title: Asymptotic properties of adaptive group Lasso for sparse reduced rank regression

abstract: This paper studies the asymptotic properties of the penalized least squares estimator using an adaptive group Lasso penalty for the reduced rank regression. The group Lasso penalty is defined in the way that the regression coefficients corresponding to each predictor are treated as one group. It is shown that under certain regularity conditions, the estimator can achieve the minimax optimal rate of convergence. Moreover, the variable selection consistency can also be achieved, that is, the relevant predictors can be identified with probability approaching one. In the asymptotic theory, the number of response variables, the number of predictors, and the rank number are allowed to grow to infinity with the sample size.


title: Correctness Guarantees for the Composition of Lane Keeping and Adaptive Cruise Control

abstract: This paper develops a control approach with correctness guarantees for the simultaneous operation of lane keeping and adaptive cruise control. The safety specifications for these driver assistance modules are expressed in terms of set invariance. Control barrier functions are used to design a family of control solutions that guarantee the forward invariance of a set, which implies satisfaction of the safety specifications. The control barrier functions are synthesized through a combination of sum-of-squares program and physics-based modeling and optimization. A real-time quadratic program is posed to combine the control barrier functions with the performance-based controllers, which can be either expressed as control Lyapunov function conditions or as black-box legacy controllers. In both cases, the resulting feedback control guarantees the safety of the composed driver assistance modules in a formally correct manner. Importantly, the quadratic program admits a closed-form solution that can be easily implemented. The effectiveness of the control approach is demonstrated by simulations in the industry-standard vehicle simulator Carsim.


title: Neumann problem for p-Laplace equation in metric spaces using a variational approach: existence, boundedness, and boundary regularity

abstract: We employ a variational approach to study the Neumann boundary value problem for the $p$-Laplacian on bounded smooth-enough domains in the metric setting, and show that solutions exist and are bounded. The boundary data considered are Borel measurable bounded functions. We also study boundary continuity properties of the solutions. One of the key tools utilized is the trace theorem for Newton-Sobolev functions, and another is an analog of the De Giorgi inequality adapted to the Neumann problem.


title: An infinite family of vertex-primitive 2-arc-transitive digraphs

abstract: We solve the long-standing existence problem of vertex-primitive 2-arc-transitive digraphs by constructing an infinite family of such digraphs.


abstract: In this paper we confirm several conjectures of Z.-W. Sun on Hankel-type determinants for some combinatorial sequences including Franel numbers, Domb numbers and Ap\'ery numbers. For any nonnegative integer $n$, define \begin{gather*}f_n:=\sum_{k=0}^n\binom nk^3,\ D_n:=\sum_{k=0}^n\binom nk^2\binom{2k}k\binom{2(n-k)}{n-k}, b_n:=\sum_{k=0}^n\binom nk^2\binom{n+k}k,\ A_n:=\sum_{k=0}^n\binom nk^2\binom{n+k}k^2. \end{gather*} For $n=0,1,2,\ldots$, we show that $6^{-n}|f_{i+j}|_{0\leq i,j\leq n}$ and $12^{-n}|D_{i+j}|_{0\le i,j\le n}$ are positive odd integers, and $10^{-n}|b_{i+j}|_{0\leq i,j\leq n}$ and $24^{-n}|A_{i+j}|_{0\leq i,j\leq n}$ are always integers.


title: A lower bound of the distortion of the Torelli group in the mapping class group with boundary components

abstract: We prove that each Torelli group of an orientable surface with any number of boundary components is at least exponentially distorted in the mapping class group by using Broaddus-Farb-Putman's techniques. Further we show that the distortion of each Torelli group in the level $d$ mapping class group is the same as that of in the mapping class group.


title: Bottom crossing probability for symmetric jump processes (full version)

abstract: We determine the decay rate of the bottom crossing probability for symmetric jump processes under the condition on heat kernel estimates. Our results are applicable to symmetric stable-like processes and stable-subordinated diffusion processes on a class of (unbounded) fractals and fractal-like spaces.


title: Complexity of virtual 3-manifolds

abstract: Virtual $3$-manifolds were introduced by S.V. Matveev in 2009 as natural generalizations of the classical $3$-manifolds. In this paper, we introduce a notion of complexity of a virtual $3$-manifold. We investigate the values of the complexity for virtual 3-manifolds presented by special polyhedra with one or two $2$-components. On the basis of these results, we establish the exact values of the complexity for a wide class of hyperbolic $3$-manifolds with totally geodesic boundary.


title: Escape rate of the Brownian motions on hyperbolic spaces

abstract: We discuss the escape rate of the Brownian motion on a hyperbolic space. We point out that the escape rate is determined by using the Brownian expression of the radial part and a generalized Kolmogorov's test for the one dimensional Brownian motion.


title: How correlated is the FeSe/SrTiO$_3$ system ?

abstract: Recent observation of $\sim$ 10 times higher critical temperature in FeSe monolayer compared with its bulk phase has drawn a great deal of attention because the electronic structure in the monolayer phase appears to be different than bulk FeSe. Using a combination of density functional theory and dynamical mean field theory, we find electronic correlations have important effects on the predicted atomic-scale geometry and the electronic structure of the monolayer FeSe on SrTiO$_3$. The electronic correlations are dominantly controlled by the Se-Fe-Se angle either in the bulk phase or the monolayer phase. But the angle sensitivity increases and the orbital differentiation decreases in the monolayer phase compared to the bulk phase. The correlations are more dependent on Hund's J than Hubbard U. The observed orbital selective incoherence to coherence crossover with temperature confirms the Hund's metallic nature of the monolayer FeSe. We also find electron doping by oxygen vacancies in SrTiO$_3$ increases the correlation strength, especially in the $d_{xy}$ orbital by reducing the Se-Fe-Se angle.


title: Ultrafast generation of skyrmionic defects with vortex beams: printing laser profiles on magnets

abstract: Controlling electric and magnetic properties of matter by laser beams is actively explored in the broad region of condensed matter physics, including spintronics and magneto-optics. Here we theoretically propose an application of optical and electron vortex beams carrying intrinsic orbital angular momentum to chiral ferro- and antiferro- magnets. We analyze the time evolution of spins in chiral magnets under irradiation of vortex beams, by using the stochastic Landau-Lifshitz-Gilbert equation. We show that beam-driven nonuniform temperature lead to a class of ring-shaped magnetic defects, what we call skyrmion multiplex, as well as conventional skyrmions. We discuss the proper beam parameters and the optimal way of applying the beams for the creation of these topological defects. Our findings provide an ultrafast scheme of generating topological magnetic defects in a way applicable to both metallic and insulating chiral (anti-) ferromagnets.


title: Fast Evolution and Waveform Generator for Extreme-Mass-Ratio Inspirals in Equatorial-Circular Orbits

abstract: In this paper we discuss the development of a fast and accurate waveform model for the quasi-circular orbital evolution of extreme-mass-ratio-inspirals (EMRIs). This model simply employs the data of a few numerical Teukoulsky-based energy fluxes and waveforms to fit out a set of polynomials for the entire fluxes and waveforms. These obtained polynomials are accurate enough in the entire evolution domain, and much more accurate than the resummation post-Newtonian (PN) energy fluxes and waveforms, especially when the spin of a black hole becomes large. The dynamical equation we adopted for orbital revolution is the effective-one-body (EOB) formalism. Because of the simplified expressions, the efficiency of calculating the orbital evolution with our polynomials is also better than the traditional method which uses the resummed PN analytical fluxes. Our model should be useful in calculation of waveform templates of EMRIs for the gravitational wave detectors such as the evolved Laser Interferometer Space Antenna (eLISA).


title: Computing the monodromy and pole order filtration on Milnor fiber cohomology of plane curves

abstract: We describe an algorithm computing the monodromy and the pole order filtration on the Milnor fiber cohomology of any reduced projective plane curve $C$. The relation to the zero set of Bernstein-Sato polynomial of the defining homogeneous polynomial for $C$ is also discussed. When $C$ has some non weighted homogeneous singularities, then we have to assume that a conjecture holds in order to get some of our results. In all the examples computed so far this conjecture holds.


title: Preliminary Test and Shrinkage Estimations of Scale Parameters for Two Exponential Distributions based on Record Values

abstract: The exponential distribution is applied in a very wide variety of statistical procedures. Among the most prominent applications are those in the field of life testing and reliability theory. When there are two record samples available for estimating the scale parameter, a preliminary test is usually used to determine whether to pool the samples or use the individual sample. In this paper, the preliminary test estimator and shrinkage estimator are studied. The optimum level of significance for preliminary test estimation and the optimum values of shrinkage coefficient are obtained based on minimax regret criterion under the weighted square error loss function.


title: Entanglement dynamics for uniformly accelerated two-level atoms coupled with electromagnetic vacuum fluctuations

abstract: We investigate the entanglement dynamics of two uniformly accelerated atoms with the same acceleration perpendicular to their separation. The two-atom system is treated as an open system coupled with fluctuating electromagnetic fields in the Minkowski vacuum, and in the Born-Markov approximation the master equation that describes the completely positive time evolution of the two-atom system is derived. In particular, we investigate the phenomena of entanglement degradation, generation, revival and enhancement. As opposed to the scalar-field case, the entanglement dynamics is crucially dependent on the polarization directions of the atoms. For the two-atom system with certain acceleration and separation, the polarization directions of the atoms may determine whether entanglement generation, revival or enhancement happens, while for entanglement degradation, they affect the decay rate of entanglement. A comparison between the entanglement evolution of accelerated atoms and that of static ones immersed in a thermal bath at the Unruh temperature shows that they are the same only when the acceleration is extremely small.


title: An Exponential Inequality for U-Statistics under Mixing Conditions

abstract: The family of U-statistics plays a fundamental role in statistics. This paper proves a novel exponential inequality for U-statistics under the time series setting. Explicit mixing conditions are given for guaranteeing fast convergence, the bound proves to be analogous to the one under independence, and extension to non-stationary time series is straightforward. The proof relies on a novel decomposition of U-statistics via exploiting the temporal correlatedness structure. Such results are of interest in many fields where high dimensional time series data are present. In particular, applications to high dimensional time series inference are discussed.


title: Effect of adding nanometre-sized heterogeneities on the structural dynamics and the excess wing of a molecular glass former

abstract: We present the relaxation dynamics of glass-forming glycerol mixed with 1.1 nm sized polyhedral oligomeric silsesquioxane (POSS) molecules using dielectric spectroscopy (DS) and two different neutron scattering (NS) techniques. Both, the reorientational dynamics as measured by DS and the density fluctuations detected by NS reveal a broadening of the alpha relaxation when POSS molecules are added. Moreover, we find a significant slowing down of the alpha-relaxation time. These effects are in accord with the heterogeneity scenario considered for the dynamics of glasses and supercooled liquids. The addition of POSS also affects the excess wing in glycerol arising from a secondary relaxation process, which seems to exhibit a dramatic increase in relative strength compared to the alpha-relaxation.


title: Strategic Seeding of Rival Opinions

abstract: We present a network influence game that models players strategically seeding the opinions of nodes embedded in a social network. A social learning dynamic, whereby nodes repeatedly update their opinions to resemble those of their neighbors, spreads the seeded opinions through the network. After a fixed period of time, the dynamic halts and each player's utility is determined by the relative strength of the opinions held by each node in the network vis-a-vis the other players. We show that the existence of a pure Nash equilibrium cannot be guaranteed in general. However, if the dynamics are allowed to progress for a sufficient amount of time so that a consensus among all of the nodes is obtained, then the existence of a pure Nash equilibrium can be guaranteed. The computational complexity of finding a pure strategy best response is shown to be NP-complete, but can be efficiently approximated to within a (1 - 1/e) factor of optimal by a simple greedy algorithm.


title: Convex PBW-type Lyndon Bases and Restricted Two-Parameter Quantum Group of Type $F_4$

abstract: We determine convex PBW-type Lyndon bases for two-parameter quantum groups $U_{r,s}(F_4)$ with detailed commutation relations. We construct a finite-dimensional Hopf algebra $\mathfrak u_{r,s}(F_4)$, as a quotient of $U_{r,s}(F_4)$ by a Hopf ideal generated by certain central elements, which is pointed, and of a Drinfel'd double structure under a certain condition. All of Hopf isomorphisms of $\mathfrak u_{r,s}(F_4)$ are determined which are important for seeking the possible new pointed objects in low order with $(\ell, 210)\ne 1$. Finally, necessary and sufficient conditions for $\mathfrak u_{r,s}(F_4)$ to be a ribbon Hopf algebra are singled out by describing the left and right integrals.


title: Algorithmic Persuasion with No Externalities

abstract: We study the algorithmics of information structure design -- a.k.a. persuasion or signaling -- in a fundamental special case introduced by Arieli and Babichenko: multiple agents, binary actions, and no inter-agent externalities. Unlike prior work on this model, we allow many states of nature. We assume that the principal's objective is a monotone set function, and study the problem both in the public signal and private signal models, drawing a sharp contrast between the two in terms of both efficacy and computational complexity. When private signals are allowed, our results are largely positive and quite general. First, we show polynomial-time equivalence between optimal signaling and the problem of maximizing the objective function plus an additive function. This yields an efficient implementation of the optimal scheme when the objective is supermodular or anonymous. Second, we exhibit a (1-1/e)-approximation of the optimal private signaling scheme, modulo an additive loss of $\epsilon$, when the objective function is submodular. These two results simplify, unify, and generalize results of [Arieli and Babichenko, 2016] and [Babichenko and Barman, 2016], extending them from a binary state of nature to many states (modulo the additive loss in the latter result). Third, we consider the binary-state case with a submodular objective, and simplify and slightly strengthen the result of [Babichenko and Barman, 2016] to obtain a (1-1/e)-approximation via an explicitly constructed scheme. When only a public signal is allowed, our results are negative. First, we show that it is NP-hard to approximate the optimal public scheme, within any constant factor, even when the objective is additive. Second, we show that the optimal private scheme can outperform the optimal public scheme, in terms of maximizing the sender's objective, by a polynomial factor.


title: Bibliographic Analysis with the Citation Network Topic Model

abstract: Bibliographic analysis considers author's research areas, the citation network and paper content among other things. In this paper, we combine these three in a topic model that produces a bibliographic model of authors, topics and documents using a non-parametric extension of a combination of the Poisson mixed-topic link model and the author-topic model. We propose a novel and efficient inference algorithm for the model to explore subsets of research publications from CiteSeerX. Our model demonstrates improved performance in both model fitting and a clustering task compared to several baselines.


title: On embeddings of the Grassmannian $Gr(2,m)$ into the Grassmannian $Gr(2,n)$

abstract: In this {paper}, we consider holomorphic embeddings of $Gr(2,m)$ into $Gr(2,n)$. We study such embeddings by finding all possible total Chern classes of the pull-back of the universal bundles under these embeddings. Using the relations between Chern classes of the universal bundles and Schubert cycles, together with properties of holomorphic vector bundles of rank $2$ on complex Grassmannians, we find conditions on $m$ and $n$ for which any holomorphic embedding of $Gr(2,m)$ into $Gr(2,n)$ is linear.


title: Six-wave mixing induced by free-carrier plasma in silicon nanowire waveguides

abstract: Nonlinear wave mixing in mesoscopic silicon structures is a fundamental nonlinear process with broad impact and applications. Silicon nanowire waveguides, in particular, have large third-order Kerr nonlinearity, enabling salient and abundant four-wave-mixing dynamics and functionalities. Besides the Kerr effect, in silicon waveguides two-photon absorption generates high free-carrier densities, with corresponding fifth-order nonlinearity in the forms of free-carrier dispersion and free-carrier absorption. However, whether these fifth-order free-carrier nonlinear effects can lead to six-wave-mixing dynamics still remains an open question until now. Here we report the demonstration of free-carrier-induced six-wave mixing in silicon nanowires. Unique features, including inverse detuning dependence of six-wave-mixing efficiency and its higher sensitivity to pump power, are originally observed and verfied by analytical prediction and numerical modeling. Additionally, asymmetric sideband generation is observed for different laser detunings, resulting from the phase-sensitive interactions between free-carrier six-wave-mixing and Kerr four-wave-mixing dynamics. These discoveries provide a new path for nonlinear multi-wave interactions in nanoscale platforms.


title: Summation identities and transformations for hypergeometric series

abstract: We find summation identities and transformations for the McCarthy's $p$-adic hypergeometric series by evaluating certain Gauss sums which appear while counting points on the family $$Z_{\lambda}: x_1^d+x_2^d=d\lambda x_1x_2^{d-1}$$ over a finite field $\mathbb{F}_p$. A. Salerno expresses the number of points over a finite field $\mathbb{F}_p$ on the family $Z_{\lambda}$ in terms of quotients of $p$-adic gamma function under the condition that $d|p-1$. In this paper, we first express the number of points over a finite field $\mathbb{F}_p$ on the family $Z_{\lambda}$ in terms of McCarthy's $p$-adic hypergeometric series for any odd prime $p$ not dividing $d(d-1)$, and then deduce two summation identities for the $p$-adic hypergeometric series. We also find certain transformations and special values of the $p$-adic hypergeometric series. We finally find a summation identity for the Greene's finite field hypergeometric series.


title: Nonparametric Density Estimation for Spatial Data with Wavelets

abstract: Nonparametric density estimators are studied for $d$-dimensional, strongly spatial mixing data which is defined on a general $N$-dimensional lattice structure. We consider linear and nonlinear hard thresholded wavelet estimators which are derived from a $d$-dimensional multiresolution analysis. We give sufficient criteria for the consistency of these estimators and derive rates of convergence in $L^{p'}$ for $p'\in [1,\infty)$. For this reason, we study density functions which are elements of a $d$-dimensional Besov space $B^s_{p,q}(\mathbb{R}^d)$. We also verify the analytic correctness of our results in numerical simulations.


title: Hawkes Processes with Stochastic Excitations

abstract: We propose an extension to Hawkes processes by treating the levels of self-excitation as a stochastic differential equation. Our new point process allows better approximation in application domains where events and intensities accelerate each other with correlated levels of contagion. We generalize a recent algorithm for simulating draws from Hawkes processes whose levels of excitation are stochastic processes, and propose a hybrid Markov chain Monte Carlo approach for model fitting. Our sampling procedure scales linearly with the number of required events and does not require stationarity of the point process. A modular inference procedure consisting of a combination between Gibbs and Metropolis Hastings steps is put forward. We recover expectation maximization as a special case. Our general approach is illustrated for contagion following geometric Brownian motion and exponential Langevin dynamics.


title: Pre-adjunctions and the Ramsey property

abstract: Showing that the Ramsey property holds for a class of finite structures can be an extremely challenging task and a slew of sophisticated methods have been proposed in literature. In this paper we propose a new strategy to show that a class of structures has the Ramsey property. The strategy is based on a relatively simple result in category theory and consists of establishing a pre-adjunction between the category of structures which is known to have the Ramsey property, and the category of structures we are interested in. This strategy was implicitly used already in 1981 by H.J. Pr\"omel and B. Voigt in their proof of the Ramsey property for the class of finite linearly ordered graphs. We demonstrate the applicability of this strategy by providing short proofs of three important well known results: we show the Ramsey property for the category of all finite linearly ordered posets with embeddings, for the category of finite convexly ordered ultrametric spaces with embeddings, and for the category of all finite linearly ordered metric spaces (rational metric spaces) with embeddings.


title: Collective Coordinates in One-Dimensional Soliton Models Revisited

abstract: We compare numerical solutions to the full field equations to simplified approaches based on implementing three collective coordinates for kink-antikink interactions within the $\varphi^4$ and $\phi^6$ models in one time and one space dimensions. We particularly pursue the question whether the collective coordinate approximation substantiates the conjecture that vibrational modes are important for resonance structures to occur in kink-antikink scattering.


title: Revisiting gravitino dark matter in thermal leptogenesis

abstract: In this paper, we revisit the gravitino dark matter scenario in the presence of the bilinear $R$-parity violating interaction. In particular, we discuss a consistency with the thermal leptogenesis. For a high reheating temperature required for the thermal leptogenesis, the gravitino dark matter tends to be overproduced, which puts a severe upper limit on the gluino mass. As we will show, a large portion of parameter space of the gravitino dark matter scenario has been excluded by combining the constraints from the gravitino abundance and the null results of the searches for the superparticles at the LHC experiments. In particular, the models with the stau (and other charged slepton) NLSP has been almost excluded by the searches for the long-lived charged particles at the LHC unless the required reheating temperature is somewhat lowered by assuming, for example, a degenerated right-handed neutrino mass spectrum.


title: Experimental time-optimal universal control of spin qubits in solids

abstract: Quantum control of systems plays important roles in modern science and technology. The ultimate goal of quantum control is to achieve high fidelity universal control in the time-optimal way. Although high fidelity universal control has been reported in various quantum systems, experimental implementation of time-optimal universal control remains elusive. Here we report the experimental realization of time-optimal universal control of spin qubits in diamond. By generalizing a recent method for solving quantum brachistochrone equations [X. Wang, et al., Phys. Rev. Lett. 114, 170501 (2015)], we obtained accurate minimum time protocols for multiple qubits with fixed qubits' interactions and constrained control field. Single- and two-qubit time-optimal gates are experimentally implemented with fidelities of 99% obtained via quantum process tomography. Our work provides a time-optimal route to achieve accurate quantum control, and unlocks new capabilities for emerging field of time-optimal control in general quantum systems.


title: Stripe and short range order in the charge density wave of 1T-Cu$_x$TiSe$_2$

abstract: We study the impact of Cu intercalation on the charge density wave (CDW) in 1T-Cu$_{\text{x}}$TiSe$_{\text{2}}$ by scanning tunneling microscopy and spectroscopy. Cu atoms, identified through density functional theory modeling, are found to intercalate randomly on the octahedral site in the van der Waals gap and to dope delocalized electrons near the Fermi level. While the CDW modulation period does not depend on Cu content, we observe the formation of charge stripe domains at low Cu content (x$<$0.02) and a breaking up of the commensurate order into 2$\times$2 domains at higher Cu content. The latter shrink with increasing Cu concentration and tend to be phase-shifted. These findings invalidate a proposed excitonic pairing as the primary CDW formation mechanism in this material.


title: High-temperature expansion for frustrated magnets: Application to the J1-J2 model on the BCC lattice

abstract: We present the high-temperature expansion up to 11th order for the specific heat $C$ and the uniform susceptibility $\chi_0$ and up to 9th order for the structure factor $S_{\bf Q}$ of the frustrated spin-half $J_1$-$J_2$ Heisenberg model on the BCC lattice. We consider ferromagnetic as well as antiferromagnetic nearest-neighbor exchange $J_1$ and frustrating antiferromagnetic next-nearest-neighbor exchange $J_2$. We discuss the influence of frustration on the temperature dependence of these quantities. Furthermore, we use the HTE series to determine the critical temperature $T_c$ as a function of the frustration parameter $J_2$.


title: Deep-Learned Collision Avoidance Policy for Distributed Multi-Agent Navigation

abstract: High-speed, low-latency obstacle avoidance that is insensitive to sensor noise is essential for enabling multiple decentralized robots to function reliably in cluttered and dynamic environments. While other distributed multi-agent collision avoidance systems exist, these systems require online geometric optimization where tedious parameter tuning and perfect sensing are necessary. We present a novel end-to-end framework to generate reactive collision avoidance policy for efficient distributed multi-agent navigation. Our method formulates an agent's navigation strategy as a deep neural network mapping from the observed noisy sensor measurements to the agent's steering commands in terms of movement velocity. We train the network on a large number of frames of collision avoidance data collected by repeatedly running a multi-agent simulator with different parameter settings. We validate the learned deep neural network policy in a set of simulated and real scenarios with noisy measurements and demonstrate that our method is able to generate a robust navigation strategy that is insensitive to imperfect sensing and works reliably in all situations. We also show that our method can be well generalized to scenarios that do not appear in our training data, including scenes with static obstacles and agents with different sizes. Videos are available at https://sites.google.com/view/deepmaca.


title: A Spectral Projection Preconditioner for Solving Ill Conditioned Linear Systems

abstract: We present a preconditioner based on spectral projection that is combined with a deflated Krylov subspace method for solving ill conditioned linear systems of equations. Our results show that the proposed algorithm requires many fewer iterations to achieve the convergence criterion for solving an ill conditioned problem than a Krylov subspace solver. In our numerical experiments, the solution obtained by the proposed algorithm is more accurate in terms of the norm of the distance to the exact solution of the linear system of equations.


title: Exact Sampling from Determinantal Point Processes

abstract: Determinantal point processes (DPPs) are an important concept in random matrix theory and combinatorics. They have also recently attracted interest in the study of numerical methods for machine learning, as they offer an elegant "missing link" between independent Monte Carlo sampling and deterministic evaluation on regular grids, applicable to a general set of spaces. This is helpful whenever an algorithm explores to reduce uncertainty, such as in active learning, Bayesian optimization, reinforcement learning, and marginalization in graphical models. To draw samples from a DPP in practice, existing literature focuses on approximate schemes of low cost, or comparably inefficient exact algorithms like rejection sampling. We point out that, for many settings of relevance to machine learning, it is also possible to draw exact samples from DPPs on continuous domains. We start from an intuitive example on the real line, which is then generalized to multivariate real vector spaces. We also compare to previously studied approximations, showing that exact sampling, despite higher cost, can be preferable where precision is needed.


title: PSVR - Self-stabilizing Publish/Subscribe Communication for Ad-hoc Networks

abstract: This paper presents the novel routing algorithm PSVR for pub/sub systems in ad-hoc networks. Its focus is on scenarios where communications links are unstable and nodes frequently change subscriptions. PSVR presents a compromise of size and maintenance effort for routing tables due to sub- and unsubscriptions and the length of routing paths. Designed in a self-stabilizing manner it scales well with network size. The evaluation reveals that PSVR only needs slightly more messages than a close to optimal routing structure for publication delivery, and creates shorter routing paths than an existing self-stabilizing algorithm. A real world deployment shows the usability of the approach


title: Efficient big data assimilation through sparse representation: A 3D benchmark case study in seismic history matching

abstract: In a previous work \citep{luo2016sparse2d_spej}, the authors proposed an ensemble-based 4D seismic history matching (SHM) framework, which has some relatively new ingredients, in terms of the type of seismic data in choice, the way to handle big seismic data and related data noise estimation, and the use of a recently developed iterative ensemble history matching algorithm. In seismic history matching, it is customary to use inverted seismic attributes, such as acoustic impedance, as the observed data. In doing so, extra uncertainties may arise during the inversion processes. The proposed SHM framework avoids such intermediate inversion processes by adopting amplitude versus angle (AVA) data. In addition, SHM typically involves assimilating a large amount of observed seismic attributes into reservoir models. To handle the big-data problem in SHM, the proposed framework adopts the following wavelet-based sparse representation procedure: First, a discrete wavelet transform is applied to observed seismic attributes. Then, uncertainty analysis is conducted in the wavelet domain to estimate noise in the resulting wavelet coefficients, and to calculate a corresponding threshold value. Wavelet coefficients above the threshold value, called leading wavelet coefficients hereafter, are used as the data for history matching. The retained leading wavelet coefficients preserve the most salient features of the observed seismic attributes, whereas rendering a substantially smaller data size. Finally, an iterative ensemble smoother is adopted to update reservoir models, in such a way that the leading wavelet coefficients of simulated seismic attributes better match those of observed seismic attributes. (The rest of the abstract was omitted for the length restriction.)

categories: physics.data-an math.NA math.OC physics.geo-ph stat.CO


title: Fluctuations in protein aggregation: Design of preclinical screening for early diagnosis of neurodegenerative disease

abstract: Autocatalytic fibril nucleation has recently been proposed to be a determining factor for the spread of neurodegenerative diseases, but the same process could also be exploited to amplify minute quantities of protein aggregates in a diagnostic context. Recent advances in microfluidic technology allow analysis of protein aggregation in micron-scale samples potentially enabling such diagnostic approaches, but the theoretical foundations for the analysis and interpretation of such data are so far lacking. Here we study computationally the onset of protein aggregation in small volumes and show that the process is ruled by intrinsic fluctuations whose volume dependent distribution we also estimate theoretically. Based on these results, we develop a strategy to quantify in silico the statistical errors associated with the detection of aggregate containing samples. Our work opens a new perspective on the forecasting of protein aggregation in asymptomatic subjects.


title: Posted Pricing sans Discrimination

abstract: In the quest for market mechanisms that are easy to implement, yet close to optimal, few seem as viable as posted pricing. Despite the growing body of impressive results, the performance of most posted price mechanisms however, rely crucially on price discrimination when multiple copies of a good are available. For the more general case with non-linear production costs on each good, hardly anything is known for general multi-good markets. With this in mind, we study a Bayesian setting where the seller can produce any number of copies of a good but faces convex production costs for the same, and buyers arrive sequentially. Our main contribution is a framework for non-discriminatory pricing in the presence of production costs: the framework yields posted price mechanisms with O(1)-approximation factors for fractionally subadditive (XoS) buyers, logarithmic approximations for subadditive buyers, and also extends to settings where the seller is oblivious to buyer valuations. Our work presents the first known results for Bayesian settings with production costs and is among the few posted price mechanisms that do not charge buyers differently for the same good.


title: On the usability of deep networks for object-based image analysis

abstract: As computer vision before, remote sensing has been radically changed by the introduction of Convolution Neural Networks. Land cover use, object detection and scene understanding in aerial images rely more and more on deep learning to achieve new state-of-the-art results. Recent architectures such as Fully Convolutional Networks (Long et al., 2015) can even produce pixel level annotations for semantic mapping. In this work, we show how to use such deep networks to detect, segment and classify different varieties of wheeled vehicles in aerial images from the ISPRS Potsdam dataset. This allows us to tackle object detection and classification on a complex dataset made up of visually similar classes, and to demonstrate the relevance of such a subclass modeling approach. Especially, we want to show that deep learning is also suitable for object-oriented analysis of Earth Observation data. First, we train a FCN variant on the ISPRS Potsdam dataset and show how the learnt semantic maps can be used to extract precise segmentation of vehicles, which allow us studying the repartition of vehicles in the city. Second, we train a CNN to perform vehicle classification on the VEDAI (Razakarivony and Jurie, 2016) dataset, and transfer its knowledge to classify candidate segmented vehicles on the Potsdam dataset.


title: Semantic Segmentation of Earth Observation Data Using Multimodal and Multi-scale Deep Networks

abstract: This work investigates the use of deep fully convolutional neural networks (DFCNN) for pixel-wise scene labeling of Earth Observation images. Especially, we train a variant of the SegNet architecture on remote sensing data over an urban area and study different strategies for performing accurate semantic segmentation. Our contributions are the following: 1) we transfer efficiently a DFCNN from generic everyday images to remote sensing images; 2) we introduce a multi-kernel convolutional layer for fast aggregation of predictions at multiple scales; 3) we perform data fusion from heterogeneous sensors (optical and laser) using residual correction. Our framework improves state-of-the-art accuracy on the ISPRS Vaihingen 2D Semantic Labeling dataset.


title: Deformation and fracture of echinoderm collagen networks

abstract: Collagen networks provide the main structural component of most tissues and represent an important ingredient for bio-mimetic materials for bio-medical applications. Here we study the mechanical properties of stiff collagen networks derived from three different echinoderms and show that they exhibit non-linear stiffening followed by brittle fracture. The disordered nature of the network leads to strong sample-to-sample fluctuations in elasticity and fracture strength. We perform numerical simulations of a three dimensional model for the deformation of a cross-linked elastic fibril network which is able to reproduce the macroscopic features of the experimental results and provide insights into the internal mechanics of stiff collagen networks. Our numerical model provides an avenue for the design of collagen membranes with tunable mechanical properties.


title: Production-Driven Patch Generation and Validation

abstract: We envision a world where the developer would receive each morning in her GitHub dashboard a list of potential patches that fix certain production failures. For this, we propose a novel program repair scheme, with the unique feature of being applicable to production directly. We present the design and implementation of a prototype system for Java, called Itzal, that performs patch generation for uncaught exceptions in production. We have performed two empirical experiments to validate our system: the first one on 34 failures from 14 different software applications, the second one on 16 seeded failures in 3 real open-source e-commerce applications for which we have set up a realistic user traffic. This validates the novel and disruptive idea of using program repair directly in production.


title: Existence of solutions for a class of fourth order cross-diffusion systems of gradient flow type

abstract: This article is concerned with the existence and the long time behavior of weak solutions to certain coupled systems of fourth-order degenerate parabolic equations of gradient flow type. The underlying metric is a Wasserstein-like transportation distance for vector-valued functions, with nonlinear mobilities in each component. Under the hypothesis of (flat) convexity of the driving free energy functional, weak solutions are constructed by means of the variational minimizing movement scheme for metric gradient flows. The essential regularity estimates are derived by variational methods.


title: Resistivity Minimum in Granular Composites and Thin Metallic Films

abstract: We analyze the temperature dependence of conductivity in thick granular ferromagnetic compounds NiSiO2 and in thin weakly coupled films of Fe, Ni and Py in vicinity of metal-insulator transition. Development of resistivity minimum followed by a logarithmic variation of conductivity at lower temperatures is attributed to granular structure of compounds and thin films fabricated by conventional deposition techniques. Resistivity minimum is identified as a transition between temperature dependent intra-granular metallic conductance and thermally activated inter-granular tunneling.


title: Long cycles in locally expanding graphs, with applications

abstract: We provide sufficient conditions for the existence of long cycles in locally expanding graphs, and present applications of our conditions and techniques to Ramsey theory, random graphs and positional games.


title: Criticality in Two-Dimensional Quantum Systems: Tensor Network Approach

abstract: Determination and characterization of criticality in two-dimensional (2D) quantum many-body systems belong to the most important challenges and problems of quantum physics. In this paper we propose an efficient scheme to solve this problem by utilizing the infinite projected entangled pair state (iPEPS), and tensor network (TN) representations. We show that the criticality of a 2D state is faithfully reproduced by the ground state (dubbed as boundary state) of a one-dimensional effective Hamiltonian constructed from its iPEPS representation. We demonstrate that for a critical state the correlation length and the entanglement spectrum of the boundary state are essentially different from those of a gapped iPEPS. This provides a solid indicator that allows to identify the criticality of the 2D state. Our scheme is verified on the resonating valence bond (RVB) states on kagom\'e and square lattices, where the boundary state of the honeycomb RVB is found to be described by a $c=1$ conformal field theory. We apply our scheme also to the ground state of the spin-1/2 XXZ model on honeycomb lattice, illustrating the difficulties of standard variational TN approaches to study the critical ground states. Our scheme is of high versatility and flexibility, and can be further applied to investigate the quantum criticality in many other phenomena, such as finite-temperature and topological phase transitions.


title: Slowly synchronizing automata with fixed alphabet size

abstract: It was conjectured by \v{C}ern\'y in 1964 that a synchronizing DFA on $n$ states always has a shortest synchronizing word of length at most $(n-1)^2$, and he gave a sequence of DFAs for which this bound is reached. In this paper, we investigate the role of the alphabet size. For each possible alphabet size, we count DFAs on $n \le 6$ states which synchronize in $(n-1)^2 - e$ steps, for all $e < 2\lceil n/2 \rceil$. Furthermore, we give constructions of automata with any number of states, and $3$, $4$, or $5$ symbols, which synchronize slowly, namely in $n^2 - 3n + O(1)$ steps. In addition, our results prove \v{C}ern\'y's conjecture for $n \le 6$. Our computation has led to $27$ DFAs on $3$, $4$, $5$ or $6$ states, which synchronize in $(n-1)^2$ steps, but do not belong to \v{C}ern\'y's sequence. Of these $27$ DFA's, $19$ are new, and the remaining $8$ which were already known are exactly the \emph{minimal} ones: they will not synchronize any more after removing a symbol. So the $19$ new DFAs are extensions of automata which were already known, including the \v{C}ern\'y automaton on $3$ states. But for $n > 3$, we prove that the \v{C}ern\'y automaton on $n$ states does not admit non-trivial extensions with the same smallest synchronizing word length $(n-1)^2$.


title: On the joint distribution of first-passage time and first-passage area of drifted Brownian motion

abstract: For drifted Brownian motion $X(t)= x - \mu t + B_t \ (\mu >0)$ starting from $x>0,$ we study the joint distribution of the first-passage time below zero, $\tau(x),$ and the first-passage area, $A(x),$ swept out by $X$ till the time $\tau(x).$ In particular, we establish differential equations with boundary conditions for the joint moments $E[\tau(x)^m A(x)^n],$ and we present an algorithm to find recursively them, for any $m$ and $n.$ Finally, the expected value of the time average of $X$ till the time $\tau(x)$ is obtained.


title: On the applicability of density dependent effective interactions in cluster-forming systems

abstract: We systematically studied the validity and transferability of effective, coarse-grained, pair potentials in ultrasoft colloidal systems. We focused on amphiphilic dendrimers, macromolecules which can aggregate into clusters of overlapping particles to minimize the contact area with the surrounding (implicit) solvent. Simulations were performed for both the monomeric and coarse-grained model in the liquid phase at densities ranging from infinite dilution up to values close to the freezing point. For every state point, each macromolecule was mapped onto a single interaction site and the effective pair potential was computed using a coarse-graining technique based on force-matching. We found excellent agreement between the spatial dendrimer distributions obtained from the coarse-grained and microscopically detailed simulations at low densities, where the macromolecules were distributed homogeneously in the system. However, the agreement deteriorated significantly when the density was increased further and the cluster occupation became more polydisperse. Under these conditions, the effective pair potential of the coarse-grained model can no longer be computed by averaging over the whole system, but the local density needs to be taken into account instead.


title: Space-time large deviations in capacity-constrained relay networks

abstract: We consider a single-cell network of random transmitters and fixed relays in a bounded domain of Euclidean space. The transmitters arrive over time and select one relay according to a spatially inhomogeneous preference kernel. Once a transmitter is connected to a relay, the connection remains and the relay is occupied. If an occupied relay is selected by another transmitters with later arrival time, this transmitter becomes frustrated. We derive a large deviation principle for the space-time evolution of frustrated transmitters in the high-density regime.


title: Carbon and hydrogen radio recombination lines from the cold clouds towards Cassiopeia A

abstract: We use the Low Frequency Array to perform a systematic high spectral resolution investigation of the low-frequency 33-78 MHz spectrum along the line of sight to Cassiopeia A. We complement this with a 304-386 MHz Westerbork Synthesis Radio telescope observation. In this first paper we focus on the carbon radio recombination lines. We detect Cn$\alpha$ lines at -47 and -38 km s$^{-1}$ in absorption for quantum numbers n=438-584 and in emission for n=257-278 with high signal to noise. These lines are associated with cold clouds in the Perseus spiral arm component. Hn$\alpha$ lines are detected in emission for n=257-278. In addition, we also detect Cn$\alpha$ lines at 0 km s$^{-1}$ associated with the Orion arm. We analyze the optical depth of these transitions and their line width. Our models show that the carbon line components in the Perseus arm are best fit with an electron temperature 85 K and an electron density 0.04 cm$^{-3}$ and can be constrained to within 15\%. The electron pressure is constrained to within 20\%. We argue that much of these carbon radio recombination lines arise in the CO-dark surface layers of molecular clouds where most of the carbon is ionized but hydrogen has made the transition from atomic to molecular. The hydrogen lines are clearly associated with the carbon line emitting clouds, but the low-frequency upperlimits indicate that they likely do not trace the same gas. Combining the hydrogen and carbon results we arrive at a firm lower limit to the cosmic ray ionization rate of 2.5$\times$10$^{-18}$ s$^{-1}$, but the actual value is likely much larger.


title: Triplet $p$-wave pairing correlation in low doped zigzag graphene nanoribbons

abstract: We reveal an edge spin triplet $p-$wave superconducting pairing correlation in slightly doped zigzag graphene nanoribbons. By employing a method that combines random-phase approximation, the finite-temperature determinant quantum Monte Carlo approach, and the ground-state constrained-path quantum Monte Carlo method, it is shown that such a spin-triplet pairing is mediated by the ferromagnetic fluctuations caused by the flat band at the edge. The spin susceptibility and effective pairing interactions at the edge strongly increase as the on-site Coulomb interaction increases, indicating the importance of electron-electron correlations. It is also found that the doping-dependent ground-state $p$-wave pairing correlation bears some similarity to the famous superconducting dome in the phase diagram of a high-temperature superconductor, while the spin correlation at the edge is weakened as the system is doped away from half filling.


title: On the difference between gamma-ray-detected and non-gamma-ray-detected pulsars

abstract: We compare radio profile widths of young, energetic gamma-ray-detected and non-gamma-ray-detected pulsars. We find that the latter typically have wider radio profiles, with the boundary between the two samples exhibiting a dependence on the rate of rotational energy loss. We also find that within the sample of gamma-ray-detected pulsars, radio profile width is correlated with both the separation of the main gamma-ray peaks and the presence of narrow gamma-ray components. These findings lead us to propose that these pulsars form a single population where the main factors determining gamma ray detectability are the rate of rotational energy loss and the proximity of the line of sight to the rotation axis. The expected magnetic inclination angle distribution will be different for radio pulsars with and without detectable gamma rays, naturally leading to the observed differences. Our results also suggest that the geometry of existing radio and outer-magnetosphere gamma-ray emission models are at least qualitatively realistic, implying that information about the viewing geometry can be extracted from profile properties of pulsars.


title: Equivalence between minimal time and minimal norm control problems for the heat equation

abstract: This paper presents the equivalence between minimal time and minimal norm control problems for internally controlled heat equations. The target is an arbitrarily fixed bounded, closed and convex set with a nonempty interior in the state space. This study differs from [G. Wang and E. Zuazua, \textit{On the equivalence of minimal time and minimal norm controls for internally controlled heat equations}, SIAM J. Control Optim., 50 (2012), pp. 2938-2958] where the target set is the origin in the state space. When the target set is the origin or a ball, centered at the origin, the minimal norm and the minimal time functions are continuous and strictly decreasing, and they are inverses of each other. However, when the target is located in other place of the state space, the minimal norm function may be no longer monotonous and the range of the minimal time function may not be connected. These cause the main difficulty in our study. We overcome this difficulty by borrowing some idea from the classical raising sun lemma (see, for instance, Lemma 3.5 and Figure 5 on Pages 121-122 in [E. M. Stein and R. Shakarchi, \textit{Real Analysis: Measure Theory, Integration, and Hilbert Spaces}, Princeton University Press, 2005]).


title: How Useful is Region-based Classification of Remote Sensing Images in a Deep Learning Framework?

abstract: In this paper, we investigate the impact of segmentation algorithms as a preprocessing step for classification of remote sensing images in a deep learning framework. Especially, we address the issue of segmenting the image into regions to be classified using pre-trained deep neural networks as feature extractors for an SVM-based classifier. An efficient segmentation as a preprocessing step helps learning by adding a spatially-coherent structure to the data. Therefore, we compare algorithms producing superpixels with more traditional remote sensing segmentation algorithms and measure the variation in terms of classification accuracy. We establish that superpixel algorithms allow for a better classification accuracy as a homogenous and compact segmentation favors better generalization of the training samples.


title: Total Order Reliable Convergecast in WBAN

abstract: This paper is the first extensive work on total order reliable convergecast in multi-hop Wireless Body Area Networks (WBAN). Convergecast is a many-to-one cooperative scheme where each node of the network transmits data towards the same sink. Our contribution is threefold. First, we stress existing WBAN convergecast strategies with respect to their capacity to be reliable and to ensure the total order delivery at sink. That is, packets sent in a specific order should be received in the same order by the sink. When stressed with transmission rates up to 500 packets per second the performances of these strategies decrease dramatically (more than 90% of packets lost). Secondly, we propose a new posture-centric model for WBAN. This model offers a good characterization of the path availability which is further used to fine tune the retransmission rate thresholds. Third, based on our model we propose a new mechanism for reliability and a new converge-cast strategy that outperforms WBAN dedicated strategies but also strategies adapted from DTN and WSN areas. Our extensive performance evaluations use essential parameters for WBAN: packet lost, total order reliability (messages sent in a specific order should be delivered in that specific order) and various human body postures. In particular, our strategy ensures zero packet order inversions for various transmission rates and mobility postures. Interestingly, our strategy respects this property without the need of additional energy-guzzler mechanisms.


title: A brief review of E theory

abstract: I begin with some memories of Abdus Salam who was my PhD supervisor. After reviewing the theory of non-linear realisations and Kac-Moody algebras, I explain how to construct the non-linear realisation based on the Kac-Moody algebra E11 and its vector representation. I explain how this field theory leads to dynamical equations which contain an infinite number of fields defined on a space time with an infinite number of coordinates. I then show that these unique dynamical equations, when truncated to low level fields and the usual coordinates of space time, lead to precisely the equations of motion of eleven dimensional supergravity theory. By taking different group decompositions of E11 we find all the maximal supergravity theories, including the gauged maximal supergravities, and as a result the non-linear realisation should be thought of as a unified theory that is the low energy effective action for type II strings and branes. These results essentially confirm the E11 conjecture given many years ago.


title: A probabilistic network for the diagnosis of acute cardiopulmonary diseases

abstract: In this paper, the development of a probabilistic network for the diagnosis of acute cardiopulmonary diseases is presented. This paper is a draft version of the article published after peer review in 2018 (https://doi.org/10.1002/bimj.201600206). A panel of expert physicians collaborated to specify the qualitative part, that is a directed acyclic graph defining a factorization of the joint probability distribution of domain variables. The quantitative part, that is the set of all conditional probability distributions defined by each factor, was estimated in the Bayesian paradigm: we applied a special formal representation, characterized by a low number of parameters and a parameterization intelligible for physicians, elicited the joint prior distribution of parameters from medical experts, and updated it by conditioning on a dataset of hospital patient records using Markov Chain Monte Carlo simulation. Refinement was cyclically performed until the probabilistic network provided satisfactory Concordance Index values for a selection of acute diseases and reasonable inference on six fictitious patient cases. The probabilistic network can be employed to perform medical diagnosis on a total of 63 diseases (38 acute and 25 chronic) on the basis of up to 167 patient findings.


title: Orthogonal Series Estimates on Strong Spatial Mixing Data

abstract: We study a nonparametric regression model for sample data which is defined on an $N$-dimensional lattice structure and which is assumed to be strong spatial mixing: we use design adapted multidimensional Haar wavelets which form an orthonormal system w.r.t. the empirical measure of the sample data. For such orthonormal systems, we consider a nonparametric hard thresholding estimator. We give sufficient criteria for the consistency of this estimator and we derive rates of convergence. The theorems reveal that our estimator is able to adapt to the local smoothness of the underlying regression function and the design distribution. We illustrate our results with simulated examples.


title: Transparent numerical boundary conditions for evolution equations: Derivation and stability analysis

abstract: The aim of this article is to propose a systematic study of transparent boundary conditions for finite difference approximations of evolution equations. We try to keep the discussion at the highest level of generality in order to apply the theory to the broadest class of problems. We deal with two main issues. We first derive transparent numerical boundary conditions, that is, we exhibit the relations satisfied by the solution to the pure Cauchy problem when the initial condition vanishes outside of some domain. Our derivation encompasses discretized transport, diffusion and dispersive equations with arbitrarily wide stencils. The second issue is to prove sharp stability estimates for the initial boundary value problem obtained by enforcing the boundary conditions derived in the first step. We focus here on discretized transport equations. Under the assumption that the numerical boundary is non-characteristic, our main result characterizes the class of numerical schemes for which the corresponding transparent boundary conditions satisfy the so-called Uniform Kreiss-Lopatinskii Condition introduced in [GKS72]. Adapting some previous works to the non-local boundary conditions considered here, our analysis culminates in the derivation of trace and semigroup estimates for such transparent numerical boundary conditions. Several examples and possible extensions are given.


title: Characterizations of projective spaces and quadrics by strictly nef bundles

abstract: In this paper, we show that if the tangent bundle of a smooth projective variety is strictly nef, then it is isomorphic to a projective space; if a projective variety $X^n$ $(n>4)$ has strictly nef $\Lambda^2 TX$, then it is isomorphic to $\mathbb{P}^n$ or quadric $\mathbb{Q}^n$. We also prove that on elliptic curves, strictly nef vector bundles are ample, whereas there exist Hermitian flat and strictly nef vector bundles on any smooth curve with genus $g\geq 2$.


title: An Analysis of Technical Debt Management Through Resources Allocation Policies in Software Maintenance Process

abstract: This paper presents an analysis of technical debt management through resources allocation policies in software maintenance process during its operation to demonstrate how different strategies leads to the emergence of different behaviors along the evolution path. To achieve this objective, this work used the System Dynamic approach for building a computational simulation model based on extensive literature review and secondary data. Most of the works that applied the System Dynamics on software projects research, focused on initial phases of its life cycle, leaving a gap to be explored regarding the long-term behaviors of the operation and maintenance phases. The results demonstrated that the excessive focus on the perfective maintenance activities could be more costly than performing regular preventive maintenance to reduce the technical debt incurred, ending up with fewer functionalities deployed, higher backlog, lower productivity, lower maintainability and higher technical principal.


title: On the denseness of minimum attaining operators

abstract: Let $H_1,H_2$ be complex Hilbert spaces and $T$ be a densely defined closed linear operator (not necessarily bounded). It is proved that for each $\epsilon>0$, there exists a bounded operator $S$ with $\|S\|\leq \epsilon$ such that $T+S$ is minimum attaining. Further, if $T$ is bounded below, then $S$ can be chosen to be rank one.


title: Distributed Training of Deep Neural Networks: Theoretical and Practical Limits of Parallel Scalability

abstract: This paper presents a theoretical analysis and practical evaluation of the main bottlenecks towards a scalable distributed solution for the training of Deep Neuronal Networks (DNNs). The presented results show, that the current state of the art approach, using data-parallelized Stochastic Gradient Descent (SGD), is quickly turning into a vastly communication bound problem. In addition, we present simple but fixed theoretic constraints, preventing effective scaling of DNN training beyond only a few dozen nodes. This leads to poor scalability of DNN training in most practical scenarios.


title: Interplay between optical, viscous and elastic forces on an optically trapped Brownian particle immersed in a viscoelastic fluid

abstract: We provide a detailed study of the interplay between the different interactions which appear in the Brownian motion of a micronsized sphere immersed in a viscoelastic fluid measured with optical trapping interferometry. To explore a wide range of viscous, elastic and optical forces, we analyze two different viscoelastic solutions at various concentrations, which provide a dynamic polymeric structure surrounding the Brownian sphere. Our experiments show that, depending of the fluid, optical forces, even if small, slightly modify the complex modulus at low frequencies. Based on our findings, we propose an alternative methodology to calibrate this kind of experimental set-up when non-Newtonian fluids are used. Understanding the influence of the optical potential is essential for a correct interpretation of the mechanical properties obtained by optically-trapped probe-based studies of biomaterials and living matter.


title: Transformation of the frequency-modulated continuous-wave field into a train of short pulses by resonant filters

abstract: The resonant filtering method transforming frequency modulated radiation field into a train of short pulses is proposed to apply in optical domain. Effective frequency modulation can be achieved by electro-optic modulator or by resonant frequency modulation of the filter with a narrow absorption line. Due to frequency modulation narrow-spectrum CW radiation field is seen by the resonant filter as a comb of equidistant spectral components separated by the modulation frequency. Tuning narrow-bandwidth filter in resonance with $n$-th spectral component of the comb transforms the radiation field into bunches of pulses with $n$ pulses in each bunch. The transformation is explained by the interference of the coherently scattered resonant component of the field with the whole comb. Constructive interference results in formation of pulses, while destructive interference is seen as dark windows between pulses. It is found that the optimal thickness of the resonant filter is several orders of magnitude smaller than the necessary thickness of the dispersive filters used before in optical domain to produce short pulses from the frequency modulated field.


title: On a delay differential equation arising from a car-following model: Wavefront solutions with constant-speed and their stability

abstract: This work is concerned with the study of a scalar delay differential equation \begin{equation*} z^{\prime\prime}(t)=h^2\,V(z(t-1)-z(t))+h\,z^\prime(t) \end{equation*} motivated by a simple car-following model on an unbounded straight line. Here, the positive real $h$ denotes some parameter, and $V$ is a so-called \textit{optimal velocity function} of the traffic model involved. We analyze the existence and local stability properties of solutions $z(t)=c\,t+d$, $t\in\mathbb{R}$, with $c,d\in\mathbb{R}$. In the case $c\not=0$, such a solution of the differential equation forms a wavefront solution of the car-following model where all cars are uniformly spaced on the line and move with the same constant velocity. In particular, it is shown that all but one of these wavefront solutions are located on two branches parametrized by $h$. Furthermore, we prove that along the one branch all solutions are unstable due to the principle of linearized instability, whereas along the other branch some of the solutions may be stable. The last point is done by carrying out a center manifold reduction as the linearization does always have a zero eigenvalue. Finally, we provide some numerical examples demonstrating the obtained analytical results.


title: EEG reconstruction and skull conductivity estimation using a Bayesian model promoting structured sparsity

abstract: M/EEG source localization is an open research issue. To solve it, it is important to have good knowledge of several physical parameters to build a reliable head operator. Amongst them, the value of the conductivity of the human skull has remained controversial. This report introduces a novel hierarchical Bayesian framework to estimate the skull conductivity jointly with the brain activity from the M/EEG measurements to improve the reconstruction quality. A partially collapsed Gibbs sampler is used to draw samples asymptotically distributed according to the associated posterior. The generated samples are then used to estimate the brain activity and the model hyperparameters jointly in a completely unsupervised framework. We use synthetic and real data to illustrate the improvement of the reconstruction. The performance of our method is also compared with two optimization algorithms introduced by Vallagh\'e \textit{et al.} and Gutierrez \textit{et al.} respectively, showing that our method is able to provide results of similar or better quality while remaining applicable in a wider array of situations.


title: Log-convexity and the cycle index polynomials with relation to compound Poisson distributions

abstract: We extend the exponential formula by Bender and Canfield (1996), which relates log-concavity and the cycle index polynomials. The extension clarifies the log-convexity relation. The proof is by noticing the property of a compound Poisson distribution together with its moment generating function. We also give a combinatorial proof of extended "log-convex part" referring Bender and Canfield's approach, where the formula by Bruijn and Erd\"os (1953) is additionally exploited. The combinatorial approach yields richer structural results more than log-convexity. Furthermore, we consider normal and binomial convolutions of sequences which satisfy the exponential formula. The operations generate interesting examples which are not included in well known laws about log-concavity/convexity.


title: Finite frequency Seebeck coefficient of metals: A memory function approach

abstract: We study the dynamical thermoelectric transport in metals subjected to the electron-impurity and the electron-phonon interactions using the memory function formalism. We introduce a generalized Drude form for the Seebeck coefficient in terms of thermoelectric memory function and calculate the later in various temperature and frequency limits. In the zero frequency and high temperature limit, we find that our results are consistent with the experimental findings and with the traditional Boltzmann equation approach. In the low temperature limit, we find that the Seebeck coefficient is quadratic in temperature. In the finite frequency regime, we report new results: In the electron-phonon interaction case, we find that the Seebeck coefficient shows frequency independent behavior both in the high frequency regime ($\omega \gg \omega_{D}$, where $\omega_{D}$ is the Debye frequency) and in the low frequency regime ($\omega \ll \omega_{D}$), whereas in the intermediate frequencies, it is a monotonically increasing function of frequency. In the case of the electron-impurity interaction, first it decays and then after passing through a minimum it increases with the increase in frequency and saturates at high frequencies.


title: Comparison of Channels: Criteria for Domination by a Symmetric Channel

abstract: This paper studies the basic question of whether a given channel $V$ can be dominated (in the precise sense of being more noisy) by a $q$-ary symmetric channel. The concept of "less noisy" relation between channels originated in network information theory (broadcast channels) and is defined in terms of mutual information or Kullback-Leibler divergence. We provide an equivalent characterization in terms of $\chi^2$-divergence. Furthermore, we develop a simple criterion for domination by a $q$-ary symmetric channel in terms of the minimum entry of the stochastic matrix defining the channel $V$. The criterion is strengthened for the special case of additive noise channels over finite Abelian groups. Finally, it is shown that domination by a symmetric channel implies (via comparison of Dirichlet forms) a logarithmic Sobolev inequality for the original channel.