title: Optimization of a Dynamic Profit Function using Euclidean Path Integral

abstract: A Euclidean path integral is used to find an optimal strategy for a firm under a Walrasian system, Pareto optimality and a non-cooperative feedback Nash Equilibrium. We define dynamic optimal strategies and develop a Feynman type path integration method to capture all non-additive convex strategies. We also show that the method can solve the non-linear case, for example Merton-Garman-Hamiltonian system, which the traditional Pontryagin maximum principle cannot solve in closed form. Furthermore, under Walrasian system we are able to solve for the optimal strategy under a linear constraint with a linear objective function with respect to strategy.


title: Distinguishing the nature of ''ambiguous'' merging systems hosting a neutron star: GW190425 in low-latency

abstract: GW190425 is the newly discovered gravitational wave (GW) source consistent with a neutron star-neutron star merger with chirp mass of $1.44\pm0.02M_\odot.$ This value falls in the $ambiguous$ interval as from the GW signal alone we can not rule out the presence of a black hole in the binary. In this case, the system would host a neutron star and a very light stellar black hole, with mass close to the maximum value for neutron stars, filling the $mass$ $gap$. No electromagnetic counterpart is firmly associated with this event, due to the poorly informative sky localisation and larger distance, compared to GW/GRB170817. We construct here kilonova light curve models for GW190425, for both double neutron star and black hole-neutron star systems, considering two equations of state consistent with current constraints from the signals of GW170817/GW190425 and the NICER results, including black hole spin effects and assuming a new formula for the mass of the ejecta. The putative presence of a light black hole in GW190425 would have produced a brighter kilonova emission compared to the double neutron star case, letting us to distinguish the nature of the companion to the neutron star. Concerning candidate counterparts of GW190425, classified later on as supernovae, our models could have discarded two transients detected in their early $r$-band evolution. Combining the chirp mass and luminosity distance information from the GW signal with a library of kilonova light curves helps identifying the electromagnetic counterpart early on. We remark that the release in low latency of the chirp mass in this interval of $ambiguous$ values appears to be vital for successful electromagnetic follow-ups.


title: Hilbert space average of transition probabilities

abstract: The typicality approach and the Hilbert space averaging method as its technical manifestation are important concepts of quantum statistical mechanics. Extensively used for expectation values we extend them in this paper to transition probabilities. In this context we also find that the transition probability of two random uniformly distributed states is connected to the spectral statistics of the considered operator. Furthermore, within our approach we are capable to consider distributions of matrix elements between states, that are not orthogonal. We will demonstrate our quite general result numerically for a kicked spin chain in the integrable resp. chaotic regime.


title: Relaxation and vibrational properties in metal alloys and other disordered systems

abstract: The relaxation dynamics and the vibrational spectra of amorphous solids, such as metal alloys, have been intensely investigated as well separated topics in the past. The aim of this review is to summarize recent results in both these areas in an attempt to establish, or unveil, deeper connections between the two phenomena of relaxation and vibration. Theoretical progress in the area of slow relaxation dynamics of liquid and glassy systems and in the area of vibrational spectra of glasses and liquids is reviewed. After laying down a generic modelling framework to connect vibration and relaxation, the physics of metal alloys is considered where the emergence of power-law exponents has been identified both in the vibrational density of states (VDOS) as well as in density correlations. Also, theoretical frameworks which connect the VDOS to the relaxation behaviour and mechanical viscoelastic response in metallic glasses are reviewed. The same generic interpretative framework is then applied to the case of molecular glass formers where the emergence of stretched-exponential relaxation in dielectric relaxation can be put in quantitative relation with the VDOS by means of memory-function approaches. Further connections between relaxation and vibration are provided by the study of phonon linewidths in liquids and glasses, where a natural starting point is given by hydrodynamic theories. Finally, an agenda of outstanding issues including the appearance of compressed exponential relaxation in the intermediate scattering function of experimental and simulated systems (metal alloys, colloidal gels, jammed packings) is presented in light of available (or yet to be developed) mathematical models, and compared to non-exponential behaviour measured with macroscopic means such as mechanical spectroscopy/rheology.

categories: cond-mat.mtrl-sci cond-mat.dis-nn cond-mat.soft cond-mat.stat-mech


title: It's Not What Machines Can Learn, It's What We Cannot Teach

abstract: Can deep neural networks learn to solve any task, and in particular problems of high complexity? This question attracts a lot of interest, with recent works tackling computationally hard tasks such as the traveling salesman problem and satisfiability. In this work we offer a different perspective on this question. Given the common assumption that $\textit{NP} \neq \textit{coNP}$ we prove that any polynomial-time sample generator for an $\textit{NP}$-hard problem samples, in fact, from an easier sub-problem. We empirically explore a case study, Conjunctive Query Containment, and show how common data generation techniques generate biased datasets that lead practitioners to over-estimate model accuracy. Our results suggest that machine learning approaches that require training on a dense uniform sampling from the target distribution cannot be used to solve computationally hard problems, the reason being the difficulty of generating sufficiently large and unbiased training sets.


title: Plasmonic Quantum Yield Enhancement of a Single Molecule Near a Nanoegg

abstract: We investigate the impact of the dipole-active modes formed via the mode-mixing of the dipole mode with higher-order surface plasmon modes of a nanoegg on the radiative decay rate and quantum yield of an excited molecule near the nanoegg. The Purcell factor, rate of power dissipation by the emitter, antenna efficiency of the nanoegg, as well as quantum yield enhancement of the emitter, were studied using the quasistatic approximation and the semiclassical theory of radiation, following the Gersten-Nitzan and Ford-Weber approaches. Compared to the concentric nanoshell, we show that the dielectric core-metallic shell nanoegg is a more efficient plasmonic nanoantenna for radiative decay rate enhancement of single emitters. The quantum yield of the emitter was found to be more enhanced near the nanoshell, while its emission rate was found to be more enhanced near the nanoegg.


title: Chern-Simons Supergravity on Supergroup Manifolds

abstract: We construct N=1 d=3 AdS supergravity within the group manifold approach and compare it with Achucarro-Townsend Chern-Simons formulation of the same theory. We clarify the relation between the off-shell super gauge transformations of the Chern-Simons theory and the off-shell worldvolume supersymmetry transformations of the group manifold action. We formulate the Achucarro-Townsend model in a double supersymmetric action where the Chern-Simons theory with a supergroup gauge symmetry is constructed on a supergroup manifold. This framework is useful to establish a correspondence of degrees of freedom and auxiliary fields between the two descriptions of d=3 supergravity.


title: Quark structure of the $\chi_{\rm c}(3P)$ and $X(4274)$ resonances and their strong and radiative decays

abstract: We calculate the masses of $\chi_{\rm c}(3P)$ states with threshold corrections in a coupled-channel model. The model was recently applied to the description of the properties of $\chi_{\rm c}(2P)$ and $\chi_{\rm b}(3P)$ multiplets [Phys.\ Lett.\ B {\bf 789}, 550 (2019)]. We also compute the open-charm strong decay widths of the $\chi_{\rm c}(3P)$ states and their radiative transitions. According to our predictions, the $\chi_{\rm c}(3P)$ states should be dominated by the charmonium core, but they may also show small meson-meson components. The $X(4274)$ is interpreted as a $c \bar c$ $\chi_{\rm c1}(3P)$ state. More informations on the other members of the $\chi_{\rm c}(3P)$ multiplet, as well as a more rigorous analysis of the $X(4274)$'s decay modes, are needed to provide further indications on the quark structure of the previous resonance.


title: Addressing Some Limitations of Transformers with Feedback Memory

abstract: Transformers have been successfully applied to sequential, auto-regressive tasks despite being feedforward networks. Unlike recurrent neural networks, Transformers use attention to capture temporal relations while processing input tokens in parallel. While this parallelization makes them computationally efficient, it restricts the model from fully exploiting the sequential nature of the input. The representation at a given layer can only access representations from lower layers, rather than the higher level representations already available. In this work, we propose the Feedback Transformer architecture that exposes all previous representations to all future representations, meaning the lowest representation of the current timestep is formed from the highest-level abstract representation of the past. We demonstrate on a variety of benchmarks in language modeling, machine translation, and reinforcement learning that the increased representation capacity can create small, shallow models with much stronger performance than comparable Transformers.


title: Inexact Tensor Methods with Dynamic Accuracies

abstract: In this paper, we study inexact high-order Tensor Methods for solving convex optimization problems with composite objective. At every step of such methods, we use approximate solution of the auxiliary problem, defined by the bound for the residual in function value. We propose two dynamic strategies for choosing the inner accuracy: the first one is decreasing as $1/k^{p + 1}$, where $p \geq 1$ is the order of the method and $k$ is the iteration counter, and the second approach is using for the inner accuracy the last progress in the target objective. We show that inexact Tensor Methods with these strategies achieve the same global convergence rate as in the error-free case. For the second approach we also establish local superlinear rates (for $p \geq 2$), and propose the accelerated scheme. Lastly, we present computational results on a variety of machine learning problems for several methods and different accuracy policies.


title: Rayleigh-B\'{e}nard Convection in Strong Vertical Magnetic Field: Flow Structure and Verification of Numerical Method

abstract: Direct numerical simulations are performed to study turbulent Rayleigh-B\'{e}nard convection in a vertical cylindrical cavity with uniform axial magnetic field. Flows at high Hartmann and Rayleigh numbers are considered. The calculations reveal that, similarly to the behavior observed in Rayleigh-B\'{e}nard convection with strong rotation, flows at strong magnetic field develop a central vortex, while the heat transfer is suppressed.


title: Learning to Simulate Complex Physics with Graph Networks

abstract: Here we present a machine learning framework and model implementation that can learn to simulate a wide variety of challenging physical domains, involving fluids, rigid solids, and deformable materials interacting with one another. Our framework---which we term "Graph Network-based Simulators" (GNS)---represents the state of a physical system with particles, expressed as nodes in a graph, and computes dynamics via learned message-passing. Our results show that our model can generalize from single-timestep predictions with thousands of particles during training, to different initial conditions, thousands of timesteps, and at least an order of magnitude more particles at test time. Our model was robust to hyperparameter choices across various evaluation metrics: the main determinants of long-term performance were the number of message-passing steps, and mitigating the accumulation of error by corrupting the training data with noise. Our GNS framework advances the state-of-the-art in learned physical simulation, and holds promise for solving a wide range of complex forward and inverse problems.


title: The Automated Inspection of Opaque Liquid Vaccines

abstract: In the pharmaceutical industry the screening of opaque vaccines containing suspensions is currently a manual task carried out by trained human visual inspectors. We show that deep learning can be used to effectively automate this process. A moving contrast is required to distinguish anomalies from other particles, reflections and dust resting on a vial's surface. We train 3D-ConvNets to predict the likelihood of 20-frame video samples containing anomalies. Our unaugmented dataset consists of hand-labelled samples, recorded using vials provided by the HAL Allergy Group, a pharmaceutical company. We trained ten randomly initialized 3D-ConvNets to provide a benchmark, observing mean AUROC scores of 0.94 and 0.93 for positive samples (containing anomalies) and negative (anomaly-free) samples, respectively. Using Frame-Completion Generative Adversarial Networks we: (i) introduce an algorithm for computing saliency maps, which we use to verify that the 3D-ConvNets are indeed identifying anomalies; (ii) propose a novel self-training approach using the saliency maps to determine if multiple networks agree on the location of anomalies. Our self-training approach allows us to augment our data set by labelling 217,888 additional samples. 3D-ConvNets trained with our augmented dataset improve on the results we get when we train only on the unaugmented dataset.


title: Dynamic Normalization for Compact Binary Coalescence Searches in Non-Stationary Noise

abstract: The output of gravitational-wave interferometers, such as LIGO and Virgo, can be highly non-stationary. Broadband detector noise can affect the detector sensitivity on the order of tens of seconds. Gravitational-wave transient searches, such as those for colliding black holes, estimate this noise in order to identify gravitational-wave events. During times of non-stationarity we see a higher rate of false events being reported. To accurately separate signal from noise, it is imperative to incorporate the changing detector state into gravitational-wave searches. We develop a new statistic which estimates the variation of the interferometric detector noise. We use this statistic to re-rank candidate events identified during LIGO-Virgo's second observing run by the PyCBC search pipeline. This results in a 5% improvement in the sensitivity volume for low mass binaries, particularly binary neutron stars mergers.


title: Improvement, generalization, and scheme conversion of Wilson-line operators on the lattice in the auxiliary field approach

abstract: Nonlocal quark bilinear operators connected by link paths are used for studying parton distribution functions (PDFs) and transverse momentum-dependent PDFs of hadrons using lattice QCD. The nonlocality makes it difficult to understand the renormalization and improvement of these operators using standard methods. In previous work, we showed that by introducing an auxiliary field on the lattice, one can understand an on-axis Wilson-line operator as the product of two local operators in an extended theory. In this paper, we provide details about the calculation in perturbation theory of the factor for conversion from our lattice-suitable renormalization scheme to the MS-bar scheme. Extending our work, we study Symanzik improvement of the extended theory to understand the pattern of discretization effects linear in the lattice spacing, $a$, which are present even if the lattice fermion action exactly preserves chiral symmetry. This provides a prospect for an eventual $O(a)$ improvement of lattice calculations of PDFs. We also generalize our approach to apply to Wilson lines along lattice diagonals and to piecewise-straight link paths.


title: Quantum Detection using Magnetic Avalanches in Single-Molecule Magnets

abstract: The detection of a single quantum of energy with high efficiency and low false positive rate is of considerable scientific interest, from serving as single quantum sensors of optical and infra-red photons to enabling the direct detection of low-mass dark matter. We report the first experimental demonstration of magnetic avalanches induced by scattering of quanta in single-molecule magnet (SMM) crystals made of Mn12-acetate, establishing the use of SMMs as particle detectors for the first time. While the current setup has an energy threshold in the MeV regime, our results motivate the exploration of a wide variety of SMMs whose properties could allow for detection of sub-eV energy depositions.

categories: physics.ins-det cond-mat.mes-hall hep-ex hep-ph


title: Search for CME in U+U and Au+Au collisions in STAR with different approaches of handling backgrounds

abstract: The chiral magnetic effect (CME) refers to charge separation along a strong magnetic field between left- and right-handed quarks, caused by interactions with topological gluon fields from QCD vacuum fluctuations. We present two approaches to handle the dominant elliptic flow ($v_2$) background in the three-particle correlator ($\Delta\gamma_{112}$), sensitive to CME. In the first approach, we present the $\Delta\gamma_{112}$ and $\Delta\gamma_{123}$ measurements in U+U and Au+Au collisions. While hydrodynamic simulations including resonance decays and local charge conservation predict that $\Delta\gamma_{112}$ scaled by $N_{\rm part}/v_2$ will be similar in U+U and Au+Au collisions, the projected B-field exhibits a distinct difference between the two systems and with varying $N_{\rm part}$. Therefore, U+U and Au+Au collisions provide configurations with different expectations for both CME signal and background. Moreover, the three-particle observable $\Delta\gamma_{123}$ scaled by $N_{\rm part}/v_3$ provide baseline measurement for only the background. In the second approach, we handle the $v_2$ background by measuring $\Delta\gamma_{112}$ with respect to the planes of spectators measured by Zero Degree Calorimeters and participants measured by Time Projection Chamber. These measurements contain different amounts of contributions from CME signal (along B-field, due to spectators) and $v_2$ background (determined by the participant geometry). With the two $\Delta\gamma_{112}$ measurements, the possible CME signal and the background contribution can be determined. We report such a measurement in Au+Au collisions at $\sqrt{s_{NN}}=$ 27 GeV with the newly installed event plane detector, and report the new findings in U+U system where the spectator-participant plane correlations are expected to differ from those in Au+Au collisions.


title: Multi-scale simulations of particle acceleration in astrophysical systems

abstract: This review aims at providing an up-to-date status and a general introduction to the subject of the numerical study of energetic particle acceleration and transport in turbulent astrophysical flows. The subject is also complemented by a short overview of recent progresses obtained in the domain of laser plasma experiments. We review the main physical processes at the heart of the production of a non-thermal distribution in both Newtonian and relativistic astrophysical flows, namely the first and second order Fermi acceleration processes. We also discuss shock drift and surfing acceleration, two processes important in the context of particle injection in shock acceleration. We analyze with some details the particle-in-cell (PIC) approach used to describe particle kinetics. We review the main results obtained with PIC simulations in the recent years concerning particle acceleration at shocks and in reconnection events. The review discusses the solution of Fokker-Planck problems with application to the study of particle acceleration at shocks but also in hot coronal plasmas surrounding compact objects. We continue by considering large scale physics. We describe recent developments in magnetohydrodynamic (MHD) simulations. We give a special emphasize on the way energetic particle dynamics can be coupled to MHD solutions either using a multi-fluid calculation or directly coupling kinetic and fluid calculations. This aspect is mandatory to investigate the acceleration of particles in the deep relativistic regimes to explain the highest Cosmic Ray energies.


title: Single-step exfoliation and functionalization of few-layers black phosphorus and its application for polymer composites

abstract: We demonstrate the large scale exfoliation of bulk black phosphorus into few-layers black phosphorus. We use wet-jet mill (WJM) as an industrially scalable exfoliation technique. In order to avoid the degradation by oxidation of the exfoliated flakes, the few-layers black phosphorus is protected by functionalisation (with 4-bromobenzene-diazonium tetrafluoroborate) during the exfoliation process. We used the exfoliated and functionalised flakes as fillers in a polymethyl methacrylate matrix to test the mechanical endurance, obtaining an improvement in the Young modulus of 106% compared to the bare polymer. These data support the development of black phosphorus materials as 2D-polymer fillers and performance-enhancers.


title: Configuration Entropy for Quarkonium in a Finite Density Plasma

abstract: In the recent years many examples appeared in the literature where the configuration entropy (CE), introduced by Gleiser and Stamatopoulos, plays the role of an indicator of stability of physical systems. It was observed that, comparing states of the same system, the lower is the value of the CE, the more stable is the state. In this work we investigate the behaviour of the differential configuration entropy (DCE) , that is appropriate for systems with continuous degrees of freedom, in a new context. We consider quasi-states of quarkonium (a vector meson made of a heavy quark anti-quark pair) inside a plasma at finite density. It is known that the density increases the dissociation effect for quasi-particles inside a plasma. So, increasing the density of a thermal medium corresponds to reducing the stability of the quasi-particles. In order to investigate how this situation is translated in the Configutation Entropy context, we use a recently developed holographic AdS/QCD model for heavy vector mesons. The quasi-normal modes describing the quasi-states are obtained and the corresponding DCE is calculated. We find, for bottomonium and charmonium $1 S$ quasi-states, that the DCE increases with the quark density, or quark chemical potential, of the medium. This result shows that the DCE works again as an indicator of stability, represented in this case by the dissociation effect associated with the density.


title: Directed Acyclic Graphs and causal thinking in clinical risk prediction modeling

abstract: Background: In epidemiology, causal inference and prediction modeling methodologies have been historically distinct. Directed Acyclic Graphs (DAGs) are used to model a priori causal assumptions and inform variable selection strategies for causal questions. Although tools originally designed for prediction are finding applications in causal inference, the counterpart has remained largely unexplored. The aim of this theoretical and simulation-based study is to assess the potential benefit of using DAGs in clinical risk prediction modeling. Methods and Findings: We explore how incorporating knowledge about the underlying causal structure can provide insights about the transportability of diagnostic clinical risk prediction models to different settings. A single-predictor model in the causal direction is likely to have better transportability than one in the anticausal direction. We further probe whether causal knowledge can be used to improve predictor selection. We empirically show that the Markov Blanket, the set of variables including the parents, children, and parents of the children of the outcome node in a DAG, is the optimal set of predictors for that outcome. Conclusions: Our findings challenge the generally accepted notion that a change in the distribution of the predictors does not affect diagnostic clinical risk prediction model calibration if the predictors are properly included in the model. Furthermore, using DAGs to identify Markov Blanket variables may be a useful, efficient strategy to select predictors in clinical risk prediction models if strong knowledge of the underlying causal structure exists or can be learned.


title: Interplanetary drivers of the magnetospheric disturbances: A brief review of incorrect approaches

abstract: One of the most promising areas of research in solar-terrestrial physics is the comparison of the responses of the magnetosphere-ionosphere-atmosphere system to various interplanetary disturbances (the so-called <<interplanetary drivers>>). Numerous studies show that different types of drivers cause a different reaction of the system for identical IMF variations. At the same time, the number of incorrect approaches in this direction of research has increased. These errors can be attributed to 4 large classes. (1) The first class includes works whose authors uncritically reacted to previously published works with incorrect driver identification and use incorrect results in their work. (2) Some authors used the wrong criteria and incorrectly determined the types of drivers. (3) Very often, authors associate the disturbance of the magnetosphere-ionosphere-atmosphere system caused by a complex driver (by a sequence of single drivers) with one of the drivers, ignoring the complex nature. For example, magnetic storm are often caused by compression region Sheath in front of the interplanetary CME (ICME), but the authors consider this events as so-called <<CME-induced>> storm, not <<Sheath-induced>> storm. (4) Finally, there is a <<lost driver>> of magnetospheric disturbances: some authors simply do not consider the compression region Sheath before ICME if there is no interplanetary shock (IS) before Sheath, although this type of driver, <<Sheath without IS>>, generates about 10% of moderate and strong magnetic storms.


title: Generic measures for translation surface flows

abstract: We consider straight line flows on a translation surface that are minimal but not uniquely ergodic. We give bounds for the number of generic invariant probability measures.


title: Layer-dependence of dielectric response and water-enhanced ambient degradation of highly-anisotropic black As

abstract: Black arsenic (BAs) is a van der Waals layered material with a puckered honeycomb structure and has received increased interest due to its anisotropic properties and promising performance in devices. Here, crystalline structure, thickness-dependent dielectric responses, and ambient stability of BAs nanosheets are investigated using STEM imaging and spectroscopy. Atomic-resolution HAADF-STEM images directly visualize the three-dimensional structure and evaluate the degree of anisotropy. STEM-EELS is used to measure the dielectric response of BAs as a function of the number of layers. Finally, BAs degradation under different ambient environments is studied highlighting high sensitivity to moisture in the air.


title: Excitation Transport in Molecular Aggregates with thermal motion

abstract: Molecular aggregates can under certain conditions transport electronic excitation energy over large distances due to dipole-dipole interactions. Here, we explore to what extent thermal motion of entire monomers can guide or enhance this excitation transport. The motion induces changes of aggregate geometry and hence modifies exciton states. Under certain conditions, excitation energy can thus be transported by the aggregate adiabatically, following a certain exciton eigenstate. While such transport is always slower than direct migration through dipole-dipole interactions, we show that transport through motion can yield higher transport efficiencies in the presence of on-site energy disorder than the static counterpart. For this we consider two simple models of molecular motion: (i) longitudinal vibrations of the monomers along the aggregation direction within their inter-molecular binding potential and (ii) torsional motion of planar monomers in a plane orthogonal to the aggregation direction. The parameters and potential shapes used are relevant to dye-molecule aggregates. We employ a quantum-classical method, in which molecules move through simplified classical molecular dynamics, while the excitation transport is treated quantum mechanically using Schr\"odinger's equation. For both models we find parameter regimes in which the motion enhances excitation transport, however these are more realistic for the torsional scenario, due to the limited motional range in a typical Morse type inter-molecular potential. We finally show that the transport enhancement can be linked to adiabatic quantum dynamics. This transport enhancement through adiabatic motion appears a useful resource to combat exciton trapping by disorder.


title: Guider l'attention dans les modeles de sequence a sequence pour la prediction des actes de dialogue

abstract: The task of predicting dialog acts (DA) based on conversational dialog is a key component in the development of conversational agents. Accurately predicting DAs requires a precise modeling of both the conversation and the global tag dependencies. We leverage seq2seq approaches widely adopted in Neural Machine Translation (NMT) to improve the modelling of tag sequentiality. Seq2seq models are known to learn complex global dependencies while currently proposed approaches using linear conditional random fields (CRF) only model local tag dependencies. In this work, we introduce a seq2seq model tailored for DA classification using: a hierarchical encoder, a novel guided attention mechanism and beam search applied to both training and inference. Compared to the state of the art our model does not require handcrafted features and is trained end-to-end. Furthermore, the proposed approach achieves an unmatched accuracy score of 85% on SwDA, and state-of-the-art accuracy score of 91.6% on MRDA.


title: A Multiclass Classification Approach to Label Ranking

abstract: In multiclass classification, the goal is to learn how to predict a random label $Y$, valued in $\mathcal{Y}=\{1,\; \ldots,\; K \}$ with $K\geq 3$, based upon observing a r.v. $X$, taking its values in $\mathbb{R}^q$ with $q\geq 1$ say, by means of a classification rule $g:\mathbb{R}^q\to \mathcal{Y}$ with minimum probability of error $\mathbb{P}\{Y\neq g(X) \}$. However, in a wide variety of situations, the task targeted may be more ambitious, consisting in sorting all the possible label values $y$ that may be assigned to $X$ by decreasing order of the posterior probability $\eta_y(X)=\mathbb{P}\{Y=y \mid X \}$. This article is devoted to the analysis of this statistical learning problem, halfway between multiclass classification and posterior probability estimation (regression) and referred to as label ranking here. We highlight the fact that it can be viewed as a specific variant of ranking median regression (RMR), where, rather than observing a random permutation $\Sigma$ assigned to the input vector $X$ and drawn from a Bradley-Terry-Luce-Plackett model with conditional preference vector $(\eta_1(X),\; \ldots,\; \eta_K(X))$, the sole information available for training a label ranking rule is the label $Y$ ranked on top, namely $\Sigma^{-1}(1)$. Inspired by recent results in RMR, we prove that under appropriate noise conditions, the One-Versus-One (OVO) approach to multiclassification yields, as a by-product, an optimal ranking of the labels with overwhelming probability. Beyond theoretical guarantees, the relevance of the approach to label ranking promoted in this article is supported by experimental results.


title: Recursive, parameter-free, explicitly defined interpolation nodes for simplices

abstract: A rule for constructing interpolation nodes for $n$th degree polynomials on the simplex is presented. These nodes are simple to define recursively from families of 1D node sets, such as the Lobatto-Gauss-Legendre (LGL) nodes. The resulting nodes have attractive properties: they are fully symmetric, they match the 1D family used in construction on the edges of the simplex, and the nodes constructed for the $(d-1)$-simplex are the boundary traces of the nodes constructed for the $d$-simplex. When compared using the Lebesgue constant to other explicit rules for defining interpolation nodes, the nodes recursively constructed from LGL nodes are nearly as good as the "warp & blend" nodes of Warburton in 2D (which, though defined differently, are very similar), and in 3D are better than other known explicit rules by increasing margins for $n > 6$. By that same measure, these recursively defined nodes are not as good as implicitly defined nodes found by optimizing the Lebesgue constant or related functions, but such optimal node sets have yet to be computed for the tetrahedron. A reference python implementation has been distributed as the `recursivenodes` package, but the simplicity of the recursive construction makes them easy to implement.


title: Robustness from Simple Classifiers

abstract: Despite the vast success of Deep Neural Networks in numerous application domains, it has been shown that such models are not robust i.e., they are vulnerable to small adversarial perturbations of the input. While extensive work has been done on why such perturbations occur or how to successfully defend against them, we still do not have a complete understanding of robustness. In this work, we investigate the connection between robustness and simplicity. We find that simpler classifiers, formed by reducing the number of output classes, are less susceptible to adversarial perturbations. Consequently, we demonstrate that decomposing a complex multiclass model into an aggregation of binary models enhances robustness. This behavior is consistent across different datasets and model architectures and can be combined with known defense techniques such as adversarial training. Moreover, we provide further evidence of a disconnect between standard and robust learning regimes. In particular, we show that elaborate label information can help standard accuracy but harm robustness.


title: Human Action Recognition using Local Two-Stream Convolution Neural Network Features and Support Vector Machines

abstract: This paper proposes a simple yet effective method for human action recognition in video. The proposed method separately extracts local appearance and motion features using state-of-the-art three-dimensional convolutional neural networks from sampled snippets of a video. These local features are then concatenated to form global representations which are then used to train a linear SVM to perform the action classification using full context of the video, as partial context as used in previous works. The videos undergo two simple proposed preprocessing techniques, optical flow scaling and crop filling. We perform an extensive evaluation on three common benchmark dataset to empirically show the benefit of the SVM, and the two preprocessing steps.


title: SummaryNet: A Multi-Stage Deep Learning Model for Automatic Video Summarisation

abstract: Video summarisation can be posed as the task of extracting important parts of a video in order to create an informative summary of what occurred in the video. In this paper we introduce SummaryNet as a supervised learning framework for automated video summarisation. SummaryNet employs a two-stream convolutional network to learn spatial (appearance) and temporal (motion) representations. It utilizes an encoder-decoder model to extract the most salient features from the learned video representations. Lastly, it uses a sigmoid regression network with bidirectional long short-term memory cells to predict the probability of a frame being a summary frame. Experimental results on benchmark datasets show that the proposed method achieves comparable or significantly better results than the state-of-the-art video summarisation methods.


title: Upset Recovery Control for Quadrotors Subjected to a Complete Rotor Failure from Large Initial Disturbances

abstract: This study has developed a fault-tolerant controller that is able to recover a quadrotor from arbitrary initial orientations and angular velocities, despite the complete failure of a rotor. This cascaded control method includes a position/altitude controller, an almost-global convergence attitude controller, and a control allocation method based on quadratic programming. As a major novelty, a constraint of undesirable angular velocity is derived and fused into the control allocator, which significantly improves the recovery performance. For validation, we have conducted a set of Monte-Carlo simulation to test the reliability of the proposed method of recovering the quadrotor from arbitrary initial attitude/rate conditions. In addition, real-life flight tests have been performed. The results demonstrate that the post-failure quadrotor can recover after being casually tossed into the air.


title: Whittle estimation for stationary state space models with finite second moments

abstract: In this paper, we consider the Whittle estimator for the parameters of a stationary solution of a continuous-time linear state space model sampled at low frequencies. In our context the driving process is a L\'evy process which allows flexible margins of the underlying model. The L\'evy process is supposed to have finite second moments. It is well known that then the class of stationary solutions of linear state space models and the class of multivariate CARMA processes coincides. We prove that the Whittle estimator, which is based on the periodogram, is strongly consistent and asymptotically normally distributed. A comparison with the classical setting of discrete-time ARMA models shows that in the continuous-time setting the limit covariance matrix of the Whittle estimator has an additional correction term for non-Gaussian models. For the proof, we investigate as well the asymptotic normality of the integrated periodogram which is interesting for its own. It can be used to construct goodness of fit tests. Furthermore, for univariate state space processes, which are CARMA processes, we introduce an adjusted version of the Whittle estimator and derive as well the asymptotic properties of this estimator. The practical applicability of our estimators is demonstrated through a simulation study.


title: Central limit theorems for Markov chains based on their convergence rates in Wasserstein distance

abstract: Many tools are available to bound the convergence rate of Markov chains in total variation (TV) distance. Such results can be used to establish central limit theorems (CLT) that enable error evaluations of Monte Carlo estimates in practice. However, convergence analysis based on TV distance is often non-scalable to high-dimensional Markov chains (Qin and Hobert (2018); Rajaratnam and Sparks (2015)). Alternatively, robust bounds in Wasserstein distance are often easier to obtain, thanks to a coupling argument. Our work is concerned with the implication of such convergence results, in particular, do they lead to CLTs of the corresponding Markov chains? One indirect and typically non-trivial way is to first convert Wasserstein bounds into total variation bounds. Alternatively, we provide two CLTs that directly depend on (sub-geometric) convergence rates in Wasserstein distance. Our CLTs hold for Lipschitz functions under certain moment conditions. Finally, we apply these CLTs to four sets of Markov chain examples including a class of nonlinear autoregressive processes, an exponential integrator version of the metropolis adjusted Langevin algorithm (EI-MALA), an unadjusted Langevin algorithm (ULA), and a special autoregressive model that generates reducible chains.


title: Effect of Transport Mechanism on the Coarsening of Bicontinuous Structures: A Comparison between Bulk and Surface Diffusion

abstract: Coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to $t^{1/4}$ for surface diffusion and $t^{1/3}$ for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.


title: Utilizing aLIGO Glitch Classifications to Validate Gravitational-Wave Candidates

abstract: Advanced LIGO data contains numerous noise transients, or "glitches", that have been shown to reduce the sensitivity of matched filter searches for gravitational waves from compact binaries by increasing the rate at which random coincidences occur. The presence of these transients has precipitated extensive work to establish that observed gravitational wave events are astrophysical in nature. We discuss the response of the PyCBC search for gravitational waves from stellar mass binaries to various common glitches that were observed during Advanced LIGO's first and second observing runs. We show how these transients can mimic waveforms from compact binary coalescences and quantify the likelihood that a given class of glitches will create a trigger in the search pipeline. We explore the specific waveform parameters that are most similar to different glitch classes and demonstrate how knowledge of these similarities can be used when evaluating the significance of gravitational-wave candidates.


title: Shared optical wireless cells for in-cabin aircraft links

abstract: The design of a wireless communication system that can support multiple users at high data rates inside an aircraft is a key requirement of aircraft manufacturers. This paper examines the design of an on-board visible light communication (VLC) system for transmitting data on board Boeing 747-400 aircraft. The reading light unit of each seat is utilised as an optical transmitter. A red, yellow, green, and blue (RYGB) laser diode (LD) is used in each reading light unit for transmitting data. An angle diversity receiver (ADR), which is an optical receiver that is composed of four branches (in this work), is evaluated. The signal-to-interference-plus-noise ratio (SINR) and data rate are determined. Three scenarios have been examined where, in the first scenario, one device is used, in the second scenario two devices are used and in the third scenario three devices are used by each passenger. The proposed system can offer high SINRs that support high data rates for each passenger by using simple on-off-keying (OOK).


title: Clusters of Galaxies Masquerading as X-Ray Quasars

abstract: Inspired by the discovery of the Phoenix cluster by the South Pole Telescope team, we initiated a search for other massive clusters of galaxies missing from the standard X-ray catalogs. We began by identifying 25 cluster candidates not included in the Meta-Catalog of X-ray Clusters of galaxies cluster compilation through cross-identification of the central galaxies of optically identified clusters in the Sloan Digital Sky Survey GMBCG catalog with bright X-ray sources in the ROSAT Bright Source Catalog. Those candidates were mostly unidentified or previously classified as X-ray active galactic nucleus (AGN). We analyzed brief Chandra X-ray Observatory observations of 14 of these X-ray sources and found that eight are X-ray luminous clusters of galaxies, only one showing evidence for a central X-ray point source. The remaining six candidates turned out to be point-source dominated, with faint detections or upper limits on any extended emission. We were not able to rule out the presence of extended X-ray emission from any of the point sources. The levels of extended emission around the six point sources are consistent with expectations based on optical richness, but could also be contaminated by scattered X-ray light from the central point source or extended nonthermal emission from possible radio lobes. We characterize the extended components of each of the well-detected cluster sources, finding that six of the eight X-ray clusters are consistent with being compact cool-core clusters. One of the newly identified low-luminosity X-ray clusters may have had an X-ray-luminous AGN 20 yr prior to the recent Chandra observations, based on the 4{\sigma} difference between its Chandra and ROSAT fluxes.


title: Inverse M-matrix, a new characterization

abstract: In this article we present a new characterization of inverse M-matrices, inverse row diagonally dominant M-matrices and inverse row and column diagonally dominant M-matrices, based on the positivity of certain inner products.


title: Embeddings using universal words in the free group of rank 2

abstract: For a countable group G = <A | R> presented by its generators A and defining relations R we discuss a simple method to embed G into such a 2-generator group T that the images of generators from A are explicitly given in T, and the defining relations for T can automatically be deduced from the relations R. The obtained method is applied on particular examples of groups, and references of its application for embeddings of recursive groups into finitely presented groups are given.


title: Few-Shot Learning via Learning the Representation, Provably

abstract: This paper studies few-shot learning via representation learning, where one uses $T$ source tasks with $n_1$ data per task to learn a representation in order to reduce the sample complexity of a target task for which there is only $n_2 (\ll n_1)$ data. Specifically, we focus on the setting where there exists a good \emph{common representation} between source and target, and our goal is to understand how much of a sample size reduction is possible. First, we study the setting where this common representation is low-dimensional and provide a fast rate of $O\left(\frac{\mathcal{C}\left(\Phi\right)}{n_1T} + \frac{k}{n_2}\right)$; here, $\Phi$ is the representation function class, $\mathcal{C}\left(\Phi\right)$ is its complexity measure, and $k$ is the dimension of the representation. When specialized to linear representation functions, this rate becomes $O\left(\frac{dk}{n_1T} + \frac{k}{n_2}\right)$ where $d (\gg k)$ is the ambient input dimension, which is a substantial improvement over the rate without using representation learning, i.e. over the rate of $O\left(\frac{d}{n_2}\right)$. This result bypasses the $\Omega(\frac{1}{T})$ barrier under the i.i.d. task assumption, and can capture the desired property that all $n_1T$ samples from source tasks can be \emph{pooled} together for representation learning. Next, we consider the setting where the common representation may be high-dimensional but is capacity-constrained (say in norm); here, we again demonstrate the advantage of representation learning in both high-dimensional linear regression and neural network learning. Our results demonstrate representation learning can fully utilize all $n_1T$ samples from source tasks.


title: Multi-Agent Learning in Double-side Auctions forPeer-to-peer Energy Trading

abstract: Distributed energy resources (DERs), such as rooftop solar panels, are growing rapidly and are reshaping power systems. To promote DERs, feed-in-tariff (FIT) is usually adopted by utilities to pay DER owners certain fixed rates for supplying energy to the grid. An alternative to FIT is a market based approach; i.e., consumers and DER owners trade energy in an auction-based peer-to-peer (P2P) market, and the rates are determined by a market clearing process. However, the complexities in sucha market and agents' bounded rationality may invalidate many well-established theories on auction design and hinder market development. To address this issue, we propose an automated bidding framework in a repeated auction based on multi-armed bandit learning, which aims to minimize each bidder's cumulative regret. Numerical results indicate convergence of such a multi-agent learning game to a steady-state. For comparison purpose, we apply the framework to three different auction designs to realize a P2P market.


title: Likelihood-free inference of experimental Neutrino Oscillations using Neural Spline Flows

abstract: In machine learning, likelihood-free inference refers to the task of performing an analysis driven by data instead of an analytical expression. We discuss the application of Neural Spline Flows, a neural density estimation algorithm, to the likelihood-free inference problem of the measurement of neutrino oscillation parameters in Long Baseline neutrino experiments. A method adapted to physics parameter inference is developed and applied to the case of the disappearance muon neutrino analysis at the T2K experiment.


title: Calibrating Deep Neural Networks using Focal Loss

abstract: Miscalibration - a mismatch between a model's confidence and its correctness - of Deep Neural Networks (DNNs) makes their predictions hard to rely on. Ideally, we want networks to be accurate, calibrated and confident. We show that, as opposed to the standard cross-entropy loss, focal loss [Lin et. al., 2017] allows us to learn models that are already very well calibrated. When combined with temperature scaling, whilst preserving accuracy, it yields state-of-the-art calibrated models. We provide a thorough analysis of the factors causing miscalibration, and use the insights we glean from this to justify the empirically excellent performance of focal loss. To facilitate the use of focal loss in practice, we also provide a principled approach to automatically select the hyperparameter involved in the loss function. We perform extensive experiments on a variety of computer vision and NLP datasets, and with a wide variety of network architectures, and show that our approach achieves state-of-the-art calibration without compromising on accuracy in almost all cases. Code is available at https://github.com/torrvision/focal_calibration.


title: Online Batch Decision-Making with High-Dimensional Covariates

abstract: We propose and investigate a class of new algorithms for sequential decision making that interacts with \textit{a batch of users} simultaneously instead of \textit{a user} at each decision epoch. This type of batch models is motivated by interactive marketing and clinical trial, where a group of people are treated simultaneously and the outcomes of the whole group are collected before the next stage of decision. In such a scenario, our goal is to allocate a batch of treatments to maximize treatment efficacy based on observed high-dimensional user covariates. We deliver a solution, named \textit{Teamwork LASSO Bandit algorithm}, that resolves a batch version of explore-exploit dilemma via switching between teamwork stage and selfish stage during the whole decision process. This is made possible based on statistical properties of LASSO estimate of treatment efficacy that adapts to a sequence of batch observations. In general, a rate of optimal allocation condition is proposed to delineate the exploration and exploitation trade-off on the data collection scheme, which is sufficient for LASSO to identify the optimal treatment for observed user covariates. An upper bound on expected cumulative regret of the proposed algorithm is provided.


title: Tearing and related field distortions in deep-depletion CCDs

abstract: Tearing patterns affecting flat field frames in CCDs are a visually striking obstacle to performing Pixel Response Non-Uniformity corrections. These patterns can be explained by lateral field distortions, caused by the non-uniform distribution of holes in the channel stops between sensor columns. Over the course of LSST camera development, a number of practical fixes have been suggested to get rid of tearing. But applying these fixes to our 16-channel Teledyne-e2v sensors leaves at best a distortion pattern at the vertical edges of every segment. Our working hypothesis is that the origin of the tearing is the parallel clocking itself, which moves the holes that are present in the channel stops regions. The efficiency of these transfers depends strongly on the details of the clocking operations, resulting in the observed variety of distortion patterns. Removal of most of the distortion patterns can therefore be achieved by executing a purge operation, which flattens back the hole distribution in the channel stops, immediately before acquiring a frame. A more effective solution is to switch all clocking operations to use a bipolar voltage set.


title: A Toolkit for Generating Code Knowledge Graphs

abstract: Knowledge graphs have been proven extremely useful in powering diverse applications in semantic search and natural language understanding. In this paper, we present GraphGen4Code, a toolkit to build code knowledge graphs that can similarly power various applications such as program search, code understanding, bug detection, and code automation. GraphGen4Code uses generic techniques to capture code semantics with the key nodes in the graph representing classes, functions, and methods. Edges indicate function usage (e.g., how data flows through function calls, as derived from program analysis of real code), and documentation about functions (e.g., code documentation, usage documentation, or forum discussions such as StackOverflow). Our toolkit uses named graphs in RDF to model graphs per program, or can output graphs as JSON. We show the scalability of the toolkit by applying it to 1.3 million Python files drawn from GitHub, 2,300 Python modules, and 47 million forum posts. This results in an integrated code graph with over 2 billion triples. We make the toolkit to build such graphs as well as the sample extraction of the 2 billion triples graph publicly available to the community for use.


title: Minimizing Localized Ratio Cut Objectives in Hypergraphs

abstract: Hypergraphs are a useful abstraction for modeling multiway relationships in data, and hypergraph clustering is the task of detecting groups of closely related nodes in such data. Graph clustering has been studied extensively, and there are numerous methods for detecting small, localized clusters without having to explore an entire input graph. However, there are only a few specialized approaches for localized clustering in hypergraphs. Here we present a framework for local hypergraph clustering based on minimizing localized ratio cut objectives. Our framework takes an input set of reference nodes in a hypergraph and solves a sequence of hypergraph minimum $s$-$t$ cut problems in order to identify a nearby well-connected cluster of nodes that overlaps substantially with the input set. Our methods extend graph-based techniques but are significantly more general and have new output quality guarantees. First, our methods can minimize new generalized notions of hypergraph cuts, which depend on specific configurations of nodes within each hyperedge, rather than just on the number of cut hyperedges. Second, our framework has several attractive theoretical properties in terms of output cluster quality. Most importantly, our algorithm is strongly-local, meaning that its runtime depends only on the size of the input set, and does not need to explore the entire hypergraph to find good local clusters. We use our methodology to effectively identify clusters in hypergraphs of real-world data with millions of nodes, millions of hyperedges, and large average hyperedge size with runtimes ranging between a few seconds and a few minutes.


title: 3He-rich solar energetic particles: Solar sources

abstract: 3He-rich solar energetic particles (SEPs), showing up to a 10,000-fold abundance enhancement of rare elements like 3He or ultra-heavy nuclei, have been a puzzle for more than 50 years. One reason for the current lack of understanding of 3He-rich SEPs is the difficulty resolving the source regions of these commonly occurring events. Since their discovery, there has been strong evidence that 3He-rich SEP production is associated with flares on the Sun. Anomalous abundances of 3He-rich SEPs have been attributed to a unique acceleration mechanism that must routinely operate at flare sites. Flares associated with 3He-rich SEPs have been often observed in jet-like forms indicating an acceleration in magnetic reconnection involving field lines open to interplanetary space. Owing to a fleet of spacecraft around the Sun, providing a greatly improved resolution of solar imaging observations, 3He-rich SEP sources are now explored in unprecedented detail. This paper outlines the current understanding of 3He-rich SEPs, mainly focusing on their solar sources.


title: A family of operators generating domino tableaux of a fixed shape and a decomposition of left cells into isotypic components

abstract: We exhibit a set of operators on pairs of domino tableaux of the same shape sending them to other such pairs with the same right tableau, in such a way that any two pairs with the same right tableau are conjugate by some composition of the operators. Using these operators we give explicit bases for the isotypic components of a classical Kazhdan-Lusztig left cell in terms of Kazhdan-Lusztig basis vectors.


title: A Novel Exact Analytical Expression for the Magnetic Field of a Solenoid

abstract: In this paper we present the analytical calculations to derive the magnetic field of a solenoid by solving exactly a fractional integral with the use of a novel method. Starting from the Biot-Savart law, we consider a coil of negligible thickness with a stationary electric current. We derive the expressions of the on and off-axes magnetic field components. The results have been compared to some simplified and known analytical formulae as well as to a commercial numerical code showing a good agreement.


title: Stability of the indirect utility process

abstract: We investigate the dynamic stability of the indirect utility process associated with a (possibly suboptimal) trading strategy under perturbations of the market. Establishing the reverse conjugacy characterizations first, we prove continuity and first-order convergence of the indirect-utility process under simultaneous perturbations of the finite variation and martingale parts of the return of the risky asset.


title: Unified description of dark energy and dark matter within the generalized hybrid metric-Palatini theory of gravity

abstract: The generalized hybrid metric-Palatini theory of gravity admits a scalar-tensor representation in terms of two interacting scalar fields. We show that, upon an appropriate choice of the interaction potential, one of the scalar fields behaves like dark energy, inducing a late-time accelerated expansion of the universe, while the other scalar field behaves like pressureless dark matter that, together with ordinary baryonic matter, dominates the intermediate phases of cosmic evolution. This unified description of dark energy and dark matter gives rise to viable cosmological solutions, which reproduce the main features of the evolution of the universe.


title: Quantum Coherence in Ergodic and Many-Body Localized Systems

abstract: Quantum coherence quantifies the amount of superposition a quantum state can have in a given basis. Since there is a difference in the structure of eigenstates of the ergodic and many-body localized systems, we expect them also to differ in terms of their coherences in a given basis. Here, we numerically calculate different measures of quantum coherence in the excited eigenstates of an interacting disordered Hamiltonian as a function of the disorder. We show that quantum coherence can be used as an order parameter to detect the well-studied ergodic to many-body-localized phase transition. We also perform quantum quench studies to distinguish the behavior of coherence in thermalized and localized phases. We then present a protocol to calculate measurement-based localizable coherence to investigate the thermal and many-body localized phases. The protocol allows one to investigate quantum correlations experimentally in a non-destructive way, in contrast to measures that require tracing out a subsystem, which always destroys coherence and correlation.


title: Hadamard parametrix of the Feynman Green's function of a five-dimensional charged scalar field

abstract: The Hadamard parametrix is a representation of the short-distance singularity structure of the Feynman Green's function for a quantum field on a curved space-time background. Subtracting these divergent terms regularizes the Feynman Green's function and enables the computation of renormalized expectation values of observables. We study the Hadamard parametrix for a charged, massive, complex scalar field in five space-time dimensions. Even in Minkowski space-time, it is not possible to write the Feynman Green's function for a charged scalar field exactly in closed form. We therefore present covariant Taylor series expansions for the biscalars arising in the Hadamard parametrix. On a general space-time background, we explicitly state the expansion coefficients up to the order required for the computation of the renormalized scalar field current. These coefficients become increasingly lengthy as the order of the expansion increases, so we give the higher-order terms required for the calculation of the renormalized stress-energy tensor in Minkowski space-time only.


title: Snel: SQL Native Execution for LLVM

abstract: Snel is a relational database engine featuring Just-In-Time (JIT) compilation of queries and columnar data representation. Snel is designed for fast on-line analytics by leveraging the LLVM compiler infrastructure. It also has custom special methods like resolving histograms as extensions to the SQL language. "Snel" means "SQL Native Execution for LLVM". Unlike traditional database engines, it does not provide a client-server interface. Instead, it exposes its interface as an extension to SQLite, for a simple interactive usage from command line and for embedding in applications. Since Snel tables are read-only, it does not provide features like transactions or updates. This allows queries to be very fast since they don't have the overhead of table locking or ensuring consistency. At its core, Snel is simply a dynamic library that can be used by client applications. It has an SQLite extension for seamless integration with a traditional SQL environment and simple interactive usage from command line.


title: Entire theta operators at unramified primes

abstract: Starting with work of Serre, Katz, and Swinnerton-Dyer, theta operators have played a key role in the study of $p$-adic and $\bmod p$ modular forms and Galois representations. This paper achieves two main results for theta operators on automorphic forms on PEL-type Shimura varieties: 1) the analytic continuation at unramified primes $p$ to the whole Shimura variety of the $\bmod p$ reduction of $p$-adic Maass--Shimura operators {\it a priori} defined only over the $\mu$-ordinary locus, and 2) the construction of new $\bmod p$ theta operators that do not arise as the $\bmod p$ reduction of Maass--Shimura operators. While the main accomplishments of this paper concern the geometry of Shimura varieties and consequences for differential operators, we conclude with applications to Galois representations. Our approach involves a careful analysis of the behavior of Shimura varieties and enables us to obtain more general results than allowed by prior techniques, including for arbitrary signature, vector weights, and unramified primes in CM fields of arbitrary degree.


title: On the Descriptive Complexity of Temporal Constraint Satisfaction Problems

abstract: Finite-domain constraint satisfaction problems are either solvable by Datalog, or not even expressible in fixed-point logic with counting. The border between the two regimes coincides with an important dichotomy in universal algebra; in particular, the border can be described by a strong height-one Maltsev condition. For infinite-domain CSPs, the situation is more complicated even if the template structure of the CSP is model-theoretically tame. We prove that there is no Maltsev condition that characterizes Datalog already for the CSPs of first-order reducts of (Q;<); such CSPs are called temporal CSPs and are of fundamental importance in infinite-domain constraint satisfaction. Our main result is a complete classification of temporal CSPs that can be expressed in one of the following logical formalisms: Datalog, fixed-point logic (with or without counting), or fixed-point logic with the Boolean rank operator. The classification shows that many of the equivalent conditions in the finite fail to capture expressibility in Datalog or fixed-point logic already for temporal CSPs.


title: Massive MIMO Channel Measurements and Achievable Rates in a Residential Area

abstract: In this paper we present a measurement set-up for massive MIMO channel sounding that shows very good long-term phase stability. Initial measurements were performed in a residential area to evaluate different conventional precoding schemes such as maximum ratio transmission and phase only precoding. A massive amount of data points was collected, with 924 times 64 complex channel weights per data point. Each data point is position-tagged using differential GPS with real-time kinematik, achieving better than 35cm position accuracy in more than 90% of the collected data points, making this dataset a rich resource for, e.g., further studying machine learning based, data-driven approaches in wireless communications.


title: BER Analysis of Reconfigurable Intelligent Surface Assisted Downlink Power Domain NOMA

abstract: The use of software controlled passive Reconfigurable Intelligent Surface (RIS) in wireless communications has attracted many researchers in recent years. RIS has a certain degree of control over the scattering and reflection characteristics of the electromagnetic waves, compared to the conventional communications in which the received signal is degraded due to the uncontrollable scattering of the transmitted signal and its interaction with the objects in propagating medium. Further, in RIS assisted communications, the phases of the multiple incoming signals can be controlled to enable constructive addition of multiple signals from different channel paths to improve Signal to Noise Ratio (SNR). On the other hand, Non-Orthogonal Multiple Access (NOMA) provides massive connectivity and low latency. The power domain variant NOMA uses superposition coded symbols with different powers for different user symbols. In this paper, a novel RIS assisted downlink NOMA system is proposed by combining the merits of both RIS and NOMA to improve the reliability of the system. Analytical expressions are derived for the Bit Error Rate (BER) performance of the proposed RIS assisted power domain NOMA system. The BER performance of the proposed system is analyzed using the numerical simulation results. It is observed that the proposed system has better performance than the conventional NOMA system.


title: On the Interaction of Galactic Cosmic Rays with Heliospheric Shocks During Forbush Decreases

abstract: Forbush decreases (FDs) are depletions in the galactic cosmic ray (GCR) count rate that last typically for about a week and can be caused by coronal mass ejections (CMEs) or corotating interacting regions (CIRs). Fast CMEs that drive shocks cause large FDs that often show a two-step decrease where the first step is attributed to the shock/sheath region, while the second step is attributed to the closed magnetic structure. Since the difference in size of shock and sheath region is significant, and since there are observed effects that can be related to shocks and not necessarily to the sheath region we expect that the physical mechanisms governing the interaction with GCRs in these two regions are different. We therefore aim to analyse interaction of GCRs with heliospheric shocks only. We approximate the shock by a structure where the magnetic field linearly changes with position within this structure. We assume protons of different energy, different pitch angle and different incoming direction. We also vary the shock parameters such as the magnetic field strength and orientation, as well as the shock thickness. The results demonstrate that protons with higher energies are less likely to be reflected. Also, thicker shocks and shocks with stronger field reflect protons more efficiently.


title: Hybrid Symbolic-Numeric Framework for Power System Modeling and Analysis

abstract: With the recent proliferation of open-source packages for computing, power system differential-algebraic equation (DAE) modeling and simulation are being revisited to reduce the programming efforts. Existing open-source tools require manual efforts to develop code for numerical equations, sparse Jacobians, and discontinuous components. This paper proposes a hybrid symbolic-numeric framework, exemplified by an open-source Python-based library ANDES, which consists of a symbolic layer for descriptive modeling and a numeric layer for vector-based numerical computation. This method enables the implementation of DAE models by mixing and matching modeling components, through which models are described. In the framework, a rich set of discontinuous components and standard transfer function blocks are provided besides essential modeling elements for rapid modeling. ANDES can automatically generate robust and fast numerical simulation code, as well as and high-quality documentation. Case studies present a) two implementations of turbine governor model TGOV1, b) power flow computation time break down for MATPOWER systems, c) validation of time-domain simulation with commercial software using three test systems with a variety of models, and d) the full eigenvalue analysis for Kundur's system. Validation shows that ANDES closely matches the commercial tool DSATools for power flow, time-domain simulation, and eigenvalue analysis.


title: Lower bounds for Cauchy data on curves in a negatively curved surface

abstract: We prove a uniform lower bound on Cauchy data on an arbitrary curve on a negatively curved surface using the Dyatlov-Jin(-Nonnenmacher) observability estimate on the global surface. In the process, we prove some further results about defect measures of restrictions of eigenfunctions to a hypersurface.


title: Tight paths in convex geometric hypergraphs

abstract: In this paper, we prove a theorem on tight paths in convex geometric hypergraphs, which is asymptotically sharp in infinitely many cases. Our geometric theorem is a common generalization of early results of Hopf and Pannwitz [12], Sutherland [19], Kupitz and Perles [16] for convex geometric graphs, as well as the classical Erd\H{o}s-Gallai Theorem [6] for graphs. As a consequence, we obtain the first substantial improvement on the Tur\'an problem for tight paths in uniform hypergraphs.


title: Sequential Submodular Maximization and Applications to Ranking an Assortment of Products

abstract: We study a submodular maximization problem motivated by applications in online retail. A platform displays a list of products to a user in response to a search query. The user inspects the first $k$ items in the list for a $k$ chosen at random from a given distribution, and decides whether to purchase an item from that set based on a choice model. The goal of the platform is to maximize the engagement of the shopper defined as the probability of purchase. This problem gives rise to a less-studied variation of submodular maximization in which we are asked to choose an $\textit{ordering}$ of a set of elements to maximize a linear combination of different submodular functions. First, using a reduction to maximizing submodular functions over matroids, we give an optimal $\left(1-1/e\right)$-approximation for this problem. We then consider a variant in which the platform cares not only about user engagement, but also about diversification across various groups of users, that is, guaranteeing a certain probability of purchase in each group. We characterize the polytope of feasible solutions and give a bi-criteria $((1-1/e)^2,(1-1/e)^2)$-approximation for this problem by rounding an approximate solution of a linear programming relaxation. For rounding, we rely on our reduction and the particular rounding techniques for matroid polytopes. For the special case in which underlying submodular functions are coverage functions -- which is practically relevant in online retail -- we propose an alternative LP relaxation and a simpler randomized rounding for the problem. This approach yields to an optimal bi-criteria $(1-1/e,1-1/e)$-approximation algorithm for the special case of the problem with coverage functions.


title: Hidden invariance of last passage percolation and directed polymers

abstract: Last passage percolation and directed polymer models on $\mathbb Z^2$ are invariant under translation and certain reflections. When these models have an integrable structure coming from either the RSK correspondence or the geometric RSK correspondence (e.g. geometric last passage percolation or the log-gamma polymer), we show that these basic invariances can be combined with a decoupling property to yield a rich new set of symmetries. Among other results, we prove shift and rearrangement invariance statements for last passage times, geodesic locations, disjointness probabilities, polymer partition functions, and quenched polymer measures. We also use our framework to find `scrambled' versions of the classical RSK correspondence, and to find an RSK correspondence for moon polyominoes. The results extend to limiting models, including the KPZ equation and the Airy sheet.


title: Parameterized Correlation Clustering in Hypergraphs and Bipartite Graphs

abstract: Motivated by applications in community detection and dense subgraph discovery, we consider new clustering objectives in hypergraphs and bipartite graphs. These objectives are parameterized by one or more resolution parameters in order to enable diverse knowledge discovery in complex data. For both hypergraph and bipartite objectives, we identify parameter regimes that are equivalent to existing objectives and share their (polynomial-time) approximation algorithms. We first show that our parameterized hypergraph correlation clustering objective is related to higher-order notions of normalized cut and modularity in hypergraphs. It is further amenable to approximation algorithms via hyperedge expansion techniques. Our parameterized bipartite correlation clustering objective generalizes standard unweighted bipartite correlation clustering, as well as bicluster deletion. For a certain choice of parameters it is also related to our hypergraph objective. Although in general it is NP-hard, we highlight a parameter regime for the bipartite objective where the problem reduces to the bipartite matching problem and thus can be solved in polynomial time. For other parameter settings, we present approximation algorithms using linear program rounding techniques. These results allow us to introduce the first constant-factor approximation for bicluster deletion, the task of removing a minimum number of edges to partition a bipartite graph into disjoint bi-cliques. In several experimental results, we highlight the flexibility of our framework and the diversity of results that can be obtained in different parameter settings. This includes clustering bipartite graphs across a range of parameters, detecting motif-rich clusters in an email network and a food web, and forming clusters of retail products in a product review hypergraph, that are highly correlated with known product categories.


title: Fine-Grained Instance-Level Sketch-Based Video Retrieval

abstract: Existing sketch-analysis work studies sketches depicting static objects or scenes. In this work, we propose a novel cross-modal retrieval problem of fine-grained instance-level sketch-based video retrieval (FG-SBVR), where a sketch sequence is used as a query to retrieve a specific target video instance. Compared with sketch-based still image retrieval, and coarse-grained category-level video retrieval, this is more challenging as both visual appearance and motion need to be simultaneously matched at a fine-grained level. We contribute the first FG-SBVR dataset with rich annotations. We then introduce a novel multi-stream multi-modality deep network to perform FG-SBVR under both strong and weakly supervised settings. The key component of the network is a relation module, designed to prevent model over-fitting given scarce training data. We show that this model significantly outperforms a number of existing state-of-the-art models designed for video analysis.


title: CP odd weak basis invariants in minimal see-saw model and Leptogenesis

abstract: In this paper, we derive the relationship between the weak basis invariants (WB) related to CP violation responsible for leptogenesis and CP violation relevant at low energy. We examine all the experimental viable cases of Frampton-Glashow and Yanagida (FGY) model, in order to construct the WB invariants in terms of left handed Majorana neutrino mass matrix elements, and thus finding the necessary and sufficient condition for CP conservation. Further for all the viable FGY texture zeros, we have shown the explicit dependence of WB invariants on Dirac type and Majorana type CP violating phases. In the end, we discuss the implication of such interrelationships on leptogenesis.


title: Privately Learning Markov Random Fields

abstract: We consider the problem of learning Markov Random Fields (including the prototypical example, the Ising model) under the constraint of differential privacy. Our learning goals include both structure learning, where we try to estimate the underlying graph structure of the model, as well as the harder goal of parameter learning, in which we additionally estimate the parameter on each edge. We provide algorithms and lower bounds for both problems under a variety of privacy constraints -- namely pure, concentrated, and approximate differential privacy. While non-privately, both learning goals enjoy roughly the same complexity, we show that this is not the case under differential privacy. In particular, only structure learning under approximate differential privacy maintains the non-private logarithmic dependence on the dimensionality of the data, while a change in either the learning goal or the privacy notion would necessitate a polynomial dependence. As a result, we show that the privacy constraint imposes a strong separation between these two learning problems in the high-dimensional data regime.


title: Private Mean Estimation of Heavy-Tailed Distributions

abstract: We give new upper and lower bounds on the minimax sample complexity of differentially private mean estimation of distributions with bounded $k$-th moments. Roughly speaking, in the univariate case, we show that $n = \Theta\left(\frac{1}{\alpha^2} + \frac{1}{\alpha^{\frac{k}{k-1}}\varepsilon}\right)$ samples are necessary and sufficient to estimate the mean to $\alpha$-accuracy under $\varepsilon$-differential privacy, or any of its common relaxations. This result demonstrates a qualitatively different behavior compared to estimation absent privacy constraints, for which the sample complexity is identical for all $k \geq 2$. We also give algorithms for the multivariate setting whose sample complexity is a factor of $O(d)$ larger than the univariate case.


title: Locally Private Hypothesis Selection

abstract: We initiate the study of hypothesis selection under local differential privacy. Given samples from an unknown probability distribution $p$ and a set of $k$ probability distributions $\mathcal{Q}$, we aim to output, under the constraints of $\varepsilon$-local differential privacy, a distribution from $\mathcal{Q}$ whose total variation distance to $p$ is comparable to the best such distribution. This is a generalization of the classic problem of $k$-wise simple hypothesis testing, which corresponds to when $p \in \mathcal{Q}$, and we wish to identify $p$. Absent privacy constraints, this problem requires $O(\log k)$ samples from $p$, and it was recently shown that the same complexity is achievable under (central) differential privacy. However, the naive approach to this problem under local differential privacy would require $\tilde O(k^2)$ samples. We first show that the constraint of local differential privacy incurs an exponential increase in cost: any algorithm for this problem requires at least $\Omega(k)$ samples. Second, for the special case of $k$-wise simple hypothesis testing, we provide a non-interactive algorithm which nearly matches this bound, requiring $\tilde O(k)$ samples. Finally, we provide sequentially interactive algorithms for the general case, requiring $\tilde O(k)$ samples and only $O(\log \log k)$ rounds of interactivity. Our algorithms are achieved through a reduction to maximum selection with adversarial comparators, a problem of independent interest for which we initiate study in the parallel setting. For this problem, we provide a family of algorithms for each number of allowed rounds of interaction $t$, as well as lower bounds showing that they are near-optimal for every $t$. Notably, our algorithms result in exponential improvements on the round complexity of previous methods.


title: Averages of long Dirichlet polynomials

abstract: We consider the asymptotic behavior of the mean square of truncations of the Dirichlet series of $\zeta(s)^k$. We discuss the connections of this problem with that of the variance of the divisor function in short intervals and in arithmetic progressions, reviewing the recent results on this topic. Finally, we show how these results can all be proved assuming a suitable version of the moments conjecture.


title: Hydrodynamics of a dense flock of sheep: edge motion and long-range correlations

abstract: Sheep are gregarious animals, and they often aggregate into dense, cohesive flocks, especially under stress. In this paper, we use image processing tools to analyze a publicly available aerial video showing a dense sheep flock moving under the stimulus of a shepherding dog. Inspired by the fluidity of the motion, we implement a hydrodynamics approach, extracting velocity fields, and measuring their propagation and correlations in space and time. We find that while the flock overall is stationary, significant dynamics happens at the edges, notably in the form of fluctuations propagating like waves, and large-scale correlations spanning the entire flock. These observations highlight the importance of incorporating interfacial dynamics, for instance in the form of line tension, when using a hydrodynamics framework to model the dynamics of dense, non-polarized swarms.


title: Solar Contamination in Extreme Precision Radial Velocity Measurements: Deleterious Effects and Prospects for Mitigation

abstract: Solar contamination, due to moonlight and atmospheric scattering of sunlight, can cause systematic errors in stellar radial velocity (RV) measurements that significantly detract from the ~10cm/s sensitivity required for the detection and characterization of terrestrial exoplanets in or near Habitable Zones of Sun-like stars. The addition of low-level spectral contamination at variable effective velocity offsets introduces systematic noise when measuring velocities using classical mask-based or template-based cross-correlation techniques. Here we present simulations estimating the range of RV measurement error induced by uncorrected scattered sunlight contamination. We explore potential correction techniques, using both simultaneous spectrometer sky fibers and broadband imaging via coherent fiber imaging bundles, that could reliably reduce this source of error to below the photon-noise limit of typical stellar observations. We discuss the limitations of these simulations, the underlying assumptions, and mitigation mechanisms. We also present and discuss the components designed and built into the NEID precision RV instrument for the WIYN 3.5m telescope, to serve as an ongoing resource for the community to explore and evaluate correction techniques. We emphasize that while "bright time" has been traditionally adequate for RV science, the goal of 10cm/s precision on the most interesting exoplanetary systems may necessitate access to darker skies for these next-generation instruments.


title: An end-to-end approach for the verification problem: learning the right distance

abstract: In this contribution, we augment the metric learning setting by introducing a parametric pseudo-distance, trained jointly with the encoder. Several interpretations are thus drawn for the learned distance-like model's output. We first show it approximates a likelihood ratio which can be used for hypothesis tests, and that it further induces a large divergence across the joint distributions of pairs of examples from the same and from different classes. Evaluation is performed under the verification setting consisting of determining whether sets of examples belong to the same class, even if such classes are novel and were never presented to the model during training. Empirical evaluation shows such method defines an end-to-end approach for the verification problem, able to attain better performance than simple scorers such as those based on cosine similarity and further outperforming widely used downstream classifiers. We further observe training is much simplified under the proposed approach compared to metric learning with actual distances, requiring no complex scheme to harvest pairs of examples.


title: Score-based likelihood ratios to evaluate forensic pattern evidence

abstract: In 2016, the European Network of Forensic Science Institutes (ENFSI) published guidelines for the evaluation, interpretation and reporting of scientific evidence. In the guidelines, ENFSI endorsed the use of the likelihood ratio (LR) as a means to represent the probative value of most types of evidence. While computing the value of a LR is practical in several forensic disciplines, calculating an LR for pattern evidence such as fingerprints, firearm and other toolmarks is particularly challenging because standard statistical approaches are not applicable. Recent research suggests that machine learning algorithms can summarize a potentially large set of features into a single score which can then be used to quantify the similarity between pattern samples. It is then possible to compute a score-based likelihood ratio (SLR) and obtain an approximation to the value of the evidence, but research has shown that the SLR can be quite different from the LR not only in size but also in direction. We provide theoretical and empirical arguments that under reasonable assumptions, the SLR can be a practical tool for forensic evaluations.


title: Learning Fairness-aware Relational Structures

abstract: The development of fair machine learning models that effectively avert bias and discrimination is an important problem that has garnered attention in recent years. The necessity of encoding complex relational dependencies among the features and variables for competent predictions require the development of fair, yet expressive relational models. In this work, we introduce Fair-A3SL, a fairness-aware structure learning algorithm for learning relational structures, which incorporates fairness measures while learning relational graphical model structures. Our approach is versatile in being able to encode a wide range of fairness metrics such as statistical parity difference, overestimation, equalized odds, and equal opportunity, including recently proposed relational fairness measures. While existing approaches employ the fairness measures on pre-determined model structures post prediction, Fair-A3SL directly learns the structure while optimizing for the fairness measures and hence is able to remove any structural bias in the model. We demonstrate the effectiveness of our learned model structures when compared with the state-of-the-art fairness models quantitatively and qualitatively on datasets representing three different modeling scenarios: i) a relational dataset, ii) a recidivism prediction dataset widely used in studying discrimination, and iii) a recommender systems dataset. Our results show that Fair-A3SL can learn fair, yet interpretable and expressive structures capable of making accurate predictions.


title: Geometric rank of tensors and subrank of matrix multiplication

abstract: Motivated by problems in algebraic complexity theory (e.g., matrix multiplication) and extremal combinatorics (e.g., the cap set problem and the sunflower problem), we introduce the geometric rank as a new tool in the study of tensors and hypergraphs. We prove that the geometric rank is an upper bound on the subrank of tensors and the independence number of hypergraphs. We prove that the geometric rank is smaller than the slice rank of Tao, and relate geometric rank to the analytic rank of Gowers and Wolf in an asymptotic fashion. As a first application, we use geometric rank to prove a tight upper bound on the (border) subrank of the matrix multiplication tensors, matching Strassen's well-known lower bound from 1987.


title: Clustering as an Evaluation Protocol for Knowledge Embedding Representation of Categorised Multi-relational Data in the Clinical Domain

abstract: Learning knowledge representation is an increasingly important technology applicable in many domain-specific machine learning problems. We discuss the effectiveness of traditional Link Prediction or Knowledge Graph Completion evaluation protocol when embedding knowledge representation for categorised multi-relational data in the clinical domain. Link prediction uses to split the data into training and evaluation subsets, leading to loss of information along training and harming the knowledge representation model accuracy. We propose a Clustering Evaluation Protocol as a replacement alternative to the traditionally used evaluation tasks. We used embedding models trained by a knowledge embedding approach which has been evaluated with clinical datasets. Experimental results with Pearson and Spearman correlations show strong evidence that the novel proposed evaluation protocol is pottentially able to replace link prediction.


title: Fast Implementation of Morphological Filtering Using ARM NEON Extension

abstract: In this paper we consider speedup potential of morphological image filtering on ARM processors. Morphological operations are widely used in image analysis and recognition and their speedup in some cases can significantly reduce overall execution time of recognition. More specifically, we propose fast implementation of erosion and dilation using ARM SIMD extension NEON. These operations with the rectangular structuring element are separable. They were implemented using the advantages of separability as sequential horizontal and vertical passes. Each pass was implemented using van Herk/Gil-Werman algorithm for large windows and low-constant linear complexity algorithm for small windows. Final implementation was improved with SIMD and used a combination of these methods. We also considered fast transpose implementation of 8x8 and 16x16 matrices using ARM NEON to get additional computational gain for morphological operations. Experiments showed 3 times efficiency increase for final implementation of erosion and dilation compared to van Herk/Gil-Werman algorithm without SIMD, 5.7 times speedup for 8x8 matrix transpose and 12 times speedup for 16x16 matrix transpose compared to transpose without SIMD.


title: Dirty higher-order Dirac semimetal: Quantum criticality and bulk-boundary correspondence

abstract: We analyze the stability of time-reversal (${\mathcal T}$) and lattice four-fold ($C_4$) symmetry breaking three-dimensional higher-order topological (HOT) Dirac semimetals (DSMs) and the associated one-dimensional hinge modes in the presence of random pointlike charge impurities. Complementary real space numerical and momentum space renormalization group (RG) analyses suggest that a HOTDSM, while being a stable phase of matter for sufficiently weak disorder, undergoes a continuous quantum phase transition into a trivial metal at finite disorder. However, the corresponding critical exponents (numerically obtained from the scaling of the density of states) are extremely close to the ones found in a dirty, but first-order DSM that on the other hand preserves ${\mathcal T}$ and $C_4$ symmetries, and support two Fermi arc surface states. This observation suggests an emergent \emph{superuniversality} (insensitive to symmetries) in the entire family of dirty DSMs, as also predicted by a leading-order RG analysis. As a direct consequence of the bulk-boundary correspondence, the hinge modes in a system with open boundaries gradually fade away with increasing randomness, and completely dissolve in the trivial metallic phase at strong disorder.


title: Hydrogen Peroxide Formation Rates in a PEMFC Anode and Cathode: Effect of Humidity and Temperature

abstract: Hydrogen peroxide (H2O2) formation rates in a proton exchange membrane (PEM) fuel cell anode and cathode were estimated as a function of humidity and temperature by studying the oxygen reduction reaction (ORR) on a rotating ring disc electrode (RRDE). Fuel cell conditions were replicated by depositing a film of Pt/Vulcan XC-72 catalyst onto the disk and by varying the temperature, dissolved O2 concentration and the acidity levels in hydrochloric acid (HClO4). The HClO4 acidity was correlated to ionomer water activity and hence fuel cell humidity. The H2O2 formation rates showed a linear dependence on oxygen concentration and square dependence on water activity. The H2O2 selectivity in ORR was independent of oxygen concentration but increased with decrease in water activity (i.e., decreased humidity). Potential dependent activation energy for the H2O2 formation reaction was estimated from data obtained at different temperatures.


title: Graph Computing based Distributed State Estimation with PMUs

abstract: Power system state estimation plays a fundamental and critical role in the energy management system (EMS). To achieve a high performance and accurate system states estimation, a graph computing based distributed state estimation approach is proposed in this paper. Firstly, a power system network is divided into multiple areas. Reference buses are selected with PMUs being installed at these buses for each area. Then, the system network is converted into multiple independent areas. In this way, the power system state estimation could be conducted in parallel for each area and the estimated system states are obtained without compromise of accuracy. IEEE 118-bus system and MP 10790-bus system are employed to verify the results accuracy and present the promising computation performance.


title: On the Search for Feedback in Reinforcement Learning

abstract: The problem of Reinforcement Learning (RL) in an unknown nonlinear dynamical system is equivalent to the search for an optimal feedback law utilizing the simulations/ rollouts of the dynamical system. Most RL techniques search over a complex global nonlinear feedback parametrization making them suffer from high training times as well as variance. Instead, we advocate searching over a local feedback representation consisting of an open-loop sequence, and an associated optimal linear feedback law completely determined by the open-loop. We show that this alternate approach results in highly efficient training, the answers obtained are repeatable and hence reliable, and the resulting closed performance is superior to global state-of-the-art RL techniques. Finally, if we replan, whenever required, which is feasible due to the fast and reliable local solution, it allows us to recover global optimality of the resulting feedback law.