title: The Variable Density Model for the Rayleigh-Taylor Instability and its transformation to the diffusive, inhomogeneous, incompressible Navier-Stokes equations

abstract: It is shown how the variable density model (VDM) that governs the Rayleigh-Taylor instability (RTI) for the miscible mixing of two incompressible fluids can be transformed into a diffusive version of the inhomogeneous, incompressible Navier-Stokes equations forced by gradients of the composition density $\rho$ of the mixing layer. This demonstrates how buoyancy-driven flows drive and enhance Navier-Stokes turbulence. The role of the potential vorticity $q = \omega\cdot\nabla\rho$ is also discussed.


title: The role of the diffuse ionized gas in metallicity calibrations

abstract: Estimates of gas-phase abundances based on strong-line methods have been calibrated for H~{\scshape ii} regions. Those methods ignore any contribution from the diffuse ionized gas (DIG), which shows enhanced collisional-to-recombination line ratios in comparison to H~{\scshape ii} regions of the same metallicity. Applying strong line methods whilst ignoring the role of the DIG thus systematically overestimates metallicities. Using integral field spectroscopy data, we show how to correct for the DIG contribution and how it biases the mass--metallicity--star formation rate relation.


title: Evaluating the impact of increasing temperatures on changes in Soil Organic Carbon stocks: sensitivity analysis and non-standard discrete approximation

abstract: A novel model is here introduced for the SOC change index defined as the normalized difference between the actual Soil Organic Carbon and the value assumed at an initial reference year. It is tailored on the RothC carbon model dynamics and assumes as baseline the value of the SOC equilibrium under constant environmental conditions. A sensitivity analysis is performed to evaluate the response of the model to changes of temperature, Net Primary Production (NPP), and land use soil class (forest, grassland, arable). A non-standard monthly time-stepping procedure has been proposed to approximate the SOC change index in the Alta Murgia National Park, a protected area in the Italian Apulia region, selected as test site. In the case of arable class, the SOC change index exhibits a negative trend which can be inverted by a suitable organic fertilization program here proposed.


title: Technical Report for HW2VEC -- A Graph Learning Tool for Automating Hardware Security

abstract: In this technical report, we present HW2VEC [11], an open-source graph learning tool for hardware security, and its implementation details (Figure 1). HW2VEC provides toolboxes for graph representation extraction in the form of Data Flow Graphs (DFGs) or Abstract Syntax Trees (ASTs) from hardware designs at RTL and GLN levels. Besides, HW2VEC also offers graph learning tools for representing hardware designs in vectors that preserve both structural features and behavioral features. To the best of our knowledge, HW2VEC is the first open-source research tool that supports applying graph learning methods to hardware designs in different abstraction levels for hardware security. We organize the remainder of this technical report as follows: Section 2 introduces the architecture of HW2VEC; Section 3 gives information about the use-case implementations; Section 4 provides the experimental results and demonstrates the performance of HW2VEC for two hardware security applications: HT detection and IP piracy detection; finally, Section 5 will conclude this report.


title: Zooming Into the Darknet: Characterizing Internet Background Radiation and its Structural Changes

abstract: Network telescopes or "Darknets" provide a unique window into Internet-wide malicious activities associated with malware propagation, denial of service attacks, scanning performed for network reconnaissance, and others. Analyses of the resulting data can provide actionable insights to security analysts that can be used to prevent or mitigate cyber-threats. Large Darknets, however, observe millions of nefarious events on a daily basis which makes the transformation of the captured information into meaningful insights challenging. We present a novel framework for characterizing Darknet behavior and its temporal evolution aiming to address this challenge. The proposed framework: (i) Extracts a high dimensional representation of Darknet events composed of features distilled from Darknet data and other external sources; (ii) Learns, in an unsupervised fashion, an information-preserving low-dimensional representation of these events (using deep representation learning) that is amenable to clustering; (iv) Performs clustering of the scanner data in the resulting representation space and provides interpretable insights using optimal decision trees; and (v) Utilizes the clustering outcomes as "signatures" that can be used to detect structural changes in the Darknet activities. We evaluate the proposed system on a large operational Network Telescope and demonstrate its ability to detect real-world, high-impact cybersecurity incidents.


title: A New Semi-supervised Learning Benchmark for Classifying View and Diagnosing Aortic Stenosis from Echocardiograms

abstract: Semi-supervised image classification has shown substantial progress in learning from limited labeled data, but recent advances remain largely untested for clinical applications. Motivated by the urgent need to improve timely diagnosis of life-threatening heart conditions, especially aortic stenosis, we develop a benchmark dataset to assess semi-supervised approaches to two tasks relevant to cardiac ultrasound (echocardiogram) interpretation: view classification and disease severity classification. We find that a state-of-the-art method called MixMatch achieves promising gains in heldout accuracy on both tasks, learning from a large volume of truly unlabeled images as well as a labeled set collected at great expense to achieve better performance than is possible with the labeled set alone. We further pursue patient-level diagnosis prediction, which requires aggregating across hundreds of images of diverse view types, most of which are irrelevant, to make a coherent prediction. The best patient-level performance is achieved by new methods that prioritize diagnosis predictions from images that are predicted to be clinically-relevant views and transfer knowledge from the view task to the diagnosis task. We hope our released Tufts Medical Echocardiogram Dataset and evaluation framework inspire further improvements in multi-task semi-supervised learning for clinical applications.


title: Computational Complexity of Synchronization under Sparse Regular Constraints

abstract: The constrained synchronization problem (CSP) asks for a synchronizing word of a given input automaton contained in a regular set of constraints. It could be viewed as a special case of synchronization of a discrete event system under supervisory control. Here, we study the computational complexity of this problem for the class of sparse regular constraint languages. We give a new characterization of sparse regular sets, which equal the bounded regular sets, and derive a full classification of the computational complexity of CSP for letter-bounded regular constraint languages, which properly contain the strictly bounded regular languages. Then, we introduce strongly self-synchronizing codes and investigate CSP for bounded languages induced by these codes. With our previous result, we deduce a full classification for these languages as well. In both cases, depending on the constraint language, our problem becomes NP-complete or polynomial time solvable.


title: Towards Continual Entity Learning in Language Models for Conversational Agents

abstract: Neural language models (LM) trained on diverse corpora are known to work well on previously seen entities, however, updating these models with dynamically changing entities such as place names, song titles and shopping items requires re-training from scratch and collecting full sentences containing these entities. We aim to address this issue, by introducing entity-aware language models (EALM), where we integrate entity models trained on catalogues of entities into the pre-trained LMs. Our combined language model adaptively adds information from the entity models into the pre-trained LM depending on the sentence context. Our entity models can be updated independently of the pre-trained LM, enabling us to influence the distribution of entities output by the final LM, without any further training of the pre-trained LM. We show significant perplexity improvements on task-oriented dialogue datasets, especially on long-tailed utterances, with an ability to continually adapt to new entities (to an extent).


title: Majorization Minimization Methods for Distributed Pose Graph Optimization

abstract: We consider the problem of distributed pose graph optimization (PGO) that has important applications in multi-robot simultaneous localization and mapping (SLAM). We propose the majorization minimization (MM) method for distributed PGO ($\mathsf{MM-PGO}$) that applies to a broad class of robust loss kernels. The $\mathsf{MM-PGO}$ method is guaranteed to converge to first-order critical points under mild conditions. Furthermore, noting that the $\mathsf{MM-PGO}$ method is reminiscent of proximal methods, we leverage Nesterov's method and adopt adaptive restarts to accelerate convergence. The resulting accelerated MM methods for distributed PGO -- both with a master node in the network ($\mathsf{AMM-PGO}^*$) and without ($\mathsf{AMM-PGO}^{\#}$) -- have faster convergence in contrast to the $\mathsf{AMM-PGO}$ method without sacrificing theoretical guarantees. In particular, the $\mathsf{AMM-PGO}^{\#}$ method, which needs no master node and is fully decentralized, features a novel adaptive restart scheme and has a rate of convergence comparable to that of the $\mathsf{AMM-PGO}^*$ method using a master node to aggregate information from all the other nodes. The efficacy of this work is validated through extensive applications to 2D and 3D SLAM benchmark datasets and comprehensive comparisons against existing state-of-the-art methods, indicating that our MM methods converge faster and result in better solutions to distributed PGO.


title: The History of Speech Recognition to the Year 2030

abstract: The decade from 2010 to 2020 saw remarkable improvements in automatic speech recognition. Many people now use speech recognition on a daily basis, for example to perform voice search queries, send text messages, and interact with voice assistants like Amazon Alexa and Siri by Apple. Before 2010 most people rarely used speech recognition. Given the remarkable changes in the state of speech recognition over the previous decade, what can we expect over the coming decade? I attempt to forecast the state of speech recognition research and applications by the year 2030. While the changes to general speech recognition accuracy will not be as dramatic as in the previous decade, I suggest we have an exciting decade of progress in speech technology ahead of us.


title: Transverse Single Spin Asymmetries of Heavy Flavor Electrons and Charged Pions in 200 GeV $p+p^{\uparrow}$ Collisions at Midrapidity

abstract: Transverse single spin asymmetries of particles produced in $p+p^{\uparrow}$ collisions provide insight on the partonic spin and momentum structure of hadrons; heavy flavor electrons provide access to initial state spin-momentum correlations of gluons in the proton, while charged pions provide access to initial and final state transverse spin effects. Charged particles are measured at midrapidity at PHENIX using the silicon vertex detector and central arm spectrometer, made of an electromagnetic calorimeter, a ring-imaging Cherenkov detector, and drift and pad chambers. Recent results for both electron and charged pion measurements from the 2015 running period will be presented.


title: A generalized mean-field game model for the dynamics of pedestrians with limited predictive abilities

abstract: This paper investigates the model for pedestrian flow firstly proposed in [Cristiani et al., DOI:10.1137/140962413]. The model assumes that each individual in the crowd moves in a known domain, aiming at minimizing a given cost functional. Both the pedestrian dynamics and the cost functional itself depend on the position of the whole crowd. In addition, pedestrians are assumed to have predictive abilities, but limited in time, extending only up to $\theta$ time units into the future, where $\theta\in[0,\infty)$ is a model parameter. 1) For $\theta=0$ (no predictive abilities), we recover the modeling assumptions of the Hughes's model, where people take decisions on the basis of the current position of the crowd only. 2) For $\theta\to\infty$, instead, we recover the standard mean field game (MFG) setting, where people are able to forecast the behavior of the others at any future time and take decisions on the basis of the current and future position of the whole crowd. 3) For intermediate values of $\theta$ we obtain something different: as in the Hughes's model, the numerical procedure to solve the problem requires to run an off-line procedure at any fixed time $t$, which returns the current optimal velocity field at time $t$ by solving an associated backward-in-time Hamilton-Jacobi-Bellman equation; but, differently from the Hughes's model, here the procedure involves a prediction of the crowd behavior in the sliding time window $[t,t+\theta)$, therefore the optimal velocity field is given by the solution to a forward-backward system which joins a Fokker-Planck equation with a Hamilton-Jacobi-Bellman equation as in the MFG approach. The fact that a different forward-backward system must be solved at any time $t$ arises new interesting theoretical questions. Numerical tests will give some clues about the well-posedness of the problem.


title: Fluctuation and dissipation in memoryless open quantum evolutions

abstract: Von Neumann entropy rate for open quantum systems is, in general, written in terms of entropy production and entropy flow rates, encompassing the second law of thermodynamics. When the open-quantum-system evolution corresponds to a quantum dynamical semigroup, we find a decomposition of the infinitesimal generator of the dynamics, that allows to relate the von Neumann entropy rate with the divergence-based quantum Fisher information, at any time. Applied to quantum Gaussian channels that are dynamical semigroups, our decomposition leads to the quantum analog of the generalized classical de Bruijn identity, thus expressing the quantum fluctuation-dissipation relation in that kind of channels. Finally, from this perspective, we analyze how stationarity arises.


title: A Convergent Numerical Method for the Reflector Antenna Problem via Optimal Transport on the Sphere

abstract: We consider a PDE approach to numerically solving the reflector antenna problem by solving an Optimal Transport problem on the unit sphere with cost function $c(x,y) = -2\log \left\Vert x - y \right\Vert$. At each point on the sphere, we replace the surface PDE with a generalized Monge-Amp\`ere type equation posed on the local tangent plane. We then utilize a provably convergent finite difference scheme to approximate the solution and construct the reflector. The method is easily adapted to take into account highly nonsmooth data and solutions, which makes it particularly well adapted to real-world optics problems. Computational examples demonstrate the success of this method in computing reflectors for a range of challenging problems including discontinuous intensities and intensities supported on complicated geometries.


title: Tensor-Train Density Estimation

abstract: Estimation of probability density function from samples is one of the central problems in statistics and machine learning. Modern neural network-based models can learn high dimensional distributions but have problems with hyperparameter selection and are often prone to instabilities during training and inference. We propose a new efficient tensor train-based model for density estimation (TTDE). Such density parametrization allows exact sampling, calculation of cumulative and marginal density functions, and partition function. It also has very intuitive hyperparameters. We develop an efficient non-adversarial training procedure for TTDE based on the Riemannian optimization. Experimental results demonstrate the competitive performance of the proposed method in density estimation and sampling tasks, while TTDE significantly outperforms competitors in training speed.


title: Reactive Synthesis Modulo Theories Using Abstraction Refinement

abstract: Reactive synthesis builds a system from a specification given as a temporal logic formula. Traditionally, reactive synthesis is defined for systems with Boolean input and output variables. Recently, new theories and techniques have been proposed to extend reactive synthesis to data domains, which are required for more sophisticated programs. In particular, Temporal stream logic(TSL) (Finkbeiner et al. 2019) extends LTL with state variables, updates, and uninterpreted functions and was created for use in synthesis. We present a synthesis procedure for TSL(T), an extension of TSL with theories. Synthesis is performed using a counter-example guided synthesis loop and an LTL synthesis procedure. Our method translates TSL(T) specifications to LTL and extracts a system if synthesis is successful. Otherwise, it analyzes the counterstrategy for inconsistencies with the theory. If the counterstrategy is theory-consistent, it proves that the specification is unrealizable. Otherwise, we add temporal assumptions and Boolean predicates to the TSL(T) specification and start the next iteration of the the loop. We show that the synthesis problem for TSL (T) is undecidable. Nevertheless our method can successfully synthesize or show unrealizability of several non-Boolean examples.


title: Effect of intravalley and intervalley electron-hole exchange on the nonlinear optical response of monolayer MoSe$_2$

abstract: The coherent third-order nonlinear response of monolayer transition-metal dichalcogenide (TMD) semiconductors, such as MoSe$_2$ is dominated by the nonlinear exciton response, as well as biexciton and trion resonances. The fact that these resonances may be spectrally close together makes identification of the signatures, for example in differential transmission (DT), challenging. Instead of focusing on explaining a given set of experimental data, a systematic study aimed at elucidating the roles of intravalley and intervalley electron-hole (e-h) exchange on the DT spectra is presented. Previous works have shown that the e-h exchange introduces a linear leading-order term in the exciton dispersion. Based on a generalized Lippmann-Schwinger equation, we show that the presence of this linear dispersion term can reduce the biexciton binding energy to zero, contrary to the conventional situation of quadratic dispersion where an arbitrarily weak (well-behaved) attractive interaction always supports bound state(s). The effects of spin-scattering and the spin-orbit interaction caused by e-h exchange is also clarified, and the DT lineshape at the exciton and trion resonance are studied as a function of e-h exchange strength. In particular, as the exciton lineshape is determined by the interplay of linear exciton susceptibility and the bound-state two-exciton resonance in the T-matrix, the lineshape at the trion is similarly determined by the interplay of the linear trion susceptibility and the bound-state exciton-trion resonance in the T-matrix.


title: Maximal fluctuation exploitation in Gaussian information engines

abstract: Understanding the connections between information and thermodynamics has been among the most visible applications of stochastic thermodynamics. While recent theoretical advances have established that the second law of thermodynamics sets limits on information-to-energy conversion, it is currently unclear to what extent real systems can achieve the predicted theoretical limits. Using a simple model of an information engine that has recently been experimentally implemented, we explore the limits of information-to-energy conversion when an information engine's benefit is limited to output energy that can be stored. We find that restricting the engine's output in this way can limit its ability to convert information to energy. Nevertheless, a feedback control that inputs work can allow the engine to store energy at the highest achievable rate. These results sharpen our theoretical understanding of the limits of real systems that convert information to energy.


title: Rigidity for general semiconvex entire solutions to the sigma-2 equation

abstract: We show that every general semiconvex entire solution to the sigma-2 equation is a quadratic polynomial. A decade ago, this result was shown for almost convex solutions.


title: Thermal Image Super-Resolution Using Second-Order Channel Attention with Varying Receptive Fields

abstract: Thermal images model the long-infrared range of the electromagnetic spectrum and provide meaningful information even when there is no visible illumination. Yet, unlike imagery that represents radiation from the visible continuum, infrared images are inherently low-resolution due to hardware constraints. The restoration of thermal images is critical for applications that involve safety, search and rescue, and military operations. In this paper, we introduce a system to efficiently reconstruct thermal images. Specifically, we explore how to effectively attend to contrasting receptive fields (RFs) where increasing the RFs of a network can be computationally expensive. For this purpose, we introduce a deep attention to varying receptive fields network (AVRFN). We supply a gated convolutional layer with higher-order information extracted from disparate RFs, whereby an RF is parameterized by a dilation rate. In this way, the dilation rate can be tuned to use fewer parameters thus increasing the efficacy of AVRFN. Our experimental results show an improvement over the state of the art when compared against competing thermal image super-resolution methods.


title: Nonreciprocal model swimmer at intermediate Reynolds numbers

abstract: Metachronal swimming, the sequential beating of limbs with a small phase lag, is observed in many organisms at various scales, but has been studied mostly in the limits of high or low Reynolds numbers. Motivated by the swimming of brine shrimp, a mesoscale organism that operates at intermediate Reynolds numbers, we computationally studied a simple nonreciprocal 2D model that performs metachronal swimming. Our swimmer is composed of two pairs of paddles beating with a phase difference that are symmetrically attached to the sides of a flat body. We numerically solved the Navier-Stokes equations and used the immersed boundary method to model the interactions between the fluid and swimmer. To investigate the effect of inertia and geometry, we performed simulations varying the paddle spacing and the Reynolds numbers in the range $Re = 0.05 - 100$. In all cases, we observed back-and-forth motion during the cycle and a finite cycle-averaged swim speed at steady state. Interestingly, we found that the swim speed of the swimmer has nonmonotonic dependence on $Re$, with a maximum around $Re\approx1$, a flat minimum between $Re=20-30$ and an eventual increase for $Re>35$. To get more insight into the mechanism behind this relationship, we first decomposed the swim stroke of the swimmer and each paddle into power and recovery strokes and characterized the forward and backward motion. Unlike for reciprocal swimmers, here, for parts of the cycle there is competition between the two strokes of the paddles - the balance of which leads to net motion. We then studied the cycle-averaged, as well as, power-stroke-averaged and recovery-stroke-averaged fluid flows, and related differences in the fluid field to the nonmonotonic behavior of the swim speed. Our results suggest the existence of distinct motility mechanisms that develop as inertia increases within the intermediate-Re range.


title: Tradeoffs inherent in solar geoengineering peak-shaving strategies

abstract: Some have suggested solar geoengineering, cooling the planet by reflecting sunlight, could be deployed temporarily to shave the peak of global temperatures. This would use limited deployment to keep global mean temperatures below a given threshold, buying time for emissions reductions and negative emissions to reduce greenhouse gas concentrations. However, we suggest the trade-offs inherent in such scenarios are not widely appreciated. In this study, we explore general features of such scenarios, using a combination of simple analytical relationships between emissions and temperatures and modelling to illustrate the strong trade-off between the duration of geoengineering deployment and the temperature overshoot that is suppressed by geoengineering, as well as the scale of negative emissions that are required. Our results indicate that, unless climate sensitivity is at the high end of estimates, shaving a substantial peak, say, of 2K to 1.5K, would require at least a century of deployment of solar geoengineering and large-scale negative emissions, under a broad range of emissions trajectories. If negative emissions are not available at the scale of tens of gigatonnes of C02 per year, substantial peak-shaving will require deployment of geoengineering for multiple centuries. Our analysis emphasise the importance of factoring in the scale, feasibility, and cost of the negative emissions required to avoid indefinite commitment to solar geoengineering.


title: Method and portable bench for tests of the laser optical calibration system components for the Baikal-GVD underwater neutrino Cherenkov telescope

abstract: The large-scale deep underwater Cherenkov neutrino telescopes like Baikal-GVD, ANTARES or KM3NeT, require calibration and testing methods of their optical modules. These methods usually include laser-based systems which allow to check the telescope responses to the light and for real-time monitoring of the optical parameters of water such as absorption and scattering lengths, which show seasonal changes in natural reservoirs of water. We will present a testing method of a laser calibration system and a set of dedicated tools developed for Baikal- GVD, which includes a specially designed and built, compact, portable, and reconfigurable scanning station. This station is adapted to perform fast quality tests of the underwater laser sets just before their deployment in the telescope structure, even on ice, without darkroom. The testing procedure includes the energy stability test of the laser device, 3D scan of the light emission from the diffuser and attenuation test of the optical elements of the laser calibration system. The test bench consists primarily of an automatic mechanical scanner with a movable Si detector, beam splitter with a reference Si detector and, optionally, Q-switched diode-pumped solid-state laser used for laboratory scans of the diffusers. The presented test bench enables a three-dimensional scan of the light emission from diffusers, which are designed to obtain the isotropic distribution of photons around the point of emission. The results of the measurement can be easily shown on a 3D plot immediately after the test and may be also implemented to a dedicated program simulating photons propagation in water, which allows to check the quality of the diffuser in the scale of the Baikal-GVD telescope geometry.


title: Internet of things: a multiprotocol gateway as solution of the interoperability problem

abstract: One of the main challenges of the Internet of Things is the interoperability of highly heterogeneous devices, mainly in terms of the communication capabilities and network protocols used. As consequence, the interconnection model of the different devices involves an intermediary device, known as gateway. This gateway is a centralized element for the management of the devices that make up an IoT application. In addition, it is essential for the transmission of information to the Internet, especially when many IoT devices are not IP-based. This chapter describes a proposed model for an IoT gateway that allows the exchange of data through different wireless technologies and forwarding of such data to the Internet. The proposed gateway has important advantages such as: supporting for multiprotocol interconnectivity, the remote configuration of wireless nodes for sensor and actuators management; a flexible algorithm totranslate the data obtained by sensors into a uniform format for transmission to a cloud server; low energy consumption due to efficient data transfer over MQTT protocol. In order to demonstrate the usefulness of the developed gateway, a proof of concept test was implemented. The implemented scenario consists of 2 wireless nodes responsible for sensing environmental variables and transmitting data to the gateway node through different communication protocols. The obtained results show the feasibility for simultaneous data transmission from the remote wireless nodes to the gateway. Metrics on energy consumption in the devices are also presented.


title: A Deep Learning Approach to Predict Blood Pressure from PPG Signals

abstract: Blood Pressure (BP) is one of the four primary vital signs indicating the status of the body's vital (life-sustaining) functions. BP is difficult to continuously monitor using a sphygmomanometer (i.e. a blood pressure cuff), especially in everyday-setting. However, other health signals which can be easily and continuously acquired, such as photoplethysmography (PPG), show some similarities with the Aortic Pressure waveform. Based on these similarities, in recent years several methods were proposed to predict BP from the PPG signal. Building on these results, we propose an advanced personalized data-driven approach that uses a three-layer deep neural network to estimate BP based on PPG signals. Different from previous work, the proposed model analyzes the PPG signal in time-domain and automatically extracts the most critical features for this specific application, then uses a variation of recurrent neural networks called Long-Short-Term-Memory (LSTM) to map the extracted features to the BP value associated with that time window. Experimental results on two separate standard hospital datasets, yielded absolute errors mean and absolute error standard deviation for systolic and diastolic BP values outperforming prior works.


title: Quantum collision finding for homomorphic hash functions

abstract: Hash functions are a basic cryptographic primitive. Certain hash functions try to prove security against collision and preimage attacks by reductions to known hard problems. These hash functions usually have some additional properties that allow for that reduction. Hash functions which are additive or multiplicative are vulnerable to a quantum attack using the hidden subgroup problem algorithm for quantum computers. Using a quantum oracle to the hash, we can reconstruct the kernel of the hash function, which is enough to find collisions and second preimages. When the hash functions are additive with respect to the group operation in an Abelian group, there is always an efficient implementation of this attack. We present concrete attack examples to provable hash functions, including a preimage attack to $\oplus$-linear hash functions and for certain multiplicative homomorphic hash schemes.


title: The Pyro game: a slow intelligent fire

abstract: In the Firefighter problem, a fire breaks out at a vertex of a graph and at each subsequent time step, the firefighter chooses a vertex to protect and then the fire spreads from each burned vertex to every unprotected neighbour. The problem can be thought of as a simplified model for the spread of gossip or disease in a network. We introduce a new two-player variation called the Pyro game, in which at each step, the fire spreads from one burned vertex to all unprotected neighbours of that vertex. The fire is no longer automated and aims to maximize the number of burned vertices. We show, that unlike the Firefighter problem, one firefighter can contain a fire on the Cartesian grid in the Pyro game. We also study both the Pyro Game and the Firefighter Problem on the infinite strong grid and the complexity of the Pyro game.


title: Near-Optimal Spanners for General Graphs in (Nearly) Linear Time

abstract: Let $G = (V,E,w)$ be a weighted undirected graph on $|V| = n$ vertices and $|E| = m$ edges, let $k \ge 1$ be any integer, and let $\epsilon < 1$ be any parameter. We present the following results on fast constructions of spanners with near-optimal sparsity and lightness, which culminate a long line of work in this area. (By near-optimal we mean optimal under Erd\H{o}s' girth conjecture and disregarding the $\epsilon$-dependencies.) - There are (deterministic) algorithms for constructing $(2k-1)(1+\epsilon)$-spanners for $G$ with a near-optimal sparsity of $O(n^{1/k} \log(1/\epsilon)/\epsilon))$. The first algorithm can be implemented in the pointer-machine model within time $O(m\alpha(m,n) \log(1/\epsilon)/\epsilon) + SORT(m))$, where $\alpha( , )$ is the two-parameter inverse-Ackermann function and $SORT(m)$ is the time needed to sort $m$ integers. The second algorithm can be implemented in the WORD RAM model within time $O(m \log(1/\epsilon)/\epsilon))$. - There is a (deterministic) algorithm for constructing a $(2k-1)(1+\epsilon)$-spanner for $G$ that achieves a near-optimal bound of $O(n^{1/k}\mathrm{poly}(1/\epsilon))$ on both sparsity and lightness. This algorithm can be implemented in the pointer-machine model within time $O(m\alpha(m,n) \mathrm{poly}(1/\epsilon) + SORT(m))$ and in the WORD RAM model within time $O(m \alpha(m,n) \mathrm{poly}(1/\epsilon))$. The previous fastest constructions of $(2k-1)(1+\epsilon)$-spanners with near-optimal sparsity incur a runtime of is $O(\min\{m(n^{1+1/k}) + n\log n,k n^{2+1/k}\})$, even regardless of the lightness. Importantly, the greedy spanner for stretch $2k-1$ has sparsity $O(n^{1/k})$ -- with no $\epsilon$-dependence whatsoever, but its runtime is $O(m(n^{1+1/k} + n\log n))$. Moreover, the state-of-the-art lightness bound of any $(2k-1)$-spanner is poor, even regardless of the sparsity and runtime.


title: Anomaly Detection with Neural Parsers That Never Reject

abstract: Reinforcement learning has recently shown promise as a technique for training an artificial neural network to parse sentences in some unknown format, through a body of work known as RL-GRIT. A key aspect of the RL-GRIT approach is that rather than explicitly inferring a grammar that describes the format, the neural network learns to perform various parsing actions (such as merging two tokens) over a corpus of sentences, with the goal of maximizing the estimated frequency of the resulting parse structures. This can allow the learning process to more easily explore different action choices, since a given choice may change the optimality of the parse (as expressed by the total reward), but will not result in the failure to parse a sentence. However, this also presents a limitation: because the trained neural network can successfully parse any sentence, it cannot be directly used to identify sentences that deviate from the format of the training sentences, i.e., that are anomalous. In this paper, we address this limitation by presenting procedures for extracting production rules from the neural network, and for using these rules to determine whether a given sentence is nominal or anomalous. When a sentence is anomalous, an attempt is made to identify the location of the anomaly. We empirically demonstrate that our approach is capable of grammatical inference and anomaly detection for both non-regular formats and those containing regions of high randomness/entropy. While a format with high randomness typically requires large sets of production rules, we propose a two pass grammatical inference method to generate parsimonious rule sets for such formats. By further improving parser learning, and leveraging the presented rule extraction and anomaly detection algorithms, one might begin to understand common errors, either benign or malicious, in practical formats.


title: Structural Guidance for Transformer Language Models

abstract: Transformer-based language models pre-trained on large amounts of text data have proven remarkably successful in learning generic transferable linguistic representations. Here we study whether structural guidance leads to more human-like systematic linguistic generalization in Transformer language models without resorting to pre-training on very large amounts of data. We explore two general ideas. The "Generative Parsing" idea jointly models the incremental parse and word sequence as part of the same sequence modeling task. The "Structural Scaffold" idea guides the language model's representation via additional structure loss that separately predicts the incremental constituency parse. We train the proposed models along with a vanilla Transformer language model baseline on a 14 million-token and a 46 million-token subset of the BLLIP dataset, and evaluate models' syntactic generalization performances on SG Test Suites and sized BLiMP. Experiment results across two benchmarks suggest converging evidence that generative structural supervisions can induce more robust and humanlike linguistic generalization in Transformer language models without the need for data intensive pre-training.


title: Deep Feature Tracker: A Novel Application for Deep Convolutional Neural Networks

abstract: Feature tracking is the building block of many applications such as visual odometry, augmented reality, and target tracking. Unfortunately, the state-of-the-art vision-based tracking algorithms fail in surgical images due to the challenges imposed by the nature of such environments. In this paper, we proposed a novel and unified deep learning-based approach that can learn how to track features reliably as well as learn how to detect such reliable features for tracking purposes. The proposed network dubbed as Deep-PT, consists of a tracker network which is a convolutional neural network simulating cross-correlation in terms of deep learning and two fully connected networks that operate on the output of intermediate layers of the tracker to detect features and predict trackability of the detected points. The ability to detect features based on the capabilities of the tracker distinguishes the proposed method from previous algorithms used in this area and improves the robustness of the algorithms against dynamics of the scene. The network is trained using multiple datasets due to the lack of specialized dataset for feature tracking datasets and extensive comparisons are conducted to compare the accuracy of Deep-PT against recent pixel tracking algorithms. As the experiments suggest, the proposed deep architecture deliberately learns what to track and how to track and outperforms the state-of-the-art methods.


title: Soft Calibration Objectives for Neural Networks

abstract: Optimal decision making requires that classifiers produce uncertainty estimates consistent with their empirical accuracy. However, deep neural networks are often under- or over-confident in their predictions. Consequently, methods have been developed to improve the calibration of their predictive uncertainty both during training and post-hoc. In this work, we propose differentiable losses to improve calibration based on a soft (continuous) version of the binning operation underlying popular calibration-error estimators. When incorporated into training, these soft calibration losses achieve state-of-the-art single-model ECE across multiple datasets with less than 1% decrease in accuracy. For instance, we observe an 82% reduction in ECE (70% relative to the post-hoc rescaled ECE) in exchange for a 0.7% relative decrease in accuracy relative to the cross entropy baseline on CIFAR-100. When incorporated post-training, the soft-binning-based calibration error objective improves upon temperature scaling, a popular recalibration method. Overall, experiments across losses and datasets demonstrate that using calibration-sensitive procedures yield better uncertainty estimates under dataset shift than the standard practice of using a cross entropy loss and post-hoc recalibration methods.


title: Comparing object recognition in humans and deep convolutional neural networks -- An eye tracking study

abstract: Deep convolutional neural networks (DCNNs) and the ventral visual pathway share vast architectural and functional similarities in visual challenges such as object recognition. Recent insights have demonstrated that both hierarchical cascades can be compared in terms of both exerted behavior and underlying activation. However, these approaches ignore key differences in spatial priorities of information processing. In this proof-of-concept study, we demonstrate a comparison of human observers (N = 45) and three feedforward DCNNs through eye tracking and saliency maps. The results reveal fundamentally different resolutions in both visualization methods that need to be considered for an insightful comparison. Moreover, we provide evidence that a DCNN with biologically plausible receptive field sizes called vNet reveals higher agreement with human viewing behavior as contrasted with a standard ResNet architecture. We find that image-specific factors such as category, animacy, arousal, and valence have a direct link to the agreement of spatial object recognition priorities in humans and DCNNs, while other measures such as difficulty and general image properties do not. With this approach, we try to open up new perspectives at the intersection of biological and computer vision research.


title: Absolute frequency measurement of the 171Yb optical lattice clock at KRISS using TAI for over a year

abstract: We report a measurement of the absolute frequency of the 1S0-3P0 transition in the 171Yb optical lattice clock at KRISS (KRISS-Yb1) for 14 months, which was referenced to the SI second by primary and secondary standards worldwide via TAI (International Atomic Time). The determined absolute frequency is 518 295 836 590 863.75(14) Hz with the relative frequency uncertainty of 2.6x10^-16, which agrees well with other reports. This result is expected to contribute to the future update of the CIPM recommendation frequency of the secondary frequency standards.


title: A Machine-learning Based Initialization for Joint Statistical Iterative Dual-energy CT with Application to Proton Therapy

abstract: Dual-energy CT (DECT) has been widely investigated to generate more informative and more accurate images in the past decades. For example, Dual-Energy Alternating Minimization (DEAM) algorithm achieves sub-percentage uncertainty in estimating proton stopping-power mappings from experimental 3-mm collimated phantom data. However, elapsed time of iterative DECT algorithms is not clinically acceptable, due to their low convergence rate and the tremendous geometry of modern helical CT scanners. A CNN-based initialization method is introduced to reduce the computational time of iterative DECT algorithms. DEAM is used as an example of iterative DECT algorithms in this work. The simulation results show that our method generates denoised images with greatly improved estimation accuracy for adipose, tonsils, and muscle tissue. Also, it reduces elapsed time by approximately 5-fold for DEAM to reach the same objective function value for both simulated and real data.


title: Design of microring resonators integrated with 2D graphene oxide films for four-wave mixing

abstract: We theoretically investigate and optimize the performance of four-wave mixing (FWM) in microring resonators (MRRs) integrated with two-dimensional (2D) layered graphene oxide (GO) films. Owing to the interaction between the MRRs and the highly nonlinear GO films as well as to the resonant enhancement effect, the FWM efficiency in GO-coated MRRs can be significantly improved. Based on previous experiments, we perform detailed analysis for the influence of the GO film parameters and MRR coupling strength on the FWM conversion efficiency (CE) of the hybrid MRRs. By optimizing the device parameters to balance the trade-off between the Kerr nonlinearity and loss, we achieve a high CE enhancement of ~18.6 dB relative to the uncoated MRR, which is ~8.3 dB higher than previous experimental results. The influence of photo-thermal changes in the GO films as well as variations in the MRR parameters such as the ring radius and waveguide dispersion on the FWM performance is also discussed. These results highlight the significantly improved FWM performance that can be achieved in MRRs incorporating GO films and provide a guide for optimizing their FWM performance.


title: Student ownership and understanding of multi-week final projects

abstract: National calls to transform laboratory courses by making them more discovery-based can be met by engaging students in multi-week final projects. One plausible outcome of this approach is that students may feel ownership of their projects. We define ownership as a dynamic relationship between students and their projects characterized by three student-project interactions that evolve over three project phases. Student-project interactions include students' contributions to, emotional responses to, and new understanding of the project. Phases include choosing the topic and team, carrying out the research, and creating and presenting end-of-project deliverables. Drawing on interviews with students collected as part of a multi-year, multi-institutional study, this paper will elaborate on the evolution of students' own new knowledge about the project across different project phases. Throughout our work, we compare the type of ownership that manifests when students work on projects in a lab course to the type of ownership gardeners feel when tending to plots in a community garden. We end with recommendations for instructors wishing to foster project ownership in their lab courses.


title: Use of Machine Learning for gamma/hadron separation with HAWC

abstract: Background showers triggered by hadrons represent over 99.9% of all particles arriving at ground-based gamma-ray observatories. An important stage in the data analysis of these observatories, therefore, is the removal of hadron-triggered showers. Currently, the High-Altitude Water Cherenkov (HAWC) gamma-ray observatory employs an algorithm based on a single cut in two variables, unlike other ground-based gamma-ray observatories (e.g. H.E.S.S., VERITAS), which employ a large number of variables to separate the primary particles. In this work, we explore machine learning techniques (Boosted Decision Trees and Neural Networks) to identify the primary particles detected by HAWC. Our new gamma/hadron separation techniques were tested on data from the Crab nebula, the standard reference in Very High Energy astronomy, showing an improvement compared to the standard HAWC background rejection method.


title: Monitoring the radio galaxy M87 with HAWC

abstract: Studies of radio galaxies at TeV energies are fascinating because their jets are misaligned concerning our sightline. Thus, it provides us with a unique opportunity to study the structure of their jets, the radiative processes, and the acceleration mechanisms involved in them. In addition, some radio galaxies have presented variability in their emission, like the giant radio galaxy M87, which has reported several activity periods. Due to its duty cycle >95% and instantaneous field of view of 2 sr, HAWC provides daily monitoring of variable sources visible from the Northern Hemisphere. In this work, we show the results of monitoring M87 between January 2015 and December 2018. HAWC's observations are consistent with the low activity state reported by other instruments (like H.E.S.S and MAGIC). However, after September 2017 (~MJD 58000), the HAWC measurements of M87 show hints of higher activity.


title: On The State of Data In Computer Vision: Human Annotations Remain Indispensable for Developing Deep Learning Models

abstract: High-quality labeled datasets play a crucial role in fueling the development of machine learning (ML), and in particular the development of deep learning (DL). However, since the emergence of the ImageNet dataset and the AlexNet model in 2012, the size of new open-source labeled vision datasets has remained roughly constant. Consequently, only a minority of publications in the computer vision community tackle supervised learning on datasets that are orders of magnitude larger than Imagenet. In this paper, we survey computer vision research domains that study the effects of such large datasets on model performance across different vision tasks. We summarize the community's current understanding of those effects, and highlight some open questions related to training with massive datasets. In particular, we tackle: (a) The largest datasets currently used in computer vision research and the interesting takeaways from training on such datasets; (b) The effectiveness of pre-training on large datasets; (c) Recent advancements and hurdles facing synthetic datasets; (d) An overview of double descent and sample non-monotonicity phenomena; and finally, (e) A brief discussion of lifelong/continual learning and how it fares compared to learning from huge labeled datasets in an offline setting. Overall, our findings are that research on optimization for deep learning focuses on perfecting the training routine and thus making DL models less data hungry, while research on synthetic datasets aims to offset the cost of data labeling. However, for the time being, acquiring non-synthetic labeled data remains indispensable to boost performance.


title: Searching for the hydrogen 21 cm line in cosmos

abstract: It is widely agreed that studying the 21 cm emission line from neutral hydrogen may be our best hope for understanding the creation of the first structures during the dark ages and the Epoch of Reionization (EoR). This hyperfine transition occurs as a result of the spin-spin interaction between the electron and the proton in hydrogen: The energy of the parallel spins (triplet) state is greater than that of the anti-parallel spins (singlet). This transition is strongly banned, with a transition probability of $2.6 \times 10^{-15} \mathrm{s}^{-1}$ being exceedingly low. Despite its low chance, the 21 cm hyperfine transition is one of the most important tools in observational astronomy due to the abundance of hydrogen in the Universe. The relative occupation number of the ground state and the occupation number of the excited state determine the strength of the 21 cm emission or absorption. Collisions or the so-called Wouthuysen-Field effect can both excite the 21 cm line. In this article, we give a brief review of the basic physics of 21 cm emission and highlight the current main worldwide 21 cm signal experiments. We also discuss several related concepts in astronomy and cosmology with these observations.


title: Role of ambient humidity underestimated in research on correlation between radioactive decay rates and space weather

abstract: In recent work, Milian-Sanchez et al. observed fluctuations in radioactive decay rate measurement series, and after excluding environmental influences (measured indoors) as root causes, they looked for possible correlations with astrophysical variables. In spite of the authors efforts to investigate possible influences of environmental parameters (such as ambient temperature, pressure and humidity) on the detectors stability, it turns out that the influence of ambient humidity on the instrumentation has been underestimated.


title: Margin-Aware Intra-Class Novelty Identification for Medical Images

abstract: Traditional anomaly detection methods focus on detecting inter-class variations while medical image novelty identification is inherently an intra-class detection problem. For example, a machine learning model trained with normal chest X-ray and common lung abnormalities, is expected to discover and flag idiopathic pulmonary fibrosis which a rare lung disease and unseen by the model during training. The nuances from intra-class variations and lack of relevant training data in medical image analysis pose great challenges for existing anomaly detection methods. To tackle the challenges, we propose a hybrid model - Transformation-based Embedding learning for Novelty Detection (TEND) which without any out-of-distribution training data, performs novelty identification by combining both autoencoder-based and classifier-based method. With a pre-trained autoencoder as image feature extractor, TEND learns to discriminate the feature embeddings of in-distribution data from the transformed counterparts as fake out-of-distribution inputs. To enhance the separation, a distance objective is optimized to enforce a margin between the two classes. Extensive experimental results on both natural image datasets and medical image datasets are presented and our method out-performs state-of-the-art approaches.


title: Finite-size scaling around the critical point in the heavy quark region of QCD

abstract: Finite-size scaling is investigated in detail around the critical point in the heavy-quark region of nonzero temperature QCD. Numerical simulations are performed with large spatial volumes up to the aspect ratio $N_s/N_t=12$ at a fixed lattice spacing with $N_t=4$. We show that the Binder cumulant and the distribution function of the Polyakov loop follow the finite-size scaling in the $Z(2)$ universality class for large spatial volumes with $N_s/N_t \ge 9$, while, for $N_s/N_t \le 8$, the Binder cumulant becomes inconsistent with the $Z(2)$ scaling. To realize the large-volume simulations in the heavy-quark region, we adopt the hopping parameter expansion for the quark determinant: We generate gauge configurations using the leading order action including the Polyakov loop term for $N_t=4$, and incorporate the next-to-leading order effects in the measurements by the multipoint reweighting method. We find that the use of the leading-order configurations is crucially effective in suppressing the overlapping problem in the reweighting and thus reducing the statistical errors.


title: Measurement of Transverse Spin Dependent Azimuthal Correlations of Charged Pion(s) in $p^{\uparrow} p$ Collisions at $\sqrt s = 200$ GeV at STAR

abstract: At the leading twist, the transversity distribution function, $h^{q}_{1}(x)$, where $x$ is the longitudinal momentum fraction of the proton carried by quark $q$, encodes the transverse spin structure of the nucleon. Extraction of it is difficult because of its chiral-odd nature. In transversely polarized proton-proton collisions ($p^\uparrow p$), $h_{1}^{q}(x)$ can be coupled with another chiral-odd partner, a spin-dependent fragmentation function (FF). The resulting asymmetries in hadron(s) azimuthal correlations directly probe $h_{1}^{q}(x)$. We report the measurement of correlation asymmetries for charged pion(s) in $p^\uparrow p$, through the Collins and the Interference FF channel.


title: Critical threshold for global regularity of Euler-Monge-Amp\`ere system with radial symmetry

abstract: We study the global wellposedness of the Euler-Monge-Amp\`ere (EMA) system. We obtain a sharp, explicit critical threshold in the space of initial configurations which guarantees the global regularity of EMA system with radially symmetric initial data. The result is obtained using two independent approaches -- one using spectral dynamics of Liu & Tadmor [Comm. Math. Physics 228(3):435-466, 2002] and another based on the geometric approach of Brenier & Loeper [Geom. Funct. Analysis 14(6):1182--1218, 2004]. The results are extended to 2D radial EMA with swirl.


title: Entropy production related properties of first passage process

abstract: With nontrivial entropy production, first passage process is one of the most common nonequilibrium process in stochastic thermodynamics. Using one dimensional birth and death precess as a model framework, approximated expressions of mean first passage time (FPT), mean total number of jumps (TNJ), and their coefficients of variation (CV), are obtained for the case far from equilibrium. Consequently, uncertainty relations for FPT and TNJ are presented. Generally, mean FPT decreases exponentially with entropy production, while mean TNJ decreases exponentially first and then tends to a starting site dependent limit. For forward biased process, the CV of TNJ decreases exponentially with entropy production, while that of FPT decreases exponentially first and then tends to a starting site dependent limit. For backward biased process, both CVs of FPT and TNJ tend to one for large absolute values of entropy production. Related properties about the case of equilibrium are also addressed briefly for comparison.


title: Nonlinear Schrodinger Equation Solitons on Quantum Droplets

abstract: Irrotational ow of a spherical thin liquid layer surrounding a rigid core is described using the defocusing nonlinear Schrodinger equation. Accordingly, azimuthal moving nonlinear waves are modeled by periodic dark solitons expressed by elliptic functions. In the quantum regime the algebraic Bethe ansatz is used in order to capture the energy levels of such motions, which we expect to be relevant for the dynamics of the nuclear clusters in deformed heavy nuclei surface modeled by quantum liquid drops. In order to validate the model we match our theoretical energy spectra with experimental results on energy, angular momentum and parity for alpha particle clustering nuclei.


title: Modelling the Earth's magnetic field with magnetometer data from a Raspberry Pi on board the ISS

abstract: We used data from a magnetometer on board the ISS to simulate the Earth's magnetic field. For that purpose, we generated two models from the data and compared them: an off-centered dipole model and a centered multipolar model. We found out that measuring the magnetic field along only two orbits is not enough to provide a good estimate of the Earth's magnetic field, as they chart a too small portion of the Earth's surface. This was our original work within the "Astro Pi Challenge", a project ran by ESA and the Raspberry Pi Foundation, and the magnetometer used was the one incorporated in the Raspberry Pi Sense Hat. Here, we improved by performing the same analysis using data from ten orbits and learned that the multipolar model already provides a good approximation to the Earth's magnetic field with this number of orbits. We also noticed a deviation of the data collected from what would be expected from the IGRF model, in both sets of data.


abstract: This work considers an extension of the Kuramoto model with run-and-tumble dynamics -- a type of self-propelled motion. The difference between the extended and the original model is that in the extended version angular velocity of individual particles is no longer fixed but can change sporadically with a new velocity drawn from a distribution g(w). Because the Kuramoto model undergoes a phase transition, it offers a simple case study for investigating phase transition for a system with self-propelled particles.


title: Proximal Quasi-Newton Methods for Multiobjective Optimization Problems

abstract: We introduce some new proximal quasi-Newton methods for unconstrained multiobjective optimization problems (in short, UMOP), where each objective function is the sum of a twice continuously differentiable strongly convex function and a proper lower semicontinuous convex but not necessarily differentiable function. We propose proximal BFGS method, proximal self-scaling BFGS method, and proximal Huang BFGS method for (UMOP) with both line searches and without line searches cases. Under mild assumputions, we show that each accumulation point of the sequence generated by these algorithms, if exists, is a Pareto stationary point of the (UMOP). Moreover, we present their applications in both constrained multiobjective optimization problems and robust multiobjective optimization problems. In particular, for robust multiobjective optimization problems, we show that the subproblems of proximal quasi-Newton algorithms can be regarded as quadratic minimization problems with quadratic inequality constraints. Numerical experiments are also carried out to verify the effectiveness of the proposed proximal quasi-Newton methods.


title: A Moving Embedded Boundary Approach For The Compressible Navier-Stokes Equations In A Block-Structured Adaptive Refinement Framework

abstract: A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement framework of AMReX. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds $-$ acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries $-$ an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating the conservative nature of the moving boundary algorithm. Viscous test cases involve that of a horizontally moving cylinder at $Re=40$, an inline oscillating cylinder at $Re=100$, and a transversely oscillating cylinder at $Re=185$. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.


title: Structure Amplification on Multi-layer Stochastic Block Models

abstract: Much of the complexity of social, biological, and engineered systems arises from a network of complex interactions connecting many basic components. Network analysis tools have been successful at uncovering latent structure termed communities in such networks. However, some of the most interesting structure can be difficult to uncover because it is obscured by the more dominant structure. Our previous work proposes a general structure amplification technique called HICODE that uncovers many layers of functional hidden structure in complex networks. HICODE incrementally weakens dominant structure through randomization allowing the hidden functionality to emerge, and uncovers these hidden structure in real-world networks that previous methods rarely uncover. In this work, we conduct a comprehensive and systematic theoretical analysis on the hidden community structure. In what follows, we define multi-layer stochastic block model, and provide theoretical support using the model on why the existence of hidden structure will make the detection of dominant structure harder compared with equivalent random noise. We then provide theoretical proofs that the iterative reducing methods could help promote the uncovering of hidden structure as well as boosting the detection quality of dominant structure.


title: Physics-informed Dyna-Style Model-Based Deep Reinforcement Learning for Dynamic Control

abstract: Model-based reinforcement learning (MBRL) is believed to have much higher sample efficiency compared to model-free algorithms by learning a predictive model of the environment. However, the performance of MBRL highly relies on the quality of the learned model, which is usually built in a black-box manner and may have poor predictive accuracy outside of the data distribution. The deficiencies of the learned model may prevent the policy from being fully optimized. Although some uncertainty analysis-based remedies have been proposed to alleviate this issue, model bias still poses a great challenge for MBRL. In this work, we propose to leverage the prior knowledge of underlying physics of the environment, where the governing laws are (partially) known. In particular, we developed a physics-informed MBRL framework, where governing equations and physical constraints are utilized to inform the model learning and policy search. By incorporating the prior information of the environment, the quality of the learned model can be notably improved, while the required interactions with the environment are significantly reduced, leading to better sample efficiency and learning performance. The effectiveness and merit have been demonstrated over a handful of classic control problems, where the environments are governed by canonical ordinary/partial differential equations.


title: Point-wise posteriori phase estimation in high-precision fringe projection profilometry

abstract: In fringe projection profilometry, the high-order harmonics information of non-sinusoidal fringes will lead to errors in the phase estimation. In order to solve this problem, a point-wise posterior phase estimation (PWPPE) method based on deep learning technique is proposed in this paper. The complex nonlinear mapping relationship between the multiple gray values and the sine / cosine value of the phase is constructed by using the feedforward neural network model. After the model training, it can estimate the phase values of each pixel location, and the accuracy is higher than the point-wise least-square (PWLS) method. To further verify the effectiveness of this method, a face mask is measured, the traditional PWLS method and the proposed PWPPE method are employed, respectively. The comparison results show that the traditional method is with periodic phase errors, while the proposed PWPPE method can effectively eliminate such phase errors caused by non-sinusoidal fringes.


title: Solution of Mumford's second problem

abstract: A complete solution of Mumford's second problem about representation of theta derivatives with rational characteristics in terms of theta constants with rational characteristics is found. An explicit formula for computing such an expression for theta derivative with an arbitrary rational characteristic is derived, and illustrated with examples. Expressions for theta derivatives appear to be homogeneous of degree $3$ with respect to theta constants.


title: Simple, Fast, and Flexible Framework for Matrix Completion with Infinite Width Neural Networks

abstract: Matrix completion problems arise in many applications including recommendation systems, computer vision, and genomics. Increasingly larger neural networks have been successful in many of these applications, but at considerable computational costs. Remarkably, taking the width of a neural network to infinity allows for improved computational performance. In this work, we develop an infinite width neural network framework for matrix completion that is simple, fast, and flexible. Simplicity and speed come from the connection between the infinite width limit of neural networks and kernels known as neural tangent kernels (NTK). In particular, we derive the NTK for fully connected and convolutional neural networks for matrix completion. The flexibility stems from a feature prior, which allows encoding relationships between coordinates of the target matrix, akin to semi-supervised learning. The effectiveness of our framework is demonstrated through competitive results for virtual drug screening and image inpainting/reconstruction. We also provide an implementation in Python to make our framework accessible on standard hardware to a broad audience.


title: A Unified Convergence Analysis of First Order Convex Optimization Methods via Strong Lyapunov Functions

abstract: We present a unified convergence analysis for first order convex optimization methods using the concept of strong Lyapunov conditions. Combining this with suitable time scaling factors, we are able to handle both convex and strong convex cases, and establish contraction properties of Lyapunov functions for many existing ordinary differential equation models. Then we derive prevailing first order optimization algorithms, such as proximal gradient methods, heavy ball methods (also known as momentum methods), Nesterov accelerated gradient methods, and accelerated proximal gradient methods from numerical discretizations of corresponding dynamical systems. We also apply strong Lyapunov conditions to the discrete level and provide a more systematical analysis framework. Another contribution is a novel second order dynamical system called Hessian-driven Nesterov accelerated gradient flow which can be used to design and analyze accelerated first order methods for smooth and non-smooth convex optimizations.


title: Online Spatio-temporal Calibration of Tightly-coupled Ultrawideband-aided Inertial Localization

abstract: The combination of ultrawideband (UWB) radios and inertial measurement units (IMU) can provide accurate positioning in environments where the Global Positioning System (GPS) service is either unavailable or has unsatisfactory performance. The two sensors, IMU and UWB radio, are often not co-located on a moving system. The UWB radio is typically located at the extremities of the system to ensure reliable communication, whereas the IMUs are located closer to its center of gravity. Furthermore, without hardware or software synchronization, data from heterogeneous sensors can arrive at different time instants resulting in temporal offsets. If uncalibrated, these spatial and temporal offsets can degrade the positioning performance. In this paper, using observability and identifiability criteria, we derive the conditions required for successfully calibrating the spatial and the temporal offset parameters of a tightly-coupled UWB-IMU system. We also present an online method for jointly calibrating these offsets. The results show that our calibration approach results in improved positioning accuracy while simultaneously estimating (i) the spatial offset parameters to millimeter precision and (ii) the temporal offset parameter to millisecond precision.


title: Factorially many maximum matchings close to the Erd\H{o}s-Gallai bound

abstract: A classical result of Erd\H{o}s and Gallai determines the maximum size $m(n,\nu)$ of a graph $G$ of order $n$ and matching number $\nu n$. We show that $G$ has factorially many maximum matchings provided that its size is sufficiently close to $m(n,\nu)$.


title: Revised event rates for extreme and extremely large mass-ratio inspirals

abstract: One of the main targets of the Laser Interferometer Space Antenna (LISA) is the detection of extreme mass-ratio inspirals (EMRIs) and extremely large mass-ratio inspirals (X-MRIs). Their orbits are expected to be highly eccentric and relativistic when entering the LISA band. Under these circumstances, the inspiral time-scale given by Peters' formula loses precision and the shift of the last-stable orbit (LSO) caused by the massive black hole spin could influence the event rates estimate. We re-derive EMRIs and X-MRIs event rates by implementing two different versions of a Kerr loss-cone angle that includes the shift in the LSO, and a corrected version of Peters' time-scale that accounts for eccentricity evolution, 1.5 post-Newtonian hereditary fluxes, and spin-orbit coupling. The main findings of our study are summarized as follows: (1) implementing a Kerr loss-cone changes the event rates by a factor ranging between 0.9 and 1.1; (2) the high-eccentricity limit of Peters' formula offers a reliable inspiral time-scale for EMRIs and X-MRIs, resulting in an event rate estimate that deviates by a factor of about 0.9 to 3 when compared to event rates computed with the corrected version of Peters' time-scale and the usual loss-cone definition. (3) Event rates estimates for systems with a wide range of eccentricities should be revised. Peters' formula overestimates the inspiral rates of highly eccentric systems by a factor of about 8 to 30 compared to the corrected values. Besides, for e$_0 \lesssim$0.8, implementing the corrected version of Peters' formula would be necessary to obtain accurate estimates.


title: Redox and Molecular Diffusion in 2D van der Waals Space

abstract: Understanding charge transfer (CT) between two chemical entities and subsequent change in their charge densities is essential not only for molecular species but also for various low-dimensional materials. Because of their extremely high fraction of surface atoms, two-dimensional (2-D) materials are most susceptible to charge exchange and exhibit drastically different physicochemical properties depending on their charge density. In this regard, spontaneous and uncontrollable ionization of graphene in the ambient air has caused much confusion and technical difficulty in achieving experimental reproducibility since its first report in 2004. Moreover, the same ambient hole doping was soon observed in 2-D semiconductors, which implied that a common mechanism should be operative and apply to other low-dimensional materials universally. In this Account, we review our breakthroughs in unraveling the chemical origin and mechanistic requirements of the hidden CT reactions using 2-D crystals. We developed in-situ optical methods to quantify charge density using Raman and photoluminescence (PL) spectroscopy and imaging. Using gas and temperature-controlled in-situ measurements, we revealed that the electrical holes are injected by the oxygen reduction reaction (ORR): $O_{2}$ + $4H^{+}$ + $4e^{-}$ $\rightleftharpoons$ $2H_{2}O$, which was independently verified by pH dependence in HCl solutions. In addition to oxygen and water vapor, the overall CT reaction requires hydrophilic dielectric substrates, which assist hydration of the sample-substrate interface. The interface-localized reaction allowed us to visualized and control interfacial molecular diffusion and CT by varing the 2-D gap spacing and introducing defects. The complete mechanism of the fundamental charge exchange summarized in this Account will be essential in exploring material and device properties of other low dimensional materials.


title: Two-photon sideband transition in a driven quantum Rabi model : Quantitative discussions with derived longitudinal drives and beyond the rotating wave approximation

abstract: In this work, we analytically and numerically study the sideband transition dynamics of the driven quantum Rabi model (QRM). We focus in particular on the conditions when the external transverse drive fields induce first-order sideband transitions. Inducing sideband transitions between two different systems is an essential technique for various physical models, including the QRM. However, despite its importance, a precise analytical study has not been reported yet that successfully explains the sideband transition rates in a driven QRM applicable for all system parameter configurations. In our study, we analytically derive the sideband transition rates based on second-order perturbation theory, not relying on the rotating wave approximation (RWA) \cite{RWA}. Our formula are valid for all ranges of drive frequencies and system's parameters. Our analytical derived formula agrees well with the numerical results in a regime of moderate drive amplitudes. Interestingly, we have found a non-trivial longitudinal drive effect derived from the transverse drive Hamiltonian. This accounts for significant corrections to the sideband transition rates that are expected without considering the derived longitudinal effect. Using this approach, one can precisely estimate the sideband transition rates in the QRM not confining themselves within specific parameter regimes. This provides important contributions for understanding experiments described by the driven QRM.


title: Learning to Control DC Motor for Micromobility in Real Time with Reinforcement Learning

abstract: Autonomous micromobility has been attracting the attention of researchers and practitioners in recent years. A key component of many micro-transport vehicles is the DC motor, a complex dynamical system that is continuous and non-linear. Learning to quickly control the DC motor in the presence of disturbances and uncertainties is desired for various applications that require robustness and stability. Techniques to accomplish this task usually rely on a mathematical system model, which is often insufficient to anticipate the effects of time-varying and interrelated sources of non-linearities. While some model-free approaches have been successful at the task, they rely on massive interactions with the system and are trained in specialized hardware in order to fit a highly parameterized controller. In this work, we learn to steer a DC motor via sample-efficient reinforcement learning. Using data collected from hardware interactions in the real world, we additionally build a simulator to experiment with a wide range of parameters and learning strategies. With the best parameters found, we learn an effective control policy in one minute and 53 seconds on a simulation and in 10 minutes and 35 seconds on a physical system.


title: Pose-Guided Feature Learning with Knowledge Distillation for Occluded Person Re-Identification

abstract: Occluded person re-identification (ReID) aims to match person images with occlusion. It is fundamentally challenging because of the serious occlusion which aggravates the misalignment problem between images. At the cost of incorporating a pose estimator, many works introduce pose information to alleviate the misalignment in both training and testing. To achieve high accuracy while preserving low inference complexity, we propose a network named Pose-Guided Feature Learning with Knowledge Distillation (PGFL-KD), where the pose information is exploited to regularize the learning of semantics aligned features but is discarded in testing. PGFL-KD consists of a main branch (MB), and two pose-guided branches, \ieno, a foreground-enhanced branch (FEB), and a body part semantics aligned branch (SAB). The FEB intends to emphasise the features of visible body parts while excluding the interference of obstructions and background (\ieno, foreground feature alignment). The SAB encourages different channel groups to focus on different body parts to have body part semantics aligned representation. To get rid of the dependency on pose information when testing, we regularize the MB to learn the merits of the FEB and SAB through knowledge distillation and interaction-based training. Extensive experiments on occluded, partial, and holistic ReID tasks show the effectiveness of our proposed network.


title: Model structures and recollements induced by duality pairs

abstract: We give some equivalent characterizations of $\mathcal{GP}$, the class of Gorenstein $(\mathcal{L}, \mathcal{A})$-projective modules, and construct some model structures associated to duality pairs and Frobenius pairs. Moreover, some rings are described by Frobenius pairs. Meanwhile, we investigate strongly Gorenstein $(\mathcal{L}, \mathcal{A})$-projective modules and obtain some equivalent characterizations of them. Also, some model structures and recollements associated to strongly Gorenstein $(\mathcal{L}, \mathcal{A})$-projective modules are constructed.


title: Mapping the structural, magnetic and electronic behavior of (Eu$_{1-x}$Ca$_{x}$)$_{2}$Ir$_{2}$O$_{7}$ across a metal-insulator transition

abstract: In this study, we employ bulk electronic properties characterization and x-ray scattering/spectroscopy techniques to map the structural, magnetic and electronic properties of (Eu$_{1-x}$Ca$_{x}$)$_{2}$Ir$_{2}$O$_{7}$ as a function of Ca-doping. As expected, the metal-insulator transition temperature, $T_{MIT}$, decreases with Ca-doping until a metallic state is realized down to 2 K. In contrast, $T_{AFM}$ becomes decoupled from the MIT and (likely short-range) AFM order persists into the metallic regime. This decoupling is understood as a result of the onset of an electronically phase separated state, the occurrence of which seemingly depends on both synthesis method and rare earth site magnetism. PDF analysis suggests that electronic phase separation occurs without accompanying chemical phase segregation or changes in the short-range crystallographic symmetry while synchrotron x-ray diffraction confirms that there is no change in the long-range crystallographic symmetry. X-ray absorption measurements confirm the $J_{eff}$ = 1/2 character of (Eu$_{1-x}$Ca$_{x}$)$_{2}$Ir$_{2}$O$_{7}$. Surprisingly these measurements also indicate a net electron doping, rather than the expected hole doping, indicating a compensatory mechanism. Lastly, XMCD measurements show a weak Ir magnetic polarization that is largely unaffected by Ca-doping.


title: Enabling Preserving Bisimulation Equivalence

abstract: Most fairness assumptions used for verifying liveness properties are criticised for being too strong or unrealistic. On the other hand, justness, arguably the minimal fairness assumption required for the verification of liveness properties, is not preserved by classical semantic equivalences, such as strong bisimilarity. To overcome this deficiency, we introduce a finer alternative to strong bisimilarity, called enabling preserving bisimilarity. We prove that this equivalence is justness-preserving and a congruence for all standard operators, including parallel composition.


title: Homotopy types of gauge groups over Riemann surfaces

abstract: Let $G$ be a compact connected Lie group with $\pi_1(G)\cong\mathbb{Z}$. We study the homotopy types of gauge groups of principal $G$-bundles over Riemann surfaces. This can be applied to an explicit computation of the homotopy groups of the moduli spaces of stable vector bundles over Riemann surfaces.


title: Personalized Stress Monitoring using Wearable Sensors in Everyday Settings

abstract: Since stress contributes to a broad range of mental and physical health problems, the objective assessment of stress is essential for behavioral and physiological studies. Although several studies have evaluated stress levels in controlled settings, objective stress assessment in everyday settings is still largely under-explored due to challenges arising from confounding contextual factors and limited adherence for self-reports. In this paper, we explore the objective prediction of stress levels in everyday settings based on heart rate (HR) and heart rate variability (HRV) captured via low-cost and easy-to-wear photoplethysmography (PPG) sensors that are widely available on newer smart wearable devices. We present a layered system architecture for personalized stress monitoring that supports a tunable collection of data samples for labeling, and present a method for selecting informative samples from the stream of real-time data for labeling. We captured the stress levels of fourteen volunteers through self-reported questionnaires over periods of between 1-3 months, and explored binary stress detection based on HR and HRV using Machine Learning Methods. We observe promising preliminary results given that the dataset is collected in the challenging environments of everyday settings. The binary stress detector is fairly accurate and can detect stressful vs non-stressful samples with a macro-F1 score of up to \%76. Our study lays the groundwork for more sophisticated labeling strategies that generate context-aware, personalized models that will empower health professionals to provide personalized interventions.


title: Manifold-Inspired Single Image Interpolation

abstract: Manifold models consider natural-image patches to be on a low-dimensional manifold embedded in a high dimensional state space and each patch and its similar patches to approximately lie on a linear affine subspace. Manifold models are closely related to semi-local similarity, a well-known property of natural images, referring to that for most natural-image patches, several similar patches can be found in its spatial neighborhood. Many approaches to single image interpolation use manifold models to exploit semi-local similarity by two mutually exclusive parts: i) searching each target patch's similar patches and ii) operating on the searched similar patches, the target patch and the measured input pixels to estimate the target patch. Unfortunately, aliasing in the input image makes it challenging for both parts. A very few works explicitly deal with those challenges and only ad-hoc solutions are proposed. To overcome the challenge in the first part, we propose a carefully-designed adaptive technique to remove aliasing in severely aliased regions, which cannot be removed from traditional techniques. This technique enables reliable identification of similar patches even in the presence of strong aliasing. To overcome the challenge in the second part, we propose to use the aliasing-removed image to guide the initialization of the interpolated image and develop a progressive scheme to refine the interpolated image based on manifold models. Experimental results demonstrate that our approach reconstructs edges with both smoothness along contours and sharpness across profiles, and achieves an average Peak Signal-to-Noise Ratio (PSNR) significantly higher than existing model-based approaches.


title: T$_k$ML-AP: Adversarial Attacks to Top-$k$ Multi-Label Learning

abstract: Top-$k$ multi-label learning, which returns the top-$k$ predicted labels from an input, has many practical applications such as image annotation, document analysis, and web search engine. However, the vulnerabilities of such algorithms with regards to dedicated adversarial perturbation attacks have not been extensively studied previously. In this work, we develop methods to create adversarial perturbations that can be used to attack top-$k$ multi-label learning-based image annotation systems (TkML-AP). Our methods explicitly consider the top-$k$ ranking relation and are based on novel loss functions. Experimental evaluations on large-scale benchmark datasets including PASCAL VOC and MS COCO demonstrate the effectiveness of our methods in reducing the performance of state-of-the-art top-$k$ multi-label learning methods, under both untargeted and targeted attacks.


title: Communication-avoiding micro-architecture to compute Xcorr scores for peptide identification

abstract: Database algorithms play a crucial part in systems biology studies by identifying proteins from mass spectrometry data. Many of these database search algorithms incur huge computational costs by computing similarity scores for each pair of sparse experimental spectrum and candidate theoretical spectrum vectors. Modern MS instrumentation techniques which are capable of generating high-resolution spectrometry data require comparison against an enormous search space, further emphasizing the need of efficient accelerators. Recent research has shown that the overall cost of scoring, and deducing peptides is dominated by the communication costs between different hierarchies of memory and processing units. However, these communication costs are seldom considered in accelerator-based architectures leading to inefficient DRAM accesses, and poor data-utilization due to irregular memory access patterns. In this paper, we propose a novel communication-avoiding micro-architecture to compute cross-correlation based similarity score by utilizing efficient local cache, and peptide pre-fetching to minimize DRAM accesses, and a custom-designed peptide broadcast bus to allow input reuse. An efficient bus arbitration scheme was designed, and implemented to minimize synchronization cost and exploit parallelism of processing elements. Our simulation results show that the proposed micro-architecture performs on average 24x better than a CPU implementation running on a 3.6 GHz Intel i7-4970 processor with 16GB memory.


title: Double Exotic States Productions in Pion and Kaon Induced Reactions

abstract: In the present work, we investigate the production of $P_{c}$ and $Z_c(3900)/Z_{cs}(3985)$ states in the processes of $\pi p \rightarrow Z_{c}(3900) P_{c} $/$Kp\rightarrow Z_{cs}(3985) P_{c}$ with an effective Lagrangian approach. Our estimations indicates that the cross sections depends on the model parameters and beam energies, in particular, the cross sections for $\pi p \to Z_c(3900) P_c(4312)$ and $\pi p \to Z_c(3900) P_c(4440)$ are similar and can reach up to 10 nb, while those for $\pi p \to Z_c(3900) P_c(4457)$ can be 30 nb. As for $Kp \to Z_{cs}(3985)P_c$, the cross sections are about 2.5 times smaller than those for $\pi p \to Z_c(3900) P_c$. The discussed processes in the present work, especially $\pi p \to Z_c(3900) P_c(4457)$, may be accessible by the the high-energy pion and kaon beams in the facilities of J-PARC and COMPASS.


title: Eavesdropping with Intelligent Reflective Surfaces: Threats and Defense Strategies

abstract: Intelligent reflecting surfaces (IRSs) have several prominent advantages, including improving the level of wireless communications security and privacy. In this work, we focus on this aspect and envision a strategy to counteract the presence of passive eavesdroppers overhearing transmissions from a base station towards legitimate users. Unlike most of the existing works addressing passive eavesdropping, the strategy we consider has low complexity and is suitable for scenarios where nodes are equipped with a limited number of antennas. Through our performance evaluation, we highlight the trade-off between the legitimate users' data rate and secrecy rate, and how the system parameters affect such a trade-off.


title: Deep Image-based Illumination Harmonization

abstract: Integrating a foreground object into a background scene with illumination harmonization is an important but challenging task in computer vision and augmented reality community. Existing methods mainly focus on foreground and background appearance consistency or the foreground object shadow generation, which rarely consider global appearance and illumination harmonization. In this paper, we formulate seamless illumination harmonization as an illumination exchange and aggregation problem. Specifically, we firstly apply a physically-based rendering method to construct a large-scale, high-quality dataset (named IH) for our task, which contains various types of foreground objects and background scenes with different lighting conditions. Then, we propose a deep image-based illumination harmonization GAN framework named DIH-GAN, which makes full use of a multi-scale attention mechanism and illumination exchange strategy to directly infer mapping relationship between the inserted foreground object and the corresponding background scene. Meanwhile, we also use adversarial learning strategy to further refine the illumination harmonization result. Our method can not only achieve harmonious appearance and illumination for the foreground object but also can generate compelling shadow cast by the foreground object. Comprehensive experiments on both our IH dataset and real-world images show that our proposed DIH-GAN provides a practical and effective solution for image-based object illumination harmonization editing, and validate the superiority of our method against state-of-the-art methods. Our IH dataset is available at https://github.com/zhongyunbao/Dataset.


title: An Empirical Analysis on Transparent Algorithmic Exploration in Recommender Systems

abstract: All learning algorithms for recommendations face inevitable and critical trade-off between exploiting partial knowledge of a user's preferences for short-term satisfaction and exploring additional user preferences for long-term coverage. Although exploration is indispensable for long success of a recommender system, the exploration has been considered as the risk to decrease user satisfaction. The reason for the risk is that items chosen for exploration frequently mismatch with the user's interests. To mitigate this risk, recommender systems have mixed items chosen for exploration into a recommendation list, disguising the items as recommendations to elicit feedback on the items to discover the user's additional tastes. This mix-in approach has been widely used in many recommenders, but there is rare research, evaluating the effectiveness of the mix-in approach or proposing a new approach for eliciting user feedback without deceiving users. In this work, we aim to propose a new approach for feedback elicitation without any deception and compare our approach to the conventional mix-in approach for evaluation. To this end, we designed a recommender interface that reveals which items are for exploration and conducted a within-subject study with 94 MTurk workers. Our results indicated that users left significantly more feedback on items chosen for exploration with our interface. Besides, users evaluated that our new interface is better than the conventional mix-in interface in terms of novelty, diversity, transparency, trust, and satisfaction. Finally, path analysis show that, in only our new interface, exploration caused to increase user-centric evaluation metrics. Our work paves the way for how to design an interface, which utilizes learning algorithm based on users' feedback signals, giving better user experience and gathering more feedback data.


title: To adjust or not to adjust? Estimating the average treatment effect in randomized experiments with missing covariates

abstract: Complete randomization allows for consistent estimation of the average treatment effect based on the difference in means of the outcomes without strong modeling assumptions on the outcome-generating process. Appropriate use of the pretreatment covariates can further improve the estimation efficiency. However, missingness in covariates is common in experiments and raises an important question: should we adjust for covariates subject to missingness, and if so, how? The unadjusted difference in means is always unbiased. The complete-covariate analysis adjusts for all completely observed covariates and improves the efficiency of the difference in means if at least one completely observed covariate is predictive of the outcome. Then what is the additional gain of adjusting for covariates subject to missingness? A key insight is that the missingness indicators act as fully observed pretreatment covariates as long as missingness is not affected by the treatment, and can thus be used in covariate adjustment to bring additional estimation efficiency. This motivates adding the missingness indicators to the regression adjustment, yielding the missingness-indicator method as a well-known but not so popular strategy in the literature of missing data. We recommend it due to its many advantages. We also propose modifications to the missingness-indicator method based on asymptotic and finite-sample considerations. To reconcile the conflicting recommendations in the missing data literature, we analyze and compare various strategies for analyzing randomized experiments with missing covariates under the design-based framework. This framework treats randomization as the basis for inference and does not impose any modeling assumptions on the outcome-generating process and missing-data mechanism.


title: Dynamic Virtual Power Plant: A New Concept for Grid Integration of Renewable Energy Sources

abstract: The notion of Virtual Power Plant (VPP) has been used many times in last years in power systems and for several reasons. As a general trend, the behavior of a classic synchronous generator is to be emulated for a class of conventional grid components like, e.g., renewable generators or/and power electronic units. Most of the times production of these units is of interest, as it is the case for the new AGC scheme of Spain which, from this point of view, looks like a VPP. However, dynamic aspects are of high importance, especially for increasing the actual rate of penetration of Renewable Energy Sources (RES). Indeed, to go above the actual rate of RES penetration, one should deal with full participation of RES to grid services. For that, we propose here a new concept called Dynamic VPP (DVPP) which fully integrates the dynamic aspects at all levels: locally (for each RES generator), globally (for grid ancillary services and interaction with other neighbor elements of the grid) and economically (for internal optimal dispatch and participation to electricity markets). A DVPP is a set of RES along with a set of control and operation procedures. This means methodologies for: choosing the participating RES, optimal and continuous operation as a whole (especially in case of loss of natural resources - e.g., wind, sun - on a part of the DVPP), regulation (in the dynamic sense) to ensure local objectives for each generator, participation to ancillary services of the DVPP as a unit and to diminish negative effects of interaction with neighbor dynamics elements of the power system, integration in both actual power systems scenarios (with mixed classic and power electronics based generation) and future ones with high degree of RES penetration.