title: BN-HTRd: A Benchmark Dataset for Document Level Offline Bangla Handwritten Text Recognition (HTR) and Line Segmentation

abstract: We introduce a new dataset for offline Handwritten Text Recognition (HTR) from images of Bangla scripts comprising words, lines, and document-level annotations. The BN-HTRd dataset is based on the BBC Bangla News corpus, meant to act as ground truth texts. These texts were subsequently used to generate the annotations that were filled out by people with their handwriting. Our dataset includes 788 images of handwritten pages produced by approximately 150 different writers. It can be adopted as a basis for various handwriting classification tasks such as end-to-end document recognition, word-spotting, word or line segmentation, and so on. We also propose a scheme to segment Bangla handwritten document images into corresponding lines in an unsupervised manner. Our line segmentation approach takes care of the variability involved in different writing styles, accurately segmenting complex handwritten text lines of curvilinear nature. Along with a bunch of pre-processing and morphological operations, both Hough line and circle transforms were employed to distinguish different linear components. In order to arrange those components into their corresponding lines, we followed an unsupervised clustering approach. The average success rate of our segmentation technique is 81.57% in terms of FM metrics (similar to F-measure) with a mean Average Precision (mAP) of 0.547.


title: Towards a Deep Multi-layered Dialectal Language Analysis: A Case Study of African-American English

abstract: Currently, natural language processing (NLP) models proliferate language discrimination leading to potentially harmful societal impacts as a result of biased outcomes. For example, part-of-speech taggers trained on Mainstream American English (MAE) produce non-interpretable results when applied to African American English (AAE) as a result of language features not seen during training. In this work, we incorporate a human-in-the-loop paradigm to gain a better understanding of AAE speakers' behavior and their language use, and highlight the need for dialectal language inclusivity so that native AAE speakers can extensively interact with NLP systems while reducing feelings of disenfranchisement.


title: Celestial Yang-Mills Amplitudes and D=4 Conformal Blocks

abstract: We discuss the properties of recently constructed "single-valued" celestial four-gluon amplitudes. We show that the amplitude factorizes into the "current" part and the "scalar" part. The current factor is given by the group-dependent part of the Wess-Zumino-Witten correlator of four holomorphic currents with a non-vanishing level of Ka\v{c}-Moody algebra. The scalar factor can be expressed in terms of a complex integral of the Koba-Nielsen form, similar to the integrals describing four-point correlators in Coulomb gas models and, more generally, in the infinite central charge limit of Liouville theory. The scalar part can be also obtained by a dimensional reduction of a single D=4 conformal block and the shadow block from Minkowski space to the celestial sphere.


title: Explainable Global Error Weighted on Feature Importance: The xGEWFI metric to evaluate the error of data imputation and data augmentation

abstract: Evaluating the performance of an algorithm is crucial. Evaluating the performance of data imputation and data augmentation can be similar since both generated data can be compared with an original distribution. Although, the typical evaluation metrics have the same flaw: They calculate the feature's error and the global error on the generated data without weighting the error with the feature importance. The result can be good if all of the feature's importance is similar. However, in most cases, the importance of the features is imbalanced, and it can induce an important bias on the features and global errors. This paper proposes a novel metric named "Explainable Global Error Weighted on Feature Importance"(xGEWFI). This new metric is tested in a whole preprocessing method that 1. detects the outliers and replaces them with a null value. 2. imputes the data missing, and 3. augments the data. At the end of the process, the xGEWFI error is calculated. The distribution error between the original and generated data is calculated using a Kolmogorov-Smirnov test (KS test) for each feature. Those results are multiplied by the importance of the respective features, calculated using a Random Forest (RF) algorithm. The metric result is expressed in an explainable format, aiming for an ethical AI.


title: Repair kinetics of DSB-foci induced by proton and helium ion microbeams of different energies

abstract: In this work, the induction and repair of radiation-induced 53BP1 foci were studied in human umbilical vein endothelial cells irradiated at the PTB microbeam with protons and {\alpha}-particles of different energies. The data were analyzed in terms of the mean number of 53BP1 foci induced by the different ion beams. The number of 53BP1 foci found at different times post-irradiation suggests that the disappearance of foci follows first order kinetics. The mean number of initially produced foci shows the expected increase with LET. The most interesting finding of this work is that the absolute number of persistent foci increases with LET but not their fraction. Furthermore, protons seem to produce more persistent foci as compared to {\alpha}-particles of even higher LET. This may be seen as experimental evidence that protons may be more effective in producing severe DNA lesions, as was already shown in other work, and that LET may not be the best suited parameter to characterize radiation quality


title: ck-means, a novel unsupervised learning method that combines fuzzy and crispy clustering methods to extract intersecting data

abstract: Clustering data is a popular feature in the field of unsupervised machine learning. Most algorithms aim to find the best method to extract consistent clusters of data, but very few of them intend to cluster data that share the same intersections between two features or more. This paper proposes a method to do so. The main idea of this novel method is to generate fuzzy clusters of data using a Fuzzy C-Means (FCM) algorithm. The second part involves applying a filter that selects a range of minimum and maximum membership values, emphasizing the border data. A {\mu} parameter defines the amplitude of this range. It finally applies a k-means algorithm using the membership values generated by the FCM. Naturally, the data having similar membership values will regroup in a new crispy cluster. The algorithm is also able to find the optimal number of clusters for the FCM and the k-means algorithm, according to the consistency of the clusters given by the Silhouette Index (SI). The result is a list of data and clusters that regroup data sharing the same intersection, intersecting two features or more. ck-means allows extracting the very similar data that does not naturally fall in the same cluster but at the intersection of two clusters or more. The algorithm also always finds itself the optimal number of clusters.


title: BURSTT: Bustling Universe Radio Survey Telescope in Taiwan

abstract: Fast Radio Bursts (FRBs) are bright millisecond-duration radio transients that appear about 1,000 times per day, all-sky, for a fluence threshold 5 Jy ms at 600 MHz. The FRB radio-emission physics and the compact objects involved in these events are subjects of intense active debate. To better constrain source models, the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is optimized to discover and localize a large sample of rare, high-fluence, nearby FRBs. This is the population most amenable to multi-messenger, multi-wavelength follow-up, allowing deeper understanding of source mechanisms. BURSTT will provide horizon-to-horizon sky coverage with a half power field-of-view (FoV) of $\sim$10$^{4}$ deg$^{2}$, a 400 MHz effective bandwidth between 300-800 MHz, and sub-arcsecond localization, made possible using outrigger stations hundreds to thousands of km from the main array. Initially, BURSTT will employ 256 antennas. After tests of various antenna designs and optimization of system performance we plan to expand to 2048 antennas. We estimate that BURSTT-256 will detect and localize $\sim$100 bright ($\geq$100 Jy ms) FRBs per year. Another advantage of BURSTT's large FoV and continuous operation will be greatly enhanced monitoring of FRBs for repetition. The current lack of sensitive all-sky observations likely means that many repeating FRBs are currently cataloged as single-event FRBs.


title: Multi-scale Super-resolution Magnetic Resonance Spectroscopic Imaging with Adjustable Sharpness

abstract: Magnetic Resonance Spectroscopic Imaging (MRSI) is a valuable tool for studying metabolic activities in the human body, but the current applications are limited to low spatial resolutions. The existing deep learning-based MRSI super-resolution methods require training a separate network for each upscaling factor, which is time-consuming and memory inefficient. We tackle this multi-scale super-resolution problem using a Filter Scaling strategy that modulates the convolution filters based on the upscaling factor, such that a single network can be used for various upscaling factors. Observing that each metabolite has distinct spatial characteristics, we also modulate the network based on the specific metabolite. Furthermore, our network is conditioned on the weight of adversarial loss so that the perceptual sharpness of the super-resolved metabolic maps can be adjusted within a single network. We incorporate these network conditionings using a novel Multi-Conditional Module. The experiments were carried out on a 1H-MRSI dataset from 15 high-grade glioma patients. Results indicate that the proposed network achieves the best performance among several multi-scale super-resolution methods and can provide super-resolved metabolic maps with adjustable sharpness.


title: TransResU-Net: Transformer based ResU-Net for Real-Time Colonoscopy Polyp Segmentation

abstract: Colorectal cancer (CRC) is one of the most common causes of cancer and cancer-related mortality worldwide. Performing colon cancer screening in a timely fashion is the key to early detection. Colonoscopy is the primary modality used to diagnose colon cancer. However, the miss rate of polyps, adenomas and advanced adenomas remains significantly high. Early detection of polyps at the precancerous stage can help reduce the mortality rate and the economic burden associated with colorectal cancer. Deep learning-based computer-aided diagnosis (CADx) system may help gastroenterologists to identify polyps that may otherwise be missed, thereby improving the polyp detection rate. Additionally, CADx system could prove to be a cost-effective system that improves long-term colorectal cancer prevention. In this study, we proposed a deep learning-based architecture for automatic polyp segmentation, called Transformer ResU-Net (TransResU-Net). Our proposed architecture is built upon residual blocks with ResNet-50 as the backbone and takes the advantage of transformer self-attention mechanism as well as dilated convolution(s). Our experimental results on two publicly available polyp segmentation benchmark datasets showed that TransResU-Net obtained a highly promising dice score and a real-time speed. With high efficacy in our performance metrics, we concluded that TransResU-Net could be a strong benchmark for building a real-time polyp detection system for the early diagnosis, treatment, and prevention of colorectal cancer. The source code of the proposed TransResU-Net is publicly available at https://github.com/nikhilroxtomar/TransResUNet.


title: Laboratory Demonstration of Real-Time Focal Plane Wavefront Control of Residual Atmospheric Speckles

abstract: Current and future high contrast imaging instruments aim to detect exoplanets at closer orbital separations, lower masses, and/or older ages than their predecessors. However, continually evolving speckles in the coronagraphic science image limit contrasts of state-of-the-art ground-based exoplanet imaging instruments. For ground-based adaptive optics (AO) instruments it remains challenging for most speckle suppression techniques to attenuate both the dynamic atmospheric as well as quasi-static instrumental speckles on-sky. We have proposed a focal plane wavefront sensing and control algorithm to address this challenge, called the Fast Atmospheric Self-coherent camera (SCC) Technique (FAST), which in theory enables the SCC to operate down to millisecond timescales even when only a few photons are detected per speckle. Here we present the first experimental results of FAST on the Santa Cruz Extreme AO Laboratory (SEAL) testbed. In particular, we illustrate the benefit of ``second stage'' AO-based focal plane wavefront control, demonstrating up to 5x contrast improvement with FAST closed-loop compensation of evolving residual atmospheric turbulence -- both for low and high order spatial modes -- down to 20 millisecond-timescales.


title: Norm Inequalities for Integral Operators on Cones

abstract: In this dissertation we explore the $[L^{\mathrm{p}},\ L^{q}]$-boundedness of certain integral operators on weighted spaces on cones in ${\mathbb R}^{n}.$ These integral operators are of the type $\displaystyle \int_{V}k(x,\ y)f(y)dy$ defined on a homogeneous cone $V$. The results of this dissertation are then applied to an important class of operators such as Riemann-Liouville's fractional integral operators, Weyl's fractional integral operators and Laplace's operators. As special cases of the above, we obtain an ${\mathbb R}^{n}$ -generalization of the celebrated Hardy's inequality on domains of positivity. We also prove dual results.


title: The stellar metallicity gradients of Local Group dwarf galaxies

abstract: Through a homogeneous analysis of spectroscopic literature data of red giant stars, we determine the radial metallicity profiles of 30 dwarf galaxies in the Local Group. We explore correlations between the calculated metallicity gradients and stellar mass, star formation history and environment, delivering the largest compilation to date of this type. The dwarf galaxies in our sample mostly show metallicity profiles decreasing with radius, with some exhibiting rather steep profiles. The derived metallicity gradients as a function of the half-light radius, $\nabla_{\rm [Fe/H]} (R/R_e)$, show no statistical differences when compared with the galaxies' morphological type, nor with their distance from the Milky Way or M31. No correlations are found with either stellar mass or star formation timescales. In particular, we do not find the linear relationship between $\nabla_{\rm [Fe/H]} (R/R_e)$ and the galaxies' median age $t_{50}$, as instead shown in the literature for a set of simulated systems. The presence of high angular momentum in some of our galaxies does not seem to have an impact on the gradient values. The strongest gradients in our sample are observed in systems that are likely to have experienced a past merger event. By excluding them, the analysed dwarf galaxies show mild gradients ($\sim -0.1$ dex $R_e^{-1}$) with little scatter between them, regardless of their stellar mass, dynamical state, and star formation history. These results are in good agreement with different sets of simulations presented in the literature and analysed using the same method as for the observed sample. The interplay between the multitude of factors that could drive the formation of metallicity gradients in dwarf galaxies likely combine in complex ways to produce in general comparable values.


title: A Southern Photometric Quasar Catalog from the Dark Energy Survey Data Release 2

abstract: We present a catalog of 1.4 million photometrically-selected quasar candidates in the southern hemisphere over the $\sim 5000\,{\rm deg^2}$ Dark Energy Survey (DES) wide survey area. We combine optical photometry from the DES second data release (DR2) with available near-infrared (NIR) and the all-sky unWISE mid-infrared photometry in the selection. We build models of quasars, galaxies, and stars with multivariate Skew-t distributions in the multi-dimensional space of relative fluxes as functions of redshift (or color for stars) and magnitude. Our selection algorithm assigns probabilities for quasars, galaxies, and stars, and simultaneously calculates photometric redshifts (photo-$z$) for quasar and galaxy candidates. Benchmarking on spectroscopically confirmed objects, we successfully classify (with photometry) 94.7% of quasars, 99.3% of galaxies, and 96.3% of stars when all IR bands (NIR $YJHK$ and WISE $W1W2$) are available. The classification and photo-$z$ regression success rates decrease when fewer bands are available. Our quasar (galaxy) photo-$z$ quality, defined as the fraction of objects with the difference between the photo-$z$ $z_p$ and the spectroscopic redshift $z_s$, $|\Delta z| = |z_s - z_p|/(1 + z_s)\le 0.1$, is 92.2% (98.1%) when all IR bands are available, decreasing to 72.2% (90.0%) using optical DES data only. Our photometric quasar catalog achieves estimated completeness of 89% and purity of 79% at $r<21.5$ (0.68 million quasar candidates), with reduced completeness and purity at $21.5<r\lesssim 24$. Among the 1.4 million quasar candidates, 87,857 have existing spectra and 84,978 (96.7%) of them are spectroscopically confirmed quasars. Finally, we provide quasar, galaxy, and star probabilities for all (0.69 billion) photometric sources in the DES DR2 coadded photometric catalog.


title: Shadows Shed Light on 3D Objects

abstract: 3D reconstruction is a fundamental problem in computer vision, and the task is especially challenging when the object to reconstruct is partially or fully occluded. We introduce a method that uses the shadows cast by an unobserved object in order to infer the possible 3D volumes behind the occlusion. We create a differentiable image formation model that allows us to jointly infer the 3D shape of an object, its pose, and the position of a light source. Since the approach is end-to-end differentiable, we are able to integrate learned priors of object geometry in order to generate realistic 3D shapes of different object categories. Experiments and visualizations show that the method is able to generate multiple possible solutions that are consistent with the observation of the shadow. Our approach works even when the position of the light source and object pose are both unknown. Our approach is also robust to real-world images where ground-truth shadow mask is unknown.


title: Efficient Order-Optimal Preconditioners for Implicit Runge-Kutta and Runge-Kutta-Nystr\"om Methods Applicable to a Large Class of Parabolic and Hyperbolic PDEs

abstract: We generalize previous work by Mardal, Nilssen, and Staff (2007, SIAM J. Sci. Comp. v. 29, pp. 361-375) and Rana, Howle, Long, Meek, and Milestone (2021, SIAM J. Sci. Comp. v. 43, p. 475-495) on order-optimal preconditioners for parabolic PDEs to a larger class of differential equations and methods. The problems considered are those of the forms $u_{t}=-\mathcal{K}u+g$ and $u_{tt}=-\mathcal{{K}}u+g$, where the operator $\mathcal{{K}}$ is defined by $\mathcal{{K}}u:=-\nabla\cdot\left(\alpha\nabla u\right)+\beta u$ and the functions $\alpha$ and $\beta$ are restricted so that $\alpha>0$, and $\beta\ge0$. The methods considered are A-stable implicit Runge--Kutta methods for the parabolic equation and implicit Runge--Kutta--Nystr\"om methods for the hyperbolic equation. We prove the order optimality of a class of block preconditioners for the stage equation system arising from these problems, and furthermore we show that the LD and DU preconditioners of Rana et al. are in this class. We carry out numerical experiments on several test problems in this class -- the 2D diffusion equation, Pennes bioheat equation, the wave equation, and the Klein--Gordon equation, with both constant and variable coefficients. Our experiments show that these preconditioners, particularly the LD preconditioner, are successful at reducing the condition number of the systems as well as improving the convergence rate and solve time for GMRES applied to the stage equations.


title: Discovery of CH2CCHC4H and a rigorous detection of CH2CCHC3N in TMC-1 with the QUIJOTE line survey

abstract: Using the QUIJOTE line survey in the 32.0-50.4 GHz range, we report the discovery of the molecule CH2CCHC4H towards the prestellar cold core TMC-1 in the Taurus region. We also present a rigorous detection of CH2CCHC3N, along with its detailed analysis. We identified a total of twenty rotational transitions for each one of these molecules. The rotational quantum numbers range from Ju=17 up to 24 and Ka<=3. The column density for CH2CCHC4H is N=(2.2+/-0.2)x 1E12 cm-2, while for CH2CCHC3N, we derived N=(1.2+/-0.15) x 1E11 cm-2. The rotational temperature is 9.0+/-0.5 K for both species. The abundance ratio between CH2CCHC4H and CH2CCHC3N is 18+/-4. We also compared the column densities of these species with those of their isomers CH3C6H and CH3C5N, derived from their J=20-19 up to J=30-29 rotational transitions observed with the QUIJOTE line survey. The observed abundances for all these species are reasonably well explained by state-of-the-art chemical models of TMC-1. The observed astronomical frequencies were merged with laboratory frequencies from the literature to derive improved spectroscopic parameters.


title: A bijection between $K$-Kohnert diagrams and reverse set-valued tableaux

abstract: Lascoux polynomials are $K$-theoretic analogues of the key polynomials. They both have combinatorial formulas involving tableaux: reverse set-valued tableaux ($\mathsf{RSVT}$) rule for Lascoux polynomials and reverse semistandard Young tableaux ($\mathsf{RSSYT}$) rule for key polynomials. Furthermore, key polynomials have a simple algorithmic model in terms of Kohnert diagrams, which are in bijection with $\mathsf{RSSYT}$. Ross and Yong introduced $K$-Kohnert diagrams, which are analogues of Kohnert diagrams. They conjectured a $K$-Kohnert diagram rule for Lascoux polynomials. We establish this conjecture by constructing a weight-preserving bijection between $\mathsf{RSVT}$ and $K$-Kohnert diagrams.


title: Robust Group Synchronization via Quadratic Programming

abstract: We propose a novel quadratic programming formulation for estimating the corruption levels in group synchronization, and use these estimates to solve this problem. Our objective function exploits the cycle consistency of the group and we thus refer to our method as detection and estimation of structural consistency (DESC). This general framework can be extended to other algebraic and geometric structures. Our formulation has the following advantages: it can tolerate corruption as high as the information-theoretic bound, it does not require a good initialization for the estimates of group elements, it has a simple interpretation, and under some mild conditions the global minimum of our objective function exactly recovers the corruption levels. We demonstrate the competitive accuracy of our approach on both synthetic and real data experiments of rotation averaging.


title: Space-time POD and the Hankel matrix

abstract: Time-delay embedding is an increasingly popular starting point for data-driven reduced-order modeling efforts. In particular, the singular value decomposition (SVD) of a block Hankel matrix formed from successive delay embeddings of the state of a dynamical system lies at the heart of several popular reduced-order modeling methods. In this paper, we show that the left singular vectors of this Hankel matrix are a discrete approximation of classical space-time proper orthogonal decomposition (POD) modes, and the singular values are square roots of the POD energies. This connection establishes a clear interpretation of the Hankel modes grounded in classical theory, and we gain insights into the Hankel modes by instead analyzing the equivalent discrete space-time POD modes in terms of the correlation matrix formed by multiplying the Hankel matrix by its conjugate transpose. These insights include the distinct meaning of rows and columns, the implied norm in which the modes are optimal, the impact of the time step between snapshots on the modes, and an interpretation of the embedding dimension/height of the Hankel matrix in terms of the time window on which the modes are optimal. Moreover, the connections we establish offer opportunities to improve the convergence and computation time in certain practical cases, and to improve the accuracy of the modes with the same data. Finally, popular variants of POD, namely the standard space-only POD and spectral POD, are recovered in the limits that snapshots used to form each column of the Hankel matrix represent flow evolution over short and long times, respectively.


title: Towards Consensus: Reducing Polarization by Perturbing Social Networks

abstract: This paper studies how a centralized planner can modify the structure of a social or information network to reduce polarization. First, polarization is found to be highly dependent on degree and structural properties of the network -- including the well-known isoperimetric number (i.e., Cheeger constant). We then formulate the planner's problem under full information, and motivate disagreement-seeking and coordinate descent heuristics. A novel setting for the planner in which the population's innate opinions are adversarially chosen is introduced, and shown to be equivalent to maximization of the Laplacian's spectral gap. We prove bounds for the effectiveness of a strategy that adds edges between vertices on opposite sides of the cut induced by the spectral gap's eigenvector. Finally, these strategies are evaluated on six real-world and synthetic networks. In several networks, we find that polarization can be significantly reduced through the addition of a small number of edges.


title: Excellent thermoelectric performance and impressive optoelectronic properties of Janus monolayer of ZrXY(X=O, S) (Y=S, Se)

abstract: Lower-dimensional TMDC materials are suitable for thermoelectric applications for their specific quantum confinement and being distinct in the density of states (DOS). Here we investigated thermoelectric parameters of the 2D TMDC monolayer of ZrXY ((X=O, S,) (Y=S, Se)) by using Density Functional Theory (DFT) combined with Boltzmann Transport Equation (BTE)


title: A review of machine learning concepts and methods for addressing challenges in probabilistic hydrological post-processing and forecasting

abstract: Probabilistic forecasting is receiving growing attention nowadays in a variety of applied fields, including hydrology. Several machine learning concepts and methods are notably relevant towards addressing the major challenges of formalizing and optimizing probabilistic forecasting implementations, as well as the equally important challenge of identifying the most useful ones among these implementations. Nonetheless, practically-oriented reviews focusing on such concepts and methods, and on how these can be effectively exploited in the above-outlined essential endeavour, are currently missing from the probabilistic hydrological forecasting literature. This absence holds despite the pronounced intensification in the research efforts for benefitting from machine learning in this same literature. It also holds despite the substantial relevant progress that has recently emerged, especially in the field of probabilistic hydrological post-processing, which traditionally provides the hydrologists with probabilistic hydrological forecasting implementations. Herein, we aim to fill this specific gap. In our review, we emphasize key ideas and information that can lead to effective popularizations, as such an emphasis can support successful future implementations and further scientific developments. In the same forward-looking direction, we identify open research questions and propose ideas to be explored in the future.


title: On the Importance of Well-Defined Thermal Correlation Functions in Simulating Vibronic Spectra

abstract: Two difficulties associated with the computations of thermal vibrational correlation functions are discussed. The first one is the lack of a well-behaved expression that is valid at both high-temperature and $T \to 0$ K limits. Specifically, if the partition function and the propagator are considered separately, then thermal vibrational correlation functions may have an indeterminate form 0/0 in the limit $T \to 0$ K. This difficulty is resolved when the partition function and the propagator are jointly considered in the harmonic approximation, which allows a problematic term that emanates from the zero-point energy to be cancelled out thereby producing a thermal correlation function with a determinate form in $T \to 0$ K limit. The second difficulty is related to the multivaluedness of the vibrational correlation function. We show numerically that an improper selection of branch leads to discontinuities in the computed correlation function and an incorrect vibronic spectra. We propose a phase tracking procedure that ensures continuity of both real and imaginary parts of the correlation function to recover the correct spectra. We support our findings by simulating the UV-vis absorption spectra of pentacene at 4 K and benzene at 298 K. Both are found to be in good agreement with their experimental counterparts.


title: Revisiting $f(R)$ gravity's rainbow: Inflation and primordial fluctuations

abstract: We study inflation and the generation of primordial fluctuations in $f(R)$ gravity's rainbow. We calculate the cosmological perturbations and then the scalar and tensor primordial power spectrum. We contrast the predictions of the model with the current observational data from PLANCK and BICEP/Keck. Particularly, we found new results for the scalar spectral index $n_s$ and the tensor-to-scalar ratio $r$ along with new observational constraints on the rainbow functions.


title: Asymptotic $C^{1,\gamma}$-regularity for value functions to uniformly elliptic dynamic programming principles

abstract: In this paper we prove an asymptotic $C^{1,\gamma}$-estimate for value functions of stochastic processes related to uniformly elliptic dynamic programming principles. As an application, this allows us to pass to the limit with a discrete gradient and then to obtain a $C^{1,\gamma}$-result for the corresponding limit PDE.


title: Cluster Generation via Deep Energy-Based Model

abstract: We present a new approach for the generation of stable structures of nanoclusters using deep learning methods. Our method consists in constructing an artificial potential energy surface, with local minima corresponding to the most stable structures and which is much smoother than "real" potential in the intermediate regions of the configuration space. To build the surface, graph convolutional networks are used. The method can extrapolates the potential surface to cases of structures with larger number of atoms than was used in training. Thus, having a sufficient number of low-energy structures in the training set, the method allows to generate new candidates for the ground-state structures, including ones with larger number of atoms. We applied the approach to silica clusters $(SiO_2)_n$ and for the first time found the stable structures with n=28...51. The method is universal and does not depend on the atomic composition and number of atoms.


title: Mining S-PLUS for Metal-Poor Stars in the Milky Way

abstract: This work presents the medium-resolution ($R \sim 1,500$) spectroscopic follow-up of 522 low-metallicity star candidates from the Southern Photometric Local Universe Survey (S-PLUS). The objects were selected from narrow-band photometry, taking advantage of the metallicity-sensitive S-PLUS colors. The follow-up observations were conducted with the Blanco and Gemini South telescopes, using the COSMOS and GMOS spectrographs, respectively. The stellar atmospheric parameters (T$_{\rm eff}$, log$g$, and [Fe/H]), as well as carbon and $\alpha$-element abundances, were calculated for the program stars in order to assess the efficacy of the color selection. Results show that $92^{+2}_{-3}\%$ of the observed stars have [Fe/H]$\leq -1.0$, $83^{+3}_{-3}\%$ have [Fe/H]$\leq -2.0$, and $15^{+3}_{-3}\%$ have [Fe/H]$\leq -3.0$, including two ultra metal-poor stars ([Fe/H]$\leq -4.0$). The 80th percentile for the metallicity cumulative distribution function of the observed sample is [Fe/H]$= -2.04$. The sample also includes 68 Carbon-Enhanced Metal-Poor (CEMP) stars. Based on the calculated metallicities, further S-PLUS, color cuts are proposed, which can increase the fractions of stars with [Fe/H]$\leq -1.0$ and $\leq -2.0$ to $98\%$ and $88\%$, respectively. Such high success rates enable targeted high-resolution spectroscopic follow-up efforts, as well as provide selection criteria for fiber-fed multiplex spectroscopic surveys.


title: Towards Efficient Active Learning of PDFA

abstract: We propose a new active learning algorithm for PDFA based on three main aspects: a congruence over states which takes into account next-symbol probability distributions, a quantization that copes with differences in distributions, and an efficient tree-based data structure. Experiments showed significant performance gains with respect to reference implementations.


title: Mapping renormalized coupled cluster methods to quantum computers through a compact unitary representation of non-unitary operators

abstract: Non-unitary theories are commonly seen in the classical simulations of quantum systems. Among these theories, the method of moments of coupled-cluster equations (MMCCs) and the ensuing classes of the renormalized coupled-cluster (CC) approaches have evolved into one of the most accurate approaches to describe correlation effects in various quantum systems. The MMCC formalism provides an effective way for correcting energies of approximate CC formulations (parent theories) using moments, or CC equations, that are not used to determine approximate cluster amplitudes. In this paper, we propose a quantum algorithm for computing MMCC ground-state energies that provide two main advantages over classical computing or other quantum algorithms: (i) the possibility of forming superpositions of CC moments of arbitrary ranks in the entire Hilbert space and using an arbitrary form of the parent cluster operator for MMCC expansion; and (ii) significant reduction in the number of measurements in quantum simulation through a compact unitary representation for a generally non-unitary operator. We illustrate the robustness of our approach over a broad class of test cases, including ~40 molecular systems with varying basis sets encoded using 4~40 qubits, and exhibit the detailed MMCC analysis for the 8-qubit half-filled, four-site, single impurity Anderson model and 12-qubit hydrogen fluoride molecular system from the corresponding noise-free and noisy MMCC quantum computations. We also outline the extension of MMCC formalism to the case of unitary CC wave function ansatz.


title: Sterically Driven Current Reversal in a Model Molecular Motor

abstract: Simulations can help unravel the complicated ways in which molecular structure determines function. Here, we use molecular simulations to show how slight alterations of a molecular motor's structure can cause the motor's typical dynamical behavior to reverse directions. Inspired by autonomous synthetic catenane motors, we study the molecular dynamics of a minimal motor model, consisting of a shuttling ring that moves along a track containing interspersed binding sites and catalytic sites. The binding sites attract the shuttling ring while the catalytic sites speed up a reaction between molecular species, which can be thought of as fuel and waste. When that fuel and waste are held in a nonequilibrium steady-state concentration, the free energy from the reaction drives directed motion of the shuttling ring along the track. Using this model and nonequilibrium molecular dynamics, we show that the shuttling ring's direction can be reversed by simply adjusting the spacing between binding and catalytic sites on the track. We present a steric mechanism behind the current reversal, supported by kinetic measurements from the simulations. These results demonstrate how molecular simulation can guide future development of artificial molecular motors.


title: Radiation pressure acceleration of high-quality ion beams using ultrashort laser pulses

abstract: The generation of compact, high-energy ion beams is one of the most promising applications of intense laser-matter interactions, but the control of the beam spectral quality remains an outstanding challenge. We show that in radiation pressure acceleration of a thin solid target the onset of electron heating is determined by the growth of the Rayleigh-Taylor-like instability at the front surface and must be controlled to produce ion beams with high spectral quality in the light sail regime. The growth rate of the instability imposes an upper limit on the laser pulse duration and intensity to achieve high spectral beam quality and we demonstrate that under this optimal regime, the maximum peak ion beam energy per nucleon is independent of target density, composition, and laser energy (transverse spot size). Our predictions are validated by two- and three-dimensional particle-in-cell simulations, which indicate that for recent and upcoming experimental facilities using ultrashort ($\lesssim 25$ fs) laser pulses it is possible to produce $100 - 300$ MeV proton beams with $\sim 30\%$ energy spread and high laser-to-proton energy conversion efficiency.


title: Orthogonal Rational Approximation of Transfer Functions for High-Frequency Circuits

abstract: Rational function approximations find applications in many areas including macro-modeling of high-frequency circuits, model order reduction for controller design, interpolation and extrapolation of system responses, surrogate models for high-energy physics, and approximation of elementary mathematical functions. The unknown denominator polynomial of the model results in a non-linear problem, which can be replaced with successive solutions of linearized problems following the Sanathanan-Koerner (SK) iteration. An orthogonal basis can be obtained based on Arnoldi resulting in a stabilized SK iteration. We present an extension of the stabilized SK, called Orthogonal Rational Approximation (ORA), which ensures real polynomial coefficients and stable poles for realizability of electrical networks. We also introduce an efficient implementation of ORA for multi-port networks based on a block QR decomposition.


title: Intelligent Blockchain-based Edge Computing via Deep Reinforcement Learning: Solutions and Challenges

abstract: The convergence of mobile edge computing (MEC) and blockchain is transforming the current computing services in wireless Internet-of-Things networks, by enabling task offloading with security enhancement based on blockchain mining. Yet the existing approaches for these enabling technologies are isolated, providing only tailored solutions for specific services and scenarios. To fill this gap, we propose a novel cooperative task offloading and blockchain mining (TOBM) scheme for a blockchain-based MEC system, where each edge device not only handles computation tasks but also deals with block mining for improving system utility. To address the latency issues caused by the blockchain operation in MEC, we develop a new Proof-of-Reputation consensus mechanism based on a lightweight block verification strategy. To accommodate the highly dynamic environment and high-dimensional system state space, we apply a novel distributed deep reinforcement learning-based approach by using a multi-agent deep deterministic policy gradient algorithm. Experimental results demonstrate the superior performance of the proposed TOBM scheme in terms of enhanced system reward, improved offloading utility with lower blockchain mining latency, and better system utility, compared to the existing cooperative and non-cooperative schemes. The paper concludes with key technical challenges and possible directions for future blockchain-based MEC research.


title: LIMO: Latent Inceptionism for Targeted Molecule Generation

abstract: Generation of drug-like molecules with high binding affinity to target proteins remains a difficult and resource-intensive task in drug discovery. Existing approaches primarily employ reinforcement learning, Markov sampling, or deep generative models guided by Gaussian processes, which can be prohibitively slow when generating molecules with high binding affinity calculated by computationally-expensive physics-based methods. We present Latent Inceptionism on Molecules (LIMO), which significantly accelerates molecule generation with an inceptionism-like technique. LIMO employs a variational autoencoder-generated latent space and property prediction by two neural networks in sequence to enable faster gradient-based reverse-optimization of molecular properties. Comprehensive experiments show that LIMO performs competitively on benchmark tasks and markedly outperforms state-of-the-art techniques on the novel task of generating drug-like compounds with high binding affinity, reaching nanomolar range against two protein targets. We corroborate these docking-based results with more accurate molecular dynamics-based calculations of absolute binding free energy and show that one of our generated drug-like compounds has a predicted $K_D$ (a measure of binding affinity) of $6 \cdot 10^{-14}$ M against the human estrogen receptor, well beyond the affinities of typical early-stage drug candidates and most FDA-approved drugs to their respective targets. Code is available at https://github.com/Rose-STL-Lab/LIMO.


title: Evolutionary Random Graph for Bitcoin Overlay and Blockchain Mining Networks

abstract: The world economy is experiencing the novel adoption of distributed currencies that are free from the control of central banks. Distributed currencies suffer from extreme volatility, and this can lead to catastrophic implications during future economic crisis. Understanding the dynamics of this new type of currencies is vital for empowering supervisory bodies from current reactive and manual incident responders to more proactive and well-informed planners. Bitcoin, the first and dominant distributed cryptocurrency, is still notoriously vague, especially for a financial instrument with market value exceeding 1 trillion. Modeling of bitcoin overlay network poses a number of important theoretical and methodological challenges. Current measuring approaches, for example, fail to identify the real network size of bitcoin miners. This drastically undermines the ability to predict forks, the suitable mining difficulty and most importantly the resilience of the network supporting bitcoin. In this work, we developed Evolutionary Random Graph, a theoretical model that describes the network of bitcoin miners. The correctness of this model has been validated using simulated and measure real bitcoin data. We then predicted forking, optimal mining difficulty, network size and consequently the network's inability to stand a drastic drop in bitcoin price using the current mining configuration.


title: Diffusion models as plug-and-play priors

abstract: We consider the problem of inferring high-dimensional data $\mathbf{x}$ in a model that consists of a prior $p(\mathbf{x})$ and an auxiliary differentiable constraint $c(\mathbf{x},\mathbf{y})$ on $x$ given some additional information $\mathbf{y}$. In this paper, the prior is an independently trained denoising diffusion generative model. The auxiliary constraint is expected to have a differentiable form, but can come from diverse sources. The possibility of such inference turns diffusion models into plug-and-play modules, thereby allowing a range of potential applications in adapting models to new domains and tasks, such as conditional generation or image segmentation. The structure of diffusion models allows us to perform approximate inference by iterating differentiation through the fixed denoising network enriched with different amounts of noise at each step. Considering many noised versions of $\mathbf{x}$ in evaluation of its fitness is a novel search mechanism that may lead to new algorithms for solving combinatorial optimization problems.


title: X-ray performance of critical-angle transmission grating prototypes for the Arcus mission

abstract: Arcus is a proposed soft x-ray grating spectrometer Explorer. It aims to explore cosmic feedback by mapping hot gases within and between galaxies and galaxy clusters and characterizing jets and winds from supermassive black holes, and to investigate the dynamics of protoplanetary discs and stellar accretion. Arcus features 12 m-focal-length grazing-incidence silicon pore optics (SPO) developed for the Athena mission. Critical-angle transmission (CAT) gratings efficiently disperse high diffraction orders onto CCDs. We report new and improved x-ray performance results for Arcus-like CAT gratings, including record resolving power for two co-aligned CAT gratings. Multiple Arcus prototype grating facets were illuminated by an SPO at the PANTER facility. The facets consist of $32\times32.5$ mm$^2$ patterned silicon membranes, bonded to metal frames. The bonding angle is adjusted according to the measured average tilt angle of the grating bars in the membrane. Two simultaneously illuminated facets show minor broadening of the Al-K$_{\alpha}$ doublet in 18$^{\rm th}$ and 21$^{\rm st}$ orders with a best fit record effective resolving power of $R_G \approx 1.3^{+\infty}_{-0.5}\times10^4$ ($3\sigma$), about 3-4 times the Arcus requirement. We measured the diffraction efficiency of quasi-fully illuminated gratings at O-K wavelengths in orders 4-7 in an Arcus-like configuration and compare results with synchrotron spot measurements. After corrections for geometrical effects and bremsstrahlung continuum we find agreement between full and spot illumination at the two different facilities, as well as with the models used for Arcus effective area predictions. We find that these flight-like gratings meet diffraction efficiency and greatly exceed resolving power Arcus requirements.


title: Unifying interval maps and branching systems with applications to relative graph C*-algebras

abstract: We describe Markov interval maps via branching systems and develop the theory of relative branching systems, characterizing when the associated representations of relative graph C*-algebras are faithful. When the Markov interval maps $f$ have escape sets, we use our results to characterize injectivity of the associated relative graph algebra representations, improving on previous work by the first, third, and fourth authors.


title: Galactic Winds across the Gas-Rich Merger Sequence II. Lyman Alpha Emission and Highly Ionized O VI and N V Outflows in Ultraluminous Infrared Galaxies

abstract: This paper is the second in a series aimed at examining the gaseous environments of z$\le$0.3 quasars and ultraluminous infrared galaxies (ULIRGs) as a function of AGN/host galaxy properties across the merger sequence. This second paper focuses on the Ly$\alpha$ emission and O VI and N V absorption features, tracers of highly ionized gas outflows, in ULIRGs observed with HST/COS. Ly$\alpha$ emission is detected in 15 out of 19 ULIRGs, and 12 of the 14 clear Ly$\alpha$ detections show emission with blueshifted velocity centroids and/or wings. The equivalent widths of the Ly$\alpha$ emission increase with increasing AGN luminosities and AGN bolometric fractions. The blueshifts of the Ly$\alpha$ emission correlate positively with those of the [O III] emission, where the latter traces the ionized gas outflows. The Ly$\alpha$ escape fractions tend to be slightly larger in objects with stronger AGN and larger outflow velocities, but they do not correlate with nebular line reddening. Among the 12 ULIRGs with good continuum signal-to-noise ratios, O VI and/or N V absorption features are robustly detected in 6 of them, all of which are blueshifted, indicative of outflows. In the combined ULIRG $+$ quasar sample, the outflows are more frequently detected in the X-ray weak or absorbed sources. The absorption equivalent widths, velocities and velocity dispersions of the outflows are also higher in the X-ray weak sources. No other strong correlations are visible between the properties of the outflows and those of the AGN or host galaxies.


title: Path-Gradient Estimators for Continuous Normalizing Flows

abstract: Recent work has established a path-gradient estimator for simple variational Gaussian distributions and has argued that the path-gradient is particularly beneficial in the regime in which the variational distribution approaches the exact target distribution. In many applications, this regime can however not be reached by a simple Gaussian variational distribution. In this work, we overcome this crucial limitation by proposing a path-gradient estimator for the considerably more expressive variational family of continuous normalizing flows. We outline an efficient algorithm to calculate this estimator and establish its superior performance empirically.


title: Highly spin-polarized multi-GeV electron beams generated by single-species plasma photocathodes

abstract: High-gradient and high-efficiency acceleration in plasma-based accelerators has been demonstrated, showing its potential as the building block for a future collider operating at the energy frontier of particle physics. However, generating and accelerating the required spin-polarized beams in such a collider using plasma-based accelerators has been a long-standing challenge. Here we show that the passage of a highly relativistic, high-current electron beam through a single-species (ytterbium) vapor excites a nonlinear plasma wake by primarily ionizing the two outer 6s electrons. Further photoionization of the resultant Yb2+ ions by a circularly polarized laser injects the 4f14 electrons into this wake generating a highly spin-polarized beam. Combining time-dependent Schrodinger equation simulations with particle-in-cell simulations, we show that a sub-femtosecond, high-current (4 kA) electron beam with up to 56% net spin polarization can be generated and accelerated to 15 GeV in just 41 cm. This relatively simple scheme solves the perplexing problem of producing spin-polarized relativistic electrons in plasma-based accelerators.

categories: physics.plasm-ph physics.acc-ph physics.atom-ph physics.optics


title: Bayesian Data Augmentation for Partially Observed Stochastic Compartmental Models

abstract: Deterministic compartmental models are predominantly used in the modeling of infectious diseases, though stochastic models are considered more realistic, yet are complicated to estimate due to missing data. In this paper we present a novel algorithm for estimating the stochastic SIR/SEIR epidemic model within a Bayesian framework, which can be readily extended to more complex stochastic compartmental models. Specifically, based on the infinitesimal conditional independence properties of the model, we are able to find a proposal distribution for a Metropolis algorithm which is very close to the correct posterior distribution. As a consequence, rather than perform a Metropolis step updating one missing data point at a time, as in the current benchmark Markov chain Monte Carlo (MCMC) algorithm, we are able to extend our proposal to the entire set of missing observations. This improves the MCMC methods dramatically and makes the stochastic models now a viable modeling option. A number of real data illustrations and the necessary mathematical theory supporting our results are presented.


title: A unified extension theory of Rota-Baxter algebras, dendriform algebras, and a fundamental sequence of Wells

abstract: A Rota-Baxter algebra $A_R$ is an algebra $A$ equipped with a distinguished Rota-Baxter operator $R$ on it. Rota-Baxter algebras are closely related to dendriform algebras introduced by Loday. In this paper, we first consider the non-abelian extension theory of Rota-Baxter algebras and classify them by introducing the non-abelian cohomology. Next, given a non-abelian extension $0 \rightarrow B_S \rightarrow E_U \rightarrow A_R \rightarrow 0$ of Rota-Baxter algebras, we construct the Wells type exact sequences and find their role in extending a Rota-Baxter automorphism $\beta \in \mathrm{Aut}(B_S)$ and lifting a Rota-Baxter automorphism $\alpha \in \mathrm{Aut}(A_R)$ to an automorphism in $\mathrm{Aut}(E_U)$. We end this paper by considering a similar study for dendriform algebras.


title: Uniform and Modular Sequent Systems for Description Logics

abstract: We introduce a framework that allows for the construction of sequent systems for expressive description logics extending ALC. Our framework not only covers a wide array of common description logics, but also allows for sequent systems to be obtained for extensions of description logics with special formulae that we call "role relational axioms." All sequent systems are sound, complete, and possess favorable properties such as height-preserving admissibility of common structural rules and height-preserving invertibility of rules.


title: Conditional Permutation Invariant Flows

abstract: We present a novel, conditional generative probabilistic model of set-valued data with a tractable log density. This model is a continuous normalizing flow governed by permutation equivariant dynamics. These dynamics are driven by a learnable per-set-element term and pairwise interactions, both parametrized by deep neural networks. We illustrate the utility of this model via applications including (1) complex traffic scene generation conditioned on visually specified map information, and (2) object bounding box generation conditioned directly on images. We train our model by maximizing the expected likelihood of labeled conditional data under our flow, with the aid of a penalty that ensures the dynamics are smooth and hence efficiently solvable. Our method significantly outperforms non-permutation invariant baselines in terms of log likelihood and domain-specific metrics (offroad, collision, and combined infractions), yielding realistic samples that are difficult to distinguish from real data.


title: Designing MacPherson Suspension Architectures using Bayesian Optimization

abstract: Engineering design is traditionally performed by hand: an expert makes design proposals based on past experience, and these proposals are then tested for compliance with certain target specifications. Testing for compliance is performed first by computer simulation using what is called a discipline model. Such a model can be implemented by a finite element analysis, multibody systems approach, etc. Designs passing this simulation are then considered for physical prototyping. The overall process may take months, and is a significant cost in practice. We have developed a Bayesian optimization system for partially automating this process by directly optimizing compliance with the target specification with respect to the design parameters. The proposed method is a general framework for computing a generalized inverse of a high-dimensional non-linear function that does not require e.g. gradient information, which is often unavailable from discipline models. We furthermore develop a two-tier convergence criterion based on (i) convergence to a solution optimally satisfying all specified design criteria, or (ii) convergence to a minimum-norm solution. We demonstrate the proposed approach on a vehicle chassis design problem motivated by an industry setting using a state-of-the-art commercial discipline model. We show that the proposed approach is general, scalable, and efficient, and that the novel convergence criteria can be implemented straightforwardly based on existing concepts and subroutines in popular Bayesian optimization software packages.


title: Efficient Object Manipulation Planning with Monte Carlo Tree Search

abstract: This paper presents an efficient approach to object manipulation planning using Monte Carlo Tree Search (MCTS) to find contact sequences and an efficient ADMM-based trajectory optimization algorithm to evaluate the dynamic feasibility of candidate contact sequences. To accelerate MCTS, we propose a methodology to learn a goal-conditioned policy-value network to direct the search towards promising nodes. Further, manipulation-specific heuristics enable to drastically reduce the search space. Systematic object manipulation experiments in a physics simulator and on real hardware demonstrate the efficiency of our approach. In particular, our approach scales favorably for long manipulation sequences thanks to the learned policy-value network, significantly improving planning success rate.


title: Partisan US News Media Representations of Syrian Refugees

abstract: We investigate how representations of Syrian refugees (2011-2021) differ across US partisan news outlets. We analyze 47,388 articles from the online US media about Syrian refugees to detail differences in reporting between left- and right-leaning media. We use various NLP techniques to understand these differences. Our polarization and question answering results indicated that left-leaning media tended to represent refugees as child victims, welcome in the US, and right-leaning media cast refugees as Islamic terrorists. We noted similar results with our sentiment and offensive speech scores over time, which detail possibly unfavorable representations of refugees in right-leaning media. A strength of our work is how the different techniques we have applied validate each other. Based on our results, we provide several recommendations. Stakeholders may utilize our findings to intervene around refugee representations, and design communications campaigns that improve the way society sees refugees and possibly aid refugee outcomes.


title: Disentangling the Signatures of Blended-Light Atmospheres in L/T Transition Brown Dwarfs

abstract: We present a technique to identify spectrophotometrically variable L7-T3 brown dwarfs with single-epoch, low-resolution, near-infrared SpeX spectra. We calculated spectral indices on known variable brown dwarfs and used them to select 11 index-index parameter spaces where known variables can be distinguished from the rest of the general population of brown dwarfs. We find 62 candidate variables, 12 of which show significant variability amplitude in independent photometric monitoring surveys. This technique constitutes the first formal method to identify a time-dependent effect such as variability from peculiarities in their integrated light spectra. This technique will be a useful tool to prioritize targets for future photometric and spectroscopic monitoring in the era of JWST and 30m-class telescopes.


title: Gluon mass from charmonium radiative decay and its impact on the hybrid meson properties

abstract: We study the photon spectrum in the radiative decays of the charmonium $\frac{1}{\Gamma}\frac{d\Gamma (J/\psi\rightarrow\gamma gg)}{dz}$ as an excellent scene for hunting the exotic mesons (hybrids $q\bar{q}g$ and glueballs $gg$) and an ideal frame for the theoretical and phenomenological investigations of the gluon nature in IR, by subtle estimation of its non-zero effective mass which seems required in Quark-Gluon Constituent Model (QGCM) for identifying the nature of these exotic species.


title: Landscape Learning for Neural Network Inversion

abstract: Many machine learning methods operate by inverting a neural network at inference time, which has become a popular technique for solving inverse problems in computer vision, robotics, and graphics. However, these methods often involve gradient descent through a highly non-convex loss landscape, causing the optimization process to be unstable and slow. We introduce a method that learns a loss landscape where gradient descent is efficient, bringing massive improvement and acceleration to the inversion process. We demonstrate this advantage on a number of methods for both generative and discriminative tasks, including GAN inversion, adversarial defense, and 3D human pose reconstruction.


title: The LCO Outbursting Objects Key Project: Overview and Year 1 Status

abstract: The LCO Outbursting Objects Key (LOOK) Project uses the telescopes of the Las Cumbres Observatory (LCO) Network to: (1) to systematically monitor a sample of Dynamically New Comets over the whole sky, and (2) use alerts from existing sky surveys to rapidly respond to and characterize detected outburst activity in all small bodies. The data gathered on outbursts helps to characterize each outburst's evolution with time, assess the frequency and magnitude distribution of outbursts in general, and contributes to the understanding of outburst processes and volatile distribution in the Solar System. The LOOK Project exploits the synergy between current and future wide-field surveys such as ZTF, PanSTARRS, and LSST as well as rapid-response telescope networks such as LCO, and serves as an excellent testbed for what will be needed the much larger number of objects coming from Rubin Observatory. We will describe the LOOK Project goals, the planning and target selection (including the use of NEOexchange as a Target and Observation Manager or "TOM"), and results from the first phase of observations, including the detection of activity and outbursts on the giant comet C/2014 UN271 (Bernardinelli-Bernstein) and the discovery and follow-up of outbursts on comets. Within these outburst discoveries, we present a high cadence of 7P/Pons-Winnecke with days, a large outburst on 57P/duToit-Neujmin-Delporte, and evidence that comet P/2020 X1 (ATLAS) was in outburst when discovered.


title: Binary Early-Exit Network for Adaptive Inference on Low-Resource Devices

abstract: Deep neural networks have significantly improved performance on a range of tasks with the increasing demand for computational resources, leaving deployment on low-resource devices (with limited memory and battery power) infeasible. Binary neural networks (BNNs) tackle the issue to an extent with extreme compression and speed-up gains compared to real-valued models. We propose a simple but effective method to accelerate inference through unifying BNNs with an early-exiting strategy. Our approach allows simple instances to exit early based on a decision threshold and utilizes output layers added to different intermediate layers to avoid executing the entire binary model. We extensively evaluate our method on three audio classification tasks and across four BNNs architectures. Our method demonstrates favorable quality-efficiency trade-offs while being controllable with an entropy-based threshold specified by the system user. It also results in better speed-ups (latency less than 6ms) with a single model based on existing BNN architectures without retraining for different efficiency levels. It also provides a straightforward way to estimate sample difficulty and better understanding of uncertainty around certain classes within the dataset.


title: Responsibility-associated Multi-agent Collision Avoidance with Social Preferences

abstract: This paper introduces a novel social preference-aware decentralized safe control framework to address the responsibility allocation problem in multi-agent collision avoidance. Considering that agents do not necessarily cooperate in symmetric ways, this paper focuses on semi-cooperative behavior among heterogeneous agents with varying cooperation levels. Drawing upon the idea of Social Value Orientation (SVO) for quantifying the individual selfishness, we propose a novel concept of Responsibility-associated Social Value Orientation (R-SVO) to express the intended relative social implications between pairwise agents. This is used to redefine each agent's social preferences or personalities in terms of corresponding responsibility shares in contributing to the coordination scenario, such as semi-cooperative collision avoidance where all agents interact in an asymmetric way. By incorporating such relative social implications through proposed Local Pairwise Responsibility Weights, we develop a Responsibility-associated Control Barrier Function-based safe control framework for individual agents, and multi-agent collision avoidance is achieved with formally provable safety guarantees. Simulations are provided to demonstrate the effectiveness and efficiency of the proposed framework in several multi-agent navigation tasks, such as a position-swapping game, a self-driving car highway ramp merging scenario, and a circular position swapping game.


title: Revisiting Viscous Transonic Decretion Disks of Be Stars

abstract: In the context of Be stars, we restudied the viscous transonic decretion disk model of these stars. This model is driven by a radiative force due to an ensemble of optically-thin lines and viscosity considering the Shakura Sunyaev prescription. The non-linear equation of motion presents a singularity (sonic point) and an eigenvalue, which is also the initial condition at the stellar surface. Then, to obtain this eigenvalue, we set it as a radial quantity and perform a detailed topological analysis. Thereafter, we describe a numerical method for solving either Nodal and Saddle transonic solutions. The value of the viscosity,"alpha", barely determine the location of the sonic point, but it determines the topology of the solution. We found two Nodal solutions, which are almost indistinguishable between them. Saddle solutions are founded for lower values of "alpha" than the required of the Nodal solutions. In addition, rotational velocity do not play a determine role in the velocity (and density) profile, because viscosity effects collapse all the solutions to almost a unique one in a small region above the stellar surface. A suitable combination of line force parameters and/or disk temperature, give location of the sonic point lower than 50 stellar radii, describing a truncated disk. This could explain the SED turndown observed in Be stars without needing a binary companion.


title: Conjunctive Queries with Free Access Patterns under Updates

abstract: We study the problem of answering conjunctive queries with free access patterns (CQAP) under updates. A free access pattern is a partition of the free variables of the query into input and output. The query returns tuples over the output variables given a tuple of values over the input variables. We introduce a fully dynamic evaluation approach for CQAP queries. We also give a syntactic characterisation of those CQAP queries that admit constant time per single-tuple update and whose output tuples can be enumerated with constant delay given a tuple of values over the input variables. Finally, we chart the complexity trade-off between the preprocessing time, update time and enumeration delay for CQAP queries. For a class of CQAP queries, our approach achieves optimal, albeit non-constant, update time and delay. Their optimality is predicated on the Online Matrix-Vector Multiplication conjecture. Our results recover prior work on the dynamic evaluation of conjunctive queries without access patterns. We also illustrate an application of our dynamic evaluation approach to tractable CQAP queries over probabilistic databases.


title: Magnetic quadrupole transitions in the relativistic energy density functional theory

abstract: Background: Magnetic quadrupole (M2) excitation represents a fundamental feature in atomic nucleus associated to nuclear magnetism induced by spin and orbital transition operator. So far it has only been investigated within the non-relativistic theoretical approaches, and available experimental data are rather limited. Purpose: We aim to investigate the properties of M2 transitions in closed and open-shell nuclei using the framework of relativistic nuclear energy density functional. The calculated M2 transition strengths could be used to constrain the quenching of the spin gyromagnetic factors. Methods: The M2 excitations are described using the relativistic quasiparticle random phase approximation (RQRPA) with the residual interaction extended with the isovector-pseudovector term. Results: The M2 transition strength distributions are described and analyzed for closed and open shell nuclei. The results are compared with available experimental data and the strength missing from the experiment is discussed. The evolution of M2 transition properties has been investigated within the $^{36-64} \rm Ca$ isotope chain. Conclusion: The main M2 transitions have rather rich underlying structure and their collectivity increases with the mass number due to larger number of contributing particle-hole configurations. Pairing correlations in open shell nuclei have strong effect, causing the M2 strength reduction and shifting of the centroid energies to higher values. The analysis of M2 transition strengths indicate that considerable amount of experimental strength may be missing, mainly due to limitations to rather restricted energy ranges. The calculated M2 strengths for Ca isotopes, together with the future experimental data will allow constraining the quenching of the $g$ factors in nuclear medium.


title: Towards Better Selective Classification

abstract: We tackle the problem of Selective Classification where the objective is to achieve the best performance on a predetermined ratio (coverage) of the dataset. Recent state-of-the-art selective methods come with architectural changes either via introducing a separate selection head or an extra abstention logit. In this paper, we challenge the aforementioned methods. The results suggest that the superior performance of state-of-the-art methods is owed to training a more generalizable classifier rather than their proposed selection mechanisms. We argue that the best performing selection mechanism should instead be rooted in the classifier itself. Our proposed selection strategy uses the classification scores and achieves better results by a significant margin, consistently, across all coverages and all datasets, without any added compute cost. Furthermore, inspired by semi-supervised learning, we propose an entropy-based regularizer that improves the performance of selective classification methods. Our proposed selection mechanism with the proposed entropy-based regularizer achieves new state-of-the-art results.


title: Homological Mirror Symmetry for the universal centralizers

abstract: We prove homological mirror symmetry for the universal centralizer $J_G$ (a.k.a the Toda space), associated to any complex reductive Lie group $G$. The A-side is a partially wrapped Fukaya category on $J_G$, and the B-side is the category of coherent sheaves on the categorical quotient of a dual maximal torus by the Weyl group action (with some modification if the center of $G$ is not connected).


title: A Note on Monte Carlo Integration in High Dimensions

abstract: Monte Carlo integration is a commonly used technique to compute intractable integrals and is typically thought to perform poorly for very high-dimensional integrals. To show that this is not always the case, we examine Monte Carlo integration using techniques from the high-dimensional statistics literature by allowing the dimension of the integral to increase. In doing so, we derive non-asymptotic bounds for the relative and absolute error of the approximation for some general classes of functions through concentration inequalities. We provide concrete examples in which the magnitude of the number of points sampled needed to guarantee a consistent estimate vary between polynomial to exponential, and show that in theory arbitrarily fast or slow rates are possible. This demonstrates that the behaviour of Monte Carlo integration in high dimensions is not uniform. Through our methods we also obtain non-asymptotic confidence intervals for Monte Carlo estimate which are valid regardless of the number of points sampled.


title: Assessing transportation accessibility equity via open data

abstract: We propose a methodology to assess transportation accessibility inequity in metropolitan areas. The methodology is based on the classic analysis tools of Lorenz curves and Gini indices, but the novelty resides in the fact that it can be easily applied in an automated way to several cities around the World, with no need for customized data treatment. Indeed, our equity metrics can be computed solely relying on open data, publicly available in standardized form. We showcase our method and study transp


title: Validation of Vector Data using Oblique Images

abstract: Oblique images are aerial photographs taken at oblique angles to the earth's surface. Projections of vector and other geospatial data in these images depend on camera parameters, positions of the geospatial entities, surface terrain, occlusions, and visibility. This paper presents a robust and scalable algorithm to detect inconsistencies in vector data using oblique images. The algorithm uses image descriptors to encode the local appearance of a geospatial entity in images. These image descriptors combine color, pixel-intensity gradients, texture, and steerable filter responses. A Support Vector Machine classifier is trained to detect image descriptors that are not consistent with underlying vector data, digital elevation maps, building models, and camera parameters. In this paper, we train the classifier on visible road segments and non-road data. Thereafter, the trained classifier detects inconsistencies in vectors, which include both occluded and misaligned road segments. The consistent road segments validate our vector, DEM, and 3-D model data for those areas while inconsistent segments point out errors. We further show that a search for descriptors that are consistent with visible road segments in the neighborhood of a misaligned road yields the desired road alignment that is consistent with pixels in the image.


title: Spectropolarimetry of the tidal disruption event AT 2019qiz: a quasispherical reprocessing layer

abstract: We present optical spectropolarimetry of the tidal disruption event (TDE) AT 2019qiz on days $+0$ and $+29$ relative to maximum brightness. Continuum polarization, which informs the shape of the electron-scattering surface, was found to be consistent with 0 per cent at peak brightness. On day $+29$, the continuum polarization rose to $\sim 1$ per cent, making this the first reported spectropolarimetric evolution of a TDE. These findings are incompatible with a naked eccentric disc that lacks significant mass outflow. Instead, the spectropolarimetry paints a picture wherein, at maximum brightness, high-frequency emission from the accretion disc is reprocessed into the optical band by a nearly spherical, optically thick, electron-scattering photosphere located far away from the black hole. We estimate the radius of the scattering photosphere to be $\sim 100\rm\, au$ at maximum brightness -- significantly larger than the tidal radius ($\sim 1\rm\, au$) and the thermalisation radius ($\sim 30\rm\, au$) where the optical continuum is formed. A month later, as the fallback rate drops and the scattering photosphere recedes, the continuum polarization increases, revealing a moderately aspherical interior. We also see evidence for smaller-scale density variations in the scattering photosphere, inferred from the scatter of the data in the Stokes $q-u$ plane. On day $+29$, the H$\alpha$ emission-line peak is depolarized to $\sim 0.3$ per cent (compared to $\sim 1$ per cent continuum polarization), and displays a gradual rise toward the line's redder wavelengths. This observation indicates the H$\alpha$ line formed near the electron-scattering radius.


title: Energy reconstruction for large liquid scintillator detectors with machine learning techniques: aggregated features approach

abstract: Large-scale detectors consisting of a liquid scintillator target surrounded by an array of photo-multiplier tubes (PMTs) are widely used in the modern neutrino experiments: Borexino, KamLAND, Daya Bay, Double Chooz, RENO, and the upcoming JUNO with its satellite detector TAO. Such apparatuses are able to measure neutrino energy which can be derived from the amount of light and its spatial and temporal distribution over PMT channels. However, achieving a fine energy resolution in large-scale detectors is challenging. In this work, we present machine learning methods for energy reconstruction in the JUNO detector, the most advanced of its type. We focus on positron events in the energy range of 0-10 MeV which corresponds to the main signal in JUNO -- neutrinos originated from nuclear reactor cores and detected via the inverse beta decay channel. We consider the following models: Boosted Decision Trees and Fully Connected Deep Neural Network, trained on aggregated features, calculated using the information collected by PMTs. We describe the details of our feature engineering procedure and show that machine learning models can provide the energy resolution $\sigma = 3\%$ at 1 MeV using subsets of engineered features. The dataset for model training and testing is generated by the Monte Carlo method with the official JUNO software.


title: Accelerating Machine Learning Training Time for Limit Order Book Prediction

abstract: Financial firms are interested in simulation to discover whether a given algorithm involving financial machine learning will operate profitably. While many versions of this type of algorithm have been published recently by researchers, the focus herein is on a particular machine learning training project due to the explainable nature and the availability of high frequency market data. For this task, hardware acceleration is expected to speed up the time required for the financial machine learning researcher to obtain the results. As the majority of the time can be spent in classifier training, there is interest in faster training steps. A published Limit Order Book algorithm for predicting stock market direction is our subject, and the machine learning training process can be time-intensive especially when considering the iterative nature of model development. To remedy this, we deploy Graphical Processing Units (GPUs) produced by NVIDIA available in the data center where the computer architecture is geared to parallel high-speed arithmetic operations. In the studied configuration, this leads to significantly faster training time allowing more efficient and extensive model development.


title: Riemannian CUR Decompositions for Robust Principal Component Analysis

abstract: Robust Principal Component Analysis (PCA) has received massive attention in recent years. It aims to recover a low-rank matrix and a sparse matrix from their sum. This paper proposes a novel nonconvex Robust PCA algorithm, coined Riemannian CUR (RieCUR), which utilizes the ideas of Riemannian optimization and robust CUR decompositions. This algorithm has the same computational complexity as Iterated Robust CUR, which is currently state-of-the-art, but is more robust to outliers. RieCUR is also able to tolerate a significant amount of outliers, and is comparable to Accelerated Alternating Projections, which has high outlier tolerance but worse computational complexity than the proposed method. Thus, the proposed algorithm achieves state-of-the-art performance on Robust PCA both in terms of computational complexity and outlier tolerance.


title: Black Hole Mass Measurements of Early-Type Galaxies NGC 1380 and NGC 6861 Through ALMA and HST Observations and Gas-Dynamical Modeling

abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 2 observations of CO(2-1) emission from the circumnuclear disks in two early-type galaxies, NGC 1380 and NGC 6861. The disk in each galaxy is highly inclined ($i\,{\sim}\,75^{\circ}$), and the projected velocities of the molecular gas near the galaxy centers are ${\sim}300\,\mathrm{km \, s^{-1}}$ in NGC 1380 and ${\sim}500\,\mathrm{km \, s^{-1}}$ in NGC 6861. We fit thin disk dynamical models to the ALMA data cubes to constrain the masses of the central black holes (BHs). We created host galaxy models using Hubble Space Telescope images for the extended stellar mass distributions and incorporated a range of plausible central dust extinction values. For NGC 1380, our best-fit model yields $M_{\mathrm{BH}} = 1.47 \times 10^8\,M_{\odot}$ with a ${\sim}40\%$ uncertainty. For NGC 6861, the lack of dynamical tracers within the BH's sphere of influence due to a central hole in the gas distribution precludes a precise measurement of $M_{\mathrm{BH}}$. However, our model fits require a value for $M_{\mathrm{BH}}$ in the range of $(1-3) \times 10^9\,M_{\odot}$ in NGC 6861 to reproduce the observations. The BH masses are generally consistent with predictions from local BH-host galaxy scaling relations. Systematic uncertainties associated with dust extinction of the host galaxy light and choice of host galaxy mass model dominate the error budget of both measurements. Despite these limitations, the measurements demonstrate ALMA's ability to provide constraints on BH masses in cases where the BH's projected radius of influence is marginally resolved or the gas distribution has a central hole.


title: Universal Complexity Bounds Based on Value Iteration and Application to Entropy Games

abstract: We develop value iteration-based algorithms to solve in a unified manner different classes of combinatorial zero-sum games with mean-payoff type rewards. These algorithms rely on an oracle, evaluating the dynamic programming operator up to a given precision. We show that the number of calls to the oracle needed to determine exact optimal (positional) strategies is, up to a factor polynomial in the dimension, of order R/sep, where the "separation" sep is defined as the minimal difference between distinct values arising from strategies, and R is a metric estimate, involving the norm of approximate sub and super-eigenvectors of the dynamic programming operator. We illustrate this method by two applications. The first one is a new proof, leading to improved complexity estimates, of a theorem of Boros, Elbassioni, Gurvich and Makino, showing that turn-based mean-payoff games with a fixed number of random positions can be solved in pseudo-polynomial time. The second one concerns entropy games, a model introduced by Asarin, Cervelle, Degorre, Dima, Horn and Kozyakin. The rank of an entropy game is defined as the maximal rank among all the ambiguity matrices determined by strategies of the two players. We show that entropy games with a fixed rank, in their original formulation, can be solved in polynomial time, and that an extension of entropy games incorporating weights can be solved in pseudo-polynomial time under the same fixed rank condition.


title: Enabling Undergrounding of Long-Distance Transmission Lines with Low Frequency AC Technology

abstract: With increasing prevalence of severe natural hazards and the ignition of wildfires by power lines, many power system planners are considering converting overhead lines to underground cables to mitigate risks related to these events. Systems with a large proportion of underground cables can bring challenges due to the capacitance and losses of the cables, factors which may limit the potential for hardening in critical areas. Low frequency AC (LFAC) transmission solves these problems, as lowering the frequency decreases both the effects of capacitance and the losses. This paper presents a tractable frequency-dependent model for underground cables and incorporates it into the open source optimal power flow tool VariableFrequencyOPF.jl. The model and implementation are used to demonstrate the benefits of LFAC in a case study involving two multi-terminal cable upgrades with LFAC, including a comparison with HVDC. The results demonstrate the value of LFAC systems for both power flow control and reduction of losses.


title: Beyond Rewards: a Hierarchical Perspective on Offline Multiagent Behavioral Analysis

abstract: Each year, expert-level performance is attained in increasingly-complex multiagent domains, where notable examples include Go, Poker, and StarCraft II. This rapid progression is accompanied by a commensurate need to better understand how such agents attain this performance, to enable their safe deployment, identify limitations, and reveal potential means of improving them. In this paper we take a step back from performance-focused multiagent learning, and instead turn our attention towards agent behavior analysis. We introduce a model-agnostic method for discovery of behavior clusters in multiagent domains, using variational inference to learn a hierarchy of behaviors at the joint and local agent levels. Our framework makes no assumption about agents' underlying learning algorithms, does not require access to their latent states or policies, and is trained using only offline observational data. We illustrate the effectiveness of our method for enabling the coupled understanding of behaviors at the joint and local agent level, detection of behavior changepoints throughout training, discovery of core behavioral concepts, demonstrate the approach's scalability to a high-dimensional multiagent MuJoCo control domain, and also illustrate that the approach can disentangle previously-trained policies in OpenAI's hide-and-seek domain.


title: Zeeman-Hyperfine Measurements of a Pseudo-Degenerate Quadruplet in CaF$_2$:Ho$^{3+}$

abstract: We report Zeeman infra-red spectroscopy of electronic-nuclear levels of $^5$I$_8 \rightarrow ^5$I$_7$ transitions of Ho$^{3+}$ in the C$_{\rm 4v}$(F$^-$) centre in CaF$_2$ with the magnetic field along the $\langle 111\rangle$ direction of the crystal. Transitions to the lowest $^5$I$_7$ state, an isolated electronic doublet, and the next group of states, a pseudo-quadruplet consisting of a doublet and two nearby singlets, exhibit strongly non-linear Zeeman splittings and intensity variations. Simulated spectra based upon a crystal-field analysis give an excellent approximation to the data, illustrating the strong predictive ability of the parametrised crystal-field approach. Anti-crossings in the hyperfine splittings, the basis of quantum information storage in rare-earth doped insulating dielectrics, are also predicted.


title: ICLR 2022 Challenge for Computational Geometry and Topology: Design and Results

abstract: This paper presents the computational challenge on differential geometry and topology that was hosted within the ICLR 2022 workshop ``Geometric and Topological Representation Learning". The competition asked participants to provide implementations of machine learning algorithms on manifolds that would respect the API of the open-source software Geomstats (manifold part) and Scikit-Learn (machine learning part) or PyTorch. The challenge attracted seven teams in its two month duration. This paper describes the design of the challenge and summarizes its main findings.


title: Fast Sinkhorn II: Collinear Triangular Matrix and Linear Time Accurate Computation of Optimal Transport

abstract: In our previous work [arXiv:2202.10042], the complexity of Sinkhorn iteration is reduced from $O(N^2)$ to the optimal $O(N)$ by leveraging the special structure of the kernel matrix. In this paper, we explore the special structure of kernel matrices by defining and utilizing the properties of the Lower-ColLinear Triangular Matrix (L-CoLT matrix) and Upper-ColLinear Triangular Matrix (U-CoLT matrix). We prove that (1) L/U-CoLT matrix-vector multiplications can be carried out in $O(N)$ operations; (2) both families of matrices are closed under the Hadamard product and matrix scaling. These properties help to alleviate two key difficulties for reducing the complexity of the Inexact Proximal point method (IPOT), and allow us to significantly reduce the number of iterations to $O(N)$. This yields the Fast Sinkhorn II (FS-2) algorithm for accurate computation of optimal transport with low algorithm complexity and fast convergence. Numerical experiments are presented to show the effectiveness and efficiency of our approach.


title: Multisolitons are the unique constrained minimizers of the KdV conserved quantities

abstract: We consider the following variational problem: minimize the $(n+1)$st polynomial conserved quantity of KdV over $H^n(\mathbb{R})$ with the first $n$ conserved quantities constrained. Maddocks and Sachs used that $n$-solitons are local minimizers for this problem in order to prove that $n$-solitons are orbitally stable in $H^n(\mathbb{R})$. Given $n$ constraints that are attainable by an $n$-soliton, we show that there is a unique set of $n$ amplitude parameters so that the corresponding multisolitons satisfy the constraints. Moreover, we prove that these multisolitons are the unique global constrained minimizers. We then use this variational characterization to provide a new proof of the orbital stability result of Maddocks and Sachs via concentration compactness. In the case when the constraints can be attained by functions in $H^n(\mathbb{R})$ but not by an $n$-soliton, we discover new behavior for minimizing sequences.


title: EPOC Emotiv EEG Basics

abstract: This document provides some basic guidance to start working with the EPOC Emotiv neuroheadset device and describes how to use it to perform basic Brain-Computer Interface (BCI) research. A brief tutorial on how to set up the device, from its electrophysiological point of view, as well as a description and practical code to perform some basic analysis, is explained. A basic experiment is introduced to detect one of the oldest and, indeed, quite still valuable electrophysiological correlate, visual occipital alpha waves, or Berger Rhythm. An additional experiment is expounded where the power spectrum of alpha waves is reduced when a subject is affected by background cognitive disturbances. This document also briefs about the extraction of information by using the EPOC Emotiv library and also with python Emokit package. This report presents a basic guide on how to use EEGLAB and MATLAB, as well as python stack to perform the neurophysiological analysis. Finally, a basic analysis on different feature extraction and classification methods is provided.


title: Non-affine deformation of semiflexible polymer and fiber networks

abstract: Networks of semiflexible or stiff polymers such as most biopolymers are known to deform inhomogeneously when sheared. The effects of such non-affine deformation have been shown to be much stronger than for flexible polymers. To date, our understanding of non-affinity in such systems is limited to simulations or specific 2D models of athermal fibers. Here, we present an effective medium theory (EMT) for non-affine deformation of semiflexible polymer and fiber networks, which is general to both 2D and 3D and in both thermal and athermal limits. The predictions of this model are in good agreement with both prior computational and experimental results for linear elasticity. Moreover, the framework we introduce can be extended to address nonlinear elasticity and network dynamics.


title: An imminent collision monitoring system with safe stopping interventions for autonomous aerial flights

abstract: Collision avoidance requires tradeoffs in planning time horizons. Depending on the planner, safety cannot always be guaranteed in uncertain environments given map updates. To mitigate situations where the planner leads the vehicle into a state of collision or the vehicle reaches a point where no trajectories are feasible, we propose a continuous collision checking algorithm. The imminent collision checking system continuously monitors vehicle safety, and plans a safe trajectory that leads the vehicle to a stop within the observed map. We test our proposed pipeline alongside a teleoperated navigation in real-life experiments, and in simulated random-forest and warehouse environments where we show that with our method, we are able to mitigate collisions with a success rate of at least 90\%.


title: Learning the parameters of a differential equation from its trajectory via the adjoint equation

abstract: The paper contributes to strengthening the relation between machine learning and the theory of differential equations. In this context, the inverse problem of fitting the parameters, and the initial condition of a differential equation to some measurements constitutes a key issue. The paper explores an abstraction that can be used to construct a family of loss functions with the aim of fitting the solution of an initial value problem to a set of discrete or continuous measurements. It is shown, that an extension of the adjoint equation can be used to derive the gradient of the loss function as a continuous analogue of backpropagation in machine learning. Numerical evidence is presented that under reasonably controlled circumstances the gradients obtained this way can be used in a gradient descent to fit the solution of an initial value problem to a set of continuous noisy measurements, and a set of discrete noisy measurements that are recorded at uncertain times.


title: Augmented Imagefication: A Data-driven Fault Detection Method for Aircraft Air Data Sensors

abstract: In this paper, a novel data-driven approach named Augmented Imagefication for Fault detection (FD) of aircraft air data sensors (ADS) is proposed. Exemplifying the FD problem of aircraft air data sensors, an online FD scheme on edge device based on deep neural network (DNN) is developed. First, the aircraft inertial reference unit measurements is adopted as equivalent inputs, which is scalable to different aircraft/flight cases. Data associated with 6 different aircraft/flight conditions are collected to provide diversity (scalability) in the training/testing database. Then Augmented Imagefication is proposed for the DNN-based prediction of flying conditions. The raw data are reshaped as a grayscale image for convolutional operation, and the necessity of augmentation is analyzed and pointed out. Different kinds of augmented method, i.e. Flip, Repeat, Tile and their combinations are discussed, the result shows that the All Repeat operation in both axes of image matrix leads to the best performance of DNN. The interpretability of DNN is studied based on Grad-CAM, which provide a better understanding and further solidifies the robustness of DNN. Next the DNN model, VGG-16 with augmented imagefication data is optimized for mobile hardware deployment. After pruning of DNN, a lightweight model (98.79% smaller than original VGG-16) with high accuracy (slightly up by 0.27%) and fast speed (time delay is reduced by 87.54%) is obtained. And the hyperparameters optimization of DNN based on TPE is implemented and the best combination of hyperparameters is determined (learning rate 0.001, iterative epochs 600, and batch size 100 yields the highest accuracy at 0.987). Finally, a online FD deployment based on edge device, Jetson Nano, is developed and the real time monitoring of aircraft is achieved. We believe that this method is instructive for addressing the FD problems in other similar fields.


title: Sub-recoil clock-transition laser cooling enabling shallow optical lattice clocks

abstract: Laser cooling is a key ingredient for quantum control of atomic systems in a variety of settings. In divalent atoms, two-stage Doppler cooling is typically used to bring atoms to the uK regime. Here, we implement a pulsed radial cooling scheme using the ultranarrow 1S0-3P0 clock transition in ytterbium to realize sub-recoil temperatures, down to tens of nK. Together with sideband cooling along the one-dimensional lattice axis, we efficiently prepare atoms in shallow lattices at an energy of 6 lattice recoils. Under these conditions key limits on lattice clock accuracy and instability are reduced, opening the door to dramatic improvements. Furthermore, tunneling shifts in the shallow lattice do not compromise clock accuracy at the 10-19 level.


title: Parameter Identifiability of a Multitype Pure-Birth Model of Speciation

abstract: Diversification models describe the random growth of evolutionary trees, modeling the historical relationships of species through speciation and extinction events. One class of such models allows for independently changing traits, or types, of the species within the tree, upon which speciation and extinction rates depend. Although identifiability of parameters is necessary to justify parameter estimation with a model, it has not been formally established for these models, despite their adoption for inference. This work establishes generic identifiability up to label swapping for the parameters of one of the simpler forms of such a model, a multitype pure birth model of speciation, from an asymptotic distribution derived from a single tree observation as its depth goes to infinity. Crucially for applications to available data, no observation of types is needed at any internal points in the tree, nor even at the leaves.


title: NASTAR: Noise Adaptive Speech Enhancement with Target-Conditional Resampling

abstract: For deep learning-based speech enhancement (SE) systems, the training-test acoustic mismatch can cause notable performance degradation. To address the mismatch issue, numerous noise adaptation strategies have been derived. In this paper, we propose a novel method, called noise adaptive speech enhancement with target-conditional resampling (NASTAR), which reduces mismatches with only one sample (one-shot) of noisy speech in the target environment. NASTAR uses a feedback mechanism to simulate adaptive training data via a noise extractor and a retrieval model. The noise extractor estimates the target noise from the noisy speech, called pseudo-noise. The noise retrieval model retrieves relevant noise samples from a pool of noise signals according to the noisy speech, called relevant-cohort. The pseudo-noise and the relevant-cohort set are jointly sampled and mixed with the source speech corpus to prepare simulated training data for noise adaptation. Experimental results show that NASTAR can effectively use one noisy speech sample to adapt an SE model to a target condition. Moreover, both the noise extractor and the noise retrieval model contribute to model adaptation. To our best knowledge, NASTAR is the first work to perform one-shot noise adaptation through noise extraction and retrieval.