title: Mahler's Work and Algebraic Dynamical Systems

abstract: After Furstenberg had provided a first glimpse of remarkable rigidity phenomena associated with the joint action of several commuting automorphisms (or endomorphisms) of a compact abelian group, further key examples motivated the development of an extensive theory of such actions. Two of Mahler's achievements, the recognition of the significance of Mahler measure of multivariate polynomials in relating the lengths and heights of products of polynomials in terms of the corresponding quantities for the constituent factors, and his work on additive relations in fields, have unexpectedly played important roles in the study of entropy and higher order mixing for these actions. This article briefly surveys these connections between Mahler's work and dynamics. It also sketches some of the dynamical outgrowths of his work that are very active today, including the investigation of the Fuglede-Kadison determinant of a convolution operator in a group von Neumann algebra as a noncommutative generalization of Mahler measure, as well as diophantine questions related to the growth rates of periodic points and their relation to entropy.


title: Learned Neural Iterative Decoding for Lossy Image Compression Systems

abstract: For lossy image compression systems, we develop an algorithm, iterative refinement, to improve the decoder's reconstruction compared to standard decoding techniques. Specifically, we propose a recurrent neural network approach for nonlinear, iterative decoding. Our decoder, which works with any encoder, employs self-connected memory units that make use of causal and non-causal spatial context information to progressively reduce reconstruction error over a fixed number of steps. We experiment with variants of our estimator and find that iterative refinement consistently creates lower distortion images of higher perceptual quality compared to other approaches. Specifically, on the Kodak Lossless True Color Image Suite, we observe as much as a 0.871 decibel (dB) gain over JPEG, a 1.095 dB gain over JPEG 2000, and a 0.971 dB gain over a competitive neural model.


title: The limited role of recombination energy in common envelope removal

abstract: We calculate the outward energy transport time by convection and photon diffusion in an inflated common envelope and find this time to be shorter than the envelope expansion time. We conclude therefore that most of the hydrogen recombination energy ends in radiation rather than in kinetic energy of the outflowing envelope. We use the stellar evolution code MESA and inject energy inside the envelope of an asymptotic giant branch star to mimic energy deposition by a spiraling-in stellar companion. During 1.7 years the envelope expands by a factor of more than 2. Along the entire evolution the convection can carry the energy very efficiently outwards, to the radius where radiative transfer becomes more efficient. The total energy transport time stays within several months, shorter than the dynamical time of the envelope. Had we included rapid mass loss, as is expected in the common envelope evolution, the energy transport time would have been even shorter. It seems that calculations that assume that most of the recombination energy ends in the outflowing gas might be inaccurate.


title: Kinetics of electron cooling in metal films at low temperatures and revision of the two-temperature model

abstract: The two-temperature model (2TM) introduced by Kaganov, Lifshitz, and Tanatarov is widely used to describe the energy relaxation in the electron-phonon system of a metallic film. At the same time, the accuracy of the description of the electron-phonon system in terms of 2TM, i.e. on the basis of the electron and phonon temperatures, has not been considered in detail until now. In this paper we present a microscopic theory of cooling of instantly heated electrons in metallic films. In framework of this theory the main features of electron cooling in thick and thin films were found, and an analysis of the accuracy of the 2TM in the low-temperature region was carried out. We consider a more accurate three-temperature model, which (in contrast to 2TM) explicitly takes into account phonons with angles of incidence exceeding the angle of total internal reflection. The contribution of these phonons to the kinetics of electron relaxation can be significant if the sound velocities in the film and the substrate are quite different.


title: The complexity of comparing multiply-labelled trees by extending phylogenetic-tree metrics

abstract: A multilabeled tree (or MUL-tree) is a rooted tree in which every leaf is labelled by an element from some set, but in which more than one leaf may be labelled by the same element of that set. In phylogenetics, such trees are used in biogeographical studies, to study the evolution of gene families, and also within approaches to construct phylogenetic networks. A multilabelled tree in which no leaf-labels are repeated is called a phylogenetic tree, and one in which every label is the same is also known as a tree-shape. In this paper, we consider the complexity of computing metrics on MUL-trees that are obtained by extending metrics on phylogenetic trees. In particular, by restricting our attention to tree shapes, we show that computing the metric extension on MUL-trees is NP complete for two well-known metrics on phylogenetic trees, namely, the path-difference and Robinson Foulds distances. We also show that the extension of the Robinson Foulds distance is fixed parameter tractable with respect to the distance parameter. The path distance complexity result allows us to also answer an open problem concerning the complexity of solving the quadratic assignment problem for two matrices that are a Robinson similarity and a Robinson dissimilarity, which we show to be NP-complete. We conclude by considering the maximum agreement subtree (MAST) distance on phylogenetic trees to MUL-trees. Although its extension to MUL-trees can be computed in polynomial time, we show that computing its natural generalization to more than two MUL-trees is NP-complete, although fixed-parameter tractable in the maximum degree when the number of given trees is bounded.


title: Capturing Structure Implicitly from Time-Series having Limited Data

abstract: Scientific fields such as insider-threat detection and highway-safety planning often lack sufficient amounts of time-series data to estimate statistical models for the purpose of scientific discovery. Moreover, the available limited data are quite noisy. This presents a major challenge when estimating time-series models that are robust to overfitting and have well-calibrated uncertainty estimates. Most of the current literature in these fields involve visualizing the time-series for noticeable structure and hard coding them into pre-specified parametric functions. This approach is associated with two limitations. First, given that such trends may not be easily noticeable in small data, it is difficult to explicitly incorporate expressive structure into the models during formulation. Second, it is difficult to know $\textit{a priori}$ the most appropriate functional form to use. To address these limitations, a nonparametric Bayesian approach was proposed to implicitly capture hidden structure from time series having limited data. The proposed model, a Gaussian process with a spectral mixture kernel, precludes the need to pre-specify a functional form and hard code trends, is robust to overfitting and has well-calibrated uncertainty estimates.


title: Free subgroups of $3$-manifold groups

abstract: We show that any closed hyperbolic $3$-manifold $M$ has a co-final tower of covers $M_i \to M$ of degrees $n_i$ such that any subgroup of $\pi_1(M_i)$ generated by $k_i$ elements is free, where $k_i \ge n_i^C$ and $C = C(M) > 0$. Together with this result we show that $\log k_i \geq C_1 sys_1(M_i)$, where $sys_1(M_i)$ denotes the systole of $M_i$, thus providing a large set of new examples for a conjecture of Gromov. In the second theorem $C_1> 0$ is an absolute constant. We also consider a generalization of these results to non-compact finite volume hyperbolic $3$-manifolds.


title: Magneto-elastic switching of magnetostrictive nanomagnets with in-plane shape anisotropy

abstract: We theoretically study the effect of a material defect (material void) on switching errors associated with magneto-elastic switching of magnetization in elliptical magnetostrictive nanomagnets having in-plane magnetic anisotropy. We find that the error probability increases significantly in the presence of the defect, indicating that magneto-elastic switching is particularly vulnerable to material imperfections. Curiously, there is a critical stress value that gives the lowest error probability in both defect-free and defective nanomagnets. The critical stress is much higher in defective nanomagnets than in defect-free ones. Since it is more difficult to generate the critical stress in small nanomagnets than in large nanomagnets (having the same energy barrier for thermal stability), it would be a challenge to downscale magneto-elastically switched nanomagnets in memory and other applications where reliable switching is required. This is likely to be further exacerbated by the presence of defects.


title: pH Sensing by Lipids in Membranes: The Fundamentals of pH-driven Migration, Polarization and Deformations of Lipid Bilayer Assemblies

abstract: Most biological molecules contain acido-basic groups that modulate their structure and interactions. A consequence is that pH gradients, local heterogeneities and dynamic variations are used by cells and organisms to drive or regulate specific biological functions including energetic metabolism, vesicular traffic, migration and spatial patterning of tissues in development. While the direct or regulatory role of pH in protein function is well documented, the role of hydrogen and hydroxyl ions in modulating the properties of lipid assemblies such as bilayer membranes is only beginning to be understood. Here, we review approaches using artificial lipid vesicles that have been instrumental in providing an understanding of the influence of pH gradients and local variations on membrane vectorial motional processes: migration, membrane curvature effects promoting global or local deformations, crowding generation by segregative polarization processes. In the case of pH induced local deformations, an extensive theoretical framework is given and an application to a specific biological issue, namely the structure and stability of mitochondrial cristae, is described.


title: Distributed Data Vending on Blockchain

abstract: Recent advances in blockchain technologies have provided exciting opportunities for decentralized applications. Specifically, blockchain-based smart contracts enable credible transactions without authorized third parties. The attractive properties of smart contracts facilitate distributed data vending, allowing for proprietary data to be securely exchanged on a blockchain. Distributed data vending can transform domains such as healthcare by encouraging data distribution from owners and enabling large-scale data aggregation. However, one key challenge in distributed data vending is the trade-off dilemma between the effectiveness of data retrieval, and the leakage risk from indexing the data. In this paper, we propose a framework for distributed data vending through a combination of data embedding and similarity learning. We illustrate our framework through a practical scenario of distributing and aggregating electronic medical records on a blockchain. Extensive empirical results demonstrate the effectiveness of our framework.


title: Virtual CNN Branching: Efficient Feature Ensemble for Person Re-Identification

abstract: In this paper we introduce an ensemble method for convolutional neural network (CNN), called "virtual branching," which can be implemented with nearly no additional parameters and computation on top of standard CNNs. We propose our method in the context of person re-identification (re-ID). Our CNN model consists of shared bottom layers, followed by "virtual" branches, where neurons from a block of regular convolutional and fully-connected layers are partitioned into multiple sets. Each virtual branch is trained with different data to specialize in different aspects, e.g., a specific body region or pose orientation. In this way, robust ensemble representations are obtained against human body misalignment, deformations, or variations in viewing angles, at nearly no any additional cost. The proposed method achieves competitive performance on multiple person re-ID benchmark datasets, including Market-1501, CUHK03, and DukeMTMC-reID.


title: Deep Structure Inference Network for Facial Action Unit Recognition

abstract: Facial expressions are combinations of basic components called Action Units (AU). Recognizing AUs is key for developing general facial expression analysis. In recent years, most efforts in automatic AU recognition have been dedicated to learning combinations of local features and to exploiting correlations between Action Units. In this paper, we propose a deep neural architecture that tackles both problems by combining learned local and global features in its initial stages and replicating a message passing algorithm between classes similar to a graphical model inference approach in later stages. We show that by training the model end-to-end with increased supervision we improve state-of-the-art by 5.3% and 8.2% performance on BP4D and DISFA datasets, respectively.


title: Synthesizing geocodes to facilitate access to detailed geographical information in large scale administrative data

abstract: We investigate whether generating synthetic data can be a viable strategy for providing access to detailed geocoding information for external researchers, without compromising the confidentiality of the units included in the database. Our work was motivated by a recent project at the Institute for Employment Research (IAB) in Germany that linked exact geocodes to the Integrated Employment Biographies (IEB), a large administrative database containing several million records. We evaluate the performance of three synthesizers regarding the trade-off between preserving analytical validity and limiting disclosure risks: One synthesizer employs Dirichlet Process mixtures of products of multinomials (DPMPM), while the other two use different versions of Classification and Regression Trees (CART). In terms of preserving analytical validity, our proposed synthesis strategy for geocodes based on categorical CART models outperforms the other two. If the risks of the synthetic data generated by the categorical CART synthesizer are deemed too high, we demonstrate that synthesizing additional variables is the preferred strategy to address the risk-utility trade-off in practice, compared to limiting the size of the regression trees or relying on the strategy of providing geographical information only on an aggregated level. We also propose strategies for making the synthesizers scalable for large files, present analytical validity measures and disclosure risk measures for the generated data, and provide general recommendations for statistical agencies considering the synthetic data approach for disseminating detailed geographical information.


title: Maxiset point of view for signal detection in inverse problems

abstract: This paper extends the successful maxiset paradigm from function estimation to signal detection in inverse problems. In this context, the maxisets do not have the same shape compared to the classical estimation framework. Nevertheless, we introduce a robustified version of these maxisets, allowing to exhibit tail conditions on the signals of interest. Under this novel paradigm we are able to compare direct and indirect testing procedures.


title: Predictions of Angular Observables for $\bar{B}_s\to K^{\ast}\ell\ell$ and $\bar{B}\to \rho\ell\ell$ in Standard Model

abstract: Exclusive semileptonic decays based on $b\to s$ transitions have been attracting a lot of attention as some angular observables deviate significantly from the Standard Model (SM) predictions in specific $q^2$ bins. B meson decays induced by other Flavor Changing Neutral Current (FCNC), $b\to d$, can also offer a probe to new physics with an additional sensitivity to the weak phase in Cabibbo-Kobayashi-Masakawa (CKM) matrix. We provide predictions for angular observables for $b\to d$ semileptonic transitions, namely $\bar{B}_s\to K^{\ast}\ell^+\ell^-$, $\bar{B}^0\to\rho^0 \ell^+\ell^-$, and their CP-conjugated modes. For $\bar{B}^0\to\rho^0 \ell^+\ell^-$ mode, $B^0-\bar{B}^0$ mixing effects have been included and predictions are made for Belle and LHCb separately. Study of these decay modes will be useful in its own right and to understand the pattern of deviations in $b\to s$ transitions.


title: Orientational alignment in Cavity Quantum Electrodynamics

abstract: We consider the orientational alignment of dipoles due to strong matter light coupling, for a non-vanishing density of excitations. We compare various approaches to this problem in the limit of large numbers of emitters, and show that direct Monte Carlo integration, mean-field theory, and large deviation methods match exactly in this limit. All three results show that orientational alignment develops in the presence of a macroscopically occupied polariton mode, and that the dipoles asymptotically approach perfect alignment in the limit of high density or low temperature.


title: The Laplace transform of the lognormal distribution

abstract: We study the analytical properties of the Laplace transform of the lognormal distribution. Two integral expressions for the analytic continuation of the Laplace transform of the lognormal distribution are provided, one of which takes the form of a Mellin-Barnes integral. As a corollary, we obtain an integral expression for the characteristic function; we show that the integral expression derived by Leipnik in [11] is incorrect. We present two approximations for the Laplace transform of the lognormal distribution, both valid in $\mathbb{C} \setminus(-\infty,0]$. In the last section, we discuss how one may use our results to compute the density of a sum of independent lognormal random variables.


title: Capillary stress and structural relaxation in moist granular materials

abstract: We propose a theoretical framework to calculate capillary stresses in complex mesoporous materials, such as moist sand, nanoporous hydrates, and drying colloidal films. Molecular simulations are mapped onto a phase-field model of the liquid-vapor mixture, whose inhomogeneous stress tensor is integrated over Voronoi polyhedra in order to calculate equal and opposite forces between each pair of neighboring grains. The method is illustrated by simulations of moisture-induced forces in small clusters and random packings of spherical grains using lattice-gas Density Functional Theory. For a nano-granular model of cement hydrates, this approach reproduces the hysteretic water sorption/desorption isotherms and predicts drying shrinkage strain isotherm in good agreement with experiments. We show that capillary stress is an effective mechanism for internal stress relaxation in colloidal random packings, which contributes to the extraordinary durability of cement paste.


title: GossipGraD: Scalable Deep Learning using Gossip Communication based Asynchronous Gradient Descent

abstract: In this paper, we present GossipGraD - a gossip communication protocol based Stochastic Gradient Descent (SGD) algorithm for scaling Deep Learning (DL) algorithms on large-scale systems. The salient features of GossipGraD are: 1) reduction in overall communication complexity from {\Theta}(log(p)) for p compute nodes in well-studied SGD to O(1), 2) model diffusion such that compute nodes exchange their updates (gradients) indirectly after every log(p) steps, 3) rotation of communication partners for facilitating direct diffusion of gradients, 4) asynchronous distributed shuffle of samples during the feedforward phase in SGD to prevent over-fitting, 5) asynchronous communication of gradients for further reducing the communication cost of SGD and GossipGraD. We implement GossipGraD for GPU and CPU clusters and use NVIDIA GPUs (Pascal P100) connected with InfiniBand, and Intel Knights Landing (KNL) connected with Aries network. We evaluate GossipGraD using well-studied dataset ImageNet-1K (~250GB), and widely studied neural network topologies such as GoogLeNet and ResNet50 (current winner of ImageNet Large Scale Visualization Research Challenge (ILSVRC)). Our performance evaluation using both KNL and Pascal GPUs indicates that GossipGraD can achieve perfect efficiency for these datasets and their associated neural network topologies. Specifically, for ResNet50, GossipGraD is able to achieve ~100% compute efficiency using 128 NVIDIA Pascal P100 GPUs - while matching the top-1 classification accuracy published in literature.


title: Soft-Photon Corrections to $\bar{B} \to D \tau^{-} \bar{\nu}_{\tau}$ Relative to $\bar{B} \to D \mu^{-} \bar{\nu}_{\mu}$

abstract: We evaluate long-distance electromagnetic (QED) contributions to $\bar{B}{}^0 \to D^+ \tau^{-} \bar{\nu}_{\tau}$ and $B^- \to D^0 \tau^{-} \bar{\nu}_{\tau}$ relative to $\bar{B}{}^0 \to D^+ \mu^{-} \bar{\nu}_{\mu}$ and $B^- \to D^0 \mu^{-} \bar{\nu}_{\mu}$, respectively, in the standard model. We point out that the QED corrections to the ratios $R(D^{+})$ and $R(D^{0})$ are not negligible, contrary to the expectation that radiative corrections are almost canceled out in the ratio of the two branching fractions. The reason is that long-distance QED corrections depend on the masses and relative velocities of the daughter particles. We find that theoretical predictions for $R(D^{+})^{\tau/\mu}$ and $R(D^{0})^{\tau/\mu}$ can be amplified by $\sim4\%$ and $\sim3\%$, respectively, for the soft-photon energy cut in range $20$-$40$ MeV.


title: Index one minimal surfaces in spherical space forms

abstract: We prove that orientable index one minimal surfaces in spherical space forms with large fundamental group have genus at most two. This confirms a conjecture of R. Schoen for an infinite class of 3-manifolds.


title: Monodromy of elliptic curve convolution, seven-point sheaves of $G_2$-type and motives of Beauville type

abstract: We study the Tannakian properties of the category of perverse sheaves on elliptic curves endowed with the convolution product. We establish that for certain sheaves with unipotent local monodromy over seven points the corresponding Tannaka group is isomorphic to $G_2$. This monodromy approach generalizes a result of Katz on the existence of $G_2$-motives in the middle cohomology of deformations of Beauville surfaces.


title: Lineshape-asymmetry elimination in weak atomic transitions driven by an intense standing wave field

abstract: Owing to the ac-Stark effect, the lineshape of a weak optical transition in an atomic beam can become significantly distorted, when driven by an intense standing wave field. We use an Yb atomic beam to study the lineshape of the 6s2 1S0 -> 5d6s 3D1 transition, which is excited with light circulating in a Fabry-Perot resonator. We demonstrate two methods to avoid the distortion of the transition profile. Of these, one relies on the operation of the resonator in multiple longitudinal modes, and the other in multiple transverse modes.


title: Identifiability of dynamical networks with partial node measurements

abstract: Much recent research has dealt with the identifiability of a dynamical network in which the node signals are connected by causal linear transfer functions and are excited by known external excitation signals and/or unknown noise signals. A major research question concerns the identifiability of the whole network - topology and all transfer functions - from the measured node signals and external excitation signals. So far all results on this topic have assumed that all node signals are measured. This paper presents the first results for the situation where not all node signals are measurable, under the assumptions that (1) the topology of the network is known, and (2) each node is excited by a known external excitation. Using graph theoretical properties, we show that the transfer functions that can be identified depend essentially on the topology of the paths linking the corresponding vertices to the measured nodes. A practical outcome is that, under those assumptions, a network can often be identified using only a small subset of node measurements.


title: Universal quantum computing with parafermions assisted by a half fluxon

abstract: Braiding of anyons such as Majoranas or parafermions provides only Clifford gates which do not form a universal set of quantum gates. We propose a robust and resource-efficient scheme to perform a non-Clifford gate on a logical qudit encoded in parafermionic zero modes via the Aharonov-Casher effect. This gate can be implemented by moving a half flux quantum around the pair of parafermionic zero modes. The parafermion modes can be realized in a two-dimensional set-up using existing proposals and a half fluxon carrying half flux quantum can be created as a part of a half fluxon/anti-half fluxon pair in a spin-triplet Josephson junction with a dipole defect. With an appropriate bias current pulse, the half fluxon can be braided around the parafermions. Supplementing this gate with the braiding of parafermions provides the avenue for universal quantum computing with parafermions without magic state distillation.


title: Coulomb-gas electrostatics controls large fluctuations of the KPZ equation

abstract: We establish a large deviation principle for the Kardar-Parisi-Zhang (KPZ) equation, providing precise control over the left tail of the height distribution for narrow wedge initial condition. Our analysis exploits an exact connection between the KPZ one-point distribution and the Airy point process -- an infinite particle Coulomb-gas which arises at the spectral edge in random matrix theory. We develop the large deviation principle for the Airy point process and use it to compute, in a straight-forward and assumption-free manner, the KPZ large deviation rate function in terms of an electrostatic problem (whose solution we evaluate). This method also applies to the half-space KPZ equation, showing that its rate function is half of the full-space rate function. In addition to these long-time estimates, we provide rigorous proof of finite-time tail bounds on the KPZ distribution which demonstrate a crossover between exponential decay with exponent $3$ (in the shallow left tail) to exponent $5/2$ (in the deep left tail). The full-space KPZ rate function agrees with the one computed in Sasorov et al. [J. Stat. Mech, 063203 (2017)] via a WKB approximation analysis of a non-local, non-linear integro-differential equation generalizing Painlev\'e II which Amir et al. [Comm. Pure Appl. Math. 64, 466 (2011)] related to the KPZ one-point distribution.

categories: cond-mat.stat-mech math-ph math.MP math.PR


title: A solvable problem in statistical mechanics: the dipole-type Hamiltonian mean field model

abstract: The present study regards the zeroth order mean field approximation of a dipole-type interaction model, which is analytically solved in the canonical and microcanonical ensembles. After writing the canonical partition function, the free and internal energies, magnetization and the specific heat are derived and graphically represented. A crucial derivation is the calculation of the free energy, which is variationally evaluated, and it is shown that the exact result coincides with the approximate trend when $N$ tends to infinity. In the microcanonical ensemble, the entropy as other thermodynamic properties are calculated. We notice that both schemes coincide in equilibrium.


title: Using Automatic Refactoring to Improve Energy Efficiency of Android Apps

abstract: The ever-growing popularity of mobile phones has brought additional challenges to the software development lifecycle. Mobile applications (apps, for short) ought to provide the same set of features as conventional software, with limited resources: such as, limited processing capabilities, storage, screen and, not less important, power source. Although energy efficiency is a valuable requirement, developers often lack knowledge of best practices. In this paper, we study whether or not automatic refactoring can aid developers ship energy efficient apps. We leverage a tool, Leafactor, with five energy code smells that tend to go unnoticed. We use Leafactor to analyze code smells in 140 free and open source apps. As a result, we detected and fixed code smells in 45 apps, from which 40% have successfully merged our changes into the official repository.


title: Blow-up results for space-time fractional stochastic partial differential equations

abstract: Consider non-linear time-fractional stochastic reaction-diffusion equations of the following type, $$\partial^\beta_tu_t(x)=-\nu(-\Delta)^{\alpha/2} u_t(x)+I^{1-\beta}_t[b(u)+ \sigma(u)\stackrel{\cdot}{F}(t,x)]$$ in $(d+1)$ dimensions, where $\nu>0, \beta\in (0,1)$, $\alpha\in (0,2]$. The operator $\partial^\beta_t$ is the Caputo fractional derivative while $-(-\Delta)^{\alpha/2} $ is the generator of an isotropic $\alpha$-stable L\'evy process and $I^{1-\beta}_t$ is the Riesz fractional integral operator. The forcing noise denoted by $\stackrel{\cdot}{F}(t,x)$ is a Gaussian noise. These equations might be used as a model for materials with random thermal memory. We derive non-existence (blow-up) of global random field solutions under some additional conditions, most notably on $b$, $\sigma$ and the initial condition. Our results complement those of P. Chow in \cite{chow2}, \cite{chow1}, and Foondun et al. in \cite{Foondun-liu-nane}, \cite{foondun-parshad} among others.


title: Restoring Heisenberg scaling in noisy quantum metrology by monitoring the environment

abstract: We study quantum frequency estimation for $N$ qubits subjected to independent Markovian noise, via strategies based on time-continuous monitoring of the environment. Both physical intuition and an extended convexity property of the quantum Fisher information (QFI) suggest that these strategies are more effective than the standard ones based on the measurement of the unconditional state after the noisy evolution. Here we focus on initial GHZ states and on parallel or transverse noise. For parallel noise, i.e. dephasing, we show that perfectly efficient time-continuous photo-detection allows to recover the unitary (noiseless) QFI, and thus to obtain a Heisenberg scaling for every value of the monitoring time. For finite detection efficiency, one falls back to the noisy standard quantum limit scaling, but with a constant enhancement due to an effective reduced dephasing. Also in the transverse noise case we obtain that the Heisenberg scaling is recovered for perfectly efficient detectors, and we find that both homodyne and photo-detection based strategies are optimal. For finite detectors efficiency, our numerical simulations show that, as expected, an enhancement can be observed, but we cannot give any conclusive statement regarding the scaling. We finally describe in detail the stable and compact numerical algorithm that we have developed in order to evaluate the precision of such time-continuous estimation strategies, and that may find application in other quantum metrology schemes.


title: Afterglow Imaging and Polarization of Misaligned Structured GRB Jets and Cocoons: Breaking the Degeneracy in GRB 170817A

abstract: The X-ray to radio afterglow emission of GRB 170817A/GW 170817 so far scales as $F_\nu\propto\nu^{-0.6}t^{0.8}$ with observed frequency and time, consistent with a single power-law segment of the synchrotron spectrum from the external shock going into the ambient medium. This requires the effective isotropic equivalent afterglow shock energy in the visible region to increase as $\sim t^{1.7}$. The two main channels for such an energy increase are (i) \emph{radial}: more energy carried by slower material (in the visible region) gradually catches up with the afterglow shock and energizes it, and (ii) \emph{angular}: more energy in relativistic outflow moving at different angles to our line of sight, whose radiation is initially beamed away from us but its beaming cone gradually reaches our line of sight as it decelerates. One cannot distinguish between these explanations (or combinations of them) using only the X-ray to radio $F_\nu(t)$. Here we demonstrate that the most promising way to break this degeneracy is through afterglow imaging and polarization, by calculating the predicted evolution of the afterglow image (its size, shape and flux centroid) and linear polarization $\Pi(t)$ for different angular and/or radial outflow structures that fit $F_\nu(t)$. We consider two angular profiles -- a Gaussian and a narrow core with power-law wings in energy per solid angle, as well as a (cocoon motivated) (quasi-) spherical flow with radial velocity profile. For a jet viewed off-axis (and a magnetic field produced in the afterglow shock) $\Pi(t)$ peaks when the jet's core becomes visible, at $\approx2t_p$ where the lightcurve peaks at $t_p$, and the image can be elongated with aspect ratios$\;\gtrsim2$. A quasi-spherical flow has an almost circular image and a much lower $\Pi(t)$ (peaking at $\approx t_p$) and flux centroid displacement.


title: Altitude Terrain Guarding and Guarding Uni-Monotone Polygons

abstract: We present an optimal, linear-time algorithm for the following version of terrain guarding: given a 1.5D terrain and a horizontal line, place the minimum number of guards on the line to see all of the terrain. We prove that the cardinality of the minimum guard set coincides with the cardinality of a maximum number of ``witnesses'', i.e., terrain points, no two of which can be seen by a single guard. We show that our results also apply to the Art Gallery problem in ``monotone mountains'', i.e., $x$-monotone polygons with a single edge as one of the boundary chains. This means that any monotone mountain is ``perfect'' (its guarding number is the same as its witness number); we thus establish the first non-trivial class of perfect polygons.


title: The diagrammatic coaction and the algebraic structure of cut Feynman integrals

abstract: We present a new formula for the coaction of a large class of integrals. When applied to one-loop (cut) Feynman integrals, it can be given a diagrammatic representation purely in terms of pinches and cuts of the edges of the graph. The coaction encodes the algebraic structure of these integrals, and offers ways to extract important properties of complicated integrals from simpler functions. In particular, it gives direct access to discontinuities of Feynman integrals and facilitates a straightforward derivation of the differential equations they satisfy, which we illustrate in the case of the pentagon.


title: Weak Modular Zilber-Pink with Derivatives

abstract: In unpublished notes Pila proposed a Modular Zilber-Pink with Derivatives (MZPD) conjecture, which is a Zilber-Pink type statement for the modular $j$-function and its derivatives. In this article we define D-special varieties, then state and prove two functional (differential) analogues of the MZPD conjecture for those varieties. In particular, we prove a weak version of MZPD. As a special case of our results, we obtain a functional Modular Andr\'e-Oort with Derivatives statement. The main tools used in the paper come from (model theoretic) differential algebra and complex analytic geometry, and the Ax-Schanuel theorem for the $j$-function and its derivatives (established by Pila and Tsimerman) plays a crucial role in our proofs. In the proof of the second Zilber-Pink type theorem we also use an Existential Closedness statement for the differential equation of the $j$-function.


title: Length in the Cremona group

abstract: The Cremona group is the group of birational transformations of the plane. A birational transformation for which there exists a pencil of lines which is sent onto another pencil of lines is called a Jonqui\`eres transformation. By the famous Noether-Castelnuovo theorem, every birational transformation $f$ is a product of Jonqui\`eres transformations. The minimal number of factors in such a product will be called the length, and written $\mathrm{lgth}(f)$. Even if this length is rather unpredictable, we provide an explicit algorithm to compute it, which only depends on the multiplicities of the linear system of $f$. As an application of this computation, we give a few properties of the dynamical length of $f$ defined as the limit of the sequence $n \mapsto \mathrm{lgth} (f^n) / n$. It follows for example that an element of the Cremona group is distorted if and only if it is algebraic. The computation of the length may also be applied to the so called Wright complex associated with the Cremona group: This has been done recently by Lonjou. Moreover, we show that the restriction of the length to the automorphism group of the affine plane is the classical length of this latter group (the length coming from its amalgamated structure). In another direction, we compute the lengths and dynamical lengths of all monomial transformations, and of some Halphen transformations. Finally, we show that the length is a lower semicontinuous map on the Cremona group endowed with its Zariski topology.


title: Contrasting information theoretic decompositions of modulatory and arithmetic interactions in neural information processing systems

abstract: Biological and artificial neural systems are composed of many local processors, and their capabilities depend upon the transfer function that relates each local processor's outputs to its inputs. This paper uses a recent advance in the foundations of information theory to study the properties of local processors that use contextual input to amplify or attenuate transmission of information about their driving inputs. This advance enables the information transmitted by processors with two distinct inputs to be decomposed into those components unique to each input, that shared between the two inputs, and that which depends on both though it is in neither, i.e. synergy. The decompositions that we report here show that contextual modulation has information processing properties that contrast with those of all four simple arithmetic operators, that it can take various forms, and that the form used in our previous studies of artificial neural nets composed of local processors with both driving and contextual inputs is particularly well-suited to provide the distinctive capabilities of contextual modulation under a wide range of conditions. We argue that the decompositions reported here could be compared with those obtained from empirical neurobiological and psychophysical data under conditions thought to reflect contextual modulation. That would then shed new light on the underlying processes involved. Finally, we suggest that such decompositions could aid the design of context-sensitive machine learning algorithms.


title: Effects of a subadiabatic layer on convection and dynamos in spherical wedge simulations

abstract: We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We find that the convective heat transport is strongly concentrated toward the equatorial and polar regions in the cases without a substantial radiative layer below the convection zone. The presence of a stable layer below the convection zone significantly reduces the anisotropy of radial enthalpy transport. Furthermore, the dynamo solutions are sensitive to subtle changes in the convection zone structure. We find that the kinetic helicity changes sign in the deeper parts of the convection zone at high latitudes in all runs. This region expands progressively toward the equator in runs with a thicker stably stratified layer.


title: Moduli spaces of $G_{2}$ and $Spin(7)-$instantons on product manifolds

abstract: Let $X$ be a closed $6-$dimensional manifold with a half-closed $SU(3)-$structure. On the product manifold $X\times S^{1}$, with respect to the product $G_{2}-$structure and on a pullback vector bundle from $X$, we show that any $G_{2}-$instanton is equivalent to a Hermitian Yang-Mills connection on $X$ via a "broken gauge". This result reveals the topological type of the moduli of $G_{2}-$instantons on $X\times S^{1}$. In dimension $8$, similar result holds for moduli of $Spin(7)-$instantons. A generalization and an example are given.


title: Toolflows for Mapping Convolutional Neural Networks on FPGAs: A Survey and Future Directions

abstract: In the past decade, Convolutional Neural Networks (CNNs) have demonstrated state-of-the-art performance in various Artificial Intelligence tasks. To accelerate the experimentation and development of CNNs, several software frameworks have been released, primarily targeting power-hungry CPUs and GPUs. In this context, reconfigurable hardware in the form of FPGAs constitutes a potential alternative platform that can be integrated in the existing deep learning ecosystem to provide a tunable balance between performance, power consumption and programmability. In this paper, a survey of the existing CNN-to-FPGA toolflows is presented, comprising a comparative study of their key characteristics which include the supported applications, architectural choices, design space exploration methods and achieved performance. Moreover, major challenges and objectives introduced by the latest trends in CNN algorithmic research are identified and presented. Finally, a uniform evaluation methodology is proposed, aiming at the comprehensive, complete and in-depth evaluation of CNN-to-FPGA toolflows.


title: Neutrino tomography of the Earth

abstract: Cosmic-ray interactions with the nuclei of the Earth's atmosphere produce a flux of neutrinos in all directions with energies extending above the TeV scale. However, the Earth is not a fully transparent medium for neutrinos with energies above a few TeV. At these energies, the charged-current neutrino-nucleon cross section is large enough so that the neutrino mean-free path in a medium with the Earth's density is comparable to the Earth's diameter. Therefore, when neutrinos of these energies cross the Earth, there is a non-negligible probability for them to be absorbed. Since this effect depends on the distance traveled by neutrinos and on their energy, studying the zenith and energy distributions of TeV atmospheric neutrinos passing through the Earth offers an opportunity to infer the Earth's density profile. Here we perform an Earth tomography with neutrinos using actual data, the publicly available one-year through-going muon sample of the atmospheric neutrino data of the IceCube neutrino telescope. We are able to determine the mass of the Earth, its moment of inertia, the mass of the Earth's core and to establish the core is denser than the mantle, using weak interactions only, in a way completely independent from gravitational measurements. Our results confirm that this can be achieved with current neutrino detectors. This method to study the Earth's internal structure, complementary to the traditional one from geophysics based on seismological data, is starting to provide useful information and it could become competitive as soon as more statistics is available thanks to the current and larger future neutrino detectors.


title: Velocity-dependent Lyapunov exponents in many-body quantum, semi-classical, and classical chaos

abstract: The exponential growth or decay with time of the out-of-time-order commutator (OTOC) is one widely used diagnostic of many-body chaos in spatially-extended systems. In studies of many-body classical chaos, it has been noted that one can define a velocity-dependent Lyapunov exponent, $\lambda({\bf v})$, which is the growth or decay rate along "rays" at that velocity. We examine the behavior of $\lambda({\bf v})$ for a variety of many-body systems, both chaotic and integrable. The so-called light cone for the spreading of operators is defined by $\lambda({\bf \hat n}v_B({\bf \hat n}))=0$, with a generally direction-dependent "butterfly speed" $v_B({\bf \hat n})$. In spatially local systems, $\lambda(v)$ is negative outside the light cone where it takes the form $\lambda(v) \sim -(v-v_B)^{\alpha}$ near $v_b$, with the exponent $\alpha$ taking on various values over the range of systems we examine. The regime inside the light cone with positive Lyapunov exponents may only exist for classical, semi-classical or large-$N$ systems, but not for "fully quantum" chaotic systems with strong short-range interactions and local Hilbert space dimensions of order one.

categories: cond-mat.stat-mech cond-mat.str-el hep-th nlin.CD quant-ph


title: Local $\beta$-deformations and Yang-Baxter sigma model

abstract: Homogeneous Yang-Baxter (YB) deformation of AdS$_5\times$S$^5$ superstring is revisited. We calculate the YB sigma model action up to quadratic order in fermions and show that homogeneous YB deformations are equivalent to $\beta$-deformations of the AdS$_5\times$S$^5$ background when the classical $r$-matrices consist of bosonic generators. In order to make our discussion clearer, we discuss YB deformations in terms of the double-vielbein formalism of double field theory. We further provide an O(10,10)-invariant string action that reproduces the Green-Schwarz type II superstring action up to quadratic order in fermions. When an AdS background contains a non-vanishing $H$-flux, it is not straightforward to perform homogeneous YB deformations. In order to get any hint for such YB deformations, we study $\beta$-deformations of $H$-fluxed AdS backgrounds and obtain various solutions of (generalized) type II supergravity.


title: A generalization of the CHSH inequality self-testing maximally entangled states of any local dimension

abstract: For every $d \geq 2$, we present a generalization of the CHSH inequality with the property that maximal violation self-tests the maximally entangled state of local dimension $d$. This is the first example of a family of inequalities with this property. Moreover, we provide a conjecture for a family of inequalities generalizing the tilted CHSH inequalities, and we conjecture that for each pure bipartite entangled state there is an inequality in the family whose maximal violation self-tests it. All of these inequalities are inspired by the self-testing correlations of [Nat. Comm. 8, 15485 (2017)].


title: Critical behavior of the van der Waals bonded ferromagnet Fe$_{3-x}$GeTe$_2$

abstract: The critical properties of the single-crystalline van der Waals bonded ferromagnet Fe$_{3-x}$GeTe$_2$ were investigated by bulk dc magnetization around the paramagnetic (PM) to ferromagnetic (FM) phase transition. The Fe$_{3-x}$GeTe$_2$ single crystals grown by self-flux method with Fe deficiency $x \approx 0.36$ exhibit bulk FM ordering below $T_c = 152$ K. The M\"{o}ssbauer spectroscopy was used to provide information on defects and local atomic environment in such crystals. Critical exponents $\beta = 0.372(4)$ with a critical temperature $T_c = 151.25(5)$ K and $\gamma = 1.265(15)$ with $T_c = 151.17(12)$ K are obtained by the Kouvel-Fisher method whereas $\delta = 4.50(1)$ is obtained by a critical isotherm analysis at $T_c = 151$ K. These critical exponents obey the Widom scaling relation $\delta = 1+\gamma/\beta$, indicating self-consistency of the obtained values. With these critical exponents the isotherm $M(H)$ curves below and above the critical temperatures collapse into two independent universal branches, obeying the single scaling equation $m = f_\pm(h)$, where $m$ and $h$ are renormalized magnetization and field, respectively. The exponents determined in this study are close to those calculated from the results of the renormalization group approach for a heuristic model of three-dimensional Heisenberg ($d = 3, n = 3$) spins coupled with the attractive long-range interactions between spins that decay as $J(r)\approx r^{-(3+\sigma)}$ with $\sigma=1.89$.


title: A note on submanifolds of $\bar{M}^{2n+1}(f_1,f_2,f_3)$ with respect to certain connections

abstract: The present paper deals with some results of almsot semi-invariant submanifolds of generalized Sasakian-space-forms in \cite{ALEGRE3} with respect to semisymmetric metric connection, semisymmetric non-metric connection, Schouten-van Kampen connection and Tanaka-Webster connection.


title: Water transport on infinite graphs

abstract: If the nodes of a graph are considered to be identical barrels - featuring different water levels - and the edges to be (locked) water-filled pipes in between the barrels, consider the optimization problem of how much the water level in a fixed barrel can be raised with no pumps available, i.e. by opening and closing the locks in an elaborate succession. This model is related to an opinion formation process, the so-called Deffuant model. We consider i.i.d. random initial water levels and ask whether the supremum of achievable levels at a given node has a degenerate distribution, i.e. concentrates on a single value. This turns out to be the case for all infinite connected quasi-transitive graphs with exactly one exception: the two-sided infinite path.


title: Singular limit of BSDEs and Optimal control of two scale stochastic systems in infinite dimensional spaces

abstract: In this paper we study by probabilistic techniques the convergence of the value function for a two-scale, infinite-dimensional, stochastic controlled system as the ratio between the two evolution speeds diverges. The value function is represented as the solution of a \textit{backward stochastic differential equation} (BSDE) that it is shown to converge towards a \textit{reduced} BSDE. The noise is assumed to be additive both in the slow and the fast equations for the state. Some non degeneracy condition on the slow equation is required. The limit BSDE involves the solution of an \textit{ergodic} BSDE and is itself interpreted as the value function of an auxiliary stochastic control problem on a reduced state space.


title: Efficient Hardware Realization of Convolutional Neural Networks using Intra-Kernel Regular Pruning

abstract: The recent trend toward increasingly deep convolutional neural networks (CNNs) leads to a higher demand of computational power and memory storage. Consequently, the deployment of CNNs in hardware has become more challenging. In this paper, we propose an Intra-Kernel Regular (IKR) pruning scheme to reduce the size and computational complexity of the CNNs by removing redundant weights at a fine-grained level. Unlike other pruning methods such as Fine-Grained pruning, IKR pruning maintains regular kernel structures that are exploitable in a hardware accelerator. Experimental results demonstrate up to 10x parameter reduction and 7x computational reduction at a cost of less than 1% degradation in accuracy versus the un-pruned case.


title: Combining the color structures and intersection points of thick center vortices and low-lying Dirac modes

abstract: We investigate several examples of Yang-Mills gauge configurations containing center vortex structures, including intersection points between vortices and nontrivial color structures residing on the vortex world-surfaces. Various topological charge contributions of the color structures and intersection points are studied in these configurations. Low-lying eigenmodes of the (overlap) Dirac operator in the presence of these vortex backgrounds are analyzed. The results indicate characteristic properties for spontaneous chiral symmetry breaking.


title: The solution of the Einstein equation in the interior of the black holes by using arbitrary distribution functions

abstract: The aim of this paper is to obtain the solution of the Einstein equation in the interior of the black holes by using arbitrary distribution functions; corresponding to Gaussian, Rayleigh, Maxwell-Boltzmann and non-Gaussian distributions. Also we calculate the Hawking temperature, the mass and heat capacity for cosmological horizon and the black hole horizon.


title: The redshift evolution of rest-UV spectroscopic properties in Lyman Break Galaxies at z ~ 2-4

abstract: We present the first comprehensive evolutionary analysis of the rest-frame UV spectroscopic properties of star-forming galaxies at z ~ 2-4. We match samples at different redshifts in UV luminosity and stellar mass, and perform systematic measurements of spectral features and stellar population modeling. By creating composite spectra grouped according to Ly$\alpha$ equivalent width (EW), and various galaxy properties, we study the evolutionary trends among Ly$\alpha$, low- and high-ionization interstellar (LIS and HIS) absorption features, and integrated galaxy properties. We also examine the redshift evolution of Ly$\alpha$ and LIS absorption kinematics, and fine-structure emission EWs. The connections among the strengths of Ly$\alpha$, LIS lines, and dust extinction are redshift-independent, as is the decoupling of Ly$\alpha$ and HIS line strengths, and the bulk outflow kinematics as traced by LIS lines. Stronger Ly$\alpha$ emission is observed at higher redshift at fixed UV luminosity, stellar mass, SFR, and age. Much of this variation in average Ly$\alpha$ strength with redshift, and the variation in Ly$\alpha$ strength at fixed redshift, can be explained in terms of variations in neutral gas covering fraction and/or dust content in the ISM and CGM. However, based on the connection between Ly$\alpha$ and CIII] emission strengths, we additionally find evidence for variations in the intrinsic production rate of Ly$\alpha$ photons at the highest Ly$\alpha$ EWs. The challenge now is to understand the observed evolution in neutral gas covering fraction and dust extinction within a coherent model for galaxy formation, and make robust predictions for the escape of ionizing radiation at z > 6.


title: Linear stability analysis of collective neutrino oscillations without spurious modes

abstract: Collective neutrino oscillations are induced by the presence of neutrinos themselves. As such they are intrinsically nonlinear phenomena and are much more complex than linear counterparts such as the vacuum or MSW oscillations. They obey integro-differential equations, numerical solutions of which are also very challenging. If one focuses on the onset of the collective oscillations, on the other hand, the equations can be linearized and the technique of linear analysis can be employed. Unfortunately, however, it is well known that such an analysis, when applied with discretizations of continuous angular distributions, suffers from the appearance of so-called spurious modes, unphysical eigenmodes of the discretized linear equations. In this paper, we analyze in detail the origin of these unphysical modes and present a simple solution to this annoying problem. We have found that the spurious modes originate from the artificial production of pole singularities instead of a branch cut in the Riemann surface by the discretizations. The branching point singularities in the Riemann surface for the original undiscretized equations can be recovered by approximating the angular distributions with polynomials and then performing the integrals analytically. We demonstrate for some examples that this simple prescription removes the spurious modes indeed. We also propose an even simpler method: a piecewise linear approximation to the angular distribution. It is shown that the same methodology is applicable to the multi-energy case as well as to the dispersion relation approach that was proposed very recently.


title: Molecular Emission from a Galaxy Associated with a z~2.2 Damped Lyman-alpha Absorber

abstract: Using the Atacama Large Millimeter/sub-millimeter Array, we have detected CO(3-2) line and far-infrared continuum emission from a galaxy associated with a high-metallicity ([M/H] = -0.27) damped Ly-alpha absorber (DLA) at z =2.19289. The galaxy is located 3.5" away from the quasar sightline, corresponding to a large impact parameter of 30 kpc at the DLA redshift. We use archival Very Large Telescope-SINFONI data to detect Halpha emission from the associated galaxy, and find that the object is dusty, with a dust-corrected star formation rate of 110 +60 -30 Msun/yr. The galaxy's molecular mass is large, Mmol = (1.4 +- 0.2) x 10^11 x (\alpha_CO/4.3) x (0.57/r_31) Msun, supporting the hypothesis that high-metallicity DLAs arise predominantly near massive galaxies. The excellent agreement in redshift between the CO(3-2) line emission and low-ion metal absorption (~40 km/s) disfavors scenarios whereby the gas probed by the DLA shows bulk motion around the galaxy. We use Giant Metrewave Radio Telescope HI 21cm absorption spectroscopy to find that the HI along the DLA sightline must be warm, with a stringent lower limit on the spin temperature of T_s > 1895 x (f/0.93) K. The detection of CI absorption in the DLA, however, also indicates the presence of cold neutral gas. To reconcile these results requires that the cold components in the DLA contribute little to the HI column density, yet contain roughly 50% of the metals of the absorber, underlining the complex multi-phase nature of the gas surrounding high-z galaxies.


title: Characterization of the Infrared/X-ray sub-second variability for the black-hole transient GX 339-4

abstract: We present a detailed analysis of the X-ray/IR fast variability of the Black-Hole Transient GX 339-4 during its low/hard state in August 2008. Thanks to simultaneous high time-resolution observations made with the VLT and RXTE, we performed the first characterisation of the sub-second variability in the near-infrared band - and of its correlation with the X-rays - for a low-mass X-ray binary, using both time and frequency-domain techniques. We found a power-law correlation between the X-ray and infrared fluxes when measured on timescales of 16 seconds, with a marginally variable slope, steeper than the one found on timescales of days at similar flux levels. We suggest the variable slope - if confirmed - could be due to the infrared flux being a non-constant combination of both optically thin and optically thick synchrotron emission from the jet, as a result of a variable self-absorption break. From cross spectral analysis we found an approximately constant infrared time lag of $\approx$ 0.1s, and a very high coherence of $\approx$ 90 per cent on timescales of tens of seconds, slowly decreasing toward higher frequencies. Finally, we report on the first detection of a linear rms-flux relation in the emission from a low-mass X-ray binary jet, on timescales where little correlation is found between the X-rays and the jet emission itself. This suggests that either the inflow variations and jet IR emission are coupled by a non-linear or time-variable transform, or that the IR rms-flux relation is not transferred from the inflow to the jet, but is an intrinsic property of emission processes in the jet.


title: Extended Gauge Theory Deformations From~Flux~Backgrounds

abstract: We consider supersymmetric deformations of gauge theories in various dimensions obtained from a String Theory realisation of branes embedded in flux backgrounds. In particular we obtain deformations which take the form of Wilson line defects, where the R-symmetry is twisted into the gauge symmetry. Furthermore we construct higher-order generalisations, also expressed a twisting of the R-symmetry, that have symmetries associated to co-dimension two and three defects.


title: The Telltale Heartbeat: Detection and Characterization of Eccentric Orbiting Planets via Tides on their Host Star

abstract: We present an analytic description of tides raised on a star by a small orbiting body. In particular, we highlight the disproportionate effect of eccentricity and thus the scope for using these tides to detect and characterise the orbits of exoplanets and brown dwarfs. The tidal distortions of the star produced by an eccentric orbit are, in comparison to a circular orbit, much richer in detail, and potentially visible from any viewing angle. The magnitude of these variations is much larger than that in a circular orbit of the same semi-major axis. These variations are visible in both photometric and spectroscopic data, and dominate other regular sources of phase variability (e.g reflection and Doppler beaming) over a particularly interesting portion of parameter space. These tidal signatures will be a useful tool for planet detection on their own, and used in concert with other methods provide powerful constraints on planetary and stellar properties.


title: Spectroscopic Follow-Up of the Hercules Aquila Cloud

abstract: We designed a follow-up program to find the spectroscopic properties of the Hercules-Aquila Cloud (HAC) and test scenarios for its formation. We measured the radial velocities (RVs) of 45 RR Lyrae in the southern portion of the HAC using the facilities at the MDM observatory, producing the first large sample of velocities in the HAC. We found a double-peaked distribution in RVs, skewed slightly to negative velocities. We compared both the morphology of HAC projected onto the plane of the sky and the distribution of velocities in this structure outlined by RR Lyrae and other tracer populations at different distances to N-body simulations. We found that the behaviour is characteristic of an old, well-mixed accretion event with small apo-galactic radius. We cannot yet rule out other formation mechanisms for the HAC. However, if our interpretation is correct, HAC represents just a small portion of a much larger debris structure spread throughout the inner Galaxy whose distinct kinematic structure should be apparent in RV studies along many lines of sight.


title: Capturing near-equilibrium solutions: a comparison between high-order discontinuous Galerkin methods and well-balanced schemes

abstract: Equilibrium or stationary solutions usually proceed through the exact balance between hyperbolic transport terms and source terms. Such equilibrium solutions are affected by truncation errors that prevent any classical numerical scheme from capturing the evolution of small amplitude waves of physical significance. In order to overcome this problem, we compare two commonly adopted strategies: going to very high order and reduce drastically the truncation errors on the equilibrium solution, or design a specific scheme that preserves by construction the equilibrium exactly, the so-called well-balanced approach. We present a modern numerical implementation of these two strategies and compare them in details, using hydrostatic but also dynamical equilibrium solutions of several simple test cases. Finally, we apply our methodology to the simulation of a protoplanetary disc in centrifugal equilibrium around its star and model its interaction with an embedded planet, illustrating in a realistic application the strength of both methods.


title: On the quantum entropy function in 4d gauged supergravity

abstract: We analyze BPS black hole attractors in the conformal 4d gauged supergravity formalism and apply the technique known as supergravity localization in order to evaluate Sen's quantum entropy function in the $\mathrm{AdS}_2\!\times\!\mathrm{S}^2$ near-horizon geometry. Under certain assumptions, we reduce the exact expression of the functional integral to a finite-dimensional integral for a number of supersymmetric black holes in gauged supergravity with AdS asymptotics subject to a holographic description via a dual field theory. Examples include the asymptotically $\mathrm{AdS}_4\!\times\!\mathrm{S}^7$ Cacciatori-Klemm black holes in M-theory and the asymptotically $\mathrm{AdS}_5\!\times\!\mathrm{S}^5$ generalizations of Gutowski-Reall black holes and Benini-Bobev black strings in type IIB, as well as the recently constructed asymptotically $\mathrm{AdS}_4\!\times\!\mathrm{S}^6$ solutions in massive type IIA. Our results provide an important first step towards a gravitational counterpart to the exact evaluation of supersymmetric partition functions at finite $N$ for the holographically dual field theories in these examples.


title: Thermodynamics of bosonic systems in adS spacetime

abstract: We analyze the thermodynamics of massless bosonic systems in D-dimensional anti-de Sitter spacetime, considering scalar, electromagnetic, and gravitational fields. Their dynamics are described by Poschl-Teller effective potentials and quantized in a unified framework, with the determination of the associated energy spectra. From the microscopic description developed, a macroscopic thermodynamic treatment is proposed, where an effective volume in anti-de Sitter geometry is defined and a suitable thermodynamic limit is considered. Partition functions are constructed for the bosonic gases, allowing the determination of several thermodynamic quantities of interest. With the obtained results, general aspects of the thermodynamics are explored.


title: Lithium abundance patterns of late-F stars: an in-depth analysis of the lithium desert

abstract: We address the existence and origin of the lithium (Li) desert, a region in the Li - T_eff plane sparsely populated by stars. Here we analyze some of the explanations that have been suggested for this region, including mixing in the late main sequence, a Li dip origin for stars with low Li abundances in the region, and a possible relation with the presence of planets. To study the Li desert, we measured the atmospheric parameters and Li abundance of 227 late-F dwarfs and subgiants, chosen to be in the T_eff range of the desert and without previous Li abundance measurements. Subsequently, we complemented those with literature data to obtain a homogeneous catalog of 2318 stars, for which we compute masses and ages. We characterize stars surrounding the region of the Li desert. We conclude that stars with low Li abundances below the desert are more massive and more evolved than stars above the desert. Given the unexpected presence of low Li abundance stars in this effective temperature range, we concentrate on finding their origin. We conclude that these stars with low Li abundance do not evolve from stars above the desert: at a given mass, stars with low Li (i.e., below the desert) are more metal-poor. Instead, we suggest that stars below the Li desert are consistent with having evolved from the Li dip, discarding the need to invoke additional mixing to explain this feature. Thus, stars below the Li desert are not peculiar and are only distinguished from other subgiants evolved from the Li dip in that their combination of atmospheric parameters locates them in a range of effective temperatures where otherwise only high Li abundance stars would be found (i.e., stars above the desert).


title: ALMA Survey of Lupus Protoplanetary Disks II: Gas Disk Radii

abstract: We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply "Keplerian masking" to enhance the signal-to-noise ratios of our 12CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than mm dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission as well as the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, alpha_visc, finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 micron continuum observations, we also calculate the mm spectral index, alpha_mm, for 70 Lupus disks; we find an anti-correlation between alpha_mm and mm flux for low-mass disks (M_dust < 5), followed by a flattening as disks approach alpha_mm = 2, which could indicate faster grain growth in higher-mass disks, but may also reflect their larger optically thick components. In sum, this work demonstrates the continuous stream of new insights into disk evolution and planet formation that can be gleaned from unbiased ALMA disk surveys.


title: Investigating the Metallicity-Mixing Length Relation

abstract: Stellar models typically use the mixing length approximation as a way to implement convection in a simplified manner. While conventionally the value of the mixing length parameter, $\alpha$, used is the solar calibrated value, many studies have shown that other values of $\alpha$ are needed to properly model stars. This uncertainty in the value of the mixing length parameter is a major source of error in stellar models and isochrones. Using asteroseismic data, we determine the value of the mixing length parameter required to properly model a set of about 450 stars ranging in $\log g$, $T_{\mathrm{eff}}$, and $\mathrm{[Fe/H]}$. The relationship between the value of $\alpha$ required and the properties of the star is then investigated. For Eddington atmosphere, non-diffusion models, we find that the value of $\alpha$ can be approximated by a linear model, in the form of $\alpha/\alpha_{\odot}=5.426 -0.101 \log (g) -1.071 \log (T_{\mathrm{eff}}) + 0.437 (\mathrm{[Fe/H]})$. This process is repeated using a variety of model physics as well as compared to previous studies and results from 3D convective simulations.


title: Update on the X-ray variability plane for active galactic nuclei: The role of the obscuration

abstract: Scaling relations are the most powerful astrophysical tools to set constraints to the physical mechanisms of astro- nomical sources and to infer properties that cannot be accessed directly. We re-investigate here one of these scaling relations in active galactic nuclei (AGN); the so-called X-ray variability plane (or mass-luminosity-timescale relation, McHardy et al. 2006). This relation links the power-spectral density (PSD) break frequency with the super-massive black hole (SMBH) mass and the bolometric luminosity. We used available XMM -Newton observations of a sample of 22 AGN to study the PSD and spectra in short segments within each observation. This allows us to report for the first time that the PSD break frequency varies for each object, showing variations in 19 out of the 22 AGN analyzed. Our analysis of the variability plane confirms the relation between the break frequency and the SMBH mass and finds that the obscuration along the line of sight NH (or the variations on the obscuration using its standard deviation, $\rm{\Delta}$(NH)) is also a required parameter, at least for the range of frequencies analyzed here (3x10E-5-5x10E-2Hz). We constrain a new variability plane of the form: log($\nu_{Break}$) = (-0.589$\rm{\pm}$0.005) log(MBH ) + (0.10$\rm{\pm}$0.01) log(NH ) - (1.5$\rm{\pm}$0.3) (or log($\nu_{Break}$)=(-0.549$\rm{\pm}$0.009)log(MBH)+(0.56$\rm{\pm}$0.06)$\rm{\Delta}$(NH)+(0.19$\rm{\pm}$0.08)). The X-ray variability plane found by McHardy et al. (2006) is roughly recovered when we use unobscured segments. We speculate that this behavior is well explained if most of the reported frequencies are related to inner clouds (within 1pc), following Kepler orbits under the gravitational field of the SMBH.


title: Learning meets Assessment: On the relation between Item Response Theory and Bayesian Knowledge Tracing

abstract: Few models have been more ubiquitous in their respective fields than Bayesian knowledge tracing and item response theory. Both of these models were developed to analyze data on learners. However, the study designs that these models are designed for differ; Bayesian knowledge tracing is designed to analyze longitudinal data while item response theory is built for cross-sectional data. This paper illustrates a fundamental connection between these two models. Specifically, the stationary distribution of the latent variable and the observed response variable in Bayesian knowledge Tracing are related to an item response theory model. This connection between these two models highlights a key missing component: the role of education in these models. A research agenda is outlined which answers how to move forward with modeling learner data. %Furthermore, recent advances in network psychometrics demonstrate how this relationship can be exploited and generalized to a network model.


title: Galactic Rotation from Cepheids with Gaia DR2 and Effects of Non-Axisymmetry

abstract: We apply a simple axisymmetric disc model to 218 Galactic Cepheids whose accurate measurements of the distance and velocities are obtained by cross-matching an existing Cepheids catalogue with the Gaia DR2 data. Our model fit determines the "local centrifugal speed", $V_\mathrm{c}$ $-$ defined as the rotation speed required to balance the local radial gravitational force $-$ at the Sun's location to be $V_{c}(R_0)=236\pm 3$ km s$^{-1}$ and the Sun's azimuthal and radial peculiar motions to be $V_{\odot}=12.4\pm0.7$ km s$^{-1}$ and $U_{\odot}=7.7\pm0.9$ km s$^{-1}$, respectively. These results are obtained with strong priors on the solar radius, $R_0=8.2\pm0.1$ kpc, and Sun's angular rotation velocity, $\Omega_{\odot}=30.24\pm0.12$ km s$^{-1}$ kpc$^{-1}$. We also applied the axisymmetric model to mock data from an N-body/hydrodynamic simulation of a Milky Way-like galaxy with a bar and spiral arms. We find that our axisymmetric model fit to the young stars recovers the local centrifugal speed reasonably well, even in the face of significant non-axisymmetry. However, the local centrifugal speed determined from our Cepheid sample could suffer from systematic uncertainty as large as 6 km s$^{-1}$.


title: RankME: Reliable Human Ratings for Natural Language Generation

abstract: Human evaluation for natural language generation (NLG) often suffers from inconsistent user ratings. While previous research tends to attribute this problem to individual user preferences, we show that the quality of human judgements can also be improved by experimental design. We present a novel rank-based magnitude estimation method (RankME), which combines the use of continuous scales and relative assessments. We show that RankME significantly improves the reliability and consistency of human ratings compared to traditional evaluation methods. In addition, we show that it is possible to evaluate NLG systems according to multiple, distinct criteria, which is important for error analysis. Finally, we demonstrate that RankME, in combination with Bayesian estimation of system quality, is a cost-effective alternative for ranking multiple NLG systems.


title: Remarks on Donaldson's symplectic submanifolds

abstract: This paper presents a few remarks about the topology of symplectic hyperplane sections and the geometry of their complements. In particular, it contains a detailed proof of the following result already stated with hints in [Gi]: for sufficiently large degrees, the complements of Donaldson's symplectic hyperplane sections are naturally Weinstein --- and less naturally Stein --- manifolds.


title: Research and Development Workstation Environment: the new class of Current Research Information Systems

abstract: Against the backdrop of the development of modern technologies in the field of scientific research the new class of Current Research Information Systems (CRIS) and related intelligent information technologies has arisen. It was called - Research and Development Workstation Environment (RDWE) - the comprehensive problem-oriented information systems for scientific research and development lifecycle support. The given paper describes design and development fundamentals of the RDWE class systems. The RDWE class system's generalized information model is represented in the article as a three-tuple composite web service that include: a set of atomic web services, each of them can be designed and developed as a microservice or a desktop application, that allows them to be used as an independent software separately; a set of functions, the functional filling-up of the Research and Development Workstation Environment; a subset of atomic web services that are required to implement function of composite web service. In accordance with the fundamental information model of the RDWE class the system for supporting research in the field of ontology engineering - the automated building of applied ontology in an arbitrary domain area, scientific and technical creativity - the automated preparation of application documents for patenting inventions in Ukraine was developed. It was called - Personal Research Information System. A distinctive feature of such systems is the possibility of their problematic orientation to various types of scientific activities by combining on a variety of functional services and adding new ones within the cloud integrated environment. The main results of our work are focused on enhancing the effectiveness of the scientist's research and development lifecycle in the arbitrary domain area.


title: The role of surface tension gradient in determining microscopic dynamic contact angle

abstract: Following Gibb's interpretation of an interface as a dividing surface, we derive a model for the microscopic dynamic contact angle by writing a force balance for a control volume encompassing the interfaces and the contact line. In doing so we identify that, in addition to the surface tension of respective interfaces, the gradient of surface tension plays an important role in determining the dynamic contact angle. This is because not only does it contribute towards an additional force, but it also accounts for the deviation of local surface tension from its static equilibrium value. It is shown that this gradient in surface tension can be attributed to the convective acceleration in the vicinity of the contact line, which in turn is a direct result of varying degree of slip in that region. In addition, we provide evidence that this gradient in surface tension is one of the key factors responsible for the difference in contact angle at the leading and trailing edge, of a steadily moving contact line. These findings are validated using molecular dynamics simulations.


title: Multilevel Monte Carlo Method for Ergodic SDEs without Contractivity

abstract: This paper proposes a new multilevel Monte Carlo (MLMC) method for the ergodic SDEs which do not satisfy the contractivity condition. By introducing the change of measure technique, we simulate the path with contractivity and add the Radon-Nykodim derivative to the estimator. We can show the strong error of the path is uniformly bounded with respect to $T.$ Moreover, the variance of the new level estimators increase linearly in $T,$ which is a great reduction compared with the exponential increase in standard MLMC. Then the total computational cost is reduced to $O(\varepsilon^{-2}|\log \varepsilon|^{2})$ from $O(\varepsilon^{-3}|\log \varepsilon|)$ of the standard Monte Carlo method. Numerical experiments support our analysis.


title: Some Closure Results for Polynomial Factorization and Applications

abstract: In a sequence of seminal results in the 80's, Kaltofen showed that the complexity class VP is closed under taking factors. A natural question in this context is to understand if other natural classes of multivariate polynomials, for instance, arithmetic formulas, algebraic branching programs, bounded depth arithmetic circuits or the class VNP, are closed under taking factors. In this paper, we show that all factors of degree at most $\log^a n$ of polynomials with poly(n) size depth $k$ circuits have poly(n) size circuits of depth at most $O(k + a)$. This partially answers a question of Shpilka-Yehudayoff and has applications to hardness-randomness tradeoffs for bounded depth arithmetic circuits. More precisely, this shows that a superpolynomial lower bound for bounded depth arithmetic circuits, for a family of explicit polynomials of degree poly$(\log n)$ implies deterministic sub-exponential time algorithms for polynomial identity testing (PIT) for bounded depth arithmetic circuits. This is incomparable to a beautiful result of Dvir et al., where they showed that super-polynomial lower bounds for constant depth arithmetic circuits for any explicit family of polynomials (of potentially high degree) implies sub-exponential time deterministic PIT for bounded depth circuits of bounded individual degree. Thus, we remove the "bounded individual degree" condition in [DSY09] at the cost of strengthening the hardness assumption to hold for polynomials of low degree. As direct applications of our techniques, we also show that the complexity class VNP is closed under taking factors, thereby confirming a conjecture of B\"urgisser and get an alternate proof of the fact (first shown by Dutta et al.) that if a polynomial $Q$ of degree at most $d$ divides a polynomial $P$ computable by a formula of size $s$, then $Q$ has a formula of size at most poly$(s, d^{\log d}, deg(P))$.


title: Dual solutions for opposing mixed convection in porous media

abstract: The problem of steady mixed convection boundary-layer flow on a cooled vertical permeable circular cylinder embedded in a fluid-saturated porous medium is studied. Here, we evaluate the flow and heat transfer characteristics numerically for various values of the governing parameters and demonstrate the existence of dual solutions beyond a critical point.


title: CIM/E Oriented Graph Database Model Architecture and Parallel Network Topology Processing

abstract: CIM/E is an easy and efficient electric power model exchange standard between different Energy Management System vendors. With the rapid growth of data size and system complexity, the traditional relational database is not the best option to store and process the data. In contrast, the graph database and graph computation show their potential advantages to handle the power system data and perform real-time data analytics and computation. The graph concept fits power grid data naturally because of the fundamental structure similarity. Vertex and edge in the graph database can act as both a parallel storage unit and a computation unit. In this paper, the CIM/E data is modeled into the graph database. Based on this model, the parallel network topology processing algorithm is established and conducted by applying graph computation. The modeling and parallel network topology processing have been demonstrated in the modified IEEE test cases and practical Sichuan power network. The processing efficiency is greatly improved using the proposed method.


title: Mapping Disorder in Entropically Ordered Crystals

abstract: Systems of hard shapes crystallize due to entropy. How is entropy distributed among translational and rotational microscopic contributions? We answer this question by decomposing thermal fluctuation of crystals of hard hexagons into collective modes, a generalization and quantification of the Onsager picture of hard rod liquid crystals. We show that at densities both near densest packing and near the solid-hexatic melting transition, solids of hard regular hexagons hold most of their entropy in translational degrees of freedom.


title: Definable decompositions for graphs of bounded linear cliquewidth

abstract: We prove that for every positive integer k, there exists an MSO_1-transduction that given a graph of linear cliquewidth at most k outputs, nondeterministically, some cliquewidth decomposition of the graph of width bounded by a function of k. A direct corollary of this result is the equivalence of the notions of CMSO_1-definability and recognizability on graphs of bounded linear cliquewidth.


title: Low Rank plus Sparse Decomposition of ODFs for Improved Detection of Group-level Differences and Variable Correlations in White Matter

abstract: A novel approach is presented for group statistical analysis of diffusion weighted MRI datasets through voxelwise Orientation Distribution Functions (ODF). Recent advances in MRI acquisition make it possible to use high quality diffusion weighted protocols (multi-shell, large number of gradient directions) for routine in vivo study of white matter architecture. The dimensionality of these data sets is however often reduced to simplify statistical analysis. While these approaches may detect large group differences, they do not fully capitalize on all acquired image volumes. Incorporation of all available diffusion information in the analysis however risks biasing the outcome by outliers. Here we propose a statistical analysis method operating on the ODF, either the diffusion ODF or fiber ODF. To avoid outlier bias and reliably detect voxelwise group differences and correlations with demographic or behavioral variables, we apply the Low-Rank plus Sparse (L + S) matrix decomposition on the voxelwise ODFs which separates the sparse individual variability in the sparse matrix S whilst recovering the essential ODF features in the low-rank matrix L. We demonstrate the performance of this ODF L + S approach by replicating the established negative association between global white matter integrity and physical obesity in the Human Connectome dataset. The volume of positive findings agrees with and expands on the volume found by TBSS, Connectivity based fixel enhancement and Connectometry. In the same dataset we further localize the correlations of brain structure with neurocognitive measures such as fluid intelligence and episodic memory. The presented ODF L + S approach will aid in the full utilization of all acquired diffusion weightings leading to the detection of smaller group differences in clinically relevant settings as well as in neuroscience applications.


title: Noise-induced rectification in out-of-equilibrium structures

abstract: We consider the motion of overdamped particles on random potentials subjected to a Gaussian white noise and a time-dependent periodic external forcing. The random potential is modeled as the potential resulting from the interaction of a point particle with a random polymer. The random polymer is made up, by means of some stochastic process, from a finite set of possible monomer types. The process is assumed to reach a non-equilibrium stationary state, which means that every realization of a random polymer can be considered as an out-of-equilibrium structure. We show that the net flux of particles on this random medium is non-vanishing when the potential profile on every monomer is symmetric. We prove that this ratchet-like phenomenon is a consequence of the irreversibility of the stochastic process generating the polymer. On the contrary, when the process generating the polymer is at equilibrium (thus fulfilling the detailed balance condition) the system is unable to rectify the motion. We calculate the net flux of the particles in the adiabatic limit for a simple model and we test our theoretical predictions by means of Langevin dynamics simulations. We also show that, out of the adiabatic limit, the system also exhibits current reversals as well as non-monotonic dependence of the diffusion coefficient as a function of forcing amplitude.


title: Smartphone picture organization: A hierarchical approach

abstract: We live in a society where the large majority of the population has a camera-equipped smartphone. In addition, hard drives and cloud storage are getting cheaper and cheaper, leading to a tremendous growth in stored personal photos. Unlike photo collections captured by a digital camera, which typically are pre-processed by the user who organizes them into event-related folders, smartphone pictures are automatically stored in the cloud. As a consequence, photo collections captured by a smartphone are highly unstructured and because smartphones are ubiquitous, they present a larger variability compared to pictures captured by a digital camera. To solve the need of organizing large smartphone photo collections automatically, we propose here a new methodology for hierarchical photo organization into topics and topic-related categories. Our approach successfully estimates latent topics in the pictures by applying probabilistic Latent Semantic Analysis, and automatically assigns a name to each topic by relying on a lexical database. Topic-related categories are then estimated by using a set of topic-specific Convolutional Neuronal Networks. To validate our approach, we ensemble and make public a large dataset of more than 8,000 smartphone pictures from 10 persons. Experimental results demonstrate better user satisfaction with respect to state of the art solutions in terms of organization.


title: Excited States of Methylene, Polyenes, and Ozone from Heat-Bath Configuration Interaction

abstract: The electronically excited states of methylene (CH$_2$), ethylene (C$_2$H$_4$), butadiene (C$_4$H$_6$), hexatriene (C$_6$H$_8$), and ozone (O$_3$) have long proven challenging due to their complex mixtures of static and dynamic correlations. Semistochastic heat-bath configuration interaction (SHCI), which efficiently and systematically approaches the full configuration interaction (FCI) limit, is used to provide close approximations to the FCI energies in these systems. This article presents the largest FCI-level calculation to date -- on hexatriene using a polarized double-zeta basis (ANO-L-pVDZ), which gives rise to a Hilbert space containing more than $10^{38}$ determinants. These calculations give vertical excitation energies of 5.58 and 5.59 eV respectively for the $2^1{\rm A}_{\rm g}$ and $1^1{\rm B}_{\rm u}$ states, showing that they are nearly degenerate. The same excitation energies in butadiene/ANO-L-pVDZ were found to be 6.58 and 6.45 eV. In addition to these benchmarks, our calculations strongly support the presence of a previously hypothesized ring-minimum species of ozone that lies 1.3 eV higher than the open-ring minimum energy structure and is separated from it by a barrier of 1.11 eV.


title: Adaptive Tube-based Nonlinear MPC for Ecological Autonomous Cruise Control of Plug-in Hybrid Electric Vehicles

abstract: This paper proposes an adaptive tube-based nonlinear model predictive control (AT-NMPC) approach to the design of autonomous cruise control (ACC) systems. The proposed method utilizes two separate models to define the constrained receding horizon optimal control problem. A fixed nominal model is used to handle the problem constraints based on a robust tube-based approach. A separate adaptive model is used to define the objective function, which utilizes least square online parameter estimators for adaption. By having two separate models, this method takes into account uncertainties, modeling errors and delayed data in the design of the controller and guaranties robust constraint handling, while adapting to them to improve control performance. Furthermore, to be able implement the designed AT-NMPC in real-time, a Newton/GMRES fast solver is employed to solve the optimization problem. Simulations performed on a high-fidelity model of the baseline vehicle, the Toyota plug-in Prius, which is a plug-in hybrid electric vehicle (PHEV), show that the proposed controller is able to handle the defined constraints in the presence of uncertainty, while improving the energy cost of the trip. Moreover, the result of the hardware-in-loop experiment demonstrates the performance of the proposed controller in real time application.


title: Closed form expressions for Appell polynomials

abstract: We show that any Appell sequence can be written in closed form as a forward difference transformation of the identity. Such transformations are actually multipliers in the abelian group of the Appell polynomials endowed with the operation of binomial convolution. As a consequence, we obtain explicit expressions for higher order convolution identities referring to various kinds of Appell polynomials in terms of the Stirling numbers. Applications of the preceding results to generalized Bernoulli and Apostol-Euler polynomials of real order are discussed in detail.


title: $L^2$ concentration of blow-up solutions for the mass-critical NLS with inverse-square potential

abstract: In this paper, we prove a refined version of a compactness lemma and we use it to establish mass-concentration for the focusing nonlinear Schr\"{o}dinger equation with an inverse-square potential.


title: Analog simulator of integro-differential equations with classical memristors

abstract: An analog computer makes use of continuously changeable quantities of a system, such as its electrical, mechanical, or hydraulic properties, to solve a given problem. While these devices are usually computationally more powerful than their digital counterparts, they suffer from analog noise which does not allow for error control. We will focus on analog computers based on active electrical networks comprised of resistors, capacitors, and operational amplifiers which are capable of simulating any linear ordinary differential equation. However, the class of nonlinear dynamics they can solve is limited. In this work, by adding memristors to the electrical network, we show that the analog computer can simulate a large variety of linear and nonlinear integro-differential equations by carefully choosing the conductance and the dynamics of the memristor state variable. To the best of our knowledge, this is the first time that circuits based on memristors are proposed for simulations. We study the performance of these analog computers by simulating integro-differential models related to fluid dynamics, nonlinear Volterra equations for population growth, and quantum models describing non-Markovian memory effects, among others. Finally, we perform stability tests by considering imperfect analog components, obtaining robust solutions with up to $13\%$ relative error for relevant timescales.


title: An Optimal Algorithm to Compute the Inverse Beacon Attraction Region

abstract: The beacon model is a recent paradigm for guiding the trajectory of messages or small robotic agents in complex environments. A beacon is a fixed point with an attraction pull that can move points within a given polygon. Points move greedily towards a beacon: if unobstructed, they move along a straight line to the beacon, and otherwise they slide on the edges of the polygon. The Euclidean distance from a moving point to a beacon is monotonically decreasing. A given beacon attracts a point if the point eventually reaches the beacon. The problem of attracting all points within a polygon with a set of beacons can be viewed as a variation of the art gallery problem. Unlike most variations, the beacon attraction has the intriguing property of being asymmetric, leading to separate definitions of attraction region and inverse attraction region. The attraction region of a beacon is the set of points that it attracts. It is connected and can be computed in linear time for simple polygons. By contrast, it is known that the inverse attraction region of a point---the set of beacon positions that attract it---could have $\Omega(n)$ disjoint connected components. In this paper, we prove that, in spite of this, the total complexity of the inverse attraction region of a point in a simple polygon is linear, and present a $O(n \log n)$ time algorithm to construct it. This improves upon the best previous algorithm which required $O(n^3)$ time and $O(n^2)$ space. Furthermore we prove a matching $\Omega(n\log n)$ lower bound for this task in the algebraic computation tree model of computation, even if the polygon is monotone.


title: Control Inversion: A Clustering-Based Method for Distributed Wide-Area Control of Power Systems

abstract: Wide-area control (WAC) has been shown to be an effective tool for damping low-frequency oscillations in power systems. In the current state of art, WAC is challenged by two main factors - namely, scalability of design and complexity of implementation. In this paper we present a control design called control inversion that bypasses both of these challenges using the idea of clustering. The basic philosophy behind this method is to project the original power system model into a lower-dimensional state-space through clustering and aggregation of generator states, and then designing an LQR controller for the lower-dimensional model. This controller is finally projected back to the original coordinates for wide-area implementation. The main problem is, therefore, posed as finding the projection which best matches the closed-loop performance of the WAC controller with that of a reference LQR controller for damping low-frequency oscillations. We verify the effectiveness of the proposed design using the NPCC 48-machine power system model.


title: Average Cost of QuickXsort with Pivot Sampling

abstract: QuickXsort is a strategy to combine Quicksort with another sorting method X, so that the result has essentially the same comparison cost as X in isolation, but sorts in place even when X requires a linear-size buffer. We solve the recurrence for QuickXsort precisely up to the linear term including the optimization to choose pivots from a sample of k elements. This allows to immediately obtain overall average costs using only the average costs of sorting method X (as if run in isolation). We thereby extend and greatly simplify the analysis of QuickHeapsort and QuickMergesort with practically efficient pivot selection, and give the first tight upper bounds including the linear term for such methods.