title: Hierarchical Integrated Spatial Process Modeling of Monotone West Antarctic Snow Density Curves

abstract: Snow density estimates below the surface, used with airplane-acquired ice-penetrating radar measurements, give a site-specific history of snow water accumulation. Because it is infeasible to drill snow cores across all of Antarctica to measure snow density and because it is critical to understand how climatic changes are affecting the world's largest freshwater reservoir, we develop methods that enable snow density estimation with uncertainty in regions where snow cores have not been drilled. In inland West Antarctica, snow density increases monotonically as a function of depth, except for possible micro-scale variability or measurement error, and it cannot exceed the density of ice. We present a novel class of integrated spatial process models that allow interpolation of monotone snow density curves. For computational feasibility, we construct the space-depth process through kernel convolutions of log-Gaussian spatial processes. We discuss model comparison, model fitting, and prediction. Using this model, we extend estimates of snow density beyond the depth of the original core and estimate snow density curves where snow cores have not been drilled. Along flight lines with ice-penetrating radar, we use interpolated snow density curves to estimate recent water accumulation and find predominantly decreasing water accumulation over recent decades.


title: Single Molecule Nonlinearity in a Plasmonic Waveguide

abstract: Plasmonic waveguides offer the unique possibility to confine light far below the diffraction limit. Past room temperature experiments of single quantum emitters coupled to such waveguides have mainly focused on efficient channelling of the spontaneous emission into waveguide modes after a linear optical excitation. However, only the simultaneous interaction of the emitter with multiple plasmonic fields would lead to functionality in a plasmonic circuit. Here, we demonstrate the nonlinear optical interaction of a single molecule and propagating plasmons. An individual terrylene diimide (TDI) molecule is placed in the nanogap between two single-crystalline silver nanowires. A visible wavelength pump pulse and a red-shifted depletion pulse travel along the waveguide, leading to stimulated emission depletion (STED) in the observed fluorescence. The efficiency increases by up to a factor of 50 compared to far-field excitation. Our study thus demonstrates remote nonlinear four-wave mixing at a single molecule with propagating plasmons. It paves the way towards functional quantum plasmonic circuits and improved nonlinear single-molecule spectroscopy.


title: Minimality of the action on the universal circle of uniform foliations

abstract: Given a uniform foliation by Gromov hyperbolic leaves on a $3$-manifold, we show that the action of the fundamental group on the universal circle is minimal and transitive on pairs of different points. We also prove two other results: we prove that general uniform Reebless foliations are $\mathbb{R}$-covered and we give a new description of the universal circle of $\mathbb{R}$-covered foliations with Gromov hyperbolic leaves in terms of the JSJ decomposition of $M$.


title: Fast large-scale boundary element algorithms

abstract: Boundary element methods (BEM) reduce a partial differential equation in a domain to an integral equation on the domain's boundary. They are particularly attractive for solving problems on unbounded domains, but handling the dense matrices corresponding to the integral operators requires efficient algorithms. This article describes two approaches that allow us to solve boundary element equations on surface meshes consisting of several millions of triangles while preserving the optimal convergence rates of the Galerkin discretization.


title: Flat histogram method comparison on 2D Ising Model

abstract: We compare the convergence of several flat-histogram methods applied to the 2D Ising model, including the recently introduced stochastic approximation with a dynamic update factor (SAD) method. We compare this method with the Wang-Landau (WL) method, the $1/t$ variant of the WL method, and standard stochastic approximation Monte Carlo (SAMC). In addition, we consider a procedure WL followed by a "production run" with fixed weights that refines the estimation of the entropy. To our knowledge, this work is the first to test this approach against other methods. We find that WL followed by a production run \emph{does} converge to the true density of states, in contrast to pure WL. Three of the methods converge robustly: SAD, $1/t$-WL, and WL followed by a production run. Of these, SAD does not require \emph{a priori} knowledge of the energy range. This work also shows that WL followed by a production run performs superior to other forms of WL while ensuring both ergodicity and detailed balance.


title: Scaling Blockchains to Support Electronic Health Records for Hospital Systems

abstract: Electronic Health Records (EHRs) have improved many aspects of healthcare and allowed for easier patient management for medical providers. Blockchains have been proposed as a promising solution for supporting Electronic Health Records (EHRs), but have also been linked to scalability concerns about supporting real-world healthcare systems. This paper quantifies the scalability issues and bottlenecks related to current blockchains and puts into perspective the limitations blockchains have with supporting healthcare systems. Particularly we show that well known blockchains such as Bitcoin, Ethereum, and IOTA cannot support transactions of a large scale hospital system such as the University of Kentucky HealthCare system and leave over 7.5M unsealed transactions per day. We then discuss how bottlenecks of blockchains can be relieved with sidechains, enabling well-known blockchains to support even larger hospital systems of over 30M transactions per day. We then introduce the Patient-Healthchain architecture to provide future direction on how scaling blockchains for EHR systems with sidechains can be achieved.


title: Unboundedness of potential dependent Riesz transforms for totally irregular measures

abstract: We prove that, for totally irregular measures $\mu$ on $\mathbb{R}^{d}$ with $d\geq3$, the $(d-1)$-dimensional Riesz transform $$ T_{A,\mu}^{V}f(x) = \int_{\mathbb{R}^d} \nabla_{1}\mathcal{E}_{A}^{V}(x,y) f(y) \, d \mu(y) $$ adapted to the Schr\"{o}dinger operator $L_{A}^{V} = -\mathrm{div} A \nabla + V$ with fundamental solution $\mathcal{E}_{A}^{V}$ is not bounded on $L^{2}(\mu)$. This generalises recent results obtained by Conde-Alonso, Mourgoglou and Tolsa for free-space elliptic operators with H\"older continuous coefficients $A$ since it allows for the presence of potentials $V$ in the reverse H\"{o}lder class $RH_{d}$. We achieve this by obtaining new exponential decay estimates for the kernel $\nabla_{1} \mathcal{E}_{A}^{V}$ as well as H\"older regularity estimates at local scales determined by the potential's critical radius function.


title: Robust preconditioning of monolithically coupled multiphysics problems

abstract: In many applications, one wants to model physical systems consisting of two different physical processes in two different domains that are coupled across a common interface. A crucial challenge is then that the solutions of the two different domains often depend critically on the interaction at the interface and therefore the problem cannot be easily decoupled into its subproblems. Here, we present a framework for finding robust preconditioners for a fairly general class of such problems by exploiting operators representing fractional and weighted Laplacians at the interface. Furthermore, we show feasibility of the framework for two common multiphysics problems; namely the Darcy-Stokes problem and a fluid--structure interaction problem. Numerical experiments that demonstrate the effectiveness of the approach are included.


title: Cell-Type-Specific Differences in KDEL Receptor Clustering in Mammalian Cells

abstract: In eukaryotic cells, KDEL receptors (KDELRs) facilitate the retrieval of endoplasmic reticulum (ER) luminal proteins from the Golgi compartment back to the ER. Apart from the well-documented retention function, recent findings reveal that the cellular KDELRs have more complex roles, e.g. in cell signalling, protein secretion, cell adhesion and tumorigenesis. Furthermore, several studies suggest that a sub-population of KDELRs is located at the cell surface, where they could form and internalize KDELR/cargo clusters after K/HDEL-ligand binding. However, so far it has been unclear whether there are cell-type- or species-specific differences in KDELR clustering. By comparing ligand-induced KDELR clustering in different mouse and human cell lines via live cell imaging, we show that macrophage cell lines from both species do not develop any clusters. Using RT-qPCR experiments and numerical analysis, we address the role of KDELR expression as well as endocytosis and exocytosis rates on the receptor clustering at the plasma membrane and discuss how the efficiency of directed transport to preferred docking sites on the membrane influences the exponent of the power-law distribution of the cluster size.


title: Robust preconditioning for coupled Stokes-Darcy problems with the Darcy problem in primal form

abstract: The coupled Darcy-Stokes problem is widely used for modeling fluid transport in physical systems consisting of a porous part and a free part. In this work we consider preconditioners for monolitic solution algorithms of the coupled Darcy-Stokes problem, where the Darcy problem is in primal form. We employ the operator preconditioning framework and utilize a fractional solver at the interface between the problems to obtain order optimal schemes that are robust with respect to the material parameters, i.e. the permeability, viscosity and Beavers-Joseph-Saffman parameter. Our approach is similar to our earlier work, but since the Darcy problem is in primal form, the mass conservation at the interface introduces some challenges. These challenges will be specifically addressed in this paper. Numerical experiments illustrating the performance are provided. The preconditioner is posed in non-standard Sobolev spaces which may be perceived as an obstacle for its use in applications. However, we detail the implementational aspects and show that the preconditioner is quite feasible to realize in practice.


title: Biorthogonal Greedy Algorithms in Convex Optimization

abstract: The study of greedy approximation in the context of convex optimization is becoming a promising research direction as greedy algorithms are actively being employed to construct sparse minimizers for convex functions with respect to given sets of elements. In this paper we propose a unified way of analyzing a certain kind of greedy-type algorithms for the minimization of convex functions on Banach spaces. Specifically, we define the class of Weak Biorthogonal Greedy Algorithms for convex optimization that contains a wide range of greedy algorithms. We analyze the introduced class of algorithms and establish the properties of convergence, rate of convergence, and numerical stability, which is understood in the sense that the steps of the algorithm are allowed to be performed not precisely but with controlled computational inaccuracies. We show that the following well-known algorithms for convex optimization -- the Weak Chebyshev Greedy Algorithm (co) and the Weak Greedy Algorithm with Free Relaxation (co) -- belong to this class, and introduce a new algorithm -- the Rescaled Weak Relaxed Greedy Algorithm (co). Presented numerical experiments demonstrate the practical performance of the aforementioned greedy algorithms in the setting of convex minimization as compared to optimization with regularization, which is the conventional approach of constructing sparse minimizers.


title: Random walk algorithm for the Dirichlet problem for parabolic integro-differential equation

abstract: We consider stochastic differential equations driven by a general L\'evy processes (SDEs) with infinite activity and the related, via the Feynman-Kac formula, Dirichlet problem for parabolic integro-differential equation (PIDE). We approximate the solution of PIDE using a numerical method for the SDEs. The method is based on three ingredients: (i) we approximate small jumps by a diffusion; (ii) we use restricted jump-adaptive time-stepping; and (iii) between the jumps we exploit a weak Euler approximation. We prove weak convergence of the considered algorithm and present an in-depth analysis of how its error and computational cost depend on the jump activity level. Results of some numerical experiments, including pricing of barrier basket currency options, are presented.


title: Improving GANs for Speech Enhancement

abstract: Generative adversarial networks (GAN) have recently been shown to be efficient for speech enhancement. However, most, if not all, existing speech enhancement GANs (SEGAN) make use of a single generator to perform one-stage enhancement mapping. In this work, we propose to use multiple generators that are chained to perform multi-stage enhancement mapping, which gradually refines the noisy input signals in a stage-wise fashion. Furthermore, we study two scenarios: (1) the generators share their parameters and (2) the generators' parameters are independent. The former constrains the generators to learn a common mapping that is iteratively applied at all enhancement stages and results in a small model footprint. On the contrary, the latter allows the generators to flexibly learn different enhancement mappings at different stages of the network at the cost of an increased model size. We demonstrate that the proposed multi-stage enhancement approach outperforms the one-stage SEGAN baseline, where the independent generators lead to more favorable results than the tied generators. The source code is available at http://github.com/pquochuy/idsegan.


title: Quantum critical point in the itinerant ferromagnet Ni$_{1-x}$Rh$_x$

abstract: We report a chemical substitution-induced ferromagnetic quantum critical point in polycrystalline Ni$_{1-x}$Rh$_x$ alloys. Through magnetization and muon spin relaxation measurements, we show that the ferromagnetic ordering temperature is suppressed continuously to zero at $x_{crit} = 0.375$ while the magnetic volume fraction remains 100% up to $x_{crit}$, pointing to a second order transition. Non-Fermi liquid behavior is observed close to $x_{crit}$, where the electronic specific heat $C_{el}/T$ diverges logarithmically, while immediately above $x_{crit}$ the volume thermal expansion coefficient $\alpha_{V}/T$ and the Gr\"uneisen ratio $\Gamma = \alpha_{V}/C_{el}$ both diverge logarithmically in the low temperature limit, further indication of a ferromagnetic quantum critical point in Ni$_{1-x}$Rh$_x$.


title: An evaluation of machine learning techniques to predict the outcome of children treated for Hodgkin-Lymphoma on the AHOD0031 trial: A report from the Children's Oncology Group

abstract: In this manuscript we analyze a data set containing information on children with Hodgkin Lymphoma (HL) enrolled on a clinical trial. Treatments received and survival status were collected together with other covariates such as demographics and clinical measurements. Our main task is to explore the potential of machine learning (ML) algorithms in a survival analysis context in order to improve over the Cox Proportional Hazard (CoxPH) model. We discuss the weaknesses of the CoxPH model we would like to improve upon and then we introduce multiple algorithms, from well-established ones to state-of-the-art models, that solve these issues. We then compare every model according to the concordance index and the brier score. Finally, we produce a series of recommendations, based on our experience, for practitioners that would like to benefit from the recent advances in artificial intelligence.


title: The Bhargava greedoid as a Gaussian elimination greedoid

abstract: Inspired by Manjul Bhargava's theory of generalized factorials, Fedor Petrov and the author have defined the "Bhargava greedoid" -- a greedoid (a matroid-like set system on a finite set) assigned to any "ultra triple" (a somewhat extended variant of a finite ultrametric space). Here we show that the Bhargava greedoid of a finite ultra triple is always a "Gaussian elimination greedoid" over any sufficiently large (e.g., infinite) field; this is a greedoid analogue of a representable matroid. We find necessary and sufficient conditions on the size of the field to ensure this.


title: VVV-WIT-01: highly obscured classical nova or protostellar collision?

abstract: A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $\mu$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. The source is projected against the Infrared Dark Cloud (IRDC) SDC G331.062$-$0.294. The chance projection probability is small for any single event ($p \approx 0.01$ to 0.02) which suggests a physical association, e.g. a collision between low mass protostars. However, black body emission at $T \sim 1000$ K is common in classical novae (especially CO novae) at the infrared peak in the light curve, due to condensation of dust $\sim$30--60 days after the explosion. Radio follow up with the Australia Telescope Compact Array (ATCA) detected a fading continuum source with properties consistent with a classical nova but probably inconsistent with colliding protostars. Considering all VVV transients that could have been projected against a catalogued IRDC raises the probability of a chance association to $p=0.13$ to 0.24. After weighing several options, it appears likely that VVV-WIT-01 was a classical nova event located behind an IRDC.


title: Accelerated Dual-Averaging Primal-Dual Method for Composite Convex Minimization

abstract: Dual averaging-type methods are widely used in industrial machine learning applications due to their ability to promoting solution structure (e.g., sparsity) efficiently. In this paper, we propose a novel accelerated dual-averaging primal-dual algorithm for minimizing a composite convex function. We also derive a stochastic version of the proposed method which solves empirical risk minimization, and its advantages on handling sparse data are demonstrated both theoretically and empirically.


title: Functional Calculus on Non-Homogeneous Operators on Nilpotent Groups

abstract: We study the functional calculus associated with a hypoelliptic left-invariant differential operator $\mathcal{L}$ on a connected and simply connected nilpotent Lie group $G$ with the aid of the corresponding \emph{Rockland} operator $\mathcal{L}_0$ on the `local' contraction $G_0$ of $G$, as well as of the corresponding Rockland operator $\mathcal{L}_\infty$ on the `global' contraction $G_\infty$ of $G$. We provide asymptotic estimates of the Riesz potentials associated with $\mathcal{L}$ at $0$ and at $\infty$, as well as of the kernels associated with functions of $\mathcal{L}$ satisfying Mihlin conditions of every order. We also prove some Mihlin-H\"ormander multiplier theorems for $\mathcal{L}$ which generalize analogous results to the non-homogeneous case. Finally, we extend the asymptotic study of the density of the `Plancherel measure' associated with $\mathcal{L}$ from the case of a quasi-homogeneous sub-Laplacian to the case of a quasi-homogeneous sum of even powers.


title: Raman Response and Transport Properties of One-Dimensional van der Waals Tellurium Nanowires

abstract: Tellurium can form nanowires of helical atomic chains. Given their unique one-dimensional van der Waals structure, these nanowires are expected to show remarkably different physical and electronic properties than bulk tellurium. Here we show that few-chain and single-chain van der Waals tellurium nanowires can be isolated using carbon nanotube and boron nitride nanotube encapsulation. With the approach, the number of atomic chains can be controlled by the inner diameter of the nanotube. The Raman response of the structures suggests that the interaction between a single-atomic tellurium chain and a carbon nanotube is weak, and that the inter-chain interaction becomes stronger as the number of chains increases. Compared with bare tellurium nanowires on SiO2, nanowires encapsulated in boron nitride nanotubes exhibit a dramatically enhanced current-carrying capacity, with a current density of 1.5*10^8 A cm-2, which exceeds that of most semiconducting nanowires. We also use our tellurium nanowires encapsulated in boron nitride nanotubes to create field-effect transistors that have a diameter of only 2 nm.


title: Insertion-Deletion Transformer

abstract: We propose the Insertion-Deletion Transformer, a novel transformer-based neural architecture and training method for sequence generation. The model consists of two phases that are executed iteratively, 1) an insertion phase and 2) a deletion phase. The insertion phase parameterizes a distribution of insertions on the current output hypothesis, while the deletion phase parameterizes a distribution of deletions over the current output hypothesis. The training method is a principled and simple algorithm, where the deletion model obtains its signal directly on-policy from the insertion model output. We demonstrate the effectiveness of our Insertion-Deletion Transformer on synthetic translation tasks, obtaining significant BLEU score improvement over an insertion-only model.


title: Ubiquitous impact of localised impurity states on the exchange coupling mechanism in magnetic topological insulators

abstract: Since the discovery of the quantum anomalous Hall effect in the magnetically doped topological insulators (MTI) Cr:(Bi,Sb)$_2$Te$_3$ and V:(Bi,Sb)$_2$Te$_3$, the search for the exchange coupling mechanisms underlying the onset of ferromagnetism has been a central issue, and a variety of different scenarios have been put forward. By combining resonant photoemission, X-ray magnetic dichroism and multiplet ligand field theory, we determine the local electronic and magnetic configurations of V and Cr impurities in (Bi,Sb)$_2$Te$_3$. While strong pd hybridisation is found for both dopant types, their 3d densities of states show pronounced differences. State-of-the-art first-principles calculations show how these impurity states mediate characteristic short-range pd exchange interactions, whose strength sensitively varies with the position of the 3d states relative to the Fermi level. Measurements on films with varying host stoichiometry support this trend. Our results establish the essential role of impurity-state mediated exchange interactions in the magnetism of MTI.


title: Physics-Constrained Bayesian Neural Network for Fluid Flow Reconstruction with Sparse and Noisy Data

abstract: In many applications, flow measurements are usually sparse and possibly noisy. The reconstruction of a high-resolution flow field from limited and imperfect flow information is significant yet challenging. In this work, we propose an innovative physics-constrained Bayesian deep learning approach to reconstruct flow fields from sparse, noisy velocity data, where equation-based constraints are imposed through the likelihood function and uncertainty of the reconstructed flow can be estimated. Specifically, a Bayesian deep neural network is trained on sparse measurement data to capture the flow field. In the meantime, the violation of physical laws will be penalized on a large number of spatiotemporal points where measurements are not available. A non-parametric variational inference approach is applied to enable efficient physics-constrained Bayesian learning. Several test cases on idealized vascular flows with synthetic measurement data are studied to demonstrate the merit of the proposed method.


title: A parabolic local problem with exponential decay of the resonance error for numerical homogenization

abstract: This paper aims at an accurate and efficient computation of effective quantities, e.g., the homogenized coefficients for approximating the solutions to partial differential equations with oscillatory coefficients. Typical multiscale methods are based on a micro-macro coupling, where the macro model describes the coarse scale behaviour, and the micro model is solved only locally to upscale the effective quantities, which are missing in the macro model. The fact that the micro problems are solved over small domains within the entire macroscopic domain, implies imposing artificial boundary conditions on the boundary of the microscopic domains. A naive treatment of these artificial boundary conditions leads to a first order error in $\varepsilon/\delta$, where $\varepsilon < \delta$ represents the characteristic length of the small scale oscillations and $\delta^d$ is the size of micro domain. This error dominates all other errors originating from the discretization of the macro and the micro problems, and its reduction is a main issue in today's engineering multiscale computations. The objective of the present work is to analyze a parabolic approach, first announced in [A. Abdulle, D. Arjmand, E. Paganoni, C. R. Acad. Sci. Paris, Ser. I, 2019], for computing the homogenized coefficients with arbitrarily high convergence rates in $\varepsilon/\delta$. The analysis covers the setting of periodic micro structure, and numerical simulations are provided to verify the theoretical findings for more general settings, e.g. random stationary micro structures.


title: Signal to Noise Ratio optimization for extended sources with a new kind of MURA masks

abstract: Using coded aperture, for localization of radioactive hot-spots, results in enhanced efficiency and under certain configurations wide Field of View (FOV). We present a coded aperture assembly technique which can be restructured easily, as well as the reduction of the intrinsic noise of coded apertures constructed with this technique, when they localize spatially extended {\gamma}-emitters. Specifically, Modified-Uniformly-Redundant-Array (MURA) coded apertures are structured by embedding lead spheres in a matrix of positions machined inside a transparent medium such as acrylic glass, resulting in an advantageous transparent to opaque area ratio and thus an improved detection efficiency. This configuration also induces a systematic, element-wise, noise on the Point-Spread-Function (PSF) of the correlation matrix. When imaging with these apertures extended hot-spots, a penumbra phenomenon occurrs and reduces this intrinsic noise, in the way a kernel filter does. Fast-Fourier-Transform (FFT) is used to analyze the effect of this phenomenon on the correlation matrix and to explain the maximization of its Signal-to-Noise Ratio (SNR) for certain extent of the hot-spots. Simulations have been used for the detailed study of the SNR dependence on the dimensions of the hot-spot, while experiments with two 99mTc cylindrical sources with 11mm and 24mm diameter, respectively and 1.5 MBq activity each, confirm the reduction of the intrinsic noise. The results define the way of optimization of the imaging setup for the detection of extended hot-spots. Such an optimization could be useful for example in the case of lymph nodes or thyroid remnant imaging in nuclear medicine. Finally, we propose a kernel filter, derived by the Auto-Correlation-Function (ACF), to be applied on PSFs with high intrinsic noise, in order to eliminate it.


title: A "Network Pruning Network" Approach to Deep Model Compression

abstract: We present a filter pruning approach for deep model compression, using a multitask network. Our approach is based on learning a a pruner network to prune a pre-trained target network. The pruner is essentially a multitask deep neural network with binary outputs that help identify the filters from each layer of the original network that do not have any significant contribution to the model and can therefore be pruned. The pruner network has the same architecture as the original network except that it has a multitask/multi-output last layer containing binary-valued outputs (one per filter), which indicate which filters have to be pruned. The pruner's goal is to minimize the number of filters from the original network by assigning zero weights to the corresponding output feature-maps. In contrast to most of the existing methods, instead of relying on iterative pruning, our approach can prune the network (original network) in one go and, moreover, does not require specifying the degree of pruning for each layer (and can learn it instead). The compressed model produced by our approach is generic and does not need any special hardware/software support. Moreover, augmenting with other methods such as knowledge distillation, quantization, and connection pruning can increase the degree of compression for the proposed approach. We show the efficacy of our proposed approach for classification and object detection tasks.


title: A new recursion for Bressoud's polynomials

abstract: A new recursion in only one variable allows very simple verifications of Bressoud's polynomial identities, which lead to the Rogers-Ramanujan identities. This approach might be compared with an earlier approach due to Chapman. Applying the $q$-Chu-Vandermonde convolution, as suggested by Cigler, makes the computations particularly simple and elementary. The same treatment is also applied to the Santos polynomials and perhaps more polynomials from a list of Rogers-Ramanujan like polynomials \cite{Sills}.


title: Nonassociativity of the Norton Algebras of some distance regular graphs

abstract: A Norton algebra is an eigenspace of a distance regular graph endowed with a commutative nonassociative product called the Norton product, which is defined as the projection of the entrywise product onto this eigenspace. The Norton algebras are useful in finite group theory as they have interesting automorphism groups. We provide a precise quantitative measurement for the nonassociativity of the Norton product on the eigenspace of the second largest eigenvalue of the Johnson graphs, Grassman graphs, Hamming graphs, and dual polar graphs, based on the formulas for this product established in previous work of Levstein, Maldonado and Penazzi. Our result shows that this product is as nonassociative as possible except for two cases, one being the trivial vanishing case while the other having connections with the integer sequence A000975 on OEIS and the so-called double minus operation studied recently by Huang, Mickey, and Xu.


title: Segmentation with Residual Attention U-Net and an Edge-Enhancement Approach Preserves Cell Shape Features

abstract: The ability to extrapolate gene expression dynamics in living single cells requires robust cell segmentation, and one of the challenges is the amorphous or irregularly shaped cell boundaries. To address this issue, we modified the U-Net architecture to segment cells in fluorescence widefield microscopy images and quantitatively evaluated its performance. We also proposed a novel loss function approach that emphasizes the segmentation accuracy on cell boundaries and encourages shape feature preservation. With a 97% sensitivity, 93% specificity, 91% Jaccard similarity, and 95% Dice coefficient, our proposed method called Residual Attention U-Net with edge-enhancement surpassed the state-of-the-art U-Net in segmentation performance as evaluated by the traditional metrics. More remarkably, the same proposed candidate also performed the best in terms of the preservation of valuable shape features, namely area, eccentricity, major axis length, solidity and orientation. These improvements on shape feature preservation can serve as useful assets for downstream cell tracking and quantification of changes in cell statistics or features over time.


title: Supervised Segmentation of Retinal Vessel Structures Using ANN

abstract: In this study, a supervised retina blood vessel segmentation process was performed on the green channel of the RGB image using artificial neural network (ANN). The green channel is preferred because the retinal vessel structures can be distinguished most clearly from the green channel of the RGB image. The study was performed using 20 images in the DRIVE data set which is one of the most common retina data sets known. The images went through some preprocessing stages like contrastlimited adaptive histogram equalization (CLAHE), color intensity adjustment, morphological operations and median and Gaussian filtering to obtain a good segmentation. Retinal vessel structures were highlighted with top-hat and bot-hat morphological operations and converted to binary image by using global thresholding. Then, the network was trained by the binary version of the images specified as training images in the dataset and the targets are the images segmented manually by a specialist. The average segmentation accuracy for 20 images was found as 0.9492.


title: Asymptotic expansion of the variance of random zeros on complex manifolds

abstract: Linear statistics of random zero sets are integrals of smooth differential forms over the zero set and as such are smooth analogues of the volume of the random zero set inside a fixed domain. We derive an asymptotic expansion for the variance of linear statistics of the zero divisors of random holomorphic sections of powers of a positive line bundle on a compact K\"ahler manifold. This expansion sharpens the leading-order asymptotics (in the codimension one case) given by Shiffman--Zelditch in 2010.


title: Substituting Gadolinium in Brain MRI Using DeepContrast

abstract: Cerebral blood volume (CBV) is a hemodynamic correlate of oxygen metabolism and reflects brain activity and function. High-resolution CBV maps can be generated using the steady-state gadolinium-enhanced MRI technique. Such a technique requires an intravenous injection of exogenous gadolinium based contrast agent (GBCA) and recent studies suggest that the GBCA can accumulate in the brain after frequent use. We hypothesize that endogenous sources of contrast might exist within the most conventional and commonly acquired structural MRI, potentially obviating the need for exogenous contrast. Here, we test this hypothesis by developing and optimizing a deep learning algorithm, which we call DeepContrast, in mice. We find that DeepContrast performs equally well as exogenous GBCA in mapping CBV of the normal brain tissue and enhancing glioblastoma. Together, these studies validate our hypothesis that a deep learning approach can potentially replace the need for GBCAs in brain MRI.


title: Resolution of spin Hall and anisotropic magnetoresistance in Pt/EuO$_{1-x}$

abstract: We report on the angular and field dependence of the magnetoresistance (MR) in bilayers of Pt/EuO_{1-x} thin films, measured in both in-plane and out-of-plane geometries at different temperatures (T). Presence of oxygen vacancies manifested by a metal-insulator transition as well as a high-T ferromagnet to paramagnet transition (T_P) were observed in the bilayers. The Anisotropic Magnetoresistance (AMR) could be extracted in the entire T-range, even above T_P, exhibiting two sign crossovers. We attribute its T-evolution to the rotation of easy axis direction from a high-T out-of-plane to a low-T in-plane orientation. In addition, considering MR contributions from the films' (111) texture and interface, we identify a T-window wherein the spin Hall effect induced spin Hall magnetoresistance (SMR) could be extracted.


title: Exponential quantum communication reductions from generalizations of the Boolean Hidden Matching problem

abstract: In this work we revisit the Boolean Hidden Matching communication problem, which was the first communication problem in the one-way model to demonstrate an exponential classical-quantum communication separation. In this problem, Alice's bits are matched into pairs according to a partition that Bob holds. These pairs are compressed using a Parity function and it is promised that the final bit-string is equal either to another bit-string Bob holds, or its complement. The problem is to decide which case is the correct one. Here we generalize the Boolean Hidden Matching problem by replacing the parity function with an arbitrary Boolean function $f$. Efficient communication protocols are presented depending on the sign-degree of $f$. If its sign-degree is less than or equal to 1, we show an efficient classical protocol. If its sign-degree is less than or equal to $2$, we show an efficient quantum protocol. We then completely characterize the classical hardness of all symmetric functions $f$ of sign-degree greater than or equal to $2$, except for one family of specific cases. We also prove, via Fourier analysis, a classical lower bound for any function $f$ whose pure high degree is greater than or equal to $2$. Similarly, we prove, also via Fourier analysis, a quantum lower bound for any function $f$ whose pure high degree is greater than or equal to $3$. These results give a large family of new exponential classical-quantum communication separations.


title: When is $\overline{M}_{0,n}(\mathbb{P}^1,1)$ a Mori dream space?

abstract: We prove that the moduli space of $n$-pointed stable maps $\overline{M}_{0,n}(\mathbb{P}^1,1)$ is a Mori dream space whenever the moduli space $\overline{M}_{0,n+3}$ of $(n+3)$-pointed rational curves is. We also show that $\overline{M}_{0,n}(\mathbb{P}^1,1)$ is a log Fano variety for $n \le 5$.


title: Probing the Influence of a Tachocline in Simulated M-Dwarf Dynamos

abstract: In recent years, observers have found that the fraction of M-stars demonstrating significant magnetic activity transitions sharply from roughly $10\%$ for main-sequence stars earlier (more massive) than spectral type M3.5 (0.35 M$_\odot$) to nearly $90\%$ for stars later than M3.5. The latter are typically rotating quite rapidly, suggesting differing spin-down histories. Tantalizingly, it is also later than M3.5 at which main-sequence stars become fully convective, and may no longer contain a tachocline, a layer of rotational shear revealed by helioseismology to separate the radiative zone (RZ) and convection zone (CZ). We turn here to the more massive M-stars to study the impact such a layer may have on their internal dynamics. Using the spherical 3D MHD simulation code Rayleigh, we compare the properties of convective dynamos generated within quickly rotating (1, 2, and 4 $\Omega_\odot$) M2 (0.4 M$_\odot$) stars, with the computational domain either terminating at the base of the convection zone or permitting overshoot into the underlying stable region. We find that a tachocline is not necessary for the organization of strong toroidal wreaths of magnetism in these stars, though its presence can increase the coupling of mean field amplitudes to the stellar rotation rate. Additionally, in stars that undergo periodic cycles, we find that the presence of a tachocline tends to make the cycles both longer and more regular than they would have otherwise been. Finally, we find that the tachocline helps to enhance the surface poloidal fields and organize them into larger spatial scales, both of which provide favorable conditions for more rapid angular momentum loss through a magnetized stellar wind.


title: Thermophoresis in Liquids and its Connection to Equilibrium Quantities

abstract: Thermophoresis is the process of particles moving along a temperature gradient in thermodynamic systems. Even though it has been studied for over 150 years, there is neither a complete theoretical description of thermophoresis in liquids nor of the Soret coefficient, the quantifiying measure of the Soret equilibrium. Recent studies connect its nature to equilibrium properties of the system, namely the excess enthalpy and the excess entropy, while there is still a debate over which of both describes the Soret coefficient more accurately and whether it can even be represented using those quantities. In this work I present a theoretical description for both cases based on density analysis by means of Brownian motion and dynamical density functional theory, where I assume local equilibrium. The interpretation of the Brownian stochastic differential equation (SDE) as an Ito or Stratonovich SDE is shown to have an influence on the outcome of the density. I argue that a Soret coefficient proportional to the excess enthalpy is connected to a system in a thermal gradient following its equation of state. Furthermore, I derive that the Boltzmann distribution law for external potentials does not hold for systems in thermal gradients but has to be replaced by a more general quantity. The theoretical predictions are consequently tested by means of BD simulations, where several systems are shown to follow their equation of state. The Soret coefficient of a solute in a toy model solvent is shown to be proportional to the solute's equilibrium solvation enthalpy. I further attempt to verify the theoretic derivations by means of MD simulations of noble gas solutes in SPC/E water. While the Soret equilibrium SPC/E water density follows its equation of state, neither Soret coefficient nor Soret equilibrium solute densitie entail a coherent agreement with any of the theoretical predictions.


title: An infinite product on the Teichm\"{u}ller space of the once-punctured torus

abstract: We prove the identity $$ \prod_{\gamma}\left(\frac{e^{l(\gamma)}+1}{e^{l(\gamma)}-1}\right)^{2h}=\exp\left(\frac{l_1+l_2+l_3}{2}\right), $$ (or $$ \prod_{\gamma}\left(\frac{t(\gamma)^2}{t(\gamma)^2-4}\right)^h=\frac{t_1+\sqrt{t_1^2-4}}{2}\cdot\frac{t_2+\sqrt{t_2^2-4}}{2}\cdot\frac{t_3+\sqrt{t_3^2-4}}{2} $$ in trace coordinates), where the product is over all simple closed geodesics on the once-punctured torus, $l(\gamma)=2\operatorname{arccosh}(t(\gamma)/2)$ is the length of the geodesic, and $l_i$ ($t_i$) are the lengths (traces) of any triple of simple geodesics $\{\gamma_i\}$ intersecting at a single point. The exponent $h=h(\gamma;\{\gamma_i\})$ is a positive integer "height" which increases as we move away from the chosen triple $\{\gamma_i\}$ in its orbit under $SL_2(\mathbb{Z})$ (see Figure 1 for the "definition by picture"). For comparison, a short proof of McShane's identity $$ \sum_{\gamma}\frac{1}{1+e^{l(\gamma)}}=\frac{1}{2}=\sum_{\gamma}\frac{1-\sqrt{1-4/t(\gamma)^2}}{2} $$ in the same spirit is given in an appendix. Both proofs are elementary and proceed by "integrating" around the chosen triple $\{\gamma_i\}$ in its Teichm\"{u}ller orbit.


title: Efficient Discrete-Event Based Particle Tracking Simulation for High Energy Physics

abstract: This work presents novel discrete event-based simulation algorithms based on the Quantized State System (QSS) numerical methods. QSS provides attractive features for particle transportation processes, in particular a very efficient handling of discontinuities in the simulation of continuous systems. We focus on High Energy Physics (HEP) particle tracking applications that typically rely on discrete time-based methods, and study the advantages of adopting a discrete event-based numerical approach that resolves efficiently the crossing of geometry boundaries by a traveling particle. For this purpose we follow two complementary strategies. First, a new co-simulation technique connects the Geant4 simulation toolkit with a standalone QSS solver. Second, a new native QSS numerical stepper is embedded into Geant4. We compare both approaches against the latest Geant4 default steppers in different HEP setups, including a complex real scenario (the CMS particle detector at CERN). Our techniques achieve relevant simulation speedups in a wide range of scenarios, particularly when the intensity of discrete-event handling dominates performance in the solving of the continuous laws of particle motion.


title: Mode-Assisted Unsupervised Learning of Restricted Boltzmann Machines

abstract: Restricted Boltzmann machines (RBMs) are a powerful class of generative models, but their training requires computing a gradient that, unlike supervised backpropagation on typical loss functions, is notoriously difficult even to approximate. Here, we show that properly combining standard gradient updates with an off-gradient direction, constructed from samples of the RBM ground state (mode), improves their training dramatically over traditional gradient methods. This approach, which we call mode training, promotes faster training and stability, in addition to lower converged relative entropy (KL divergence). Along with the proofs of stability and convergence of this method, we also demonstrate its efficacy on synthetic datasets where we can compute KL divergences exactly, as well as on a larger machine learning standard, MNIST. The mode training we suggest is quite versatile, as it can be applied in conjunction with any given gradient method, and is easily extended to more general energy-based neural network structures such as deep, convolutional and unrestricted Boltzmann machines.


title: An Attack-Resilient Pulse-Based Synchronization Strategy for General Connected Topologies

abstract: Synchronization of pulse-coupled oscillators (PCOs) has gained significant attention recently due to increased applications in sensor networks and wireless communications. Given the distributed and unattended nature of wireless sensor networks, it is imperative to enhance the resilience of pulse-based synchronization against malicious attacks. However, most existing results on resilient PCO synchronization are obtained for all-to-all networks. We propose a new pulse-based synchronization mechanism to improve the resilience of PCO synchronization that is applicable under general connected topologies. Under the proposed synchronization mechanism, we rigorously characterize the condition for stealthy Byzantine attacks and prove that perfect synchronization of legitimate oscillators can be guaranteed in the presence of multiple stealthy Byzantine attackers, irrespective of whether the attackers collude with each other or not. The new mechanism can guarantee resilient synchronization even when the initial phases of legitimate oscillators are widely distributed in a half circle, which is in distinct difference from most existing attack-resilient synchronization algorithms (including the seminal paper from Lamport and Melliar-Smith [1]) that require a priori (almost) synchronization among legitimate oscillators. Numerical simulation results are given to confirm the theoretical results.


title: Partial Coherence in modern optics: Emil Wolf's legacy in the 21st century

abstract: We highlight the impact of Emil Wolf's work on coherence and polarization on an ever increasing amount of applications in the 21 st century. We present a brief review of how partial coherence at the level of increasing order of coherence functions is leading to evolution in the better methods for microscopy, imaging, optical coherence tomography; speckle imaging; propagation through random media. This evolution in our capabilities is expected to have wide ramifications in Science and Engineering.


title: 1/8-BPS Couplings and Exceptional Automorphic Functions

abstract: Unlike the $\mathcal{R}^4$ and $\nabla^4\mathcal{R}^4$ couplings, whose coefficients are Langlands-Eisenstein series of the U-duality group, the coefficient $\mathcal{E}_{(0,1)}^{(d)}$ of the $\nabla^6\mathcal{R}^4$ interaction in the low-energy effective action of type II strings compactified on a torus $T^d$ belongs to a more general class of automorphic functions, which satisfy Poisson rather than Laplace-type equations. In earlier work, it was proposed that the exact coefficient is given by a two-loop integral in exceptional field theory, with the full spectrum of mutually 1/2-BPS states running in the loops, up to the addition of a particular Langlands-Eisenstein series. Here we compute the weak coupling and large radius expansions of these automorphic functions for any $d$. We find perfect agreement with perturbative string theory up to genus three, along with non-perturbative corrections which have the expected form for 1/8-BPS instantons and bound states of 1/2-BPS instantons and anti-instantons. The additional Langlands-Eisenstein series arises from a subtle cancellation between the two-loop amplitude with 1/4-BPS states running in the loops, and the three-loop amplitude with mutually 1/2-BPS states in the loops. For $d=4$, the result is shown to coincide with an alternative proposal in terms of a covariantised genus-two string amplitude, due to interesting identities between the Kawazumi-Zhang invariant of genus-two curves and its tropical limit, and between double lattice sums for the particle and string multiplets, which may be of independent mathematical interest.


title: The equivariant parametrized $h$-cobordism theorem, the non-manifold part

abstract: We construct a map from the suspension $G$-spectrum $\Sigma_G^\infty M$ of a smooth compact $G$-manifold to the equivariant $A$-theory spectrum $A_G(M)$, and we show that its fiber is, on fixed points, a wedge of stable $h$-cobordism spectra. This map is constructed as a map of spectral Mackey functors, which is compatible with tom Dieck style splitting formulas on fixed points. In order to synthesize different definitions of the suspension $G$-spectrum as a spectral Mackey functor, we present a new perspective on spectral Mackey functors, viewing them as multifunctors on indexing categories for "rings on many objects" and modules over such. This perspective should be of independent interest.


title: Simplification of Indoor Space Footprints

abstract: Simplification is one of the fundamental operations used in geoinformation science (GIS) to reduce size or representation complexity of geometric objects. Although different simplification methods can be applied depending on one's purpose, a simplification that many applications employ is designed to preserve their spatial properties after simplification. This article addresses one of the 2D simplification methods, especially working well on human-made structures such as 2D footprints of buildings and indoor spaces. The method simplifies polygons in an iterative manner. The simplification is segment-wise and takes account of intrusion, extrusion, offset, and corner portions of 2D structures preserving its dominant frame.


title: Fractional Orlicz-Sobolev embeddings

abstract: The optimal Orlicz target space is exhibited for embeddings of fractional-order Orlicz-Sobolev spaces in $\mathbb R^n$. An improved embedding with an Orlicz-Lorentz target space, which is optimal in the broader class of all rearrangement-invariant spaces, is also established. Both spaces of order $s\in (0,1)$, and higher-order spaces are considered. Related Hardy type inequalities are proposed as well. An extension theorem is proved, that enables us to derive embeddings for spaces defined in Lipschitz domains. Necessary and sufficient conditions for the compactness of fractional Orlicz-Sobolev embeddings are provided.


title: Image Segmentation Using Deep Learning: A Survey

abstract: Image segmentation is a key topic in image processing and computer vision with applications such as scene understanding, medical image analysis, robotic perception, video surveillance, augmented reality, and image compression, among many others. Various algorithms for image segmentation have been developed in the literature. Recently, due to the success of deep learning models in a wide range of vision applications, there has been a substantial amount of works aimed at developing image segmentation approaches using deep learning models. In this survey, we provide a comprehensive review of the literature at the time of this writing, covering a broad spectrum of pioneering works for semantic and instance-level segmentation, including fully convolutional pixel-labeling networks, encoder-decoder architectures, multi-scale and pyramid based approaches, recurrent networks, visual attention models, and generative models in adversarial settings. We investigate the similarity, strengths and challenges of these deep learning models, examine the most widely used datasets, report performances, and discuss promising future research directions in this area.


title: Newtonian Monte Carlo: single-site MCMC meets second-order gradient methods

abstract: Single-site Markov Chain Monte Carlo (MCMC) is a variant of MCMC in which a single coordinate in the state space is modified in each step. Structured relational models are a good candidate for this style of inference. In the single-site context, second order methods become feasible because the typical cubic costs associated with these methods is now restricted to the dimension of each coordinate. Our work, which we call Newtonian Monte Carlo (NMC), is a method to improve MCMC convergence by analyzing the first and second order gradients of the target density to determine a suitable proposal density at each point. Existing first order gradient-based methods suffer from the problem of determining an appropriate step size. Too small a step size and it will take a large number of steps to converge, while a very large step size will cause it to overshoot the high density region. NMC is similar to the Newton-Raphson update in optimization where the second order gradient is used to automatically scale the step size in each dimension. However, our objective is to find a parameterized proposal density rather than the maxima. As a further improvement on existing first and second order methods, we show that random variables with constrained supports don't need to be transformed before taking a gradient step. We demonstrate the efficiency of NMC on a number of different domains. For statistical models where the prior is conjugate to the likelihood, our method recovers the posterior quite trivially in one step. However, we also show results on fairly large non-conjugate models, where NMC performs better than adaptive first order methods such as NUTS or other inexact scalable inference methods such as Stochastic Variational Inference or bootstrapping.


title: Notes on Communication and Computation in Secure Distributed Matrix Multiplication

abstract: We consider the problem of secure distributed matrix multiplication in which a user wishes to compute the product of two matrices with the assistance of honest but curious servers. In this paper, we answer the following question: Is it beneficial to offload the computations if security is a concern? We answer this question in the affirmative by showing that by adjusting the parameters in a polynomial code we can obtain a trade-off between the user's and the servers' computational time. Indeed, we show that if the computational time complexity of an operation in $\mathbb{F}_q$ is at most $\mathcal{Z}_q$ and the computational time complexity of multiplying two $n\times n$ matrices is $\mathcal{O}(n^\omega \mathcal{Z}_q)$ then, by optimizing the trade-off, the user together with the servers can compute the multiplication in $\mathcal{O}(n^{4-\frac{6}{\omega+1}} \mathcal{Z}_q)$ time. We also show that if the user is only concerned in optimizing the download rate, a common assumption in the literature, then the problem can be converted into a simple private information retrieval problem by means of a scheme we call Private Oracle Querying. However, this comes at large upload and computational costs for both the user and the servers.


title: Foundations of Quantum Mechanics according to my teachers

abstract: These are author's recollections of informal discussions on foundations of quantum mechanics, which happened in his presence in 1970-80s


title: A survey for variable young stars with small telescopes: II -- Mapping a protoplanetary disk with stable structures at 0.15 AU

abstract: The HOYS citizen science project conducts long term, multifilter, high cadence monitoring of large YSO samples with a wide variety of professional and amateur telescopes. We present the analysis of the light curve of V1490Cyg in the Pelican Nebula. We show that colour terms in the diverse photometric data can be calibrated out to achieve a median photometric accuracy of 0.02mag in broadband filters, allowing detailed investigations into a variety of variability amplitudes over timescales from hours to several years. Using GaiaDR2 we estimate the distance to the Pelican Nebula to be 870$^{+70}_{-55}$pc. V1490Cyg is a quasi-periodic dipper with a period of 31.447$\pm$0.011d. The obscuring dust has homogeneous properties, and grains larger than those typical in the ISM. Larger variability on short timescales is observed in U and R$_c-$H$\alpha$, with U-amplitudes reaching 3mag on timescales of hours, indicating the source is accreting. The H$\alpha$ equivalent width and NIR/MIR colours place V1490Cyg between CTTS/WTTS and transition disk objects. The material responsible for the dipping is located in a warped inner disk, about 0.15AU from the star. This mass reservoir can be filled and emptied on time scales shorter than the period at a rate of up to 10$^{-10}$M$_\odot$/yr, consistent with low levels of accretion in other T Tauri stars. Most likely the warp at this separation from the star is induced by a protoplanet in the inner accretion disk. However, we cannot fully rule out the possibility of an AA Tau-like warp, or occultations by the Hill sphere around a forming planet.


title: On Model Evaluation under Non-constant Class Imbalance

abstract: Many real-world classification problems are significantly class-imbalanced to detriment of the class of interest. The standard set of proper evaluation metrics is well-known but the usual assumption is that the test dataset imbalance equals the real-world imbalance. In practice, this assumption is often broken for various reasons. The reported results are then often too optimistic and may lead to wrong conclusions about industrial impact and suitability of proposed techniques. We introduce methods focusing on evaluation under non-constant class imbalance. We show that not only the absolute values of commonly used metrics, but even the order of classifiers in relation to the evaluation metric used is affected by the change of the imbalance rate. Finally, we demonstrate that using subsampling in order to get a test dataset with class imbalance equal to the one observed in the wild is not necessary, and eventually can lead to significant errors in classifier's performance estimate.


title: A C Code Generator for Fast Inference and Simple Deployment of Convolutional Neural Networks on Resource Constrained Systems

abstract: Inference of Convolutional Neural Networks in time critical applications usually requires a GPU. In robotics or embedded devices these are often not available due to energy, space and cost constraints. Furthermore, installation of a deep learning framework or even a native compiler on the target platform is not possible. This paper presents a neural network code generator (NNCG) that generates from a trained CNN a plain ANSI C code file that encapsulates the inference in single a function. It can easily be included in existing projects and due to lack of dependencies, cross compilation is usually possible. Additionally, the code generation is optimized based on the known trained CNN and target platform following four design principles. The system is evaluated utilizing small CNN designed for this application. Compared to TensorFlow XLA and Glow speed-ups of up to 11.81 can be shown and even GPUs are outperformed regarding latency.


title: Consumer-Driven Explanations for Machine Learning Decisions: An Empirical Study of Robustness

abstract: Many proposed methods for explaining machine learning predictions are in fact challenging to understand for nontechnical consumers. This paper builds upon an alternative consumer-driven approach called TED that asks for explanations to be provided in training data, along with target labels. Using semi-synthetic data from credit approval and employee retention applications, experiments are conducted to investigate some practical considerations with TED, including its performance with different classification algorithms, varying numbers of explanations, and variability in explanations. A new algorithm is proposed to handle the case where some training examples do not have explanations. Our results show that TED is robust to increasing numbers of explanations, noisy explanations, and large fractions of missing explanations, thus making advances toward its practical deployment.


title: Advbox: a toolbox to generate adversarial examples that fool neural networks

abstract: In recent years, neural networks have been extensively deployed for computer vision tasks, particularly visual classification problems, where new algorithms reported to achieve or even surpass the human performance. Recent studies have shown that they are all vulnerable to the attack of adversarial examples. Small and often imperceptible perturbations to the input images are sufficient to fool the most powerful neural networks. \emph{Advbox} is a toolbox to generate adversarial examples that fool neural networks in PaddlePaddle, PyTorch, Caffe2, MxNet, Keras, TensorFlow and it can benchmark the robustness of machine learning models. Compared to previous work, our platform supports black box attacks on Machine-Learning-as-a-service, as well as more attack scenarios, such as Face Recognition Attack, Stealth T-shirt, and DeepFake Face Detect. The code is licensed under the Apache 2.0 and is openly available at https://github.com/advboxes/AdvBox. Advbox now supports Python 3.


title: Ownership Structure Variation and Firm Efficiency

abstract: Firms with different ownership structures could be argued to have different levels of efficiency.Highly concentrated firms are expected to be more efficient as this type of ownership structure may alleviate the conflict of interest between managers and shareholders.In Malaysia, public-listed firms have been found to have highly concentrated ownership structure.However, whether this evidence holds for every industry has not been established.Hence, the objective of this paper is to investigate whether there are variations in ownership structure and firm's efficiency across sectors.To achieve this objective, the frequency distributions of ownership structure were calculated and firms efficiency scores for consumer products, industrial products, construction and trading/services sectors were measured.Data Envelopment Analysis(DEA) under the assumptions of constant returns to scale(CRS) and variable returns to scale(VRS) was employed to estimate firms efficiency scores.A sample of 156 firms listed on the Kuala Lumpur Stock Exchange(KLSE) was selected using the stratified random sampling method. The findings have shown that there are variations in firm ownership structure and efficiency across sectors.


title: Novel Approaches to High-Power Proton Beams

abstract: The potentially realizable beam power at the Fermilab long-baseline neutrino program has motivated a reinvigorated design and optimization effort for a rapid-cycling synchrotron (RCS) intensity upgrade of the Fermilab proton complex. We examine areas of technological development with the potential for high-impact on the Fermilab RCS design - low-loss slip-stacking, advanced neutrino-target R\&D, laser-stripping H$^{-}$ injection, fast-ramping super-ferric magnets, nonlinear integrable optics, electron lens devices, and next-generation halo monitors. A brief overview is given on the Fermilab Accelerator Science \& Technology (FAST) facility, where the latter three technologies will undergo comprehensive beam tests.


title: Stable homotopy hypothesis in the Tamsamani model

abstract: We prove that symmetric monoidal weak n-groupoids in the Tamsamani model provide a model for stable n-types. Moreover, we recover the classical statement that Picard categories model stable 1-types.


title: VSEC-LDA: Boosting Topic Modeling with Embedded Vocabulary Selection

abstract: Topic modeling has found wide application in many problems where latent structures of the data are crucial for typical inference tasks. When applying a topic model, a relatively standard pre-processing step is to first build a vocabulary of frequent words. Such a general pre-processing step is often independent of the topic modeling stage, and thus there is no guarantee that the pre-generated vocabulary can support the inference of some optimal (or even meaningful) topic models appropriate for a given task, especially for computer vision applications involving "visual words". In this paper, we propose a new approach to topic modeling, termed Vocabulary-Selection-Embedded Correspondence-LDA (VSEC-LDA), which learns the latent model while simultaneously selecting most relevant words. The selection of words is driven by an entropy-based metric that measures the relative contribution of the words to the underlying model, and is done dynamically while the model is learned. We present three variants of VSEC-LDA and evaluate the proposed approach with experiments on both synthetic and real databases from different applications. The results demonstrate the effectiveness of built-in vocabulary selection and its importance in improving the performance of topic modeling.


title: Three-level rate equations in cold, disordered Rydberg gases

abstract: We have investigated formation of structures of Rydberg atoms excited from a disordered gas of ultracold atoms, using rate equations for two-photon Rydberg excitation in a single atom without eliminating the intermediate state. We have explored the validity range of these rate equations and defined a simple measure to determine, whether our model is applicable for a given set of laser parameters. We have applied these rate equations in Monte Carlo simulations of ultracold gases, for different laser beam profiles, and compared these simulations to experimental observations and find a general agreement.


title: Power Cycling Test Bench for Accelerated Life Testing for Reliability Assessment of SiC-MOSFET in Extreme Offshore Environment

abstract: The reliability of power semiconductor switches is important when considering their vital role in power electronic converters for downhole subsea applications. Respect to technology advancements in material sciences, power MOSFETs with wide band gap materials have been proposed such as silicon carbide (SiC) and gallium nitride (GaN) as an alternative to existing silicon (Si) based MOSFETs and IGBTs. However, reliability analysis should be performed before substituting SiC-MOSFETs in the place of existing Si-MOSFETs and IGBTs. Due to costly equipment of experimental test setup for accelerated life test, a good reliable and precise simulation-based test bench should be used to test the life test procedure before implementing actual hardware. Therefore, this paper introduces a power cycle (PC) test bench for accelerated life testing for reliability assessment of SiC-MOSFET in harsh offshore environment. The introduced test bench is a simulation-based of power switch in SimScape and LTspice and has been validated with datasheet of 1.2 kV SiC-MOSFET, CAS300M12BM2 by CREE. Preliminary hardware circuits are also shown for further experimental tests. The captured data from the Device-Under-Test (DUT) in different ambient temperatures are envisioned and provide critical information about the failure mechanisms and lifetime characteristics of power devices. The provided lifetime characteristics data of SiC-MOSFET can be used to statistically estimate the Remaining-Useful-Lifetime (RUL) of component in a real application such as downhole motor drives.


title: Complete and Sufficient Spatial Domination of Multidimensional Rectangles

abstract: Rectangles are used to approximate objects, or sets of objects, in a plethora of applications, systems and index structures. Many tasks, such as nearest neighbor search and similarity ranking, require to decide if objects in one rectangle A may, must, or must not be closer to objects in a second rectangle B, than objects in a third rectangle R. To decide this relation of "Spatial Domination" it can be shown that using minimum and maximum distances it is often impossible to detect spatial domination. This spatial gem provides a necessary and sufficient decision criterion for spatial domination that can be computed efficiently even in higher dimensional space. In addition, this spatial gem provides an example, pseudocode and an implementation in Python.


title: The Error Probability of Maximum-Likelihood Decoding over Two Deletion Channels

abstract: This paper studies the problem of reconstructing a word given several of its noisy copies. This setup is motivated by several applications, among them is reconstructing strands in DNA-based storage systems. Under this paradigm, a word is transmitted over some fixed number of identical independent channels and the goal of the decoder is to output the transmitted word or some close approximation. The main focus of this paper is the case of two deletion channels and studying the error probability of the maximum-likelihood (ML) decoder under this setup. First, it is discussed how the ML decoder operates. Then, we observe that the dominant error patterns are deletions in the same run or errors resulting from alternating sequences. Based on these observations, it is derived that the error probability of the ML decoder is roughly $\frac{3q-1}{q-1}p^2$, when the transmitted word is any $q$-ary sequence and $p$ is the channel's deletion probability. We also study the cases when the transmitted word belongs to the Varshamov Tenengolts (VT) code or the shifted VT code. Lastly, the insertion channel is studied as well. These theoretical results are verified by corresponding simulations.


title: Compressing Permutation Groups into Grammars and Polytopes. A Graph Embedding Approach

abstract: It can be shown that each permutation group $G \sqsubseteq S_n$ can be embedded, in a well defined sense, in a connected graph with $O(n+|G|)$ vertices. Some groups, however, require much fewer vertices. For instance, $S_n$ itself can be embedded in the $n$-clique $K_n$, a connected graph with n vertices. In this work, we show that the minimum size of a context-free grammar generating a finite permutation group $G \sqsubseteq S_n$ can be upper bounded by three structural parameters of connected graphs embedding $G$: the number of vertices, the treewidth, and the maximum degree. More precisely, we show that any permutation group $G \sqsubseteq S_n$ that can be embedded into a connected graph with $m$ vertices, treewidth k, and maximum degree $\Delta$, can also be generated by a context-free grammar of size $2^{O(k\Delta\log\Delta)}\cdot m^{O(k)}$. By combining our upper bound with a connection between the extension complexity of a permutation group and the grammar complexity of a formal language, we also get that these permutation groups can be represented by polytopes of extension complexity $2^{O(k \Delta\log \Delta)}\cdot m^{O(k)}$. The above upper bounds can be used to provide trade-offs between the index of permutation groups, and the number of vertices, treewidth and maximum degree of connected graphs embedding these groups. In particular, by combining our main result with a celebrated $2^{\Omega(n)}$ lower bound on the grammar complexity of the symmetric group $S_n$ we have that connected graphs of treewidth $o(n/\log n)$ and maximum degree $o(n/\log n)$ embedding subgroups of $S_n$ of index $2^{cn}$ for some small constant $c$ must have $n^{\omega(1)}$ vertices. This lower bound can be improved to exponential on graphs of treewidth $n^{\varepsilon}$ for $\varepsilon<1$ and maximum degree $o(n/\log n)$.


title: A higher rank rigidity theorem for convex real projective manifolds

abstract: For convex real projective manifolds we prove an analogue of the higher rank rigidity theorem of Ballmann and Burns-Spatzier.


title: GPU Tensor Cores for fast Arithmetic Reductions

abstract: This work proposes a GPU tensor core approach that encodes the arithmetic reduction of $n$ numbers as a set of chained $m \times m$ matrix multiply accumulate (MMA) operations executed in parallel by GPU tensor cores. The asymptotic running time of the proposed chained tensor core approach is $T(n)=5 log_{m^2}{n}$ and its speedup is $S=\dfrac{4}{5} log_{2}{m^2}$ over the classic $O(n \log n)$ parallel reduction algorithm. Experimental performance results show that the proposed reduction method is $\sim 3.2 \times$ faster than a conventional GPU reduction implementation, and preserves the numerical precision because the sub-results of each chain of $R$ MMAs is kept as a 32-bit floating point value, before being all reduced into as a final 32-bit result. The chained MMA design allows a flexible configuration of thread-blocks; small thread-blocks of 32 or 128 threads can still achieve maximum performance using a chain of $R=4,5$ MMAs per block, while large thread-blocks work best with $R=1$. The results obtained in this work show that tensor cores can indeed provide a significant performance improvement to non-Machine Learning applications such as the arithmetic reduction, which is an integration tool for studying many scientific phenomena.


title: Queue-Aware Beam Scheduling for Half-Duplex mmWave Relay Networks

abstract: Millimeter wave (mmWave) bands are considered a powerful key enabler for next generation (5G) mobile networks by providing multi-Gbps data rates. However, their severe pathloss and sensitivity to blockage present challenges in practical implementation. One effective way to mitigate these effects and to increase communication range is beamforming in combination with relaying. In this paper, we focus on two typical mmWave relay networks and for each network, we propose three beam scheduling methods to approach the network information theoretic capacity. The proposed beam scheduling methods include the deterministic horizontal continuous edge coloring (HC-EC) scheduler, the adaptive back pressure (BP) scheduler and the adaptive low-delay new back pressure (newBP) scheduler. With the aid of computer simulations, we show that within the network capacity range, the proposed schedulers provide good guarantees for the network stability, meanwhile achieve very low packet end-to-end delay.


title: Expected performances of the Characterising Exoplanet Satellite (CHEOPS) II. The CHEOPS simulator

abstract: The CHaracterising ExOPlanet Satellite (CHEOPS) is a mission dedicated to the search for exoplanetary transits through high precision photometry of bright stars already known to host planets. The telescope will provide the unique capability of determining accurate radii for planets whose masses have already been measured from ground-based spectroscopic surveys. This will allow a first-order characterisation of the planets' internal structure through the determination of the bulk density, providing direct insight into their composition. The CHEOPS simulator has been developed to perform detailed simulations of the data which is to be received from the CHEOPS satellite. It generates accurately simulated images that can be used to explore design options and to test the on-ground data processing, in particular, the pipeline producing the photometric time series. It is, thus, a critical tool for estimating the photometric performance expected in flight and to guide photometric analysis. It can be used to prepare observations, consolidate the noise budget, and asses the performance of CHEOPS in realistic astrophysical fields that are difficult to reproduce in the laboratory. Images generated by CHEOPSim take account of many detailed effects, including variations of the incident signal flux and backgrounds, and detailed modelling of the satellite orbit, pointing jitter and telescope optics, as well as the CCD response, noise and readout. The simulator results presented in this paper have been used in the context of validating the data reduction processing chain, in which image time series generated by CHEOPSim were used to generate light curves for simulated planetary transits across real and simulated targets. Independent analysts were successfully able to detect the planets and measure their radii to an accuracy within the science requirements of the mission.


title: Preferential concentration in the particle-induced convective instability

abstract: Heavy particles in turbulent flows have been shown to accumulate in regions of high strain rate or low vorticity, a process otherwise known as preferential concentration. This can be observed in geophysical flows, and is inferred to occur in astrophysical environments, often resulting in rapid particle growth which is critical to physical processes such as rain or planet formation. Here we study the effects of preferential concentration in a two-way coupled system in the context of the particle-driven convective instability. To do so, we use Direct Numerical Simulations and adopt the two-fluid approximation. We focus on a particle size range for which the latter is valid, namely when the Stokes number is $\lesssim O(0.1)$. For Stokes number above $\sim 0.01$, we find that the maximum particle concentration enhancement over the mean scales with the rms fluid velocity $u_{\rm{rms}}$, the particle stopping time $\tau_p$, and the particle diffusivity $\kappa_p$, as $u_{\rm{rms}}^2 \tau_p / \kappa_p$. We show that this scaling can be understood from simple arguments of dominant balance. We also show that the typical particle concentration enhancement over the mean scales as $(u_{\rm{rms}}^2 \tau_p/ \kappa_p)^{1/2}$. We finally find that the probability distribution function of the particle concentration enhancement over the mean has an exponential tail whose slope scales as $(u_{\rm{rms}}^2 \tau_p / \kappa_p)^{-1/2}$. We apply our model to geophysical and astrophysical examples, and discuss its limitations.


title: Variations on twins in permutations

abstract: Let $\pi$ be a permutation of the set $[n]=\{1,2,\dots, n\}$. Two disjoint order-isomorphic subsequences of $\pi$ are called twins. How long twins are contained in every permutation? The well known Erd\H{o}s-Szekeres theorem implies that there is always a pair of twins of length $\Omega(\sqrt{n})$. On the other hand, by a simple probabilistic argument Gawron proved that for every $n\geqslant 1$ there exist permutations with all twins having length $O(n^{2/3})$. He conjectured that the latter bound is the correct size of the longest twins guaranteed in every permutation. We support this conjecture by showing that almost all permutations contain twins of length $\Omega(n^{2/3}/\log n^{1/3})$. Recently, Bukh and Rudenko have tweaked our proof and removed the log-factor. For completeness, we also present our version of their proof (see Remark 1.2 below on the interrelation between the two proofs). In addition, we study several variants of the problem with diverse restrictions imposed on the twins. For instance, if we restrict attention to twins avoiding a fixed permutation $\tau$, then the corresponding extremal function equals $\Theta(\sqrt{n})$, provided that $\tau$ is not monotone. In case of block twins (each twin occupies a segment) we prove that it is $(1+o(1))\frac{\log n}{\log\log n}$, while for random permutations it is twice as large. For twins that jointly occupy a segment (tight twins), we prove that for every $n$ there are permutations avoiding them on all segments of length greater than $24$.


title: Improvement of heavy-heavy and heavy-light currents with the Oktay-Kronfeld action

abstract: The CKM matrix elements $V_{cb}$ and $V_{ub}$ can be obtained by combining data from the experiments with lattice QCD results for the semi-leptonic form factors for the $\bar{B} \to D^\ast \ell \bar{\nu}$ and $\bar{B} \to \pi \ell \bar{\nu}$ decays. It is highly desirable to use the Oktay-Kronfeld (OK) action for the form factor calculation on the lattice, since the OK action is designed to reduce the heavy quark discretization error down to the $\mathcal{O}(\lambda^4)$ level in the power counting rules of the heavy quark effective theory (HQET). Here, we present a matching calculation to improve heavy-heavy and heavy-light currents up to the $\lambda^3$ order in HQET, the same level of improvement as the OK action. Our final results for the improved currents are being used in a lattice QCD calculation of the semi-leptonic form factors for the $\bar{B} \to D^\ast \ell \bar{\nu}$ and $\bar{B} \to D \ell \bar{\nu}$ decays.


title: Distributed, partially collapsed MCMC for Bayesian Nonparametrics

abstract: Bayesian nonparametric (BNP) models provide elegant methods for discovering underlying latent features within a data set, but inference in such models can be slow. We exploit the fact that completely random measures, which commonly used models like the Dirichlet process and the beta-Bernoulli process can be expressed as, are decomposable into independent sub-measures. We use this decomposition to partition the latent measure into a finite measure containing only instantiated components, and an infinite measure containing all other components. We then select different inference algorithms for the two components: uncollapsed samplers mix well on the finite measure, while collapsed samplers mix well on the infinite, sparsely occupied tail. The resulting hybrid algorithm can be applied to a wide class of models, and can be easily distributed to allow scalable inference without sacrificing asymptotic convergence guarantees.


title: Boost-invariant description of polarization within hydrodynamics with spin

abstract: We briefly review a recently proposed formalism of perfect-fluid hydrodynamics with spin, which is a generalization of the standard hydrodynamic framework and provides a natural tool for describing the evolution of spin-polarized systems of particles with spin 1/2. It is based on the de Groot - van Leeuwen - van Weert forms of energy-momentum and spin tensors and conservation laws. Using Bjorken symmetry we show how this formalism may be used to determine observables describing the polarization of particles measured in the experiment.


title: The centre-quotient property and weak centrality for $C^*$-algebras

abstract: We give a number of equivalent conditions (including weak centrality) for a general $C^*$-algebra to have the centre-quotient property. We show that every $C^*$-algebra $A$ has a largest weakly central ideal $J_{wc}(A)$. For an ideal $I$ of a unital $C^*$-algebra $A$, we find a necessary and sufficient condition for a central element of $A/I$ to lift to a central element of $A$. This leads to a characterisation of the set $V_A$ of elements of an arbitrary $C^*$-algebra $A$ which prevent $A$ from having the centre-quotient property. The complement $\mathrm{CQ}(A):= A \setminus V_A$ always contains $Z(A)+J_{wc}(A)$ (where $Z(A)$ is the centre of $A$), with equality if and only if $A/J_{wc}(A)$ is abelian. Otherwise, $\mathrm{CQ}(A)$ fails spectacularly to be a $C^*$-subalgebra of $A$.


title: Thermal Emission Spectroscopy of Single, Isolated Carbon Nanoparticles: Effects of Particle Size, Material, Charge, Excitation Wavelength, and Thermal History

abstract: Results are presented for thermal emission from individually trapped carbon nanoparticles (NPs) in the temperature range from 1000 to 2100 K. We explore the effects on the magnitude and wavelength dependency of the emissivity, epsilon (lambda), of the NP size and charge and of the type of carbon material, including graphite, graphene, diamond, carbon black, and carbon dots. In addition, it is found that heating the NPs, particularly to temperatures above 1900 K, results in significant changes in the emission properties, which is attributed to changes in the distribution of surface and defect sites caused by annealing and sublimation.


title: A Cameron-Storvick type theorem on $C_{a,b}^2[0,T]$ with applications

abstract: The purpose of this paper is to establish a very general Cameron-Storvick theorem involving the generalized analytic Feynman integral of functionals on the product function space $C_{a,b}^2[0,T]$. The function space $C_{a,b}[0,T]$ can be induced by the generalized Brownian motion process associated with continuous functions $a$ and $b$. To do this we first introduce the class $\mathcal F_{A_1,A_2}^{\,\,a,b}$ of functionals on $C_{a,b}^2[0,T]$ which is a generalization of the Kallianpur and Bromley Fresnel class $\mathcal F_{A_1,A_2}$. We then proceed to establish a Cameron-Storvick type theorem on the product function space $C_{a,b}^2[0,T]$. Finally we use our Cameron--Storvick type theorem to obtain several meaningful results and examples.


title: Windows for cdgas

abstract: We study a Fourier-Mukai kernel associated to a GIT wall-crossing for arbitrarily singular (not necessarily reduced or irreducible) affine varieties over any field. This kernel is closely related to a derived fiber product diagram for the wall-crossing and simple to understand from the viewpoint of commutative differential graded algebras. However, from the perspective of algebraic varieties, the kernel can be quite complicated, corresponding to a complex with multiple homology sheaves. Under mild assumptions in the Calabi-Yau case, we prove that this kernel provides an equivalence between the category of perfect complexes on the two GIT quotients. More generally, we obtain semi-orthogonal decompositions which show that these categories differ by a certain number of copies of the derived category of the derived fixed locus. The derived equivalence for the Mukai flop is recovered as a very special case.


title: Stellar Characterization of M-dwarfs from the APOGEE Survey: A Calibrator Sample for the M-dwarf Metallicities

abstract: We present spectroscopic determinations of the effective temperatures, surface gravities and metallicities for 21 M-dwarfs observed at high-resolution (R $\sim$ 22,500) in the \textit{H}-band as part of the SDSS-IV APOGEE survey. The atmospheric parameters and metallicities are derived from spectral syntheses with 1-D LTE plane parallel MARCS models and the APOGEE atomic/molecular line list, together with up-to-date H$_{2}$O and FeH molecular line lists. Our sample range in $T_{\rm eff}$ from $\sim$ 3200 to 3800K, where eleven stars are in binary systems with a warmer (FGK) primary, while the other 10 M-dwarfs have interferometric radii in the literature. We define an $M_{K_{S}}$--Radius calibration based on our M-dwarf radii derived from the detailed analysis of APOGEE spectra and Gaia DR2 distances, as well as a mass-radius relation using the spectroscopically-derived surface gravities. A comparison of the derived radii with interferometric values from the literature finds that the spectroscopic radii are slightly offset towards smaller values, with $\Delta$ = -0.01 $\pm$ 0.02 $R{\star}$/$R_{\odot}$. In addition, the derived M-dwarf masses based upon the radii and surface gravities tend to be slightly smaller (by $\sim$5-10\%) than masses derived for M-dwarf members of eclipsing binary systems for a given stellar radius. The metallicities derived for the 11 M-dwarfs in binary systems, compared to metallicities obtained for their hotter FGK main-sequence primary stars from the literature, shows excellent agreement, with a mean difference of [Fe/H](M-dwarf - FGK primary) = +0.04 $\pm$ 0.18 dex, confirming the APOGEE metallicity scale derived here for M-dwarfs.


title: Banjo Break Angle Tension Modulation as Parametric Oscillation

abstract: The motion of the floating bridge of the banjo, in conjunction with the break angle of the strings over that bridge, produces string tension modulation that is first order in the amplitude of the string motion. This note refines a previous suggestion regarding the impact on the frequencies of the strings' and bridge's motion. For a given mode frequency pair of string and bridge, the resulting tension modulation produces a new, additional motion characterized by the sum and difference of the original ones. Strictly speaking, this corresponds to canonical "frequency modulation" only in the limit of modulation slow compared to the string frequency. The more general result is precisely an example of what is known as "parametric oscillation," first analyzed by Rayleigh. The qualitative impact of tension modulation on banjo timbre remains as suggested previously. It is only the precise math and physics that warrants this correction.


title: Semisimple flat F-manifolds in higher genus

abstract: In this paper, we generalize the Givental theory for Frobenius manifolds and cohomological field theories to flat F-manifolds and F-cohomological field theories. In particular, we define a notion of Givental cone for flat F-manifolds, and we provide a generalization of the Givental group as a matrix loop group acting on them. We show that this action is transitive on semisimple flat F-manifolds. We then extend this action to F-cohomological field theories in all genera. We show that, given a semisimple flat F-manifold and a Givental group element connecting it to the constant flat F-manifold at its origin, one can construct a family of F-CohFTs in all genera, parameterized by a vector in the associative algebra at the origin, whose genus $0$ part is the given flat F-manifold. If the flat F-manifold is homogeneous, then the associated family of F-CohFTs contains a subfamily of homogeneous F-CohFTs. However, unlike in the case of Frobenius manifolds and CohFTs, these homogeneous F-CohFTs can have different conformal dimensions, which are determined by the properties of a certain metric associated to the flat F-manifold.


title: Population synthesis of helium white dwarf--red giant star mergers and the formation of lithium-rich giants and carbon stars

abstract: The formation histories of lithium-rich and carbon-rich red giants are not yet understood. It has been proposed that the merger of a helium-core white dwarf with a red giant branch star might provide a solution. We have computed an extended grid of post-merger evolution models and combined these with predictions of binary star population synthesis. The results strongly support the proposal that the merger of a helium white dwarf with a red giant branch star can provide the progenitors of both lithium-rich red clump stars and early-R carbon stars. The distribution of post-merger models in $T_{\rm eff}$, $\log g$, $\log L$, the surface abundances of lithium and carbon, and the predicted space densities agree well with the observed distributions of these parameters for Li-rich and early-R stars in the Galaxy.


title: Revealing the Double Nucleus of NGC 4490

abstract: NGC 4490/85 (UGC 7651/48) or Arp 269 is well known for being one of the closest interacting/merging galactic systems. NGC 4490 has a high star formation rate (SFR) and is surrounded by an enormous HI feature stretching about 60 kpc north and south of the optically visible galaxies. Both the driver for the high SFR in NGC 4490 and the formation mechanism of the HI structure are puzzling aspects of this system. We have used mid-infrared Spitzer data to show that NGC 4490 has a double nucleus morphology. One nucleus is visible in the optical, while the other is only visible at infrared and radio wavelengths. We find the optical nucleus and the potential infrared visible nucleus have similar sizes, masses, and luminosities. Both are comparable in mass and luminosity to other nuclei found in interacting galaxy pairs and much more massive and luminous compared with typical non-nuclear star-forming complexes. We examine possible origin scenarios for the infrared feature, and conclude that it is likely that NGC 4490 is itself a merger remnant, which is now interacting with NGC 4485. This earlier encounter provides both a possible driver for extended star formation in NGC 4490, and multiple pathways for the formation of the extended HI plume.


title: Sequential Selection for Accelerated Life Testing via Approximate Bayesian Inference

abstract: Accelerated life testing (ALT) is typically used to assess the reliability of material's lifetime under desired stress levels. Recent advances in material engineering have made a variety of material alternatives readily available. To identify the most reliable material setting with efficient experimental design, a sequential test planning strategy is preferred. To guarantee a tractable statistical mechanism for information collection and update, we develop explicit model parameter update formulas via approximate Bayesian inference. Theories show that our explicit update formulas give consistent parameter estimates. Simulation study and a case study show that the proposed sequential selection approach can significantly improve the probability of identifying the material alternative with best reliability performance over other design approaches.


title: Resilient quantum electron microscopy

abstract: We investigate the fundamental limit of biological quantum electron microscopy, which is designed to go beyond the shot noise limit. Inelastic scattering is expected to be the main obstacle in this setting, especially for thick specimens of actual biological interest. Here we describe a measurement procedure that, in principle, significantly neutralizes the effect of inelastic scattering.


title: ALMA band 8 observations of DLA2233+131 at z=3.150

abstract: We present our ALMA Band 8 observations of a damped Ly${\alpha}$ absorption (DLA) system at $z$=3.150 observed in the spectrum of the quasar Q2233+131 at $z$=3.295. The optical counterpart of this DLA has been identified and it shows a double-peaked Ly${\alpha}$ emission line. Since one possible origin of DLAs at high redshift is an outflowing gas from star-forming galaxies, DLA2233+131 provides a good laboratory to investigate the nature of high-$z$ DLAs. Motivated by this, we have carried out ALMA band 8 observations to study the [C II] line in this system. However, we do not detect any significant emission line in the observed pass bands. Instead, we have serendipitously found three submm continuum sources in the observed sky area. One appears to be the quasar Q2233+131 itself while the other two sources are newly identified submm galaxies (SMGs), called SMG1 and SMG2 in this paper. They are located at a separation of 400.7 and 800.1 from Q2233+131, respectively. Their 646 ${\mu}$m fluxes are 6.35 mJy and 6.43 mJy, respectively, being higher than that of Q2233+131, 3.62 mJy. Since these two SMGs are not detected in the optical images obtained with the Hubble Space Telescope and the Subaru Telescope, they have a very red spectral energy distribution. It is, therefore,suggested that they are high-redshift galaxies or very dusty galaxies at intermediate redshift although we cannot rule out the possibility that they are optically very faint SMG analogs at low redshift. Follow up observations will be necessary to explore the nature of this interesting region.


title: Revisiting the detection rate for axion haloscopes

abstract: The cavity haloscope has been employed to detect microwave photons resonantly converted from invisible cosmic axions under a strong magnetic field. In this scheme, the axion-photon conversion power has been formulated to be valid for certain conditions, either $Q_{cavity}\ll Q_{\rm axion}$ or $Q_{cavity} \gg Q_{axion}$. This remedy, however, fails when these two quantities are comparable to each other. Furthermore, the noise power flow has been treated independently of the impedance mismatch of the system, which could give rise to misleading estimates of the experimental sensitivity. We revisit the analytical approaches to derive a general description of the signal and noise power. We also optimize the coupling strength of a receiver to yield the maximal sensitivity for axion search experiments.


title: Full field-of-view multi-targets imaging through scattering beyond 3D optical memory effect

abstract: A robust method and strategy for efficient full field-ofview and depth separation optical imaging through scattering media regardless of the three-dimensional (3D) optical memory effect are proposed. In this method, the problem of imaging de-aliasing, decomposition, and separation of speckle patterns are solved taking advantages of the spatial decorrelation characteristics of speckles by employing randomly modulated illumination strategy and independent component analysis methods. Full field-of-view imaging of multi-targets locate at diverse spatial positions behind a scattering layer are realized and observed experimentally, for the first time, to the best of our knowledge. The method and strategy provide a potentially useful means for incoherent imaging through scattering in a wide class of fields such as optical microscopy, biomedical imaging, and astronomical imaging.


title: NNSC-Cobordism of Bartnik Data in High Dimensions

abstract: In this short note, we formulate three problems relating to nonnegative scalar curvature (NNSC) fill-ins. Loosely speaking, the first two problems focus on: When are $(n-1)$-dimensional Bartnik data $\big(\Sigma_i ^{n-1}, \gamma_i, H_i\big)$, $i=1,2$, NNSC-cobordant? (i.e., there is an $n$-dimensional compact Riemannian manifold $\big(\Omega^n, g\big)$ with scalar curvature $R(g)\geq 0$ and the boundary $\partial \Omega=\Sigma_{1} \cup \Sigma_{2}$ such that $\gamma_i$ is the metric on $\Sigma_i ^{n-1}$ induced by $g$, and $H_i$ is the mean curvature of $\Sigma_i$ in $\big(\Omega^n, g\big)$). If $\big(\mathbb{S}^{n-1},\gamma_{\rm std},0\big)$ is positive scalar curvature (PSC) cobordant to $\big(\Sigma_1 ^{n-1}, \gamma_1, H_1\big)$, where $\big(\mathbb{S}^{n-1}, \gamma_{\rm std}\big)$ denotes the standard round unit sphere then $\big(\Sigma_1 ^{n-1}, \gamma_1, H_1\big)$ admits an NNSC fill-in. Just as Gromov's conjecture is connected with positive mass theorem, our problems are connected with Penrose inequality, at least in the case of $n=3$. Our third problem is on $\Lambda\big(\Sigma^{n-1}, \gamma\big)$ defined below.