title: Pairing Conceptual Modeling with Machine Learning

abstract: Both conceptual modeling and machine learning have long been recognized as important areas of research. With the increasing emphasis on digitizing and processing large amounts of data for business and other applications, it would be helpful to consider how these areas of research can complement each other. To understand how they can be paired, we provide an overview of machine learning foundations and development cycle. We then examine how conceptual modeling can be applied to machine learning and propose a framework for incorporating conceptual modeling into data science projects. The framework is illustrated by applying it to a healthcare application. For the inverse pairing, machine learning can impact conceptual modeling through text and rule mining, as well as knowledge graphs. The pairing of conceptual modeling and machine learning in this this way should help lay the foundations for future research.


title: Soliton solutions of the shifted nonlocal NLS and MKdV equations

abstract: We find one- and two-soliton solutions of shifted nonlocal NLS and MKdV equations. We discuss the singular structures of these soliton solutions and present some of the graphs of them.


title: An efficient and secure scheme of verifiable computation for Intel SGX

abstract: Cloud computing offers resource-constrained users big-volume data storage and energy-consuming complicated computation. However, owing to the lack of full trust in the cloud, the cloud users prefer privacy-preserving outsourced data computation with correctness verification. However, cryptography-based schemes introduce high computational costs to both the cloud and its users for verifiable computation with privacy preservation, which makes it difficult to support complicated computations in practice. Intel Software Guard Extensions (SGX) as a trusted execution environment is widely researched in various fields (such as secure data analytics and computation), and is regarded as a promising way to achieve efficient outsourced data computation with privacy preservation over the cloud. But we find two types of threats towards the computation with SGX: Disarranging Data-Related Code threat and Output Tampering and Misrouting threat. In this paper, we depict these threats using formal methods and successfully conduct the two threats on the enclave program constructed by Rust SGX SDK to demonstrate their impacts on the correctness of computations over SGX enclaves. In order to provide countermeasures, we propose an efficient and secure scheme to resist the threats and realize verifiable computation for Intel SGX. We prove the security and show the efficiency and correctness of our proposed scheme through theoretic analysis and extensive experiments. Furthermore, we compare the performance of our scheme with that of some cryptography-based schemes to show its high efficiency.


title: Pluri-subharmonic functions on complex tori, Ricci curvature and convexity

abstract: We show that, in toric manifolds, one can characterize the sign of the Ricci curvature in terms of the convexity of the volume functional. More generally we discuss relationships between (i) Ricci curvature and volume, (ii) totally real and Lagrangian submanifolds, (iii) pluri-subharmonic functions and convexity.


title: On Graphical Models and Convex Geometry

abstract: We introduce a mixture-model of beta distributions to identify significant correlations among $P$ predictors when $P$ is large. The method relies on theorems in convex geometry, which we use to show how to control the error rate of edge detection in graphical models. Our `betaMix' method does not require any assumptions about the network structure, nor does it assume that the network is sparse. The results in this article hold for a wide class of data generating distributions that include light-tailed and heavy-tailed spherically symmetric distributions.


title: Using deep learning to detect patients at risk for prostate cancer despite benign biopsies

abstract: Background: Transrectal ultrasound guided systematic biopsies of the prostate is a routine procedure to establish a prostate cancer diagnosis. However, the 10-12 prostate core biopsies only sample a relatively small volume of the prostate, and tumour lesions in regions between biopsy cores can be missed, leading to a well-known low sensitivity to detect clinically relevant cancer. As a proof-of-principle, we developed and validated a deep convolutional neural network model to distinguish between morphological patterns in benign prostate biopsy whole slide images from men with and without established cancer. Methods: This study included 14,354 hematoxylin and eosin stained whole slide images from benign prostate biopsies from 1,508 men in two groups: men without an established prostate cancer (PCa) diagnosis and men with at least one core biopsy diagnosed with PCa. 80% of the participants were assigned as training data and used for model optimization (1,211 men), and the remaining 20% (297 men) as a held-out test set used to evaluate model performance. An ensemble of 10 deep convolutional neural network models was optimized for classification of biopsies from men with and without established cancer. Hyperparameter optimization and model selection was performed by cross-validation in the training data . Results: Area under the receiver operating characteristic curve (ROC-AUC) was estimated as 0.727 (bootstrap 95% CI: 0.708-0.745) on biopsy level and 0.738 (bootstrap 95% CI: 0.682 - 0.796) on man level. At a specificity of 0.9 the model had an estimated sensitivity of 0.348. Conclusion: The developed model has the ability to detect men with risk of missed PCa due to under-sampling of the prostate. The proposed model has the potential to reduce the number of false negative cases in routine systematic prostate biopsies and to indicate men who could benefit from MRI-guided re-biopsy.


title: Use of Machine Learning Technique to maximize the signal over background for $H \rightarrow \tau \tau$

abstract: In recent years, artificial neural networks (ANNs) have won numerous contests in pattern recognition and machine learning. ANNS have been applied to problems ranging from speech recognition to prediction of protein secondary structure, classification of cancers, and gene prediction. Here, we intend to maximize the chances of finding the Higgs boson decays to two $\tau$ leptons in the pseudo dataset using a Machine Learning technique to classify the recorded events as signal or background.


title: Sparse Logistic Tensor Decomposition for Binary Data

abstract: Tensor data are increasingly available in many application domains. We develop several tensor decomposition methods for binary tensor data. Different from classical tensor decompositions for continuous-valued data with squared error loss, we formulate logistic tensor decompositions for binary data with a Bernoulli likelihood. To enhance the interpretability of estimated factors and improve their stability further, we propose sparse formulations of logistic tensor decomposition by considering $\ell_{1}$-norm and $\ell_{0}$-norm regularized likelihood. To handle the resulting optimization problems, we develop computational algorithms which combine the strengths of tensor power method and majorization-minimization (MM) algorithm. Through simulation studies, we demonstrate the utility of our methods in analysis of binary tensor data. To illustrate the effectiveness of the proposed methods, we analyze a dataset concerning nations and their political relations and perform co-clustering of estimated factors to find associations between the nations and political relations.


title: SDOF-Tracker: Fast and Accurate Multiple Human Tracking by Skipped-Detection and Optical-Flow

abstract: Multiple human tracking is a fundamental problem for scene understanding. Although both accuracy and speed are required in real-world applications, recent tracking methods based on deep learning have focused on accuracy and require substantial running time. This study aims to improve running speed by performing human detection at a certain frame interval because it accounts for most of the running time. The question is how to maintain accuracy while skipping human detection. In this paper, we propose a method that complements the detection results with optical flow, based on the fact that someone's appearance does not change much between adjacent frames. To maintain the tracking accuracy, we introduce robust interest point selection within human regions and a tracking termination metric calculated by the distribution of the interest points. On the MOT20 dataset in the MOTChallenge, the proposed SDOF-Tracker achieved the best performance in terms of the total running speed while maintaining the MOTA metric. Our code is available at https://github.com/hitottiez/sdof-tracker.


title: Nondegenerate bright solitons in coupled nonlinear Schr\"{o}dinger systems: Recent developments on optical vector solitons

abstract: Nonlinear dynamics of an optical pulse or a beam continue to be one of the active areas of research in the field of optical solitons. Especially, in multi-mode fibers or fiber arrays and photorefractive materials, the vector solitons display rich nonlinear phenomena. Due to their fascinating and intriguing novel properties, the theory of optical vector solitons has been developed considerably both from theoretical and experimental points of view leading to soliton based promising potential applications. In the recent past, many types of vector solitons have been identified both in the integrable and non-integrable coupled nonlinear Schr\"{o}dinger (CNLS) equations framework. In this article, we review some of the recent progress in understanding the dynamics of the so called nondegenerate vector bright solitons in nonlinear optics, where the fundamental soliton can have more than one propagation constant. We address this theme by considering the integrable two CNLS family of equations, namely Manakov system, mixed 2-CNLS system, coherently CNLS system, generalized CNLS system and two-component long-wave short-wave resonance interaction (LSRI) system. In these models, we discuss the existence of nondegenerate vector solitons and their associated novel multi-hump geometrical profile nature by deriving their analytical forms through the Hirota bilinear method. Then we reveal the novel collision properties of the nondegenerate solitons in the Manakov system as an example. The asymptotic analysis shows that the nondegenerate solitons, in general, undergo three types of elastic collisions without any energy redistribution among the modes. Further, we show that the energy sharing collision exhibiting vector solitons arises as a special case of the newly reported nondegenerate vector solitons. Finally, we point out the possible further developments in this subject and potential applications.


title: Pulsating iron spectral features in the emission of X-ray pulsar V0332+53

abstract: We present results of phase- and time-resolved study of iron spectral features in the emission of the Be/X-ray transient pulsar V0332+53, during its type II outburst in 2004 using archival RXTE/PCA data. Coherent pulsations of both fluorescent iron line at {\approx}6.4 keV and neutral iron K-edge at {\approx}7.1 keV have been detected throughout the entire outburst. The pulsating iron K-edge is reported for the first time for this object. Near the peak of the outburst, the 3-12 keV pulse profile shows two deep, F_max / F_min ~ 2, and narrow dips of nearly identical shape, separated by exactly {\Delta}{\phi}=0.5 in phase. The dip spectra are nearly identical to each other and very similar in shape to the spectra outside the dips. The iron K-edge peaks at the phase intervals corresponding to the dips, although its optical depth {\tau}_K ~ 0.05 is by far insufficient to explain the dips. The iron line shows pulsations with a complex pulse profile without any obvious correlation with the total flux or optical depth of the K-edge. Accounting for the component associated with reprocessing of the pulsar emission by the surface of the donor star and circumstellar material, we find a very high pulsation amplitude of the iron line flux, F_max / F_min ~ 10. We demonstrate that these properties of V0332+53 can not be explained by contemporary emission models for accreting X-ray pulsars and speculate about the origin of the observed iron spectral features.


title: Edge-Unfolding Prismatoids: Tall or Rectangular Base

abstract: We show how to edge-unfold a new class of convex polyhedra, specifically a new class of prismatoids (the convex hull of two parallel convex polygons, called the top and base), by constructing a nonoverlapping "petal unfolding" in two new cases: (1) when the top and base are sufficiently far from each other; and (2) when the base is a rectangle and all other faces are nonobtuse triangles. The latter result extends a previous result by O'Rourke that the petal unfolding of a prismatoid avoids overlap when the base is a triangle (possibly obtuse) and all other faces are nonobtuse triangles. We also illustrate the difficulty of extending this result to a general quadrilateral base by giving a counterexample to our technique.


title: The Role of Shift Vector in High-Harmonic Generation from Non-Centrosymmetric Topological Insulators under Strong Laser Fields

abstract: As a promising avenue to obtain new extreme ultraviolet light source and detect electronic properties, high-harmonic generation (HHG) has been actively developed in both theory and experiment. In solids lacking inversion symmetry, when electrons undergo a nonadiabatic transition, a directional charge shift occurs and is characterized by shift vector, which measures the real-space shift of the photoexcited electron and hole. For the first time, we have revealed that shift vector plays prominent roles in the real-space tunneling mechanism of three-step model for electrons under strong laser fields. Since shift vector is determined by the topological properties of related wave functions, we expect HHG with its contribution can provide direct knowledge on the band topology in noncentrosymmetric topological insulators (TIs). In both Kane-Mele model and realistic material BiTeI, we have found that the shift vector reverses when band inversion happens during the topological phase transition between normal and topological insulators. Under oscillating strong laser fields, the reversal of shift vector leads to completely opposite radiation time of high-order harmonics. This makes HHG a feasible all-optical strong-field method to directly identify the band inversion in non-centrosymmetric TIs.


title: Geometric fields and new enantio-sensitive observables in photoionization of chiral molecules

abstract: Chiral molecules are instrumental for molecular recognition in living organisms. Distinguishing between two opposite enantiomers, the mirror twins of the same chiral molecule, is both vital and challenging. Photoelectron circular dichroism (PECD), an extremely sensitive probe of molecular chirality via photoionization, outperforms standard optical methods by many orders of magnitude. Here we show that the physical origin of PECD in chiral molecules is linked to the concept of geometric magnetism, which enables a broad class of phenomena in solids including the anomalous electron velocity, the Hall effect, and related topological phenomena. We uncover the geometric field in molecular photoionization, which leads to a new class of enantio-sensitive observables emerging due to ultrafast excitation of chiral electronic or vibronic currents prior to ionization. Next, we introduce the first member of this new class: enantio-sensitive orientation of chiral molecules via photoionization. This effect opens new routes to both enantio-separation and imaging of chiral dynamics on ultrafast time scales. Our work suggests that geometric fields in photoionization provide the bridge between the two geometrical properties, chirality and topology.


title: Reward-Based 1-bit Compressed Federated Distillation on Blockchain

abstract: The recent advent of various forms of Federated Knowledge Distillation (FD) paves the way for a new generation of robust and communication-efficient Federated Learning (FL), where mere soft-labels are aggregated, rather than whole gradients of Deep Neural Networks (DNN) as done in previous FL schemes. This security-per-design approach in combination with increasingly performant Internet of Things (IoT) and mobile devices opens up a new realm of possibilities to utilize private data from industries as well as from individuals as input for artificial intelligence model training. Yet in previous FL systems, lack of trust due to the imbalance of power between workers and a central authority, the assumption of altruistic worker participation and the inability to correctly measure and compare contributions of workers hinder this technology from scaling beyond small groups of already entrusted entities towards mass adoption. This work aims to mitigate the aforementioned issues by introducing a novel decentralized federated learning framework where heavily compressed 1-bit soft-labels, resembling 1-hot label predictions, are aggregated on a smart contract. In a context where workers' contributions are now easily comparable, we modify the Peer Truth Serum for Crowdsourcing mechanism (PTSC) for FD to reward honest participation based on peer consistency in an incentive compatible fashion. Due to heavy reductions of both computational complexity and storage, our framework is a fully on-blockchain FL system that is feasible on simple smart contracts and therefore blockchain agnostic. We experimentally test our new framework and validate its theoretical properties.


title: Metric description of flexible octahedra

abstract: For flexibility of an octahedron we find necessary metric conditions in terms of edge lengths. These conditions yield a new description of Bricard's octahedra, suitable for solving some problems in metric geometry of octahedra, in particular, for searching the proof of I.\,Hh~Sabitov hypothesis that all non-leading coefficients of the volume polynomial for an octahedron of third type are zero.


title: Revisiting the distance to radio Loops I and IV using Gaia and radio/optical polarization data

abstract: Galactic synchrotron emission exhibits large-angular-scale features known as radio spurs and loops. Determining the physical size of these structures is important for understanding the local interstellar structure and for modeling the Galactic magnetic field. However, the distance to these structures is either under debate or entirely unknown. We revisit a classical method of finding the location of radio spurs by comparing optical polarization angles with those of synchrotron emission as a function of distance. We consider three tracers of the magnetic field: stellar polarization, polarized synchrotron radio emission, and polarized thermal dust emission. We employ archival measurements of optical starlight polarization and Gaia distances, and construct a new map of polarized synchrotron emission from WMAP and Planck data. We confirm that synchrotron, dust emission, and stellar polarization angles all show a statistically significant alignment at high Galactic latitude. We obtain distance limits to three regions towards Loop I of 112$\pm$17 pc, 135$\pm$20 pc, and $<105$ pc. Our results strongly suggest that the polarized synchrotron emission towards the North Polar Spur at $b > 30^\circ$ is local. This is consistent with the conclusions of earlier work based on stellar polarization and extinction, but in stark contrast with the Galactic center origin recently revisited on the basis of X-ray data. We also obtain a distance measurement towards part of Loop IV (180$\pm$15 pc) and find evidence that its synchrotron emission arises from chance overlap of structures located at different distances. Future optical polarization surveys will allow the expansion of this analysis to other radio spurs.


title: Brauer Groups of Algebraic Stacks and GIT-Quotients

abstract: In this paper we consider the Brauer groups of algebraic stacks and GIT quotients: the only algebraic stacks we consider in this paper are quotient stacks [X/G], where X is a smooth scheme of finite type over a field k, and G is a linear algebraic group over k and acting on X, as well as various moduli stacks of principal G-bundles on a smooth projective curve X, associated to a reductive group G. We also consider the Brauer groups of the corresponding coarse moduli spaces, which most often identify with the corresponding GIT-quotients. One conclusion that we seem to draw then is that the Brauer groups (or their $\ell$-primary torsion parts, for a fixed prime $\ell$ different from char(k)) of the corresponding stacks and coarse moduli spaces depend strongly on the Brauer group of the given scheme X.


title: Deep Learning for Technical Document Classification

abstract: In large technology companies, the requirements for managing and organizing technical documents created by engineers and managers have increased dramatically in recent years, which has led to a higher demand for more scalable, accurate, and automated document classification. Prior studies have only focused on processing text for classification, whereas technical documents often contain multimodal information. To leverage multimodal information for document classification to improve the model performance, this paper presents a novel multimodal deep learning architecture, TechDoc, which utilizes three types of information, including natural language texts and descriptive images within documents and the associations among the documents. The architecture synthesizes the convolutional neural network, recurrent neural network, and graph neural network through an integrated training process. We applied the architecture to a large multimodal technical document database and trained the model for classifying documents based on the hierarchical International Patent Classification system. Our results show that TechDoc presents a greater classification accuracy than the unimodal methods and other state-of-the-art benchmarks. The trained model can potentially be scaled to millions of real-world multimodal technical documents, which is useful for data and knowledge management in large technology companies and organizations.


title: Relativistic correction of the Coulomb interaction in the local density approximation for energies and radii in doubly-magic nuclei

abstract: Effects of the relativistic correction of the Coulomb interaction on doubly-magic nuclei are discussed with Skyrme Hartree--Fock calculations. The relativistic correction is treated by using the local density approximation. It is found that the correction to the total energy is about $ 2.4 \, \mathrm{MeV} $ for $ {}^{208} \mathrm{Pb} $, while proton and neutron radii do not change significantly. This difference is larger than the difference of the Coulomb exchange (Fock) energy calculated with the local density (Hartree--Fock--Slater) approximation and that with the exact treatment and the neutron finite-size effect. Effects of the correction are also compared to the correction due to the vacuum polarization. It is shown that the two contributions to the total energy are comparable in light nuclei, but the latter dominates in heavy nuclei, while the contribution of the relativistic correction to the total energy is non-negligible compared to the target accuracy of the DFT calculation.


title: Relaxation in Ordered Assembly of Magnetic Nanoparticles

abstract: We study the relaxation characteristics in the two-dimensional ($l^{}_x \times l^{}_y$) array of magnetic nanoparticles (MNPs) as a function of aspect ratio $A^{}_r=l^{}_y/l^{}_x$, dipolar interaction strength $h^{}_d$ and anisotropy axis orientation using computer simulation. The anisotropy axes of all the MNPs are assumed to have the same direction, $\alpha$ being the orientational angle. Irrespective of $\alpha$ and $A^{}_r$, the functional form of the magnetization-decay curve is perfectly exponentially decaying with $h^{}_d\leq0.2$. There exists a transition in relaxation behaviour at $h^{}_d\approx0.4$; magnetization relaxes slowly for $\alpha\leq45^\circ$; it relaxes rapildy with $\alpha>45^\circ$. Interestingly, it decays rapidly for $h^{}_d>0.6$, irrespective of $\alpha$. It is because the dipolar interaction promotes antiferromagnetic coupling in such cases. There is a strong effect of $\alpha$ on the magnetic relaxation in the highly anisotropic system ($A^{}_r\geq25$). Interesting physics unfolds in the case of a huge aspect ratio $A^{}_r=400$. There is a rapid decay of magnetization with $\alpha$, even for weakly interacting MNPs. Remarkably, magnetization does not relax even with a moderate value of $h^{}_d=0.4$ and $\alpha=0^\circ$ because of ferromagnetic coupling dominance. Surprisingly, there is a complete magnetization reversal from saturation (+1) to $-1$ state with $\alpha>60^\circ$. The dipolar field and anisotropy axis tend to get aligned antiparallel to each other in such a case. The effective N\'eel relaxation time $\tau^{}_N$ depends weakly on $\alpha$ for small $h^{}_d$ and $A^{}_r\leq25.0$. For large $A^{}_r$, there is a rapid fall in $\tau^{}_N$ as $\alpha$ is incremented from 0 to $90^\circ$. These results benefit applications in data and energy storages where such controlled magnetization alignment and desired structural anisotropy are desirable.


title: Anomalously Strong Near-Neighbor Attraction in Doped 1D Cuprate Chains

abstract: In the cuprates, one-dimensional chain compounds provide a unique opportunity to understand the microscopic physics due to the availability of reliable theories. However, progress has been limited by the inability to controllably dope these materials. Here, we report the synthesis and spectroscopic analysis of the one-dimensional cuprate Ba$_{2-x}$Sr$_x$CuO$_{3+\delta}$ over a wide range of hole doping. Our angle-resolved photoemission experiments reveal the doping evolution of the holon and spinon branches. We identify a prominent folding branch whose intensity fails to match predictions of the simple Hubbard model. An additional strong near-neighbor attraction, which may arise from coupling to phonons, quantitatively explains experiments for all accessible doping levels. Considering structural and quantum chemistry similarities among cuprates, this attraction will play a similarly crucial role in the high-$T_C$ superconducting counterparts


title: A Systematic Review of Bio-Cyber Interface Technologies and Security Issues for Internet of Bio-Nano Things

abstract: Advances in synthetic biology and nanotechnology have contributed to the design of tools that can be used to control, reuse, modify, and re-engineer cells' structure, as well as enabling engineers to effectively use biological cells as programmable substrates to realize Bio-Nano Things (biological embedded computing devices). Bio-NanoThings are generally tiny, non-intrusive, and concealable devices that can be used for in-vivo applications such as intra-body sensing and actuation networks, where the use of artificial devices can be detrimental. Such (nano-scale) devices can be used in various healthcare settings such as continuous health monitoring, targeted drug delivery, and nano-surgeries. These services can also be grouped to form a collaborative network (i.e., nanonetwork), whose performance can potentially be improved when connected to higher bandwidth external networks such as the Internet, say via 5G. However, to realize the IoBNT paradigm, it is also important to seamlessly connect the biological environment with the technological landscape by having a dynamic interface design to convert biochemical signals from the human body into an equivalent electromagnetic signal (and vice versa). This, unfortunately, risks the exposure of internal biological mechanisms to cyber-based sensing and medical actuation, with potential security and privacy implications. This paper comprehensively reviews bio-cyber interface for IoBNT architecture, focusing on bio-cyber interfacing options for IoBNT like biologically inspired bio-electronic devices, RFID enabled implantable chips, and electronic tattoos. This study also identifies known and potential security and privacy vulnerabilities and mitigation strategies for consideration in future IoBNT designs and implementations.


title: Learning Mesh Representations via Binary Space Partitioning Tree Networks

abstract: Polygonal meshes are ubiquitous, but have only played a relatively minor role in the deep learning revolution. State-of-the-art neural generative models for 3D shapes learn implicit functions and generate meshes via expensive iso-surfacing. We overcome these challenges by employing a classical spatial data structure from computer graphics, Binary Space Partitioning (BSP), to facilitate 3D learning. The core operation of BSP involves recursive subdivision of 3D space to obtain convex sets. By exploiting this property, we devise BSP-Net, a network that learns to represent a 3D shape via convex decomposition without supervision. The network is trained to reconstruct a shape using a set of convexes obtained from a BSP-tree built over a set of planes, where the planes and convexes are both defined by learned network weights. BSP-Net directly outputs polygonal meshes from the inferred convexes. The generated meshes are watertight, compact (i.e., low-poly), and well suited to represent sharp geometry. We show that the reconstruction quality by BSP-Net is competitive with those from state-of-the-art methods while using much fewer primitives. We also explore variations to BSP-Net including using a more generic decoder for reconstruction, more general primitives than planes, as well as training a generative model with variational auto-encoders. Code is available at https://github.com/czq142857/BSP-NET-original.


title: Learning without Forgetting for 3D Point Cloud Objects

abstract: When we fine-tune a well-trained deep learning model for a new set of classes, the network learns new concepts but gradually forgets the knowledge of old training. In some real-life applications, we may be interested in learning new classes without forgetting the capability of previous experience. Such learning without forgetting problem is often investigated using 2D image recognition tasks. In this paper, considering the growth of depth camera technology, we address the same problem for the 3D point cloud object data. This problem becomes more challenging in the 3D domain than 2D because of the unavailability of large datasets and powerful pretrained backbone models. We investigate knowledge distillation techniques on 3D data to reduce catastrophic forgetting of the previous training. Moreover, we improve the distillation process by using semantic word vectors of object classes. We observe that exploring the interrelation of old and new knowledge during training helps to learn new concepts without forgetting old ones. Experimenting on three 3D point cloud recognition backbones (PointNet, DGCNN, and PointConv) and synthetic (ModelNet40, ModelNet10) and real scanned (ScanObjectNN) datasets, we establish new baseline results on learning without forgetting for 3D data. This research will instigate many future works in this area.


title: Disconnectivity-Aware Energy-Efficient Cargo-UAV Trajectory Planning with Minimum Handoffs

abstract: On-board battery consumption, cellular disconnectivity, and frequent handoff are key challenges for unmanned aerial vehicle (UAV) based delivery missions, a.k.a., cargo-UAV. Indeed, with the introduction of UAV technology into cargo shipping and logistics, designing energy-efficient paths becomes a serious issue for the next retail industry transformation. Typically, the latter has to guarantee uninterrupted or slightly interrupted cellular connectivity for the UAV's command and control through a small number of handoffs. In this paper, we formulate the trajectory planning as a multi-objective problem aiming to minimize both the UAV's energy consumption and the handoff rate, constrained by the UAV battery size and disconnectivity rate. Due to the problem's complexity, we propose a dynamic programming based solution. Through simulations, we demonstrate the efficiency of our approach in providing optimized UAV trajectories. Also, the impact of several parameters, such as the cargo-UAV altitude, disconnectivity rate, and type of environment, are investigated. The obtained results allow to draw recommendations and guidelines for cargo-UAV operations.


title: How many moments does MMD compare?

abstract: We present a new way of study of Mercer kernels, by corresponding to a special kernel $K$ a pseudo-differential operator $p({\mathbf x}, D)$ such that $\mathcal{F} p({\mathbf x}, D)^\dag p({\mathbf x}, D) \mathcal{F}^{-1}$ acts on smooth functions in the same way as an integral operator associated with $K$ (where $\mathcal{F}$ is the Fourier transform). We show that kernels defined by pseudo-differential operators are able to approximate uniformly any continuous Mercer kernel on a compact set. The symbol $p({\mathbf x}, {\mathbf y})$ encapsulates a lot of useful information about the structure of the Maximum Mean Discrepancy distance defined by the kernel $K$. We approximate $p({\mathbf x}, {\mathbf y})$ with the sum of the first $r$ terms of the Singular Value Decomposition of $p$, denoted by $p_r({\mathbf x}, {\mathbf y})$. If ordered singular values of the integral operator associated with $p({\mathbf x}, {\mathbf y})$ die down rapidly, the MMD distance defined by the new symbol $p_r$ differs from the initial one only slightly. Moreover, the new MMD distance can be interpreted as an aggregated result of comparing $r$ local moments of two probability distributions. The latter results holds under the condition that right singular vectors of the integral operator associated with $p$ are uniformly bounded. But even if this is not satisfied we can still hold that the Hilbert-Schmidt distance between $p$ and $p_r$ vanishes. Thus, we report an interesting phenomenon: the MMD distance measures the difference of two probability distributions with respect to a certain number of local moments, $r^\ast$, and this number $r^\ast$ depends on the speed with which singular values of $p$ die down.


title: Degrees of Faithful Irreducible Representations of Certain Metabelian Groups and a Question of Sim

abstract: In this paper, we answer affirmatively a question of H S Sim on representations in characteristic $0$, for a class of metabelian groups. Moreover, we provide examples to point out that the analogous answer is no longer valid if the solvable group has derived length larger than 2. Let $F$ be a field of characteristic $0$ and $\overline{F}$ be its algebraic closure. We prove that if $G$ is a finite metabelian group containing a maximal abelian normal subgroup which is a p-group with abelian quotient, all possible faithful irreducible representations over $F$ have the same degree and that the Schur index of any faithful irreducible $\overline{F}$-representation with respect to $F$ is always $1$ or $2$. H S Sim had proven such a result for metacyclic groups when the characteristic of $F$ is positive and posed the question in characteristic $0$. Our result answers this question for the above class of metabelian groups affirmatively. We also determine explicitly the Wedderburn component corresponding to any faithful irreducible $\overline{F}$-representation in the group algebra $F[G]$.


title: Numerical dispersion effects on the energy cascade in large-eddy simulation

abstract: Implicitly filtered large-eddy simulation (LES) is by nature numerically under-resolved. With the sole exception of Fourier-spectral methods, discrete numerical derivative operators cannot accurately represent the dynamics of all of the represented scales. Since the resolution scale in an LES usually lies in the inertial range, these poorly represented scales are dynamically significant and errors in their dynamics can affect all resolved scales. This Letter is focused on characterizing the effects of numerical dispersion error by studying the energy cascade in LES of convecting homogeneous isotropic turbulence. Numerical energy and transfer spectra reveal that energy is not transferred at the appropriate rate to wavemodes where significant dispersion error is present. This leads to a deficiency of energy in highly dispersive modes and an accompanying pile up of energy in the well resolved modes, since dissipation by the subgrid model is diminished. An asymptotic analysis indicates that dispersion error causes a phase decoherence between triad interacting wavemodes, leading to a reduction in the mean energy transfer rate for these scales. These findings are relevant to a wide range of LES, since turbulence commonly convects through the grid in practical simulations. Further, these results indicate that the resolved scales should be defined to not include the dispersive modes.


title: Algorithmic Randomness and Kolmogorov Complexity for Qubits

abstract: Nies and Scholz defined quantum Martin-L\"of randomness (q-MLR) for states (infinite qubitstrings). We define a notion of quantum Solovay randomness and show it to be equivalent to q-MLR using purely linear algebraic methods. Quantum Schnorr randomness is then introduced. A quantum analogue of the law of large numbers is shown to hold for quantum Schnorr random states. We introduce quantum-K, ($QK$) a measure of the descriptive complexity of density matrices using classical prefix-free Turing machines and show that the initial segments of weak Solovay random and quantum Schnorr random states are incompressible in the sense of $QK$. Several connections between Solovay randomness and $K$ carry over to those between weak Solovay randomness and $QK$. We then define $QK_C$, using computable measure machines and connect it to quantum Schnorr randomness. We then explore a notion of `measuring' a state. We formalize how `measurement' of a state induces a probability measure on the space of infinite bitstrings. A state is `measurement random' ($mR$) if the measure induced by it, under any computable basis, assigns probability one to the set of Martin-L\"of randoms. I.e., measuring a $mR$ state produces a Martin-L\"of random bitstring almost surely. While quantum-Martin-L\"of random states are $mR$, the converse fails: there is a $mR$ state, $\rho$ which is not quantum-Martin-L\"of random. In fact, something stronger is true. While $\rho$ is computable and can be easily constructed, measuring it in any computable basis yields an arithmetically random sequence with probability one. So, classical randomness can be generated from a computable state which is not quantum random. We conclude by studying the asymptotic von Neumann entropy of computable states.


title: The Einstein ring GAL-CLUS-022058s: a Lensed Ultrabright Submillimeter Galaxy at z=1.4796

abstract: We report an ultra-bright lensed submillimeter galaxy at $z_{spec}=1.4796$, identified as a result of a full-sky cross-correlation of the AllWISE and Planck compact source catalogs aimed to search for bright submillimeter galaxies at $z \sim 1.5-2.8$. APEX/LABOCA observations of the candidate galaxy reveal a source with flux (S$_{870 \mu m}= 54\pm 8$ mJy). The position of the APEX source coincides with the position of the AllWISE mid-IR source, and with the Einstein ring GAL-CLUS-022058s, observed with the HST. Archival VLT/FORS observations reveal the redshift of this Einstein ring, $z_{spec}=1.4796$, and detection of the CO(5-4) line at $z_{spec} = 1.4802$ with APEX/nFLASH230 confirms the redshift of the submillimeter emission. The lensed source appears to be gravitationally magnified by a massive foreground galaxy cluster lens at $z = 0.36$. We use Lenstool to model the gravitational lensing, which is near to a "fold arc" configuration for an elliptical mass distribution of the central halo, where four images of the lensed galaxy are seen; the mean magnification is $\mu_{\rm L} =18\pm 4$. We have determined an intrinsic rest-frame infrared luminosity of $L_{IR} \approx 10^{12} L_\odot $ and a likely star formation rate of $\sim 70-170$ $M_\odot\ yr^{-1}$. The molecular gas mass is $M_{mol} \sim 2.6 \times 10^{10} M_\odot$ and the gas fraction is $f = 0.34\pm 0.07$. We also obtain a stellar mass log$(M_\ast/M_\odot) = 10.7 \pm 0.1$ and a specific star formation rate log$(sSFR/Gyr^{-1})=0.15 \pm 0.03$. This galaxy lies on the so-called main sequence of star-forming galaxies at this redshift.


title: A Closer Look at How Fine-tuning Changes BERT

abstract: Given the prevalence of pre-trained contextualized representations in today's NLP, there have been many efforts to understand what information they contain, and why they seem to be universally successful. The most common approach to use these representations involves fine-tuning them for an end task. Yet, how fine-tuning changes the underlying embedding space is less studied. In this work, we study the English BERT family and use two probing techniques to analyze how fine-tuning changes the space. We hypothesize that fine-tuning affects classification performance by increasing the distances between examples associated with different labels. We confirm this hypothesis with carefully designed experiments on five different NLP tasks. Via these experiments, we also discover an exception to the prevailing wisdom that "fine-tuning always improves performance". Finally, by comparing the representations before and after fine-tuning, we discover that fine-tuning does not introduce arbitrary changes to representations; instead, it adjusts the representations to downstream tasks while largely preserving the original spatial structure of the data points.


title: Bilateral Birth and death process in quantum calculus

abstract: In this paper I shall give the complete solution of the equations governing the bilateral birth and death process on path set $\mathbb{R}_q=\{q^n,\quad n\in\mathbb{Z}\}$ in which the birth and death rates $\lambda_n=q^{2\nu-2n}$ and $\mu_n=q^{-2n}$ where $0<q<1$ and $\nu>-1$ . The mathematical methods employed here are based on $q$-Bessel Fourier analysis.


title: Effective grey-box testing with partial FSM models

abstract: For partial, nondeterministic, finite state machines, a new conformance relation called strong reduction is presented. It complements other existing conformance relations in the sense that the new relation is well-suited for model-based testing of systems whose inputs are enabled or disabled, depending on the actual system state. Examples of such systems are graphical user interfaces and systems with interfaces that can be enabled or disabled in a mechanical way. We present a new test generation algorithm producing complete test suites for strong reduction. The suites are executed according to the grey-box testing paradigm: it is assumed that the state-dependent sets of enabled inputs can be identified during test execution, while the implementation states remain hidden, as in black-box testing. It is shown that this grey-box information is exploited by the generation algorithm in such a way that the resulting best-case test suite size is only linear in the state space size of the reference model. Moreover, examples show that this may lead to significant reductions of test suite size in comparison to true black-box testing for strong reduction.


title: Twin Vertices in Fault-Tolerant Metric Sets and Fault-Tolerant Metric Dimension of Multistage Interconnection Networks

abstract: A set of vertices $S\subseteq V(G)$ is a basis or resolving set of a graph $G$ if for each $x,y\in V(G)$ there is a vertex $u\in S$ such that $d(x,u)\neq d(y,u)$. A basis $S$ is a fault-tolerant basis if $S\setminus \{x\}$ is a basis for every $x \in S$. The fault-tolerant metric dimension (FTMD) $\beta'(G)$ of $G$ is the minimum cardinality of a fault-tolerant basis. It is shown that each twin vertex of $G$ belongs to every fault-tolerant basis of $G$. As a consequence, $\beta'(G) = n(G)$ iff each vertex of $G$ is a twin vertex, which corrects a wrong characterization of graphs $G$ with $\beta'(G) = n(G)$ from [Mathematics 7(1) (2019) 78]. This FTMD problem is reinvestigated for Butterfly networks, Benes networks, and silicate networks. This extends partial results from [IEEE Access 8 (2020) 145435--145445], and at the same time, disproves related conjectures from the same paper.


title: Exact Simulation of Pigment-Protein Complexes Unveils Vibronic Renormalization of Electronic Parameters in Ultrafast Spectroscopy

abstract: The primary steps of photosynthesis rely on the generation, transport, and trapping of excitons in pigment-protein complexes (PPCs). Generically, PPCs possess highly structured vibrational spectra, combining many discrete intra-pigment modes and a quasi-continuous of protein modes, with vibrational and electronic couplings of comparable strength. The intricacy of the resulting vibronic dynamics poses significant challenges in establishing a quantitative connection between spectroscopic data and underlying microscopic models. Here we show how to address this challenge using numerically exact simulation methods by considering two model systems, namely the water-soluble chlorophyll-binding protein of cauliflower and the special pair of bacterial reaction centers. We demonstrate that the inclusion of the full multi-mode vibronic dynamics in numerical calculations of linear spectra leads to systematic and quantitatively significant corrections to electronic parameter estimation. These multi-mode vibronic effects are shown to be relevant in the longstanding discussion regarding the origin of long-lived oscillations in multidimensional nonlinear spectra.


title: "Smashing more than two": Deuteron production in relativistic heavy ion collisions via stochastic multi-particle reactions

abstract: We study the deuteron production via the deuteron pion and nucleon catalysis reactions, $\pi p n \leftrightarrow \pi d$ and $N p n \leftrightarrow N d$, by employing stochastic multi-particle reactions in the hadronic transport approach SMASH for the first time. This is an improvement compared to previous studies, which introduced an artificial fake resonance $d'$ to simulate these $3 \leftrightarrow 2$ reactions as a chain of $2\leftrightarrow 2$ reactions. The derivation of the stochastic criterion for multi-particle reactions is presented in a comprehensive fashion and its implementation is tested against an analytic expression for the scattering rate and the equilibrating particle yields in box calculations. We then study Au + Au collisions at $\sqrt{s_{\mathrm{NN}}} = 7.7$ GeV, where we find that multi-particle collisions substantially reduce the time required for deuterons to reach partial chemical equilibrium with nucleons. Subsequently, the final yield of $d$ is practically independent from the number of $d$ at particlization, confirming the results of previous studies. The mean transverse momentum and the integrated elliptic flow as a function of centrality are rather insensitive to the exact realization of the $2\leftrightarrow 3$ reactions.


title: Outage Performance Analysis of Widely Linear Receivers in Uplink Multi-user MIMO Systems

abstract: This paper considers the application of widely linear (WL) receivers in an uplink multi-user system using real-valued modulation schemes, where the cellular base station (BS) with multiple antennas provides connectivity for randomly deployed single-antenna users. The targeted use case is massive machine type communication (mMTC) with grant-free access in the uplink, where the network is required to host a large number of low data rate devices transmitting in an uncoordinated fashion. Four types of WL receivers are investigated, namely the WL zero-forcing (ZF) and the WL minimum meansquared error (MMSE) receivers, along with their enhanced versions employing successive interference cancellation (SIC) with channel-dependent ordering, i.e., the WL-ZF-SIC and WL-MMSE-SIC receivers. The outage performances of these receivers are analytically characterized in the high signal-to-noise ratio (SNR) regime and compared to those of conventional linear (CL) receivers using complex-valued modulation schemes. For the non-SIC receivers, we show that, when compared to the CL counterparts, the WL receivers yield a higher diversity gain when decoding the same number of users and have the same diversity gain but a decreased coding gain when the number of users is nearly doubled. The outage performance analysis of WL-SIC receivers is facilitated by the marginal distribution of ordered eigenvalues of a real-valued Wishart matrix. It is shown that the SIC operation with channel-dependent ordering brings no additional diversity gain to the WL receivers but instead increases the coding gain. Moreover, the coding gain of WL-SIC receivers grows as the number of users increases and even exceeds that of CL-SIC receivers under suitable conditions.


title: Global Convergence of Gradient Descent for Asymmetric Low-Rank Matrix Factorization

abstract: We study the asymmetric low-rank factorization problem: \[\min_{\mathbf{U} \in \mathbb{R}^{m \times d}, \mathbf{V} \in \mathbb{R}^{n \times d}} \frac{1}{2}\|\mathbf{U}\mathbf{V}^\top -\mathbf{\Sigma}\|_F^2\] where $\mathbf{\Sigma}$ is a given matrix of size $m \times n$ and rank $d$. This is a canonical problem that admits two difficulties in optimization: 1) non-convexity and 2) non-smoothness (due to unbalancedness of $\mathbf{U}$ and $\mathbf{V}$). This is also a prototype for more complex problems such as asymmetric matrix sensing and matrix completion. Despite being non-convex and non-smooth, it has been observed empirically that the randomly initialized gradient descent algorithm can solve this problem in polynomial time. Existing theories to explain this phenomenon all require artificial modifications of the algorithm, such as adding noise in each iteration and adding a balancing regularizer to balance the $\mathbf{U}$ and $\mathbf{V}$. This paper presents the first proof that shows randomly initialized gradient descent converges to a global minimum of the asymmetric low-rank factorization problem with a polynomial rate. For the proof, we develop 1) a new symmetrization technique to capture the magnitudes of the symmetry and asymmetry, and 2) a quantitative perturbation analysis to approximate matrix derivatives. We believe both are useful for other related non-convex problems.


title: Darker than Black-Box: Face Reconstruction from Similarity Queries

abstract: Several methods for inversion of face recognition models were recently presented, attempting to reconstruct a face from deep templates. Although some of these approaches work in a black-box setup using only face embeddings, usually, on the end-user side, only similarity scores are provided. Therefore, these algorithms are inapplicable in such scenarios. We propose a novel approach that allows reconstructing the face querying only similarity scores of the black-box model. While our algorithm operates in a more general setup, experiments show that it is query efficient and outperforms the existing methods.


title: A new low-beta regime for unstable proton firehose modes in bi-Kappa distributed plasmas

abstract: In the solar wind plasma an excess of kinetic temperature along the background magnetic field stimulates proton firehose modes to grow if the parallel plasma beta parameter is sufficiently high, i.e., $\beta_{p \parallel}\gtrsim 1$. This instability can prevent the expansion-driven anisotropy from increasing indefinitely, and explain the observations. Moreover, such kinetic instabilities are expected to be even more effective in the presence of suprathermal Kappa-distributed populations, which are ubiquitous in the solar wind, are less affected by collisions than the core population, but contribute with an additional free energy. In this work we use both linear and extended quasi-linear (QL) frameworks to characterize the unstable periodic proton firehose modes (propagating parallel to the magnetic field) under the influence of suprathermal protons. Linear theory predicts a systematic stimulation of the instability, suprathermals amplifying the growth rates and decreasing the instability thresholds to lower anisotropies and lower plasma betas ($\beta_{p \parallel}<1$). In perfect agreement with these results, the QL approach reveals a significant enhancement of the resulting electromagnetic fluctuations up to the saturation with a stronger back reaction on protons, leading also to a faster and more efficient relaxation of the temperature anisotropy.


title: Knee Osteoarthritis Severity Prediction using an Attentive Multi-Scale Deep Convolutional Neural Network

abstract: Knee Osteoarthritis (OA) is a destructive joint disease identified by joint stiffness, pain, and functional disability concerning millions of lives across the globe. It is generally assessed by evaluating physical symptoms, medical history, and other joint screening tests like radiographs, Magnetic Resonance Imaging (MRI), and Computed Tomography (CT) scans. Unfortunately, the conventional methods are very subjective, which forms a barrier in detecting the disease progression at an early stage. This paper presents a deep learning-based framework, namely OsteoHRNet, that automatically assesses the Knee OA severity in terms of Kellgren and Lawrence (KL) grade classification from X-rays. As a primary novelty, the proposed approach is built upon one of the most recent deep models, called the High-Resolution Network (HRNet), to capture the multi-scale features of knee X-rays. In addition, we have also incorporated an attention mechanism to filter out the counterproductive features and boost the performance further. Our proposed model has achieved the best multiclass accuracy of 71.74% and MAE of 0.311 on the baseline cohort of the OAI dataset, which is a remarkable gain over the existing best-published works. We have also employed the Gradient-based Class Activation Maps (Grad-CAMs) visualization to justify the proposed network learning.


title: On determining the homological Conley index of Poincar\'e maps in autonomous systems

abstract: A theorem on computation of the homological Conley index of an isolated invariant set of the Poincar\'e map associated to a section in a rotating local dynamical system $\phi$ is proved. Let $(N,L)$ be an index pair for a discretization $\phi^h$ of $\phi$, where $h>0$, and let $S$ denote the invariant part of $N\setminus L$; it follows that the section $S_0$ of $S$ is an isolated invariant set of the Poincar\'e map. The theorem asserts that if the sections $N_0$ of $N$ and $L_0$ of $L$ are ANRs, the homology classes $[u_j]$ of some cycles $u_j$ form a basis of $H(N_0,L_0)$, and for some scalars $a_{ij}$, the cycles $u_j$ and $\sum a_{ij}u_i$ are homologous in the covering pair $(\widetilde N,\widetilde L)$ of $(N,L)$ and the homology relation is preserved in $(\widetilde N,\widetilde L)$ under the transformation induced by $\phi^t$ for $t\in [0,h]$ then the homological Conley index of $S_0$ is equal to the Leray reduction of the matrix $[a_{ij}]$. In particular, no information on the values of the Poincar\'e map or its approximations is required. In a special case of the system generated by a $T$-periodic non-autonomous ordinary differential equation with rational $T/h>1$, the theorem was proved in the paper M.\,Mrozek, R.\,Srzednicki, and F.\,Weilandt, SIAM J. Appl. Dyn. Syst. 14 (2015), 1348-1386, and it motivated a construction of an algorithm for determining the index.


title: Social Inequality Analysis of Fiber Bundle Model Statistics and Prediction of Materials Failure

abstract: Inequalities are abundant in a society with a number of agents competing for a limited amount of resource. Statistics of such social inequalities are usually represented by the Lorenz function $L(p)$, where $p$ fraction of the population possesses $L(p)$ fraction of the total wealth, when the population is arranged in the ascending order of their wealth. Similarly, in scientometrics, such inequalities can be represented by a plot of the citation count against the respective number of papers by a scientist, again arranged in the ascending order of their citation counts. Quantitatively, these inequalities are captured by the corresponding inequality indices, namely the Kolkata $k$ and the Hirsch $h$ indices, given by the fixed points of these nonlinear functions. In statistical physics of criticality, the fixed points of the Renormalization Group generator functions are studied in their self-similar limit, where its (fractal) structure converges to a unique form. The statistical indices in the social science, however, correspond to the fixed points where the values of the generator function (wealth or citation sizes) are commensurately abundant in fractions or numbers (of persons or papers). We introduce and study these indices for the inequalities of (pre-failure) avalanches, given by their size distributions in the Fiber Bundle Models (FBM) of non-brittle materials. We show how a prior knowledge of the terminal and almost universal value of the $k$ index for a wide range of disorder parameter, can help in predicting an imminent catastrophic breakdown in the model. This observation has also been complemented by noting a similar (but not identical) behavior of the Hirsch index ($h$), redefined for such avalanche statistics.


title: A Reinforcement Learning Approach for Sequential Spatial Transformer Networks

abstract: Spatial Transformer Networks (STN) can generate geometric transformations which modify input images to improve the classifier's performance. In this work, we combine the idea of STN with Reinforcement Learning (RL). To this end, we break the affine transformation down into a sequence of simple and discrete transformations. We formulate the task as a Markovian Decision Process (MDP) and use RL to solve this sequential decision-making problem. STN architectures learn the transformation parameters by minimizing the classification error and backpropagating the gradients through a sub-differentiable sampling module. In our method, we are not bound to the differentiability of the sampling modules. Moreover, we have freedom in designing the objective rather than only minimizing the error; e.g., we can directly set the target as maximizing the accuracy. We design multiple experiments to verify the effectiveness of our method using cluttered MNIST and Fashion-MNIST datasets and show that our method outperforms STN with a proper definition of MDP components.


title: Ferromagnetic Cr4PtGa17: A Novel Half-Heusler-Type Compound with a Breathing Pyrochlore Lattice

abstract: We describe the crystal structure and elementary magnetic properties of a previously unreported ternary intermetallic compound, Cr4PtGa17, which crystallizes in a rhombohedral unit cell in the noncentrosymmetric space group R3m. The crystal structure is closely related to those of XYZ half-Heusler compounds, where X, Y and Z are reported to be single elements only, occupying three different face-centered cubic sublattices. The new material, Cr4PtGa17, can be most straightforwardly illustrated by writing the formula as (PtGa2)(Cr4Ga14)Ga (X=PtGa2, Y = Cr4Ga14, Z = Ga), that is, the X and Y sites are occupied by clusters instead of single elements. The magnetic Cr occupies a breathing pyrochlore lattice. Ferromagnetic ordering is found below TC ~61 K, by both neutron diffraction and magnetometer studies, with a small, saturated moment of ~0.25 muB/Cr observed at 2 K, making Cr4PtGa17 the first ferromagnetically ordered material with a breathing pyrochlore lattice. A magnetoresistance of ~140% was observed at 2 K. DFT calculations suggest that the material has a nearly-half-metallic electronic structure. The new material, Cr4PtGa17, the first realization of both a half-Heusler-type structure and a breathing pyrochlore lattice, might pave a new way to achieve novel types of half-Heusler compounds.


title: Hadron gas in the presence of a magnetic field using non-extensive statistics: A transition from diamagnetic to paramagnetic system

abstract: Non-central heavy-ion collisions at ultra-relativistic energies are unique in producing magnetic fields of the largest strength in the laboratory. Such fields being produced at the early stages of the collision could affect the properties of Quantum Chromodynamics (QCD) matter formed in the relativistic heavy-ion collisions. The transient magnetic field leaves its reminiscence, which in principle, can affect the thermodynamic and transport properties of the final state dynamics of the system. In this work, we study the thermodynamic properties of a hadron gas in the presence of an external static magnetic field using a thermodynamically consistent non-extensive Tsallis distribution function. Various thermodynamical observables such as energy density ($\epsilon$), entropy density ($s$), pressure ($P$) and speed of sound ($c_{\rm s}$) are studied. Investigation of magnetization ($M$) is also performed and this analysis reveals an interplay of diamagnetic and paramagnetic nature of the system in the presence of a magnetic field of varying strength. Further, to understand the system dynamics under equilibrium and non-equilibrium conditions, the effect of the non-extensive parameter ($q$) on the above observables is also studied.


title: A2A: 21,000 bulge stars from the ARGOS survey with stellar parameters on the APOGEE scale

abstract: We use the data-driven method, The Cannon, to bring 21,000 stars from the ARGOS bulge survey, including 10,000 red clump stars, onto the parameter and abundance scales of the cross-Galactic survey, APOGEE, obtaining rms precisions of 0.10 dex, 0.07 dex, 74 K, and 0.18 dex for [Fe/H], [Mg/Fe], Teff, and log(g), respectively. The re-calibrated ARGOS survey - which we refer to as the A2A survey - is combined with the APOGEE survey to investigate the abundance structure of the Galactic bulge. We find X-shaped [Fe/H] and [Mg/Fe] distributions in the bulge that are more pinched than the bulge density, a signature of its disk origin. The mean abundance along the major axis of the bar varies such that the stars are more [Fe/H]-poor and [Mg/Fe]-rich near the Galactic center than in the long bar/outer bulge region. The vertical [Fe/H] and [Mg/Fe] gradients vary between the inner bulge and long bar with the inner bulge showing a flattening near the plane that is absent in the long bar. The [Fe/H]-[Mg/Fe] distribution shows two main maxima, an ``[Fe/H]-poor [Mg/Fe]- rich'' maximum and an ``[Fe/H]-rich [Mg/Fe]-poor'' maximum, that vary in strength with position in the bulge. In particular, the outer long bar close to the Galactic plane is dominated by super-solar [Fe/H], [Mg/Fe]-normal stars. Stars composing the [Fe/H]-rich maximum show little kinematic dependence on [Fe/H], but for lower [Fe/H] the rotation and dispersion of the bulge increase slowly. Stars with [Fe/H]<-1 dex have a very different kinematic structure than stars with higher [Fe/H]. Comparing with recent models for the Galactic boxy-peanut bulge, the abundance gradients and distribution, and the relation between [Fe/H] and kinematics suggest that the stars comprising each maximum have separate disk origins with the ``[Fe/H]-poor [Mg/Fe]-rich'' stars originating from a thicker disk than the ``[Fe/H]-rich [Mg/Fe]-poor'' stars.


title: Quantum and classical cosmology in the Brans-Dicke theory

abstract: In this paper we discuss classical and quantum aspects of cosmological models in Brans-Dicke theory. First, we review cosmological bounce solution in Brans-Dicke theory that obeys energy conditions (without ghost) for a universe filled with radiative fluid. Then we quantize this classical model in a canonical way, establishing the corresponding Wheeler-DeWitt equation in the minisuperspace, and analyze the quantum solutions. When the energy conditions are violated, corresponding to the case $\omega<-3/2$, the energy is bounded from below and singularity-free solutions are found. However, in the case $\omega> -3/2$ we cannot compute the evolution of the scale factor by evaluating the expectation values because the wave function is not finite (energy spectrum is not bounded from below). But we can analyze this case using Bohmian mechanics and the de Broglie-Bohm interpretation of quantum mechanics. Using this approach, the classical and quantum results can be compared for any value of $\omega$.


title: ASK: Adversarial Soft k-Nearest Neighbor Attack and Defense

abstract: K-Nearest Neighbor (kNN)-based deep learning methods have been applied to many applications due to their simplicity and geometric interpretability. However, the robustness of kNN-based classification models has not been thoroughly explored and kNN attack strategies are underdeveloped. In this paper, we propose an Adversarial Soft kNN (ASK) loss to both design more effective kNN attack strategies and to develop better defenses against them. Our ASK loss approach has two advantages. First, ASK loss can better approximate the kNN's probability of classification error than objectives proposed in previous works. Second, the ASK loss is interpretable: it preserves the mutual information between the perturbed input and the in-class-reference data. We use the ASK loss to generate a novel attack method called the ASK-Attack (ASK-Atk), which shows superior attack efficiency and accuracy degradation relative to previous kNN attacks. Based on the ASK-Atk, we then derive an ASK-\underline{Def}ense (ASK-Def) method that optimizes the worst-case training loss induced by ASK-Atk. Experiments on CIFAR-10 (ImageNet) show that (i) ASK-Atk achieves $\geq 13\%$ ($\geq 13\%$) improvement in attack success rate over previous kNN attacks, and (ii) ASK-Def outperforms the conventional adversarial training method by $\geq 6.9\%$ ($\geq 3.5\%$) in terms of robustness improvement.


title: How one can repair non-integrable Kahan discretizations. II. A planar system with invariant curves of degree 6

abstract: We find a novel one-parameter family of integrable quadratic Cremona maps of the plane preserving a pencil of curves of degree 6 and of genus 1. They turn out to serve as Kahan-type discretizations of a novel family of quadratic vector fields possessing a polynomial integral of degree 6 whose level curves are of genus 1, as well. These vector fields are non-homogeneous generalizations of reduced Nahm systems for magnetic monopoles with icosahedral symmetry, introduced by Hitchin, Manton and Murray. The straightforward Kahan discretization of these novel non-homogeneous systems is non-integrable. However, this drawback is repaired by introducing adjustments of order $O(\epsilon^2)$ in the coefficients of the discretization, where $\epsilon$ is the stepsize.


title: Non-destructive interrogation of nuclear waste barrels through muon tomography: A Monte Carlo study based on dual-parameter analysis via GEANT4 simulations

abstract: In this study, we employ the Monte Carlo simulations by using the GEANT4 code to demonstrate the capability of muon tomography based on the dual-parameter analysis in the examination of the nuclear waste barrels. Our current hodoscope setup consists of three top and three bottom plastic scintillators made of polyvinyl toluene with the thickness of 0.4 cm, and the composite target material is a cylindrical nuclear waste drum with the height of 96 cm and the radius of 29.6 cm where the outermost layer is stainless steel with the lateral thickness of 3.2 cm and the filling material is ordinary concrete that encapsulates the nuclear materials of dimensions 20$\times$20$\times$20 cm$^{3}$. By simulating with a narrow planar muon beam of 1$\times$1 cm$^{2}$ over the uniform energy interval between 0.1 and 8 GeV, we determine the variation of the average scattering angle together with the standard deviation by utilizing a 0.5-GeV bin length, the counts of the scattering angle by using a 1-mrad step, and the number of the absorption events for the five prevalent nuclear materials starting from cobalt and ending in plutonium. Via the duo-parametric analysis that is founded on the scattering angle as well as the absorption in the present study, we show that the presence of the nuclear materials in the waste barrels is numerically visible in comparison with the concrete-filled waste drum without any nuclear material, and the muon tomography is capable of distinguishing these nuclear materials by coupling the information about the scattering angle and the number of absorption in the cases where one of these two parameters yields strong similarity for certain nuclear materials.


title: A measurement of the scale of homogeneity in the Early Universe

abstract: We present the first measurement of the homogeneity index, $\mathcal{H}$, a fractal or Hausdorff dimension of the early Universe from the Planck CMB temperature variations $\delta T$ in the sky. This characterization of the isotropy scale is model-free and purely geometrical, independent of the amplitude of $\delta T$. We find evidence of homogeneity ($\mathcal{H}=0$) for scales larger than $\theta_{\mathcal{H}} = 65.9 \pm 9.2 \deg $ on the CMB sky. This finding is at odds with the $\Lambda$CDM prediction, which assumes a scale invariant infinite universe. Such anomaly is consistent with the well known low quadrupule amplitude in the angular $\delta T$ spectrum, but quantified in a direct and model independent way. We estimate the significance of our finding for $\mathcal{H}=0$ using a principal component analysis from the sampling variations of the observed sky. This analysis is validated with theoretical prediction of the covariance matrix \textcolor{black}{ and simulations, booth base purely on data or in the $\Lambda$CDM prediction.} Assuming translation invariance (and flat geometry) we can convert the isotropy scale $\theta_\mathcal{H}$ into a (comoving) homogeneity scale which is very close to the trapped surface generated by the observed cosmological constant $\Lambda$.


title: The consequences of hesitation: Axelrod model with intrinsic noise

abstract: We study the influence of hesitating agents in the Axelrod model by introducing an intrinsic noise, which is proportional to the disagreement between the interacting agents, and thus coupled to the dynamics. Our results show that, unlike the effect of the cultural drift where the noise is just controlled by an external parameter, the intrinsic noise never allows the system to reach a frozen state in the thermodynamic limit. Moreover, we show that when the intrinsic noise affects the agents' variables, the system's behaviour is also different from the case when it affects the network of their interactions.


title: Unsupervised Skill Discovery with Bottleneck Option Learning

abstract: Having the ability to acquire inherent skills from environments without any external rewards or supervision like humans is an important problem. We propose a novel unsupervised skill discovery method named Information Bottleneck Option Learning (IBOL). On top of the linearization of environments that promotes more various and distant state transitions, IBOL enables the discovery of diverse skills. It provides the abstraction of the skills learned with the information bottleneck framework for the options with improved stability and encouraged disentanglement. We empirically demonstrate that IBOL outperforms multiple state-of-the-art unsupervised skill discovery methods on the information-theoretic evaluations and downstream tasks in MuJoCo environments, including Ant, HalfCheetah, Hopper and D'Kitty.


title: 3D Reconstruction through Fusion of Cross-View Images

abstract: 3D recovery from multi-stereo and stereo images, as an important application of the image-based perspective geometry, serves many applications in computer vision, remote sensing and Geomatics. In this chapter, the authors utilize the imaging geometry and present approaches that perform 3D reconstruction from cross-view images that are drastically different in their viewpoints. We introduce our framework that takes ground-view images and satellite images for full 3D recovery, which includes necessary methods in satellite and ground-based point cloud generation from images, 3D data co-registration, fusion and mesh generation. We demonstrate our proposed framework on a dataset consisting of twelve satellite images and 150k video frames acquired through a vehicle-mounted Go-pro camera and demonstrate the reconstruction results. We have also compared our results with results generated from an intuitive processing pipeline that involves typical geo-registration and meshing methods.


title: Geometric Processing for Image-based 3D Object Modeling

abstract: Image-based 3D object modeling refers to the process of converting raw optical images to 3D digital representations of the objects. Very often, such models are desired to be dimensionally true, semantically labeled with photorealistic appearance (reality-based modeling). Laser scanning was deemed as the standard (and direct) way to obtaining highly accurate 3D measurements of objects, while one would have to abide the high acquisition cost and its unavailability on some of the platforms. Nowadays the image-based methods backboned by the recently developed advanced dense image matching algorithms and geo-referencing paradigms, are becoming the dominant approaches, due to its high flexibility, availability and low cost. The largely automated geometric processing of images in a 3D object reconstruction workflow, from ordered/unordered raw imagery to textured meshes, is becoming a critical part of the reality-based 3D modeling. This article summarizes the overall geometric processing workflow, with focuses on introducing the state-of-the-art methods of three major components of geometric processing: 1) geo-referencing; 2) Image dense matching 3) texture mapping. Finally, we will draw conclusions and share our outlooks of the topics discussed in this article.


title: Concentration of Contractive Stochastic Approximation and Reinforcement Learning

abstract: Using a martingale concentration inequality, concentration bounds `from time $n_0$ on' are derived for stochastic approximation algorithms with contractive maps and both martingale difference and Markov noises. These are applied to reinforcement learning algorithms, in particular to asynchronous Q-learning and TD(0).


title: Change Detection for Geodatabase Updating

abstract: The geodatabase (vectorized data) nowadays becomes a rather standard digital city infrastructure; however, updating geodatabase efficiently and economically remains a fundamental and practical issue in the geospatial industry. The cost of building a geodatabase is extremely high and labor intensive, and very often the maps we use have several months and even years of latency. One solution is to develop more automated methods for (vectorized) geospatial data generation, which has been proven a difficult task in the past decades. An alternative solution is to first detect the differences between the new data and the existing geospatial data, and then only update the area identified as changes. The second approach is becoming more favored due to its high practicality and flexibility. A highly relevant technique is change detection. This article aims to provide an overview the state-of-the-art change detection methods in the field of Remote Sensing and Geomatics to support the task of updating geodatabases. Data used for change detection are highly disparate, we therefore structure our review intuitively based on the dimension of the data, being 1) change detection with 2D data; 2) change detection with 3D data. Conclusions will be drawn based on the reviewed efforts in the field, and we will share our outlooks of the topic of updating geodatabases.


title: Optimal Control of Closed Quantum Systems via B-Splines with Carrier Waves

abstract: We consider the optimal control problem of determining electromagnetic pulses for implementing logical gates in a closed quantum system, where the Hamiltonian models the dynamics of coupled superconducting qudits. The quantum state is governed by Schr\"odinger's equation, which we formulate in terms of the real and imaginary parts of the state vector and solve by the St\"ormer-Verlet scheme, which is a symplectic partitioned Runge-Kutta method. A novel parameterization of the control functions based on B-splines with carrier waves is introduced. The carrier waves are used to trigger the resonant frequencies in the system Hamiltonian, and the B-spline functions specify their amplitude and phase. This approach allows the number of control parameters to be independent of, and significantly smaller than, the number of time steps for integrating Schr\"odinger's equation. We present numerical examples of how the proposed technique can be combined with an interior point L-BFGS algorithm for realizing quantum gates, and generalize our approach to calculate risk-neutral controls that are resilient to noise in the Hamiltonian model. The proposed method is also shown to compare favorably with QuTiP/pulse\_optim and Grape-Tensorflow.


title: Using low dose X-ray Speckle Visibility Spectroscopy to study dynamics of soft matter samples

abstract: We demonstrate the successful application of X-ray Speckle Visibility Spectroscopy (XSVS) experiments to study the dynamics of radiation sensitive, biological samples with unprecedentedly small X-ray doses of 45 Gy and below. Using XSVS, we track the dynamics of casein micelles in native, concentrated, and acidified solution conditions, while substantially reducing the deposited dose as compared to alternative techniques like sequential X-ray photon correlation spectroscopy (XPCS). The Brownian motion in a skim milk sample yields the hydrodynamic radius of the casein micelles while deviations from Brownian motion with a characteristic $q$-dependent diffusion coefficient $D(q)$ can be observed in more concentrated solution conditions. The low dose applied in our experiments allows the observation of static, frozen speckle patterns from gelled acidic milk. We show that the XSVS technique is especially suitable for tracking dynamics of radiation sensitive samples in combination with the improved coherent properties of new generation X-ray sources, emphasizing the great potential for further investigations of protein dynamics using fourth generation synchrotrons and free electron lasers.


title: An Analysis of Multi-hop Iterative Approximate Byzantine Consensus with Local Communication

abstract: Iterative Approximate Byzantine Consensus (IABC) is a fundamental problem of fault-tolerant distributed computing where machines seek to achieve approximate consensus to arbitrary exactness in the presence of Byzantine failures. We present a novel algorithm for this problem, named Relay-IABC, which relies on the usage of a multi-hop relayed messaging system and crytographically secure message signatures. The use of signatures and relays allows the strict necessary network conditions of traditional IABC algorithms to be circumvented. In addition, we show evidence that Relay-IABC achieves faster convergence than traditional algorithms even under these strict network conditions with both theoretical analysis and experimental results.


title: AniVis: Generating Animated Transitions Between Statistical Charts with a Tree Model

abstract: Animated transitions help viewers understand changes between related visualizations. To clearly present the underlying relations between statistical charts, animation authors need to have a high level of expertise and a considerable amount of time to describe the relations with reasonable animation stages. We present AniVis, an automated approach for generating animated transitions to demonstrate the changes between two statistical charts. AniVis models each statistical chart into a tree-based structure. Given an input chart pair, the differences of data and visual properties of the chart pair are formalized as tree edit operations. The edit operations can be mapped to atomic transition units. Through this approach, the animated transition between two charts can be expressed as a set of transition units. Then, we conduct a formative study to understand people's preferences for animation sequences. Based on the study, we propose a set of principles and a sequence composition algorithm to compose the transition units into a meaningful animation sequence. Finally, we synthesize these units together to deliver a smooth and intuitive animated transition between charts. To test our approach, we present a prototype system and its generated results to illustrate the usage of our framework. We perform a comparative study to assess the transition sequence derived from the tree model. We further collect qualitative feedback to evaluate the effectiveness and usefulness of our method.


title: Truncated Metric Dimension for Finite Graphs

abstract: A graph $G=(V,E)$ with geodesic distance $d(\cdot,\cdot)$ is said to be resolved by a non-empty subset $R$ of its vertices when, for all vertices $u$ and $v$, if $d(u,r)=d(v,r)$ for each $r\in R$, then $u=v$. The metric dimension of $G$ is the cardinality of its smallest resolving set. In this manuscript, we present and investigate the notions of resolvability and metric dimension when the geodesic distance is truncated with a certain threshold $k$; namely, we measure distances in $G$ using the metric $d_k(u,v):=\min\{d(u,v),k+1\}$. We denote the metric dimension of $G$ with respect to $d_k$ as $\beta_k(G)$. We study the behavior of this quantity with respect to $k$ as well as the diameter of $G$. We also characterize the truncated metric dimension of paths and cycles as well as graphs with extreme metric dimension, including graphs of order $n$ such that $\beta_k(G)=n-2$ and $\beta_k(G)=n-1$. We conclude with a study of various problems related to the truncated metric dimension of trees.


title: Point to Point Ethernet TransmissionWireless Backhaul Links Clustering

abstract: Arxiv is acting weird and throwing error: "Bad character(s) in field Abstract." for no reason. Please refer to the manuscript.


title: PYInfer: Deep Learning Semantic Type Inference for Python Variables

abstract: Python type inference is challenging in practice. Due to its dynamic properties and extensive dependencies on third-party libraries without type annotations, the performance of traditional static analysis techniques is limited. Although semantics in source code can help manifest intended usage for variables (thus help infer types), they are usually ignored by existing tools. In this paper, we propose PYInfer, an end-to-end learning-based type inference tool that automatically generates type annotations for Python variables. The key insight is that contextual code semantics is critical in inferring the type for a variable. For each use of a variable, we collect a few tokens within its contextual scope, and design a neural network to predict its type. One challenge is that it is difficult to collect a high-quality human-labeled training dataset for this purpose. To address this issue, we apply an existing static analyzer to generate the ground truth for variables in source code. Our main contribution is a novel approach to statically infer variable types effectively and efficiently. Formulating the type inference as a classification problem, we can handle user-defined types and predict type probabilities for each variable. Our model achieves 91.2% accuracy on classifying 11 basic types in Python and 81.2% accuracy on classifying 500 most common types. Our results substantially outperform the state-of-the-art type annotators. Moreover, PYInfer achieves 5.2X more code coverage and is 187X faster than a state-of-the-art learning-based tool. With similar time consumption, our model annotates 5X more variables than a state-of-the-art static analysis tool. Our model also outperforms a learning-based function-level annotator on annotating types for variables and function arguments. All our tools and datasets are publicly available to facilitate future research in this direction.


title: Necessary and sufficient conditions for regularity of interval parametric matrices

abstract: Matrix regularity is a key to various problems in applied mathematics. The sufficient conditions, used for checking regularity of interval parametric matrices, usually fail in case of large parameter intervals. We present necessary and sufficient conditions for regularity of interval parametric matrices in terms of boundary parametric hypersurfaces, parametric solution sets, determinants, real spectral radiuses. The initial n-dimensional problem involving K interval parameters is replaced by numerous problems involving 1<= t <= min(n-1, K) interval parameters, in particular t=1 is most attractive. The advantages of the proposed methodology are discussed along with its application for finding the interval hull solution to interval parametric linear system and for determining the regularity radius of an interval parametric matrix.


title: Effects of water currents on fish migration through a Feynman-type path integral approach under $\sqrt{8/3}$ Liouville-like quantum gravity surfaces

abstract: A stochastic differential game theoretic model has been proposed to determine optimal behavior of a fish while migrating against water currents both in rivers and oceans. Then, a dynamic objective function is maximized subject to two stochastic dynamics, one represents its location and another its relative velocity against water currents. In relative velocity stochastic dynamics, a Cucker-Smale type stochastic differential equation is introduced under white noise. As the information regarding hydrodynamic environment is incomplete and imperfect, a Feynman type path integral under $\sqrt{8/3}$ Liouville-like quantum gravity surface has been introduced to obtain a Wick-rotated Schr\"odinger type equation to determine an optimal strategy of a fish during its migration. The advantage of having Feynman type path integral is that, it can be used in more generalized nonlinear stochastic differential equations where constructing a Hamiltonian-Jacobi-Bellman (HJB) equation is impossible. The mathematical analytic results show exact expression of an optimal strategy of a fish under imperfect information and uncertainty.


title: A Fast Response Mission to Rendezvous with an Interstellar Object

abstract: A solar sail propelled small satellite mission concept to intercept and potentially rendezvous with newly discovered transient interstellar objects (ISOs) is described. The mission concept derives from the proposal for a technology demonstration mission for exiting the solar system at high velocity, eventually to reach the focal region of the solar gravitational lens. The ISO mission concept is to fly a solar sail toward a holding orbit around the Sun and when the ISO orbit is confirmed, target the sailcraft to reach an escape velocity of over 6 AU/year. This would permit rapid response to a new ISO discovery and an intercept within 10 AU from the Sun. Two new proven interplanetary technologies are utilized to enable such a mission: i) interplanetary smallsats, such as those demonstrated by the MarCO mission, and ii) solar sails, such as demonstrated by LightSail and IKAROS missions and developed for NEA Scout and Solar Cruiser missions. Current technology work suggests that already within this decade such a mission could fly and reach an ISO moving through the solar system. It might enable the first encounter with an ISO to allow for imaging and spectroscopy, measurements of size and mass, potentially giving a unique information about the object's origin and composition. A similar approach may be used to allow for a sample return.


title: Legendre Deep Neural Network (LDNN) and its application for approximation of nonlinear Volterra Fredholm Hammerstein integral equations

abstract: Various phenomena in biology, physics, and engineering are modeled by differential equations. These differential equations including partial differential equations and ordinary differential equations can be converted and represented as integral equations. In particular, Volterra Fredholm Hammerstein integral equations are the main type of these integral equations and researchers are interested in investigating and solving these equations. In this paper, we propose Legendre Deep Neural Network (LDNN) for solving nonlinear Volterra Fredholm Hammerstein integral equations (VFHIEs). LDNN utilizes Legendre orthogonal polynomials as activation functions of the Deep structure. We present how LDNN can be used to solve nonlinear VFHIEs. We show using the Gaussian quadrature collocation method in combination with LDNN results in a novel numerical solution for nonlinear VFHIEs. Several examples are given to verify the performance and accuracy of LDNN.


title: Draw Me a Flower: Processing and Grounding Abstraction in Natural Language

abstract: Abstraction is a core tenet of human cognition and communication. When composing natural language instructions, humans naturally evoke abstraction to convey complex procedures in an efficient and concise way. Yet, interpreting and grounding abstraction expressed in NL has not yet been systematically studied in NLP, with no accepted benchmarks specifically eliciting abstraction in NL. In this work, we set the foundation for a systematic study of processing and grounding abstraction in NLP. First, we deliver a novel abstraction elicitation method and present Hexagons, a 2D instruction-following game. Using Hexagons we collected over 4k naturally-occurring visually-grounded instructions rich with diverse types of abstractions. From these data, we derive an instruction-to-execution task and assess different types of neural models. Our results show that contemporary models and modeling practices are substantially inferior to human performance, and that models' performance is inversely correlated with the level of abstraction, showing less satisfying performance on higher levels of abstraction. These findings are consistent across models and setups, confirming that abstraction is a challenging phenomenon deserving further attention and study in NLP/AI research.


title: The Transient Ultra-luminous X-ray Source, ULX-4, in M51

abstract: We present the results of a temporal and spectral analysis of the transient source ULX-4 in the galaxy M51. The data used were drawn from Chandra, XMM-Newton and Swift-XRT archives, spanning the years 2000-2019. The X-ray flux of the source is seen to vary by two orders of magnitudes within a month but a short-term variability was not observed over the time intervals of 100-2000 second in the 0.3-10 keV energy band. We find some evidence for the existence of bi-modality feature in the flux distribution of ULX-4. We identified two optical sources as possible counterparts within an error radius of 0."18 at 95% confidence level for ULX-4 based on the archival HST/ACS and HST/WFC3 data. Blackbody fits of the spectral energy distributions (SEDs) indicate the spectral type to be B-type stars. One of these counterparts exhibits a low-amplitude optical periodicity of 264 days in the F606W filter; if we assume this apparent periodicity is associated with the orbital motion of the donor, then it is more likely that the donor is a red supergiant satisfying the long periodicity and accretion via Roche-lobe overflow. Consequently, the SED would then have to be interpreted as a superposition of emissions from a cold donor and a hot flow component, most likely from an accretion disk. If, on the other hand, the periodicity is super orbital in nature i.e., due to possible interactions of the compact object with a circumstellar disk, the donor could then be a Be/X star hosting a neutron star.


title: Use of Variational Inference in Music Emotion Recognition

abstract: This work was developed aiming to employ Statistical techniques to the field of Music Emotion Recognition, a well-recognized area within the Signal Processing world, but hardly explored from the statistical point of view. Here, we opened several possibilities within the field, applying modern Bayesian Statistics techniques and developing efficient algorithms, focusing on the applicability of the results obtained. Although the motivation for this project was the development of a emotion-based music recommendation system, its main contribution is a highly adaptable multivariate model that can be useful interpreting any database where there is an interest in applying regularization in an efficient manner. Broadly speaking, we will explore what role a sound theoretical statistical analysis can play in the modeling of an algorithm that is able to understand a well-known database and what can be gained with this kind of approach.


title: Learning stochastic object models from medical imaging measurements by use of advanced ambient generative adversarial networks

abstract: Purpose: To objectively assess new medical imaging technologies via computer-simulations, it is important to account for the variability in the ensemble of objects to be imaged. This source of variability can be described by stochastic object models (SOMs). It is generally desirable to establish SOMs from experimental imaging measurements acquired by use of a well-characterized imaging system, but this task has remained challenging. Approach: A generative adversarial network (GAN)-based method that employs AmbientGANs with modern progressive or multiresolution training approaches is proposed. AmbientGANs established using the proposed training procedure are systematically validated in a controlled way using computer-simulated magnetic resonance imaging (MRI) data corresponding to a stylized imaging system. Emulated single-coil experimental MRI data are also employed to demonstrate the methods under less stylized conditions. Results: The proposed AmbientGAN method can generate clean images when the imaging measurements are contaminated by measurement noise. When the imaging measurement data are incomplete, the proposed AmbientGAN can reliably learn the distribution of the measurement components of the objects. Conclusions: Both visual examinations and quantitative analyses, including task-specific validations using the Hotelling observer, demonstrated that the proposed AmbientGAN method holds promise to establish realistic SOMs from imaging measurements.


title: Photonic Topological Baths for Quantum Simulation

abstract: Quantum simulation involves engineering devices to implement different Hamiltonians and measuring their quantized spectra to study quantum many-body systems. Recent developments in topological photonics have shown the possibility of studying novel quantum phenomena by controlling the topological properties of such devices. Here, using coupled arrays of upto 16 high Q nano-cavities we experimentally realize quantum photonic baths which are analogs of the Su-Schrieffer-Heeger model. We investigate the effect of fabrication induced disorder on these baths by probing individual super-modes and demonstrate the design mitigation steps required to overcome the disorder effects on the quantum phenomena.


title: Last-iterate Convergence in Extensive-Form Games

abstract: Regret-based algorithms are highly efficient at finding approximate Nash equilibria in sequential games such as poker games. However, most regret-based algorithms, including counterfactual regret minimization (CFR) and its variants, rely on iterate averaging to achieve convergence. Inspired by recent advances on last-iterate convergence of optimistic algorithms in zero-sum normal-form games, we study this phenomenon in sequential games, and provide a comprehensive study of last-iterate convergence for zero-sum extensive-form games with perfect recall (EFGs), using various optimistic regret-minimization algorithms over treeplexes. This includes algorithms using the vanilla entropy or squared Euclidean norm regularizers, as well as their dilated versions which admit more efficient implementation. In contrast to CFR, we show that all of these algorithms enjoy last-iterate convergence, with some of them even converging exponentially fast. We also provide experiments to further support our theoretical results.


title: Waves propagating parallel to the magnetic field in relativistically hot plasmas: A hydrodynamic model

abstract: The high-frequency part of spectrum of electromagnetic waves propagating parallel to the external magnetic field is considered for the macroscopically motionless plasmas with the relativistic temperatures $T\sim m_{e}c^{2}$, where $m_{e}$ is the mass of electron, $c$ is the speed of light. The analysis is based on the novel hydrodynamic model based on four equations for the material fields which can be combined in two four vectors. These material fields are the concentration and the velocity field \emph{and} the average reverse relativistic $\gamma$ functor and the flux of the reverse relativistic $\gamma$ functor. In the nonrelativistic regime we have three waves (the ions are assumed to be motionless). Strong thermal effects lead to a coefficient in front of cyclotron frequency which decreases the effective contribution of the cyclotron frequency. At $T=0.1m_{e}c^{2}$ we have a decrease of area of existence of fast magneto-sound wave from the area of the large frequencies. While the area of existence of extraordinary waves becomes larger towards smaller frequencies. The strong magnetic field limit $\mid\Omega_{e}\mid > \omega_{Le}$ additional wave appears with frequency below thermally decreased cyclotron frequency, where $\mid\Omega_{e}\mid$ is the electron cyclotron frequency, and $\omega_{Le}$ is the Langmuir frequency. Further increase of temperature leads to the disappearance of fast magneto-sound wave and to the considerable increase of area of existence of extraordinary towards smaller frequencies.