title: Entanglement phase transitions in measurement-only dynamics

abstract: Unitary circuits subject to repeated projective measurements can undergo an entanglement phase transition (EPT) as a function of the measurement rate. This transition is generally understood in terms of a competition between the scrambling effects of unitary dynamics and the disentangling effects of measurements. We find that, surprisingly, EPTs are possible even in the absence of scrambling unitary dynamics, where they are best understood as arising from measurements alone. This motivates us to introduce \emph{measurement-only models}, in which the "scrambling" and "un-scrambling" effects driving the EPT are fundamentally intertwined and cannot be attributed to physically distinct processes. This represents a novel form of an EPT, conceptually distinct from that in hybrid unitary-projective circuits. We explore the entanglement phase diagrams, critical points, and quantum code properties of some of these measurement-only models. We find that the principle driving the EPTs in these models is \emph{frustration}, or mutual incompatibility, of the measurements. Suprisingly, an entangling (volume-law) phase is the generic outcome when measuring sufficiently long but still local ($\gtrsim 3$-body) operators. We identify a class of exceptions to this behavior ("bipartite ensembles") which cannot sustain an entangling phase, but display dual area-law phases, possibly with different kinds of quantum order, separated by self-dual critical points. Finally, we introduce a measure of information spreading in dynamics with measurements and use it to demonstrate the emergence of a statistical light-cone, despite the non-locality inherent to quantum measurements.


title: Stabilizing Quasi-Time-Optimal Nonlinear Model Predictive Control with Variable Discretization

abstract: This paper deals with the development and analysis of novel time-optimal point-to-point model predictive control concepts for nonlinear systems. Recent approaches in the literature apply a time transformation, however, which do not maintain recursive feasibility for piecewise constant control parameterization. The key idea in this paper is to introduce uniform grids with variable discretization. A shrinking-horizon grid adaptation scheme ensures convergence to a specific region around the target state and recursive feasibility. The size of the region is configurable by design parameters. This facilitates the systematic dual-mode design for quasi-time-optimal control to restore asymptotic stability and establish a smooth stabilization. Two nonlinear program formulations with different sparsity patterns are introduced to realize and implement the underlying optimal control problem. For a class of numerical integration schemes, even nominal asymptotic stability and true time-optimality are achieved without dual-mode. A comparative analysis as well as experimental results demonstrate the effectiveness of the proposed techniques.


title: Fermi-liquid behaviors in three-dimensional polaron system with tunable dispersion

abstract: We investigate the electronic properties and electron correlations in a three-dimensional (3D) polaron system with tunable dispersion parameter in Fermi-liquid picture. The polaronic coupling is considered as a weakly attractive bare coupling. Both the single and many-polaron, zero and finite-temperature cases are considered in this paper. The physical quantities, including the self-energy, spectral function, optical conductivity (both the planar and 3D one), free energy, and momentum distribution function, are studied and are been verified exhibit Fermi-liquid features. Both the analytic derivation and numerical calculation are done in this paper to make our results reliable.


title: Learning Entangled Single-Sample Distributions via Iterative Trimming

abstract: In the setting of entangled single-sample distributions, the goal is to estimate some common parameter shared by a family of distributions, given one \emph{single} sample from each distribution. We study mean estimation and linear regression under general conditions, and analyze a simple and computationally efficient method based on iteratively trimming samples and re-estimating the parameter on the trimmed sample set. We show that the method in logarithmic iterations outputs an estimation whose error only depends on the noise level of the $\lceil \alpha n \rceil$-th noisiest data point where $\alpha$ is a constant and $n$ is the sample size. This means it can tolerate a constant fraction of high-noise points. These are the first such results for the method under our general conditions. It also justifies the wide application and empirical success of iterative trimming in practice. Our theoretical results are complemented by experiments on synthetic data.


title: Driven and active colloids at fluid interfaces

abstract: We derive expressions for the leading-order far-field flows generated by externally driven and active (swimming) colloids at planar fluid-fluid interfaces. We consider colloids adjacent to the interface or adhered to the interface with a pinned contact line. The Reynolds and capillary numbers are assumed much less than unity, in line with typical micron-scale colloids involving air- or alkane-aqueous interfaces. For driven colloids, the leading-order flow is given by the point-force (and/or torque) response of this system. For active colloids, the force-dipole (stresslet) response occurs at leading order. At clean (surfactant-free) interfaces, these hydrodynamic modes are essentially a restricted set of the usual Stokes multipoles in a bulk fluid. To leading order, driven colloids exert Stokeslets parallel to the interface, while active colloids drive differently oriented stresslets depending on the colloid's orientation. We then consider how these modes are altered by the presence of an incompressible interface, a typical circumstance for colloidal systems at small capillary numbers in the presence of surfactant. The leading-order modes for driven and active colloids are restructured dramatically. For driven colloids, interfacial incompressibility substantially weakens the far-field flow normal to the interface; the point-force response drives flow only parallel to the interface. However, Marangoni stresses induce a new dipolar mode, which lacks an analogue on a clean interface. Surface-viscous stresses, if present, potentially generate very long-ranged flow on the interface and the surrounding fluids. Our results have important implications for colloid assembly and advective mass transport enhancement near fluid boundaries.


title: Sparse aNETT for Solving Inverse Problems with Deep Learning

abstract: We propose a sparse reconstruction framework (aNETT) for solving inverse problems. Opposed to existing sparse reconstruction techniques that are based on linear sparsifying transforms, we train an autoencoder network $D \circ E$ with $E$ acting as a nonlinear sparsifying transform and minimize a Tikhonov functional with learned regularizer formed by the $\ell^q$-norm of the encoder coefficients and a penalty for the distance to the data manifold. We propose a strategy for training an autoencoder based on a sample set of the underlying image class such that the autoencoder is independent of the forward operator and is subsequently adapted to the specific forward model. Numerical results are presented for sparse view CT, which clearly demonstrate the feasibility, robustness and the improved generalization capability and stability of aNETT over post-processing networks.


title: Toward a unified equation of state for multi-messenger astronomy

abstract: We present a first step in developing a benchmark equation-of-state (EoS) model for multi-messenger astronomy that unifies the thermodynamics of quark and hadronic degrees of freedom. A Lagrangian approach to the thermodynamic potential of quark-meson-nucleon (QMN) matter was used. In this approach, dynamical chiral-symmetry breaking is described by the scalar mean-field dynamics coupled to quarks and nucleons and their chiral partners, whereby its restoration occurs in the hadronic phase by parity doubling, as well as in the quark phase. Quark confinement was achieved by an auxiliary scalar field that parametrizes a dynamical infrared cutoff in the quark sector, serving as an ultraviolet cutoff for the nucleonic phase space. The gap equations were solved for the isospin-symmetric case, as well as for neutron star (NS) conditions. We also calculated the mass-radius (MR) relation of NSs and their tidal deformability parameter. The obtained EoS is in accordance with nuclear matter properties at saturation density and with the flow constraint from heavy ion collision experiments. For isospin-asymmetric matter, a sequential occurrence of light quark flavors is obtained, allowing for a mixed phase of chirally-symmetric nucleonic matter with deconfined down quarks. The MR relations and TDs for compact stars fulfill the constraints from the latest astrophysical observations for PSR J0740+6620, PSR J0030+0451, and the NS merger GW170817, whereby the tension between the maximum mass and compactness constraints rather uniquely fixes the model parameters. The model predicts the existence of stars with a core of chirally restored but purely hadronic (confined) matter for masses beyond $1.8~M_\odot$. Stars with pure-quark matter cores are found to be unstable against the gravitational collapse. This instability is shifted to even higher densities if repulsive interactions between quarks are included.


title: Turbulence and Energetic Particles in Radiative Shock Waves in the Cygnus Loop I: Shock Properties

abstract: We have obtained a contiguous set of long-slit spectra of a shock wave in the Cygnus Loop to investigate its structure, which is far from the morphology predicted by 1D models. Proper motions from Hubble Space Telescope images combined with the known distance to the Cygnus Loop provide an accurate shock speed. Earlier analyses of shock spectra estimated the shock speed, postshock density, temperature, and elemental abundances. In this paper we determine several more shock parameters: a more accurate shock speed, ram pressure, density, compression ratio, dust destruction efficiency, magnetic field strength, and vorticity in the cooling region. From the derived shock properties we estimate the emissivities of synchrotron emission in the radio and pion decay emission in the gamma rays. Both are consistent with the observations if we assume simple adiabatic compression of ambient cosmic rays as in the van der Laan mechanism. We also find that, although the morphology is far from that predicted by 1D models and the line ratios vary dramatically from point to point, the average spectrum is matched reasonably well by 1D shock models with the shock speed derived from the measured proper motion. A subsequent paper will analyze the development of turbulence in the cooling zone behind the shock.


title: Renormalized entanglement entropy and curvature invariants

abstract: Renormalized entanglement entropy can be defined using the replica trick for any choice of renormalization scheme; renormalized entanglement entropy in holographic settings is expressed in terms of renormalized areas of extremal surfaces. In this paper we show how holographic renormalized entanglement entropy can be expressed in terms of the Euler invariant of the surface and renormalized curvature invariants. For a spherical entangling region in an odd-dimensional CFT, the renormalized entanglement entropy is proportional to the Euler invariant of the holographic entangling surface, with the coefficient of proportionality capturing the (renormalized) F quantity. Variations of the entanglement entropy can be expressed elegantly in terms of renormalized curvature invariants, facilitating general proofs of the first law of entanglement.


title: Neural network compression via learnable wavelet transforms

abstract: Wavelets are well known for data compression, yet have rarely been applied to the compression of neural networks. This paper shows how the fast wavelet transform can be used to compress linear layers in neural networks. Linear layers still occupy a significant portion of the parameters in recurrent neural networks (RNNs). Through our method, we can learn both the wavelet bases and corresponding coefficients to efficiently represent the linear layers of RNNs. Our wavelet compressed RNNs have significantly fewer parameters yet still perform competitively with the state-of-the-art on synthetic and real-world RNN benchmarks. Wavelet optimization adds basis flexibility, without large numbers of extra weights. Source code is available at https://github.com/v0lta/Wavelet-network-compression.


title: Optical cycling, radiative deflection and laser cooling of barium monohydride (BaH)

abstract: We present the first experimental demonstration of radiation pressure force deflection and direct laser cooling for barium monohydride (BaH) molecules resulting from multiple photon scattering. Despite a small recoil velocity (2.7 mm/s) and a long excited state lifetime (137 ns), we use 1060 nm laser light exciting the $X\rightarrow A$ electronic transition of BaH to deflect a cryogenic buffer-gas beam and reduce its transverse velocity spread. Multiple experimental methods are employed to characterize the optical cycling dynamics and benchmark theoretical estimates based on rate equation models as well as solutions of the Lindblad master equation for the complete multilevel system. Broader implications for laser cooling and magneto-optical trapping of heavy-metal-containing molecules with narrow transition linewidths are presented. Our results pave the way for producing a new class of ultracold molecules -- alkaline earth monohydrides -- via direct laser cooling and trapping, opening the door to realizing a new method for delivering ultracold hydrogen atoms (Lane 2015 \textit{Phys. Rev. A} 92, 022511).


title: Language-agnostic Multilingual Modeling

abstract: Multilingual Automated Speech Recognition (ASR) systems allow for the joint training of data-rich and data-scarce languages in a single model. This enables data and parameter sharing across languages, which is especially beneficial for the data-scarce languages. However, most state-of-the-art multilingual models require the encoding of language information and therefore are not as flexible or scalable when expanding to newer languages. Language-independent multilingual models help to address this issue, and are also better suited for multicultural societies where several languages are frequently used together (but often rendered with different writing systems). In this paper, we propose a new approach to building a language-agnostic multilingual ASR system which transforms all languages to one writing system through a many-to-one transliteration transducer. Thus, similar sounding acoustics are mapped to a single, canonical target sequence of graphemes, effectively separating the modeling and rendering problems. We show with four Indic languages, namely, Hindi, Bengali, Tamil and Kannada, that the language-agnostic multilingual model achieves up to 10% relative reduction in Word Error Rate (WER) over a language-dependent multilingual model.


title: Dark Matter properties through cosmic history

abstract: We perform the first test of dark matter (DM) stress-energy evolution through cosmic history, using cosmic microwave background measurements supplemented with baryon acoustic oscillation data and the Hubble Space Telescope key project data. We constrain the DM equation of state (EoS) in 8 redshift bins, and its sound speed and (shear) viscosity in 9 redshift bins, finding no convincing evidence for non-$\Lambda$CDM values in any of the redshift bins. Despite this enlarged parameter space, the sound speed and viscosity are constrained relatively well at late times (due to the inclusion of CMB lensing), whereas the EoS is most strongly constrained around recombination. These results constrain for the first time the level of "coldness" required of DM across various cosmological epochs at both the background and perturbative levels. We show that simultaneously allowing time dependence for both the EoS and sound speed parameters shifts the posterior of the DM abundance before recombination to a higher value, while keeping the present day DM abundance similar to the $\Lambda$CDM value. This shifts the posterior for the present day Hubble constant compared to $\Lambda$CDM, suggesting that DM with time-dependent parameters is well-suited to explore possible solutions to persistent tensions within the $\Lambda$CDM model. We perform a detailed comparison with our previous study involving a vanishing sound speed and viscosity using the same datasets in order to explain the physical mechanism behind these shifts.


title: Utilizing Mask R-CNN for Waterline Detection in Canoe Sprint Video Analysis

abstract: Determining a waterline in images recorded in canoe sprint training is an important component for the kinematic parameter analysis to assess an athlete's performance. Here, we propose an approach for the automated waterline detection. First, we utilized a pre-trained Mask R-CNN by means of transfer learning for canoe segmentation. Second, we developed a multi-stage approach to estimate a waterline from the outline of the segments. It consists of two linear regression stages and the systematic selection of canoe parts. We then introduced a parameterization of the waterline as a basis for further evaluations. Next, we conducted a study among several experts to estimate the ground truth waterlines. This not only included an average waterline drawn from the individual experts annotations but, more importantly, a measure for the uncertainty between individual results. Finally, we assessed our method with respect to the question whether the predicted waterlines are in accordance with the experts annotations. Our method demonstrated a high performance and provides opportunities for new applications in the field of automated video analysis in canoe sprint.


title: A Note on Load Balancing in Many-Server Heavy-Traffic Regime

abstract: In this note, we apply Stein's method to analyze the performance of general load balancing schemes in the many-server heavy-traffic regime. In particular, consider a load balancing system of $N$ servers and the distance of arrival rate to the capacity region is given by $N^{1-\alpha}$ with $\alpha > 1$. We are interested in the performance as $N$ goes to infinity under a large class of policies. We establish different asymptotics under different scalings and conditions. Specifically, (i) If the second moments linearly increase with $N$ with coefficients $\sigma_a^2$ and $\nu_s^2$, then for any $\alpha > 4$, the distribution of the sum queue length scaled by $N^{-\alpha}$ converges to an exponential random variable with mean $\frac{\sigma_a^2 + \nu_s^2}{2}$. (3) If the second moments quadratically increase with $N$ with coefficients $\tilde{\sigma}_a^2$ and $\tilde{\nu}_s^2$, then for any $\alpha > 3$, the distribution of the sum queue length scaled by $N^{-\alpha-1}$ converges to an exponential random variable with mean $\frac{\tilde{\sigma}_a^2 + \tilde{\nu}_s^2}{2}$. Both results are simple applications of our previously developed framework of Stein's method for heavy-traffic analysis in \cite{zhou2020note}.


title: Predicting nucleation near the spinodal in the Ising model using machine learning

abstract: We use a Convolutional Neural Network (CNN) and two logistic regression models to predict the probability of nucleation in the two-dimensional Ising model. The three models successfully predict the probability for the Nearest Neighbor Ising model for which classical nucleation is observed. The CNN outperforms the logistic regression models near the spinodal of the Long Range Ising model, but the accuracy of its predictions decreases as the quenches approach the spinodal. Occlusion analysis suggests that this decrease is due to the vanishing difference between the density of the nucleating droplet and the background. Our results are consistent with the general conclusion that predictability decreases near a critical point.


title: LSQ+: Improving low-bit quantization through learnable offsets and better initialization

abstract: Unlike ReLU, newer activation functions (like Swish, H-swish, Mish) that are frequently employed in popular efficient architectures can also result in negative activation values, with skewed positive and negative ranges. Typical learnable quantization schemes [PACT, LSQ] assume unsigned quantization for activations and quantize all negative activations to zero which leads to significant loss in performance. Naively using signed quantization to accommodate these negative values requires an extra sign bit which is expensive for low-bit (2-, 3-, 4-bit) quantization. To solve this problem, we propose LSQ+, a natural extension of LSQ, wherein we introduce a general asymmetric quantization scheme with trainable scale and offset parameters that can learn to accommodate the negative activations. Gradient-based learnable quantization schemes also commonly suffer from high instability or variance in the final training performance, hence requiring a great deal of hyper-parameter tuning to reach a satisfactory performance. LSQ+ alleviates this problem by using an MSE-based initialization scheme for the quantization parameters. We show that this initialization leads to significantly lower variance in final performance across multiple training runs. Overall, LSQ+ shows state-of-the-art results for EfficientNet and MixNet and also significantly outperforms LSQ for low-bit quantization of neural nets with Swish activations (e.g.: 1.8% gain with W4A4 quantization and upto 5.6% gain with W2A2 quantization of EfficientNet-B0 on ImageNet dataset). To the best of our knowledge, ours is the first work to quantize such architectures to extremely low bit-widths.


title: Emergent conformal symmetry in non-unitary random dynamics of free fermions

abstract: We present random quantum circuit models for non-unitary quantum dynamics of free fermions in one spatial dimension. Numerical simulations reveal that the dynamics tends towards steady states with logarithmic violations of the entanglement area law and power law correlation functions. Moreover, starting with a short-range entangled many-body state, the dynamical evolution of entanglement and correlations quantitatively agrees with the predictions of two-dimensional conformal field theory with a space-like time direction. We argue that this behavior is generic in non-unitary free quantum dynamics with time-dependent randomness, and show that the emergent conformal dynamics of two-point functions arises out of a simple "nonlinear master equation".

categories: quant-ph cond-mat.dis-nn cond-mat.stat-mech nlin.AO


title: CLOPS: Continual Learning of Physiological Signals

abstract: Deep learning algorithms are known to experience destructive interference when instances violate the assumption of being independent and identically distributed (i.i.d). This violation, however, is ubiquitous in clinical settings where data are streamed temporally and from a multitude of physiological sensors. To overcome this obstacle, we propose CLOPS, a replay-based continual learning strategy. In three continual learning scenarios based on three publically-available datasets, we show that CLOPS can outperform the state-of-the-art methods, GEM and MIR. Moreover, we propose end-to-end trainable parameters, which we term task-instance parameters, that can be used to quantify task difficulty and similarity. This quantification yields insights into both network interpretability and clinical applications, where task difficulty is poorly quantified.


title: Sparse Oblique Decision Tree for Power System Security Rules Extraction and Embedding

abstract: Increasing the penetration of variable generation has a substantial effect on the operational reliability of power systems. The higher level of uncertainty that stems from this variability makes it more difficult to determine whether a given operating condition will be secure or insecure. Data-driven techniques provide a promising way to identify security rules that can be embedded in economic dispatch model to keep power system operating states secure. This paper proposes using a sparse weighted oblique decision tree to learn accurate, understandable, and embeddable security rules that are linear and can be extracted as sparse matrices using a recursive algorithm. These matrices can then be easily embedded as security constraints in power system economic dispatch calculations using the Big-M method. Tests on several large datasets with high renewable energy penetration demonstrate the effectiveness of the proposed method. In particular, the sparse weighted oblique decision tree outperforms the state-of-art weighted oblique decision tree while keeping the security rules simple. When embedded in the economic dispatch, these rules significantly increase the percentage of secure states and reduce the average solution time.


title: Parameter estimation of path-dependent McKean-Vlasov stochastic differential equations

abstract: The work concerns a class of path-dependent McKean-Vlasov stochastic differential equations with unknown parameters. First, we prove the existence and uniqueness of these equations under non-Lipschitz conditions. Second, we construct maximum likelihood estimators of these parameters and then discuss their strong consistency. Third, a numerical simulation method for the class of path-dependent McKean-Vlasov stochastic differential equations is offered. Moreover, we estimate the errors between solutions of these equations and that of their numerical equations. Finally, we give an example to explain our result.


title: Human-Collective Collaborative Site Selection

abstract: Robotic collectives are large groups (at least 50) of locally sensing and communicating robots that encompass characteristics of swarms and colonies, whose emergent behaviors accomplish complex tasks. Future human-collective teams will extend the ability of operators to monitor, respond, and make decisions in disaster response, search and rescue, and environmental monitoring problems. This manuscript evaluates two collective best-of-n decision models for enabling collectives to identify and choose the highest valued target from a finite set of n targets. Two challenges impede the future use of human-collective shared decisions: 1) environmental bias reduces collective decision accuracy when poorer targets are easier to evaluate than higher quality targets, and 2) little is understood about shared human-collective decision making interaction strategies. The two evaluated collective best-of-n models include an existing insect colony decision model and an extended bias-reducing model that attempts to reduce environmental bias in order to improve accuracy. Collectives using these two strategies are compared independently and as members of human-collective teams. Independently, the extended model is slower than the original model, but the extended algorithm is 57% more accurate in decisions where the optimal option is more difficult to evaluate. Human-collective teams using the bias-reducing model require less operator influence and achieve 25% higher accuracy with difficult decisions, than the human-collective teams using the original model. Further, a novel human-collective interaction strategy enables operators to adjust collective autonomy while making multiple simultaneous decisions.


title: Finite size effects in critical fiber networks

abstract: Fibrous networks such as collagen are common in physiological systems. One important function of these networks is to provide mechanical stability for cells and tissues. At physiological levels of connectivity, such networks would be mechanically unstable with only central-force interactions. While networks can be stabilized by bending interactions, it has also been shown that they exhibit a critical transition from floppy to rigid as a function of applied strain. Beyond a certain strain threshold, it is predicted that underconstrained networks with only central-force interactions exhibit a discontinuity in the shear modulus. We study the finite-size scaling behavior of this transition and identify both the mechanical discontinuity and critical exponents in the thermodynamic limit. We find both non-mean-field behavior and evidence for a hyperscaling relation for the critical exponents, for which the network stiffness is analogous to the heat capacity for thermal phase transitions. Further evidence for this is also found in the self-averaging properties of fiber networks.


title: FlashFlow: A Secure Speed Test for Tor

abstract: The Tor network uses a measurement system to estimate its relays' forwarding capacity and to balance traffic among them. This system has been shown to be vulnerable to adversarial manipulation. Moreover, its accuracy and effectiveness in benign circumstances has never been fully quantified. We first obtain such a quantification by analyzing Tor metrics data and performing experiments on the live network. Our results show that Tor currently underestimates its true capacity by about 50% and improperly balances its traffic by 15-25%. Then, to solve the problems with security and accuracy, we present FlashFlow, a system to measure the capacity of Tor relays. Our analysis shows that FlashFlow limits a malicious relay to obtaining a capacity estimate at most 1.33 times its true capacity. Through realistic Internet experiments, we find that FlashFlow measures relay capacity with at least 89% accuracy 95% of the time. Through simulation, we find that FlashFlow can measure the entire Tor network in less than 5 hours using 3 measurers with 1 Gbit/s of bandwidth each. Finally, simulations using FlashFlow for load balancing shows that, compared to TorFlow, network weight error decreases by 86%, while the median of 50 KiB, 1 MiB, and 5 MiB transfer times decreases by 15%, 29%, and 37%, respectively. Moreover, FlashFlow yields more consistent client performance: the median rate of transfer timeouts decreases by 100%, while the standard deviation of 50 KiB, 1 MiB, and 5 MiB transfer times decreases by 55%, 61%, and 41%, respectively. We also find that the performance improvements increase relative to TorFlow as the total client-traffic load increases, demonstrating that FlashFlow is better suited to supporting network growth.


title: ViSQOL v3: An Open Source Production Ready Objective Speech and Audio Metric

abstract: Estimation of perceptual quality in audio and speech is possible using a variety of methods. The combined v3 release of ViSQOL and ViSQOLAudio (for speech and audio, respectively,) provides improvements upon previous versions, in terms of both design and usage. As an open source C++ library or binary with permissive licensing, ViSQOL can now be deployed beyond the research context into production usage. The feedback from internal production teams at Google has helped to improve this new release, and serves to show cases where it is most applicable, as well as to highlight limitations. The new model is benchmarked against real-world data for evaluation purposes. The trends and direction of future work is discussed.


title: Natural gas shortages during the "coal-to-gas" transition in China have caused a large redistribution of air pollution

abstract: The Chinese "coal-to-gas" strategy aims at reducing coal consumption and related air pollution by promoting the use of clean and low carbon fuels in northern China. Here we show that on top of meteorological influences, these measures achieved an average decrease of fine particulate matter (PM2.5) concentrations of ~14% during winter 2017 in Beijing and surrounding areas (the "2+26" pilot cities). However, the localized air quality improvement was accompanied by a contemporaneous ~15% upsurge of PM2.5 concentrations over large areas in southern China. We find that the pollution transfer that resulted from a shift in emissions was caused by a natural gas shortage in the south due to the "coal-to-gas" transition in the north. The overall shortage of natural gas greatly jeopardized the air quality benefits of the "coal-to-gas" strategy in winter 2017 and reflects structural challenges and potential threats in China's clean energy transition.


title: Unusually High CO Abundance of the First Active Interstellar Comet

abstract: Comets spend most of their lives at large distances from any star, during which time their interior compositions remain relatively unaltered. Cometary observations can therefore provide direct insight into the chemistry that occurred during their birth at the time of planet formation. To-date, there have been no confirmed observations of parent volatiles (gases released directly from the nucleus) of a comet from any planetary system other than our own. Here we present high-resolution, interferometric observations of 2I/Borisov, the first confirmed interstellar comet, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) on 15th-16th December 2019. Our observations reveal emission from hydrogen cyanide (HCN), and carbon monoxide (CO), coincident with the expected position of 2I/Borisov's nucleus, with production rates Q(HCN)=$(7.0\pm1.1)\times10^{23}$ s$^{-1}$ and Q(CO)=$(4.4\pm0.7)\times10^{26}$ s$^{-1}$. While the HCN abundance relative to water (0.06-0.16%) appears similar to that of typical, previously observed comets in our Solar System, the abundance of CO (35-105%) is among the highest observed in any comet within 2 au of the Sun. This shows that 2I/Borisov must have formed in a relatively CO-rich environment - probably beyond the CO ice-line in the very cold, outer regions of a distant protoplanetary accretion disk, as part of a population of small, icy bodies analogous to our Solar System's own proto-Kuiper Belt.


title: On the Symmetries of Cosmological Perturbations

abstract: The space of inflationary models is vast, containing wide varieties of mechanisms, symmetries, and spectra of particles. Consequently, the space of observational signatures is similarly complex. Hence, it is natural to look for boundaries of the space of models and their signatures. In this paper, we explore the possible symmetries associated with the primordial cosmological perturbations and their correlators in the asymptotic future. Assuming the observed homogeneity, isotropy and (approximate) scale invariance, we prove three main results. First, correlation functions of scalar metric fluctuations are uniquely characterized by soft theorems and are free from ambiguity under field redefinitions. Second, whatever the particle content and interactions, when the standard soft theorems apply, invariance under de Sitter boosts (linearly realized conformal invariance) is only possible if all connected correlators vanish identically, i.e. if the theory is free. Third, conformal invariance is the largest set of linearly realized (bosonic) symmetries of the correlators of any single scalar, irrespectively of any soft theorems or particle content.


title: On The Problem of Relevance in Statistical Inference

abstract: This paper is dedicated to the "50 Years of the Relevance Problem" - a long-neglected topic that begs attention from practical statisticians who are concerned with the problem of drawing inference from large-scale heterogeneous data.


title: Weighted Cheeger and Buser Inequalities, with Applications to Clustering and Cutting Probability Densities

abstract: In this paper, we show how sparse or isoperimetric cuts of a probability density function relate to Cheeger cuts of its principal eigenfunction, for appropriate definitions of `sparse cut' and `principal eigenfunction'. We construct these appropriate definitions of sparse cut and principal eigenfunction in the probability density setting. Then, we prove Cheeger and Buser type inequalities similar to those for the normalized graph Laplacian of Alon-Milman. We demonstrate that no such inequalities hold for most prior definitions of sparse cut and principal eigenfunction. We apply this result to generate novel algorithms for cutting probability densities and clustering data, including a principled variant of spectral clustering.


title: Almost-Reed--Muller Codes Achieve Constant Rates for Random Errors

abstract: This paper considers '$\delta$-almost Reed-Muller codes', i.e., linear codes spanned by evaluations of all but a $\delta$ fraction of monomials of degree at most $d$. It is shown that for any $\delta > 0$ and any $\varepsilon>0$, there exists a family of $\delta$-almost Reed-Muller codes of constant rate that correct $1/2-\varepsilon$ fraction of random errors with high probability. For exact Reed-Muller codes, the analogous result is not known and represents a weaker version of the longstanding conjecture that Reed-Muller codes achieve capacity for random errors (Abbe-Shpilka-Wigderson STOC '15). Our approach is based on the recent polarization result for Reed-Muller codes, combined with a combinatorial approach to establishing inequalities between the Reed-Muller code entropies.


title: Semi-closed form prices of barrier options in the Hull-White model

abstract: In this paper we derive semi-closed form prices of barrier (perhaps, time-dependent) options for the Hull-White model, ie., where the underlying follows a time-dependent OU process with a mean-reverting drift. Our approach is similar to that in (Carr and Itkin, 2020) where the method of generalized integral transform is applied to pricing barrier options in the time-dependent OU model, but extends it to an infinite domain (which is an unsolved problem yet). Alternatively, we use the method of heat potentials for solving the same problems. By semi-closed solution we mean that first, we need to solve numerically a linear Volterra equation of the first kind, and then the option price is represented as a one-dimensional integral. Our analysis shows that computationally our method is more efficient than the backward and even forward finite difference methods (if one uses them to solve those problems), while providing better accuracy and stability.


title: Black-Box Strategies and Equilibrium for Games with Cumulative Prospect Theoretic Players

abstract: The betweenness property of preference relations states that a probability mixture of two lotteries should lie between them in preference. It is a weakened form of the independence property and hence satisfied in expected utility theory (EUT). Experimental violations of betweenness are well-documented and several preference theories, notably cumulative prospect theory (CPT), do not satisfy betweenness. We prove that CPT preferences satisfy betweenness if and only if they conform with EUT preferences. In game theory, lack of betweenness in the players' preference relations makes it essential to distinguish between the two interpretations of a mixed action by a player - conscious randomizations by the player and the uncertainty in the beliefs of the opponents. We elaborate on this distinction and study its implication for the definition of Nash equilibrium. This results in four different notions of equilibrium, with pure and mixed action Nash equilibrium being two of them. We dub the other two pure and mixed black-box strategy Nash equilibrium respectively. We resolve the issue of existence of such equilibria and examine how these different notions of equilibrium compare with each other.


title: The Conformational Space of a Flexible Amino Acid at Metallic Surfaces

abstract: In interfaces between inorganic and biological materials relevant for technological applications, the general challenge of structure determination is exacerbated by the high flexibility of bioorganic components, chemical bonding, and charge rearrangement at the interface. In this paper, we investigate a chemically complex building block, namely, the arginine (Arg) amino-acid interfaced with Cu, Ag and Au (111) surfaces. We investigate how the environment changes the accessible conformational space of this amino acid, by building and analyzing a database of thousands of structures optimized with the PBE functional including screened pairwise van der Waals interactions. When in contact with metallic surfaces, the accessible space for Arg is dramatically reduced, while the one for Arg-H$^+$ is instead increased if compared to the gas-phase. This is explained by the formation of strong bonds between Arg and the surfaces and by their absence and charge screening on Arg-H$^+$ upon adsorption. We also observe protonation-dependent stereoselective binding of the amino acid to the metal surfaces: Arg adsorbs with its chiral C$_\alpha$H center pointing H away from the surfaces while Arg-H$^+$ adsorbs with H pointing toward the surface.


title: Lindemann unjamming of emulsions

abstract: We study the bulk and shear elastic properties of barely-compressed, "athermal" emulsions and find that the rigidity of the jammed solid fails at remarkably large critical osmotic pressures. The minuscule yield strain and similarly small Brownian particle displacement of solid emulsions close to this transition suggests that this catastrophic failure corresponds to a plastic-entropic instability: the solid becomes too soft and weak to resist the thermal agitation of the droplets that compose it and fails. We propose a modified Lindemann stability criterion to describe this transition and derive a scaling law for the critical osmotic pressure that agrees quantitatively with experimental observations.


title: Moreau-Yosida regularization for optimal control of fractional PDEs with state constraints: parabolic case

abstract: This paper considers optimal control of fractional parabolic PDEs with both state and control constraints. The key challenge is how to handle the state constraints. Similarly, to the elliptic case, in this paper, we establish several new mathematical tools in the parabolic setting that are of wider interest. For example, existence of solution to the fractional parabolic equation with measure data on the right-hand-side. We employ the Moreau-Yosida regularization to handle the state constraints. We establish convergence, with rate, of the regularized optimal control problem to the original one. Numerical experiments confirm what we have proven theoretically.


title: On-the-fly Detection of User Engagement Decrease in Spontaneous Human-Robot Interaction, International Journal of Social Robotics, 2019

abstract: In this paper, we consider the detection of a decrease of engagement by users spontaneously interacting with a socially assistive robot in a public space. We first describe the UE-HRI dataset that collects spontaneous Human-Robot Interactions following the guidelines provided by the Affective Computing research community to collect data "in-the-wild". We then analyze the users' behaviors, focusing on proxemics, gaze, head motion, facial expressions and speech during interactions with the robot. Finally, we investigate the use of deep learning techniques (Recurrent and Deep Neural Networks) to detect user engagement decrease in realtime. The results of this work highlight, in particular, the relevance of taking into account the temporal dynamics of a user's behavior. Allowing 1 to 2 seconds as buffer delay improves the performance of taking a decision on user engagement.


title: Keck/NIRC2 $L$'-Band Imaging of Jovian-Mass Accreting Protoplanets around PDS 70

abstract: We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.


title: Slowly rotating black holes in Einsteinian cubic gravity

abstract: We construct slowly rotating black-hole solutions of Einsteinian cubic gravity (ECG) in four dimensions with flat and AdS asymptotes. At leading order in the rotation parameter, the only modification with respect to the static case is the appearance of a non-vanishing $g_{t\phi}$ component. Similarly to the static case, the order of the equation determining such component can be reduced twice, giving rise to a second-order differential equation which can be easily solved numerically as a function of the ECG coupling. We study how various physical properties of the solutions are modified with respect to the Einstein gravity case, including its angular velocity, photon sphere, photon rings, shadow, and innermost stable circular orbits (in the case of timelike geodesics).


title: Distributed Resources for the Earth System Grid Advanced Management (DREAM)

abstract: The DREAM project was funded more than 3 years ago to design and implement a next-generation ESGF (Earth System Grid Federation [1]) architecture which would be suitable for managing and accessing data and services resources on a distributed and scalable environment. In particular, the project intended to focus on the computing and visualization capabilities of the stack, which at the time were rather primitive. At the beginning, the team had the general notion that a better ESGF architecture could be built by modularizing each component, and redefining its interaction with other components by defining and exposing a well defined API. Although this was still the high level principle that guided the work, the DREAM project was able to accomplish its goals by leveraging new practices in IT that started just about 3 or 4 years ago: the advent of containerization technologies (specifically, Docker), the development of frameworks to manage containers at scale (Docker Swarm and Kubernetes), and their application to the commercial Cloud. Thanks to these new technologies, DREAM was able to improve the ESGF architecture (including its computing and visualization services) to a level of deployability and scalability beyond the original expectations.


title: Why do People Share Misinformation during the COVID-19 Pandemic?

abstract: The World Health Organization have emphasised that misinformation - spreading rapidly through social media - poses a serious threat to the COVID-19 response. Drawing from theories of health perception and cognitive load, we develop and test a research model hypothesizing why people share unverified COVID-19 information through social media. Our findings suggest a person's trust in online information and perceived information overload are strong predictors of unverified information sharing. Furthermore, these factors, along with a person's perceived COVID-19 severity and vulnerability influence cyberchondria. Females were significantly more likely to suffer from cyberchondria, however, males were more likely to share news without fact checking their source. Our findings suggest that to mitigate the spread of COVID-19 misinformation and cyberchondria, measures should be taken to enhance a healthy skepticism of health news while simultaneously guarding against information overload.


title: An Automated Pipeline for Character and Relationship Extraction from Readers' Literary Book Reviews on Goodreads.com

abstract: Reader reviews of literary fiction on social media, especially those in persistent, dedicated forums, create and are in turn driven by underlying narrative frameworks. In their comments about a novel, readers generally include only a subset of characters and their relationships, thus offering a limited perspective on that work. Yet in aggregate, these reviews capture an underlying narrative framework comprised of different actants (people, places, things), their roles, and interactions that we label the "consensus narrative framework". We represent this framework in the form of an actant-relationship story graph. Extracting this graph is a challenging computational problem, which we pose as a latent graphical model estimation problem. Posts and reviews are viewed as samples of sub graphs/networks of the hidden narrative framework. Inspired by the qualitative narrative theory of Greimas, we formulate a graphical generative Machine Learning (ML) model where nodes represent actants, and multi-edges and self-loops among nodes capture context-specific relationships. We develop a pipeline of interlocking automated methods to extract key actants and their relationships, and apply it to thousands of reviews and comments posted on Goodreads.com. We manually derive the ground truth narrative framework from SparkNotes, and then use word embedding tools to compare relationships in ground truth networks with our extracted networks. We find that our automated methodology generates highly accurate consensus narrative frameworks: for our four target novels, with approximately 2900 reviews per novel, we report average coverage/recall of important relationships of > 80% and an average edge detection rate of >89\%. These extracted narrative frameworks can generate insight into how people (or classes of people) read and how they recount what they have read to others.


title: Integer Quantization for Deep Learning Inference: Principles and Empirical Evaluation

abstract: Quantization techniques can reduce the size of Deep Neural Networks and improve inference latency and throughput by taking advantage of high throughput integer instructions. In this paper we review the mathematical aspects of quantization parameters and evaluate their choices on a wide range of neural network models for different application domains, including vision, speech, and language. We focus on quantization techniques that are amenable to acceleration by processors with high-throughput integer math pipelines. We also present a workflow for 8-bit quantization that is able to maintain accuracy within 1% of the floating-point baseline on all networks studied, including models that are more difficult to quantize, such as MobileNets and BERT-large.


title: Echo Chambers on Social Media: A comparative analysis

abstract: Recent studies have shown that online users tend to select information adhering to their system of beliefs, ignore information that does not, and join groups - i.e., echo chambers - around a shared narrative. Although a quantitative methodology for their identification is still missing, the phenomenon of echo chambers is widely debated both at scientific and political level. To shed light on this issue, we introduce an operational definition of echo chambers and perform a massive comparative analysis on more than 1B pieces of contents produced by 1M users on four social media platforms: Facebook, Twitter, Reddit, and Gab. We infer the leaning of users about controversial topics - ranging from vaccines to abortion - and reconstruct their interaction networks by analyzing different features, such as shared links domain, followed pages, follower relationship and commented posts. Our method quantifies the existence of echo-chambers along two main dimensions: homophily in the interaction networks and bias in the information diffusion toward likely-minded peers. We find peculiar differences across social media. Indeed, while Facebook and Twitter present clear-cut echo chambers in all the observed dataset, Reddit and Gab do not. Finally, we test the role of the social media platform on news consumption by comparing Reddit and Facebook. Again, we find support for the hypothesis that platforms implementing news feed algorithms like Facebook may elicit the emergence of echo-chambers.


title: Frequency Control Analysis based on Unit Commitment Schemes with High Wind Power Integration: a Spanish Isolated Power System Case Study

abstract: The relevant integration of wind power into the grid has involved a remarkable impact on power system operation, mainly in terms of security and reliability due to the inherent loss of the rotational inertia as a consequence of such new generation units decoupled from the grid. In these weak scenarios, the contribution of wind turbines to frequency control is considered as a suitable solution to improve system stability. With regard to frequency response analysis and grid stability, most contributions introduce wind control discuss generation tripping for isolated power systems under arbitrary power imbalance conditions. Frequency response is then analyzed for hypothetical imbalances usually ranged between 5% and 20%, and assuming averaged energy schedule scenarios. In this paper, a more realistic framework is proposed to evaluate frequency deviations by including high wind power integration. With this aim, unit commitment schemes and frequency load shedding are considered in this work for frequency response analysis under high wind power penetration. The Gran Canaria Island's isolated power system (Spain) is used for evaluation purposes. Results provide a variety of influences from wind frequency control depending not only on the wind power integration, but also the generation units under operation, the rotational inertia reductions as well as the available reserves from each resource, aspects that have not been addressed previously in the specific literature to evaluate frequency excursions under high wind power integration.


title: Cliques and constructors in "Hats'"game

abstract: The following general variant of deterministic Hats game is analyzed. Several sages wearing colored hats occupy the vertices of a graph, the $k$-th sage can have hats of one of $h(k)$ colors. Each sage tries to guess the color of his own hat merely on the basis of observing the hats of his neighbors without exchanging any information. A predetermined guessing strategy is winning if it guarantees at least one correct individual guess for every assignment of colors. For complete graphs and for cycles we solve the problem of describing functions $h(k)$ for which the sages win. We demonstrate here winning strategies for the sages on complete graphs, and analyze the Hats game on almost complete graphs. We develop "theory of constructors", that is a collection of theorems demonstrating how one can construct new graphs for which the sages win. We define also new game "Check by rook" which is equivalent to Hats game on 4-cycle and give complete analysis of this game.


title: H$\alpha$ Distances to the Leading Arm of the Magellanic Stream

abstract: The Leading Arm is a tidal feature that is in front of the Magellanic Clouds on their orbit through the Galaxy's halo. Many physical properties of the Leading Arm, such as its mass and size, are poorly constrained because it has few distance measurements. While H$\alpha$ measurements have been used to estimate the distances to halo clouds, many studies have been unsuccessful in detecting H$\alpha$ from the Leading Arm. In this study, we explore a group of H I clouds which lie $75^{\circ} - 90^{\circ}$ from the Magellanic Clouds. Through ultraviolet and 21-cm radio spectroscopy, this region, dubbed the Leading Arm Extension, was found to have chemical and kinematic similarities to the Leading Arm. Using the Wisconsin H$\alpha$ Mapper, we detect H$\alpha$ emission in four out of seven of our targets. Assuming that this region is predominantly photoionized, we use a radiation model that incorporates the contributions of the Galaxy, Magellanic Clouds, and the extragalactic background at $\rm z = 0$ to derive a heliocentric distance of $d_{\odot}\ge13.4~kpc$. We also use this model to rederive H$\alpha$ distances of $d_{\odot} \geq 5.0$ kpc and $d_{\odot} \geq 22.9~kpc$ to two clouds in the literature that might also be associated with the Leading Arm. Using these new measurements, and others in the literature, we provide a general trend of the variation of Leading Arm heliocentric distance as a function of Magellanic Stream longitude, and explore its implications for the origin and closest point of approach of the Leading Arm.


title: Data Processing for Optimizing Naturalness of Vietnamese Text-to-speech System

abstract: Abstract End-to-end text-to-speech (TTS) systems has proved its great success in the presence of a large amount of high-quality training data recorded in anechoic room with high-quality microphone. Another approach is to use available source of found data like radio broadcast news. We aim to optimize the naturalness of TTS system on the found data using a novel data processing method. The data processing method includes 1) utterance selection and 2) prosodic punctuation insertion to prepare training data which can optimize the naturalness of TTS systems. We showed that using the processing data method, an end-to-end TTS achieved a mean opinion score (MOS) of 4.1 compared to 4.3 of natural speech. We showed that the punctuation insertion contributed the most to the result. To facilitate the research and development of TTS systems, we distributed the processed data of one speaker at https://forms.gle/6Hk5YkqgDxAaC2BU6.


title: Flow-based Algorithms for Improving Clusters: A Unifying Framework, Software, and Performance

abstract: Clustering points in a vector space or nodes in a graph is a ubiquitous primitive in statistical data analysis, and it is commonly used for exploratory data analysis. In practice, it is often of interest to "refine" or "improve" a given cluster that has been obtained by some other method. In this survey, we focus on principled algorithms for this cluster improvement problem. Many such cluster improvement algorithms are flow-based methods, by which we mean that operationally they require the solution of a sequence of maximum flow problems on a (typically implicitly) modified data graph. These cluster improvement algorithms are powerful, both in theory and in practice, but they have not been widely adopted for problems such as community detection, local graph clustering, semi-supervised learning, etc. Possible reasons for this are: the steep learning curve for these algorithms; the lack of efficient and easy to use software; and the lack of detailed numerical experiments on real-world data that demonstrate their usefulness. Our objective here is to address these issues. To do so, we guide the reader through the whole process of understanding how to implement and apply these powerful algorithms. We present a unifying fractional programming optimization framework that permits us to distill, in a simple way, the crucial components of all these algorithms. It also makes apparent similarities and differences between related methods. Viewing these cluster improvement algorithms via a fractional programming framework suggests directions for future algorithm development. Finally, we develop efficient implementations of these algorithms in our LocalGraphClustering Python package, and we perform extensive numerical experiments to demonstrate the performance of these methods on social networks and image-based data graphs.


title: Solar slow magneto-acoustic-gravity waves: an erratum correction and a revisited scenario

abstract: Slow waves are commonly observed on the entire solar atmosphere. Assuming a thin flux tube approximation, the cut-off periods of slow-mode magneto-acoustic-gravity waves that travel from the photosphere to the corona were obtained in Costa et al. (2018). In that paper, however, a typo in the specific heat coefficient at constant pressure $c_{\mathrm{p}}$ value led to an inconsistency in the cut-off calculation, which is only significant at the transition region. Due to the abrupt temperature change in the region, a change of the mean atomic weight (by a factor of approximately two) also occurs, but is often overlooked in analytical models for simplicity purposes. In this paper, we revisit the calculation of the cut-off periods of magneto-acoustic-gravity waves in Costa et al. (2018) by considering an atmosphere in hydrostatic equilibrium with a temperature profile, with the inclusion of the variation of the mean atomic weight and the correction of the inconsistency aforementioned. In addition, we show that the cut-off periods obtained analytically are consistent with the corresponding periods measured in observations of a particular active region.


title: Deep variational network for rapid 4D flow MRI reconstruction

abstract: Phase-contrast magnetic resonance imaging (MRI) provides time-resolved quantification of blood flow dynamics that can aid clinical diagnosis. Long in vivo scan times due to repeated three-dimensional (3D) volume sampling over cardiac phases and breathing cycles necessitate accelerated imaging techniques that leverage data correlations. Standard compressed sensing reconstruction methods require tuning of hyperparameters and are computationally expensive, which diminishes the potential reduction of examination times. We propose an efficient model-based deep neural reconstruction network and evaluate its performance on clinical aortic flow data. The network is shown to reconstruct undersampled 4D flow MRI data in under a minute on standard consumer hardware. Remarkably, the relatively low amounts of tunable parameters allowed the network to be trained on images from 11 reference scans while generalizing well to retrospective and prospective undersampled data for various acceleration factors and anatomies.


title: Reduced Tensor Product on the Drinfeld Center

abstract: The annulus comes with a "stacking" operation which glues two annuli into one. This provides a tensor product structure on the category of boundary values $Z_{\text{CY}}(\text{Ann})$ associated to the annulus in an extended Crane-Yetter TQFT. It is known that $Z_{\text{CY}}(\text{Ann}) \simeq \mathcal{Z}(\mathcal{A})$, the Drinfeld center of the premodular category $\mathcal{A}$ from which $Z_{\text{CY}}$ is constructed. We give an explicit formula for the tensor product on $\mathcal{Z}(\mathcal{A})$ that corresponds to the stacking operation on $Z_{\text{CY}}(\text{Ann})$.


title: A Critical Overview of Privacy-Preserving Approaches for Collaborative Forecasting

abstract: Cooperation between different data owners may lead to an improvement in forecast quality - for instance by benefiting from spatial-temporal dependencies in geographically distributed time series. Due to business competitive factors and personal data protection questions, said data owners might be unwilling to share their data, which increases the interest in collaborative privacy-preserving forecasting. This paper analyses the state-of-the-art and unveils several shortcomings of existing methods in guaranteeing data privacy when employing Vector Autoregressive (VAR) models. The paper also provides mathematical proofs and numerical analysis to evaluate existing privacy-preserving methods, dividing them into three groups: data transformation, secure multi-party computations, and decomposition methods. The analysis shows that state-of-the-art techniques have limitations in preserving data privacy, such as a trade-off between privacy and forecasting accuracy, while the original data in iterative model fitting processes, in which intermediate results are shared, can be inferred after some iterations.


title: Fixed-Target Runtime Analysis

abstract: Runtime analysis aims at contributing to our understanding of evolutionary algorithms through mathematical analyses of their runtimes. In the context of discrete optimization problems, runtime analysis classically studies the time needed to find an optimal solution. However, both from a practical and from a theoretical viewpoint, more fine-grained performance measures are needed to gain a more detailed understanding of the main working principles and their resulting performance implications. Two complementary approaches have been suggested: fixed-budget analyses and fixed-target analyses. In this work, we conduct an in-depth study on the advantages and the limitations of fixed-target analyses. We show that, different from fixed-budget analyses, many classical methods from the runtime analysis of discrete evolutionary algorithms yield fixed-target results without greater effort. We use this to conduct a number of new fixed-target analyses. However, we also point out examples where an extension of existing runtime results to fixed-target results is highly non-trivial.


title: 100,000 frames-per-second compressive imaging with a conventional rolling-shutter camera by random point-spread-function engineering

abstract: We demonstrate an approach that allows taking videos at very high-speeds of over 100,000 frames per second (fps) by exploiting the fast sampling rate of the standard rolling-shutter readout mechanism, common to most conventional sensors, and a compressive-sampling acquisition scheme. Our approach is directly applied to a conventional imaging system by the simple addition of a diffuser to the pupil plane, randomly encoding the entire field-of-view to each camera row, while maintaining diffraction-limited resolution. A short video is reconstructed from a single camera frame via a compressed-sensing reconstruction algorithm, exploiting inherent sparsity of the imaged scene.


title: Sampling and Inference of Networked Dynamics using Log-Koopman Nonlinear Graph Fourier Transform

abstract: Networked nonlinear dynamics underpin the complex functionality of many engineering, social, biological, and ecological systems. Monitoring the networked dynamics via the minimum subset of nodes is essential for a variety of scientific and operational purposes. When there is a lack of a explicit model and networked signal space, traditional evolution analysis and non-convex methods are insufficient. An important data-driven state-of-the-art method use the Koopman operator to generate a linear evolution model for a vector-valued observable of original state-space. As a result, one can derive a sampling strategy via the linear evolution property of observable. However, current polynomial Koopman operators result in a large sampling space due to: (i) the large size of polynomial based observables ($O(N^2)$, $N$ number of nodes in network), and (ii) not factoring in the nonlinear dependency between observables. In this work, to achieve linear scaling ($O(N)$) and a small set of sampling nodes, we propose to combine a novel Log-Koopman operator and nonlinear Graph Fourier Transform (NL-GFT) scheme. First, the Log-Koopman operator is able to reduce the size of observables by transforming multiplicative poly-observable to logarithm summation. Second, a nonlinear GFT concept and sampling theory are provided to exploit the nonlinear dependence of observables for Koopman linearized evolution analysis. Combined, the sampling and reconstruction algorithms are designed and demonstrated on two established application areas. The results demonstrate that the proposed Log-Koopman NL-GFT scheme can (i) linearize unknown nonlinear dynamics using $O(N)$ observables, and (ii) achieve lower number of sampling nodes, compared with the state-of-the art polynomial Koopman linear evolution analysis.


title: Global existence of weak solutions to a diffuse interface model for magnetic fluids

abstract: This article is devoted to the derivation and analysis of a system of partial differential equations modeling a diffuse interface flow of two Newtonian incompressible magnetic fluids. The system consists of the incompressible Navier-Stokes equations coupled with an evolutionary equation for the magnetization vector and the Cahn-Hilliard equations. We show global in time existence of weak solutions to the system using the time discretization method.


title: A geometric characterization of VES and Kadiyala-type production functions

abstract: The basic concepts of the differential geometry are shortly reviewed and applied to the study of VES production function in the spirit of the works of V\^ilcu and collaborators. A similar characterization is given for a more general production function, namely the Kadiyala production function, in the case of developable surfaces.


title: Global well-posedness for the defocusing, cubic nonlinear Schr{\"o}dinger equation with initial data in a critical space

abstract: In this note we prove global well-posedness for the defocusing, cubic nonlinear Schr{\"o}dinger equation with initial data lying in a critical Sobolev space.


title: Vilamb: Low Overhead Asynchronous Redundancy for Direct Access NVM

abstract: Vilamb provides efficient asynchronous systemredundancy for direct access (DAX) non-volatile memory (NVM) storage. Production storage deployments often use system-redundancy in form of page checksums and cross-page parity. State-of-the-art solutions for maintaining system-redundancy for DAX NVM either incur a high performance overhead or require specialized hardware. The Vilamb user-space library maintains system-redundancy with low overhead by delaying and amortizing the system-redundancy updates over multiple data writes. As a result, Vilamb provides 3--5x the throughput of the state-of-the-art software solution at high operation rates. For applications that need system-redundancy with high performance, and can tolerate some delaying of data redundancy, Vilamb provides a tunable knob between performance and quicker redundancy. Even with the delayed coverage, Vilamb increases the mean time to data loss due to firmware-induced corruptions by up to two orders of magnitude in comparison to maintaining no system-redundancy.


title: Symplectic duality and implosions

abstract: We discuss symplectic and hyperk\"ahler implosion and present candidates for the symplectic duals of the universal hyperk\"ahler implosion for various groups.


title: Characterizing consensus in the Heard-Of model

abstract: The Heard-Of model is a simple and relatively expressive model of distributed computation. Because of this, it has gained a considerable attention of the verification community. We give a characterization of all algorithms solving consensus in a fragment of this model. The fragment is big enough to cover many prominent consensus algorithms. The characterization is purely syntactic: it is expressed in terms of some conditions on the text of the algorithm. One of the recent methods of verification of distributed algorithms is to abstract an algorithm to the Heard-Of model and then to verify the abstract algorithm using semi-automatic procedures. Our results allow, in some cases, to avoid the second step in this methodology.


title: Ground-state properties of the $K-\Gamma$ model on a honeycomb lattice

abstract: We investigate the ground-state properies of the $K-\Gamma$ model on a honeycomb lattice using series expansions and numerical exact diagonalizations, where the model includes Kitaev ($K$) and symmetric off-diagonal ($\Gamma$) interactions. Starting from the weakly interacting dimers on the specific bond, we strengthen the interdimer interactions to the isotropically interacting system. We show that depending on $\Gamma$ and $K$, the dimer state survives up to the isotropically interacting system, where the phase transition occurs, or obeys a phase transition to a magnetically ordered state at an anisotropic interaction. The results are summarized in the phase diagram. We also show that the Kekul\'{e} dimerized state is unstable in the isotropic $K-\Gamma$ model.


title: Fast Multivariate Probit Estimation via a Two-Stage Composite Likelihood

abstract: The multivariate probit is popular for modeling correlated binary data, with an attractive balance of flexibility and simplicity. However, considerable challenges remain in computation and in devising a clear statistical framework. Interest in the multivariate probit has increased in recent years. Current applications include genomics and precision medicine, where simultaneous modeling of multiple traits may be of interest, and computational efficiency is an important consideration. We propose a fast method for multivariate probit estimation via a two-stage composite likelihood. We explore computational and statistical efficiency, and note that the approach sets the stage for extensions beyond the purely binary setting.


title: Well-posedness for the Majda-Biello system on the half line

abstract: We study the initial-boundary value problem for the Majda-Biello system posed on the right half line. We prove local well-posedness on the half line, matching the local theory on the real line established by Oh (2008). The approach combines the Laplace transform method of Bona-Sun-Zhang with adapted estimates from the work of Colliander and Kenig on the KdV half line initial-boundary value problem.


title: A New Community Definition For MultiLayer Networks And A Novel Approach For Its Efficient Computation

abstract: As the use of MultiLayer Networks (or MLNs) for modeling and analysis is gaining popularity, it is becoming increasingly important to propose a community definition that encompasses the multiple features represented by MLNs and develop algorithms for efficiently computing communities on MLNs. Currently, communities for MLNs, are based on aggregating the networks into single graphs using different techniques (type independent, projection-based, etc.) and applying single graph community detection algorithms, such as Louvain and Infomap on these graphs. This process results in different types of information loss (semantics and structure). To the best of our knowledge, in this paper we propose, for the first time, a definition of community for heterogeneous MLNs (or HeMLNs) which preserves semantics as well as the structure. Additionally, our basic definition can be extended to appropriately match the analysis objectives as needed. In this paper, we present a structure and semantics preserving community definition for HeMLNs that is compatible with and is an extension of the traditional definition for single graphs. We also present a framework for its efficient computation using a newly proposed decoupling approach. First, we define a k-community for connected k layers of a HeMLN. Then we propose a family of algorithms for its computation using the concept of bipartite graph pairings. Further, for a broader analysis, we introduce several pairing algorithms and weight metrics for composing binary HeMLN communities using participating community characteristics. Essentially, this results in an extensible family of community computations. We provide extensive experimental results for showcasing the efficiency and analysis flexibility of the proposed computation using popular IMDb and DBLP data sets.


title: An exposition of the equivalence of Heegaard Floer homology and embedded contact homology

abstract: This article is a survey on the authors' proof of the isomorphism between Heegaard Floer homology and embedded contact homology appeared in This article is a survey on the authors' proof of the isomorphism between Heegaard Floer homology and embedded contact homology appeared in arXiv:1208.1074, arXiv:1208.1077 and arXiv:1208.1526.


title: Inference by Stochastic Optimization: A Free-Lunch Bootstrap

abstract: Assessing sampling uncertainty in extremum estimation can be challenging when the asymptotic variance is not analytically tractable. Bootstrap inference offers a feasible solution but can be computationally costly especially when the model is complex. This paper uses iterates of a specially designed stochastic optimization algorithm as draws from which both point estimates and bootstrap standard errors can be computed in a single run. The draws are generated by the gradient and Hessian computed from batches of data that are resampled at each iteration. We show that these draws yield consistent estimates and asymptotically valid frequentist inference for a large class of regular problems. The algorithm provides accurate standard errors in simulation examples and empirical applications at low computational costs. The draws from the algorithm also provide a convenient way to detect data irregularities.


title: Two-Level Lattice Neural Network Architectures for Control of Nonlinear Systems

abstract: In this paper, we consider the problem of automatically designing a Rectified Linear Unit (ReLU) Neural Network (NN) architecture (number of layers and number of neurons per layer) with the guarantee that it is sufficiently parametrized to control a nonlinear system. Whereas current state-of-the-art techniques are based on hand-picked architectures or heuristic based search to find such NN architectures, our approach exploits the given model of the system to design an architecture; as a result, we provide a guarantee that the resulting NN architecture is sufficient to implement a controller that satisfies an achievable specification. Our approach exploits two basic ideas. First, assuming that the system can be controlled by an unknown Lipschitz-continuous state-feedback controller with some Lipschitz constant upper-bounded by $K_\text{cont}$, we bound the number of affine functions needed to construct a Continuous Piecewise Affine (CPWA) function that can approximate the unknown Lipschitz-continuous controller. Second, we utilize the authors' recent results on a novel NN architecture named as the Two-Level Lattice (TLL) NN architecture, which was shown to be capable of implementing any CPWA function just from the knowledge of the number of affine functions that compromises this CPWA function.


title: Self-Supervised Feature Extraction for 3D Axon Segmentation

abstract: Existing learning-based methods to automatically trace axons in 3D brain imagery often rely on manually annotated segmentation labels. Labeling is a labor-intensive process and is not scalable to whole-brain analysis, which is needed for improved understanding of brain function. We propose a self-supervised auxiliary task that utilizes the tube-like structure of axons to build a feature extractor from unlabeled data. The proposed auxiliary task constrains a 3D convolutional neural network (CNN) to predict the order of permuted slices in an input 3D volume. By solving this task, the 3D CNN is able to learn features without ground-truth labels that are useful for downstream segmentation with the 3D U-Net model. To the best of our knowledge, our model is the first to perform automated segmentation of axons imaged at subcellular resolution with the SHIELD technique. We demonstrate improved segmentation performance over the 3D U-Net model on both the SHIELD PVGPe dataset and the BigNeuron Project, single neuron Janelia dataset.


title: Optimization of Chaos-based Coded Modulations for Compensation of Amplifier Nonlinearities

abstract: In this work we expand the possibilities of a known class of chaos-based coded modulation (CCM) systems to address the problem of high-power amplifier (HPA) nonlinearity. The hypothesis is that the CCM nonlinear nature can be exploited to counteract its effects on the system performance, avoiding the need of a predistorter. We propose an optimization method for the design of the CCM to prove this hypothesis. The results show that, for a given back-off level, a nonlinear mapping function for the chaos-based encoder and decoder can be designed to compensate the HPA effect on the error rate.


title: Non-Precipitating Shallow Cumulus Clouds: Theory and Direct Numerical Simulation

abstract: We model cumulus clouds as transient diabatic plumes, and present a single-fluid formulation for the study of the dynamics of developing shallow cumulus clouds. The fluid, air, carries water vapour and liquid water as scalars that are advected with the flow. We employ the Boussinesq approximation and a simplified form of the Clausius-Clapeyron law for phase changes. We use this formulation to study the ambient conditions required for the formation of cumulus clouds. We show that the temperature lapse rate in the ambient and the relative humidity of the ambient are crucial in deciding when cumulus clouds form. In the phase-space of these two parameters, we determine the boundary that separates the regions where cumulus clouds can form and where they cannot.


title: Intelligent Querying for Target Tracking in Camera Networks using Deep Q-Learning with n-Step Bootstrapping

abstract: Surveillance camera networks are a useful infrastructure for various visual analytics applications, where high-level inferences and predictions could be made based on target tracking across the network. Most multi-camera tracking works focus on target re-identification and trajectory association problems to track the target. However, since camera networks can generate enormous amount of video data, inefficient schemes for making re-identification or trajectory association queries can incur prohibitively large computational requirements. In this paper, we address the problem of intelligent scheduling of re-identification queries in a multi-camera tracking setting. To this end, we formulate the target tracking problem in a camera network as an MDP and learn a reinforcement learning based policy that selects a camera for making a re-identification query. The proposed approach to camera selection does not assume the knowledge of the camera network topology but the resulting policy implicitly learns it. We have also shown that such a policy can be learnt directly from data. Using the NLPR MCT and the Duke MTMC multi-camera multi-target tracking benchmarks, we empirically show that the proposed approach substantially reduces the number of frames queried.


title: CHY Theory for Several Fields

abstract: The Cachazo-He-Yuan (CHY) formula was originally proposed to describe on-shell scattering of particles from a single massless field. We present a method to modify it to include several interacting scalar fields, all possessing different masses and possibly off-shell momenta. The method is applied to Yukawa interactions between a number of scalar nucleons and pions, and to the $\f_1\f_2\f_3$ coupling of three different scalar fields. Composite models constructed from existing theories can be used to broaden the scope of the method. The modification is applied to describe Compton scattering from a massive particle, and to photon bremsstrahlung. It is also employed to generalize the disk function $Z$ and the sphere function $J$.


title: Radiation reaction friction: Resistive material medium

abstract: We explore a novel method of describing the radiation friction of particles traveling through a mechanically resistive medium. We introduce a particle motion induced matter warping along the path in a manner assuring that charged particle dynamics occurs subject to radiative energy loss described by the Larmor formula. We compare our description with the Landau-Lifshitz-like model for the radiation friction and show that the established model exhibits non-physical behavior. Our approach predicts in the presence of large mechanical friction an upper limit on radiative energy loss being equal to the energy loss due to the mechanical medium resistance. We demonstrate that mechanical friction due to strong interactions, for example of quarks in quark-gluon plasma, can induce significant soft photon radiation.


title: Twisted conjugacy in linear algebraic groups

abstract: Let $k$ be an algebraically closed field, $G$ a linear algebraic group over $k$ and $\varphi\in Aut(G)$, the group of all algebraic group automorphisms of $G$. Two elements $x, y$ of $G$ are said to be $\varphi$-twisted conjugate if $y=gx\varphi(g)^{-1}$ for some $g\in G$. In this paper we prove that for a connected non-solvable linear algebraic group $G$ over $k$, the number of its $\varphi$-twisted conjugacy classes is infinite for every $\varphi\in Aut(G)$.


title: Searching for globular cluster chemical anomalies on the main sequence of a young massive cluster

abstract: The spectroscopic and photometric signals of the star-to-star abundance variations found in globular clusters seem to be correlated with global parameters like the cluster's metallicity, mass and age. Understanding this behaviour could bring us closer to the origin of these intriguing abundance spreads. In this work we use deep HST photometry to look for evidence of abundance variations in the main sequence of a young massive cluster NGC 419 ($\sim10^5$ M$_{\odot}$, $\sim1.4$ Gyr). Unlike previous studies, here we focus on stars in the same mass range found in old globulars ($\sim0.75-1$ M$_{\odot}$), where light elements variations are detected. We find no evidence for N abundance variations among these stars in the $Un-B$ and $U-B$ CMD of NGC 419. This is at odds with the N-variations found in old globulars like 47 Tuc, NGC 6352 and NGC 6637 with similar metallicity to NGC 419. Although the signature of the abundance variations characteristic of old globulars appears to be significantly smaller or absent in this young cluster, we cannot conclude if this effect is mainly driven by its age or its mass.


title: Grassmannian stochastic analysis and the stochastic quantization of Euclidean Fermions

abstract: We introduce a stochastic analysis of Grassmann random variables suitable for the stochastic quantization of Euclidean fermionic quantum field theories. Analysis on Grassmann algebras is developed here from the point of view of quantum probability: a Grassmann random variable is an homomorphism of an abstract Grassmann algebra into a quantum probability space, i.e. a $C^{\ast}$-algebra endowed with a suitable state. We define the notion of Gaussian processes, Brownian motion and stochastic (partial) differential equations taking values in Grassmann algebras. We use them to study the long time behavior of finite and infinite dimensional Langevin Grassmann stochastic differential equations driven by Gaussian space-time white noise and to describe their invariant measures. As an application we give a proof of the stochastic quantization and of the removal of the space cut-off for the Euclidean Yukawa model.


title: Stability of steady state solutions in a predator-prey model with density-dependent diffusion

abstract: In this note we present a study of the solutions associated to a particular spatial extension of the Rosenzweig-MacArthur model for predator and prey. The analysis presented here shows that positive steady state solutions emerge via a transcritical bifurcation mechanism, in accordance with the insight obtained from previous numerical and analytical results. In the model under discussion, prey is assumed to move avoiding crowds via a density-dependent diffusion and also incorporates the existence of a refuge zone, where predators cannot consume prey. Saturation in prey consumption is also included through a Holling type II functional response.


title: Impairment of insulin-stimulated glucose utilization is associated with burn-induced insulin resistance in mouse muscle by hyperinsulinemic-isoglycemic clamp

abstract: Burn-induced insulin resistance is associated with increased morbidity and mortality; however, the impact of burn injury on tissue-specific insulin sensitivity and its molecular mechanisms with consideration of insulin state remains unknown in rodent models. This study was designed to characterize a burn mouse model with tissue-specific insulin resistance under insulin clamp conditions. C57BL6/J mice were subjected to 30% full-thickness burn injury and underwent the combination of hyperinsulinemic isoglycemicclamp (HIC) and positron emission tomography (PET). Hepatic glucose production (HGP) and peripheral glucose disappearance rate (Rd) were measured at different time points up to 7 days post injury. Burned mice showed a significant fasting hypoglycemia and hypoinsulinemia (P < 0.01) on post-burn day (PBD) 3 and 7 along with significantly higher energy expenditure (P < 0.01). HICon PBD 3 demonstrated that burn injury induced systemic insulin resistance, resulting from a significant decrease in insulin-stimulated Rd (33.0 +/- 10.2 vs 68.3 +/- 5.9 mg/kg/min; P < 0.05). In contrast, HGP of burned and sham mice was comparable both in the basal and clamp period. PET on PBD 3 showed a lower insulin-stimulated 18F-labeled 2-fluoro-2-deoxy-D-glucose uptake in the quadriceps of burned mice compared with sham-burned mice. Gastrocnemius muscle harvested from burned mice on PBD 3 showed decreased insulin-stimulated tyrosine phosphorylation of insulin receptor substrate-1 to 34.7% of that in sham-burn mice by immunoblotting analysis (P < 0.05). These findings suggest that impaired insulin-stimulated Rd in skeletal muscle, not elevated HGP, plays a role in the development of burn-induced insulin resistance in a mouse model.


title: Mechanism Design for Online Resource Allocation: A Unified Approach

abstract: This paper concerns the mechanism design for online resource allocation in a strategic setting. In this setting, a single supplier allocates capacity-limited resources to requests that arrive in a sequential and arbitrary manner. Each request is associated with an agent who may act selfishly to misreport the requirement and valuation of her request. The supplier charges payment from agents whose requests are satisfied, but incurs a load-dependent supply cost. The goal is to design an incentive compatible online mechanism, which determines not only the resource allocation of each request, but also the payment of each agent, so as to (approximately) maximize the social welfare (i.e., aggregate valuations minus supply cost). We study this problem under the framework of competitive analysis. The major contribution of this paper is the development of a unified approach that achieves the best-possible competitive ratios for setups with different supply costs. Specifically, we show that when there is no supply cost or the supply cost function is linear, our model is essentially a standard 0-1 knapsack problem, for which our approach achieves logarithmic competitive ratios that match the state-of-the-art (which is optimal). For the more challenging setup when the supply cost is strictly-convex, we provide online mechanisms, for the first time, that lead to the optimal competitive ratios as well. To the best of our knowledge, this is the first approach that unifies the characterization of optimal competitive ratios in online resource allocation for different setups including zero, linear and strictly-convex supply costs.


title: Symmetry Adapted Gram Spectrahedra

abstract: This paper explores the geometric structure of the spectrahedral cone, called the symmetry adapted PSD cone, and the symmetry adapted Gram spectrahedron of a symmetric polynomial. In particular, we determine the dimension of the symmetry adapted PSD cone, describe its extreme rays, and discuss the structure of its matrix representations. We also consider the symmetry adapted Gram spectrahedra for specific families of symmetric polynomials including binary symmetric polynomials, quadratics, and ternary quartics and sextics which give us further insight into these symmetric SOS polynomials. Finally, we discuss applications of the theory of sums of squares and symmetric polynomials which arise from symmetric function inequalities.


title: Kibble-Zurek Dynamics in a Trapped Ultracold Bose Gas

abstract: The dynamical evolution of an inhomogeneous ultracold atomic gas quenched at different controllable rates through the Bose-Einstein condensation phase transition is studied numerically in the premise of a recent experiment in an anisotropic harmonic trap. Our findings based on the stochastic (projected) Gross-Pitaevskii equation are shown to be consistent at early times with the predictions of the homogeneous Kibble-Zurek mechanism. This is demonstrated by collapsing the early dynamical evolution of densities, spectral functions and correlation lengths for different quench rates, based on an appropriate characterization of the distance to criticality felt by the quenched system. The subsequent long-time evolution, beyond the identified dynamical critical region, is also investigated by looking at the behaviour of the density wavefront evolution and the corresponding phase ordering dynamics.


title: How exponentially ill-conditioned are contiguous submatrices of the Fourier matrix?

abstract: We show that the condition number of any cyclically contiguous $p\times q$ submatrix of the $N\times N$ discrete Fourier transform (DFT) matrix is at least $$ \exp \left( \frac{\pi}{2} \left[\min(p,q)- \frac{pq}{N}\right] \right)~, $$ up to algebraic prefactors. That is, fixing any shape parameters $(\alpha,\beta):=(p/N,q/N)\in(0,1)^2$, the growth is $e^{\rho N}$ as $N\to\infty$ with rate $\rho = \frac{\pi}{2}[\min(\alpha,\beta)- \alpha\beta]$. Such Vandermonde system matrices arise in many applications, such as Fourier continuation, super-resolution, and diffraction imaging. Our proof uses the Kaiser-Bessel transform pair (of which we give a self-contained proof), and estimates on sums over distorted sinc functions, to construct a localized trial vector whose DFT is also localized. We warm up with an elementary proof of the above but with half the rate, via a periodized Gaussian trial vector. Using low-rank approximation of the kernel $e^{ixt}$, we also prove another lower bound $(4/e\pi \alpha)^q$, up to algebraic prefactors, which is stronger than the above for small $\alpha, \beta$. When combined, the bounds are within a factor of two of the numerically-measured empirical asymptotic rate, uniformly over $(0,1)^2$, and they become sharp in certain regions. However, the results are not asymptotic: they apply to essentially all $N$, $p$, and $q$, and with all constants explicit.


title: Site-bond percolation solution to preventing the propagation of \textit{Phytophthora} zoospores on plantations

abstract: We propose a strategy based on the site-bond percolation to minimize the propagation of \textit{Phytophthora} zoospores on plantations, consisting in introducing physical barriers between neighboring plants. Two clustering processes are distinguished: i) one of cells with the presence of the pathogen, detected on soil analysis; and ii) that of diseased plants, revealed from a visual inspection of the plantation. The former is well described by the standard site-bond percolation. In the latter, the percolation threshold is fitted by a Tsallis distribution when no barriers are introduced. We provide, for both cases, the formulae for the minimal barrier density to prevent the emergence of the spanning cluster. Though this work is focused on a specific pathogen, the model presented here can also be applied to prevent the spreading of other pathogens that disseminate, by other means, from one plant to the neighboring ones. Finally, the application of this strategy to three types of commercialy important Mexican chili plants is also shown.


title: Flattening the Curve: Insights From Queueing Theory

abstract: The worldwide outbreak of the coronavirus was first identified in 2019 in Wuhan, China. Since then, the disease has spread worldwide. As it currently spreading in the United States, policy makers, public health officials and citizens are racing to understand the impact of this virus on the United States healthcare system. They fear that the rapid influx of patients will overwhelm the healthcare system leading to unnecessary fatalities. Most countries and states in America have introduced mitigation strategies, such as social distancing, to decrease the rate of newly infected people, i.e. flattening the curve.In this paper, we analyze the time evolution of the number of people hospitalized due to the coronavirus using the methods of queueing theory. Given that the rate of new infections varies over time as the pandemic evolves, we model the number of coronavirus patients as a dynamical system based on the theory of infinite server queues with non-stationary Poisson arrival rates. With this model we are able to quantify how flattening the curve affects the peak demand for hospital resources. This allows us to characterize how aggressively society must flatten the curve in order to avoid overwhelming the capacity of healthcare system. We also demonstrate how flattening the curve impacts the elapsed time between the peak rate of hospitalizations and the time of the peak demand for the hospital resources. Finally, we present empirical evidence from China, South Korea, Italy and the United States that supports the insights from the model.


title: Causal network learning with non-invertible functional relationships

abstract: Discovery of causal relationships from observational data is an important problem in many areas. Several recent results have established the identifiability of causal DAGs with non-Gaussian and/or nonlinear structural equation models (SEMs). In this paper, we focus on nonlinear SEMs defined by non-invertible functions, which exist in many data domains, and propose a novel test for non-invertible bivariate causal models. We further develop a method to incorporate this test in structure learning of DAGs that contain both linear and nonlinear causal relations. By extensive numerical comparisons, we show that our algorithms outperform existing DAG learning methods in identifying causal graphical structures. We illustrate the practical application of our method in learning causal networks for combinatorial binding of transcription factors from ChIP-Seq data.