title: Scaling Data Association for Hypothesis-Oriented MHT

abstract: Multi-hypothesis tracking is a flexible and intuitive approach to tracking multiple nearby objects. However, the original formulation of its data association step is widely thought to scale poorly with the number of tracked objects. We propose enhancements including handling undetected objects and false measurements without inflating the size of the problem, early stopping during solution calculation, and providing for sparse or gated input. These changes collectively improve the computational time and space requirements of data association so that hundreds or thousands of hypotheses over hundreds of objects may be considered in real time. A multi-sensor simulation demonstrates that scaling up the hypothesis count can significantly improve performance in some applications.


title: Cloud-cloud collision in the DR 21 cloud as a trigger of massive star formation

abstract: We report on a possible cloud-cloud collision in the DR 21 region, which we found through molecular observations with the Nobeyama 45-m telescope. We mapped an area of 8'x12' around the region with twenty molecular lines including the 12CO(J=1-0) and 13CO(J=1-0) emission lines, and sixteen of them were significantly detected. Based on the 12CO and 13CO data, we found five distinct velocity components in the observed region, and we call molecular gas associated with these components -42, -22, -3, 9, and 17 km/s clouds taking after their typical radial velocities. The -3 km/s cloud is the main filamentary cloud (31,000 Mo) associated with young massive stars such as DR21 and DR21(OH), and the 9 km/s cloud is a smaller cloud (3,400 Mo) which may be an extension of the W75 region in the north. The other clouds are much smaller. We found a clear anticorrelation in the distributions of the -3 and 9 km/s clouds, and detected faint 12CO emission having intermediate velocities bridging the two clouds at their intersection. These facts strongly indicate that the two clouds are colliding against each other. In addition, we found that DR21 and DR21(OH) are located in the periphery of the densest part of the 9 km/s cloud, which is consistent with results of recent numerical simulations of cloud-cloud collisions. We therefore suggest that the -3 and 9 km/s clouds are colliding, and that the collision induced the massive star formation in the DR21 cloud. The interaction of the -3 and 9 km/s clouds was previously suggested by Dickel et al. (1978), and our results strongly support their hypothesis of the interaction.


title: Lasserre Integrality Gaps for Graph Spanners and Related Problems

abstract: There has been significant recent progress on algorithms for approximating graph spanners, i.e., algorithms which approximate the best spanner for a given input graph. Essentially all of these algorithms use the same basic LP relaxation, so a variety of papers have studied the limitations of this approach and proved integrality gaps for this LP in a variety of settings. We extend these results by showing that even the strongest lift-and-project methods cannot help significantly, by proving polynomial integrality gaps even for $n^{\Omega(\epsilon)}$ levels of the Lasserre hierarchy, for both the directed and undirected spanner problems. We also extend these integrality gaps to related problems, notably Directed Steiner Network and Shallow-Light Steiner Network.


title: 4D Seismic History Matching Incorporating Unsupervised Learning

abstract: The work discussed and presented in this paper focuses on the history matching of reservoirs by integrating 4D seismic data into the inversion process using machine learning techniques. A new integrated scheme for the reconstruction of petrophysical properties with a modified Ensemble Smoother with Multiple Data Assimilation (ES-MDA) in a synthetic reservoir is proposed. The permeability field inside the reservoir is parametrised with an unsupervised learning approach, namely K-means with Singular Value Decomposition (K-SVD). This is combined with the Orthogonal Matching Pursuit (OMP) technique which is very typical for sparsity promoting regularisation schemes. Moreover, seismic attributes, in particular, acoustic impedance, are parametrised with the Discrete Cosine Transform (DCT). This novel combination of techniques from machine learning, sparsity regularisation, seismic imaging and history matching aims to address the ill-posedness of the inversion of historical production data efficiently using ES-MDA. In the numerical experiments provided, I demonstrate that these sparse representations of the petrophysical properties and the seismic attributes enables to obtain better production data matches to the true production data and to quantify the propagating waterfront better compared to more traditional methods that do not use comparable parametrisation techniques.


title: A semantic-aided particle filter approach for AUV localization

abstract: This paper presents a novel approach to AUV localization, based on a semantic-aided particle filter. Particle filters have been used successfully for robotics localization since many years. Most of the approaches are however based on geometric measurements and geometric information and simulations. In the past years more and more efforts from research goes towards cognitive robotics and the marine domain is not exception. Moving from signal to symbol becomes therefore paramount for more complex applications. This paper presents a contribution in the well-known area of underwater localization, incorporating semantic information. An extension to the standard particle filter approach is presented, based on semantic information of the environment. A comparison with the geometric approach shows the advantages of a semantic layer to successfully perform self-localization.


title: Learning Open Information Extraction of Implicit Relations from Reading Comprehension Datasets

abstract: The relationship between two entities in a sentence is often implied by word order and common sense, rather than an explicit predicate. For example, it is evident that "Fed chair Powell indicates rate hike" implies (Powell, is a, Fed chair) and (Powell, works for, Fed). These tuples are just as significant as the explicit-predicate tuple (Powell, indicates, rate hike), but have much lower recall under traditional Open Information Extraction (OpenIE) systems. Implicit tuples are our term for this type of extraction where the relation is not present in the input sentence. There is very little OpenIE training data available relative to other NLP tasks and none focused on implicit relations. We develop an open source, parse-based tool for converting large reading comprehension datasets to OpenIE datasets and release a dataset 35x larger than previously available by sentence count. A baseline neural model trained on this data outperforms previous methods on the implicit extraction task.


title: An Informational Approach to Cosmological Parameter Estimation

abstract: We introduce a new approach for cosmological parameter estimation based on the information-theoretical Jensen-Shannon divergence (${\cal D}_{\rm JS}$), calculating it for models in the restricted parameter space $\{H_0, w_0, w_a\}$, where $H_0$ is the value of the Hubble constant today, and $w_0$ and $w_a$ are dark energy parameters, with the other parameters held fixed at their best-fit values from the Planck 2018 data. As an application, we investigate the $H_0$ tension between the Planck temperature power spectrum data (TT) and the local astronomical data by comparing the $\Lambda$CDM model with the $w$CDM and the $w_0w_a$CDM dynamic dark energy models. We find agreement with other works using the standard Bayesian inference for parameter estimation; in addition, we show that while the ${\cal D}_{\rm JS}$ is equally minimized for both values of $H_0$ along the $(w_0,w_a)$ plane, the lines of degeneracy are different for each value of $H_0$. This allows for distinguishing between the two, once the value of either $w_0$ or $w_a$ is known.


title: Adaptively Truncating Backpropagation Through Time to Control Gradient Bias

abstract: Truncated backpropagation through time (TBPTT) is a popular method for learning in recurrent neural networks (RNNs) that saves computation and memory at the cost of bias by truncating backpropagation after a fixed number of lags. In practice, choosing the optimal truncation length is difficult: TBPTT will not converge if the truncation length is too small, or will converge slowly if it is too large. We propose an adaptive TBPTT scheme that converts the problem from choosing a temporal lag to one of choosing a tolerable amount of gradient bias. For many realistic RNNs, the TBPTT gradients decay geometrically in expectation for large lags; under this condition, we can control the bias by varying the truncation length adaptively. For RNNs with smooth activation functions, we prove that this bias controls the convergence rate of SGD with biased gradients for our non-convex loss. Using this theory, we develop a practical method for adaptively estimating the truncation length during training. We evaluate our adaptive TBPTT method on synthetic data and language modeling tasks and find that our adaptive TBPTT ameliorates the computational pitfalls of fixed TBPTT.


title: Decision-Making in a Social Multi-Armed Bandit Task: Behavior, Electrophysiology and Pupillometry

abstract: Understanding, predicting, and learning from other people's actions are fundamental human social-cognitive skills. Little is known about how and when we consider other's actions and outcomes when making our own decisions. We developed a novel task to study social influence in decision-making: the social multi-armed bandit task. This task assesses how people learn policies for optimal choices based on their own outcomes and another player's (observed) outcomes. The majority of participants integrated information gained through observation of their partner similarly as information gained through their own actions. This lead to a suboptimal decision-making strategy. Interestingly, event-related potentials time-locked to stimulus onset qualitatively similar but the amplitudes are attenuated in the solo compared to the dyadic version. This might indicate that arousal and attention after receiving a reward are sustained when a second agent is present but not when playing alone.


title: Automated 3D recovery from very high resolution multi-view satellite images

abstract: This paper presents an automated pipeline for processing multi-view satellite images to 3D digital surface models (DSM). The proposed pipeline performs automated geo-referencing and generates high-quality densely matched point clouds. In particular, a novel approach is developed that fuses multiple depth maps derived by stereo matching to generate high-quality 3D maps. By learning critical configurations of stereo pairs from sample LiDAR data, we rank the image pairs based on the proximity of the results to the sample data. Multiple depth maps derived from individual image pairs are fused with an adaptive 3D median filter that considers the image spectral similarities. We demonstrate that the proposed adaptive median filter generally delivers better results in general as compared to normal median filter, and achieved an accuracy of improvement of 0.36 meters RMSE in the best case. Results and analysis are introduced in detail.


title: Analysis of critical parameters of satellite stereo image for 3D reconstruction and mapping

abstract: Although nowadays advanced dense image matching (DIM) algorithms are able to produce LiDAR (Light Detection And Ranging) comparable dense point clouds from satellite stereo images, the accuracy and completeness of such point clouds heavily depend on the geometric parameters of the satellite stereo images. The intersection angle between two images are normally seen as the most important one in stereo data acquisition, as the state-of-the-art DIM algorithms work best on narrow baseline (smaller intersection angle) stereos (E.g. Semi-Global Matching regards 15-25 degrees as good intersection angle). This factor is in line with the traditional aerial photogrammetry configuration, as the intersection angle directly relates to the base-high ratio and texture distortion in the parallax direction, thus both affecting the horizontal and vertical accuracy. However, our experiments found that even with very similar (and good) intersection angles, the same DIM algorithm applied on different stereo pairs (of the same area) produced point clouds with dramatically different accuracy as compared to the ground truth LiDAR data. This raises a very practical question that is often asked by practitioners: what factors constitute a good satellite stereo pair, such that it produces accurate and optimal results for mapping purpose? In this work, we provide a comprehensive analysis on this matter by performing stereo matching over 1,000 satellite stereo pairs with different acquisition parameters including their intersection angles, off-nadir angles, sun elevation & azimuth angles, as well as time differences, thus to offer a thorough answer to this question. This work will potentially provide a valuable reference to researchers working on multi-view satellite image reconstruction, as well as industrial practitioners minimizing costs for high-quality large-scale mapping.


title: Full Transport Model of GW170817-Like Disk Produces a Blue Kilonova

abstract: The 2017 detection of the inspiral and merger of two neutron stars in gravitational waves and gamma rays was accompanied by a quickly-reddening transient. Such a transient was predicted to occur following a rapid neutron capture (r-process) nucleosynthesis event, which synthesizes neutron-rich, radioactive nuclei and can take place in both dynamical ejecta and in the wind driven off the accretion torus formed after a neutron star merger. We present the first three-dimensional general relativistic, full transport neutrino radiation magnetohydrodynamics (GRRMHD) simulations of the black hole-accretion disk-wind system produced by the GW170817 merger. We show that the small but non-negligible optical depths lead to neutrino transport globally coupling the disk electron fraction, which we capture by solving the transport equation with a Monte Carlo method. The resulting absorption drives up the electron fraction in a structured, continuous outflow, with electron fraction as high as $Y_e\sim 0.4$ in the extreme polar region. We show via nuclear reaction network and radiative transfer calculations that nucleosynthesis in the disk wind will produce a blue kilonova.


title: Dueling Decoders: Regularizing Variational Autoencoder Latent Spaces

abstract: Variational autoencoders learn unsupervised data representations, but these models frequently converge to minima that fail to preserve meaningful semantic information. For example, variational autoencoders with autoregressive decoders often collapse into autodecoders, where they learn to ignore the encoder input. In this work, we demonstrate that adding an auxiliary decoder to regularize the latent space can prevent this collapse, but successful auxiliary decoding tasks are domain dependent. Auxiliary decoders can increase the amount of semantic information encoded in the latent space and visible in the reconstructions. The semantic information in the variational autoencoder's representation is only weakly correlated with its rate, distortion, or evidence lower bound. Compared to other popular strategies that modify the training objective, our regularization of the latent space generally increased the semantic information content.


title: Incentive Design for Efficient Federated Learning in Mobile Networks: A Contract Theory Approach

abstract: To strengthen data privacy and security, federated learning as an emerging machine learning technique is proposed to enable large-scale nodes, e.g., mobile devices, to distributedly train and globally share models without revealing their local data. This technique can not only significantly improve privacy protection for mobile devices, but also ensure good performance of the trained results collectively. Currently, most the existing studies focus on optimizing federated learning algorithms to improve model training performance. However, incentive mechanisms to motivate the mobile devices to join model training have been largely overlooked. The mobile devices suffer from considerable overhead in terms of computation and communication during the federated model training process. Without well-designed incentive, self-interested mobile devices will be unwilling to join federated learning tasks, which hinders the adoption of federated learning. To bridge this gap, in this paper, we adopt the contract theory to design an effective incentive mechanism for simulating the mobile devices with high-quality (i.e., high-accuracy) data to participate in federated learning. Numerical results demonstrate that the proposed mechanism is efficient for federated learning with improved learning accuracy.


title: Elastic Weyl points and surface arc states in 3D structures

abstract: The study of Weyl points in electronic systems has inspired many recent researches in classical systems such as photonic and acoustic lattices. Here we show how the Weyl physics can also inspire the design of novel elastic structures. We construct a single-phase 3D structure, an analogue of the AA-stacked honeycomb lattice, and predict the existence of Weyl points with opposite topological charges (${\pm 1}$), elastic Fermi arcs, and the associated gapless topologically protected surface states. We employ full-scale numerical simulations on the elastic 3D structure, and present a clear visualization of topological surface states that are directional and robust. Such designed lattices can pave the way for novel vibration control and energy harvesting on structures that are ubiquitous in many engineering applications.


title: Multilinear Compressive Learning

abstract: Compressive Learning is an emerging topic that combines signal acquisition via compressive sensing and machine learning to perform inference tasks directly on a small number of measurements. Many data modalities naturally have a multi-dimensional or tensorial format, with each dimension or tensor mode representing different features such as the spatial and temporal information in video sequences or the spatial and spectral information in hyperspectral images. However, in existing compressive learning frameworks, the compressive sensing component utilizes either random or learned linear projection on the vectorized signal to perform signal acquisition, thus discarding the multi-dimensional structure of the signals. In this paper, we propose Multilinear Compressive Learning, a framework that takes into account the tensorial nature of multi-dimensional signals in the acquisition step and builds the subsequent inference model on the structurally sensed measurements. Our theoretical complexity analysis shows that the proposed framework is more efficient compared to its vector-based counterpart in both memory and computation requirement. With extensive experiments, we also empirically show that our Multilinear Compressive Learning framework outperforms the vector-based framework in object classification and face recognition tasks, and scales favorably when the dimensionalities of the original signals increase, making it highly efficient for high-dimensional multi-dimensional signals.


title: Learning to Reconstruct 3D Manhattan Wireframes from a Single Image

abstract: In this paper, we propose a method to obtain a compact and accurate 3D wireframe representation from a single image by effectively exploiting global structural regularities. Our method trains a convolutional neural network to simultaneously detect salient junctions and straight lines, as well as predict their 3D depth and vanishing points. Compared with the state-of-the-art learning-based wireframe detection methods, our network is simpler and more unified, leading to better 2D wireframe detection. With global structural priors from parallelism, our method further reconstructs a full 3D wireframe model, a compact vector representation suitable for a variety of high-level vision tasks such as AR and CAD. We conduct extensive evaluations on a large synthetic dataset of urban scenes as well as real images. Our code and datasets have been made public at https://github.com/zhou13/shapeunity.


title: Deterministic Combinatorial Replacement Paths and Distance Sensitivity Oracles

abstract: In this work we derandomize two central results in graph algorithms, replacement paths and distance sensitivity oracles (DSOs) matching in both cases the running time of the randomized algorithms. For the replacement paths problem, let G = (V,E) be a directed unweighted graph with n vertices and m edges and let P be a shortest path from s to t in G. The {\sl replacement paths} problem is to find for every edge e \in P the shortest path from s to t avoiding e. Roditty and Zwick [ICALP 2005] obtained a randomized algorithm with running time of ~O(m \sqrt{n}). Here we provide the first deterministic algorithm for this problem, with the same ~O(m \sqrt{n}) time. For the problem of distance sensitivity oracles, let G = (V,E) be a directed graph with real-edge weights. An f-Sensitivity Distance Oracle (f-DSO) gets as input the graph G=(V,E) and a parameter f, preprocesses it into a data-structure, such that given a query (s,t,F) with s,t \in V and F \subseteq E \cup V, |F| \le f being a set of at most f edges or vertices (failures), the query algorithm efficiently computes the distance from s to t in the graph G \setminus F ({\sl i.e.}, the distance from s to t in the graph G after removing from it the failing edges and vertices F). For weighted graphs with real edge weights, Weimann and Yuster [FOCS 2010] presented a combinatorial randomized f-DSO with ~O(mn^{4-\alpha}) preprocessing time and subquadratic ~O(n^{2-2(1-\alpha)/f}) query time for every value of 0 < \alpha < 1. We derandomize this result and present a combinatorial deterministic f-DSO with the same asymptotic preprocessing and query time.


title: Exact and approximate mechanisms for pure bending of sheets

abstract: Direct measurement of the moment-curvature response of sheets or wires up to high curvatures can aid in modeling creasing, pleating, and other forming operations. We consider theoretical issues related to the geometry of pure bending. We present a linkage design that, for homogeneous deformation of the sample, results in an exact pure bending state up to arbitrarily high curvature, by combining a cochleoidal trajectory with an angle-doubling mechanism. The full mechanism is single degree of freedom and position controlled, a desirable feature for measurement of soft materials. We also present an optimal approximation to the cochleoid using a more easily implemented circular trajectory, and compare this with an existing commercial system for fabric testing. We quantify the error of the approximate test by calculating the deformation of an Euler \emph{elastica}, a structure with linear moment-curvature response. While the circular mechanism can approximate the exact boundary conditions quite well up to moderately large curvatures, the resulting curvature of the test sample still deviates significantly from homogeneous response. We also briefly discuss expectations of localization behavior in pure bending tests.


title: On the convergence of nuclear effective field theory with perturbative pions

abstract: The classic paper by Fleming, Mehen and Stewart cast doubts on the convergence of spin-triplet nucleon-nucleon partial wave scattering amplitudes when following the proposal of Kaplan, Savage and Wise to construct nuclear effective field theory around the unitary fermion limit with perturbative pion exchange. FMS identified the subclass of iterated one-pion exchange potential graphs as the cause of this poor convergence, which they showed persisted in the chiral limit. Theoretical tools are developed here to compute these Feynman graphs analytically to high order in all angular momentum channels simultaneously, examining the amplitudes computed to seven loops in the $L=J$ channels, and three loops in the coupled $L=J\pm1$ channels. One finds that there is nothing pathological about the perturbative expansion of a $1/r^3$ potential in general, and that the expansion converges satisfactorily in all partial waves except those with the lowest angular momentum, particularly the ${}^3P_0$ and the coupled ${}^3S_1-{}^3D_1$ channels. The results corroborate work by Birse, which suggests possible avenues to explore for improving the range of validity of the EFT expansion.


title: Aeronautical Ad Hoc Networking for the Internet-Above-The-Clouds

abstract: The engineering vision of relying on the ``smart sky" for supporting air traffic and the ``Internet above the clouds" for in-flight entertainment has become imperative for the future aircraft industry. Aeronautical ad hoc Networking (AANET) constitutes a compelling concept for providing broadband communications above clouds by extending the coverage of Air-to-Ground (A2G) networks to oceanic and remote airspace via autonomous and self-configured wireless networking amongst commercial passenger airplanes. The AANET concept may be viewed as a new member of the family of Mobile ad hoc Networks (MANETs) in action above the clouds. However, AANETs have more dynamic topologies, larger and more variable geographical network size, stricter security requirements and more hostile transmission conditions. These specific characteristics lead to more grave challenges in aircraft mobility modeling, aeronautical channel modeling and interference mitigation as well as in network scheduling and routing. This paper provides an overview of AANET solutions by characterizing the associated scenarios, requirements and challenges. Explicitly, the research addressing the key techniques of AANETs, such as their mobility models, network scheduling and routing, security and interference are reviewed. Furthermore, we also identify the remaining challenges associated with developing AANETs and present their prospective solutions as well as open issues. The design framework of AANETs and the key technical issues are investigated along with some recent research results. Furthermore, a range of performance metrics optimized in designing AANETs and a number of representative multi-objective optimization algorithms are outlined.


title: Designing of a magnetodielectric system in hybrid organic-inorganic framework, a perovskite layered phosphonate MnO3PC6H4-m-Br.H2O

abstract: The research on multiferrocity and magnetoelectric coupling in metal-organic system is rare. Very few hybrid organic-inorganic frameworks (HOIF) exhibit direct magnetoelectric coupling (coupling between spins and dipoles) and also restricted to particular COOH-based system. We show how one can design a hybrid system to obtain such coupling based on the rational design of the organic ligands. The layered phosphonate, MnO3PC6H5.H2O, consisting of perovskite layers stacked with organic phenyl layers, is used as a starting potential candidate. To introduce dipole moment, a closely related metal-phosphonate, MnO3PC6H4-m-Br.H2O is designed. For this purpose, this phosphonate is prepared from 3-bromophenylphosphonic acid that features one electronegative bromine atom directly attached on the aromatic ring in meta position, lowering the symmetry of precursor itself. Thus, bromobenzene moieties in MnO3PC6H4-m-Br.H2O induce a finite dipole moment. This new designed compound exhibits complex magnetism, as observed in layered alkyl chains MnO3PCnH2n+1.H2O materials, namely, 2D magnetic ordering around 20 K followed by weak ferromagnetic ordering below 12 K(T1) with a magnetic field (H)-induced transition around 25 kOe below T1. All these magnetic features are exactly captured in T and H-dependent dielectric constant, epsilon(T) and epsilon(H). This demonstrates direct magnetodielectric coupling in this designed hybrid and yields a new path to tune multiferroic ordering and magnetodielectric coupling.


title: Automatic Posterior Transformation for Likelihood-Free Inference

abstract: How can one perform Bayesian inference on stochastic simulators with intractable likelihoods? A recent approach is to learn the posterior from adaptively proposed simulations using neural network-based conditional density estimators. However, existing methods are limited to a narrow range of proposal distributions or require importance weighting that can limit performance in practice. Here we present automatic posterior transformation (APT), a new sequential neural posterior estimation method for simulation-based inference. APT can modify the posterior estimate using arbitrary, dynamically updated proposals, and is compatible with powerful flow-based density estimators. It is more flexible, scalable and efficient than previous simulation-based inference techniques. APT can operate directly on high-dimensional time series and image data, opening up new applications for likelihood-free inference.


title: TASI Lectures on Future Colliders

abstract: These lectures review the main motivations for future high-energy colliders, focusing on the understanding of electroweak symmetry breaking and on the search for physics beyond the Standard Model. The open questions and the challenges are common to all future projects; for concreteness, I will use studies of the potential of $\rm e^+e^-$ and pp circular colliders to provide examples of the anticipated physics reach.


title: Sequential training algorithm for neural networks

abstract: A sequential training method for large-scale feedforward neural networks is presented. Each layer of the neural network is decoupled and trained separately. After the training is completed for each layer, they are combined together. The performance of the network would be sub-optimal compared to the full network training if the optimal solution would be achieved. However, achieving the optimal solution for the full network would be infeasible or require long computing time. The proposed sequential approach reduces the required computer resources significantly and would have better convergences as a single layer is optimised for each optimisation step. The required modifications of existing algorithms to implement the sequential training are minimal. The performance is verified by a simple example.


title: Positioning aiding using LiDAR in GPS signal loss scenarios

abstract: In the presented scenario, an autonomous surface vehicle (ASV) equipped with a laser scanner navigates on a inland pathway surrounded and crossed by man-made structures such as bridges and locks. {GPS} receiver present on board experiences signal loss and multipath reflections in situation when the view of the sky is obscured by a bridge or tall buildings. In both cases, a potentially dangerous situation is provoked as the robot has no or inaccurate positioning data. A sensor data processing scheme is proposed where these gaps are smoothly filled in by positioning data generated from scan matching and registration of the laser data. This article shows preliminary results of positioning data improvement during trials in harbor-river environment.


title: Manifestation of the Purcell effect in current transport through a dot-cavity-QED system

abstract: We study the transport properties of a wire-dot system coupled to a cavity and a photon reservoir. Tuning the photon energy, Rabi-resonant states emerge and in turn resonant current peaks are observed. We demonstrate the effects of the cavity-photon reservoir coupling, the mean photon number in the reservoir, the electron-photon coupling and the photon polarization on the intraband transitions occurring between the Rabi-resonant states, and on the corresponding resonant current peaks. The Rabi-splitting can be controlled by the photon polarization and the electron-photon coupling strength. In the selected range of parameters, we observe the results of the Purcell effect enhancing the current peaks through the cavity by increasing the cavity-reservoir coupling, while they decrease with increasing the electron-photon coupling. In addition, the resonant current peaks are also sensitive to the mean number of photons in the reservoir.


title: The exponential-nonexponential transition in quantum decay as a phenomenon of interference in time domain of the decaying particle with itself

abstract: By using an exact analytical non-Hermitian approach in terms of resonance (quasinormal) states we express the decaying wave function as the sum of exponential and nonexponential decaying solutions to the time-dependent Schr\"odinger equation. We show that the exponential-nonexponential transition of decay at long times represents physically a process of interference in time domain of the decaying particle with itself.


title: On groups in which Engel sinks are cyclic

abstract: For an element $g$ of a group $G$, an Engel sink is a subset $\mathcal{E}(g)$ such that for every $ x\in G $ all sufficiently long commutators $ [x,g,g,\ldots,g] $ belong to $\mathcal{E}(g)$. We conjecture that if $G$ is a profinite group in which every element admits a sink that is a procyclic subgroup, then $G$ is procyclic-by-(locally nilpotent). We prove the conjecture in two cases -- when $G$ is a finite group, or a soluble pro-$p$ group.


title: Multi-field Inflation in High-Slope Potentials

abstract: We present two families of multi-field potentials that support inflation while satisfying the refined de Sitter and the distance swampland conjectures. Both families feature Planck-compatible phenomenology. The first is a helix-type potential, in a flat field-space metric, that satisfies the conjectures via a high turning rate. This model has a tensor-toscalar ratio close to, but below, the current experimental limits and small non-gaussianities. The second family, an example of orbital inflation, utilizes a negatively curved field metric to achieve prolonged inflation with nontrivial turning in the presence of a tachyonic direction. Although perturbations in this model undergo an exponential growth before horizon exit, it is always possible to match the measured amplitude of the power spectrum by lowering the scale of inflation if the turning rate is low enough. We identify a Planck-compatible region of parameter space in which the scale of inflation is above that of nucleosynthesis. Due to the rapid growth, this model predicts an exponentially suppressed value for the tensor-to-scalar ratio.


title: Complex polynomial Bohnenblust--Hille inequality with polynomial bounds

abstract: The Bohnenblust-Hille inequality and its variants have found applications in several areas of Mathematics and related fields. The control of the constants for the variant for complex $m$-homogeneous polynomials is of special interest for applications in Harmonic Analysis and Number Theory. Up to now, the best known estimates for its constants are dominated by $\kappa\left(1+\varepsilon\right) ^{m}$, where $\varepsilon>0$ is arbitrary and $\kappa>0$ depends on the choice of $\varepsilon$. For the special cases in which the number of variables in each monomial is bounded by some fixed number $M$, it has been shown that the optimal constant is dominated by a constant depending solely on $M$. In this note, based on a deep result of Bayart, we prove an inequality for any subset of the indices, showing how summability of arbitrary restrictions on monomials can be related to the combinatorial dimension associated with them.


title: A comprehensive study of speech separation: spectrogram vs waveform separation

abstract: Speech separation has been studied widely for single-channel close-talk microphone recordings over the past few years; developed solutions are mostly in frequency-domain. Recently, a raw audio waveform separation network (TasNet) is introduced for single-channel data, with achieving high Si-SNR (scale-invariant source-to-noise ratio) and SDR (source-to-distortion ratio) comparing against the state-of-the-art solution in frequency-domain. In this study, we incorporate effective components of the TasNet into a frequency-domain separation method. We compare both for alternative scenarios. We introduce a solution for directly optimizing the separation criterion in frequency-domain networks. In addition to speech separation objective and subjective measurements, we evaluate the separation performance on a speech recognition task as well. We study the speech separation problem for far-field data (more similar to naturalistic audio streams) and develop multi-channel solutions for both frequency and time-domain separators with utilizing spectral, spatial and speaker location information. For our experiments, we simulated multi-channel spatialized reverberate WSJ0-2mix dataset. Our experimental results show that spectrogram separation can achieve competitive performance with better network design. Multi-channel framework as well is shown to improve the single-channel performance relatively up to +35.5% and +46% in terms of WER and SDR, respectively.


title: Rethinking Atmospheric Turbulence Mitigation

abstract: State-of-the-art atmospheric turbulence image restoration methods utilize standard image processing tools such as optical flow, lucky region and blind deconvolution to restore the images. While promising results have been reported over the past decade, many of the methods are agnostic to the physical model that generates the distortion. In this paper, we revisit the turbulence restoration problem by analyzing the reference frame generation and the blind deconvolution steps in a typical restoration pipeline. By leveraging tools in large deviation theory, we rigorously prove the minimum number of frames required to generate a reliable reference for both static and dynamic scenes. We discuss how a turbulence agnostic model can lead to potential flaws, and how to configure a simple spatial-temporal non-local weighted averaging method to generate references. For blind deconvolution, we present a new data-driven prior by analyzing the distributions of the point spread functions. We demonstrate how a simple prior can outperform state-of-the-art blind deconvolution methods.


title: LR-GLM: High-Dimensional Bayesian Inference Using Low-Rank Data Approximations

abstract: Due to the ease of modern data collection, applied statisticians often have access to a large set of covariates that they wish to relate to some observed outcome. Generalized linear models (GLMs) offer a particularly interpretable framework for such an analysis. In these high-dimensional problems, the number of covariates is often large relative to the number of observations, so we face non-trivial inferential uncertainty; a Bayesian approach allows coherent quantification of this uncertainty. Unfortunately, existing methods for Bayesian inference in GLMs require running times roughly cubic in parameter dimension, and so are limited to settings with at most tens of thousand parameters. We propose to reduce time and memory costs with a low-rank approximation of the data in an approach we call LR-GLM. When used with the Laplace approximation or Markov chain Monte Carlo, LR-GLM provides a full Bayesian posterior approximation and admits running times reduced by a full factor of the parameter dimension. We rigorously establish the quality of our approximation and show how the choice of rank allows a tunable computational-statistical trade-off. Experiments support our theory and demonstrate the efficacy of LR-GLM on real large-scale datasets.


title: Classification of some vertex operator algebras of rank 3

abstract: We discuss the classification of strongly regular vertex operator algebras (VOAs) with exactly three simple modules whose character vector satisfies a monic modular linear differential equation with irreducible monodromy. Our Main Theorem provides a classification of all such VOAs in the form of one infinite family of affine VOAs, one individual affine algebra and two Virasoro algebras, together with a family of eleven exceptional character vectors and associated data that we call the $U$-series. We prove that there are at least $15$ VOAs in the $U$-series occurring as commutants in a Schellekens list holomorphic VOA. These include the affine algebra $E_{8,2}$ and H\"ohn's Baby Monster VOA $\mathbf{VB}^\natural_{(0)}$ but the other $13$ seem to be new. The idea in the proof of our Main Theorem is to exploit properties of a family of vector-valued modular forms with rational functions as Fourier coefficients, which solves a family of modular linear differential equations in terms of generalized hypergeometric series.


title: Enforcing constraints for time series prediction in supervised, unsupervised and reinforcement learning

abstract: We assume that we are given a time series of data from a dynamical system and our task is to learn the flow map of the dynamical system. We present a collection of results on how to enforce constraints coming from the dynamical system in order to accelerate the training of deep neural networks to represent the flow map of the system as well as increase their predictive ability. In particular, we provide ways to enforce constraints during training for all three major modes of learning, namely supervised, unsupervised and reinforcement learning. In general, the dynamic constraints need to include terms which are analogous to memory terms in model reduction formalisms. Such memory terms act as a restoring force which corrects the errors committed by the learned flow map during prediction. For supervised learning, the constraints are added to the objective function. For the case of unsupervised learning, in particular generative adversarial networks, the constraints are introduced by augmenting the input of the discriminator. Finally, for the case of reinforcement learning and in particular actor-critic methods, the constraints are added to the reward function. In addition, for the reinforcement learning case, we present a novel approach based on homotopy of the action-value function in order to stabilize and accelerate training. We use numerical results for the Lorenz system to illustrate the various constructions.


title: ACE of Space: Estimating Genetic Components of High-Dimensional Imaging Data

abstract: It is of great interest to quantify the contributions of genetic variation to brain structure and function, which are usually measured by high-dimensional imaging data (e.g., magnetic resonance imaging). In addition to the variance, the covariance patterns in the genetic effects of a functional phenotype are of biological importance, and covariance patterns have been linked to psychiatric disorders. The aim of this paper is to develop a scalable method to estimate heritability and the non-stationary covariance components in high-dimensional imaging data from twin studies. Our motivating example is from the Human Connectome Project (HCP). Several major big-data challenges arise from estimating the genetic and environmental covariance functions of functional phenotypes extracted from imaging data, such as cortical thickness with 60,000 vertices. Notably, truncating to positive eigenvalues and their eigenfunctions from unconstrained estimators can result in large bias. This motivated our development of a novel estimator ensuring positive semidefiniteness. Simulation studies demonstrate large improvements over existing approaches, both with respect to heritability estimates and covariance estimation. We applied the proposed method to cortical thickness data from the HCP. Our analysis suggests fine-scale differences in covariance patterns, identifying locations in which genetic control is correlated with large areas of the brain and locations where it is highly localized.


title: 3DViewGraph: Learning Global Features for 3D Shapes from A Graph of Unordered Views with Attention

abstract: Learning global features by aggregating information over multiple views has been shown to be effective for 3D shape analysis. For view aggregation in deep learning models, pooling has been applied extensively. However, pooling leads to a loss of the content within views, and the spatial relationship among views, which limits the discriminability of learned features. We propose 3DViewGraph to resolve this issue, which learns 3D global features by more effectively aggregating unordered views with attention. Specifically, unordered views taken around a shape are regarded as view nodes on a view graph. 3DViewGraph first learns a novel latent semantic mapping to project low-level view features into meaningful latent semantic embeddings in a lower dimensional space, which is spanned by latent semantic patterns. Then, the content and spatial information of each pair of view nodes are encoded by a novel spatial pattern correlation, where the correlation is computed among latent semantic patterns. Finally, all spatial pattern correlations are integrated with attention weights learned by a novel attention mechanism. This further increases the discriminability of learned features by highlighting the unordered view nodes with distinctive characteristics and depressing the ones with appearance ambiguity. We show that 3DViewGraph outperforms state-of-the-art methods under three large-scale benchmarks.


title: Story Ending Prediction by Transferable BERT

abstract: Recent advances, such as GPT and BERT, have shown success in incorporating a pre-trained transformer language model and fine-tuning operation to improve downstream NLP systems. However, this framework still has some fundamental problems in effectively incorporating supervised knowledge from other related tasks. In this study, we investigate a transferable BERT (TransBERT) training framework, which can transfer not only general language knowledge from large-scale unlabeled data but also specific kinds of knowledge from various semantically related supervised tasks, for a target task. Particularly, we propose utilizing three kinds of transfer tasks, including natural language inference, sentiment classification, and next action prediction, to further train BERT based on a pre-trained model. This enables the model to get a better initialization for the target task. We take story ending prediction as the target task to conduct experiments. The final result, an accuracy of 91.8%, dramatically outperforms previous state-of-the-art baseline methods. Several comparative experiments give some helpful suggestions on how to select transfer tasks. Error analysis shows what are the strength and weakness of BERT-based models for story ending prediction.


title: Object Rearrangement with Nested Nonprehensile Manipulation Actions

abstract: This paper considers the problem of rearrangement planning, i.e finding a sequence of manipulation actions that displace multiple objects from an initial configuration to a given goal configuration. Rearrangement is a critical skill for robots so that they can effectively operate in confined spaces that contain clutter. Examples of tasks that require rearrangement include packing objects inside a bin, wherein objects need to lay according to a predefined pattern. In tight bins, collision-free grasps are often unavailable. Nonprehensile actions, such as pushing and sliding, are preferred because they can be performed using minimalistic end-effectors that can easily be inserted in the bin. Rearrangement with nonprehensile actions is a challenging problem as it requires reasoning about object interactions in a combinatorially large configuration space of multiple objects. This work revisits several existing rearrangement planning techniques and introduces a new one that exploits nested nonprehensile actions by pushing several similar objects simultaneously along the same path, which removes the need to rearrange each object individually. Experiments in simulation and using a real Kuka robotic arm show the ability of the proposed approach to solve difficult rearrangement tasks while reducing the length of the end-effector's trajectories.


title: Parts4Feature: Learning 3D Global Features from Generally Semantic Parts in Multiple Views

abstract: Deep learning has achieved remarkable results in 3D shape analysis by learning global shape features from the pixel-level over multiple views. Previous methods, however, compute low-level features for entire views without considering part-level information. In contrast, we propose a deep neural network, called Parts4Feature, to learn 3D global features from part-level information in multiple views. We introduce a novel definition of generally semantic parts, which Parts4Feature learns to detect in multiple views from different 3D shape segmentation benchmarks. A key idea of our architecture is that it transfers the ability to detect semantically meaningful parts in multiple views to learn 3D global features. Parts4Feature achieves this by combining a local part detection branch and a global feature learning branch with a shared region proposal module. The global feature learning branch aggregates the detected parts in terms of learned part patterns with a novel multi-attention mechanism, while the region proposal module enables locally and globally discriminative information to be promoted by each other. We demonstrate that Parts4Feature outperforms the state-of-the-art under three large-scale 3D shape benchmarks.


title: Enveloping algebras with just infinite Gelfand-Kirillov dimension

abstract: Let $\mf g$ be the Witt algebra or the positive Witt algebra. It is well known that the enveloping algebra $U(\mf g )$ has intermediate growth and thus infinite Gelfand-Kirillov (GK-) dimension. We prove that the GK-dimension of $U(\mf g)$ is {\em just infinite} in the sense that any proper quotient of $U(\mf g)$ has polynomial growth. This proves a conjecture of Petukhov and the second named author for the positive Witt algebra. We also establish the corresponding results for quotients of the symmetric algebra $S(\mf g)$ by proper Poisson ideals. In fact, we prove more generally that any central quotient of the universal enveloping algebra of the Virasoro algebra has just infinite GK-dimension. We give several applications. In particular, we easily compute the annihilators of Verma modules over the Virasoro algebra.

categories: math.RA math-ph math.AG math.MP math.QA math.RT


title: Cross-referencing using Fine-grained Topic Modeling

abstract: Cross-referencing, which links passages of text to other related passages, can be a valuable study aid for facilitating comprehension of a text. However, cross-referencing requires first, a comprehensive thematic knowledge of the entire corpus, and second, a focused search through the corpus specifically to find such useful connections. Due to this, cross-reference resources are prohibitively expensive and exist only for the most well-studied texts (e.g. religious texts). We develop a topic-based system for automatically producing candidate cross-references which can be easily verified by human annotators. Our system utilizes fine-grained topic modeling with thousands of highly nuanced and specific topics to identify verse pairs which are topically related. We demonstrate that our system can be cost effective compared to having annotators acquire the expertise necessary to produce cross-reference resources unaided.


title: Supersymmetric generalized power functions

abstract: Complex-valued functions defined on a finite interval $[a,b]$ generalizing power functions of the type $(x-x_0)^n$ for $n\geq 0$ are studied. These functions called $\Phi$-generalized powers, $\Phi$ being a given nonzero complex-valued function on the interval, were considered to contruct a general solution representation of the Sturm-Liouville equation in terms of the spectral parameter \cite{kravchenko2008, kravporter2010}. The $\Phi$-generalized powers can be considered as a natural basis functions for the one-dimensional supersymmetric quantum mechanics systems taking $\Phi=\psi_0^2$, where the function $\psi_0(x)$ is the ground state wave function of one of the supersymmetric scalar Hamiltonians. Several properties are obtained such as $\Phi$-symmetric conjugate and antisymmetry of the $\Phi$-generalized powers, a supersymmetric binomial identity for these functions, a supersymmetric Pythagorean elliptic (hyperbolic) identity involving four $\Phi$-trigonometric ($\Phi$-hyperbolic) functions as well as a supersymmetric Taylor series expressed in terms of the $\Phi$-derivatives. We show that the first $n$ $\Phi$-generalized powers are a fundamental set of solutions associated with a nonconstant coefficients homogeneous linear ordinary differential equations of order $n+1$. Finally, we present a general solution representation of the stationary Schr\"odinger equation in terms of geometric series where the Volterra compositions of the first type is considered.


title: Discovery of Algebraic Reynolds-Stress Models Using Sparse Symbolic Regression

abstract: ** This article is published (open-access). ** A novel deterministic symbolic regression method SpaRTA (Sparse Regression of Turbulent Stress Anisotropy) is introduced to infer algebraic stress models for the closure of RANS equations directly from high-fidelity LES or DNS data. The models are written as tensor polynomials and are built from a library of candidate functions. The machine-learning method is based on elastic net regularisation which promotes sparsity of the inferred models. By being data-driven the method relaxes assumptions commonly made in the process of model development. Model-discovery and cross-validation is performed for three cases of separating flows, i.e. periodic hills ($Re$=10595), converging-diverging channel ($Re$=12600) and curved backward-facing step ($Re$=13700). The predictions of the discovered models are significantly improved over the $k$-$\omega$ SST also for a true prediction of the flow over periodic hills at $Re$=37000. This study shows a systematic assessment of SpaRTA for rapid machine-learning of robust corrections for standard RANS turbulence models.


title: HALLS: An Energy-Efficient Highly Adaptable Last Level STT-RAM Cache for Multicore Systems

abstract: Spin-Transfer Torque RAM (STT-RAM) is widely considered a promising alternative to SRAM in the memory hierarchy due to STT-RAM's non-volatility, low leakage power, high density, and fast read speed. The STT-RAM's small feature size is particularly desirable for the last-level cache (LLC), which typically consumes a large area of silicon die. However, long write latency and high write energy still remain challenges of implementing STT-RAMs in the CPU cache. An increasingly popular method for addressing this challenge involves trading off the non-volatility for reduced write speed and write energy by relaxing the STT-RAM's data retention time. However, in order to maximize energy saving potential, the cache configurations, including STT-RAM's retention time, must be dynamically adapted to executing applications' variable memory needs. In this paper, we propose a highly adaptable last level STT-RAM cache (HALLS) that allows the LLC configurations and retention time to be adapted to applications' runtime execution requirements. We also propose low-overhead runtime tuning algorithms to dynamically determine the best (lowest energy) cache configurations and retention times for executing applications. Compared to prior work, HALLS reduced the average energy consumption by 60.57% in a quad-core system, while introducing marginal latency overhead.


title: SplitNet: Sim2Sim and Task2Task Transfer for Embodied Visual Navigation

abstract: We propose SplitNet, a method for decoupling visual perception and policy learning. By incorporating auxiliary tasks and selective learning of portions of the model, we explicitly decompose the learning objectives for visual navigation into perceiving the world and acting on that perception. We show dramatic improvements over baseline models on transferring between simulators, an encouraging step towards Sim2Real. Additionally, SplitNet generalizes better to unseen environments from the same simulator and transfers faster and more effectively to novel embodied navigation tasks. Further, given only a small sample from a target domain, SplitNet can match the performance of traditional end-to-end pipelines which receive the entire dataset. Code is available https://github.com/facebookresearch/splitnet


title: Phase diagrams of spin-1 Ising model with bilinear and quadrupolar interactions under magnetic field. Two-particle cluster approximation

abstract: The spin-1 Ising model with bilinear and quadrupolar short-range interactions under magnetic field is investigated within the two-particle cluster approximation. It is shown that for those values of the quadrupolar interaction when at zero magnetic field the system undergoes a temperature phase transition between quadrupolar and paramagnetic phases, a triple point may exist in the temperature vs magnetic field phase diagrams, necessarily along with a critical point. It is also shown that the critical points in the temperature vs magnetic field phase diagrams of the investigated model can be of three types.


title: Bondi-Metzner-Sachs invariance and electric-magnetic duality

abstract: We exhibit a Hamiltonian formulation, both for electromagnetism and gravitation, in which it is not required that the Bondi "news" vanish, but only that the incoming news be equal to the outgoing ones. This requirement is implemented by defining the fields on a two-sheeted hyperbolic surface, which we term "the hourglass". It is a spacelike deformation of the complete lightcone. On it one approaches asymptotically (null) past and future infinity while remaining at a fixed (hyperbolic) time, by going to large spatial distances on its two sheets. The Hamiltonian formulation and -- in particular -- a conserved angular momentum, can only be constructed if one brings in both, the electric and magnetic BMS charges, together with their canonically conjugate "memories". This reveals a close interplay between the BMS and electric-magnetic duality symmetries.


title: Learning Perspective Undistortion of Portraits

abstract: Near-range portrait photographs often contain perspective distortion artifacts that bias human perception and challenge both facial recognition and reconstruction techniques. We present the first deep learning based approach to remove such artifacts from unconstrained portraits. In contrast to the previous state-of-the-art approach, our method handles even portraits with extreme perspective distortion, as we avoid the inaccurate and error-prone step of first fitting a 3D face model. Instead, we predict a distortion correction flow map that encodes a per-pixel displacement that removes distortion artifacts when applied to the input image. Our method also automatically infers missing facial features, i.e. occluded ears caused by strong perspective distortion, with coherent details. We demonstrate that our approach significantly outperforms the previous state-of-the-art both qualitatively and quantitatively, particularly for portraits with extreme perspective distortion or facial expressions. We further show that our technique benefits a number of fundamental tasks, significantly improving the accuracy of both face recognition and 3D reconstruction and enables a novel camera calibration technique from a single portrait. Moreover, we also build the first perspective portrait database with a large diversity in identities, expression and poses, which will benefit the related research in this area.


title: Neutrinoless double beta decay in effective field theory

abstract: We discuss the contributions of lepton-number-violating sources to neutrinoless double beta decay ($0\nu\beta\beta$). Assuming that these sources arise at scales well above the electroweak scale, they can be described within an effective field theory. Here, we outline the steps required to express the $0\nu\beta\beta$ half-life in terms of the effective interactions, focusing on the dimension-five operator that induces a Majorana mass for the neutrinos. This process involves the evolution of the operators down to scales of a few GeV where they can be matched onto Chiral Perturbation Theory. The resulting Chiral Lagrangian can then used be to derive the lepton-number violating potential, which, in combination with many-body methods, gives the $0\nu\beta\beta$ half-life. We will show that consistent renormalization requires the inclusion of a new contact interaction at leading order in this potential. We also briefly comment on the constraints that can be set on the operators appearing beyond dimension five.


title: Degree bounds for Gr\"{o}bner bases of modules

abstract: Let $F$ be a non-negatively graded free module over a polynomial ring $\mathbb{K}[x_1,\dots,x_n]$ generated by $m$ basis elements. Let $M$ be a submodule of $F$ generated by elements in $F$ with degrees bounded by $D$ and dim $F/M$=$r$. We prove that if $M$ is graded, the degree of the reduced Gr\"{o}bner basis of $M$ for any term order is bounded by $2\left[1/2((Dm)^{n-r}m+D) \right]^{2^{r-1}}$. If $M$ is not graded, the bound is $2\left[1/2((Dm)^{(n-r)^2}m+D) \right]^{2^{r}}$. This is a generalization of Dub\'{e}(1990) and Mayr-Ritscher(2013)'s bounds for ideals in a polynomial ring.


title: Smooth quasi-developable surfaces bounded by smooth curves

abstract: Computing a quasi-developable strip surface bounded by design curves finds wide industrial applications. Existing methods compute discrete surfaces composed of developable lines connecting sampling points on input curves which are not adequate for generating smooth quasi-developable surfaces. We propose the first method which is capable of exploring the full solution space of continuous input curves to compute a smooth quasi-developable ruled surface with as large developability as possible. The resulting surface is exactly bounded by the input smooth curves and is guaranteed to have no self-intersections. The main contribution is a variational approach to compute a continuous mapping of parameters of input curves by minimizing a function evaluating surface developability. Moreover, we also present an algorithm to represent a resulting surface as a B-spline surface when input curves are B-spline curves.


title: An Adaptive Global-Local Approach for Phase-Field Modeling of Anisotropic Brittle Fracture

abstract: This work addresses an efficient Global-Local approach supplemented with predictor-corrector adaptivity applied to anisotropic phase-field brittle fracture. The phase-field formulation is used to resolve the sharp crack surface topology on the anisotropic/non-uniform local state in the regularized concept. To resolve the crack phase-field by a given single preferred direction, second-order structural tensors are imposed to both the bulk and crack surface density functions. Accordingly, a split in tension and compression modes in anisotropic materials is considered. A Global-Local formulation is proposed, in which the full displacement/phase-field problem is solved on a lower (local) scale, while dealing with a purely linear elastic problem on an upper (global) scale. Robin-type boundary conditions are introduced to relax the stiff local response at the global scale and enhancing its stabilization. Another important aspect of this contribution is the development of an adaptive Global-Local approach, where a predictor-corrector scheme is designed in which the local domains are dynamically updated during the computation. To cope with different finite element discretizations at the interface between the two nested scales, a non-matching dual mortar method is formulated. Hence, more regularity is achieved on the interface. Several numerical results substantiate our developments.


title: Geometric Measures of Information for Quantum State Characterization

abstract: We analyze the geometry of a joint distribution over a set of discrete random variables. We briefly review Shannon's entropy, conditional entropy, mutual information and conditional mutual information. We review the entropic information distance formula of Rokhlin and Rajski. We then define an analogous information area. We motivate this definition and discuss its properties. We extend this definition to higher-dimensional volumes. We briefly discuss the potential utility for these geometric measures in quantum information processing.


title: Fluctuations of Electric Fields in the Active Site of the Enzyme Ketosteroid Isomerase

abstract: We report the effect of conformational dynamics on the fluctuations of electric fields in the active site of the enzyme Ketosteroid Isomerase (KSI). While KSI is considered rigid with little conformational variation to support different stages of the catalytic cycle, we show that KSI utilizes cooperative side chain motions of the entire protein scaffold outside the active site, which contribute negligibly to the electric fields on the substrate, by progressively stabilizing electric field contributions by particular active site residues at different timescales. The design of synthetic enzymes could benefit from strategies that can take advantage of the dynamics by using electric fields fluctuations as a guide.


title: Properties of Quantum Reactivity for a Multipartite State

abstract: We discuss the properties of quantum state reactivity as a measure for quantum correlation. This information geometry-based definition is a generalization of the two qubit construction of Schumacher to multipartite quantum states. It requires a generalization of information distance to information areas as well as to higher-dimensional volumes. The reactivity is defined in the usual chemistry way as a ratio of surface area to volume. The reactivity is an average over all detector settings. We show that this measure posses the key features required for a measure of quantum correlation. We show that it is invariant under local unitary transformations, non-increasing under local operations and classical communication, and monotonic. Its maximum bound can't be obtained using only classical correlation. Furthermore, reactivity is an analytic function of measurement probabilities and easily extendable to higher multipartite states.


title: Pore-scale direct numerical simulation of Haines jumps in a porous media model

abstract: Direct numerical simulations are presented for a porous media model consisting of two immiscible fluids, an invading and defending phase, in a two-dimensional micro-geometry filled with randomly sized and randomly distributed cylinders. First, interface instability and penetration modes are studied when varying the wetting features of a single pore in the porous medium. It is found that the displacement patterns not only change with the capillary number, as previously observed, but also are a function of the contact angle, even for a viscosity ratio of unity. This is an important conclusion suggesting that capillary number and viscosity ratio alone cannot completely describe the pore-scale displacement. Second, rapid pore-scale displacement is considered, where the displacements are accompanied by sudden interface jumps from one site to another, known as Haines jumps. The characteristic time and length scales of a Haines jump are examined to better understand the transient dynamics of the jump. We then focus on analyzing the Haines jump in a simple pore configuration where cylinders of equal size are placed at the vertices of equilateral triangles. We use this geometry to provide more insight into the effect of the contact angle at which the Haines jump is predicted.


title: Global large solutions for the Navier-Stokes equations with the Coriolis force

abstract: In this paper, we construct a class of global large solution to the three-dimensional Navier-Stokes equations with the Coriolis force in critical Fourier-Besov space $\dot{FB}^{2-\frac{3}{p}}_{p,r}(\mathbb{R}^3)$. In fact, our choice of special initial data $u_0$ can be arbitrarily large in $\dot{FB}^{s}_{p,r}(\mathbb{R}^3)$ for any $s\in\mathbb{R}$ and $1\leq p,r\leq \infty$.


title: Geometric model of the fracture as a manifold immersed in porous media

abstract: In this work, we analyze the flow filtration process of slightly compressible fluids in porous media containing man made fractures with complex geometries. We model the coupled fracture-porous media system where the linear Darcy flow is considered in porous media and the nonlinear Forchheimer equation is used inside the fracture. We develop a model to examine the flow inside fractures with complex geometries and variable thickness, on a Riemannian manifold. The fracture is represented as the normal variation of a surface immersed in $\mathbb{R}^3$. Using operators of Laplace Beltrami type and geometric identities, we model an equation that describes the flow in the fracture. A reduced model is obtained as a low dimensional BVP. We then couple the model with the porous media. Theoretical and numerical analysis have been performed to compare the solutions between the original geometric model and the reduced model in reservoirs containing fractures with complex geometries. We prove that the two solutions are close, and therefore, the reduced model can be effectively used in large scale simulators for long and thin fractures with complicated geometry.


title: Distribution Electricity Pricing under Uncertainty

abstract: Distribution locational marginal prices (DLMPs) facilitate the efficient operation of low-voltage electric power distribution systems. We propose an approach to internalize the stochasticity of renewable distributed energy resources (DERs) and risk tolerance of the distribution system operator in DLMP computations. This is achieved by means of applying conic duality to a chance-constrained AC optimal power flow. We show that the resulting DLMPs consist of the terms that allow to itemize the prices for the active and reactive power production, balancing regulation, and voltage support provided. Finally, we prove the proposed DLMPs constitute a competitive equilibrium, which can be leveraged for designing a distribution electricity market, and show that imposing chance constraints on voltage limits distorts the equilibrium.


title: Brief Notes and History Computing in Mexico during 50 years

abstract: The history of computing in Mexico can not be thought without the name of Prof. Harold V. McIntosh (1929-2015). For almost 50 years, in Mexico he contributed to the development of computer science with wide international recognition. Approximately in 1964, McIntosh began working in the Physics Department of the Advanced Studies Center (CIEA) of the National Polytechnic Institute (IPN), now called CINVESTAV. In 1965, at the National Center of Calculus (CeNaC), he was a founding member of the Master in Computing, first in Latin America. With the support of Mario Baez Camargo and Enrique Melrose, McIntosh continues his research of Martin-Baltimore Computer Center and University of Florida at IBM 709.


title: When Social Sensing Meets Edge Computing: Vision and Challenges

abstract: This paper overviews the state of the art, research challenges, and future opportunities in an emerging research direction: Social Sensing based Edge Computing (SSEC). Social sensing has emerged as a new sensing application paradigm where measurements about the physical world are collected from humans or from devices on their behalf. The advent of edge computing pushes the frontier of computation, service, and data along the cloud-to-things continuum. The merging of these two technical trends generates a set of new research challenges that need to be addressed. In this paper, we first define the new SSEC paradigm that is motivated by a few underlying technology trends. We then present a few representative real-world case studies of SSEC applications and several key research challenges that exist in those applications. Finally, we envision a few exciting research directions in future SSEC. We hope this paper will stimulate discussions of this emerging research direction in the community.


title: Multinomial Distribution Learning for Effective Neural Architecture Search

abstract: Architectures obtained by Neural Architecture Search (NAS) have achieved highly competitive performance in various computer vision tasks. However, the prohibitive computation demand of forward-backward propagation in deep neural networks and searching algorithms makes it difficult to apply NAS in practice. In this paper, we propose a Multinomial Distribution Learning for extremely effective NAS,which considers the search space as a joint multinomial distribution, i.e., the operation between two nodes is sampled from this distribution, and the optimal network structure is obtained by the operations with the most likely probability in this distribution. Therefore, NAS can be transformed to a multinomial distribution learning problem, i.e., the distribution is optimized to have a high expectation of the performance. Besides, a hypothesis that the performance ranking is consistent in every training epoch is proposed and demonstrated to further accelerate the learning process. Experiments on CIFAR10 and ImageNet demonstrate the effectiveness of our method. On CIFAR-10, the structure searched by our method achieves 2.55% test error, while being 6.0x (only 4 GPU hours on GTX1080Ti) faster compared with state-of-the-art NAS algorithms. On ImageNet, our model achieves 75.2% top1 accuracy under MobileNet settings (MobileNet V1/V2), while being 1.2x faster with measured GPU latency. Test code with pre-trained models are available at https://github.com/tanglang96/MDENAS


title: Factor Models for High-Dimensional Tensor Time Series

abstract: Large tensor (multi-dimensional array) data are now routinely collected in a wide range of applications, due to modern data collection capabilities. Often such observations are taken over time, forming tensor time series. In this paper we present a factor model approach for analyzing high-dimensional dynamic tensor time series and multi-category dynamic transport networks. Two estimation procedures along with their theoretical properties and simulation results are presented. Two applications are used to illustrate the model and its interpretations.


title: Lyapunov Exponent of Rank One Matrices: Ergodic Formula and Inapproximability of the Optimal Distribution

abstract: The Lyapunov exponent corresponding to a set of square matrices $\mathcal{A} = \{A_1, \dots, A_n \}$ and a probability distribution $p$ over $\{1, \dots, n\}$ is $\lambda(\mathcal{A},p) := \lim_{k \to \infty} \frac{1}{k} \,\mathbb{E} \log \|A_{\sigma_k} \cdots A_{\sigma_2}A_{\sigma_1}\|$, where $\sigma_i$ are i.i.d. according to $p$. This quantity is of fundamental importance to control theory since it determines the asymptotic convergence rate $e^{\lambda(\mathcal{A},p)}$ of the stochastic linear dynamical system $x_{k+1} = A_{\sigma_k} x_k$. This paper investigates the following "design problem": given $\mathcal{A}$, compute the distribution $p$ minimizing $\lambda(\mathcal{A},p)$. Our main result is that it is NP-hard to decide whether there exists a distribution $p$ for which $\lambda(\mathcal{A},p)< 0$, i.e. it is NP-hard to decide whether this dynamical system can be stabilized. This hardness result holds even in the "simple"' case where $\mathcal{A}$ contains only rank-one matrices. Somewhat surprisingly, this is in stark contrast to the Joint Spectral Radius -- the deterministic kindred of the Lyapunov exponent -- for which the analogous optimization problem for rank-one matrices is known to be exactly computable in polynomial time. To prove this hardness result, we first observe via Birkhoff's Ergodic Theorem that the Lyapunov exponent of rank-one matrices admits a simple formula and in fact is a quadratic form in $p$. Hardness of the design problem is shown through a reduction from the Independent Set problem. Along the way, simple examples are given illustrating that $p \mapsto \lambda(\mathcal{A},p)$ is neither convex nor concave in general. We conclude with extensions to continuous distributions, exchangeable processes, Markov processes, and stationary ergodic processes.


title: A Structure-Tensor Approach to Integer Matrix Completion with Applications to Differentiated Energy Services

abstract: Efficient resource allocation is one of the main driving forces of human civilizations. Of the many existing approaches to resource allocation, matrix completion is one that is frequently applied. In this paper, we investigate a special type of matrix completion problem concerning the class of $(0,1)$-matrices with given row/column sums and certain zeros prespecified. We provide a necessary and sufficient condition under which such a class is nonempty. The condition is stated in the form of the nonnegativity of a structure tensor constructed from the information regarding the given row/column sums and fixed zeros. Moreover, we show that a more general matrix completion problem can be studied in a similar manner, namely that involving the class of nonnegative integer matrices with prescribed row/column sums, predetermined zeros, and different bounds across the rows. To illustrate the utility of our results, we apply them to demand response applications in smart grids. Specifically, we address two adequacy problems in differentiated energy services, namely, the problems of supply/demand matching and minimum purchase profile.


title: Massively Parallel Computation via Remote Memory Access

abstract: We introduce the Adaptive Massively Parallel Computation (AMPC) model, which is an extension of the Massively Parallel Computation (MPC) model. At a high level, the AMPC model strengthens the MPC model by storing all messages sent within a round in a distributed data store. In the following round, all machines are provided with random read access to the data store, subject to the same constraints on the total amount of communication as in the MPC model. Our model is inspired by the previous empirical studies of distributed graph algorithms using MapReduce and a distributed hash table service. This extension allows us to give new graph algorithms with much lower round complexities compared to the best known solutions in the MPC model. In particular, in the AMPC model we show how to solve maximal independent set in $O(1)$ rounds and connectivity/minimum spanning tree in $O(\log\log_{m/n} n)$ rounds both using $O(n^\delta)$ space per machine for constant $\delta < 1$. In the same memory regime for MPC, the best known algorithms for these problems require polylog $n$ rounds. Our results imply that the 2-Cycle conjecture, which is widely believed to hold in the MPC model, does not hold in the AMPC model.


title: Adiabatic and Nonadiabatic Energy Dissipation during Scattering of Vibrationally Excited CO from Au(111)

abstract: A high-dimensional potential energy surface (PES) for CO interaction with the Au(111) surface is developed using a machine-learning algorithm. Including both molecular and surface coordinates, this PES enables the simulation of the recent experiment on scattering of vibrationally excited CO from Au(111). Trapping in a physisorption well is observed to increase with decreasing incidence energy. While energy dissipation of physisorbed CO is slow, due to weak coupling with both the phonons and electron-hole pairs, its access to the chemisorption well facilitates fast vibrational relaxation of CO through nonadiabatic coupling with surface electron-hole pairs.


title: Perfect 1-factorisations of $K_{16}$

abstract: We report the results of a computer enumeration that found that there are 3155 perfect 1-factorisations (P1Fs) of the complete graph $K_{16}$. Of these, 89 have a non-trivial automorphism group (correcting an earlier claim of 88 by Meszka and Rosa). We also (i) describe a new invariant which distinguishes between the P1Fs of $K_{16}$, (ii) observe that the new P1Fs produce no atomic Latin squares of order 15 and (iii) record P1Fs for a number of large orders that exceed prime powers by one.


title: A Fixed-Flexible BESS Allocation Scheme for Transmission Networks Considering Uncertainties

abstract: Battery energy storage systems (BESSs) can play a key role in mitigating the intermittency and uncertainty associated with adding large amounts of renewable energy to the bulk power system (BPS). Two BESS technologies that have gained prominence in this regard are Lithium-ion (LI) BESS and Vanadium redox flow (VRF) BESS. This paper proposes a fixed-flexible BESS allocation scheme that exploits the complementary characteristics of LI and VRF BESSs to attain optimal techno-economic benefits in a wind-integrated BPS. Studies carried out on relatively large transmission networks demonstrate that benefits such as reduction in system operation cost, wind spillage, voltage fluctuations, and discounted payback period, can be realized by using the proposed scheme.


title: Simultaneous Generation of Arbitrary Assembly of Polarization States with Geometrical-Scaling-Induced Phase Modulation

abstract: Manipulating the polarization of light on the microscale or nanoscale is essential for integrated photonics and quantum optical devices. Nowadays, the metasurface allows one to build on-chip devices that efficiently manipulate the polarization states. However, it remains challenging to generate different types of polarization states simultaneously, which is required to encode information for quantum computing and quantum cryptography applications. By introducing geometrical-scaling-induced (GSI) phase modulations, we demonstrate that an assembly of circularly polarized (CP) and linearly polarized (LP) states can be simultaneously generated by a single metasurface made of L-shaped resonators with different geometrical features. Upon illumination, each resonator diffracts the CP state with a certain GSI phase. The interaction of these diffractions leads to the desired output beams, where the polarization state and the propagation direction can be accurately tuned by selecting the geometrical shape, size, and spatial sequence of each resonator in the unit cell. As an example of potential applications, we show that an image can be encoded with different polarization profiles at different diffraction orders and decoded with a polarization analyzer. This approach resolves a challenging problem in integrated optics and is inspiring for on-chip quantum information processing.


title: Frame measures for infinitely many measures

abstract: For every frame spectral measure $ \mu $, there exists a discrete measure $ \nu $ as a frame measure. Since if $ \mu $ is not a frame spectral measure, then there is not any general statement about the existence of frame measures $ \nu $ for $ \mu $, we were motivated to examine Bessel and frame measures. We construct infinitely many measures $ \mu $ which admit frame measures $ \nu $, and we show that there exist infinitely many frame spectral measures $ \mu $ such that besides having a discrete frame measure, they admit continuous frame measures too.


title: Forecasts for warm dark matter from weakly lensed Type Ia supernovae

abstract: We investigate the possibility to have a constraint on the mass of thermal warm dark matter (WDM) particle from the expected data of the Wide Field Infrared Survey Telescope (WFIRST) survey if all the dark matter is warm. For this purpose we consider the lensing effect of large scale structure based on the warm dark matter scenario on the apparent magnitude of SNe Ia. We use HALOFIT as non-linear matter power spectrum and the variance of PDF. We preform a Fisher matrix analysis and obtain the lower bound of $m_{\rm WDM}>0.167$keV.


title: Practical Bayesian Optimization with Threshold-Guided Marginal Likelihood Maximization

abstract: We propose a practical Bayesian optimization method using Gaussian process regression, of which the marginal likelihood is maximized where the number of model selection steps is guided by a pre-defined threshold. Since Bayesian optimization consumes a large portion of its execution time in finding the optimal free parameters for Gaussian process regression, our simple, but straightforward method is able to mitigate the time complexity and speed up the overall Bayesian optimization procedure. Finally, the experimental results show that our method is effective to reduce the execution time in most of cases, with less loss of optimization quality.


title: Rejuvenation in $z\sim0.8$ quiescent galaxies in LEGA-C

abstract: We use reconstructed star-formation histories (SFHs) of quiescent galaxies at $z=0.6-1$ in the LEGA-C survey to identify secondary star-formation episodes that, after an initial period of quiescence, moved the galaxies back to the star-forming main sequence (blue cloud). $16\pm3$\% of the $z\sim0.8$ quiescent population has experienced such rejuvenation events in the redshift range $0.7<z<1.5$ after reaching quiescence at some earlier time. On average, these galaxies first became quiescent at $z=1.2$, and those that rejuvenated, remained quiescent for $\sim1$Gyr before their secondary SF episode which lasted $\sim0.7$Gyr. The stellar mass attributed to rejuvenation is on average 10\% of the galaxy stellar mass, with rare instances of an increase of more than a factor 2. Overall, rejuvenation events only contribute $\sim2$\% of the total stellar mass in $z\sim0.8$ quiescent galaxies and we conclude that rejuvenation is not an important evolutionary channel when considering the growth of the red sequence. However, our results complicate the interpretation of galaxy demographics in color space: the galaxies with rejuvenation events tend to lie in the so-called `green valley', yet their progenitors were quiescent at $z\sim2$.


title: Semi-Supervised Monocular Depth Estimation with Left-Right Consistency Using Deep Neural Network

abstract: There has been tremendous research progress in estimating the depth of a scene from a monocular camera image. Existing methods for single-image depth prediction are exclusively based on deep neural networks, and their training can be unsupervised using stereo image pairs, supervised using LiDAR point clouds, or semi-supervised using both stereo and LiDAR. In general, semi-supervised training is preferred as it does not suffer from the weaknesses of either supervised training, resulting from the difference in the cameras and the LiDARs field of view, or unsupervised training, resulting from the poor depth accuracy that can be recovered from a stereo pair. In this paper, we present our research in single image depth prediction using semi-supervised training that outperforms the state-of-the-art. We achieve this through a loss function that explicitly exploits left-right consistency in a stereo reconstruction, which has not been adopted in previous semi-supervised training. In addition, we describe the correct use of ground truth depth derived from LiDAR that can significantly reduce prediction error. The performance of our depth prediction model is evaluated on popular datasets, and the importance of each aspect of our semi-supervised training approach is demonstrated through experimental results. Our deep neural network model has been made publicly available.


title: High-SNR snapshot multiplex spectrometer with sub-Hadamard-S matrix coding

abstract: We present a robust high signal-to-noise ratio (SNR) snapshot multiplex spectrometer with sub-Hadamard-S matrix coding. We demonstrated for the first time that the sub-Hadamard-S matrix coding could provide comparable SNR improvement with Hadamard-S matrix in Hadamard transform spectrometer (HTS). Normally, HTS should change the coding mask to obtain a reasonable spectrum result, causing unexpected time-consuming. An extra imaging path to collect the light intensity of the aperture is added in this paper. Both light intensity of the aperture and overlapped spectra are captured within one shot, turning Hadamard-S matrix coding into sub-Hadamard-S matrix coding. Simulations and experiments show that the proposed method could obtain comparable SNR improvement with the traditional HTS, maintaining snapshot.


title: Driver Surge Pricing

abstract: Ride-hailing marketplaces like Uber and Lyft use dynamic pricing, often called surge, to balance the supply of available drivers with the demand for rides. We study driver-side payment mechanisms for such marketplaces, presenting the theoretical foundation that has informed the design of Uber's new additive driver surge mechanism. We present a dynamic stochastic model to capture the impact of surge pricing on driver earnings and their strategies to maximize such earnings. In this setting, some time periods (surge) are more valuable than others (non-surge), and so trips of different time lengths vary in the induced driver opportunity cost. First, we show that multiplicative surge, historically the standard on ride-hailing platforms, is not incentive compatible in a dynamic setting. We then propose a structured, incentive-compatible pricing mechanism. This closed-form mechanism has a simple form and is well-approximated by Uber's new additive surge mechanism. Finally, through both numerical analysis and real data from a ride-hailing marketplace, we show that additive surge is more incentive compatible in practice than is multiplicative surge.


title: A sharp $L_p$-regularity result for second-order stochastic partial differential equations with unbounded and fully degenerate leading coefficients

abstract: We present existence, uniqueness, and sharp regularity results of solution to the stochastic partial differential equation (SPDE) \begin{align} \label{abs eqn} du=(a^{ij}(\omega,t)u_{x^ix^j}+f)dt + (\sigma^{ik}(\omega,t)u_{x^i}+g^k)dw^k_t, \quad u(0,x)=u_0, \end{align} where $\{w^k_t:k=1,2,\cdots\}$ is a sequence of independent Brownian motions. The coefficients are merely measurable in $(\omega,t)$ and can be unbounded and fully degenerate, that is, coefficients $a^{ij}$, $\sigma^{ik}$ merely satisfy \begin{align} \label{abs only} \left(\alpha^{ij}(\omega,t)\right)_{d\times d}:= \left(a^{ij}(\omega,t)-\frac{1}{2}\sum_{k=1}^{\infty} \sigma^{ik}(\omega,t)\sigma^{jk}(\omega,t)\right) \geq 0. \end{align} In this article, we prove that there exists a unique solution $u$ to \eqref{abs eqn}, and \begin{align} \notag \|u_{xx}\|_{\mathbb{H}^\gamma_p(\tau,\delta)} &\leq N(d,p) \bigg( \|u_0\|_{\mathbb{B}_p^{\gamma+2 \left(1-1/ p \right)}} + \| f\|_{\mathbb{H}^\gamma_p( \tau,\delta^{1-p} )} \label{abs est} &\qquad \qquad+\|g_x\|^p_{\mathbb{H}^\gamma_p( \tau, |\sigma|^p \delta^{1-p},l_2)}+ \| g_x\|_{\mathbb{H}^\gamma_p( \tau,\delta^{1-p/2},l_2)} \bigg), \end{align} where $p\geq 2$, $\gamma\in \mathbf{R}$, $\tau$ is an arbitrary stopping time, $\delta(\omega, t)$ is the smallest eigenvalue of $\alpha^{ij}(\omega, t)$, $\mathbb{H}_p^\gamma(\tau, \delta)$ is a weighted stochastic Sobolev space, and $\mathbb{B}_p^{\gamma+2 \left(1-1/ p \right)}$ is a stochastic Besov space.


title: Hedging crop yields against weather uncertainties -- a weather derivative perspective

abstract: The effects of weather on agriculture in recent years have become a major global concern. Hence, the need for an effective weather risk management tool (i.e., weather derivatives) that can hedge crop yields against weather uncertainties. However, most smallholder farmers and agricultural stakeholders are unwilling to pay for the price of weather derivatives (WD) because of the presence of basis risks (product-design and geographical) in the pricing models. To eliminate product-design basis risks, a machine learning ensemble technique was used to determine the relationship between maize yield and weather variables. The results revealed that the most significant weather variable that affected the yield of maize was average temperature. A mean-reverting model with a time-varying speed of mean reversion, seasonal mean, and local volatility that depended on the local average temperature was then proposed. The model was extended to a multi-dimensional model for different but correlated locations. Based on these average temperature models, pricing models for futures, options on futures, and basket futures for cumulative average temperature and growing degree-days are presented. Pricing futures on baskets reduces geographical basis risk, as buyers have the opportunity to select the most appropriate weather stations with their desired weight preference. With these pricing models, farmers and agricultural stakeholders can hedge their crops against the perils of extreme weather.


title: Kantorovich distance on a finite metric space

abstract: Kantorovich distance (or 1-Wasserstein distance) on the probability simplex of a finite metric space is the value of a Linear Programming problem for which a closed-form expression is known in some cases. When the ground distance is defined by a graph, a few examples have already been studied. In the present paper, after re-deriving, with different tools, the result for trees, we prove that, for an arbitrary weighted graph, the K-distance is the minimum of the K-distances over all the spanning trees associated with the graph. We work in the dual LP-problem by using Arens-Eells norm associated with the metric space. Finally, we introduce new norms that are naturally related to $\ell_1$-embeddable distances and allows for a partial extension of our results to this new setting.


title: Topological Phase Transition Independent of System Non-Hermiticity

abstract: Non-Hermiticity can vary the topology of system, induce topological phase transition, and even invalidate the conventional bulk-boundary correspondence. Here, we show the introducing of non-Hermiticity without affecting the topological properties of the original chiral symmetric Hermitian systems. Conventional bulk-boundary correspondence holds, topological phase transition and the (non)existence of edge states are unchanged even though the energy bands are inseparable due to non-Hermitian phase transition. Chern number for energy bands of the generalized non-Hermitian system in two dimension is proved to be unchanged and favorably coincides with the simulated topological charge pumping. Our findings provide insights into the interplay between non-Hermiticity and topology. Topological phase transition independent of non-Hermitian phase transition is a unique feature that beneficial for future applications of non-Hermitian topological materials.

categories: cond-mat.mes-hall cond-mat.other cond-mat.quant-gas cond-mat.str-el quant-ph


title: Multi-Armed Bandits for Boolean Connectives in Hybrid System Falsification (Extended Version)

abstract: Hybrid system falsification is an actively studied topic, as a scalable quality assurance methodology for real-world cyber-physical systems. In falsification, one employs stochastic hill-climbing optimization to quickly find a counterexample input to a black-box system model. Quantitative robust semantics is the technical key that enables use of such optimization. In this paper, we tackle the so-called scale problem regarding Boolean connectives that is widely recognized in the community: quantities of different scales (such as speed [km/h] vs. RPM, or worse, RPH) can mask each other's contribution to robustness. Our solution consists of integration of the multi-armed bandit algorithms in hill climbing-guided falsification frameworks, with a technical novelty of a new reward notion that we call hill-climbing gain. Our experiments show our approach's robustness under the change of scales, and that it outperforms a state-of-the-art falsification tool.


title: Strong equivalence for $\rm LP^{MLN}$ programs

abstract: Strong equivalence is a well-studied and important concept in answer set programming (ASP). $\rm LP^{MLN}$ is a probabilistic extension of answer set programs with the weight scheme adapted from Markov Logic. Because of the semantic differences, strong equivalence for ASP does not simply carry over to $\rm LP^{MLN}$. I study the concept of strong equivalence in $\rm LP^{MLN}$ with the goal of extending strong equivalence to $\rm LP^{MLN}$ programs. My study shows that the verification of strong equivalence in $\rm LP^{MLN}$ can be reduced to equivalence checking in classical logic plus weight consideration.The result allows us to leverage an answer set solver for checking strong equivalence in $\rm LP^{MLN}$. Furthermore, this study also suggests us a few reformulations of the $\rm LP^{MLN}$ semantics using choice rules, logic of here and there, and the second-order logic. I will present my work result of strong equivalence for $\rm LP^{MLN}$ and talk about my next steps for research: one is approximately strong equivalence, and another is the integration of fuzzy logic with neural network.


title: Ordinal Sums of Fuzzy Negations: Main Classes and Natural Negations

abstract: In the context of fuzzy logic, ordinal sums provide a method for constructing new functions from existing functions, which can be triangular norms, triangular conorms, fuzzy negations, copulas, overlaps, uninorms, fuzzy implications, among others. As our main contribution, we establish conditions for the ordinal sum of a family of fuzzy negations to be a fuzzy negation of a specific class, such as strong, strict, continuous, invertible and frontier. Also, we relate the natural negation of the ordinal sum on families of t-norms, t-conorms and fuzzy implications with the ordinal sum of the natural negations of the respective families of t-norms, t- conorms and fuzzy implications. This motivated us to introduces a new kind of ordinal sum for families of fuzzy implications.


title: Phase noise accumulation in recirculating frequency shifting loop based programmable optical frequency comb

abstract: The phenomenon of linewidth continuously broadening along with the recirculation number in recirculating frequency shifting loop is observed. In this paper, a novel method of measuring the phase noise accumulation induced by EDFA in recirculating frequency shifting loop is proposed and the experiment results support the viewpoint of the laser linewidth will be broadening by EDFA. An empirical formula has been extracted to estimate the linewidth deterioration of the RFS output. We demonstrate the relationship between the linewidth and the recirculation number of the RFS based optical frequency comb in both theoretically and experimentally. By employing the recirculation frequency shifting loop, the phase noise accumulation induced by EDFA could be measured obviously and the linewidth of each tone can be measured precisely.