title: Multi-Particle Simulation Techniques

abstract: The nonlinear space-charge effects play an important role in high intensity/high brightness accelerators. These effects can be self-consistently studied using multi-particle simulations. In this lecture, we will discuss the particle-in-cell method and the symplectic tracking model for self-consistent multi-particle simulations.


title: Dynamical response and competing orders in two-band Hubbard model

abstract: We present a dynamical mean-field study of two-particle dynamical response functions in two-band Hubbard model across several phase transitions. We observe that the transition between theexcitonic condensate and spin-state ordered state is continuous with a narrow strip of supersolidphase separating the two. Approaching transition from the excitonic condensate is announced bysoftening of the excitonic mode at theMpoint of the Brillouin zone. Inside the spin-state orderedphase there is a magnetically ordered state with 2x2 periodicity, which has no precursor in thenormal phase.


title: Explicit formulas for the joint third and fourth central moments of the multinomial distribution

abstract: We give the first explicit formulas for the joint third and fourth central moments of the multinomial distribution, by differentiating the moment generating function. A general formula for the joint factorial moments was previously given in Mosimann (1962).


title: Remarks on Ornstein's non-inequality in $\mathbb{R}^{2\times 2}$

abstract: We give a very concise proof of Ornstein's $L^1$ non-inequality for first- and second-order operators in two dimensions. The proof just needs a two-dimensional laminate supported on three points.


title: Bayesian Approximations to Hidden Semi-Markov Models for Telemetric Monitoring of Physical Activity

abstract: We propose a Bayesian hidden Markov model for analyzing time series and sequential data where a special structure of the transition probability matrix is embedded to model explicit-duration semi-Markovian dynamics. Our formulation allows for the development of highly flexible and interpretable models that can integrate available prior information on state durations while keeping a moderate computational cost to perform efficient posterior inference. We show the benefits of choosing a Bayesian approach for HSMM estimation over its frequentist counterpart, in terms of model selection and out-of-sample forecasting, also highlighting the computational feasibility of our inference procedure whilst incurring negligible statistical error. The use of our methodology is illustrated in an application relevant to e-Health, where we investigate rest-activity rhythms using telemetric activity data collected via a wearable sensing device. This analysis considers for the first time Bayesian model selection for the form of the explicit state dwell distribution. We further investigate the inclusion of a circadian covariate into the emission density and estimate this in a data-driven manner.


title: Thermodynamics of heavy quarkonium in magnetic field background

abstract: We study the effect of magnetic field on heavy quark-antiquark pair in both Einstein-Maxwell(EM) and Einstein-Maxwell-Dilaton(EMD) model. The interquark distance, free energy, entropy, binding energy and internal energy of the heavy quarkonium are calculated. It is found that the free energy suppresses and the entropy increases quickly with the increase of the magnetic field $B$. The binding energy vanishes at smaller distance when increasing the magnetic field, which indicates the quark-antiquark pair dissociates at smaller distance. The internal energy which consists of free energy and entropy will increase at large separating distance for non-vanishing magnetic field. These conclusions are consistent both in the EM and EMD model. Moreover, we also find that the quarkonium will dissociate easier in the parallel direction than that in the transverse direction for EMD model, but the conclusion is opposite in EM model. Lattice results are in favor of EMD model. Besides, a Coulomb-plus-linear potential(Cornell potential) can be realized only in EMD model. Thus, a dilaton field is proved to be important in holographic model. Finally, we also show that the free energy, entropy and internal energy of a single quark in EMD model with the presence of magnetic field.


title: Numerical evidence towards a positive answer to Morrey's problem

abstract: We report on numerical experiments suggesting that rank-one convexity imples quasiconvexity in the planar case. We give a simple heuristic explanation of our findings.


title: Entanglement marginal problems

abstract: We consider the entanglement marginal problem, which consists of deciding whether a number of reduced density matrices are compatible with an overall separable quantum state. To tackle this problem, we propose hierarchies of semidefinite programming relaxations of the set of quantum state marginals admitting a fully separable extension. We connect the completeness of each hierarchy to the resolution of an analog classical marginal problem and thus identify relevant experimental situations where the hierarchies are complete. For finitely many parties on a star configuration or a chain, we find that we can achieve an arbitrarily good approximation to the set of nearest-neighbour marginals of separable states with a time (space) complexity polynomial (linear) on the system size. Our results even extend to infinite systems, such as translation-invariant systems in 1D, as well as higher spatial dimensions with extra symmetries.


title: The limits of min-max optimization algorithms: convergence to spurious non-critical sets

abstract: Compared to ordinary function minimization problems, min-max optimization algorithms encounter far greater challenges because of the existence of periodic cycles and similar phenomena. Even though some of these behaviors can be overcome in the convex-concave regime, the general case is considerably more difficult. On that account, we take an in-depth look at a comprehensive class of state-of-the art algorithms and prevalent heuristics in non-convex / non-concave problems, and we establish the following general results: a) generically, the algorithms' limit points are contained in the ICT sets of a common, mean-field system; b) the attractors of this system also attract the algorithms in question with arbitrarily high probability; and c) all algorithms avoid the system's unstable sets with probability 1. On the surface, this provides a highly optimistic outlook for min-max algorithms; however, we show that there exist spurious attractors that do not contain any stationary points of the problem under study. In this regard, our work suggests that existing min-max algorithms may be subject to inescapable convergence failures. We complement our theoretical analysis by illustrating such attractors in simple, two-dimensional, almost bilinear problems.


title: Multi-Agent Reinforcement Learning for Adaptive User Association in Dynamic mmWave Networks

abstract: Network densification and millimeter-wave technologies are key enablers to fulfill the capacity and data rate requirements of the fifth generation (5G) of mobile networks. In this context, designing low-complexity policies with local observations, yet able to adapt the user association with respect to the global network state and to the network dynamics is a challenge. In fact, the frameworks proposed in literature require continuous access to global network information and to recompute the association when the radio environment changes. With the complexity associated to such an approach, these solutions are not well suited to dense 5G networks. In this paper, we address this issue by designing a scalable and flexible algorithm for user association based on multi-agent reinforcement learning. In this approach, users act as independent agents that, based on their local observations only, learn to autonomously coordinate their actions in order to optimize the network sum-rate. Since there is no direct information exchange among the agents, we also limit the signaling overhead. Simulation results show that the proposed algorithm is able to adapt to (fast) changes of radio environment, thus providing large sum-rate gain in comparison to state-of-the-art solutions.


title: A specifically designed machine learning algorithm for GNSS position time series prediction and its applications in outlier and anomaly detection and earthquake prediction

abstract: We present a simple yet efficient supervised machine learning algorithm that is designed for the GNSS position time series prediction. This algorithm has four steps. First, the mean value of the time series is subtracted from it. Second, the trends in the time series are removed. Third, wavelets are used to separate the high and low frequencies. And fourth, a number of frequencies are derived and used for finding the weights between the hidden and the output layers, using the product of the identity and sine and cosine functions. The role of the observation precision is taken into account in this algorithm. A large-scale study of three thousand position times series of GNSS stations across the globe is presented. Seventeen different machine learning algorithms are examined. The accuracy levels of these algorithms are checked against the rigorous statistical method of Theta. It is shown that the most accurate machine learning algorithm is the method we present, in addition to being faster. Two applications of the algorithm are presented. In the first application, it is shown that the outliers and anomalies in a time series can be detected and removed by the proposed algorithm. In a large scale study, ten other methods of time series outlier detection are compared with the proposed algorithm. The study reveals that the proposed algorithm is approximately 3.22 percent more accurate in detecting outliers. In the second application, the suitability of the algorithm for earthquake prediction is investigated. A case study is presented for the Tohoku 2011 earthquake. It is shown that this earthquake could have been predicted approximately 2 hours before its happening, solely based on each of the 845 GEONET station time series. Comparison with four different studies show the improvement in prediction of the time of the earthquake.


title: Formation of aqua planets with water of nebular origin: Effects of water enrichment on the structure and mass of captured atmospheres of terrestrial planets

abstract: Recent detection of exoplanets with Earth-like insolation attracts growing interest in how common Earth-like aqua planets are beyond the solar system. While terrestrial planets are often assumed to capture icy or water-rich planetesimals, a primordial atmosphere of nebular origin itself can produce water through oxidation of the atmospheric hydrogen with oxidising minerals from incoming planetesimals or the magma ocean. Thermodynamically, normal oxygen buffers produce water comparable in mole number to or more than hydrogen. Thus, the primordial atmosphere would likely be highly enriched with water vapour; however, the primordial atmosphereshave been always assumed to have the solar abundances. Here we integrate the 1D structure of such an enriched atmosphere of sub-Earths embedded in a protoplanetary disc around an M dwarf of 0.3$M_\odot$ and investigate the effects of water enrichment on the atmospheric properties with focus on water amount. We find that the well-mixed, highly-enriched atmosphere is more massive by a few orders of magnitude than the solar-abundance atmosphere, and that even a Mars-mass planet can obtain water comparable to the present Earth's oceans. Although close-in Mars-mass planets likely lose the captured water via disc dispersal and photo-evaporation, these results suggest that there are more sub-Earths with Earth-like water contents than previously predicted. How much water terrestrial planets really obtain and retain against subsequent loss, however, depends on efficiencies of water production, mixing in the atmosphere and magma ocean, and photo-evaporation, detailed investigation for which should be made in the future.


title: Entanglement and symmetry resolution in two dimensional free quantum field theories

abstract: We present a thorough analysis of the entanglement entropies related to different symmetry sectors of free quantum field theories (QFT) with an internal U(1) symmetry. We provide explicit analytic computations for the charged moments of Dirac and complex scalar fields in two spacetime dimensions, both in the massive and massless cases, using two different approaches. The first one is based on the replica trick, the computation of the partition function on Riemann surfaces with the insertion of a flux $\alpha$, and the introduction of properly modified twist fields, whose two-point function directly gives the scaling limit of the charged moments. With the second method, the diagonalisation in replica space maps the problem to the computation of a partition function on a cut plane, that can be written exactly in terms of the solutions of non-linear differential equations of the Painlev\'e V type. Within this approach, we also derive an asymptotic expansion for the short and long distance behaviour of the charged moments. Finally, the Fourier transform provides the desired symmetry resolved entropies: at the leading order, they satisfy entanglement equipartition and we identify the subleading terms that break it. Our analytical findings are tested against exact numerical calculations in lattice models.


title: Topological Dirac states in a layered telluride TaPdTe$_5$ with quasi-one-dimensional PdTe$_2$ chains

abstract: We report the synthesis and systematic studies of a new layered ternary telluride TaPdTe5 with quasi-one-dimensional PdTe2 chains. This compound crystalizes in a layered orthorhombic structure with space group Cmcm. Analysis of its curved field-dependent Hall resistivity, using the two-band model, indicates the hole-dominated transport with a high mobility ${\mu}_h$ = 2.38 $\times$ 10$^3$ cm$^2$ V$^{-1}$ s$^{-1}$ at low temperatures. The in-plane magnetoresistance (MR) displays significant anisotropy with field applied along the crystallographic $b$ axis. The MR with the current applied along the $c$-axis is also measured in high magnetic fields up to 51.7 T. Remarkably, it follows a power-law dependence and reaches (9.5 $\times$ 10$^3$)% at 2.1 K without any signature of saturation. The De Haas-van Alphen oscillations show a small Fermi-surface pocket with a nontrivial Berry phase. The Shubnikov-de Haas (SdH) oscillations are detected at low temperatures and under magnetic fields above 28.5 T. Two effective masses $m^*$ (0.26$m_e$ and 0.41$m_e$) are extracted from the oscillatory SdH data. Our first-principles calculations unveil a topological Dirac cone in its surface states, and, in particular, the topological index indicates that TaPdTe$_5$ is a topologically nontrivial material.


title: Multilinear commutators in the two-weight setting

abstract: We extend the recently much-studied two-weight commutator estimates to the multilinear setting. In contrast to previous results, our result respects the multilinear nature of the problem fully and is formulated with the genuinely multilinear weights.


title: Divergence-free measures in the plane and inverse potential problems in divergence form

abstract: We show that a divergence-free measure on the plane is a continuous sum of unit tangent vector fields on rectifiable Jordan curves. This loop decomposition is more precise than the general decomposition in elementary solenoids given by S.K. Smirnov, when applied to the planar case. The proof involves extending the Fleming-Rishel formula to homogeneous BV functions (in any dimension), and establishing for such functions approximate continuity of measure theoretic connected components of suplevel sets as functions of the level. We apply these results to inverse potential problems whose source term is the divergence of some unknown (vector-valued) measure. A prototypical case is that of inverse magnetization problems when magnetizations are modeled by R3-valued Borel measures. We investigate methods for recovering a magnetization {\mu} by penalizing its measure theoretic total variation norm (TV). In particular, we prove that if a magnetization is supported in a plane, then TV-regularization schemes always have a unique minimizer, even in the presence of noise. It is further shown thatTV-norm minimization (among magnetizations generating the same field) uniquely recovers planar magnetizations in the following cases: when the magnetization is carried by a collection of sufficiently separated line segments and a set that is purely 1-unrectifiable, or when a superset of the support is tree-like. We note that such magnetizations can be recovered via TV-regularization schemes in the zero noise limit, by taking the regularization parameter to zero. This suggests definitions of sparsity in the present infinite dimensional context, that generate results akin to compressed sensing


title: Mucko: Multi-Layer Cross-Modal Knowledge Reasoning for Fact-based Visual Question Answering

abstract: Fact-based Visual Question Answering (FVQA) requires external knowledge beyond visible content to answer questions about an image, which is challenging but indispensable to achieve general VQA. One limitation of existing FVQA solutions is that they jointly embed all kinds of information without fine-grained selection, which introduces unexpected noises for reasoning the final answer. How to capture the question-oriented and information-complementary evidence remains a key challenge to solve the problem. In this paper, we depict an image by a multi-modal heterogeneous graph, which contains multiple layers of information corresponding to the visual, semantic and factual features. On top of the multi-layer graph representations, we propose a modality-aware heterogeneous graph convolutional network to capture evidence from different layers that is most relevant to the given question. Specifically, the intra-modal graph convolution selects evidence from each modality and cross-modal graph convolution aggregates relevant information across different modalities. By stacking this process multiple times, our model performs iterative reasoning and predicts the optimal answer by analyzing all question-oriented evidence. We achieve a new state-of-the-art performance on the FVQA task and demonstrate the effectiveness and interpretability of our model with extensive experiments.


title: Quantitative Group Testing and the rank of random matrices

abstract: Given a random Bernoulli matrix $ A\in \{0,1\}^{m\times n} $, an integer $ 0< k < n $ and the vector $ y:=Ax $, where $ x \in \{0,1\}^n $ is of Hamming weight $ k $, the objective in the {\em Quantitative Group Testing} (QGT) problem is to recover $ x $. This problem is more difficult the smaller $m$ is. For parameter ranges of interest to us, known polynomial time algorithms require values of $m$ that are much larger than $k$. In this work, we define a seemingly easier problem that we refer to as {\em Subset Select}. Given the same input as in QGT, the objective in Subset Select is to return a subset $ S \subseteq [n] $ of cardinality $ m $, such that for all $ i\in [n] $, if $ x_i = 1 $ then $ i\in S $. We show that if the square submatrix of $A$ defined by the columns indexed by $S$ has nearly full rank, then from the solution of the Subset Select problem we can recover in polynomial-time the solution $x$ to the QGT problem. We conjecture that for every polynomial time Subset Select algorithm, the resulting output matrix will satisfy the desired rank condition. We prove the conjecture for some classes of algorithms. Using this reduction, we provide some examples of how to improve known QGT algorithms. Using theoretical analysis and simulations, we demonstrate that the modified algorithms solve the QGT problem for values of $ m $ that are smaller than those required for the original algorithms.


title: PERL: Pivot-based Domain Adaptation for Pre-trained Deep Contextualized Embedding Models

abstract: Pivot-based neural representation models have lead to significant progress in domain adaptation for NLP. However, previous works that follow this approach utilize only labeled data from the source domain and unlabeled data from the source and target domains, but neglect to incorporate massive unlabeled corpora that are not necessarily drawn from these domains. To alleviate this, we propose PERL: A representation learning model that extends contextualized word embedding models such as BERT with pivot-based fine-tuning. PERL outperforms strong baselines across 22 sentiment classification domain adaptation setups, improves in-domain model performance, yields effective reduced-size models and increases model stability.


title: Connectors

abstract: Recently, the author and Yamamoto invented a new proof of the duality for multiple zeta values. The technique is applicable in other series identities. In this article, we exhibit such proofs for some series identities.


title: AstroSat Observation of Non-Resonant Type-C QPOs in MAXI J1535-571

abstract: Galactic transient black hole candidate (BHC) MAXI J1535-571 was discovered on 2017 September 02 simultaneously by {\it MAXI}/GSC and {\it Swift}/BAT instruments. It has also been observed by India's first multi-wavelength astronomy-mission satellite {\it AstroSat}, during the rising phase of its 2017-18 outburst. We make both the spectral and the temporal analysis of the source during 2017 September 12-17 using data of {\it AstroSat}'s Large Area X-ray Proportional Counter (LAXPC) in the energy range of $3-40$~keV to infer the accretion flow properties of the source. Spectral analysis is done with the physical two-component advective flow (TCAF) solution-based {\it fits} file. From the nature of the variation of the TCAF model fitted physical flow parameters, we conclude and confirm that the source was in the intermediate spectral state during our analysis period. We observe sharp type-C quasi-periodic oscillations (QPOs) in the frequency range of $\sim 1.75-2.81$~Hz. For a better understanding of the nature and evolution of these type-C QPOs, a dynamic study of the power density spectra is done. We also investigate the origin of these QPOs from the shock oscillation model. We find that non-satisfaction of Rankine-Hugoniot conditions for non-dissipative shocks and not their resonance oscillations is the cause of the observed type-C QPOs.


title: Fluctuation-dissipation theorem and fundamental photon commutation relations in lossy nanostructures using quasinormal modes

abstract: We provide theory and formal insight on the Green function quantization method for absorptive and dispersive spatial-inhomogeneous media in the context of dielectric media. We show that a fundamental Green function identity, which appears, e.g., in the fundamental commutation relation of the electromagnetic fields, is also valid in the limit of non-absorbing media. We also demonstrate how the zero-point field fluctuations yields a non-vanishing surface term in configurations without absorption, when using a more formal procedure of the Green function quantization method. We then apply the presented method to a recently developed theory of photon quantization using quasinormal modes [Franke et al., Phys. Rev. Lett. 122, 213901 (2019)] for finite nanostructures embedded in a lossless background medium. We discuss the strict dielectric limit of the commutation relations of the quasinormal mode operators and present different methods to obtain them, connected to the radiative loss for non-absorptive but open resonators. We show exemplary calculations of a fully three-dimensional photonic crystal beam cavity, including the lossless limit, which supports a single quasinormal mode and discuss the limits of the commutation relation for vanishing damping (no material loss and no radiative loss).


title: First detection of ozone in the mid-infrared at Mars: implications for methane detection

abstract: The ExoMars Trace Gas Orbiter (TGO) was sent to Mars in March 2016 to search for trace gases diagnostic of active geological or biogenic processes. We report the first observation of the spectral features of Martian ozone (O3) in the mid-infrared range using the Atmospheric Chemistry Suite (ACS) Mid-InfaRed (MIR) channel, a cross-dispersion spectrometer operating in solar occultation mode with the finest spectral resolution of any remote sensing mission to Mars. Observations of ozone were made at high northern latitudes (>65N) prior to the onset of the 2018 global dust storm (Ls = 163-193). During this fast transition phase between summer and winter ozone distribution, the O3 volume mixing ratio observed is 100-200 ppbv near 20 km. These amounts are consistent with past observations made at the edge of the southern polar vortex in the ultraviolet range. The observed spectral signature of ozone at 3000-3060 cm-1 directly overlaps with the spectral range of the methane (CH4) nu3 vibration-rotation band, and it, along with a newly discovered CO2 band in the same region, may interfere with measurements of methane abundance.


title: Direct Vlasov solvers

abstract: In these proceedings we will describe the theory and practical steps required to build Vlasov solvers such as those commonly used to compute coherent instabilities in synchrotrons. Thanks to a Hamiltonian formalism, we will derive a compact and general form of the linearized Vlasov equation, written using Poisson brackets. This in turn will be the basis of a procedure to build Vlasov solvers, applied to the specific example of transverse instabilities arising from beam coupling impedance.


title: Progressive Skeletonization: Trimming more fat from a network at initialization

abstract: Recent studies have shown that skeletonization (pruning parameters) of networks \textit{at initialization} provides all the practical benefits of sparsity both at inference and training time, while only marginally degrading their performance. However, we observe that beyond a certain level of sparsity (approx $95\%$), these approaches fail to preserve the network performance, and to our surprise, in many cases perform even worse than trivial random pruning. To this end, we propose an objective to find a skeletonized network with maximum {\em foresight connection sensitivity} (FORCE) whereby the trainability, in terms of connection sensitivity, of a pruned network is taken into consideration. We then propose two approximate procedures to maximize our objective (1) Iterative SNIP: allows parameters that were unimportant at earlier stages of skeletonization to become important at later stages; and (2) FORCE: iterative process that allows exploration by allowing already pruned parameters to resurrect at later stages of skeletonization. Empirical analyses on a large suite of experiments show that our approach, while providing at least as good a performance as other recent approaches on moderate pruning levels, provides remarkably improved performance on higher pruning levels (could remove up to $99.5\%$ parameters while keeping the networks trainable). Code can be found in https://github.com/naver/force.


title: A New Comprehensive Data Set of Solar Filaments of 100 yr Interval. I

abstract: Filaments are very common physical phenomena on the Sun and are often taken as important proxies of solar magnetic activities. The study of filaments has become a hot topic in the space weather research. For a more comprehensive understanding of filaments, especially for an understanding of solar activities of multiple solar cycles, it is necessary to perform a combined multifeature analysis by constructing a data set of multiple solar cycle data. To achieve this goal, we constructed a centennial data set that covers the H$\alpha$ data from five observatories around the world. During the data set construction, we encountered varieties of problems, such as data fusion, accurate determination of the solar edge, classifying data by quality, dynamic threshold, and so on, which arose mainly due to multiple sources and a large time span of data. But fortunately, these problems were well solved. The data set includes seven types of data products and eight types of feature parameters with which we can implement the functions of data searching and statistical analyses. It has the characteristics of better continuity and highly complementary to space observation data, especially in the wavelengths not covered by space observations, and covers many solar cycles (including more than 60 yr of high-cadence data). We expect that this new comprehensive data set as well as the tools will help researchers to significantly speed up their search for features or events of interest, for either statistical or case study purposes, and possibly help them get a better and more comprehensive understanding of solar filament mechanisms.


title: Reusing Trained Layers of Convolutional Neural Networks to Shorten Hyperparameters Tuning Time

abstract: Hyperparameters tuning is a time-consuming approach, particularly when the architecture of the neural network is decided as part of this process. For instance, in convolutional neural networks (CNNs), the selection of the number and the characteristics of the hidden (convolutional) layers may be decided. This implies that the search process involves the training of all these candidate network architectures. This paper describes a proposal to reuse the weights of hidden (convolutional) layers among different trainings to shorten this process. The rationale is that if a set of convolutional layers have been trained to solve a given problem, the weights calculated in this training may be useful when a new convolutional layer is added to the network architecture. This idea has been tested using the CIFAR-10 dataset, testing different CNNs architectures with up to 3 convolutional layers and up to 3 fully connected layers. The experiments compare the training time and the validation loss when reusing and not reusing convolutional layers. They confirm that this strategy reduces the training time while it even increases the accuracy of the resulting neural network. This finding opens up the future possibility of integrating this strategy in existing AutoML methods with the purpose of reducing the total search time.


title: Stochastic Unit Commitment in Electricity-Gas Coupled Integrated Energy Systems based on Modified Progressive Hedging

abstract: The increasing number of gas-fired units has significantly intensified the coupling between power and gas networks. Traditionally, the nonlinearity and nonconvexity in gas flow equations, together with renewable-induced stochasticity, result in a computationally expensive model for unit commitment in electricity-gas coupled integrated energy systems (IES). To accelerate stochastic day-ahead scheduling, we applied and modified Progressive Hedging (PH), a heuristic approach that can be computed in parallel to yield scenario-independent unit commitment. By applying a termination and enumeration technique, the modified PH algorithm saves considerable computational time, especially when the unit production prices are similar for all generators, and when the scale of IES is large. Moreover, an adapted second-order cone relaxation (SOCR) is utilized to tackle the nonconvex gas flow equation. Case studies are performed on the IEEE 24-bus system/Belgium 20-node gas system and the IEEE 118-bus system/Belgium 20-node gas system. The computational efficiency when employing PH is 188 times that of commercial software, even outperforming Benders Decomposition. Meanwhile, the gap between the PH algorithm and the benchmark is less than 0.01% in both IES systems, which proves that the solution produced by PH reaches acceptable optimality in this stochastic UC problem.


title: MCRapper: Monte-Carlo Rademacher Averages for Poset Families and Approximate Pattern Mining

abstract: We present MCRapper, an algorithm for efficient computation of Monte-Carlo Empirical Rademacher Averages (MCERA) for families of functions exhibiting poset (e.g., lattice) structure, such as those that arise in many pattern mining tasks. The MCERA allows us to compute upper bounds to the maximum deviation of sample means from their expectations, thus it can be used to find both statistically-significant functions (i.e., patterns) when the available data is seen as a sample from an unknown distribution, and approximations of collections of high-expectation functions (e.g., frequent patterns) when the available data is a small sample from a large dataset. This feature is a strong improvement over previously proposed solutions that could only achieve one of the two. MCRapper uses upper bounds to the discrepancy of the functions to efficiently explore and prune the search space, a technique borrowed from pattern mining itself. To show the practical use of MCRapper, we employ it to develop an algorithm TFP-R for the task of True Frequent Pattern (TFP) mining. TFP-R gives guarantees on the probability of including any false positives (precision) and exhibits higher statistical power (recall) than existing methods offering the same guarantees. We evaluate MCRapper and TFP-R and show that they outperform the state-of-the-art for their respective tasks.


title: Growth of Eigenfunctions and R-limits on Graphs

abstract: A characterization of the essential spectrum $\sigma_{\text ess}$ of Schr\"odinger operators on infinite graphs is derived involving the concept of $\mathcal{R}$-limits. This concept, which was introduced previously for operators on $\mathbb{N}$ and $\mathbb{Z}^d$ as "right-limits", captures the behaviour of the operator at infinity. For graphs with sub-exponential growth rate we show that each point in $\sigma_{\text ess}(H)$ corresponds to a bounded generalized eigenfunction of a corresponding $\mathcal{R}$-limit of $H$. If, additionally, the graph is of uniform sub-exponential growth, also the converse inclusion holds.


title: In situ and Operando X-ray Imaging of Directed Energy Deposition Additive Manufacturing

abstract: The mechanical performance of Directed Energy Deposition Additive Manufactured (DED-AM) components can be highly material dependent. Through in situ and operando synchrotron X-ray imaging we capture the underlying phenomena controlling build quality of stainless steel (SS316) and titanium alloy (Ti6242 or Ti-6Al-2Sn-4Zr-2Mo). We reveal three mechanisms influencing the build efficiency of titanium alloys compared to stainless steel: blown powder sintering; reduced melt-pool wetting due to the sinter; and pore pushing in the melt-pool. The former two directly increase lack of fusion porosity, while the later causes end of track porosity. Each phenomenon influences the melt-pool characteristics, wetting of the substrate and hence build efficacy and undesirable microstructural feature formation. We demonstrate that porosity is related to powder characteristics, pool flow, and solidification front morphology. Our results clarify DED-AM process dynamics, illustrating why each alloy builds differently, facilitating the wider application of additive manufacturing to new materials.


title: Lie point symmetries of near-horizon geometry equation

abstract: All the Lie point symmetries of the near extremal horizon geometry equation, in the case of 4-dimensional Einstein vacuum spacetime with cosmological constant, are the diffeomorphisms of the space of the null generators of the horizon. This result is also generalised to the Maxwell-Einstein spacetime.


title: Redundancy of triangle groups in spherical CR representations

abstract: Falbel, Koseleff and Rouillier computed a large number of boundary unipotent CR representations of fundamental groups of non compact three-manifolds. Those representations are not always discrete. By experimentally computing their limit set, one can determine that those with fractal limit sets are discrete. Many of those discrete representations can be related to (3,3,n) complex hyperbolic triangle groups. By exact computations, we verify the existence of those triangle representations, which have boundary unipotent holonomy. We also show that many representations are redundant: for n fixed, all the (3,3,n) representations encountered are conjugate and only one among them is uniformizable.


title: Response of Vulnerable Road Users to Visual Information from Autonomous Vehicles in Shared Spaces

abstract: Completely unmanned autonomous vehicles have been anticipated for a while. Initially, these are expected to drive only under certain conditions on some roads, and advanced functionality is required to cope with the ever-increasing challenges of safety. To enhance the public's perception of road safety and trust in new vehicular technologies, we investigate in this paper the effect of several interaction paradigms with vulnerable road users by developing and applying algorithms for the automatic analysis of pedestrian body language. We assess behavioral patterns and determine the impact of the coexistence of AVs and other road users on general road safety in a shared space for VRUs and vehicles. Results showed that the implementation of visual communication cues for interacting with VRUs is not necessarily required for a shared space in which informal traffic rules apply.


title: Flatness is a False Friend

abstract: Hessian based measures of flatness, such as the trace, Frobenius and spectral norms, have been argued, used and shown to relate to generalisation. In this paper we demonstrate that for feed forward neural networks under the cross entropy loss, we would expect low loss solutions with large weights to have small Hessian based measures of flatness. This implies that solutions obtained using $L2$ regularisation should in principle be sharper than those without, despite generalising better. We show this to be true for logistic regression, multi-layer perceptrons, simple convolutional, pre-activated and wide residual networks on the MNIST and CIFAR-$100$ datasets. Furthermore, we show that for adaptive optimisation algorithms using iterate averaging, on the VGG-$16$ network and CIFAR-$100$ dataset, achieve superior generalisation to SGD but are $30 \times$ sharper. This theoretical finding, along with experimental results, raises serious questions about the validity of Hessian based sharpness measures in the discussion of generalisation. We further show that the Hessian rank can be bounded by the a constant times number of neurons multiplied by the number of classes, which in practice is often a small fraction of the network parameters. This explains the curious observation that many Hessian eigenvalues are either zero or very near zero which has been reported in the literature.


title: Learning Rates as a Function of Batch Size: A Random Matrix Theory Approach to Neural Network Training

abstract: We study the effect of mini-batching on the loss landscape of deep neural networks using spiked, field-dependent random matrix theory. We demonstrate that the magnitude of the extremal values of the batch Hessian are larger than those of the empirical Hessian. We also derive similar results for the Generalised Gauss-Newton matrix approximation of the Hessian. As a consequence of our theorems we derive an analytical expressions for the maximal learning rates as a function of batch size, informing practical training regimens for both stochastic gradient descent (linear scaling) and adaptive algorithms, such as Adam (square root scaling), for smooth, non-convex deep neural networks. Whilst the linear scaling for stochastic gradient descent has been derived under more restrictive conditions, which we generalise, the square root scaling rule for adaptive optimisers is, to our knowledge, completely novel. %For stochastic second-order methods and adaptive methods, we derive that the minimal damping coefficient is proportional to the ratio of the learning rate to batch size. We validate our claims on the VGG/WideResNet architectures on the CIFAR-$100$ and ImageNet datasets. Based on our investigations of the sub-sampled Hessian we develop a stochastic Lanczos quadrature based on the fly learning rate and momentum learner, which avoids the need for expensive multiple evaluations for these key hyper-parameters and shows good preliminary results on the Pre-Residual Architecure for CIFAR-$100$.


title: Characterization of Dielectric Materials by Sparse Signal Processing with Iterative Dictionary Updates

abstract: Estimating parameters and properties of various materials without causing damage to the material under test (MUT) is important in many applications. Thus, in this letter, we address this by wireless sensing. Here, the accuracy of the estimation depends on the accurate estimation of the properties of the reflected signal from the MUT (e.g., number of reflections, their amplitudes and time delays). For a layered MUT, there are multiple reflections and, due to the limited bandwidth at the receiver, these reflections superimpose each other. Since the number of reflections coming from the MUT is limited, we propose sparse signal processing (SSP) to decompose the reflected signal. In SSP, a so called dictionary is required to obtain a sparse representation of the signal. Here, instead of a fixed dictionary, a dictionary update technique is proposed to improve the estimation of the reflected signal. To validate the proposed method, a vector network analyzer (VNA) based measurement setup is used. It turns out that the estimated dielectric constants are in close agreement with the dielectric constants of the MUTs reported in literature. Further, the proposed approach outperforms the state-of-the-art model-based curve-fitting approach in thickness estimation.


title: Another approach to non-repetitive colorings of graphs of bounded degree

abstract: We propose a new proof technique that aims to be applied to the same problems as the Lov\'asz Local Lemma or the entropy-compression method. We present this approach in the context of non-repetitive colorings and we use it to improve upper-bounds relating different non-repetitive numbers to the maximal degree of a graph. It seems that there should be other interesting applications to the presented approach. In terms of upper-bound our approach seems to be as strong as entropy-compression, but the proofs are more elementary and shorter. The application we provide in this paper are upper bounds for graphs of maximal degree at most $\Delta$: a minor improvement on the upper-bound of the non-repetitive number, a $4.25\Delta +o(\Delta)$ upper-bound on the weak total non-repetitive number and a $ \Delta^2+\frac{3}{2^\frac{1}{3}}\Delta^{\frac{5}{3}}+ o(\Delta^{\frac{5}{3}})$ upper-bound on the total non-repetitive number of graphs. This last result implies the same upper-bound for the non-repetitive index of graphs, which improves the best known bound.


title: Fatou and Julia like sets II

abstract: This paper is a continuation of authors work: Fatou and Julia like sets,Ukranian J. Math., to appear/arXiv:2006.08308[math.CV](see [4]). Here, we introduce escaping like set and generalized escaping like set for a family of holomorphic functions on an arbitrary domain, and establish some distinctive properties of these sets. The connectedness of the Julia like set is also proved.


title: The Bures Metric for Generative Adversarial Networks

abstract: Generative Adversarial Networks (GANs) are performant generative methods yielding high-quality samples. However, under certain circumstances, the training of GANs can lead to mode collapse or mode dropping, i.e. the generative models not being able to sample from the entire probability distribution. To address this problem, we use the last layer of the discriminator as a feature map to study the distribution of the real and the fake data. During training, we propose to match the real batch diversity to the fake batch diversity by using the Bures distance between covariance matrices in feature space. The computation of the Bures distance can be conveniently done in either feature space or kernel space in terms of the covariance and kernel matrix respectively. We observe that diversity matching reduces mode collapse substantially and has a positive effect on the sample quality. On the practical side, a very simple training procedure, that does not require additional hyperparameter tuning, is proposed and assessed on several datasets.


title: Accelerated Alternating Minimization and Adaptability to Strong Convexity

abstract: In the first part of the paper we consider accelerated first order optimization method for convex functions with $L$-Lipschitz-continuous gradient, that is able to automatically adapts to problems which satisfies Polyak-{\L}ojasiewicz condition or which is strongly convex with the value of parameter equal to $\mu$. In these cases method poses linear convergence with factor $\frac{\mu}{L}$, if $\mu$ is unknown. If the problem is strongly convex and $\mu$ is known, than the method possesses linear convergence with the factor $\sqrt{\frac{\mu}{L}}$. If that are not the cases, the method converges with a rate $O(1/k^2)$. The second part contains generalization of the method to the problems, that allows alternating minimization and proofs of the same asymptotic convergence rates. Also it is considered the approach called Adaptive Catalyst, which allows to increase convergence rate up to $O(1/k^2)$ and also it is provided an experimental comparison of the approach with AGM, Alternating AGM, APDAGD and Sinkhorn's algorithm for the dual problem to the discrete entropically regularized optimal transport problem. The result of the work is the attempt to explain the reason why Alternating AGM converge faster than AGM or Adaptive Catalyst despite of the asymptotic theoretical rate $O(1/k^2)$. The hypothesis relies on the fact that Alternating AGM adapts to strong convexity. The hypothesis was tested on quadratic functions, on that Alternating AGM also showed faster convergence.


title: Periodic Hamiltonian systems in shape optimization problems with Neumann boundary conditions

abstract: The recent approach based on Hamiltonian systems and the implicit parametri\-za\-tion theorem, provides a general fixed domain approximation method in shape optimization problems, using optimal control theory. In previous works, we have examined Dirichlet boundary conditions with distributed or boundary observation. Here, we discuss the case of Neumann boundary conditions, with a combined cost functional, including both distributed and boundary observation. Extensions to nonlinear state systems are possible. This new technique allows simultaneous boundary and topological variations and we also report numerical experiments confirming the theoretical results.


title: Identifying the BLE Advertising Channel for Reliable Distance Estimation on Smartphones

abstract: As a response to the global COVID-19 surge in 2020, many countries have implemented lockdown or stay-at-home policies. If, however, the contact persons of every infected patient could be identified, the number of virus transmissions could be reduced, while the more incisive measures could be softened. For this purpose, contact tracing using smartphones is being considered as a promising technique. Here, smartphones emit and scan for Bluetooth Low Energy (BLE) signals for detecting devices in range. When a device is detected, its distance is estimated by evaluating its received signal strength. The main insight that is exploited for distance estimation is that the attenuation of a signal increases with the distance along which it has traveled. However, besides distance, there are multiple additional factors that impact the attenuation and hence disturb distance estimation. Among them, frequency-selective hardware and signal propagation belong to the most significant ones. For example, a BLE device transmits beacons on three different frequencies (channels), while the transmit power and the receiver sensitivity depend on the frequency. As a result, the received signal strength varies for each channel, even when the distance remains constant. However, the information on which wireless channel a beacon has been received is not made available to a smartphone. Hence, this error cannot be compensated, e.g., by calibration. In this paper, we for the first time provide a solution to detect the wireless channel on which a packet has been received on a smartphone. We experimentally evaluate our proposed technique on multiple different smartphone models. Our results help to make contact tracing more robust by improving the accuracy of distance estimation.


title: Learning to Solve Vehicle Routing Problems with Time Windows through Joint Attention

abstract: Many real-world vehicle routing problems involve rich sets of constraints with respect to the capacities of the vehicles, time windows for customers etc. While in recent years first machine learning models have been developed to solve basic vehicle routing problems faster than optimization heuristics, complex constraints rarely are taken into consideration. Due to their general procedure to construct solutions sequentially route by route, these methods generalize unfavorably to such problems. In this paper, we develop a policy model that is able to start and extend multiple routes concurrently by using attention on the joint action space of several tours. In that way the model is able to select routes and customers and thus learns to make difficult trade-offs between routes. In comprehensive experiments on three variants of the vehicle routing problem with time windows we show that our model called JAMPR works well for different problem sizes and outperforms the existing state-of-the-art constructive model. For two of the three variants it also creates significantly better solutions than a comparable meta-heuristic solver.


title: Graph manifolds that admit arbitrarily many Anosov flows

abstract: For each natural number n, we construct an example of a graph manifold supporting at least n different Anosov flows that are not orbit equivalent. Our construction is reminiscent of the Thurston-Handel construction: we cut a geodesic flow on a surface of constant negative curvature into two pieces, modify the flow in each piece by pulling back to finite covers, and glue together compatible pairs of pullback flows along their boundary tori to get many distinct flows on the resulting graph manifold.


title: UCSG-Net -- Unsupervised Discovering of Constructive Solid Geometry Tree

abstract: Signed distance field (SDF) is a prominent implicit representation of 3D meshes. Methods that are based on such representation achieved state-of-the-art 3D shape reconstruction quality. However, these methods struggle to reconstruct non-convex shapes. One remedy is to incorporate a constructive solid geometry framework (CSG) that represents a shape as a decomposition into primitives. It allows to embody a 3D shape of high complexity and non-convexity with a simple tree representation of Boolean operations. Nevertheless, existing approaches are supervised and require the entire CSG parse tree that is given upfront during the training process. On the contrary, we propose a model that extracts a CSG parse tree without any supervision - UCSG-Net. Our model predicts parameters of primitives and binarizes their SDF representation through differentiable indicator function. It is achieved jointly with discovering the structure of a Boolean operators tree. The model selects dynamically which operator combination over primitives leads to the reconstruction of high fidelity. We evaluate our method on 2D and 3D autoencoding tasks. We show that the predicted parse tree representation is interpretable and can be used in CAD software.


title: Jackson-type inequalities and widths of functional classes in the Musielak-Orlicz type spaces

abstract: In the Musielak-Orlicz type spaces ${\mathcal S}_{\bf M}$, exact Jackson-type inequalities are obtained in terms of best approximations of functions and the averaged values of their generalized moduli of smoothness. The values of Kolmogorov, Bernstein, linear, and projective widths in ${\mathcal S}_{\bf M}$ are found for classes of periodic functions defined by certain conditions on the averaged values of the generalized moduli of smoothness.


title: New Interpretations of Normalization Methods in Deep Learning

abstract: In recent years, a variety of normalization methods have been proposed to help train neural networks, such as batch normalization (BN), layer normalization (LN), weight normalization (WN), group normalization (GN), etc. However, mathematical tools to analyze all these normalization methods are lacking. In this paper, we first propose a lemma to define some necessary tools. Then, we use these tools to make a deep analysis on popular normalization methods and obtain the following conclusions: 1) Most of the normalization methods can be interpreted in a unified framework, namely normalizing pre-activations or weights onto a sphere; 2) Since most of the existing normalization methods are scaling invariant, we can conduct optimization on a sphere with scaling symmetry removed, which can help stabilize the training of network; 3) We prove that training with these normalization methods can make the norm of weights increase, which could cause adversarial vulnerability as it amplifies the attack. Finally, a series of experiments are conducted to verify these claims.


title: Moir\'e excitons in MoSe$_2$-WSe$_2$ heterobilayers and heterotrilayers

abstract: Layered two-dimensional materials exhibit rich transport and optical phenomena in twisted or lattice-incommensurate heterostructures with spatial variations of interlayer hybridization arising from moir\'e interference effects. Here, we report experimental and theoretical studies of excitons in twisted heterobilayers and heterotrilayers of transition metal dichalcogenides. Using MoSe$_2$-WSe$_2$ stacks as representative realizations of twisted van der Waals bilayer and trilayer heterostructures, we observe contrasting optical signatures and interpret them in the theoretical framework of interlayer moir\'e excitons in different spin and valley configurations. We conclude that the photoluminescence of MoSe$_2$-WSe$_2$ heterobilayer is consistent with joint contributions from radiatively decaying valley-direct interlayer excitons and phonon-assisted emission from momentum-indirect reservoirs that reside in spatially distinct regions of moir\'e supercells, whereas the heterotrilayer emission is entirely due to momentum-dark interlayer excitons of hybrid-layer valleys. Our results highlight the profound role of interlayer hybridization for transition metal dichalcogenide heterostacks and other realizations of multi-layered semiconductor van der Waals heterostructures.


title: The distributions under two species-tree models of the number of root ancestral configurations for matching gene trees and species trees

abstract: For a pair consisting of a gene tree and a species tree, the ancestral configurations at an internal node of the species tree are the distinct sets of gene lineages that can be present at that node. Ancestral configurations appear in computations of gene tree probabilities under evolutionary models conditional on fixed species trees, and the enumeration of root ancestral configurations -- ancestral configurations at the root of the species tree -- assists in describing the complexity of these computations. In the case that the gene tree matches the species tree in topology, we study the distribution of the number of root ancestral configurations of a random labeled tree topology under each of two models.


title: Learning from Demonstration with Weakly Supervised Disentanglement

abstract: Robotic manipulation tasks, such as wiping with a soft sponge, require control from multiple rich sensory modalities. Human-robot interaction, aimed at teaching robots, is difficult in this setting as there is potential for mismatch between human and machine comprehension of the rich data streams. We treat the task of interpretable learning from demonstration as an optimisation problem over a probabilistic generative model. To account for the high-dimensionality of the data, a high-capacity neural network is chosen to represent the model. The latent variables in this model are explicitly aligned with high-level notions and concepts that are manifested in a set of demonstrations. We show that such alignment is best achieved through the use of labels from the end user, in an appropriately restricted vocabulary, in contrast to the conventional approach of the designer picking a prior over the latent variables. Our approach is evaluated in the context of two table-top robot manipulation tasks performed by a PR2 robot -- that of dabbing liquids with a sponge (forcefully pressing a sponge and moving it along a surface) and pouring between different containers. The robot provides visual information, arm joint positions and arm joint efforts. We have made videos of the tasks and data available - see supplementary materials at: https://sites.google.com/view/weak-label-lfd.


title: An STPA-based Approach for Systematic Security Analysis of In-vehicle Diagnostic and Software Update Systems

abstract: The in-vehicle diagnostic and software update system, which supports remote diagnostic and Over-The-Air (OTA) software updates, is a critical attack goal in automobiles. Adversaries can inject malicious software into vehicles or steal sensitive information through communication channels. Therefore, security analysis, which identifies potential security issues, needs to be conducted in system design. However, existing security analyses of in-vehicle systems are threat-oriented, which start with threat identification and assess risks by brainstorming. In this paper, a system-oriented approach is proposed on the basis of the System-Theoretic Process Analysis (STPA). The proposed approach extends the original STPA from the perspective of data flows and is applicable for information-flow-based systems. Besides, we propose a general model for in-vehicle diagnostic and software update systems and use it to establish a security analysis guideline. In comparison with threat-oriented approaches, the proposed approach shifts from focusing on threats to system vulnerabilities and seems to be efficient to prevent the system from known or even unknown threats. Furthermore, as an extension of the STPA, which has been proven to be applicable to high level designs, the proposed approach can be well integrated into high-level analyses and perform co-design in different disciplines within a unified STPA framework.


title: How to Probe Sentence Embeddings in Low-Resource Languages: On Structural Design Choices for Probing Task Evaluation

abstract: Sentence encoders map sentences to real valued vectors for use in downstream applications. To peek into these representations - e.g., to increase interpretability of their results - probing tasks have been designed which query them for linguistic knowledge. However, designing probing tasks for lesser-resourced languages is tricky, because these often lack large-scale annotated data or (high-quality) dependency parsers as a prerequisite of probing task design in English. To investigate how to probe sentence embeddings in such cases, we investigate sensitivity of probing task results to structural design choices, conducting the first such large scale study. We show that design choices like size of the annotated probing dataset and type of classifier used for evaluation do (sometimes substantially) influence probing outcomes. We then probe embeddings in a multilingual setup with design choices that lie in a 'stable region', as we identify for English, and find that results on English do not transfer to other languages. Fairer and more comprehensive sentence-level probing evaluation should thus be carried out on multiple languages in the future.


title: Bidirectional Teleportation using Fisher Information

abstract: In this contribution, we reformulated the bidirectional teleportation protocol suggested in [7], by means of Bloch vectors as well as the local operations are represented by using Pauli operators. Analytical and numerical calculations for the teleported state and Fisher information are introduced. It is shown that both quantities depend on the initial state settings of the teleported qubits and their triggers. The Fidelities and the Fisher information of the bidirectionally teleported states are maximized when the qubit and its trigger are polarized in the same direction. The minimum values are predicted if both initial qubits have different polarization or non-zero phase. The maximum values of the Fidelity and the quantum Fisher information are the same, but they are predicted at different polarization angles. We display that the multi-parameter form is much better than the single parameter form, where it satisfies the bounds of classical, entangled systems and the uncertainty principle.


title: Unified SVM Algorithm Based on LS-DC Loss

abstract: Over the past two decades, support vector machine (SVM) has become a popular supervised machine learning model, and plenty of distinct algorithms are designed separately based on different KKT conditions of the SVM model for classification/regression with different losses, including the convex loss or nonconvex loss. In this paper, we propose an algorithm that can train different SVM models in a \emph{unified} scheme. First, we introduce a definition of the \emph{LS-DC} (\textbf{l}east \textbf{s}quares type of \textbf{d}ifference of \textbf{c}onvex) loss and show that the most commonly used losses in the SVM community are LS-DC loss or can be approximated by LS-DC loss. Based on DCA (difference of convex algorithm), we then propose a unified algorithm, called \emph{UniSVM}, which can solve the SVM model with any convex or nonconvex LS-DC loss, in which only a vector is computed, especially by the specifically chosen loss. Particularly, for training robust SVM models with nonconvex losses, UniSVM has a dominant advantage over all existing algorithms because it has a closed-form solution per iteration, while the existing algorithms always need to solve an L1SVM/L2SVM per iteration. Furthermore, by the low-rank approximation of the kernel matrix, UniSVM can solve the large-scale nonlinear problems with efficiency. To verify the efficacy and feasibility of the proposed algorithm, we perform many experiments on some small artificial problems and some large benchmark tasks with/without outliers for classification and regression for comparison with state-of-the-art algorithms. The experimental results demonstrate that UniSVM can achieve comparable performance in less training time. The foremost advantage of UniSVM is that its core code in Matlab is less than 10 lines; hence, it can be easily grasped by users or researchers.


title: Electroweak Top Couplings, Partial Compositeness and Top Partner Searches

abstract: Partial top quark compositeness is a crucial aspect of theories with strong electroweak symmetry breaking. Together with the heavy top partners that lift the top quark mass to its observed value, these theories predict correlated modifications of the top quark's electroweak couplings. Associated measurements therefore provide direct constraints on the ultraviolet structure of the underlying hypercolour dynamics. In this paper we employ a minimal version of top compositeness to discuss how measurements related to the top's electroweak gauge interactions can inform the potential composite nature of the TeV scale. In doing so, we identify the dominant factors that limit the BSM sensitivity. Extrapolating to a future 100 TeV hadron collider, we demonstrate that top quark measurements performed at highest precision can provide complementary information to resonance search by performing a representative resonant top partner search that specifically targets the correlated resonant electroweak top partner signatures.


title: Topological defects and confinement with machine learning: the case of monopoles in compact electrodynamics

abstract: We investigate the advantages of machine learning techniques to recognize the dynamics of topological objects in quantum field theories. We consider the compact U(1) gauge theory in three spacetime dimensions as the simplest example of a theory that exhibits confinement and mass gap phenomena generated by monopoles. We train a neural network with a generated set of monopole configurations to distinguish between confinement and deconfinement phases, from which it is possible to determine the deconfinement transition point and to predict several observables. The model uses a supervised learning approach and treats the monopole configurations as three-dimensional images (holograms). We show that the model can determine the transition temperature with accuracy, which depends on the criteria implemented in the algorithm. More importantly, we train the neural network with configurations from a single lattice size before making predictions for configurations from other lattice sizes, from which a reliable estimation of the critical temperatures are obtained.


title: Adversarial representation learning for private speech generation

abstract: As more and more data is collected in various settings across organizations, companies, and countries, there has been an increase in the demand of user privacy. Developing privacy preserving methods for data analytics is thus an important area of research. In this work we present a model based on generative adversarial networks (GANs) that learns to obfuscate specific sensitive attributes in speech data. We train a model that learns to hide sensitive information in the data, while preserving the meaning in the utterance. The model is trained in two steps: first to filter sensitive information in the spectrogram domain, and then to generate new and private information independent of the filtered one. The model is based on a U-Net CNN that takes mel-spectrograms as input. A MelGAN is used to invert the spectrograms back to raw audio waveforms. We show that it is possible to hide sensitive information such as gender by generating new data, trained adversarially to maintain utility and realism.


title: Discretization of the Lamperti representation of a positive self-similar Markov process

abstract: This paper considers discretization of the L\'evy process appearing in the Lamperti representation of a strictly positive self-similar Markov process. Limit theorems for the resulting approximation are established under some regularity assumptions on the given L\'evy process. Additionally, the scaling limit of a positive self-similar Markov process at small times is provided. Finally, we present an application to simulation of self-similar L\'evy processes conditioned to stay positive.


title: 1st place solution for AVA-Kinetics Crossover in AcitivityNet Challenge 2020

abstract: This technical report introduces our winning solution to the spatio-temporal action localization track, AVA-Kinetics Crossover, in ActivityNet Challenge 2020. Our entry is mainly based on Actor-Context-Actor Relation Network. We describe technical details for the new AVA-Kinetics dataset, together with some experimental results. Without any bells and whistles, we achieved 39.62 mAP on the test set of AVA-Kinetics, which outperforms other entries by a large margin. Code will be available at: https://github.com/Siyu-C/ACAR-Net.


title: End-to-End Real-time Catheter Segmentation with Optical Flow-Guided Warping during Endovascular Intervention

abstract: Accurate real-time catheter segmentation is an important pre-requisite for robot-assisted endovascular intervention. Most of the existing learning-based methods for catheter segmentation and tracking are only trained on small-scale datasets or synthetic data due to the difficulties of ground-truth annotation. Furthermore, the temporal continuity in intraoperative imaging sequences is not fully utilised. In this paper, we present FW-Net, an end-to-end and real-time deep learning framework for endovascular intervention. The proposed FW-Net has three modules: a segmentation network with encoder-decoder architecture, a flow network to extract optical flow information, and a novel flow-guided warping function to learn the frame-to-frame temporal continuity. We show that by effectively learning temporal continuity, the network can successfully segment and track the catheters in real-time sequences using only raw ground-truth for training. Detailed validation results confirm that our FW-Net outperforms state-of-the-art techniques while achieving real-time performance.


title: $Q$-learning with Logarithmic Regret

abstract: This paper presents the first non-asymptotic result showing that a model-free algorithm can achieve a logarithmic cumulative regret for episodic tabular reinforcement learning if there exists a strictly positive sub-optimality gap in the optimal $Q$-function. We prove that the optimistic $Q$-learning studied in [Jin et al. 2018] enjoys a ${\mathcal{O}}\left(\frac{SA\cdot \mathrm{poly}\left(H\right)}{\Delta_{\min}}\log\left(SAT\right)\right)$ cumulative regret bound, where $S$ is the number of states, $A$ is the number of actions, $H$ is the planning horizon, $T$ is the total number of steps, and $\Delta_{\min}$ is the minimum sub-optimality gap. This bound matches the information theoretical lower bound in terms of $S,A,T$ up to a $\log\left(SA\right)$ factor. We further extend our analysis to the discounted setting and obtain a similar logarithmic cumulative regret bound.


title: Query Intent Detection from the SEO Perspective

abstract: Google users have different intents from their queries such as acquiring information, buying products, comparing or simulating services, looking for products, and so on. Understanding the right intention of users helps to provide i) better content on web pages from the Search Engine Optimization (SEO) perspective and ii) more user-satisfying results from the search engine perspective. In this study, we aim to identify the user query's intent by taking advantage of Google results and machine learning methods. Our proposed approach is a clustering model that exploits some features to detect query's intent. A list of keywords extracted from the clustered queries is used to identify the intent of a new given query. Comparing the clustering results with the intents predicted by filtered keywords show the efficiency of the extracted keywords for detecting intents.


title: Continuity method with movable singularities for classical Monge-Amp\`ere equations

abstract: On a compact K\"ahler manifold $(X,\omega)$, we study the strong continuity of solutions with prescribed singularities of complex Monge-Amp\`ere equations with integrable Lebesgue densities. Moreover, we give sufficient conditions for the strong continuity of solutions when the right-hand sides are modified to include all (log) K\"ahler-Einstein metrics with prescribed singularities. Our findings can be interpreted as closedness of new continuity methods in which the densities vary together with the prescribed singularities. For Monge-Amp\`ere equations of Fano type, we also prove an openness result when the singularities decrease. As an application, we deduce a strong stability result for (log-)K\"ahler Einstein metrics on semi-K\"ahler classes given as modifications of $\{\omega\}$.


title: Black-Hole evaporation and Quantum-depletion in Bose-Einstein condensates

abstract: We study the analogy between the Hawking radiation in Black-Holes and the Quantum depletion process of a Bose-Einstein condensate by using the Bogoliubov transformations method. We find that the relation between the Bogoliubov coefficients is similar in both cases (in the appropriate regimes). We then connect the condensate variables with those associated to the Black-Hole, demonstrating then that the zero temperature regime of the condensate is equivalent to the existence of an event horizon in gravity.

categories: cond-mat.quant-gas cond-mat.stat-mech gr-qc hep-th


title: Transfer of Manure as Fertilizer from Livestock Farms to Crop Fields: The Case of Catalonia

abstract: Intensive livestock production might have a negative environmental impact, by producing large amounts of animal manure, which, if not properly managed, can contaminate nearby water bodies with nutrient excess. However, if animal manure is exported to nearby crop fields, to be used as organic fertilizer, pollution can be mitigated. It is a single-objective optimization problem, in regards to finding the best solution for the logistics process of satisfying nutrient needs of crops by means of livestock manure. This paper proposes three different approaches to solve the problem: a centralized optimal algorithm (COA), a decentralized nature-inspired cooperative technique, based on the foraging behaviour of ants (AIA), as well as a naive neighbour-based method (NBS), which constitutes the existing practice used today in an ad hoc, uncoordinated manner in Catalonia. Results show that the COA approach is 8.5% more efficient than the AIA. However, the AIA approach is fairer to the farmers and more balanced in terms of average transportation distances that need to be covered by each livestock farmer, while it is 1.07 times more eefficient than the NBS. Our work constitutes the first application of a decentralized AIA to this interesting real-world problem, in a domain where swarm intelligence methods are still under-exploited.


title: Algorithms with Predictions

abstract: We introduce algorithms that use predictions from machine learning applied to the input to circumvent worst-case analysis. We aim for algorithms that have near optimal performance when these predictions are good, but recover the prediction-less worst case behavior when the predictions have large errors.


title: Dynamically induced Planck scale and inflation in the Palatini formulation

abstract: We study non-minimal Coleman-Weinberg inflation in the Palatini formulation of gravity in the presence of an $R^2$ term. The Planck scale is dynamically generated by the vacuum expectation value of the inflaton via its non-minimal coupling to the curvature scalar $R$. We show that the addition of the $R^2$ term in Palatini gravity makes non-minimal Coleman-Weinberg inflation again compatible with observational data.


title: A small molecule drug candidate targeting SARS-CoV-2 main protease

abstract: A new coronavirus identified as SARS-CoV-2 virus has brought the world to a state of crisis, causing a major pandemic, claiming more than 433,000 lives and instigating major financial damage to the global economy. Despite current efforts, developing safe and effective treatments remains a major challenge. Moreover, new strains of the virus are likely to emerge in the future. To prevent future pandemics, several drugs with various mechanisms of action are required. Drug discovery efforts against the virus fall into two main categories: (a) monoclonal antibodies targeting the spike protein of the virus and blocking it from entry; (b) small molecule inhibitors targeting key proteins of the virus, interfering with replication and translation of the virus. In this study, we are presenting a computational investigation of a potential drug candidate that targets SARS-CoV-2 protease, a viral protein critical for replication and translation of the virus.


title: Evolutionary Algorithms with Self-adjusting Asymmetric Mutation

abstract: Evolutionary Algorithms (EAs) and other randomized search heuristics are often considered as unbiased algorithms that are invariant with respect to different transformations of the underlying search space. However, if a certain amount of domain knowledge is available the use of biased search operators in EAs becomes viable. We consider a simple (1+1) EA for binary search spaces and analyze an asymmetric mutation operator that can treat zero- and one-bits differently. This operator extends previous work by Jansen and Sudholt (ECJ 18(1), 2010) by allowing the operator asymmetry to vary according to the success rate of the algorithm. Using a self-adjusting scheme that learns an appropriate degree of asymmetry, we show improved runtime results on the class of functions OneMax$_a$ describing the number of matching bits with a fixed target $a\in\{0,1\}^n$.


title: Quantum simulation and circuit design for solving multidimensional Poisson equations

abstract: Many methods solve Poisson equations by using grid techniques which discretize the problem in each dimension. Most of these algorithms are subject to the curse of dimensionality, so that they need exponential runtime. In the paper "Quantum algorithm and circuit design solving the Poisson equation" a quantum algorithm is shown running in polylog time to produce a quantum state representing the solution of the Poisson equation. In this paper a quantum simulation of an extended circuit design based on this algorithm is made on a classical computer. Our purpose is to test an efficient circuit design which can break the curse of dimensionality on a quantum computer. Due to the exponential rise of the Hilbert space this design is optimized on a small number of qubits. We use Microsoft's Quantum Development Kit and its simulator of an ideal quantum computer to validate the correctness of this algorithm.


title: Rethinking the Role of Gradient-Based Attribution Methods for Model Interpretability

abstract: Current methods for the interpretability of discriminative deep neural networks commonly rely on the model's input-gradients, i.e., the gradients of the output logits w.r.t. the inputs. The common assumption is that these input-gradients contain information regarding $p_{\theta} ( y \mid x)$, the model's discriminative capabilities, thus justifying their use for interpretability. However, in this work we show that these input-gradients can be arbitrarily manipulated as a consequence of the shift-invariance of softmax without changing the discriminative function. This leaves an open question: if input-gradients can be arbitrary, why are they highly structured and explanatory in standard models? We investigate this by re-interpreting the logits of standard softmax-based classifiers as unnormalized log-densities of the data distribution and show that input-gradients can be viewed as gradients of a class-conditional density model $p_{\theta}(x \mid y)$ implicit within the discriminative model. This leads us to hypothesize that the highly structured and explanatory nature of input-gradients may be due to the alignment of this class-conditional model $p_{\theta}(x \mid y)$ with that of the ground truth data distribution $p_{\text{data}} (x \mid y)$. We test this hypothesis by studying the effect of density alignment on gradient explanations. To achieve this alignment we use score-matching, and propose novel approximations to this algorithm to enable training large-scale models. Our experiments show that improving the alignment of the implicit density model with the data distribution enhances gradient structure and explanatory power while reducing this alignment has the opposite effect. Overall, our finding that input-gradients capture information regarding an implicit generative model implies that we need to re-think their use for interpreting discriminative models.


title: Model Explanations with Differential Privacy

abstract: Black-box machine learning models are used in critical decision-making domains, giving rise to several calls for more algorithmic transparency. The drawback is that model explanations can leak information about the training data and the explanation data used to generate them, thus undermining data privacy. To address this issue, we propose differentially private algorithms to construct feature-based model explanations. We design an adaptive differentially private gradient descent algorithm, that finds the minimal privacy budget required to produce accurate explanations. It reduces the overall privacy loss on explanation data, by adaptively reusing past differentially private explanations. It also amplifies the privacy guarantees with respect to the training data. We evaluate the implications of differentially private models and our privacy mechanisms on the quality of model explanations.


title: K\"ahler-Einstein metrics with prescribed singularities on Fano manifolds

abstract: Given a Fano manifold $(X,\omega)$ we develop a variational approach to characterize analytically the existence of K\"ahler-Einstein metrics with prescribed singularities, assuming that these singularities can be approximated algebraically. Moreover, we define a function $\alpha_{\omega}$ on the set of prescribed singularities which generalizes Tian's $\alpha$-invariant, showing that its upper level set $\{\alpha_{\omega}(\cdot)>\frac{n}{n+1}\}$ produces a subset of the K\"ahler-Einstein locus, i.e. of the locus given by all prescribed singularities that admit K\"ahler-Einstein metrics. In particular, we prove that many $K$-stable manifolds admit all possible K\"ahler-Einstein metrics with prescribed singularities. Conversely, we show that enough positivity of the $\alpha$-invariant function at non-trivial prescribed singularities (or other conditions) implies the existence of genuine K\"ahler-Einstein metrics. Finally, through a continuity method, we also prove the strong continuity of K\"ahler-Einstein metrics on curves of totally ordered prescribed singularities when the relative automorphism groups are discrete.


title: Observation of cyclotron resonance and measurement of the hole mass in optimally-doped La$_{2-x}$Sr$_{x}$CuO$_4$

abstract: Using time-domain terahertz spectroscopy in pulsed magnetic fields up to 31 T, we measure the terahertz optical conductivity in an optimally-doped thin film of the high temperature superconducting cuprate La$_{1.84}$Sr$_{0.16}$CuO$_4$. We observe systematic changes in the circularly-polarized complex optical conductivity that are consistent with cyclotron absorption of p-type charge carriers characterized by a cyclotron mass of $4.9\pm 0.8$ $m_{\rm e}$, and a scattering rate that increases with magnetic field. These results open the door to studies aimed at characterizing the degree to which electron-electron interactions influence carrier masses in cuprate superconductors.


title: On the Decidability of Reachability in Continuous Time Linear Time-Invariant Systems

abstract: We consider the decidability of state-to-state reachability in linear time-invariant control systems over continuous time. We analyse this problem with respect to the allowable control sets, which are assumed to be the image under a linear map of the unit hypercube. This naturally models bounded (sometimes called saturated) controls. Decidability of the version of the reachability problem in which control sets are affine subspaces of $\mathbb{R}^n$ is a fundamental result in control theory. Our first result is decidability in two dimensions ($n=2$) if the matrix $A$ satisfies some spectral conditions, and conditional decidablility in general. If the transformation matrix $A$ is diagonal with rational entries (or rational multiples of the same algebraic number) then the reachability problem is decidable. If the transformation matrix $A$ only has real eigenvalues, the reachability problem is conditionally decidable. The time-bounded reachability problem is conditionally decidable, and unconditionally decidable in two dimensions. Some of our decidability results are conditional in that they rely on the decidability of certain mathematical theories, namely the theory of the reals with exponential ($\mathfrak{R}_{\exp}$) and with bounded sine ($\mathfrak{R}_{\exp,\sin}$). We also obtain a hardness result for a mild generalization of the problem where the target is simple set (hypercube of dimension $n-1$ or hyperplane) instead of a point, and the control set is a convex bounded polytope. In this case, we show that the problem is at least as hard as the \emph{Continuous Positivity problem} or the \emph{Nontangential Continuous Positivity problem}.


title: Gradient formula for transition semigroup corresponding to stochastic equation driven by a system of independent L\'evy processes

abstract: Let $(P_t)$ be the transition semigroup of the Markov family $(X^x(t))$ defined by SDE $$ d X= b(X) dt + d Z, \qquad X(0)=x, $$ where $Z=\left(Z_1, \ldots, Z_d\right)^*$ is a system of independent real-valued L\'evy processes. Using the Malliavin calculus we establish the following gradient formula $$ \nabla P_tf(x)= \mathbb{E}\, f\left(X^x(t)\right) Y(t,x), \qquad f\in B_b(\mathbb{R}^d), $$ where the random field $Y$ does not depend on $f$. Sharp estimates on $\nabla P_tf(x)$ when $Z_1, \ldots , Z_d$ are $\alpha$-stable processes, $\alpha \in (0,2)$, are also given.


title: AlphaGAN: Fully Differentiable Architecture Search for Generative Adversarial Networks

abstract: Generative Adversarial Networks (GANs) are formulated as minimax game problems, whereby generators attempt to approach real data distributions by virtue of adversarial learning against discriminators. The intrinsic problem complexity poses the challenge to enhance the performance of generative networks. In this work, we aim to boost model learning from the perspective of network architectures, by incorporating recent progress on automated architecture search into GANs. To this end, we propose a fully differentiable search framework for generative adversarial networks, dubbed alphaGAN. The searching process is formalized as solving a bi-level minimax optimization problem, in which the outer-level objective aims for seeking a suitable network architecture towards pure Nash Equilibrium conditioned on the generator and the discriminator network parameters optimized with a traditional GAN loss in the inner level. The entire optimization performs a first-order method by alternately minimizing the two-level objective in a fully differentiable manner, enabling architecture search to be completed in an enormous search space. Extensive experiments on CIFAR-10 and STL-10 datasets show that our algorithm can obtain high-performing architectures only with 3-GPU hours on a single GPU in the search space comprised of approximate 2 ? 1011 possible configurations. We also provide a comprehensive analysis on the behavior of the searching process and the properties of searched architectures, which would benefit further research on architectures for generative models. Pretrained models and codes are available at https://github.com/yuesongtian/AlphaGAN.


title: Chandra spectral measurements of the O supergiant $\zeta$ Puppis indicate a surprising increase in the wind mass-loss rate over 18 years

abstract: New long Chandra grating observations of the O supergiant $\zeta$ Pup show not only a brightening of the x-ray emission line flux of 13 per cent in the 18 years since Chandra's first observing cycle, but also clear evidence - at more than four sigma significance - of increased wind absorption signatures in its Doppler-broadened x-ray emission line profiles. We demonstrate this with non-parametric analysis of the profiles as well as Gaussian fitting and then use the line-profile model fitting to derive a mass-loss rate of $2.47 \pm 0.09 \times 10^{-6}$ Msun/yr, which is a 40 per cent increase over the value obtained from the cycle 1 data. The increase in the individual emission line fluxes is greater for short-wavelength lines than long-wavelength lines, as would be expected if a uniform increase in line emission is accompanied by an increase in the wavelength-dependent absorption by the cold wind in which the shock-heated plasma is embedded.


title: When Does Self-Supervision Help Graph Convolutional Networks?

abstract: Self-supervision as an emerging technique has been employed to train convolutional neural networks (CNNs) for more transferrable, generalizable, and robust representation learning of images. Its introduction to graph convolutional networks (GCNs) operating on graph data is however rarely explored. In this study, we report the first systematic exploration and assessment of incorporating self-supervision into GCNs. We first elaborate three mechanisms to incorporate self-supervision into GCNs, analyze the limitations of pretraining & finetuning and self-training, and proceed to focus on multi-task learning. Moreover, we propose to investigate three novel self-supervised learning tasks for GCNs with theoretical rationales and numerical comparisons. Lastly, we further integrate multi-task self-supervision into graph adversarial training. Our results show that, with properly designed task forms and incorporation mechanisms, self-supervision benefits GCNs in gaining more generalizability and robustness. Our codes are available at https://github.com/Shen-Lab/SS-GCNs.


title: Multipole analysis of substrate-supported dielectric nanoresonator arrays with T-matrix method

abstract: Substrates, and layered media in general, are ubiquitous, affect the properties of whatever is in their vicinity, and their influence is, in an arbitrary framework, challenging to quantify analytically, especially for large arrays which escape explicit numerical treatment due to the computational burden. In this work, we develop a versatile T-matrix based framework in which we generalize the coupled multipole model towards arbitrarily high multipole orders and substrate-supported arrays. It allows us to study substrate-supported random/amorphous arrays of high index dielectric nanoparticles which are of wide interest due to relatively low losses and a highly tunable optical response, making them promising elements for nanophotonic devices. We discuss how multipole coupling rules evolve in the presence of a substrate in amorphous arrays for three interaction mechanisms: direct coupling between particles, substrate-mediated interparticle coupling and substrate-mediated self-coupling. We show the interplay between array density, distance from the substrate and its refractive in determining the optical response of an array. As an example, we use this framework to analyze refractometric sensing with substrate-supported arrays and demonstrate that the substrate plays a crucial role in determining the array sensitivity.


title: Multi-level Monte Carlo path integral molecular dynamics for thermal average calculation in the nonadiabatic regime

abstract: With the path integral approach, the thermal average in a multi-electronic-state quantum systems can be approximated by the ring polymer representation on an extended configuration space, where the additional degrees of freedom are associated with the surface index of each bead. The primary goal of this work is to propose a more efficient sampling algorithm for the calculation of such thermal averages. We reformulate the extended ring polymer approximation according to the configurations of the surface indexes, and by introducing a proper reference measure, the reformulation is recast as a ratio of two expectations of function expansions. By quantitatively estimating the sub-estimators, and minimizing the total variance of the sampled average, we propose a multi-level Monte Carlo path integral molecular dynamics method (MLMC-PIMD) to achieve an optimal balance of computational cost and accuracy.


title: Finite Groups with Some s-semipermutable subgroups

abstract: Suppose that $G$ is a finite group and $H$ is a subgroup of $G$. We say that $H$ is s-semipermutable in $G$ if $HG_p = G_pH$ for any Sylow $p$-subgroup $G_p$ of $G$ with $(p, |H|) = 1$. We investigate the influence of s-semipermutable subgroups on the structure of finite groups. Some recent results are generalized.


title: Perpetual Integral Functionals of Multidimensional Stochastic Processes

abstract: The paper is devoted to the existence of integral functionals $\int_0^\infty f(X(t))\,{\mathrm{d}t}$ for several classes of processes in $\mathbb{R}$ with $d\ge 3$. Some examples such as Brownian motion, fractional Brownian motion, compound Poisson process, Markov processes admitting densities of transitional probabilities are considered.


title: Improving accuracy and speeding up Document Image Classification through parallel systems

abstract: This paper presents a study showing the benefits of the EfficientNet models compared with heavier Convolutional Neural Networks (CNNs) in the Document Classification task, essential problem in the digitalization process of institutions. We show in the RVL-CDIP dataset that we can improve previous results with a much lighter model and present its transfer learning capabilities on a smaller in-domain dataset such as Tobacco3482. Moreover, we present an ensemble pipeline which is able to boost solely image input by combining image model predictions with the ones generated by BERT model on extracted text by OCR. We also show that the batch size can be effectively increased without hindering its accuracy so that the training process can be sped up by parallelizing throughout multiple GPUs, decreasing the computational time needed. Lastly, we expose the training performance differences between PyTorch and Tensorflow Deep Learning frameworks.


title: Cogradient Descent for Bilinear Optimization

abstract: Conventional learning methods simplify the bilinear model by regarding two intrinsically coupled factors independently, which degrades the optimization procedure. One reason lies in the insufficient training due to the asynchronous gradient descent, which results in vanishing gradients for the coupled variables. In this paper, we introduce a Cogradient Descent algorithm (CoGD) to address the bilinear problem, based on a theoretical framework to coordinate the gradient of hidden variables via a projection function. We solve one variable by considering its coupling relationship with the other, leading to a synchronous gradient descent to facilitate the optimization procedure. Our algorithm is applied to solve problems with one variable under the sparsity constraint, which is widely used in the learning paradigm. We validate our CoGD considering an extensive set of applications including image reconstruction, inpainting, and network pruning. Experiments show that it improves the state-of-the-art by a significant margin.


title: Ionized Gas in the NGC 5253 Supernebula:High Spatial and Spectral Resolution Observations with the JVLA and TEXES

abstract: The youngest, closest and most compact embedded massive star cluster known excites the supernebula in the nearby dwarf galaxy NGC 5253. It is a crucial target and test case for studying the birth and evolution of the most massive star clusters. We present observations of the ionized gas in this source with high spatial and spectral resolution. The data includes continuum images of free-free emission with ~0.15'' resolution made with the JVLA at 15, 22 and 33 GHz, and a full data cube of the [SIV]10.5 micron fine-structure emission line with ~4.5 km/s velocity resolution and 0.3'' beam, obtained with TEXES on Gemini North. We find that 1) the ionized gas extends out from the cluster in arms or jets, and 2) the ionized gas comprises two components offset both spatially and in velocity. We discuss mechanisms that may have created the observed velocity field; possibilities include large-scale jets or a subcluster falling onto the main source.


title: Totally Asynchronous Large-Scale Quadratic Programming: Regularization, Convergence Rates, and Parameter Selection

abstract: Quadratic programs arise in robotics, communications, smart grids, and many other applications. As these problems grow in size, finding solutions becomes more computationally demanding, and new algorithms are needed to efficiently solve them at massive scales. Targeting large-scale problems, we develop a multi-agent quadratic programming framework in which each agent updates only a small number of the total decision variables in a problem. Agents communicate their updated values to each other, though we do not impose any restrictions on the timing with which they do so, nor on the delays in these transmissions. Furthermore, we allow weak parametric coupling among agents, in the sense that they are free to independently choose their step sizes, subject to mild restrictions. We further provide the means for agents to independently regularize the problems they solve, thereby improving convergence properties while preserving agents' independence in selecting parameters and ensuring a global bound on regularization error is satisfied. Larger regularizations accelerate convergence but increase error in the solution obtained, and we quantify the trade off between convergence rates and quality of solutions. Simulation results are presented to illustrate these developments.


title: End-to-End Nanophotonic Inverse Design for Imaging and Polarimetry

abstract: By co-designing a meta-optical front end in conjunction with an image-processing back end, we demonstrate noise sensitivity and compactness substantially superior to either an optics-only or a computation-only approach, illustrated by two examples: subwavelength imaging and reconstruction of the full polarization coherence matrices of multiple light sources. Our end-to-end inverse designs couple the solution of the full Maxwell equations---exploiting all aspects of wave physics arising in subwavelength scatterers---with inverse-scattering algorithms in a single large-scale optimization involving $\gtrsim 10^4$ degrees of freedom. The resulting structures scatter light in a way that is radically different from either a conventional lens or a random microstructure, and suppress the noise sensitivity of the inverse-scattering computation by several orders of magnitude. Incorporating the full wave physics is especially crucial for detecting spectral and polarization information that is discarded by geometric optics and scalar diffraction theory.

categories: physics.optics eess.IV physics.app-ph physics.comp-ph