title: Efficient Fast Multipole Accelerated Boundary Elements via Recursive Computation of Multipole Expansions of Integrals

abstract: In boundary element methods (BEM) in $\mathbb{R}^3$, matrix elements and right hand sides are typically computed via analytical or numerical quadrature of the layer potential multiplied by some function over line, triangle and tetrahedral volume elements. When the problem size gets large, the resulting linear systems are often solved iteratively via Krylov subspace methods, with fast multipole methods (FMM) used to accelerate the matrix vector products needed. When FMM acceleration is used, most entries of the matrix never need be computed explicitly - {\em they are only needed in terms of their contribution to the multipole expansion coefficients.} We propose a new fast method - \emph{Quadrature to Expansion (Q2X)} - for the analytical generation of the multipole expansion coefficients produced by the integral expressions for single and double layers on surface triangles; charge distributions over line segments and over tetrahedra in the volume; so that the overall method is well integrated into the FMM, with controlled error. The method is based on the $O(1)$ per moment cost recursive computation of the moments. The method is developed for boundary element methods involving the Laplace Green's function in ${\mathbb R}^3$. The derived recursions are first compared against classical quadrature algorithms, and then integrated into FMM accelerated boundary element and vortex element methods. Numerical tests are presented and discussed.


title: Some Arguments for the Wave Equation in Quantum Theory

abstract: We clarify some arguments concerning Jefimenko's equations, as a way of constructing solutions to Maxwell's equations, for charge and current satisfying the continuity equation. We then isolate a condition on non-radiation in all inertial frames, which is intuitively reasonable for the stability of an atomic system, and prove that the condition is equivalent to the charge and current satisfying certain relations, including the wave equations. Finally, we prove that with these relations, the energy in the electromagnetic field is quantised and displays the properties of the Balmer series.

categories: math.AP math-ph math.MP physics.class-ph quant-ph


title: Demystifying speckle field quantitative phase microscopy

abstract: Quantitative phase microscopy (QPM) has found significant applications in the field of biomedical imaging which works on the principle of interferometry. The theory behind achieving interference in QPM with conventional light sources such as white light and lasers is very well developed. Recently, the use of dynamic speckle illumination (DSI) in QPM has attracted attention due to its advantages over conventional light sources such as high spatial phase sensitivity, single shot, scalable field of view (FOV) and resolution. However, the understanding behind obtaining interference fringes in QPM with DSI has not been convincingly covered previously. This imposes a constraint on obtaining interference fringes in QPM using DSI and limits its widespread penetration in the field of biomedical imaging. The present article provides the basic understanding of DSI through both simulation and experiments that is essential to build interference optical microscopy systems such as QPM, digital holographic microscopy and optical coherence tomography. Using the developed theory of DSI we demonstrate its capabilities of using non-identical objective lenses in both arms of the interference microscopy without degrading the interference fringe contrast and providing the flexibility to use user-defined microscope objective lens. It is also demonstrated that the interference fringes are not washed out over a large range of optical path difference (OPD) between the object and the reference arm providing competitive edge over low temporal coherence light sources. The theory and explanation developed here would enable wider penetration of DSI based QPM for applications in biology and material sciences.


title: Selective and enhanced nickel adsorption from sulfate- and calcium-rich solutions using chitosan

abstract: Nickel (Ni) is an economically important metal characterized by its mechanical strength and anticorrosion properties. With the increasing industrial demand for Ni and the depleting accessible Ni primary ores, sustainable technologies are required for the recovery of this metal from alternative resources. In this study, adsorption using chitosan was investigated as a sustainable technique to recover Ni from sulfate (SO42-) and calcium (Ca) rich secondary resources. The effects of pH, contact time, and the presence of SO42- and Ca on Ni adsorption were investigated in batch experiments. Chemical speciation modeling was performed to analyze how the predominant Ni species present under different conditions may affect the efficiency of the adsorption process. The comparison of the maximum Ni adsorption capacities of chitosan in the absence (1.00 mmol/g) and presence (1.49 mmol/g) of 500 mM SO42- and 10 mM Ca indicated the positive effect of these ions on Ni adsorption. The predominance of the neutrally charged (NiSO4)0 species in Ni/SO42- system has contributed to the enhanced Ni adsorption on chitosan as verified by X-ray photoelectron spectroscopy (XPS) analysis. However, kinetic studies confirmed that the Ni adsorption rate decreased by 4.5 times when SO42- was present. The subsequent continuous Ni adsorption from a real SO42-rich leachate in a column setup revealed that chitosan is selective for Ni over Ca and Cr with selectivity quotients of 9.6 (K Ni/Ca) and 3.0 (K Ni/Cr). Overall, this study indicated that Ni complexation with SO42- enhances the Ni adsorption capacity of chitosan, but slows down the adsorption process.


title: LMC N132D: a mature supernova remnant with a youthful gamma-ray spectrum

abstract: The supernova remnant LMC N132D is a remarkably luminous gamma-ray emitter at $\sim$50 kpc with an age of $\sim$2500 years. It belongs to the small group of oxygen-rich SNRs, which includes Cassiopeia A (Cas A) and Puppis A. N132D is interacting with a nearby molecular cloud. By adding 102 hours of new observations with the High Energy Stereoscopic System (H.E.S.S.) to the previously published data with exposure time of 150 hours, we achieve the significant detection of N132D at a 5.7$\sigma$ level in the very high energy (VHE) domain. The gamma-ray spectrum is compatible with a single power law extending above 10 TeV. We set a lower limit on an exponential cutoff energy at 8 TeV with 95% CL. The multi-wavelength study supports a hadronic origin of VHE gamma-ray emission indicating the presence of sub-PeV cosmic-ray protons. The detection of N132D is remarkable since the TeV luminosity is higher than that of Cas A by more than an order of magnitude. Its luminosity is comparable to, or even exceeding the luminosity of RX J1713.7-3946 or HESS J1640-465. Moreover, the extended power-law tail in the VHE spectrum of N132D is surprising given both the exponential cutoff at 3.5 TeV in the spectrum of its 340-year-old sibling, Cassiopeia A, and the lack of TeV emission from a Fermi- LAT 2FHL source (E > 50 GeV) associated with Puppis A. We discuss a physical scenario leading to the enhancement of TeV emission via the interaction between N132D and a near molecular cloud.


title: On Integral Theorems and their Statistical Properties

abstract: We introduce a class of integral theorems based on cyclic functions and Riemann sums approximating integrals. The Fourier integral theorem, derived as a combination of a transform and inverse transform, arises as a special case. The integral theorems provide natural estimators of density functions via Monte Carlo methods. Assessments of the quality of the density estimators can be used to obtain optimal cyclic functions, alternatives to the sin function, which minimize square integrals. Our proof techniques rely on a variational approach in ordinary differential equations and the Cauchy residue theorem in complex analysis.


title: Architecture-Guided Test Resource Allocation Via Logic

abstract: We introduce a new logic named Quantitative Confidence Logic (QCL) that quantifies the level of confidence one has in the conclusion of a proof. By translating a fault tree representing a system's architecture to a proof, we show how to use QCL to give a solution to the test resource allocation problem that takes the given architecture into account. We implemented a tool called Astrahl and compared our results to other testing resource allocation strategies.


title: Optimal configurations and "Pauli crystals" of quantum clusters

abstract: Broken rotational and translational symmetries are the hallmarks of solid state materials. In contrast, quantum liquids and gases do not exhibit such properties. However, if we regard the logarithm of the absolute square of a quantum liquid as an energy ${\cal E}= -{\rm ln}|\Psi|^2$, a geometric pattern naturally occurs at the minimum, i.e. the optimal configuration. Such geometric patterns have recently been studied for non-interacting fermions, and have been named "Pauli crystals". However, such patterns exist in all interacting gases (Bose or Fermi), independent of statistics. Here, we present an algorithm to determine the optimal configurations of quantum clusters solely from the images of their densities and without theoretical inputs. We establish its validity by recovering a number of exact results, showing that it can identify the changes in the cluster's ground state which corresponds to phase transitions in bulk systems.


title: Pre-Clustering Point Clouds of Crop Fields Using Scalable Methods

abstract: In order to apply the recent successes of machine learning and automated plant phenotyping on a large scale using agricultural robotics, efficient and general algorithms must be designed to intelligently split crop fields into small, yet actionable, portions that can then be processed by more complex algorithms. In this paper, we notice a similarity between the current state-of-the-art for separating corn plants and a commonly used density-based clustering algorithm, Quickshift. Exploiting this similarity we propose a number of novel, application-specific algorithms with the goal of producing a general and scalable field segmentation algorithm. The novel algorithms proposed in this work are shown to produce quantitatively better results than the current state-of-the-art while being less sensitive to input parameters and maintaining the same algorithmic time complexity. When incorporated into field-scale phenotyping systems, the proposed algorithms should work as a drop-in replacement that can greatly improve the accuracy of results while ensuring that performance and scalability remain undiminished.


title: A shadow of the repulsive Rutherford scattering in the laboratory frame

abstract: The paper explores the Rutherford scattering shadow in an entire class of comoving frames -- inertial frames moving along the initial projectile direction -- of which the laboratory frame, where the target is initially at rest, is a representative example. The paper is a continuation of the previous work addressing the scattering shadow in the fixed-target and the center-of-mass frame. It is shown that the transition from these frames is technically quite involved, due to the scattering shadow forming at an infinite distance from an initial target position. The central procedure involves solving the $5^\mathrm{th}$ degree polynomial as a part of an associated extremization procedure. The shadow existence itself is subject to certain conditions, dependent on a given comoving frame. A new and unexpected phenomenon is found within a certain set of comoving frames, including the laboratory frame itself. It consists of a phase transition between an entirely smooth type of shadow and the one characterized by a formation of a sharp edge.


title: Predicting the redshift of gamma-ray loud AGNs using supervised machine learning

abstract: AGNs are very powerful galaxies characterized by extremely bright emissions coming out from their central massive black holes. Knowing the redshifts of AGNs provides us with an opportunity to determine their distance to investigate important astrophysical problems such as the evolution of the early stars, their formation along with the structure of early galaxies. The redshift determination is challenging because it requires detailed follow-up of multi-wavelength observations, often involving various astronomical facilities. Here, we employ machine learning algorithms to estimate redshifts from the observed gamma-ray properties and photometric data of gamma-ray loud AGN from the Fourth Fermi-LAT Catalog. The prediction is obtained with the Superlearner algorithm, using LASSO selected set of predictors. We obtain a tight correlation, with a Pearson Correlation Coefficient of 71.3% between the inferred and the observed redshifts, an average {\Delta}z_norm = 11.6 x 10^-4. We stress that notwithstanding the small sample of gamma-ray loud AGNs, we obtain a reliable predictive model using Superlearner, which is an ensemble of several machine learning models.


title: Chance-Constrained Motion Planning using Modeled Distance-to-Collision Functions

abstract: This paper introduces Chance Constrained Gaussian Process-Motion Planning (CCGP-MP), a motion planning algorithm for robotic systems under motion and state estimate uncertainties. The paper's key idea is to capture the variations in the distance-to-collision measurements caused by the uncertainty in state estimation techniques using a Gaussian Process (GP) model. We formulate the planning problem as a chance constraint problem and propose a deterministic constraint that uses the modeled distance function to verify the chance-constraints. We apply Simplicial Homology Global Optimization (SHGO) approach to find the global minimum of the deterministic constraint function along the trajectory and use the minimum value to verify the chance-constraints. Under this formulation, we can show that the optimization function is smooth under certain conditions and that SHGO converges to the global minimum. Therefore, CCGP-MP will always guarantee that all points on a planned trajectory satisfy the given chance-constraints. The experiments in this paper show that CCGP-MP can generate paths that reduce collisions and meet optimality criteria under motion and state uncertainties. The implementation of our robot models and path planning algorithm can be found on GitHub.


title: The Rise and Fall of the Eclipsing Binary HS Hydrae

abstract: HS Hydrae is a short period eclipsing binary (P_orb=1.57 day) that belongs to a rare group of systems observed to have rapidly changing inclinations. This evolution is due to a third star on an intermediate orbit, and results in significant differences in eclipse depths and timings year-to-year. Zasche & Paschke (2012) revealed that HS Hydrae's eclipses were rapidly fading from view, predicting they would cease around 2022. Using 25 days of photometric data from Sector 009 of the Transiting Exoplanet Survey Satellite (TESS), we find that the primary eclipses for HS Hydrae were only 0.00173+/-0.00007 mag in depth in March 2019. This data from TESS likely represents the last eclipses detected from HS Hydrae. We also searched the Digitization of the Harvard Astronomical Plate Collection (DASCH) archive for historic data from the system. With a total baseline of over 125 years, this unique combination of data sets - from photographic plates to precision space-based photometry - allows us to trace the emergence and decay of eclipses from HS Hydrae, and further constrain its evolution. Recent TESS observations from Sector 035 confirm that eclipses have ceased for HS Hya, and we estimate they will begin again in 2195.


title: Learning Linear Polytree Structural Equation Models

abstract: We are interested in the problem of learning the directed acyclic graph (DAG) when data are generated from a linear structural equation model (SEM) and the causal structure can be characterized by a polytree. Under the Gaussian polytree models, we study sufficient conditions on the sample sizes for the well-known Chow-Liu algorithm to exactly recover both the skeleton and the equivalence class of the polytree, which is uniquely represented by a CPDAG. On the other hand, necessary conditions on the required sample sizes for both skeleton and CPDAG recovery are also derived in terms of information-theoretic lower bounds, which match the respective sufficient conditions and thereby give a sharp characterization of the difficulty of these tasks. We also consider the problem of inverse correlation matrix estimation under the linear polytree models, and establish the estimation error bound in terms of the dimension and the total number of v-structures. We also consider an extension of group linear polytree models, in which each node represents a group of variables. Our theoretical findings are illustrated by comprehensive numerical simulations, and experiments on benchmark data also demonstrate the robustness of polytree learning when the true graphical structures can only be approximated by polytrees.


title: Reciprocal Multi-Robot Collision Avoidance with Asymmetric State Uncertainty

abstract: We present a general decentralized formulation for a large class of collision avoidance methods and show that all collision avoidance methods of this form are guaranteed to be collision free. This class includes several existing algorithms in the literature as special cases. We then present a particular instance of this collision avoidance method, CARP (Collision Avoidance by Reciprocal Projections), that is effective even when the estimates of other agents' positions and velocities are noisy. The method's main computational step involves the solution of a small convex optimization problem, which can be quickly solved in practice, even on embedded platforms, making it practical to use on computationally-constrained robots such as quadrotors. This method can be extended to find smooth polynomial trajectories for higher dynamic systems such at quadrotors. We demonstrate this algorithm's performance in simulations and on a team of physical quadrotors. Our method finds optimal projections in a median time of 17.12ms for 285 instances of 100 randomly generated obstacles, and produces safe polynomial trajectories at over 60hz on-board quadrotors. Our paper is accompanied by an open source Julia implementation and ROS package.


title: Ego-GNNs: Exploiting Ego Structures in Graph Neural Networks

abstract: Graph neural networks (GNNs) have achieved remarkable success as a framework for deep learning on graph-structured data. However, GNNs are fundamentally limited by their tree-structured inductive bias: the WL-subtree kernel formulation bounds the representational capacity of GNNs, and polynomial-time GNNs are provably incapable of recognizing triangles in a graph. In this work, we propose to augment the GNN message-passing operations with information defined on ego graphs (i.e., the induced subgraph surrounding each node). We term these approaches Ego-GNNs and show that Ego-GNNs are provably more powerful than standard message-passing GNNs. In particular, we show that Ego-GNNs are capable of recognizing closed triangles, which is essential given the prominence of transitivity in real-world graphs. We also motivate our approach from the perspective of graph signal processing as a form of multiplex graph convolution. Experimental results on node classification using synthetic and real data highlight the achievable performance gains using this approach.


title: Random subcomplexes of finite buildings, and fibering of commutator subgroups of right-angled Coxeter groups

abstract: The main theme of this paper is higher virtual algebraic fibering properties of right-angled Coxeter groups (RACGs), with a special focus on those whose defining flag complex is a finite building. We prove for particular classes of finite buildings that their random induced subcomplexes have a number of strong properties, most prominently that they are highly connected. From this we are able to deduce that the commutator subgroup of a RACG, with defining flag complex a finite building of a certain type, admits an epimorphism to $\mathbb{Z}$ whose kernel has strong topological finiteness properties. We additionally use our techniques to present examples where the kernel is of type $\textrm{F}_2$ but not $\textrm{FP}_3$, and examples where the RACG is hyperbolic and the kernel is finitely generated and non-hyperbolic. The key tool we use is a generalization of an approach due to Jankiewicz-Norin-Wise involving Bestvina-Brady discrete Morse theory applied to the Davis complex of a RACG, together with some probabilistic arguments.


title: Inference for High Dimensional Censored Quantile Regression

abstract: With the availability of high dimensional genetic biomarkers, it is of interest to identify heterogeneous effects of these predictors on patients' survival, along with proper statistical inference. Censored quantile regression has emerged as a powerful tool for detecting heterogeneous effects of covariates on survival outcomes. To our knowledge, there is little work available to draw inference on the effects of high dimensional predictors for censored quantile regression. This paper proposes a novel procedure to draw inference on all predictors within the framework of global censored quantile regression, which investigates covariate-response associations over an interval of quantile levels, instead of a few discrete values. The proposed estimator combines a sequence of low dimensional model estimates that are based on multi-sample splittings and variable selection. We show that, under some regularity conditions, the estimator is consistent and asymptotically follows a Gaussian process indexed by the quantile level. Simulation studies indicate that our procedure can properly quantify the uncertainty of the estimates in high dimensional settings. We apply our method to analyze the heterogeneous effects of SNPs residing in lung cancer pathways on patients' survival, using the Boston Lung Cancer Survival Cohort, a cancer epidemiology study on the molecular mechanism of lung cancer.


title: Implicit Rate-Constrained Optimization of Non-decomposable Objectives

abstract: We consider a popular family of constrained optimization problems arising in machine learning that involve optimizing a non-decomposable evaluation metric with a certain thresholded form, while constraining another metric of interest. Examples of such problems include optimizing the false negative rate at a fixed false positive rate, optimizing precision at a fixed recall, optimizing the area under the precision-recall or ROC curves, etc. Our key idea is to formulate a rate-constrained optimization that expresses the threshold parameter as a function of the model parameters via the Implicit Function theorem. We show how the resulting optimization problem can be solved using standard gradient based methods. Experiments on benchmark datasets demonstrate the effectiveness of our proposed method over existing state-of-the art approaches for these problems. The code for the proposed method is available at https://github.com/google-research/google-research/tree/master/implicit_constrained_optimization .


title: Coherent coupling of mechanics to a single nuclear spin

abstract: Nuclear spins interact weakly with their environment. In particular, they are generally insensitive to mechanical vibrations. Here, we successfully demonstrate the coherent coupling of mechanics to a single nuclear spin. This coupling is mediated by a silicon vacancy (SiV) centre in diamond, taking advantage of its large strain susceptibility and hyperfine interaction with nuclear spins. Importantly, we demonstrate that the nuclear spin retains its excellent coherence properties even in the presence of this coupling. This provides a way to leverage nuclear spins as quantum memories for mechanical systems in the quantum regime.


title: Bayesian uncertainty quantification for nuclear matter incompressibility

abstract: Within a Bayesian statistical framework using the standard Skyrme-Hartree-Fock model, the maximum {\it a posteriori} (MAP) values and uncertainties of nuclear matter incompressibility and isovector interaction parameters are inferred from the experimental data of giant resonances and neutron-skin thicknesses of typical heavy nuclei. With the uncertainties of the isovector interaction parameters constrained by the data of the isovector giant dipole resonance and the neutron-skin thickness, we have obtained $K_0 = 223_{-8}^{+7}$ MeV at 68% confidence level using the data of the isoscalar giant monopole resonance in $^{208}$Pb measured at the Research Center for Nuclear Physics (RCNP), Japan, and at the Texas A&M University (TAMU), USA. Although the corresponding $^{120}$Sn data gives a MAP value for $K_0$ about 5 MeV smaller than the $^{208}$Pb data, there are significant overlaps in their posterior probability distribution functions.


title: Compositional Models: Multi-Task Learning and Knowledge Transfer with Modular Networks

abstract: Conditional computation and modular networks have been recently proposed for multitask learning and other problems as a way to decompose problem solving into multiple reusable computational blocks. We propose a new approach for learning modular networks based on the isometric version of ResNet with all residual blocks having the same configuration and the same number of parameters. This architectural choice allows adding, removing and changing the order of residual blocks. In our method, the modules can be invoked repeatedly and allow knowledge transfer to novel tasks by adjusting the order of computation. This allows soft weight sharing between tasks with only a small increase in the number of parameters. We show that our method leads to interpretable self-organization of modules in case of multi-task learning, transfer learning and domain adaptation while achieving competitive results on those tasks. From practical perspective, our approach allows to: (a) reuse existing modules for learning new task by adjusting the computation order, (b) use it for unsupervised multi-source domain adaptation to illustrate that adaptation to unseen data can be achieved by only manipulating the order of pretrained modules, (c) show how our approach can be used to increase accuracy of existing architectures for image classification tasks such as ImageNet, without any parameter increase, by reusing the same block multiple times.


title: Using smartphone photographs of the Moon to acquaint students with non-Euclidean geometry

abstract: Although they are sometimes considered problematic to grasp by students, the concepts behind non-Euclidean geometry can be taught using astronomical images. By using photographs of the Moon taken with a smartphone through a simple telescope, we were able to introduce these concepts to high-school students and college newcomers. By recognizing different Moon geological structures within the photograph, we teach students how to calculate distances of mountain ranges or areas of craters on the Moon's surface, introducing the notions of geodesics and spherical triangles. Furthermore, students can empirically see that the correct estimations for the actual values cannot be obtained using flat geometry. Instead, by using three--dimensional curved geometry, precise estimates of lengths and areas of geological elements in the Moon can be computed with less than 4\% error. These procedures help students understand, concretely, non-Euclidean geometry concepts.


title: Quantum Otto engines at relativistic energies

abstract: Relativistic quantum systems exhibit unique features not present at lower energies, such as the existence of both particles and antiparticles, and restrictions placed on the system dynamics due to the light cone. In order to understand what impact these relativistic phenomena have on the performance of quantum thermal machines we analyze a quantum Otto engine with a working medium of a relativistic particle in an oscillator potential evolving under Dirac or Klein-Gordon dynamics. We examine both the low-temperature, non-relativistic and high-temperature, relativistic limits of the dynamics and find that the relativistic engine operates with higher work output, but an effectively reduced compression ratio, leading to significantly smaller efficiency than its non-relativistic counterpart. Using the framework of endoreversible thermodynamics we determine the efficiency at maximum power of the relativistic engine, and find it to be equivalent to the Curzon-Ahlborn efficiency.


title: Dynamic variability of the phytoplankton electron requirement for carbon fixation in eastern Australian waters

abstract: Fast Repetition Rate fluorometry (FRRf) generates high-resolution measures of phytoplankton primary productivity as electron transport rates (ETRs). How ETRs scale to corresponding inorganic carbon (C) uptake rates (the so-called electron requirement for carbon fixation, e,C), inherently describes the extent and effectiveness with which absorbed light energy drives C-fixation. However, it remains unclear whether and how e,C follows predictable patterns for oceanographic datasets spanning physically dynamic, and complex, environmental gradients. We utilise a unique high-throughput approach, coupling ETRs and 14C-incubations to produce a semi-continuous dataset of e,C (n = 80), predominantly from surface waters, along the Australian coast (Brisbane to the Tasman Sea), including the East Australian Current (EAC). Environmental conditions along this transect could be generally grouped into cooler, more nutrient-rich waters dominated by larger size-fractionated Chl-a (>10 um) versus warmer nutrient-poorer waters dominated by smaller size-fractionated Chl-a (< 2 um). Whilst e,C was higher for warmer water samples, environmental conditions alone explained less than 20% variance of e,C, and changes in predominant size-fraction(s) distributions of Chl-a (biomass) failed to explain variance of e,C. Instead, NPQNSV was a better predictor of e,C, explaining 55% of observed variability. NPQNSV is a physiological descriptor that accounts for changes in both long-term driven acclimation in non-radiative decay, and quasi-instantaneous PSII downregulation, and thus may prove a useful predictor of e,C across physically-dynamic regimes, provided the slope describing their relationship is predictable.


title: Radio Continuum Sources behind the Large Magellanic Cloud

abstract: We present a comprehensive multi-frequency catalogue of radio sources behind the Large Magellanic Cloud between 0.2 and 20 GHz, gathered from a combination of new and legacy radio continuum surveys. This catalogue covers an area of $\sim$144~deg$^2$ at angular resolutions from 45 arcsec to $\sim$3 arcmin. We find 6434 discrete radio sources in total, of which 3789 are detected at two or more radio frequencies. We estimate the median spectral index ($\alpha$; where $S_{v}\sim\nu^\alpha$) of $\alpha = -0.89 $ and mean of $-0.88 \pm 0.48$ for 3636 sources detected exclusively at two frequencies (0.843 and 1.384 GHz) with similar resolution (FWHM $\sim$40-45 arcsec). The large frequency range of the surveys makes it an effective tool to investigate Gigahertz Peak Spectrum (GPS), Compact Steep Spectrum (CSS) and Infrared Faint Radio sources populations within our sample. We find 10 GPS candidates with peak frequencies near 5 GHz, from which we estimate their linear size. 1866 sources from our catalogue are (CSS) candidates with $\alpha <-0.8$. We found six candidates for High Frequency Peaker (HFP) sources, whose radio fluxes peak above 5 GHz and no sources with unconstrained peaks and $\alpha~>0.5$. We found optical counterparts for 343 of the radio continuum sources, of which 128have a redshift measurement. Finally, we investigate the population of 123 Infrared Faint Radio Sources (IFRSs) found in this study.


title: What? Who? Why? Stellify

abstract: In his 1981 article, Roberts highlights the term 'stellify' defined as "to transform (a person or thing) into a star or constellation, to place among the stars." Using the case of the Tabwa people of central Africa, not the Democratic Republic of Congo, Roberts presents among other things the sky as a mnemonic for remembering migrations and remembering culture heroes. We do not know the details of the processes of stellification, however we do know what has been stellified in many cultures by examining their names for stars and asterisms and their skylore. Of the many ideas presented in his latest book, Aveni teases out the ideas of the sky stories having connections to celestial motions, as well as being a mnemonic for remembering seasonal activities and a mnemonic for remembering locally embedded moral, ethical, and sociocultural codes, thus overlapping with Roberts' supposition of the sky serving as a mnemonic. I draw on case studies to flesh out three themes 1. celestial motions, 2. moral, ethical, and sociocultural codes, and 3. seasonal activities within African sky stories. As previously stated, though the human process of assigning names and stories to the night sky as well as stellifying aspects of their lives is not fully understood, these three themes hold promise for being foundational if not part of every culture's practice of stellification.


title: Learning Quadruped Locomotion Policies using Logical Rules

abstract: Quadruped animals are capable of exhibiting a diverse range of locomotion gaits. While progress has been made in demonstrating such gaits on robots, current methods rely on motion priors, dynamics models, or other forms of extensive manual efforts. People can use natural language to describe dance moves. Could one use a formal language to specify quadruped gaits? To this end, we aim to enable easy gait specification and efficient policy learning. Leveraging Reward Machines~(RMs) for high-level gait specification over foot contacts, our approach is called RM-based Locomotion Learning~(RMLL), and supports adjusting gait frequency at execution time. Gait specification is enabled through the use of a few logical rules per gait (e.g., alternate between moving front feet and back feet) and does not require labor-intensive motion priors. Experimental results in simulation highlight the diversity of learned gaits (including two novel gaits), their energy consumption and stability across different terrains, and the superior sample-efficiency when compared to baselines. We also demonstrate these learned policies with a real quadruped robot. Video and supplementary materials: https://sites.google.com/view/rm-locomotion-learning/home


title: The decomposition of the higher-order homology embedding constructed from the $k$-Laplacian

abstract: The null space of the $k$-th order Laplacian $\mathbf{\mathcal L}_k$, known as the {\em $k$-th homology vector space}, encodes the non-trivial topology of a manifold or a network. Understanding the structure of the homology embedding can thus disclose geometric or topological information from the data. The study of the null space embedding of the graph Laplacian $\mathbf{\mathcal L}_0$ has spurred new research and applications, such as spectral clustering algorithms with theoretical guarantees and estimators of the Stochastic Block Model. In this work, we investigate the geometry of the $k$-th homology embedding and focus on cases reminiscent of spectral clustering. Namely, we analyze the {\em connected sum} of manifolds as a perturbation to the direct sum of their homology embeddings. We propose an algorithm to factorize the homology embedding into subspaces corresponding to a manifold's simplest topological components. The proposed framework is applied to the {\em shortest homologous loop detection} problem, a problem known to be NP-hard in general. Our spectral loop detection algorithm scales better than existing methods and is effective on diverse data such as point clouds and images.


title: Adaptively Weighted Top-N Recommendation for Organ Matching

abstract: Reducing the shortage of organ donations to meet the demands of patients on the waiting list has being a major challenge in organ transplantation. Because of the shortage, organ matching decision is the most critical decision to assign the limited viable organs to the most suitable patients. Currently, organ matching decisions were only made by matching scores calculated via scoring models, which are built by the first principles. However, these models may disagree with the actual post-transplantation matching performance (e.g., patient's post-transplant quality of life (QoL) or graft failure measurements). In this paper, we formulate the organ matching decision-making as a top-N recommendation problem and propose an Adaptively Weighted Top-N Recommendation (AWTR) method. AWTR improves performance of the current scoring models by using limited actual matching performance in historical data set as well as the collected covariates from organ donors and patients. AWTR sacrifices the overall recommendation accuracy by emphasizing the recommendation and ranking accuracy for top-N matched patients. The proposed method is validated in a simulation study, where KAS [60] is used to simulate the organ-patient recommendation response. The results show that our proposed method outperforms seven state-of-the-art top-N recommendation benchmark methods.


title: Automatic Construction of Lane-level HD Maps for Urban Scenes

abstract: High definition (HD) maps have demonstrated their essential roles in enabling full autonomy, especially in complex urban scenarios. As a crucial layer of the HD map, lane-level maps are particularly useful: they contain geometrical and topological information for both lanes and intersections. However, large scale construction of HD maps is limited by tedious human labeling and high maintenance costs, especially for urban scenarios with complicated road structures and irregular markings. This paper proposes an approach based on semantic-particle filter to tackle the automatic lane-level mapping problem in urban scenes. The map skeleton is firstly structured as a directed cyclic graph from online mapping database OpenStreetMap. Our proposed method then performs semantic segmentation on 2D front-view images from ego vehicles and explores the lane semantics on a birds-eye-view domain with true topographical projection. Exploiting OpenStreetMap, we further infer lane topology and reference trajectory at intersections with the aforementioned lane semantics. The proposed algorithm has been tested in densely urbanized areas, and the results demonstrate accurate and robust reconstruction of the lane-level HD map.


title: Coexistence of topological and nontopological Fermi-superfluid phases

abstract: The two-dimensional spin-imbalanced Fermi gas subject to s-wave pairing and spin-orbit coupling is considered a promising platform for realizing a topological chiral-p-wave superfluid. In the BCS limit of s-wave pairing, i.e., when Cooper pairs are only weakly bound, the system enters the topological phase via a second-order transition driven by increasing the Zeeman spin-splitting energy. Stronger attractive two-particle interactions cause the system to undergo the BCS-BEC crossover, in the course of which the topological transition becomes first-order. As a result, topological and nontopological superfluids coexist in spatially separated domains in an extended region of phase space spanned by the strength of s-wave interactions and the Zeeman energy. Here we investigate this phase-coexistence region theoretically using a zero-temperature mean-field approach. Exact numerical results are presented to illustrate basic physical characteristics of the coexisting phases and to validate an approximate analytical description derived for weak spin-orbit coupling. Besides extending our current understanding of spin-imbalanced superfluid Fermi systems, the present approach also provides a platform for future studies of unconventional Majorana excitations that, according to topology, should be present at the internal interface between coexisting topological and nontopological superfluid parts of the system.


title: A note on sharp oracle bounds for Slope and Lasso

abstract: In this paper, we study the sharp oracle bounds for Slope and Lasso and generalize the results in Bellec et al. (2018) to allow the case that the parameter vector is not exactly sparse and obtain the optimal bounds for $\ell_q$ estimation errors with $1\leq q\leq \infty$ by using some extended Restricted Eigenvalue type conditions.


title: Call for a Framework for Reporting Evidence for Life Beyond Earth

abstract: Ours could realistically be the generation to discover evidence of life beyond Earth. With this privileged potential comes responsibility. The magnitude of the question, "are we alone?", and the public interest therein, opens the possibility that results may be taken to imply more than the observations support, or than the observers intend. As life detection objectives become increasingly prominent in space sciences, it is essential to open a community dialog about how to convey information in a subject matter that is diverse, complicated, and has high potential to be sensationalized. Establishing best practices for communicating about life detection can serve to set reasonable expectations on the early stages of a hugely challenging endeavor, attach value to incremental steps along the path, and build public trust by making clear that "false starts" and "dead ends" are an expected and potentially productive part of the scientific process. Here, we endeavor to motivate and seed the discussion with basic considerations and offer an example of how such considerations might be incorporated and applied in a proof-of-concept-level framework. Everything mentioned herein, including the name of the confidence scale, is intended not as a prescription, but simply as the beginning of an important dialogue.


title: Federated Learning Versus Classical Machine Learning: A Convergence Comparison

abstract: In the past few decades, machine learning has revolutionized data processing for large scale applications. Simultaneously, increasing privacy threats in trending applications led to the redesign of classical data training models. In particular, classical machine learning involves centralized data training, where the data is gathered, and the entire training process executes at the central server. Despite significant convergence, this training involves several privacy threats on participants' data when shared with the central cloud server. To this end, federated learning has achieved significant importance over distributed data training. In particular, the federated learning allows participants to collaboratively train the local models on local data without revealing their sensitive information to the central cloud server. In this paper, we perform a convergence comparison between classical machine learning and federated learning on two publicly available datasets, namely, logistic-regression-MNIST dataset and image-classification-CIFAR-10 dataset. The simulation results demonstrate that federated learning achieves higher convergence within limited communication rounds while maintaining participants' anonymity. We hope that this research will show the benefits and help federated learning to be implemented widely.


title: Tsformer: Time series Transformer for tourism demand forecasting

abstract: AI-based methods have been widely applied to tourism demand forecasting. However, current AI-based methods are short of the ability to process long-term dependency, and most of them lack interpretability. The Transformer used initially for machine translation shows an incredible ability to long-term dependency processing. Based on the Transformer, we proposed a time series Transformer (Tsformer) with Encoder-Decoder architecture for tourism demand forecasting. The proposed Tsformer encodes long-term dependency with encoder, captures short-term dependency with decoder, and simplifies the attention interactions under the premise of highlighting dominant attention through a series of attention masking mechanisms. These improvements make the multi-head attention mechanism process the input sequence according to the time relationship, contributing to better interpretability. What's more, the context processing ability of the Encoder-Decoder architecture allows adopting the calendar of days to be forecasted to enhance the forecasting performance. Experiments conducted on the Jiuzhaigou valley and Siguniang mountain tourism demand datasets with other nine baseline methods indicate that the proposed Tsformer outperformed all baseline models in the short-term and long-term tourism demand forecasting tasks. Moreover, ablation studies demonstrate that the adoption of the calendar of days to be forecasted contributes to the forecasting performance of the proposed Tsformer. For better interpretability, the attention weight matrix visualization is performed. It indicates that the Tsformer concentrates on seasonal features and days close to days to be forecast in short-term forecasting.


title: Kurtosis of von Neumann entanglement entropy

abstract: In this work, we study the statistical behavior of entanglement in quantum bipartite systems under the Hilbert-Schmidt ensemble as assessed by the standard measure - the von Neumann entropy. Expressions of the first three exact cumulants of von Neumann entropy are known in the literature. The main contribution of the present work is the exact formula of the corresponding fourth cumulant that controls the tail behavior of the distribution. As a key ingredient in deriving the result, we make use of newly observed unsimplifiable summation bases that lead to a complete cancellation. In addition to providing further evidence of the conjectured Gaussian limit of the von Neumann entropy, the obtained formula also provides an improved finite-size approximation to the distribution.


title: Rectifying Administrated ERC20 Tokens

abstract: The developers of Ethereum smart contracts often implement administrating patterns, such as censoring certain users, creating or destroying balances on demand, destroying smart contracts, or injecting arbitrary code. These routines turn an ERC20 token into an administrated token - the type of Ethereum smart contract that we scrutinize in this research. We discover that many smart contracts are administrated, and the owners of these tokens carry lesser social and legal responsibilities compared to the traditional centralized actors that those tokens intend to disrupt. This entails two major problems: a) the owners of the tokens have the ability to quickly steal all the funds and disappear from the market; and b) if the private key of the owner's account is stolen, all the assets might immediately turn into the property of the attacker. We develop a pattern recognition framework based on 9 syntactic features characterizing administrated ERC20 tokens, which we use to analyze existing smart contracts deployed on Ethereum Mainnet. Our analysis of 84,062 unique Ethereum smart contracts reveals that nearly 58% of them are administrated ERC20 tokens, which accounts for almost 90% of all ERC20 tokens deployed on Ethereum. To protect users from the frivolousness of unregulated token owners without depriving the ability of these owners to properly manage their tokens, we introduce SafelyAdministrated - a library that enforces a responsible ownership and management of ERC20 tokens. The library introduces three mechanisms: deferred maintenance, board of trustees and safe pause. We implement and test SafelyAdministrated in the form of Solidity abstract contract, which is ready to be used by the next generation of safely administrated ERC20 tokens.


title: Economic Recession Prediction Using Deep Neural Network

abstract: We investigate the effectiveness of different machine learning methodologies in predicting economic cycles. We identify the deep learning methodology of Bi-LSTM with Autoencoder as the most accurate model to forecast the beginning and end of economic recessions in the U.S. We adopt commonly-available macro and market-condition features to compare the ability of different machine learning models to generate good predictions both in-sample and out-of-sample. The proposed model is flexible and dynamic when both predictive variables and model coefficients vary over time. It provided good out-of-sample predictions for the past two recessions and early warning about the COVID-19 recession.


title: Score-Based Point Cloud Denoising

abstract: Point clouds acquired from scanning devices are often perturbed by noise, which affects downstream tasks such as surface reconstruction and analysis. The distribution of a noisy point cloud can be viewed as the distribution of a set of noise-free samples $p(x)$ convolved with some noise model $n$, leading to $(p * n)(x)$ whose mode is the underlying clean surface. To denoise a noisy point cloud, we propose to increase the log-likelihood of each point from $p * n$ via gradient ascent -- iteratively updating each point's position. Since $p * n$ is unknown at test-time, and we only need the score (i.e., the gradient of the log-probability function) to perform gradient ascent, we propose a neural network architecture to estimate the score of $p * n$ given only noisy point clouds as input. We derive objective functions for training the network and develop a denoising algorithm leveraging on the estimated scores. Experiments demonstrate that the proposed model outperforms state-of-the-art methods under a variety of noise models, and shows the potential to be applied in other tasks such as point cloud upsampling. The code is available at \url{https://github.com/luost26/score-denoise}.


title: Triangular Matrix Categories over path Categories and Quasi-hereditary Categories, as well as one point extensions by Projectives

abstract: In this paper, we prove that the lower triangular matrix category $\Lambda =\left [ \begin{smallmatrix} \mathcal{T}&0\\ M&\mathcal{U} \end{smallmatrix} \right ]$, where $\mathcal{T}$ and $\mathcal{U}$ are quasi-hereditary $\mathrm{Hom}$-finite Krull-Schmidt $K$-categories and $M$ is a $\mathcal U\otimes_K \mathcal T^{op}$-module that satisfies suitable conditions, is quasi-hereditary in the sense of \cite{LGOS1} and \cite{Martin}. Moreover, we solve the problem of finding quotients of path categories isomorphic to the lower triangular matrix category $\Lambda$, where $\mathcal T=K\mathcal{R/J}$ and $\mathcal U=K\mathcal{Q/I}$ are path categories of infinity quivers modulo admissible ideals. Finally, we study the case where $\Lambda$ is a path category of a quiver $Q$ with relations and $\mathcal U$ is the full additive subcategory of $\Lambda$ obtained by deleting a source vertex $*$ in $Q$ and $\mathcal T=\mathrm{add} \{*\}$. We then show that there exists an adjoint pair of functors $(\mathcal R, \mathcal E)$ between the functor categories $\mathrm{mod} \ \Lambda$ and $\mathrm{mod} \ \mathcal U$ that preserve orthogonality and exceptionality; see \cite{Assem1}. We then give some examples of how to extend classical tilting subcategories of $\mathcal U$-modules to classical tilting subcategories of $\Lambda$-modules.


title: Spontaneous chiral symmetry breaking in holographic soft wall models

abstract: We investigate non-linear extensions of the holographic soft wall model proposed by Karch, Katz, Son and Stephanov [1] including non-minimal couplings in the five-dimensional action. The non-minimal couplings bring a new parameter $a_0$ which controls the transition between spontaneous and explicit symmetry breaking near the limit of massless quarks (the chiral limit). In the physical region (positive quark mass), we show that above a critical value of the parameter $a_0$ the chiral condensate $\langle \bar{q} q \rangle$ is finite in the chiral limit, signifying spontaneous chiral symmetry breaking. This result is supported by the lightest states arising in the spectrum of the pseudoscalar mesons, which become massless in the chiral limit and are therefore intrepreted as Nambu-Goldstone bosons. Moreover, the decay constants of the pseudoscalar mesons also support this conclusion, as well as the Gell-Mann-Oakes-Renner (GOR) relation satisfied by the lightest states. We also calculate the spectrum of scalar, vector, and axial-vector mesons with their corresponding decay constants. We describe the evolution of masses and decay constants with the increasing of the quark mass and for the physical mass we compare our results against available experimental data. Finally, we do not find instabilities in our model for the physical region (positive quark mass).


title: Human Pose Transfer with Augmented Disentangled Feature Consistency

abstract: Deep generative models have made great progress in synthesizing images with arbitrary human poses and transferring poses of one person to others. Though many different methods have been proposed to generate images with high visual fidelity, the main challenge remains and comes from two fundamental issues: pose ambiguity and appearance inconsistency. To alleviate the current limitations and improve the quality of the synthesized images, we propose a pose transfer network with augmented Disentangled Feature Consistency (DFC-Net) to facilitate human pose transfer. Given a pair of images containing the source and target person, DFC-Net extracts pose and static information from the source and target respectively, then synthesizes an image of the target person with the desired pose from the source. Moreover, DFC-Net leverages disentangled feature consistency losses in the adversarial training to strengthen the transfer coherence and integrates a keypoint amplifier to enhance the pose feature extraction. With the help of the disentangled feature consistency losses, we further propose a novel data augmentation scheme that introduces unpaired support data with the augmented consistency constraints to improve the generality and robustness of DFC-Net. Extensive experimental results on Mixamo-Pose and EDN-10k have demonstrated DFC-Net achieves state-of-the-art performance on pose transfer.


title: A collection of results relating the geometry of plane domains and the exit time of planar Brownian motion

abstract: We prove a number of results relating exit times of planar Brownian with the geometric properties of the domains in question. Included are proofs of the conformal invariance of moduli of rectangles and annuli using Brownian motion; similarly probabilistic proofs of some recent results of Karafyllia on harmonic measure on starlike domains; examples of domains and their complements which are simultaneously large when measured by the moments of exit time of Brownian motion, and examples of domains and their complements which are simultaneously small; and proofs of several identities involving the Cauchy distribution using the optional stopping theorem.


title: Lower Bounds for Symmetric Circuits for the Determinant

abstract: Dawar and Wilsenach (ICALP 2020) introduce the model of symmetric arithmetic circuits and show an exponential separation between the sizes of symmetric circuits for computing the determinant and the permanent. The symmetry restriction is that the circuits which take a matrix input are unchanged by a permutation applied simultaneously to the rows and columns of the matrix. Under such restrictions we have polynomial-size circuits for computing the determinant but no subexponential size circuits for the permanent. Here, we consider a more stringent symmetry requirement, namely that the circuits are unchanged by arbitrary even permutations applied separately to rows and columns, and prove an exponential lower bound even for circuits computing the determinant. The result requires substantial new machinery. We develop a general framework for proving lower bounds for symmetric circuits with restricted symmetries, based on a new support theorem and new two-player restricted bijection games. These are applied to the determinant problem with a novel construction of matrices that are bi-adjacency matrices of graphs based on the CFI construction. Our general framework opens the way to exploring a variety of symmetry restrictions and studying trade-offs between symmetry and other resources used by arithmetic circuits.


title: Octo-Tiger's New Hydro Module and Performance Using HPX+CUDA on ORNL's Summit

abstract: Octo-Tiger is a code for modeling three-dimensional self-gravitating astrophysical fluids. It was particularly designed for the study of dynamical mass transfer between interacting binary stars. Octo-Tiger is parallelized for distributed systems using the asynchronous many-task runtime system, the C++ standard library for parallelism and concurrency (HPX) and utilizes CUDA for its gravity solver. Recently, we have remodeled Octo-Tiger's hydro solver to use a three-dimensional reconstruction scheme. In addition, we have ported the hydro solver to GPU using CUDA kernels. We present scaling results for the new hydro kernels on ORNL's Summit machine using a Sedov-Taylor blast wave problem. We also compare Octo-Tiger's new hydro scheme with its old hydro scheme, using a rotating star as a test problem.


title: Formalizing Galois Theory

abstract: We describe a project to formalize Galois theory using the Lean theorem prover, which is part of a larger effort to formalize all of the standard undergraduate mathematics curriculum in Lean. We discuss some of the challenges we faced and the decisions we made in the course of this project. The main theorems we formalized are the primitive element theorem, the fundamental theorem of Galois theory, and the equivalence of several characterizations of finite degree Galois extensions.


title: Estimating Predictive Uncertainty Under Program Data Distribution Shift

abstract: Deep learning (DL) techniques have achieved great success in predictive accuracy in a variety of tasks, but deep neural networks (DNNs) are shown to produce highly overconfident scores for even abnormal samples. Well-defined uncertainty indicates whether a model's output should (or should not) be trusted and thus becomes critical in real-world scenarios which typically involves shifted input distributions due to many factors. Existing uncertainty approaches assume that testing samples from a different data distribution would induce unreliable model predictions thus have higher uncertainty scores. They quantify model uncertainty by calibrating DL model's confidence of a given input and evaluate the effectiveness in computer vision (CV) and natural language processing (NLP)-related tasks. However, their methodologies' reliability may be compromised under programming tasks due to difference in data representations and shift patterns. In this paper, we first define three different types of distribution shift in program data and build a large-scale shifted Java dataset. We implement two common programming language tasks on our dataset to study the effect of each distribution shift on DL model performance. We also propose a large-scale benchmark of existing state-of-the-art predictive uncertainty on programming tasks and investigate their effectiveness under data distribution shift. Experiments show that program distribution shift does degrade the DL model performance to varying degrees and that existing uncertainty methods all present certain limitations in quantifying uncertainty on program dataset.


title: Detail Preserving Residual Feature Pyramid Modules for Optical Flow

abstract: Feature pyramids and iterative refinement have recently led to great progress in optical flow estimation. However, downsampling in feature pyramids can cause blending of foreground objects with the background, which will mislead subsequent decisions in the iterative processing. The results are missing details especially in the flow of thin and of small structures. We propose a novel Residual Feature Pyramid Module (RFPM) which retains important details in the feature map without changing the overall iterative refinement design of the optical flow estimation. RFPM incorporates a residual structure between multiple feature pyramids into a downsampling module that corrects the blending of objects across boundaries. We demonstrate how to integrate our module with two state-of-the-art iterative refinement architectures. Results show that our RFPM visibly reduces flow errors and improves state-of-art performance in the clean pass of Sintel, and is one of the top-performing methods in KITTI. According to the particular modular structure of RFPM, we introduce a special transfer learning approach that can dramatically decrease the training time compared to a typical full optical flow training schedule on multiple datasets.


title: A novel meta-learning initialization method for physics-informed neural networks

abstract: Physics-informed neural networks (PINNs) have been widely used to solve various scientific computing problems. However, large training costs limit PINNs for some real-time applications. Although some works have been proposed to improve the training efficiency of PINNs, few consider the influence of initialization. To this end, we propose a New Reptile initialization based Physics-Informed Neural Network (NRPINN). The original Reptile algorithm is a meta-learning initialization method based on labeled data. PINNs can be trained with less labeled data or even without any labeled data by adding partial differential equations (PDEs) as a penalty term into the loss function. Inspired by this idea, we propose the new Reptile initialization to sample more tasks from the parameterized PDEs and adapt the penalty term of the loss. The new Reptile initialization can acquire initialization parameters from related tasks by supervised, unsupervised, and semi-supervised learning. Then, PINNs with initialization parameters can efficiently solve PDEs. Besides, the new Reptile initialization can also be used for the variants of PINNs. Finally, we demonstrate and verify the NRPINN considering both forward problems, including solving Poisson, Burgers, and Schr\"odinger equations, as well as inverse problems, where unknown parameters in the PDEs are estimated. Experimental results show that the NRPINN training is much faster and achieves higher accuracy than PINNs with other initialization methods.


title: Design of the Propulsion System of Nano satellite: StudSat2

abstract: The increase in the application of the satellite has skyrocketed the number of satellites, especially in the low earth orbit. The major concern today is after the end of life, these satellites become debris which negatively affects the space environment. As per the international guidelines of the European Space Agency, it is mandatory to deorbit the satellite within 25 years of the end of life. StudSat1, which was successfully launched on 12th July 2010, is the first Pico satellite developed in India by undergraduate students from seven different engineering colleges across South India. Now, the team is developing StudSat2, which is India's first twin satellite mission having two nanosatellites whose overall mass is less than 10kg. This paper is aimed to design the propulsion system, cold gas thruster, to deorbit StudSat2 from its original orbit i.e. 600 km to lower orbit i.e. 400km. The propulsion system mainly consists of a storage tank, pipes, Convergent Divergent nozzle, and electronic actuators. The paper also gives information about the components of cold gas thruster, which have been designed in the CATIA V5, and the structural and flow analysis of the same has been done in ANSYS. The concept of Hohmann transfer has been used to deorbit the satellite and STK has been used to simulate it.


title: A 60-GHz Radar Sensor for Micron-Scale Motion Detection

abstract: A compact, continuous-wave, mmWave radar sensor is developed for non-contact detection of micron-scale motions. This board-integrated radar system consists of a pair of mmWave transmitter and receiver, two series-fed microstrip patch arrays, an IF subsystem, and a microcontroller. Working at 60-GHz frequency, this super-heterodyne, digital-IF radar sensor exhibits an agile sensitivity and a robust anti-noise performance. Assisted by a gradient-descent DC offset estimation and a phase demodulation algorithm based on extended differential and cross-multiply, a 45-{\mu}m pendulum swing can be experimentally detected with a 20.6-dB SNR, verifying its strong ability in Doppler micro-motion detections.


title: High-Level Coupled-Cluster Energetics by Merging Moment Expansions with Selected Configuration Interaction

abstract: Inspired by our earlier semi-stochastic work aimed at converging high-level coupled-cluster (CC) energetics [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017); J. Chem. Phys. 154, 124103 (2021)], we propose a novel form of the CC($P$;$Q$) theory in which the stochastic Quantum Monte Carlo propagations, used to identify dominant higher-than-doubly excited determinants, are replaced by the selected configuration interaction (CI) approach using the CIPSI algorithm. The advantages of the resulting CIPSI-driven CC($P$;$Q$) methodology are illustrated by a few molecular examples, including the dissociation of $\mathrm{F_2}$ and the automerization of cyclobutadiene, where we recover the electronic energies corresponding to the CC calculations with a full treatment of singles, doubles, and triples based on the information extracted from compact CI wave functions originating from relatively inexpensive Hamiltonian diagonalizations.


title: Weather on Other Worlds. VI. Optical Spectrophotometry of Luhman 16B Reveals Large-amplitude Variations in the Alkali Lines

abstract: Using a novel wide-slit, multi-object approach with the GMOS spectrograph on the 8-meter Gemini South telescope, we have obtained precise time-series spectrophotometry of the binary brown dwarf Luhman 16 at optical wavelengths over two full nights. The B component of this binary system is known to be variable in the red optical and near-infrared with a period of 5 hr and an amplitude of 5--20%. Our observations probe its spectrally-resolved variability in the 6000--10000 Angstrom range. At wavelengths affected by the extremely strong, broadened spectral lines of the neutral alkali metals (the potassium doublet centered near 7682 Angstroms and the sodium doublet at 5893 Angstroms), we see photometric variations that differ strikingly from the those of the 8000--10000 Angstrom `red continuum' that dominates our detected flux. On UT 2014 February 24, these variations are anticorrelated with the red continuum, while on Feb 25 they have a large relative phase shift. The extent to which the wavelength-dependent photometric behavior diverges from that of the red continuum appears to correlate with the strength of the alkali absorption. We consider but ultimately reject models in which our observations are explained by lightning or auroral activity. A more likely cause is cloud-correlated, altitude-dependent variations in the gas-phase abundances of sodium and potassium, which are in chemical equilibrium with their chlorides in brown dwarf atmospheres. Clouds could influence these chemical equilibria by changing the atmospheric temperature profile and/or through cloud particles acting as chemical catalysts.


title: Communication Efficiency in Federated Learning: Achievements and Challenges

abstract: Federated Learning (FL) is known to perform Machine Learning tasks in a distributed manner. Over the years, this has become an emerging technology especially with various data protection and privacy policies being imposed FL allows performing machine learning tasks whilst adhering to these challenges. As with the emerging of any new technology, there are going to be challenges and benefits. A challenge that exists in FL is the communication costs, as FL takes place in a distributed environment where devices connected over the network have to constantly share their updates this can create a communication bottleneck. In this paper, we present a survey of the research that is performed to overcome the communication constraints in an FL setting.


title: Resource Efficient Mountainous Skyline Extraction using Shallow Learning

abstract: Skyline plays a pivotal role in mountainous visual geo-localization and localization/navigation of planetary rovers/UAVs and virtual/augmented reality applications. We present a novel mountainous skyline detection approach where we adapt a shallow learning approach to learn a set of filters to discriminate between edges belonging to sky-mountain boundary and others coming from different regions. Unlike earlier approaches, which either rely on extraction of explicit feature descriptors and their classification, or fine-tuning general scene parsing deep networks for sky segmentation, our approach learns linear filters based on local structure analysis. At test time, for every candidate edge pixel, a single filter is chosen from the set of learned filters based on pixel's structure tensor, and then applied to the patch around it. We then employ dynamic programming to solve the shortest path problem for the resultant multistage graph to get the sky-mountain boundary. The proposed approach is computationally faster than earlier methods while providing comparable performance and is more suitable for resource constrained platforms e.g., mobile devices, planetary rovers and UAVs. We compare our proposed approach against earlier skyline detection methods using four different data sets. Our code is available at \url{https://github.com/TouqeerAhmad/skyline_detection}.


title: Pruning Ternary Quantization

abstract: Inference time, model size, and accuracy are three key factors in deep model compression. Most of the existing work addresses these three key factors separately as it is difficult to optimize them all at the same time. For example, low-bit quantization aims at obtaining a faster model; weight sharing quantization aims at improving compression ratio and accuracy; and mixed-precision quantization aims at balancing accuracy and inference time. To simultaneously optimize bit-width, model size, and accuracy, we propose pruning ternary quantization (PTQ): a simple, effective, symmetric ternary quantization method. We integrate L2 normalization, pruning, and the weight decay term to reduce the weight discrepancy in the gradient estimator during quantization, thus producing highly compressed ternary weights. Our method brings the highest test accuracy and the highest compression ratio. For example, it produces a 939kb (49$\times$) 2bit ternary ResNet-18 model with only 4\% accuracy drop on the ImageNet dataset. It compresses 170MB Mask R-CNN to 5MB (34$\times$) with only 2.8\% average precision drop. Our method is verified on image classification, object detection/segmentation tasks with different network structures such as ResNet-18, ResNet-50, and MobileNetV2.


title: Toroidal compactifications and Borel--Serre compactifications

abstract: We discuss connections of toroidal compactifications and Borel--Serre compactifications in view of the fundamental diagram of extended period domains. We give a complement to a work of Goresky--Tai.


title: Generalized Bloch model: a theory for pulsed magnetization transfer

abstract: Purpose: The paper introduces a classical model to describe the dynamics of large spin-1/2 ensembles associated with nuclei bound in large molecule structures, commonly referred to as the semi-solid spin pool, and their magnetization transfer (MT) to spins of nuclei in Theory and Methods: Like quantum-mechanical descriptions of spin dynamics and like the original Bloch equations, but unlike existing MT models, the proposed model is based on the algebra of angular momentum in the sense that it explicitly models the rotations induced by radio-frequency (RF) pulses. It generalizes the original Bloch model to non-exponential decays, which are, e.g., observed for semi-solid spin pools. The combination of rotations with non-exponential decays is facilitated by describing the latter as Green's functions, comprised in an integro-differential equation. Results: Our model describes the data of an inversion-recovery magnetization-transfer experiment with varying durations of the inversion pulse substantially better than established models. We made this observation for all measured data, but in particular for pulse durations small than 300$\mu$s. Furthermore, we provide a linear approximation of the generalized Bloch model that reduces the simulation time by approximately a factor 15,000, enabling simulation of the spin dynamics caused by a rectangular RF-pulse in roughly 2$\mu$s. Conclusion: The proposed theory unifies the original Bloch model, Henkelman's steady-state theory for magnetization transfer, and the commonly assumed rotation induced by hard pulses (i.e., strong and infinitesimally short applications of RF fields) and describes experimental data better than previous models.


title: Photon-Starved Scene Inference using Single Photon Cameras

abstract: Scene understanding under low-light conditions is a challenging problem. This is due to the small number of photons captured by the camera and the resulting low signal-to-noise ratio (SNR). Single-photon cameras (SPCs) are an emerging sensing modality that are capable of capturing images with high sensitivity. Despite having minimal read-noise, images captured by SPCs in photon-starved conditions still suffer from strong shot noise, preventing reliable scene inference. We propose photon scale-space a collection of high-SNR images spanning a wide range of photons-per-pixel (PPP) levels (but same scene content) as guides to train inference model on low photon flux images. We develop training techniques that push images with different illumination levels closer to each other in feature representation space. The key idea is that having a spectrum of different brightness levels during training enables effective guidance, and increases robustness to shot noise even in extreme noise cases. Based on the proposed approach, we demonstrate, via simulations and real experiments with a SPAD camera, high-performance on various inference tasks such as image classification and monocular depth estimation under ultra low-light, down to < 1 PPP.


title: Rate-Optimal Robust Estimation of High-Dimensional Vector Autoregressive Models

abstract: High-dimensional time series data appear in many scientific areas in the current data-rich environment. Analysis of such data poses new challenges to data analysts because of not only the complicated dynamic dependence between the series, but also the existence of aberrant observations, such as missing values, contaminated observations, and heavy-tailed distributions. For high-dimensional vector autoregressive (VAR) models, we introduce a unified estimation procedure that is robust to model misspecification, heavy-tailed noise contamination, and conditional heteroscedasticity. The proposed methodology enjoys both statistical optimality and computational efficiency, and can handle many popular high-dimensional models, such as sparse, reduced-rank, banded, and network-structured VAR models. With proper regularization and data truncation, the estimation convergence rates are shown to be almost optimal in the minimax sense under a bounded $(2+2\epsilon)$-th moment condition. When $\epsilon\geq1$, the rates of convergence match those obtained under the sub-Gaussian assumption. Consistency of the proposed estimators is also established for some $\epsilon\in(0,1)$, with minimax optimal convergence rates associated with $\epsilon$. The efficacy of the proposed estimation methods is demonstrated by simulation and a U.S. macroeconomic example.


title: Model Selection for Offline Reinforcement Learning: Practical Considerations for Healthcare Settings

abstract: Reinforcement learning (RL) can be used to learn treatment policies and aid decision making in healthcare. However, given the need for generalization over complex state/action spaces, the incorporation of function approximators (e.g., deep neural networks) requires model selection to reduce overfitting and improve policy performance at deployment. Yet a standard validation pipeline for model selection requires running a learned policy in the actual environment, which is often infeasible in a healthcare setting. In this work, we investigate a model selection pipeline for offline RL that relies on off-policy evaluation (OPE) as a proxy for validation performance. We present an in-depth analysis of popular OPE methods, highlighting the additional hyperparameters and computational requirements (fitting/inference of auxiliary models) when used to rank a set of candidate policies. We compare the utility of different OPE methods as part of the model selection pipeline in the context of learning to treat patients with sepsis. Among all the OPE methods we considered, fitted Q evaluation (FQE) consistently leads to the best validation ranking, but at a high computational cost. To balance this trade-off between accuracy of ranking and computational efficiency, we propose a simple two-stage approach to accelerate model selection by avoiding potentially unnecessary computation. Our work serves as a practical guide for offline RL model selection and can help RL practitioners select policies using real-world datasets. To facilitate reproducibility and future extensions, the code accompanying this paper is available online at https://github.com/MLD3/OfflineRL_ModelSelection.


title: Domain Adaptive Video Segmentation via Temporal Consistency Regularization

abstract: Video semantic segmentation is an essential task for the analysis and understanding of videos. Recent efforts largely focus on supervised video segmentation by learning from fully annotated data, but the learnt models often experience clear performance drop while applied to videos of a different domain. This paper presents DA-VSN, a domain adaptive video segmentation network that addresses domain gaps in videos by temporal consistency regularization (TCR) for consecutive frames of target-domain videos. DA-VSN consists of two novel and complementary designs. The first is cross-domain TCR that guides the prediction of target frames to have similar temporal consistency as that of source frames (learnt from annotated source data) via adversarial learning. The second is intra-domain TCR that guides unconfident predictions of target frames to have similar temporal consistency as confident predictions of target frames. Extensive experiments demonstrate the superiority of our proposed domain adaptive video segmentation network which outperforms multiple baselines consistently by large margins.


title: SU(2) representations and a large surgery formula

abstract: A knot $K\subset S^3$ is called $SU(2)$-abundant if for all but finitely many $r\in\mathbb{Q}\backslash\{0\}$, there is an irreducible representation $\pi_1(S^3_r(K))\to SU(2)$, and the slope $r=u/v\neq 0$ with no irreducible $SU(2)$ representation must satisfy $\Delta_K(\zeta^2)= 0$ for some $u$-th root of unity $\zeta$. We prove that a nontrivial knot $K\subset S^3$ is $SU(2)$-abundant unless it is a prime knot and the coefficients of its Alexander polynomial $\Delta_K(t)$ lie in $\{-1,0,1\}$. In particular, any hyperbolic alternating knot is $SU(2)$-abundant. The proof is based on a large surgery formula that relates instanton knot homology $KHI(S^3,K)$ and the framed instanton homology $I^\sharp(S^3_n(K))$ for any integer $n$ satisfies $|n|\ge 2g(K)+1$. By the same technique, we can calculate many examples of instanton Floer homology. First, for any Berge knot $K$, the spaces $KHI(S^3,K)$ and $\widehat{HFK}(S^3,K)$ have the same dimension. Second, for any dual knot $K_r\subset S^3_r(K)$ of a Berge knot $K$ with $r> 2g(K)-1$, we show $\dim_\mathbb{C}KHI(S^3_r(K),K_r)=|H_1(S^3_r(K);\mathbb{Z})|$. Third, for any genus-one alternating knot $K$ and any $r\in\mathbb{Q}\backslash\{0\}$, the spaces $I^\sharp(S^3_r(K))$ and $\widehat{HF}(S_r^3(K))$ have the same dimension.


title: The Fast Evolving, Tremendous and Blue Superoutburst in ASASSN-21au Reveals a Dichotomy in the Outbursts of Long-period AM CVns

abstract: ASASSN-21au is an ultracompact accreting white dwarf binary (AM CVn) with a period of approx. 58 min. Using multiwavelength observations of the system, we discovered a dichotomy in the behavior of outbursts in AM CVns. The binary showed an initial brightness increase which lasted for at least 82 days, followed by an additional increase which lasted 2 weeks. Afterwards ASASSN-21au went into superoutburst with a total duration of 19 days, showing an amplitude with respect to quiescence of approx. 7.5 mags in g, with a precursor and an echo outburst. A correlation between X-rays, UV and optical was identified for the first time in an AM CVn during this stage. The color evolution of ASASSN-21au indicates that during the superoutburst the dominant component was the accretion disk. The short duration, large amplitude and color evolution of the superoutburst agree with expectations from the disk instability model. These characteristics are opposite to the ones observed in SDSS J080710+485259 and SDSS J113732+405458, which have periods of approx. 53 min and 60 min, respectively. The initially slow brightness increase in the light curve of ASASSN-21au and the behavior after the superoutburst favors a scenario in which changes in the mass-transfer rate led to disk instabilities, while the outburst mechanism of SDSS J080710+485259 and SDSS J113732+405458 has been attributed to enhanced mass-transfer alone. Further observations are needed to understand the origin of this dichotomy.


title: Dynamic Proximal Unrolling Network for Compressive Imaging

abstract: Compressive imaging aims to recover a latent image from under-sampled measurements, suffering from a serious ill-posed inverse problem. Recently, deep neural networks have been applied to this problem with superior results, owing to the learned advanced image priors. These approaches, however, require training separate models for different imaging modalities and sampling ratios, leading to overfitting to specific settings. In this paper, a dynamic proximal unrolling network (dubbed DPUNet) was proposed, which can handle a variety of measurement matrices via one single model without retraining. Specifically, DPUNet can exploit both the embedded observation model via gradient descent and imposed image priors by learned dynamic proximal operators, achieving joint reconstruction. A key component of DPUNet is a dynamic proximal mapping module, whose parameters can be dynamically adjusted at the inference stage and make it adapt to different imaging settings. Experimental results demonstrate that the proposed DPUNet can effectively handle multiple compressive imaging modalities under varying sampling ratios and noise levels via only one trained model, and outperform the state-of-the-art approaches.


title: SuperCaustics: Real-time, open-source simulation of transparent objects for deep learning applications

abstract: Transparent objects are a very challenging problem in computer vision. They are hard to segment or classify due to their lack of precise boundaries, and there is limited data available for training deep neural networks. As such, current solutions for this problem employ rigid synthetic datasets, which lack flexibility and lead to severe performance degradation when deployed on real-world scenarios. In particular, these synthetic datasets omit features such as refraction, dispersion and caustics due to limitations in the rendering pipeline. To address this issue, we present SuperCaustics, a real-time, open-source simulation of transparent objects designed for deep learning applications. SuperCaustics features extensive modules for stochastic environment creation; uses hardware ray-tracing to support caustics, dispersion, and refraction; and enables generating massive datasets with multi-modal, pixel-perfect ground truth annotations. To validate our proposed system, we trained a deep neural network from scratch to segment transparent objects in difficult lighting scenarios. Our neural network achieved performance comparable to the state-of-the-art on a real-world dataset using only 10% of the training data and in a fraction of the training time. Further experiments show that a model trained with SuperCaustics can segment different types of caustics, even in images with multiple overlapping transparent objects. To the best of our knowledge, this is the first such result for a model trained on synthetic data. Both our open-source code and experimental data are freely available online.


title: Reductions of certain 3-dimensional crystalline representations

abstract: We compute reduction $\bar \rho $ of 3-dimensional irreducible crystalline representations $\rho$ of $G_{\mathbb Q_p}$ with Hodge-Tate weights $\{0, r , s\}$ satisfying $2 \leq r \leq p-2, \ \ 2+p \leq s \leq r + p-2.$ If $\bar \rho$ is irreducible then we show that the generic fiber of crystalline universal deformation ring of $\bar \rho$ is connected.


title: 3D Brain Reconstruction by Hierarchical Shape-Perception Network from a Single Incomplete Image

abstract: 3D shape reconstruction is essential in the navigation of minimally-invasive and auto robot-guided surgeries whose operating environments are indirect and narrow, and there have been some works that focused on reconstructing the 3D shape of the surgical organ through limited 2D information available. However, the lack and incompleteness of such information caused by intraoperative emergencies (such as bleeding) and risk control conditions have not been considered. In this paper, a novel hierarchical shape-perception network (HSPN) is proposed to reconstruct the 3D point clouds (PCs) of specific brains from one single incomplete image with low latency. A branching predictor and several hierarchical attention pipelines are constructed to generate point clouds that accurately describe the incomplete images and then complete these point clouds with high quality. Meanwhile, attention gate blocks (AGBs) are designed to efficiently aggregate geometric local features of incomplete PCs transmitted by hierarchical attention pipelines and internal features of reconstructing point clouds. With the proposed HSPN, 3D shape perception and completion can be achieved spontaneously. Comprehensive results measured by Chamfer distance and PC-to-PC error demonstrate that the performance of the proposed HSPN outperforms other competitive methods in terms of qualitative displays, quantitative experiment, and classification evaluation.


title: VisDA-2021 Competition Universal Domain Adaptation to Improve Performance on Out-of-Distribution Data

abstract: Progress in machine learning is typically measured by training and testing a model on the same distribution of data, i.e., the same domain. This over-estimates future accuracy on out-of-distribution data. The Visual Domain Adaptation (VisDA) 2021 competition tests models' ability to adapt to novel test distributions and handle distributional shift. We set up unsupervised domain adaptation challenges for image classifiers and will evaluate adaptation to novel viewpoints, backgrounds, modalities and degradation in quality. Our challenge draws on large-scale publicly available datasets but constructs the evaluation across domains, rather that the traditional in-domain bench-marking. Furthermore, we focus on the difficult "universal" setting where, in addition to input distribution drift, methods may encounter missing and/or novel classes in the target dataset. Performance will be measured using a rigorous protocol, comparing to state-of-the-art domain adaptation methods with the help of established metrics. We believe that the competition will encourage further improvement in machine learning methods' ability to handle realistic data in many deployment scenarios.


title: Finding signatures of the nuclear symmetry energy in heavy-ion collisions with deep learning

abstract: A deep convolutional neural network (CNN) is developed to study symmetry energy $E_{\rm sym}(\rho)$ effects by learning the mapping between the symmetry energy and the two-dimensional (transverse momentum and rapidity) distributions of protons and neutrons in heavy-ion collisions. Supervised training is performed with labelled data-set from the ultrarelativistic quantum molecular dynamics (UrQMD) model simulation. It is found that, by using proton spectra on event-by-event basis as input, the accuracy for classifying the soft and stiff $E_{\rm sym}(\rho)$ is about 60% due to large event-by-event fluctuations, while by setting event-summed proton spectra as input, the classification accuracy increases to 98%. The accuracy for 5-label (5 different $E_{\rm sym}(\rho)$) classification task are about 58% and 72% by using proton and neutron spectra, respectively. For the regression task, the mean absolute error (MAE) which measures the average magnitude of the absolute differences between the predicted and actual $L$ (the slope parameter of $E_{\rm sym}(\rho)$) are about 20.4 and 14.8 MeV by using proton and neutron spectra, respectively. Fingerprints of the density-dependent nuclear symmetry energy on the transverse momentum and rapidity distributions of protons and neutrons can be identified by convolutional neural network algorithm.


title: Beamforming Design and Power Allocation for Transmissive RMS-based Transmitter Architectures

abstract: This letter investigates a downlink multiple input single output (MISO) system based on transmissive reconfigurable metasurface (RMS) transmitter. Specifically, a transmitter design based on a transmissive RMS equipped with a feed antenna is first proposed. Then, in order to maximize the achievable sum-rate of the system, the beamforming design and power allocation are jointly optimized. Since the optimization variables are coupled, this formulated optimization problem is non-convex, so it is difficult to solve it directly. To solve this problem, we propose an alternating optimization (AO) technique based on difference-of-convex (DC) programming and successive convex approximation (SCA). Simulation results verify that the proposed algorithm can achieve convergence and improve the achievable sum-rate of the system.


title: Post-Treatment Confounding in Causal Mediation Studies: A Cutting-Edge Problem and A Novel Solution via Sensitivity Analysis

abstract: In causal mediation studies that decompose an average treatment effect into a natural indirect effect (NIE) and a natural direct effect (NDE), examples of post-treatment confounding are abundant. Past research has generally considered it infeasible to adjust for a post-treatment confounder of the mediator-outcome relationship due to incomplete information: it is observed under the actual treatment condition while missing under the counterfactual treatment condition. This study proposes a new sensitivity analysis strategy for handling post-treatment confounding and incorporates it into weighting-based causal mediation analysis without making extra identification assumptions. Under the sequential ignorability of the treatment assignment and of the mediator, we obtain the conditional distribution of the post-treatment confounder under the counterfactual treatment as a function of not just pretreatment covariates but also its counterpart under the actual treatment. The sensitivity analysis then generates a bound for the NIE and that for the NDE over a plausible range of the conditional correlation between the post-treatment confounder under the actual and that under the counterfactual conditions. Implemented through either imputation or integration, the strategy is suitable for binary as well as continuous measures of post-treatment confounders. Simulation results demonstrate major strengths and potential limitations of this new solution. A re-analysis of the National Evaluation of Welfare-to-Work Strategies (NEWWS) Riverside data reveals that the initial analytic results are sensitive to omitted post-treatment confounding.


title: Trajectory Design for UAV-Based Internet-of-Things Data Collection: A Deep Reinforcement Learning Approach

abstract: In this paper, we investigate an unmanned aerial vehicle (UAV)-assisted Internet-of-Things (IoT) system in a sophisticated three-dimensional (3D) environment, where the UAV's trajectory is optimized to efficiently collect data from multiple IoT ground nodes. Unlike existing approaches focusing only on a simplified two-dimensional scenario and the availability of perfect channel state information (CSI), this paper considers a practical 3D urban environment with imperfect CSI, where the UAV's trajectory is designed to minimize data collection completion time subject to practical throughput and flight movement constraints. Specifically, inspired from the state-of-the-art deep reinforcement learning approaches, we leverage the twin-delayed deep deterministic policy gradient (TD3) to design the UAV's trajectory and present a TD3-based trajectory design for completion time minimization (TD3-TDCTM) algorithm. In particular, we set an additional information, i.e., the merged pheromone, to represent the state information of UAV and environment as a reference of reward which facilitates the algorithm design. By taking the service statuses of IoT nodes, the UAV's position, and the merged pheromone as input, the proposed algorithm can continuously and adaptively learn how to adjust the UAV's movement strategy. By interacting with the external environment in the corresponding Markov decision process, the proposed algorithm can achieve a near-optimal navigation strategy. Our simulation results show the superiority of the proposed TD3-TDCTM algorithm over three conventional non-learning based baseline methods.


title: Joint Communication and Trajectory Design for Intelligent Reflecting Surface Empowered UAV SWIPT Networks

abstract: Aiming at the limited battery capacity of widely deployed low-power smart devices in the Internet-of-things (IoT), this paper proposes a novel intelligent reflecting surface (IRS) empowered unmanned aerial vehicle (UAV) simultaneous wireless information and power transfer (SWIPT) network framework, in which IRS is used to reconstruct the wireless channel to enhance the wireless energy transmission efficiency and coverage area of the UAV SWIPT networks. In this paper, we formulate an achievable sum-rate maximization problem by jointly optimizing UAV trajectory, successive interference cancellation (SIC) decoding order, UAV transmit power allocation, power splitting (PS) ratio and IRS reflection coefficient while taking account of user non-orthogonal multiple access (NOMA) and a non-linear energy harvesting model. Due to the coupling of optimization variables, this problem is a complex non-convex optimization problem, and it is challenging to solve it directly. We first transform the problem, and then apply the alternating optimization (AO) algorithm framework to divide the transformed problem into four sub-problems to solve it. Specifically, by applying successive convex approximation (SCA), penalty function method and difference-convex (DC) programming, UAV trajectory, SIC decoding order, UAV transmit power allocation, PS ratio and IRS reflection coefficient are alternately optimized until the convergence is achieved. Numerical simulation results verify the effectiveness of our proposed algorithm compared to other algorithms.


title: Treating Divergent Perturbation Theory: Lessons from Exactly Solvable 2D Models at Large $N$

abstract: We consider the operator product expansion (OPE) of correlation functions in the supersymmetric $O(N)$ non-linear sigma model at sub-leading order in the large $N$ limit in order to study the cancellation between ambiguities coming from infrared renormalons and those coming from various operators in the OPE. As has been discussed in the context of supersymmetric Yang-Mills theory in four dimensions, supersymmetry presents a challenge to this cancellation. In a bid to solve this problem we consider $O(N)$ as a toy model. A background field method inspired by Polyakov's treatment of the renormalization of the bosonic $O(N)$ model is used to identify explicit operators in the OPE of the two-point functions of bosonic and fermionic fields in the model. In order to identify the coefficient functions in the OPE, the exact two-point functions at sub-leading order in large $N$ are expanded in powers of the natural infrared length scale. The ambiguities arising from renormalons in the coefficient functions and vacuum expectation values of operators in the OPE are shown to cancel to all orders. The question of supersymmetric Yang-Mills theory without matter remains open.


title: $L_p$ John ellipsoids for log-concave functions

abstract: The aim of this paper is to develop the $L_p$ John ellipsoid for the geometry of log-concave functions. Using the results of the $L_p$ Minkowski theory for log-concave function established in \cite{fan-xin-ye-geo2020}, we characterize the $L_p$ John ellipsoid for log-concave function, and establish some inequalities of the $L_p$ John ellipsoid for log-concave function. Finally, the analog of Ball's volume ratio inequality for the $L_p$ John ellipsoid of log-concave function is established.


title: Generating Large-scale Dynamic Optimization Problem Instances Using the Generalized Moving Peaks Benchmark

abstract: This document describes the generalized moving peaks benchmark (GMPB) and how it can be used to generate problem instances for continuous large-scale dynamic optimization problems. It presents a set of 15 benchmark problems, the relevant source code, and a performance indicator, designed for comparative studies and competitions in large-scale dynamic optimization. Although its primary purpose is to provide a coherent basis for running competitions, its generality allows the interested reader to use this document as a guide to design customized problem instances to investigate issues beyond the scope of the presented benchmark suite. To this end, we explain the modular structure of the GMPB and how its constituents can be assembled to form problem instances with a variety of controllable characteristics ranging from unimodal to highly multimodal, symmetric to highly asymmetric, smooth to highly irregular, and various degrees of variable interaction and ill-conditioning.