title: Conditioning of implicit Runge-Kutta integration for finite element approximation of linear diffusion equations on anisotropic meshes

abstract: The conditioning of implicit Runge-Kutta (RK) integration for linear finite element approximation of diffusion equations on general anisotropic meshes is investigated. Bounds are established for the condition number of the resulting linear system with and without diagonal preconditioning for the implicit Euler and general implicit RK methods. Two solution strategies are considered for the linear system resulting from general implicit RK integration: the simultaneous solution (the system is solved as a whole) and a successive solution which follows the commonly used implementation of implicit RK methods to first transform the system into smaller systems using the Jordan normal form of the RK matrix and then solve them successively. For the simultaneous solution in case of a positive semidefinite symmetric part of the RK coefficient matrix and for the successive solution it is shown that . If the smallest eigenvalue of the symmetric part of the RK coefficient matrix is negative and the simultaneous solution strategy is used, an upper bound on the time step is given so that the system matrix is positive definite. The obtained bounds for the condition number have explicit geometric interpretations and take the interplay between the diffusion matrix and the mesh geometry into full consideration. They show that there are three mesh-dependent factors that can affect the conditioning: the number of elements, the mesh nonuniformity measured in the Euclidean metric, and the mesh nonuniformity with respect to the inverse of the diffusion matrix. They also reveal that the preconditioning using the diagonal of the system matrix, the mass matrix, or the lumped mass matrix can effectively eliminate the effects of the mesh nonuniformity measured in the Euclidean metric. Numerical examples are given.


title: Optimal strengthening of particle-loaded liquid foams

abstract: Foams made of complex fluids such as particle suspensions have a great potential for the development of advanced aerated materials. In this paper we study the rheological behavior of liquid foams loaded with granular suspensions. We focus on the effect of small particles, i.e. particle-to-bubble size ratio smaller than 0.1, and we measure the complex modulus as a function of particle size and particle volume fraction. With respect to previous work, the results highlight a new elastic regime characterized by unequaled modulus values as well as independence of size ratio. A careful investigation of the material microstructure reveals that particles organize through the network between the gas bubbles and form a granular skeleton structure with tightly packed particles. The latter is proven to be responsible for the reported new elastic regime. Rheological probing performed by strain sweep reveals a two-step yielding of the material: the first one occurs at small strain and is clearly attributed to yielding of the granular skeleton; the second one corresponds to the yielding of the bubble assembly, as observed for particle-free foams. Moreover the elastic modulus measured at small strain is quantitatively described by models for solid foams in assuming that the granular skeleton possesses a bulk elastic modulus of order 100 kPa. Additional rheology experiments performed on the bulk granular material indicate that this surprisingly high value can be understood as soon as the magnitude of the confinement pressure exerted by foam bubbles on packed grains is considered.


title: Existence and uniqueness of steady weak solutions to the Navier-Stokes equations in $\mathbb{R}^2$

abstract: The existence of weak solutions to the stationary Navier-Stokes equations in the whole plane $\mathbb{R}^2$ is proven. This particular geometry was the only case left open since the work of Leray in 1933. The reason is that due to the absence of boundaries the local behavior of the solutions cannot be controlled by the enstrophy in two dimensions. We overcome this difficulty by constructing approximate weak solutions having a prescribed mean velocity on some given bounded set. As a corollary, we obtain infinitely many weak solutions in $\mathbb{R}^2$ parameterized by this mean velocity, which is reminiscent of the expected convergence of the velocity field at large distances to any prescribed constant vector field. This explicit parameterization of the weak solutions allows us to prove a weak-strong uniqueness theorem for small data. The question of the asymptotic behavior of the weak solutions remains however open, when the uniqueness theorem doesn't apply.


title: Adsorption of atoms and molecules on s-triazine sheet with embedded manganese atom: First-principles calculations

abstract: The mechanical, structural, electronic and magnetic properties of s-triazine sheet (C6N6) with embedded Mn atom (Mn-C6N6) is investigated under the influence of external environment using density functional theory. Our results show that Mn-C6N6 system is structurally and mechanically stable. The binding energy of Mn embedded in C6N6 sheet can be modulated under the influence of symmetric deformation and perpendicular electric field respectively. The semiconducting property of pure C6N6 sheet is maintained upon embedment of Mn atom in the porous site. It is also found that small increment in bi-axial tensile strain enhances the band gap (from 0.630 eV at zero strain to 0.802 eV at 5% strain) while the magnetic moment of the embedded Mn atom is preserved. The electronic and magnetic properties of the Mn-C6N6 systems are maintained up to 10 V/nm in electric field strength. We also explore the geometries, electronic and magnetic properties of Mn-C6N6 with adsorbed atoms and molecules. The Mn-C6N6 with adsorbed O atom and O2 molecule systems shows half-metallic character whereas the remaining systems preserve their semiconducting property. The total magnetic moment per unit cell in most of the systems is found to reduce as compared to that of the Mn-C6N6 sheet. The reduction in magnetic moment can be related to the strong interactions among the Mn atom and the surrounding atoms which lead to the formation of low-spin configurations. Overall, our results indicate that the Mn-C6N6 systems with and without adsorbed atoms and molecules can serve as potential candidates for future spintronics and catalysis applications.


title: Il Fattore di Sylvester

abstract: Sylvester factor, an essential part of the asymptotic formula of Hardy and Littlewood which is the extended Goldbach conjecture, regarded as strongly multiplicative arithmetic function, has several remarkable properties.


title: The Carnegie-Chicago Hubble Program. II. The Distance to IC 1613: The Tip of the Red Giant Branch and RR Lyrae Period-Luminosity Relations

abstract: IC 1613 is an isolated dwarf galaxy within the Local Group. Low foreground and internal extinction, low metallicity, and low crowding make it an invaluable testbed for the calibration of the local distance ladder. We present new, high-fidelity distance estimates to IC 1613 via its Tip of the Red Giant Branch (TRGB) and its RR Lyrae (RRL) variables as part of the Carnegie-Chicago Hubble Program, which seeks an alternate local route to \ho using Population II stars. We have measured a TRGB magnitude I=20.35+/-0.01 (statistical)+/-0.01 (systematic) using wide-field observations obtained from the IMACS camera on the Magellan-Baade telescope. We have further constructed optical and near-infrared RRL light curves using archival BI- and new H- band observations from the ACS/WFC and WFC3/IR instruments aboard the Hubble Space Telescope (HST). In advance of future Gaia data releases, we set provisional values for the TRGB luminosity via the Large Magellanic Cloud and Galactic RRL zero-points via HST parallaxes. We find corresponding true distance moduli \mu(TRGB)=24.30+/-0.03 (statistical)+/-0.05 (systematic) and \mu(RRL)=24.28+/-0.04 (statistical+systematic). We compare our results to a body of recent publications on IC 1613 and find no statistically significant difference between the distances derived from stars of Population I and II.


title: On $H^2$-gradient Flows for the Willmore Energy

abstract: We show that the concept of $H^2$-gradient flow for the Willmore energy and other functionals that depend at most quadratically on the second fundamental form is well-defined in the space of immersions of Sobolev class $W^{2,p}$ from a compact, $n$-dimensional manifold into Euclidean space, provided that $p \geq 2$ and $p>n$. We also discuss why this is not true for Sobolev class $H^2=W^{2,2}$. In the case of equality constraints, we provide sufficient conditions for the existence of the projected $H^2$-gradient flow and demonstrate its usability for optimization with several numerical examples.


title: Reconfigurable re-entrant cavity for wireless coupling to an electro-optomechanical device

abstract: An electro-optomechanical device capable of microwave-to-optics conversion has recently been demonstrated, with the vision of enabling optical networks of superconducting qubits. Here we present an improved converter design that uses a three-dimensional (3D) microwave cavity for coupling between the microwave transmission line and an integrated LC resonator on the converter chip. The new design simplifies the optical assembly and decouples it from the microwave part of the setup. Experimental demonstrations show that the modular device assembly allows us to flexibly tune the microwave coupling to the converter chip while maintaining small loss. We also find that electromechanical experiments are not impacted by the additional microwave cavity. Our design is compatible with a high-finesse optical cavity and will improve optical performance.


title: Multi-Timescale, Gradient Descent, Temporal Difference Learning with Linear Options

abstract: Deliberating on large or continuous state spaces have been long standing challenges in reinforcement learning. Temporal Abstraction have somewhat made this possible, but efficiently planing using temporal abstraction still remains an issue. Moreover using spatial abstractions to learn policies for various situations at once while using temporal abstraction models is an open problem. We propose here an efficient algorithm which is convergent under linear function approximation while planning using temporally abstract actions. We show how this algorithm can be used along with randomly generated option models over multiple time scales to plan agents which need to act real time. Using these randomly generated option models over multiple time scales are shown to reduce number of decision epochs required to solve the given task, hence effectively reducing the time needed for deliberation.


title: Reconstruction of $\alpha$-attractor supergravity models of inflation

abstract: In this paper, we apply reconstruction techniques to recover the potential parameters for a particular class of single-field models, the $\alpha$-attractor (supergravity) models of inflation. This also allows to derive the inflaton vacuum expectation value at horizon crossing. We show how to use this value as one of the input variables to constrain the postaccelerated inflationary phase. We assume that the tensor-to-scalar ratio $r$ is of the order of $10^{-3}$ , a level reachable by the expected sensitivity of the next-generation CMB experiments.


title: A directional uncertainty principle for periodic functions

abstract: In this paper we introduce a notion of a directional uncertainty product for multivariate periodic functions. It measures a localization of a function along a particular direction. We study properties of the uncertainty product and give an example of well localized multivariate periodic Parseval wavelet frames.


title: Fast Diffusion of Long Guest Rods in a Lamellar Phase of Short Host Particles

abstract: We investigate the dynamic behavior of long guest rod-like particles immersed in liquid crystalline phases formed by shorter host rods, tracking both guest and host particles by fluorescence microscopy. Counter-intuitively, we evidence that long rods diffuse faster than short rods forming the one-dimensional ordered smectic-A phase. This results from the larger and non-commensurate size of the guest particles as compared to the wavelength of the energy landscape set by the lamellar stack of liquid slabs. The long guest particles are also shown to be still mobile in the crystalline smectic-B phase, as they generate their own voids in the adjacent layers.


title: Integrated cladding-pumped multi-core, few-mode erbium-doped fibre amplifier for space-division multiplexed communications

abstract: Space-division multiplexing (SDM), whereby multiple spatial channels in multimode and multicore optical fibers are used to increase the total transmission capacity per fiber, is being investigated to avert a data capacity crunch and reduce the cost per transmitted bit. With the number of channels employed in SDM transmission experiments continuing to rise, there is a requirement for integrated SDM components that are scalable. Here, we demonstrate a cladding-pumped SDM erbium-doped fiber amplifier (EDFA) that consists of six uncoupled multimode erbium-doped cores. Each core supports three spatial modes, which enables the EDFA to amplify a total of 18 spatial channels simultaneously with a single pump diode and a complexity similar to a single-mode EDFA. The amplifier delivers more than 20-dBm total output power per core and less than 7-dB noise figure over the C-band. This cladding-pumped EDFA enables combined space-division and wavelength-division multiplexed transmission over multiple multimode fiber spans.


title: Practical Coreset Constructions for Machine Learning

abstract: We investigate coresets - succinct, small summaries of large data sets - so that solutions found on the summary are provably competitive with solution found on the full data set. We provide an overview over the state-of-the-art in coreset construction for machine learning. In Section 2, we present both the intuition behind and a theoretically sound framework to construct coresets for general problems and apply it to $k$-means clustering. In Section 3 we summarize existing coreset construction algorithms for a variety of machine learning problems such as maximum likelihood estimation of mixture models, Bayesian non-parametric models, principal component analysis, regression and general empirical risk minimization.


title: A deeper view of the CoRoT-9 planetary system. A small non-zero eccentricity for CoRoT-9b likely generated by planet-planet scattering

abstract: CoRoT-9b is one of the rare long-period (P=95.3 days) transiting giant planets with a measured mass known to date. We present a new analysis of the CoRoT-9 system based on five years of radial-velocity (RV) monitoring with HARPS and three new space-based transits observed with CoRoT and Spitzer. Combining our new data with already published measurements we redetermine the CoRoT-9 system parameters and find good agreement with the published values. We uncover a higher significance for the small but non-zero eccentricity of CoRoT-9b ($e=0.133^{+0.042}_{-0.037}$) and find no evidence for additional planets in the system. We use simulations of planet-planet scattering to show that the eccentricity of CoRoT-9b may have been generated by an instability in which a $\sim 50~M_\oplus$ planet was ejected from the system. This scattering would not have produced a spin-orbit misalignment, so we predict that CoRoT-9b orbit should lie within a few degrees of the initial plane of the protoplanetary disk. As a consequence, any significant stellar obliquity would indicate that the disk was primordially tilted.


title: Breakdown of traditional many-body theories for correlated electrons

abstract: Starting from the (Hubbard) model of an atom, we demonstrate that the uniqueness of the mapping from the interacting to the noninteracting Green's function, $G\to G_0$, is strongly violated, by providing numerous explicit examples of different $G_0$ leading to the same physical $G$. We argue that there are indeed infinitely many such $G_0$, with numerous crossings with the physical solution. We show that this rich functional structure is directly related to the divergence of certain classes of (irreducible vertex) diagrams, with important consequences for traditional many-body physics based on diagrammatic expansions. Physically, we ascribe the onset of these highly non-perturbative manifestations to the progressive suppression of the charge susceptibility induced by the formation of local magnetic moments and/or RVB states in strongly correlated electron systems.


title: Lift of Frobenius and Descent to Constants

abstract: In differential algebra, a proper scheme $X$ defined over an algebraically closed field $K$ with a derivation $\partial$ on it descends to the field of constants $K^\partial$ if $X$ itself lifts the derivation $\partial$. This is a result by A. Buium. Now in the arithmetic case, the notion of a derivation is replaced by the notion of a $\pi$-derivation $\delta$ or equivalently in the flat case, a lift of Frobenius $\phi$. We will show an analogous result in the arithmetic case of equal characteristic. We show our results using the arithmetic analogue of Taylor expansion using Witt vectors.


title: The Dyck and the Preiss separation uniformly

abstract: We are concerned with two separation theorems about analytic sets by Dyck and Preiss, the former involves the positively-defined subsets of the Cantor space and the latter the Borel-convex subsets of finite dimensional Banach spaces. We show by introducing the corresponding separation trees that both of these results admit a constructive proof. This enables us to give the uniform version of the separation theorems, and derive as corollaries the results, which are analogous to the fundamental fact "HYP is effectively bi-analytic" provided by the Souslin-Kleene Theorem.


title: Electron-electron attraction in an engineered electromechanical system

abstract: Two electrons in a quantum dot repel each other: their interaction can be characterized by a positive interaction energy. From the theory of superconductivity, we also know that mechanical vibrations of the crystal lattice can make the electron-electron interaction attractive. Analogously, if a quantum dot interacts with a mechanical degree of freedom, the effective interaction energy can be negative; that is, the electron-electron interaction might be attractive. In this work, we propose and theoretically study an engineered electromechanical system that exhibits electron-electron attraction: a quantum dot suspended on a nonlinear mechanical resonator, tuned by a bottom and a top gate electrode. We focus on the example of a dot embedded in a suspended graphene ribbon, for which we identify conditions for electron-electron attraction. Our results suggest the possibility of electronic transport via tunneling of packets of multiple electrons in such devices, similar to that in superconducting nanostructures, but without the use of any superconducting elements.


title: Locating a robber with multiple probes

abstract: We consider a game in which a cop searches for a moving robber on a connected graph using distance probes, which is a slight variation on one introduced by Seager. Carragher, Choi, Delcourt, Erickson and West showed that for any $n$-vertex graph $G$ there is a winning strategy for the cop on the graph $G^{1/m}$ obtained by replacing each edge of $G$ by a path of length $m$, if $m\geq n$. The present authors showed that, for all but a few small values of $n$, this bound may be improved to $m\geq n/2$, which is best possible. In this paper we consider the natural extension in which the cop probes a set of $k$ vertices, rather than a single vertex, at each turn. We consider the relationship between the value of $k$ required to ensure victory on the original graph and the length of subdivisions required to ensure victory with $k=1$. We give an asymptotically best-possible linear bound in one direction, but show that in the other direction no subexponential bound holds. We also give a bound on the value of $k$ for which the cop has a winning strategy on any (possibly infinite) connected graph of maximum degree $\Delta$, which is best possible up to a factor of $(1-o(1))$.


title: Enhancing Physical Layer Security in Dual-Hop Multiuser Transmission

abstract: In this paper, we consider the Physical Layer Security(PLS) problem in orthogonal frequency division multiple access (OFDMA) based dual-hop system which consists of multiple users, multiple amplify and forward relays, and an eavesdropper. The aim is to enhance PLS of the entire system by maximizing sum secrecy rate of secret users through optimal resource allocation under various practical constraints. Specifically, the sub-carrier allocation to different users, the relay assignments, and the power loading over different sub-carriers at transmitting nodes are optimized. The joint optimization problem is modeled as a mixed binary integer programming problem subject to exclusive sub-carrier allocation and separate power budget constraints at each node. A joint optimization solution is obtained through Lagrangian dual decomposition where KKT conditions are exploited to find the optimal power allocation at base station. Further, to reduce the complexity, a sub-optimal scheme is presented where the optimal power allocation is derived under fixed sub-carrier-relay assignment. Simulation results are also provided to validate the performance of proposed schemes.


title: Characterization theorems for $Q$-independent random variables with values in a locally compact Abelian group

abstract: Let $X$ be a locally compact Abelian group, $Y$ be its character group. Following A. Kagan and G. Sz\'ekely we introduce a notion of $Q$-independence for random variables with values in $X$. We prove group analogues of the Cram\'er, Kac-Bernstein, Skitovich-Darmois and Heyde theorems for $Q$-independent random variables with values in $X$. The proofs of these theorems are reduced to solving some functional equations on the group $Y$.


title: Near Optimal Hamiltonian-Control and Learning via Chattering

abstract: Many applications require solving non-linear control problems that are classically not well behaved. This paper develops a simple and efficient chattering algorithm that learns near optimal decision policies through an open-loop feedback strategy. The optimal control problem reduces to a series of linear optimization programs that can be easily solved to recover a relaxed optimal trajectory. This algorithm is implemented on a real-time enterprise scheduling and control process.


title: Calculations of CR energy spectra within the NoRD model

abstract: Energy spectra of galactic cosmic rays calculated within the framework of the NoRD (nonlocal relativistic diffusion) model are presented. The model accounts for the turbulent character of the interstellar medium and the relativistic speed limit requirement. Calculations account for spiral distribution of sources, boundedness of halo, spallation of nuclei, energy dependence of the diffusion coefficient, tempered power law injection spectrum. We show that the knee in the background spectra from the ensemble of supernovae can arise due to relativistic speed limit requirement. We compare our calculations with the results obtained in frames of the local models (LoD and LoRD) and nonlocal nonrelativistic model (NoD).


title: Anisotropic triangulations via discrete Riemannian Voronoi diagrams

abstract: The construction of anisotropic triangulations is desirable for various applications, such as the numerical solving of partial differential equations and the representation of surfaces in graphics. To solve this notoriously difficult problem in a practical way, we introduce the discrete Riemannian Voronoi diagram, a discrete structure that approximates the Riemannian Voronoi diagram. This structure has been implemented and was shown to lead to good triangulations in $\mathbb{R}^2$ and on surfaces embedded in $\mathbb{R}^3$ as detailed in our experimental companion paper. In this paper, we study theoretical aspects of our structure. Given a finite set of points $\cal P$ in a domain $\Omega$ equipped with a Riemannian metric, we compare the discrete Riemannian Voronoi diagram of $\cal P$ to its Riemannian Voronoi diagram. Both diagrams have dual structures called the discrete Riemannian Delaunay and the Riemannian Delaunay complex. We provide conditions that guarantee that these dual structures are identical. It then follows from previous results that the discrete Riemannian Delaunay complex can be embedded in $\Omega$ under sufficient conditions, leading to an anisotropic triangulation with curved simplices. Furthermore, we show that, under similar conditions, the simplices of this triangulation can be straightened.


title: Eulerian Composition of Certain Franklin Squares

abstract: Several specific Franklin squares and magic squares are decomposed into their quotient and remainder squares. The results support the conjecture that Franklin used the Eulerian composition method to construct many of his squares. This method also can be used to construct new Franklin squares as illustrated herein.


title: Generalized Twisted Quantum Doubles of a Finite Group and Rational Orbifolds

abstract: In previous work the authors introduced a new class of modular quasi-Hopf algebras $D^{\omega}(G, A)$ associated to a finite group $G$, a central subgroup $A$, and a $3$-cocycle $\omega\in Z^3(G, C^x)$. In the present paper we propose a description of the class of orbifold models of rational vertex operator algebras whose module category is tensor equivalent to $D^{\omega}(G, A)$-mod. The paper includes background on quasi-Hopf algebras and a discussion of some relevant orbifolds.


title: Generating Multi-label Discrete Patient Records using Generative Adversarial Networks

abstract: Access to electronic health record (EHR) data has motivated computational advances in medical research. However, various concerns, particularly over privacy, can limit access to and collaborative use of EHR data. Sharing synthetic EHR data could mitigate risk. In this paper, we propose a new approach, medical Generative Adversarial Network (medGAN), to generate realistic synthetic patient records. Based on input real patient records, medGAN can generate high-dimensional discrete variables (e.g., binary and count features) via a combination of an autoencoder and generative adversarial networks. We also propose minibatch averaging to efficiently avoid mode collapse, and increase the learning efficiency with batch normalization and shortcut connections. To demonstrate feasibility, we showed that medGAN generates synthetic patient records that achieve comparable performance to real data on many experiments including distribution statistics, predictive modeling tasks and a medical expert review. We also empirically observe a limited privacy risk in both identity and attribute disclosure using medGAN.


title: Gestalt Phenomenon in Music? A Neurocognitive Physics Study with EEG

abstract: The term gestalt has been widely used in the field of psychology which defined the perception of human mind to group any object not in part but as a unified whole. Music in general is polytonic i.e. a combination of a number of pure tones (frequencies) mixed together in a manner that sounds harmonius. The study of human brain response due to different frequency groups of acoustic signal can give us an excellent insight regarding the neural and functional architecture of brain functions. In this work we have tried to analyze the effect of different frequency bands of music on the various frequency rhythms of human brain obtained from EEG data of 5 participants. Four (4) widely popular Rabindrasangeet clips were subjected to Wavelet Transform method for extracting five resonant frequency bands from the original music signal. These resonant frequency bands were presented to the subjects as auditory stimulus and EEG signals recorded simultaneously in 19 different locations of the brain. The recorded EEG signals were noise cleaned and subjected to Multifractal Detrended Fluctuation Analysis (MFDFA) technique on the alpha, theta and gamma frequency range. Thus, we obtained the complexity values (in the form of multifractal spectral width) in alpha, theta and gamma EEG rhythms corresponding to different frequency bands of music. We obtain frequency specific arousal based response in different lobes of brain as well as in specific EEG bands corresponding to musical stimuli. This revelation can be of immense importance when it comes to the field of cognitive music therapy.


title: VQABQ: Visual Question Answering by Basic Questions

abstract: Taking an image and question as the input of our method, it can output the text-based answer of the query question about the given image, so called Visual Question Answering (VQA). There are two main modules in our algorithm. Given a natural language question about an image, the first module takes the question as input and then outputs the basic questions of the main given question. The second module takes the main question, image and these basic questions as input and then outputs the text-based answer of the main question. We formulate the basic questions generation problem as a LASSO optimization problem, and also propose a criterion about how to exploit these basic questions to help answer main question. Our method is evaluated on the challenging VQA dataset and yields state-of-the-art accuracy, 60.34% in open-ended task.


title: Subsurface bending and reorientation of tilted vortex lattices in the bulk due to Coulomb-like repulsion at the surface

abstract: We study vortex lattices (VLs) in superconducting weak-pinning platelet-like crystals of $\beta$-Bi$_{2}$Pd in tilted fields with a Scanning Tunneling Microscope. We show that vortices exit the sample perpendicular to the surface and are thus bent beneath the surface. The structure and orientation of tilted VL in the bulk are, for large tilt angles, strongly affected by Coulomb-type intervortex repulsion at the surface due to stray fields.


title: Coupling parallel adaptive mesh refinement with a nonoverlapping domain decomposition solver

abstract: We study the effect of adaptive mesh refinement on a parallel domain decomposition solver of a linear system of algebraic equations. These concepts need to be combined within a parallel adaptive finite element software. A prototype implementation is presented for this purpose. It uses adaptive mesh refinement with one level of hanging nodes. Two and three-level versions of the Balancing Domain Decomposition based on Constraints (BDDC) method are used to solve the arising system of algebraic equations. The basic concepts are recalled and components necessary for the combination are studied in detail. Of particular interest is the effect of disconnected subdomains, a typical output of the employed mesh partitioning based on space-filling curves, on the convergence and solution time of the BDDC method. It is demonstrated using a large set of experiments that while both refined meshes and disconnected subdomains have a negative effect on the convergence of BDDC, the number of iterations remains acceptable. In addition, scalability of the three-level BDDC solver remains good on up to a few thousands of processor cores. The largest presented problem using adaptive mesh refinement has over 10^9 unknowns and is solved on 2048 cores.


title: Acceleration of Convergence of Some Infinite Sequences $\boldsymbol{\{A_n\}}$ Whose Asymptotic Expansions Involve Fractional Powers of $\boldsymbol{n}$ via the $\tilde{d}^{(m)}$ transformation

abstract: In this paper, we discuss the application of the author's $\tilde{d}^{(m)}$ transformation to accelerate the convergence of infinite series $\sum^\infty_{n=1}a_n$ when the terms $a_n$ have asymptotic expansions that can be expressed in the form $$ a_n\sim(n!)^{s/m}\exp\left[\sum^{m}_{i=0}q_in^{i/m}\right]\sum^\infty_{i=0}w_i n^{\gamma-i/m}\quad\text{as $n\to\infty$},\quad s\ \text{integer.}$$ We discuss the implementation of the $\tilde{d}^{(m)}$ transformation via the recursive W-algorithm of the author. We show how to apply this transformation and how to assess in a reliable way the accuracies of the approximations it produces, whether the series converge or they diverge. We classify the different cases that exhibit unique numerical stability issues in floating-point arithmetic. We show that the $\tilde{d}^{(m)}$ transformation can also be used efficiently to accelerate the convergence of infinite products $\prod^\infty_{n=1}(1+v_n)$, where $v_n\sim \sum^\infty_{i=0}e_in^{-t/m-i/m}$ as $n\to\infty$,\ $t\geq m+1$ an integer. Finally, we give several numerical examples that attest the high efficiency of the $\tilde{d}^{(m)}$ transformation for the different cases.


title: A Graph-Based Approach to Analyze Flux-Balanced Pathways in Metabolic Networks

abstract: An Elementary Flux Mode (EFM) is a pathway with minimum set of reactions that are functional in steady-state constrained space. Due to the high computational complexity of calculating EFMs, different approaches have been proposed to find these flux-balanced pathways. In this paper, an approach to find a subset of EFMs is proposed based on a graph data model. The given metabolic network is mapped to the graph model and decisions for reaction inclusion can be made based on metabolites and their associated reactions. This notion makes the approach more convenient to categorize the output pathways. Implications of the proposed method on metabolic networks are discussed.


title: A thermodynamic analysis of the spider silk and the importance of complexity

abstract: The spider silk is one of the most interesting bio-materials investigated in the last years. One of the main reasons that brought scientists to study this organized system is its high level of resistance if compared to other artificial materials characterized by higher density. Subsequently, researchers discovered that the spider silk is a complex system formed by different kinds of proteins, organized (or disorganized) to guarantee the required resistance, which is function of the final application and of the environmental conditions. Some spider species are able to make different silks, up to twelve, having a composition that seems to be function of the final use (i.e. dragline web, capture web, etc). The aim of this paper is to analyze the properties of the spider silk by means of a thermodynamic approach, taking advantage of the well-known theories applied to polymers, and to try to underline and develop some intriguing considerations. Moreover, this study can be taken as an example to introduce and discuss the importance of the concept of optionality and of the anti-fragile systems proposed by N. N. Thaleb in his book "Antifragile: Things that gain from disorder".


title: A Polya-Vinogradov-type inequality on $\mathbb{Z}[i]$

abstract: We establish a Polya-Vinogradov-type bound for finite periodic multipicative characters on the Gaussian integers.


title: Image denoising by median filter in wavelet domain

abstract: The details of an image with noise may be restored by removing noise through a suitable image de-noising method. In this research, a new method of image de-noising based on using median filter (MF) in the wavelet domain is proposed and tested. Various types of wavelet transform filters are used in conjunction with median filter in experimenting with the proposed approach in order to obtain better results for image de-noising process, and, consequently to select the best suited filter. Wavelet transform working on the frequencies of sub-bands split from an image is a powerful method for analysis of images. According to this experimental work, the proposed method presents better results than using only wavelet transform or median filter alone. The MSE and PSNR values are used for measuring the improvement in de-noised images.


title: P-N junctions dynamics in graphene channel induced by ferroelectric domains motion

abstract: The p-n junctions dynamics in graphene channel induced by stripe domains nucleation, motion and reversal in a ferroelectric substrate is explored using self-consistent approach based on Landau-Ginzburg-Devonshire phenomenology combined with classical electrostatics. We revealed the extrinsic size effect in the dependence of the graphene channel conductivity on its length. For the case of perfect electric contact between the ferroelectric and graphene, relatively low gate voltages are required to induce the pronounced hysteresis of ferroelectric polarization and graphene charge in response to the periodic gate voltage. Pronounced nonlinear hysteresis of graphene conductance with a wide memory window corresponds to high amplitudes of gate voltage. We predict that the considered nano-structure "top gate/dielectric layer/graphene channel/ferroelectric substrate" can be a promising candidate for the fabrication of new generation of modulators and rectifiers based on the graphene p-n junctions.


title: M\'etodos de Otimiza\c{c}\~ao Combinat\'oria Aplicados ao Problema de Compress\~ao MultiFrases

abstract: The Internet has led to a dramatic increase in the amount of available information. In this context, reading and understanding this flow of information have become costly tasks. In the last years, to assist people to understand textual data, various Natural Language Processing (NLP) applications based on Combinatorial Optimization have been devised. However, for Multi-Sentences Compression (MSC), method which reduces the sentence length without removing core information, the insertion of optimization methods requires further study to improve the performance of MSC. This article describes a method for MSC using Combinatorial Optimization and Graph Theory to generate more informative sentences while maintaining their grammaticality. An experiment led on a corpus of 40 clusters of sentences shows that our system has achieved a very good quality and is better than the state-of-the-art.


title: Energy transfer between modes in a nonlinear beam equation

abstract: We consider the nonlinear nonlocal beam evolution equation introduced by Woinowsky- Krieger. We study the existence and behavior of periodic solutions: these are called nonlinear modes. Some solutions only have two active modes and we investigate whether there is an energy transfer between them. The answer depends on the geometry of the energy function which, in turn, depends on the amount of compression compared to the spatial frequencies of the involved modes. Our results are complemented with numerical experiments, overall, they give a complete picture of the instabilities that may occur in the beam. We expect these results to hold also in more complicated dynamical system


title: CLTune: A Generic Auto-Tuner for OpenCL Kernels

abstract: This work presents CLTune, an auto-tuner for OpenCL kernels. It evaluates and tunes kernel performance of a generic, user-defined search space of possible parameter-value combinations. Example parameters include the OpenCL workgroup size, vector data-types, tile sizes, and loop unrolling factors. CLTune can be used in the following scenarios: 1) when there are too many tunable parameters to explore manually, 2) when performance portability across OpenCL devices is desired, or 3) when the optimal parameters change based on input argument values (e.g. matrix dimensions). The auto-tuner is generic, easy to use, open-source, and supports multiple search strategies including simulated annealing and particle swarm optimisation. CLTune is evaluated on two GPU case-studies inspired by the recent successes in deep learning: 2D convolution and matrix-multiplication (GEMM). For 2D convolution, we demonstrate the need for auto-tuning by optimizing for different filter sizes, achieving performance on-par or better than the state-of-the-art. For matrix-multiplication, we use CLTune to explore a parameter space of more than two-hundred thousand configurations, we show the need for device-specific tuning, and outperform the clBLAS library on NVIDIA, AMD and Intel GPUs.


title: Detection of adsorbed transition-metal porphyrins by spin-dependent conductance of graphene nanoribbon

abstract: Electronic transport in a zig-zag-edge graphene nanoribbon (GNR) and its modification by adsorbed transition metal porphyrins is studied by means of density functional theory calculations. The detachment reaction of the metal centre of the porphyrin is investigated both in the gas phase and for molecules adsorbed on the GNR. As most metal porphyrins are very stable against this reaction, it is found that these molecules bind only weakly to a perfect nanoribbon. However, interaction with a single-atom vacancy in the GNR results in chemical bonding by the transition metal centre being shared between nitrogen atoms in the porphyrin ring and the carbon atoms next to the vacancy in the GNR. For both the physisorbed and the chemisorbed geometry, the inclusion of van der Waals interaction results in a significant enlargement of the binding energy and reduction of the adsorption height. Electronic transport calculations using non-equilibrium Greens functions show that the conductivity of the GNR is altered by the chemisorbed porphyrin molecules. Since the metal centers of porphyrins carry an element-specific magnetic moment, not only the net conductance, but also the spin-dependent conductance of the GNR is affected. In particular, the adsorption of Ru-porphyrin on the single-atom vacancy results in a very large spin polarization of the current of 88% at small applied source-drain voltages. Based on our results, we suggest that a spin valve constructed from a GNR with ferromagnetic contacts could be used as a sensitive detector that could discriminate between various metal porphyrins.


title: Multiloop functional renormalization group that sums up all parquet diagrams

abstract: We present a multiloop flow equation for the four-point vertex in the functional renormalization group (fRG) framework. The multiloop flow consists of successive one-loop calculations and sums up all parquet diagrams to arbitrary order. This provides substantial improvement of fRG computations for the four-point vertex and, consequently, the self-energy. Using the X-ray-edge singularity as an example, we show that solving the multiloop fRG flow is equivalent to solving the (first-order) parquet equations and illustrate this with numerical results.


title: Electroweak Baryogenesis and Dark Matter via a Pseudoscalar vs. Scalar

abstract: We study the electroweak baryogenesis in a fermionic dark matter scenario with a (pseudo)scalar being the mediator in the Higgs portal. It is discussed that the electroweak phase transition turns to be first-order after taking into account the role of the (pseudo)scalar in the thermal effective potential in our extended standard model. Imposing the relic density constraint from the WMAP/Planck and the bounds from the direct detection experiments XENON100/LUX, we show that the dark matter scenario with a scalar mediator is hardly capable of explaining the baryogenesis while the same model with a pseudoscalar mediator is able to explain the baryon asymmetry. For the latter, we constrain more the model with {\it Fermi}-LAT upper limit on dark matter annihilation into $b\bar b$ and $\tau^+\tau^-$. The allowed dark matter mass that leads to correct relic abundance, renders the electroweak phase transition strongly first-order, and respects the {\it Fermi}-LAT limit, will be in the range $110-320$ GeV. The exotic and invisible Higgs decay bounds and the mono-jet search limit at the LHC do not affect the viable space of parameters.


title: Discontinuities in the Electromagnetic Fields of Vortex Beams in the Complex Source/Sink Model

abstract: An analytical discontinuity is reported in what was thought to be the discontinuity-free exact nonparaxial vortex beam phasor obtained within the complex source/sink model. This discontinuity appears for all odd values of the orbital angular momentum mode. Such discontinuities in the phasor lead to nonphysical discontinuities in the real electromagnetic field components. We identify the source of the discontinuities, and provide graphical evidence of the discontinuous real electric fields for the first and third orbital angular momentum modes. A simple means of avoiding these discontinuities is presented.


title: Distinct stick-slip modes in adhesive polymer interfaces

abstract: Stick-slip, manifest as intermittent tangential motion between two solids, is a well-known friction instability that occurs in a number of natural and engineering systems. In the context of adhesive polymer interfaces, this phenomenon has often been solely associated with Schallamach waves, which are termed slow waves due to their low propagation speeds. We study the dynamics of a model polymer interface using coupled force measurements and high speed \emph{in situ} imaging, to explore the occurrence of stick-slip linked to other slow wave phenomena. Two new waves---slip pulse and separation pulse---both distinct from Schallamach waves, are described. The slip pulse is a sharp stress front that propagates in the same direction as the Schallamach wave, while the separation pulse involves local interface detachment and travels in the opposite direction. Transitions between these stick-slip modes are easily effected by changing the sliding velocity or normal load. The properties of these three waves, and their relation to stick-slip is elucidated. We also demonstrate the important role of adhesion in effecting wave propagation.


title: Parametric Polynomial Preserving Recovery on Manifolds

abstract: This paper investigates gradient recovery schemes for data defined on discretized manifolds. The proposed method, parametric polynomial preserving recovery (PPPR), does not require the tangent spaces of the exact manifolds, and they have been assumed for some significant gradient recovery methods in the literature. Another advantage of PPPR is that superconvergence is guaranteed without the symmetric condition which has been asked in the existing techniques. There is also numerical evidence that the superconvergence by PPPR is high curvature stable, which distinguishes itself from the others. As an application, we show its capability of constructing an asymptotically exact \textit{a posteriori} error estimator. Several numerical examples on two-dimensional surfaces are presented to support the theoretical results and comparisons with existing methods are documented.


title: Constructive tensor field theory: The $T^{4}_{4}$ model

abstract: We continue our constructive study of tensor field theory through the next natural model, namely the rank four tensor theory with quartic melonic interactions and propagator inverse of the Laplacian on $U(1)^4$. This superrenormalizable tensor field theory has a power counting quite similar to ordinary $\phi^4_3$. We control the model via a multiscale loop vertex expansion which has to be pushed quite beyond the one of the $T^4_3$ model and we establish its Borel summability in the coupling constant. This paper is also a step to prepare the constructive treatment of just renormalizable models, such as the $T^4_5$ model with quartic melonic interactions.


title: Global surface temperature trends and the effect of World War II: a parametric analysis (long version)

abstract: We fit 12 independent surface temperature time series (zones of latitude), 6 land- and 6 ocean-based, with a Gaussian (centered around WW2) on a quadratic background. The four polar zones are unusable. Each of the remaining 8 shows evidence of a 0.2-0.3 degree C bump, although one zone is dropped from the final analysis because of an aberrant fit. Thus, we find that a surface temperature increase around WW2 is a robust, statistically significant (5.2 sigma) feature of the data, and must either be a coincidence, due to systematic error that came and went around WW2, or direct evidence of the effect of human activity (an increase in combustion of various sorts) on the global climate. The fact that, at least in that era, the temperature recovered quickly, is encouraging. On the other hand, extrapolation of the background temperature rise predicts an increase of 0.5 degree C in the next 20 years if conditions remain the same.


title: Quantum-Chaotic Cryptography

abstract: In this work, it is presented an optical scheme for quantum key distribution employing two synchronized optoelectronic oscillators (OEO) working in the chaotic regime. The produced key depends on the chaotic dynamic and the synchronization between Alice's and Bob's OEOs uses quantum states. An attack on the synchronization signals will disturb the synchronization of the chaotic systems increasing the error rate in the final key.


title: Bernoulli Rank-$1$ Bandits for Click Feedback

abstract: The probability that a user will click a search result depends both on its relevance and its position on the results page. The position based model explains this behavior by ascribing to every item an attraction probability, and to every position an examination probability. To be clicked, a result must be both attractive and examined. The probabilities of an item-position pair being clicked thus form the entries of a rank-$1$ matrix. We propose the learning problem of a Bernoulli rank-$1$ bandit where at each step, the learning agent chooses a pair of row and column arms, and receives the product of their Bernoulli-distributed values as a reward. This is a special case of the stochastic rank-$1$ bandit problem considered in recent work that proposed an elimination based algorithm Rank1Elim, and showed that Rank1Elim's regret scales linearly with the number of rows and columns on "benign" instances. These are the instances where the minimum of the average row and column rewards $\mu$ is bounded away from zero. The issue with Rank1Elim is that it fails to be competitive with straightforward bandit strategies as $\mu \rightarrow 0$. In this paper we propose Rank1ElimKL which simply replaces the (crude) confidence intervals of Rank1Elim with confidence intervals based on Kullback-Leibler (KL) divergences, and with the help of a novel result concerning the scaling of KL divergences we prove that with this change, our algorithm will be competitive no matter the value of $\mu$. Experiments with synthetic data confirm that on benign instances the performance of Rank1ElimKL is significantly better than that of even Rank1Elim, while experiments with models derived from real data confirm that the improvements are significant across the board, regardless of whether the data is benign or not.


title: Recurrent Collective Classification

abstract: We propose a new method for training iterative collective classifiers for labeling nodes in network data. The iterative classification algorithm (ICA) is a canonical method for incorporating relational information into classification. Yet, existing methods for training ICA models rely on the assumption that relational features reflect the true labels of the nodes. This unrealistic assumption introduces a bias that is inconsistent with the actual prediction algorithm. In this paper, we introduce recurrent collective classification (RCC), a variant of ICA analogous to recurrent neural network prediction. RCC accommodates any differentiable local classifier and relational feature functions. We provide gradient-based strategies for optimizing over model parameters to more directly minimize the loss function. In our experiments, this direct loss minimization translates to improved accuracy and robustness on real network data. We demonstrate the robustness of RCC in settings where local classification is very noisy, settings that are particularly challenging for ICA.


title: Fractal Weyl laws and wave decay for general trapping

abstract: We prove a Weyl upper bound on the number of scattering resonances in strips for manifolds with Euclidean infinite ends. In contrast with previous results, we do not make any strong structural assumptions on the geodesic flow on the trapped set (such as hyperbolicity) and instead use propagation statements up to the Ehrenfest time. By a similar method we prove a decay statement with high probability for linear waves with random initial data. The latter statement is related heuristically to the Weyl upper bound. For geodesic flows with positive escape rate, we obtain a power improvement over the trivial Weyl bound and exponential decay up to twice the Ehrenfest time.


title: Scale Magnetic Effect in Quantum Electrodynamics and the Wigner-Weyl Formalism

abstract: The Scale Magnetic Effect (SME) is the generation of electric current due to conformal anomaly in external magnetic field in curved spacetime. The effect appears in a vacuum with electrically charged massless particles. Similarly to the Hall effect, the direction of the induced anomalous current is perpendicular to the direction of the external magnetic field $\bf B$ and to the gradient of the conformal factor $\tau$, while the strength of the current is proportional to the beta function of the theory. In massive electrodynamics the SME remains valid, but the value of the induced current differs from the current generated in the system of massless fermions. In the present paper we use the Wigner--Weyl formalism to demonstrate that in accordance with the decoupling property of heavy fermions the corresponding anomalous conductivity vanishes in the large-mass limit with $m^2 \gg |e {\bf B}|$ and $m \gg |\nabla \tau|$.


title: Assessment of two hybrid van der Waals density functionals for covalent and non-covalent binding of molecules

abstract: Two hybrid van der Waals density functionals (vdW-DFs) are constructed using 25%, Fock exchange with i) the consistent-exchange vdW-DF-cx functional and ii) with the vdW-DF2 functional. The ability to describe covalent and non-covalent binding properties of molecules are assessed. For properties related to covalent binding, atomization energies (G2-1 set), molecular reaction energies (G2RC set), as well as ionization energies (G21IP set) are benchmarked against experimental reference values. We find that hybrid-vdW-DF-cx yields results that are rather similar to those of the standard non-empirical hybrid PBE0 [JCP 110, 6158 (1996)]. Hybrid vdW-DF2 follows somewhat different trends, showing on average significantly larger deviations from the reference energies, with a MAD of 14.5 kcal/mol for the G2-1 set. Non-covalent binding properties of molecules are assessed using the S22 benchmark set of non-covalently bonded dimers and the X40 set of dimers of small halogenated molecules, using wavefunction-based quantum chemistry results for references. For the S22 set, hybrid-vdW-DF-cx performs better than standard vdW-DF-cx for the mostly hydrogen-bonded systems. Hybrid-vdW-DF2 offers a slight improvement over standard vdW-DF2. Similar trends are found for the X40 set, with hybrid-vdW-DF-cx performing particularly well for binding involving the strongly polar hydrogen halides, but poorly for systems with tiny binding energies. Our study of the X40 set reveals both the potential of mixing Fock exchange with vdW-DF, but also highlights shortcomings of the hybrids constructed here. The solid performance of hybrid-vdW-DF-cx for covalent-bonded systems, as well as the strengths and issues uncovered for non-covalently bonded systems, makes this study a good starting point for developing even more precise hybrid vdW-DFs.


title: Optimal Quantization For Mixed Distributions

abstract: The basic goal of quantization for probability distribution is to reduce the number of values, which is typically uncountable, describing a probability distribution to some finite set and thus approximation of a continuous probability distribution by a discrete distribution. Mixed distributions are an exciting new area for optimal quantization. In this paper, we have determined the optimal sets of $n$-means, the $n$th quantization errors, and the quantization dimensions of different mixed distributions. Besides, we have discussed whether the quantization coefficients for the mixed distributions exist. The results in this paper will give a motivation and insight into more general problems in quantization for mixed distributions.


title: Convergence of Diffusion Generated Motion to Motion by Mean Curvature

abstract: We provide a new proof of convergence to motion by mean curvature (MMC) for the Merriman-Bence-Osher (MBO) thresholding algorithm. The proof is elementary and does not rely on maximum principle for the scheme. The strategy is to construct a natural ansatz of the solution and then estimate the error. The proof thus also provides a convergence rate. Only some weak integrability assumptions of the heat kernel, but not its positivity, is used. Currently the result is proved in the case when smooth and classical solution of MMC exists.


title: Detecting Oriented Text in Natural Images by Linking Segments

abstract: Most state-of-the-art text detection methods are specific to horizontal Latin text and are not fast enough for real-time applications. We introduce Segment Linking (SegLink), an oriented text detection method. The main idea is to decompose text into two locally detectable elements, namely segments and links. A segment is an oriented box covering a part of a word or text line; A link connects two adjacent segments, indicating that they belong to the same word or text line. Both elements are detected densely at multiple scales by an end-to-end trained, fully-convolutional neural network. Final detections are produced by combining segments connected by links. Compared with previous methods, SegLink improves along the dimensions of accuracy, speed, and ease of training. It achieves an f-measure of 75.0% on the standard ICDAR 2015 Incidental (Challenge 4) benchmark, outperforming the previous best by a large margin. It runs at over 20 FPS on 512x512 images. Moreover, without modification, SegLink is able to detect long lines of non-Latin text, such as Chinese.


title: Log-Canonical Coordinates for Poisson Brackets and Rational Changes of Coordinates

abstract: Goodearl and Launois have shown that for a log-canonical Poisson bracket on affine space there is no rational change of coordinates for which the Poisson bracket is constant. Our main result is that if affine space is given a log-canonical Poisson bracket, then there does not exist any rational change of coordinates for which the Poisson bracket is linear. Hence, log-canonical coordinates can be thought of as the simplest possible algebraic coordinates for affine space with a log-canonical coordinate system. In proving this conjecture we find certain invariants of log-canonical Poisson brackets on affine space which linear Poisson brackets do not have.


title: Late-time cosmological evolution of a general class of f(R, T) gravity with minimal curvature-matter coupling

abstract: In this work, we study the late-time cosmological solutions of f(R,T)=g(R)+h(-T) models assuming that the conservation of the energy-momentum tensor (EMT) is violated. We perform our analysis through constructing an autonomous dynamical system for the equations of motion. We study the stability properties of solutions via considering linear perturbations about the related equilibrium points. These parameters which are constructed out of the functions g(R) and h(-T) play the main role in finding the late time behavior of the solutions.


title: Local bulk physics from intersecting modular Hamiltonians

abstract: We show that bulk quantities localized on a minimal surface homologous to a boundary region correspond in the CFT to operators that commute with the modular Hamiltonian associated with the boundary region. If two such minimal surfaces intersect at a point in the bulk then CFT operators which commute with both extended modular Hamiltonians must be localized at the intersection point. We use this to construct local bulk operators purely from CFT considerations, without knowing the bulk metric, using intersecting modular Hamiltonians. For conformal field theories at zero and finite temperature the appropriate modular Hamiltonians are known explicitly and we recover known expressions for local bulk observables.


title: Uniform bounds for rational points on complete intersections of two quadric surfaces

abstract: We give uniform upper bounds for the number of rational points of height at most $B$ on non-singular complete intersections of two quadrics in $\mathbb{P}^3$ defined over $\mathbb{Q}$. To do this, we combine determinant methods with descent arguments.


title: Quantum Brownian motion in an analog Friedmann-Robertson-Walker geometry

abstract: In this paper we study the effects of quantum scalar field vacuum fluctuations on scalar test particles in an analog model for the Friedmann-Robertson-Walker spatially flat geometry. In this scenario, the cases with one and two perfectly reflecting plane boundaries are considered as well the case without boundary. We find that the particles can undergo Brownian motion with a nonzero mean squared velocity induced by the quantum vacuum fluctuations due to the time dependent background and the presence of the boundaries. Typical singularities which appears due to the presence of the boundaries in flat spacetime can be naturally regularized for an asymptotically bounded expanding scale function. Thus, shifts in the velocity could be, at least in principle, detectable experimentally. The possibility to implement this observation in an analog cosmological model by the use of a Bose-Einstein condensate is also discussed.


title: On Optimal 2- and 3-Planar Graphs

abstract: A graph is $k$-planar if it can be drawn in the plane such that no edge is crossed more than $k$ times. While for $k=1$, optimal $1$-planar graphs, i.e., those with $n$ vertices and exactly $4n-8$ edges, have been completely characterized, this has not been the case for $k \geq 2$. For $k=2,3$ and $4$, upper bounds on the edge density have been developed for the case of simple graphs by Pach and T\'oth, Pach et al. and Ackerman, which have been used to improve the well-known "Crossing Lemma". Recently, we proved that these bounds also apply to non-simple $2$- and $3$-planar graphs without homotopic parallel edges and self-loops. In this paper, we completely characterize optimal $2$- and $3$-planar graphs, i.e., those that achieve the aforementioned upper bounds. We prove that they have a remarkably simple regular structure, although they might be non-simple. The new characterization allows us to develop notable insights concerning new inclusion relationships with other graph classes.


title: Vision-based Real-Time Aerial Object Localization and Tracking for UAV Sensing System

abstract: The paper focuses on the problem of vision-based obstacle detection and tracking for unmanned aerial vehicle navigation. A real-time object localization and tracking strategy from monocular image sequences is developed by effectively integrating the object detection and tracking into a dynamic Kalman model. At the detection stage, the object of interest is automatically detected and localized from a saliency map computed via the image background connectivity cue at each frame; at the tracking stage, a Kalman filter is employed to provide a coarse prediction of the object state, which is further refined via a local detector incorporating the saliency map and the temporal information between two consecutive frames. Compared to existing methods, the proposed approach does not require any manual initialization for tracking, runs much faster than the state-of-the-art trackers of its kind, and achieves competitive tracking performance on a large number of image sequences. Extensive experiments demonstrate the effectiveness and superior performance of the proposed approach.


title: Universal Consistency and Robustness of Localized Support Vector Machines

abstract: The massive amount of available data potentially used to discover patters in machine learning is a challenge for kernel based algorithms with respect to runtime and storage capacities. Local approaches might help to relieve these issues. From a statistical point of view local approaches allow additionally to deal with different structures in the data in different ways. This paper analyses properties of localized kernel based, non-parametric statistical machine learning methods, in particular of support vector machines (SVMs) and methods close to them. We will show there that locally learnt kernel methods are universal consistent. Furthermore, we give an upper bound for the maxbias in order to show statistical robustness of the proposed method.


title: The Structure of Extreme Level Sets in Branching Brownian Motion

abstract: We study the structure of extreme level sets of a standard one dimensional branching Brownian motion, namely the sets of particles whose height is within a fixed distance from the order of the global maximum. It is well known that such particles congregate at large times in clusters of order-one genealogical diameter around local maxima which form a Cox process in the limit. We add to these results by finding the asymptotic size of extreme level sets and the typical height and shape of those clusters which carry such level sets. We also find the right tail decay of the distribution of the distance between the two highest particles. These results confirm two conjectures of Brunet and Derrida, 2011. The proofs rely on studying the cluster distribution and should carry over to the branching random walk and the two-dimensional discrete Gaussian free field with no conceptual difficulty.


title: A multi-Frey approach to Fermat equations of signature $(r,r,p)$

abstract: In this paper, we give a resolution of the generalized Fermat equations $$x^5 + y^5 = 3 z^n \text{ and } x^{13} + y^{13} = 3 z^n,$$ for all integers $n \ge 2$, and all integers $n \ge 2$ which are not a multiple of $7$, respectively, using the modular method with Frey elliptic curves over totally real fields. The results require a refined application of the multi-Frey technique, which we show to be effective in new ways to reduce the bounds on the exponents $n$. We also give a number of results for the equations $x^5 + y^5 = d z^n$, where $d = 1, 2$, under additional local conditions on the solutions. This includes a result which is reminiscent of the second case of Fermat's Last Theorem, and which uses a new application of level raising at $p$ modulo $p$.


title: Optimally solving the joint order batching and picker routing problem

abstract: In this work we investigate the problem of order batching and picker routing in storage areas. These are labour and capital intensive problems, often responsible for a substantial share of warehouse operating costs. In particular, we consider the case of online grocery shopping in which orders may be composed of dozens of items. We present a formulation for the problem based on an exponential number of connectivity constraints and we introduce a significant number of valid inequalities based on the standard layout of warehouses, composed of parallel aisles and two or more cross-aisles. The proposed inequalities are highly effective and greatly improve computational results. Instances involving up to 20 orders are solved to proven optimality when we jointly consider order batching and picker routing. Instances involving up to 5000 orders are considered where order batching is done heuristically, but picker routing is done optimally.


title: Stochastic fluctuations can reveal the feedback signs of gene regulatory networks at the single-molecule level

abstract: Understanding the relationship between spontaneous stochastic fluctuations and the topology of the underlying gene regulatory network is of fundamental importance for the study of single-cell stochastic gene expression. Here by solving the analytical steady-state distribution of the protein copy number in a general kinetic model of stochastic gene expression with nonlinear feedback regulation, we reveal the relationship between stochastic fluctuations and feedback topology at the single-molecule level, which provides novel insights into how and to what extent a feedback loop can enhance or suppress molecular fluctuations. Based on such relationship, we also develop an effective method to extract the topological information of a gene regulatory network from single-cell gene expression data. The theory is demonstrated by numerical simulations and, more importantly, validated quantitatively by single-cell data analysis of a synthetic gene circuit integrated in human kidney cells.

categories: q-bio.MN cond-mat.stat-mech physics.bio-ph q-bio.QM


title: Holographic coherent states from random tensor networks

abstract: Random tensor networks provide useful models that incorporate various important features of holographic duality. A tensor network is usually defined for a fixed graph geometry specified by the connection of tensors. In this paper, we generalize the random tensor network approach to allow quantum superposition of different spatial geometries. We set up a framework in which all possible bulk spatial geometries, characterized by weighted adjacent matrices of all possible graphs, are mapped to the boundary Hilbert space and form an overcomplete basis of the boundary. We name such an overcomplete basis as holographic coherent states. A generic boundary state can be expanded on this basis, which describes the state as a superposition of different spatial geometries in the bulk. We discuss how to define distinct classical geometries and small fluctuations around them. We show that small fluctuations around classical geometries define "code subspaces" which are mapped to the boundary Hilbert space isometrically with quantum error correction properties. In addition, we also show that the overlap between different geometries is suppressed exponentially as a function of the geometrical difference between the two geometries. The geometrical difference is measured in an area law fashion, which is a manifestation of the holographic nature of the states considered.


title: KIC 7385478: An eclipsing binary with a {\gamma} Doradus component

abstract: We present spectroscopic and photometric analysis of the eclipsing binary KIC\,7385478. We find that the system is formed by F1V + K4III-IV components. Combining results from analysis of spectroscopic data and $Kepler$ photometry, we calculate masses and radii of the primary and the secondary components as M$_{1}$ = 1.71 $\pm$ 0.08 \Msun, M$_{2}$ = 0.37 $\pm$ 0.04 \Msun~ and R$_{1}$ = 1.59 $\pm$ 0.03\Rsun, R$_{2}$ = 1.90 $\pm$ 0.03\Rsun, respectively. Position of the primary component in HR diagram is in the region of $\gamma$ Doradus type pulsators and residuals from light curve modeling exhibit additional light variation with a dominant period of $\sim$0.5 day. These are clear evidences of the $\gamma$ Doradus type pulsations on the primary component. We also observe occasional increase in amplitude of the residuals, where the orbital period becomes the most dominant period. These may be attributed to the cool star activity originating from the secondary component.


title: Light Bending in the Infinite Derivative Theories of Gravity

abstract: Light bending is one of the significant predictions of general relativity (GR) and it has been confirmed with great accuracy during the last one hundred years. In this paper, we semiclassically calculate the deflection angle for the photons that just grazing the Sun in the infinite derivative theories of gravity (IDG) which is a ghost and singularity free theory of gravity. From our calculations, we find that the deflection angle $\theta$ only depends on $\Lambda/E$. $\theta \rightarrow \theta_E$ when $\Lambda/E \rightarrow \infty$ and decrease to zero when $\Lambda/E \rightarrow 0$. The transition interval occurs at $10^{4}< E/\Lambda < 10^{7}$. It should be pointed out that this model can be tested by the Chandra X-ray Observatory if $0.01 eV < \Lambda < 0.1 eV$.


title: The Relationship Between Agnostic Selective Classification Active Learning and the Disagreement Coefficient

abstract: A selective classifier (f,g) comprises a classification function f and a binary selection function g, which determines if the classifier abstains from prediction, or uses f to predict. The classifier is called pointwise-competitive if it classifies each point identically to the best classifier in hindsight (from the same class), whenever it does not abstain. The quality of such a classifier is quantified by its rejection mass, defined to be the probability mass of the points it rejects. A "fast" rejection rate is achieved if the rejection mass is bounded from above by O(1/m) where m is the number of labeled examples used to train the classifier (and O hides logarithmic factors). Pointwise-competitive selective (PCS) classifiers are intimately related to disagreement-based active learning and it is known that in the realizable case, a fast rejection rate of a known PCS algorithm (called Consistent Selective Strategy) is equivalent to an exponential speedup of the well-known CAL active algorithm. We focus on the agnostic setting, for which there is a known algorithm called LESS that learns a PCS classifier and achieves a fast rejection rate (depending on Hanneke's disagreement coefficient) under strong assumptions. We present an improved PCS learning algorithm called ILESS for which we show a fast rate (depending on Hanneke's disagreement coefficient) without any assumptions. Our rejection bound smoothly interpolates the realizable and agnostic settings. The main result of this paper is an equivalence between the following three entities: (i) the existence of a fast rejection rate for any PCS learning algorithm (such as ILESS); (ii) a poly-logarithmic bound for Hanneke's disagreement coefficient; and (iii) an exponential speedup for a new disagreement-based active learner called ActiveiLESS.


title: A Controlled Set-Up Experiment to Establish Personalized Baselines for Real-Life Emotion Recognition

abstract: We design, conduct and present the results of a highly personalized baseline emotion recognition experiment, which aims to set reliable ground-truth estimates for the subject's emotional state for real-life prediction under similar conditions using a small number of physiological sensors. We also propose an adaptive stimuli-selection mechanism that would use the user's feedback as guide for future stimuli selection in the controlled-setup experiment and generate optimal ground-truth personalized sessions systematically. Initial results are very promising (85% accuracy) and variable importance analysis shows that only a few features, which are easy-to-implement in portable devices, would suffice to predict the subject's emotional state.


title: Scalable Content Delivery with Coded Caching in Multi-Antenna Fading Channels

abstract: We consider the content delivery problem in a fading multi-input single-output channel with cache-aided users. We are interested in the scalability of the equivalent content delivery rate when the number of users, $K$, is large. Analytical results show that, using coded caching and wireless multicasting, without channel state information at the transmitter (CSIT), linear scaling of the content delivery rate with respect to $K$ can be achieved in some different ways. First, if the multicast transmission spans over $L$ independent sub-channels, e.g., in quasi-static fading if $L = 1$, and in block fading or multi-carrier systems if $L>1$, linear scaling can be obtained when the product of the number of transmit antennas and the number of sub-channels scales logarithmically with $K$. Second, even with a fixed number of antennas, we can achieve the linear scaling with a threshold-based user selection requiring only one-bit feedbacks from the users. When CSIT is available, we propose a mixed strategy that combines spatial multiplexing and multicasting. Numerical results show that, by optimizing the power split between spatial multiplexing and multicasting, we can achieve a significant gain of the content delivery rate with moderate cache size.


title: Voltage Tunable Plasmon Propagation in Dual Gated Bilayer Graphene

abstract: In this paper, we theoretically investigate plasmon propagation characteristics in AB and AA stacked bilayer graphene (BLG) in the presence of energy asymmetry due to an electrostatic field oriented perpendicularly to the plane of the graphene sheet. We first derive the optical conductivity of BLG using the Kubo formalism incorporating energy asymmetry and finite electron scattering. All results are obtained for room temperature (300K) operation. By solving Maxwell's equations in a dual gate device setup, we obtain the wavevector of propagating plasmon modes in the transverse electric (TE) and transverse magnetic (TM) directions at terahertz frequencies. The plasmon wavevector allows us to compare the compression factor, propagation length, and the mode confinement of TE and TM plasmon modes in bilayer and monolayer graphene sheets and also study the impact of material parameters on plasmon characteristics. Our results show that the energy asymmetry can be harnessed to increase the propagation length of TM plasmons in BLG. AA stacked BLG shows a larger increase in propagation length than AB stacked BLG; conversely, it is very insensitive to the Fermi level variations. Additionally, the dual gate structure allows independent modulation of the energy asymmetry and the Fermi level in BLG, which is advantageous for reconfiguring plasmon characteristics post device fabrication.


title: Nonexistence of Efficient Dominating Sets in the Cayley Graphs Generated by Transposition Trees of Diameter 3

abstract: Let $d,n$ be positive integers such that $d<n$, and let $X^d_n$ be a Cayley graph generated by a transposition tree of diameter $d$. It is known that every $X^d_n$ with $d<3$ splits into efficient dominating sets. The main result of this paper is that $X^3_n$ does not have efficient dominating sets.


title: Native Language Identification using Stacked Generalization

abstract: Ensemble methods using multiple classifiers have proven to be the most successful approach for the task of Native Language Identification (NLI), achieving the current state of the art. However, a systematic examination of ensemble methods for NLI has yet to be conducted. Additionally, deeper ensemble architectures such as classifier stacking have not been closely evaluated. We present a set of experiments using three ensemble-based models, testing each with multiple configurations and algorithms. This includes a rigorous application of meta-classification models for NLI, achieving state-of-the-art results on three datasets from different languages. We also present the first use of statistical significance testing for comparing NLI systems, showing that our results are significantly better than the previous state of the art. We make available a collection of test set predictions to facilitate future statistical tests.


title: Constructing unextendible product bases from the old ones

abstract: We studied the construction problem of the unextendible product basis (UPB). We mainly give a method to construct a UPB of a quantum system through the UPBs of its subsystem. Using this method and the UPBs which are known for us, we construct different kinds of UPBs in general bipartite quantum system. Then we use these UPBs to construct a family of UPBs in multipartite quantum system. The UPBs can be used to construct the bound entangled states with different ranks.


title: Autoequivalences of tensor categories attached to quantum groups at roots of $1$

abstract: We compute the group of braided tensor autoequivalences and the Brauer-Picard group of the representation category of the small quantum group $\mathfrak{u}_q(\mathfrak{g})$, where $q$ is a root of unity.


title: Squark-mediated Higgs+jets production at the LHC

abstract: We investigate possible scenarios of light-squark production at the LHC as a new mechanism to produce Higgs bosons in association with jets. The study is motivated by the SUSY search for H+jets events, performed by the CMS collaboration on 8 and 13 TeV data using the razor variables. Two simplified models are proposed to interpret the observations in this search. The constraint from Run I and the implications for Run II and beyond are discussed.


title: Hardening Stratum, the Bitcoin Pool Mining Protocol

abstract: Stratum, the de-facto mining communication protocol used by blockchain based cryptocurrency systems, enables miners to reliably and efficiently fetch jobs from mining pool servers. In this paper we exploit Stratum's lack of encryption to develop passive and active attacks on Bitcoin's mining protocol, with important implications on the privacy, security and even safety of mining equipment owners. We introduce StraTap and ISP Log attacks, that infer miner earnings if given access to miner communications, or even their logs. We develop BiteCoin, an active attack that hijacks shares submitted by miners, and their associated payouts. We build BiteCoin on WireGhost, a tool we developed to hijack and surreptitiously maintain Stratum connections. Our attacks reveal that securing Stratum through pervasive encryption is not only undesirable (due to large overheads), but also ineffective: an adversary can predict miner earnings even when given access to only packet timestamps. Instead, we devise Bedrock, a minimalistic Stratum extension that protects the privacy and security of mining participants. We introduce and leverage the mining cookie concept, a secret that each miner shares with the pool and includes in its puzzle computations, and that prevents attackers from reconstructing or hijacking the puzzles. We have implemented our attacks and collected 138MB of Stratum protocol traffic from mining equipment in the US and Venezuela. We show that Bedrock is resilient to active attacks even when an adversary breaks the crypto constructs it uses. Bedrock imposes a daily overhead of 12.03s on a single pool server that handles mining traffic from 16,000 miners.


title: Partial actions of C*-quantum groups I: Restriction and Globalization

abstract: Partial actions of groups on C*-algebras and the closely related actions and coactions of Hopf algebras received much attention over the last decades. They arise naturally as restrictions of their global counterparts to non-invariant subalgebras, and the ambient eveloping global (co)actions have proven useful for the study of associated crossed products. In this article, we introduce the partial coactions of C*-bialgebras, focussing on C*-quantum, and prove existence of an enveloping global coaction under mild technical assumptions. The construction of the latter provides a left adjoint to the forgetful functor from coactions to partial coactions. We also show that partial coactions of the function algebra of a discrete group correspond to partial actions on direct summands of a C*-algebra, and relate partial coactions of a compact or its dual discrete C*-quantum group to partial coactions or partial actions of the dense Hopf subalgebra.


title: Truth-Telling Mechanism for Secure Two-Way Relay Communications with Energy-Harvesting Revenue

abstract: This paper brings the novel idea of paying the utility to the winning agents in terms of some physical entity in cooperative communications. Our setting is a secret two-way communication channel where two transmitters exchange information in the presence of an eavesdropper. The relays are selected from a set of interested parties such that the secrecy sum rate is maximized. In return, the selected relay nodes' energy harvesting requirements will be fulfilled up to a certain threshold through their own payoff so that they have the natural incentive to be selected and involved in the communication. However, relays may exaggerate their private information in order to improve their chance to be selected. Our objective is to develop a mechanism for relay selection that enforces them to reveal the truth since otherwise they may be penalized. We also propose a joint cooperative relay beamforming and transmit power optimization scheme based on an alternating optimization approach. Note that the problem is highly non-convex since the objective function appears as a product of three correlated Rayleigh quotients. While a common practice in the existing literature is to optimize the relay beamforming vector for given transmit power via rank relaxation, we propose a second-order cone programming (SOCP)-based approach in this paper which requires a significantly lower computational task. The performance of the incentive control mechanism and the optimization algorithm has been evaluated through numerical simulations.


title: Riesz continuity of the Atiyah-Singer Dirac operator under perturbations of local boundary conditions

abstract: On a smooth complete Riemannian spin manifold with smooth compact boundary, we demonstrate that the Atiyah-Singer Dirac operator $\mathrm{D}_{\mathcal B}$ in $\mathrm{L}^{2}$ depends Riesz continuously on $\mathrm{L}^{\infty}$ perturbations of local boundary conditions ${\mathcal B}$. The Lipschitz bound for the map ${\mathcal B} \to {\mathrm{D}}_{\mathcal B}(1 + {\mathrm{D}}_{\mathcal B}^2)^{-\frac{1}{2}}$ depends on Lipschitz smoothness and ellipticity of ${\mathcal B}$ and bounds on Ricci curvature and its first derivatives as well as a lower bound on injectivity radius. More generally, we prove perturbation estimates for functional calculi of elliptic operators on manifolds with local boundary conditions.


title: On Helmholtz free energy for finite abstract simplicial complexes

abstract: We prove a Gauss-Bonnet formula X(G) = sum_x K(x), where K(x)=(-1)^dim(x) (1-X(S(x))) is a curvature of a vertex x with unit sphere S(x) in the Barycentric refinement G1 of a simplicial complex G. K(x) is dual to (-1)^dim(x) for which Gauss-Bonnet is the definition of Euler characteristic X. Because the connection Laplacian L'=1+A' of G is unimodular, where A' is the adjacency matrix of of the connection graph G', the Green function values g(x,y) = (1+A')^-1_xy are integers and 1-X(S(x))=g(x,x). Gauss-Bonnet for K^+ reads therefore as str(g)=X(G), where str is the super trace. As g is a time-discrete heat kernel, this is a cousin to McKean-Singer str(exp(-Lt)) = X(G) for the Hodge Laplacian L=dd^* +d^*d which lives on the same Hilbert space than L'. Both formulas hold for an arbitrary finite abstract simplicial complex G. Writing V_x(y)= g(x,y) for the Newtonian potential of the connection Laplacian, we prove sum_y V_x(y) = K(x), so that by the new Gauss-Bonnet formula, the Euler characteristic of G agrees with the total potential theoretic energy sum_x,y g(x,y)=X(G) of G. The curvature K now relates to the probability measure p minimizing the internal energy U(p)=sum_x,y g(x,y) p(x) p(y) of the complex. Since both the internal energy (here linked to topology) and Shannon entropy are natural and unique in classes of functionals, we then look at critical points p the Helmholtz free energy F(p)=(1-T) U(p)-T S(p) which combines the energy functional U and the entropy functional S(p)=-sum_x p(x) log(p(x)). As the temperature T changes, we observe bifurcation phenomena. Already for G=K_3 both a saddle node bifurcation and a pitchfork bifurcation occurs. The saddle node bifurcation leads to a catastrophe: the function T -> F(p(T),T) is discontinuous if p(T) is a free energy minimizer.


title: Iwasawa Invariants of Some Non-Cyclotomic $\mathbb Z_p$-extensions

abstract: Iwasawa showed that there are non-cyclotomic $\mathbb Z_p$-extensions with positive $\mu$-invariant. We show that these $\mu$-invariants can be evaluated explicitly in many situations when $p=2$ and $p=3$.


title: The Molecular Gas Environment in the 20 km s$^{-1}$ Cloud in the Central Molecular Zone

abstract: We recently reported a population of protostellar candidates in the 20 km s$^{-1}$ cloud in the Central Molecular Zone of the Milky Way, traced by H$_2$O masers in gravitationally bound dense cores. In this paper, we report high-angular-resolution ($\sim$3'') molecular line studies of the environment of star formation in this cloud. Maps of various molecular line transitions as well as the continuum at 1.3 mm are obtained using the Submillimeter Array. Five NH$_3$ inversion lines and the 1.3 cm continuum are observed with the Karl G. Jansky Very Large Array. The interferometric observations are complemented with single-dish data. We find that the CH$_3$OH, SO, and HNCO lines, which are usually shock tracers, are better correlated spatially with the compact dust emission from dense cores among the detected lines. These lines also show enhancement in intensities with respect to SiO intensities toward the compact dust emission, suggesting the presence of slow shocks or hot cores in these regions. We find gas temperatures of $\gtrsim$100 K at 0.1-pc scales based on RADEX modelling of the H$_2$CO and NH$_3$ lines. Although no strong correlations between temperatures and linewidths/H$_2$O maser luminosities are found, in high-angular-resolution maps we notice several candidate shock heated regions offset from any dense cores, as well as signatures of localized heating by protostars in several dense cores. Our findings suggest that at 0.1-pc scales in this cloud star formation and strong turbulence may together affect the chemistry and temperature of the molecular gas.


title: Gravitational waves from a supercooled electroweak phase transition and their detection with pulsar timing arrays

abstract: We investigate the properties of a stochastic gravitational wave background produced by a first-order electroweak phase transition in the regime of extreme supercooling. We study a scenario whereby the percolation temperature that signifies the completion of the transition, $T_p$, can be as low as a few MeV (nucleosynthesis temperature), while most of the true vacuum bubbles are formed much earlier at the nucleation temperature, $T_n\sim 50$ GeV. This implies that the gravitational wave spectrum is mainly produced by the collisions of large bubbles and characterised by a large amplitude and a peak frequency as low as $f \sim 10^{-9}-10^{-7}$ Hz. We show that such a scenario can occur in (but not limited to) a model based on a non-linear realisation of the electroweak gauge group, such that the Higgs vacuum configuration is altered by a cubic coupling. In order to carefully quantify the evolution of the phase transition of this model over such a wide temperature range, we go beyond the usual fast transition approximation, taking into account the expansion of the Universe as well as the behaviour of the nucleation probability at low temperatures. Our computation shows that there exists a range of parameters for which the gravitational wave spectrum lies at the edge between the exclusion limits of current pulsar timing array experiments and the detection band of the future Square Kilometre Array observatory.


title: An entrainment model for fully-developed wind farms: effects of atmospheric stability and an ideal limit for wind farm performance

abstract: While a theoretical limit has long been established for the performance of a single turbine, no corresponding upper bound exists for the power output from a large wind farm, making it difficult to evaluate the available potential for further performance gains. Here we build a model describing the essential features of a large array of turbines with arbitrary design and layout, by considering a fully-developed wind farm whose upper edge is bounded by a self-similar boundary layer. The exchanges between the wind farm, the overlaying boundary layer, and the outer flow are parameterized by means of the classical entrainment hypothesis. We obtain a concise expression for the wind farm's power density (corresponding to power output per unit planform area), as a function of three coefficients, which represent the array thrust and the turbulent exchanges at each of the two interfaces. We assess the performance of our model by comparing the predicted power density to field data, laboratory measurements and large-eddy simulations for the fully-developed regions of wind farms, finding good agreement. We extend our model to include the effect of atmospheric stability on power output, by using an established parameterization relating entrainment coefficients to local Froude numbers. Our predictions for power variation with atmospheric stability are in broad agreement with field measurements. To the best of our knowledge, this constitutes the first quantitative comparison between an atmospheric-stability-dependent model and field data. Finally, we consider an ideal limit for array operation, whereby turbines are designed to maximize momentum exchange with the overlying boundary layer. This enables us to obtain an upper bound for the performance of large wind farms, which we determine to be an order of magnitude larger than the output of contemporary turbine arrays.


title: Object category understanding via eye fixations on freehand sketches

abstract: The study of eye gaze fixations on photographic images is an active research area. In contrast, the image subcategory of freehand sketches has not received as much attention for such studies. In this paper, we analyze the results of a free-viewing gaze fixation study conducted on 3904 freehand sketches distributed across 160 object categories. Our analysis shows that fixation sequences exhibit marked consistency within a sketch, across sketches of a category and even across suitably grouped sets of categories. This multi-level consistency is remarkable given the variability in depiction and extreme image content sparsity that characterizes hand-drawn object sketches. In our paper, we show that the multi-level consistency in the fixation data can be exploited to (a) predict a test sketch's category given only its fixation sequence and (b) build a computational model which predicts part-labels underlying fixations on objects. We hope that our findings motivate the community to deem sketch-like representations worthy of gaze-based studies vis-a-vis photographic images.


title: Substrate effects on the exciton fine structure of black phosphorus quantum dots

abstract: We study the size-dependent exciton fine structure in monolayer black phosphorus quantum dots (BPQDs) deposited on different substrates (isolated, Si and SiO$_2$) using a combination of tight-binding method to calculate the single-particle states, and the configuration interaction formalism to determine the excitonic spectrum. We demonstrate that the substrate plays a dramatic role on the excitonic gaps and excitonic spectrum of the QDs. For reasonably high dielectric constants ($\varepsilon_{sub} \sim \varepsilon_{Si} = 11.7 \varepsilon_0$), the excitonic gap can be described by a single power law $E_X(R) = E_X^{(bulk)} + C/R^{\gamma}$. For low dielectric constants $\varepsilon_{sub} \leq \varepsilon_{SiO_2} = 3.9 \varepsilon_0$, the size dependence of the excitonic gaps requires the sum of two power laws $E_X(R) = E_g^{(bulk)} + A/ R^{n} - B/R^{m}$ to describe both strong and weak quantum confinement regimes, where $A$, $B$, $C$, $\gamma$, $n$, and $m$ are substrate-dependent parameters. We also predict that the exciton lifetimes exhibit a strong temperature dependence, ranging between 2-8 ns (Si substrate) and 3-11 ns (SiO$_2$ substrate) for QDs up 10 nm in size.


title: Connected Dominating Sets in Graphs With Stability Number Three

abstract: In the special case of graphs G of independence number a(G)=3 without induced chordless cycles C7 it is shown that exists connected dominating set D of vertices with number of vertices n(D)<=4. Using the concept of connected dominating sets, we defined a new invariant h(G) that does not exceed the number of Hadwiger. For the considered graphs it is shown that h(G)>=n(G)/4.


title: Strain Mode of General Flow: Characterization and Implications for Flow Pattern Structures

abstract: Understanding the mixing capability of mixing devices based on their geometric shape is an important issue both for predicting mixing processes and for designing new mixers. The flow patterns in mixers are directly connected with the modes of the local strain rate, which is generally a combination of elongational flow and planar shear flow. We develop a measure to characterize the modes of the strain rate for general flow occurring in mixers. The spatial distribution of the volumetric strain rate (or non-planar strain rate) in connection with the flow pattern plays an essential role in understanding distributive mixing. With our measure, flows with different types of screw elements in a twin-screw extruder are numerically analyzed. The difference in flow pattern structure between conveying screws and kneading disks is successfully characterized by the distribution of the volumetric strain rate. The results suggest that the distribution of the strain rate mode offers an essential and convenient way for characterization of the relation between flow pattern structure and the mixer geometry.


title: Using Maximum Entry-Wise Deviation to Test the Goodness-of-Fit for Stochastic Block Models

abstract: The stochastic block model is widely used for detecting community structures in network data. How to test the goodness-of-fit of the model is one of the fundamental problems and has gained growing interests in recent years. In this article, we propose a novel goodness-of-fit test based on the maximum entry of the centered and re-scaled adjacency matrix for the stochastic block model. One noticeable advantage of the proposed test is that the number of communities can be allowed to grow linearly with the number of nodes ignoring a logarithmic factor. We prove that the null distribution of the test statistic converges in distribution to a Gumbel distribution, and we show that both the number of communities and the membership vector can be tested via the proposed method. Further, we show that the proposed test has asymptotic power guarantee against a class of alternatives. We also demonstrate that the proposed method can be extended to the degree-corrected stochastic block model. Both simulation studies and real-world data examples indicate that the proposed method works well.