title: Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics

abstract: Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.


title: A semi-implicit compressible model for atmospheric flows with seamless access to soundproof and hydrostatic dynamics

abstract: We introduce a second-order numerical scheme for compressible atmospheric motions at small to planetary scales. The collocated finite volume method treats the advection of mass, momentum, and mass-weighted potential temperature in conservation form while relying on Exner pressure for the pressure gradient term. It discretises the rotating compressible equations by evolving full variables rather than perturbations around a background state, and operates with time steps constrained by the advection speed only. Perturbation variables are only used as auxiliary quantities in the formulation of the elliptic problem. Borrowing ideas on forward-in-time differencing, the algorithm reframes the authors' previously proposed schemes into a sequence of implicit midpoint, advection, and implicit trapezoidal steps that allows for a time integration unconstrained by the internal gravity wave speed. Compared with existing approaches, results on a range of benchmarks of nonhydrostatic- and hydrostatic-scale dynamics are competitive. The test suite includes a new planetary-scale inertia-gravity wave test highlighting the properties of the scheme and its large time step capabilities. In the hydrostatic-scale cases the model is run in pseudo-incompressible and hydrostatic mode with simple switching within a uniform discretization framework. The differences with the compressible runs return expected relative magnitudes. By providing seamless access to soundproof and hydrostatic dynamics, the developments represent a necessary step towards an all-scale blended multimodel solver.


title: An "On The Fly" Framework for Efficiently Generating Synthetic Big Data Sets

abstract: Collecting, analyzing and gaining insight from large volumes of data is now the norm in an ever increasing number of industries. Data analytics techniques, such as machine learning, are powerful tools used to analyze these large volumes of data. Synthetic data sets are routinely relied upon to train and develop such data analytics methods for several reasons: to generate larger data sets than are available, to generate diverse data sets, to preserve anonymity in data sets with sensitive information, etc. Processing, transmitting and storing data is a key issue faced when handling large data sets. This paper presents an "On the fly" framework for generating big synthetic data sets, suitable for these data analytics methods, that is both computationally efficient and applicable to a diverse set of problems. An example application of the proposed framework is presented along with a mathematical analysis of its computational efficiency, demonstrating its effectiveness.


title: Low-frequency radio study of MACS clusters at 610 and 235 MHz using the GMRT

abstract: Studies have shown that mergers of massive galaxy clusters produce shocks and turbulence in the intra-cluster medium, the possible event that creates radio relics, as well as the radio halos. Here we present GMRT dual-band (235 and 610~MHz) radio observations of four such clusters from the MAssive Cluster Survey (MACS) catalogue. We report the discovery of a very faint, diffuse, elongated radio source with a projected size of about 0.5~Mpc in cluster MACSJ0152.5-2852. We also confirm the presence of a radio relic-like source (about 0.4~Mpc, previously reported at 325~MHz) in MACSJ0025.4-1222 cluster. Proposed relics in both these clusters are found apparently inside the virial radius instead of their usual peripheral location, while no radio halos are detected. These high-redshift clusters ($z=0.584$ and $0.413$) are among the earliest merging systems detected with cluster radio emissions. In MACSJ1931-2635 cluster, we found a radio mini-halo and an interesting highly bent pair of radio jets. Further, we present here a maiden study of low frequency (GMRT $235\;\&\;610$~MHz) spectral and morphological signatures of a previously known radio cluster MACSJ0014.3-3022 (Abell~2744). This cluster hosts a relatively flat spectrum ($\alpha^{610}_{235}\sim -1.15$), giant ($\sim 1.6$~Mpc each) halo-relic structure and a close-by high-speed ($1769\pm^{148}_{359}$~km~s$^{-1}$) merger-shock ($\mathcal{M}=2.02\pm^{0.17}_{0.41}$) originated from a possible second merger in the cluster.


title: Day-Ahead Hourly Forecasting of Power Generation from Photovoltaic Plants

abstract: The ability to accurately forecast power generation from renewable sources is nowadays recognised as a fundamental skill to improve the operation of power systems. Despite the general interest of the power community in this topic, it is not always simple to compare different forecasting methodologies, and infer the impact of single components in providing accurate predictions. In this paper we extensively compare simple forecasting methodologies with more sophisticated ones over 32 photovoltaic plants of different size and technology over a whole year. Also, we try to evaluate the impact of weather conditions and weather forecasts on the prediction of PV power generation.


title: Securing Accelerators with Dynamic Information Flow Tracking

abstract: Systems-on-chip (SoCs) are becoming heterogeneous: they combine general-purpose processor cores with application-specific hardware components, also known as accelerators, to improve performance and energy efficiency. The advantages of heterogeneity, however, come at a price of threatening security. The architectural dissimilarities of processors and accelerators require revisiting the current security techniques. With this hardware demo, we show how accelerators can break dynamic information flow tracking (DIFT), a well-known security technique that protects systems against software-based attacks. We also describe how the security guarantees of DIFT can be re-established with a hardware solution that has low performance and area penalties.


title: Workflow-Driven Distributed Machine Learning in CHASE-CI: A Cognitive Hardware and Software Ecosystem Community Infrastructure

abstract: The advances in data, computing and networking over the last two decades led to a shift in many application domains that includes machine learning on big data as a part of the scientific process, requiring new capabilities for integrated and distributed hardware and software infrastructure. This paper contributes a workflow-driven approach for dynamic data-driven application development on top of a new kind of networked Cyberinfrastructure called CHASE-CI. In particular, we present: 1) The architecture for CHASE-CI, a network of distributed fast GPU appliances for machine learning and storage managed through Kubernetes on the high-speed (10-100Gbps) Pacific Research Platform (PRP); 2) A machine learning software containerization approach and libraries required for turning such a network into a distributed computer for big data analysis; 3) An atmospheric science case study that can only be made scalable with an infrastructure like CHASE-CI; 4) Capabilities for virtual cluster management for data communication and analysis in a dynamically scalable fashion, and visualization across the network in specialized visualization facilities in near real-time; and, 5) A step-by-step workflow and performance measurement approach that enables taking advantage of the dynamic architecture of the CHASE-CI network and container management infrastructure.


title: Pion induced Drell-Yan: the transverse momentum structure of the pion

abstract: We present an analysis of unpolarized Drell-Yan pair production in pion-nucleus scattering with a particular focus into the pion dynamics. The study consists in analyzing the effect of the partonic longitudinal and, especially, transverse distributions of the pion in a Nambu--Jona-Lasinio (NJL) framework, with Pauli-Villars regularization. In order to consistently take into account the QCD evolution effects, we have estimated the hadronic scale corresponding to the NJL model's degrees of freedom through a minimization procedure at NLO: The NLO evolved pion distributions have been compared to rapidity differential Drell-Yan cross sections data. That hadronic scale so determined represents the only free parameter in our approach. The NJL transverse momentum PDF, evolved up to next-to-leading logarithmic accuracy, is then tested against the transverse momentum spectrum of dilepton pairs up to a transverse momentum of 2 GeV. We found a fair agreement with available pion-nucleus data. We find sizable evolution effects on the shape of the distributions and on the generated average transverse momentum of the dilepton pair.


title: Regularized and Distributionally Robust Data-Enabled Predictive Control

abstract: In this paper, we study a data-enabled predictive control (DeePC) algorithm applied to unknown stochastic linear time-invariant systems. The algorithm uses noise-corrupted input/output data to predict future trajectories and compute optimal control policies. To robustify against uncertainties in the input/output data, the control policies are computed to minimize a worst-case expectation of a given objective function. Using techniques from distributionally robust stochastic optimization, we prove that for certain objective functions, the worst-case optimization problem coincides with a regularized version of the DeePC algorithm. These results support the previously observed advantages of the regularized algorithm and provide probabilistic guarantees for its performance. We illustrate the robustness of the regularized algorithm through a numerical case study.


title: Performance of ePix10K, a high dynamic range, gain auto-ranging pixel detector for FELs

abstract: ePix10K is a hybrid pixel detector developed at SLAC for demanding free-electron laser (FEL) applications, providing an ultrahigh dynamic range (245 eV to 88 MeV) through gain auto-ranging. It has three gain modes (high, medium and low) and two auto-ranging modes (high-to-low and medium-to-low). The first ePix10K cameras are built around modules consisting of a sensor flip-chip bonded to 4 ASICs, resulting in 352x384 pixels of 100 $\mu$m x 100 $\mu$m each. We present results from extensive testing of three ePix10K cameras with FEL beams at LCLS, resulting in a measured noise floor of 245 eV rms, or 67 e$^-$ equivalent noise charge (ENC), and a range of 11000 photons at 8 keV. We demonstrate the linearity of the response in various gain combinations: fixed high, fixed medium, fixed low, auto-ranging high to low, and auto-ranging medium-to-low, while maintaining a low noise (well within the counting statistics), a very low cross-talk, perfect saturation response at fluxes up to 900 times the maximum range, and acquisition rates of up to 480 Hz. Finally, we present examples of high dynamic range x-ray imaging spanning more than 4 orders of magnitude dynamic range (from a single photon to 11000 photons/pixel/pulse at 8 keV). Achieving this high performance with only one auto-ranging switch leads to relatively simple calibration and reconstruction procedures. The low noise levels allow usage with long integration times at non-FEL sources. ePix10K cameras leverage the advantages of hybrid pixel detectors with high production yield and good availability, minimize development complexity through sharing the hardware, software and DAQ development with all other versions of ePix cameras, while providing an upgrade path to 5 kHz, 25 kHz and 100 kHz in three steps over the next few years, matching the LCLS-II requirements.


title: Do the Chinese Astronomical Records Dated AD 776 January 12/13 Describe an Auroral Display or a Lunar Halo? A Critical Re-examination

abstract: The enhancement of carbon-14 in tree rings around AD 774/775 has generated wide interest in solar activity at that time. The historical auroral records have been examined critically. Of particular interest was the "white vapour" observed in China on AD 776 January 12/13. Both Usoskin et al. (2013, Astron. Astrophys. 55, L3; U13) and Stephenson (2015, Adv. Sp. Res. 55, 1537; S15) interpreted this record as an auroral display. Subsequently, Neuh\"auser and Neuh\"auser (2015, Astron. Nachr. 336, 225; NN15) proposed five "criteria" for the likeliness of aurorae and on this basis rejected an auroral interpretation. Instead, they interpreted it as a lunar halo, and suggested there were no auroral records as a proxy of solar activity in the interval AD 774-785. We consider if their "lunar halo hypothesis" and their auroral criteria could be of use in future researches on historical auroral candidates. We first show a counter-example for the lunar halo hypothesis from a parallel record on 1882 November 17, which was seen as a whitish colour, in the southerly direction, and near the Moon. We then consider NN15's criteria on colour, direction, and sky brightness and investigate other counter-examples from early-modern auroral observations. We also consider the extension of the white vapour in AD 776 according to the distribution of Chinese asterisms, and show that its large extension was inconsistent with the lunar halo hypothesis. Conversely, the streaks of white vapour penetrating the eight Chinese asterisms can be reproduced if we consider auroral-ray structures at altitudes between 97 km and 170 km, along geomagnetic field lines between the L-shells L=1.55 and 1.64. Our investigations show that we should consider candidate auroral records in historical documents not on the basis of the newly suggested a priori criteria by NN15 but on all the available observational evidence.


title: Generating urban morphologies at large scales

abstract: At large scales, typologies of urban form and corresponding generating processes remain an open question with important implications regarding urban planning policies and sustainability. We propose in this paper to generate urban configurations at large scales, typically of districts, with morphogenesis models, and compare these to real configurations according to morphological indicators. Real values are computed on a large sample of districts taken in European urban areas. We calibrate each model and show their complementarity to approach the variety of real urban configurations, paving the way to multi-model approaches of urban morphogenesis.


title: SOFIA {\it FORCAST} Photometry of 12 Extended Green Objects in the Milky Way

abstract: Massive young stellar objects are known to undergo an evolutionary phase in which high mass accretion rates drive strong outflows. A class of objects believed to trace this phase accurately is the GLIMPSE Extended Green Object (EGO) sample, so named for the presence of extended 4.5 $\mu$m emission on sizescales of $\sim$0.1 pc in \textit{Spitzer} images. We have been conducting a multi-wavelength examination of a sample of 12 EGOs with distances of 1 to 5 kpc. In this paper, we present mid-infrared images and photometry of these EGOs obtained with the SOFIA telescope, and subsequently construct SEDs for these sources from the near-IR to sub-millimeter regimes using additional archival data. We compare the results from greybody models and several publicly-available software packages which produce model SEDs in the context of a single massive protostar. The models yield typical \rstar\/$\sim$10 \rsun, \tstar\/$\sim$10$^3$ to 10$^4$ K, and \lstar\/$\sim$1~$-$~40~$\times$~10$^3$ \lsun; the median $L/M$ for our sample is 24.7 \lsun/\msun. Model results rarely converge for \rstar\/ and \tstar, but do for \lstar, which we take to be an indication of the multiplicity and inherently clustered nature of these sources even though, typically, only a single source dominates in the mid-infrared. The median $L/M$ value for the sample suggests that these objects may be in a transitional stage between the commonly described ``IR-quiet'' and ``IR-bright'' stages of MYSO evolution. The median $T_{dust}$ for the sample is less conclusive, but suggests that these objects are either in this transitional stage or occupy the cooler (and presumably younger) part of the IR-bright stage.


title: Energy-corrected FEM and explicit time-stepping for parabolic problems

abstract: The presence of corners in the computational domain, in general, reduces the regularity of solutions of parabolic problems and diminishes the convergence properties of the finite element approximation introducing a so-called "pollution effect". Standard remedies based on mesh refinement around the singular corner result in very restrictive stability requirements on the time-step size when explicit time integration is applied. In this article, we introduce and analyse the energy-corrected finite element method for parabolic problems, which works on quasi-uniform meshes, and, based on it, create fast explicit time discretisation. We illustrate these results with extensive numerical investigations not only confirming the theoretical results but also showing the flexibility of the method, which can be applied in the presence of multiple singular corners and a three-dimensional setting. We also propose a fast explicit time-stepping scheme based on a piecewise cubic energy-corrected discretisation in space completed with mass-lumping techniques and numerically verify its efficiency.


title: Astro2020: Unleashing the Potential of Dust Emission as a Window onto Galaxy Evolution

abstract: We present the severe, systematic uncertainties currently facing our understanding of dust emission, which stymie our ability to truly exploit dust as a tool for studying galaxy evolution. We propose a program of study to tackle these uncertainties, describe the necessary facilities, and discuss the potential science gains that will result. This white paper was submitted to the US National Academies' Astro2020 Decadal Survey on Astronomy and Astrophysics.


title: Multi-camera calibration with pattern rigs, including for non-overlapping cameras: CALICO

abstract: This paper describes CALICO, a method for multi-camera calibration suitable for challenging contexts: stationary and mobile multi-camera systems, cameras without overlapping fields of view, and non-synchronized cameras. Recent approaches are roughly divided into infrastructure- and pattern-based. Infrastructure-based approaches use the scene's features to calibrate, while pattern-based approaches use calibration patterns. Infrastructure-based approaches are not suitable for stationary camera systems, and pattern-based approaches may constrain camera placement because shared fields of view or extremely large patterns are required. CALICO is a pattern-based approach, where the multi-calibration problem is formulated using rigidity constraints between patterns and cameras. We use a {\it pattern rig}: several patterns rigidly attached to each other or some structure. We express the calibration problem as that of algebraic and reprojection error minimization problems. Simulated and real experiments demonstrate the method in a variety of settings. CALICO compared favorably to Kalibr. Mean reconstruction accuracy error was $\le 0.71$ mm for real camera rigs, and $\le 1.11$ for simulated camera rigs. Code and data releases are available at \cite{tabb_amy_2019_3520866} and \url{https://github.com/amy-tabb/calico}.


title: Splitting Algorithms for Rare Events of Semimartingale Reflecting Brownian motions

abstract: We study rare event simulations of semimartingale reflecting Brownian motions (SRBMs) in an orthant. The rare event of interest is that a $d$-dimensional positive recurrent SRBM enters the set $B_n = \{z\in\mathbb{R}^d: \sum_{k=1}^d z_k = n\}$ before reaching a small neighborhood of the origin as $n\to\infty$. We show that under a proper scaling and some regularity conditions, the probability of interest satisfies a large deviation principle. We then construct a subsolution to the variational problem for our rare event, and based on this subsolution construct particle based simulation algorithms to estimate the probability of the rare event. It is shown that the proposed algorithm is stable and theoretically superior to standard Monte Carlo for a broad class of positive recurrent SRBMs.


title: Machine Learning the Voltage of Electrode Materials in Metal-ion Batteries

abstract: Machine learning (ML) techniques have rapidly found applications in many domains of materials chemistry and physics where large data sets are available. Aiming to accelerate the discovery of materials for battery applications, in this work, we develop a tool (http://se.cmich.edu/batteries) based on ML models to predict voltages of electrode materials for metal-ion batteries. To this end, we use deep neural network, support vector machine, and kernel ridge regression as ML algorithms in combination with data taken from the Materials Project Database, as well as feature vectors from properties of chemical compounds and elemental properties of their constituents. We show that our ML models have predictive capabilities for different reference test sets and, as an example, we utilize them to generate a voltage profile diagram and compare it to density functional theory calculations. In addition, using our models, we propose nearly 5,000 candidate electrode materials for Na- and K-ion batteries. We also make available a web-accessible tool that, within a minute, can be used to estimate the voltage of any bulk electrode material for a number of metal-ions. These results show that ML is a promising alternative for computationally demanding calculations as a first screening tool of novel materials for battery applications.


title: Generate What You Can't See - a View-dependent Image Generation

abstract: In order to operate autonomously, a robot should explore the environment and build a model of each of the surrounding objects. A common approach is to carefully scan the whole workspace. This is time-consuming. It is also often impossible to reach all the viewpoints required to acquire full knowledge about the environment. Humans can perform shape completion of occluded objects by relying on past experience. Therefore, we propose a method that generates images of an object from various viewpoints using a single input RGB image. A deep neural network is trained to imagine the object appearance from many viewpoints. We present the whole pipeline, which takes a single RGB image as input and returns a sequence of RGB and depth images of the object. The method utilizes a CNN-based object detector to extract the object from the natural scene. Then, the proposed network generates a set of RGB and depth images. We show the results both on a synthetic dataset and on real images.


title: Optimal Classification and Outlier Detection for Stripped-Envelope Core-Collapse Supernovae

abstract: In the current era of time-domain astronomy, it is increasingly important to have rigorous, data driven models for classifying transients, including supernovae. We present the first application of Principal Component Analysis to the spectra of stripped-envelope core-collapse supernovae. We use one of the largest compiled optical datasets of stripped-envelope supernovae, containing 160 SNe and 1551 spectra. We find that the first 5 principal components capture 79\% of the variance of our spectral sample, which contains the main families of stripped supernovae: Ib, IIb, Ic and broad-lined Ic. We develop a quantitative, data-driven classification method using a support vector machine, and explore stripped-envelope supernovae classification as a function of phase relative to V-band maximum light. Our classification method naturally identifies "transition" supernovae and supernovae with contested labels, which we discuss in detail. We find that the stripped-envelope supernovae types are most distinguishable in the later phase ranges of $10\pm5$ days and $15\pm5$ days relative to V-band maximum, and we discuss the implications of our findings for current and future surveys such as ZTF and LSST.


title: On the Use of Planetary Science Data for Studying Extrasolar Planets

abstract: There is an opportunity to advance both solar system and extrasolar planetary studies that does not require the construction of new telescopes or new missions but better use and access to inter-disciplinary data sets. This approach leverages significant investment from NASA and international space agencies in exploring this solar system and using those discoveries as "ground truth" for the study of extrasolar planets. This white paper illustrates the potential, using phase curves and atmospheric modeling as specific examples. A key advance required to realize this potential is to enable seamless discovery and access within and between planetary science and astronomical data sets. Further, seamless data discovery and access also expands the availability of science, allowing researchers and students at a variety of institutions, equipped only with Internet access and a decent computer to conduct cutting-edge research.


title: The Largest Prime Dividing the Maximal Order of an Element of $S_n$

abstract: We define $g(n)$ to be the maximal order of an element of the symmetric group on $n$ elements. Results about the prime factorization of $g(n)$ allow a reduction of the upper bound on the largest prime divisor of $g(n)$ to $1.328\sqrt{n\log n}$.


title: Generalization from correlated sets of patterns in the perceptron

abstract: Generalization is a central aspect of learning theory. Here, we propose a framework that explores an auxiliary task-dependent notion of generalization, and attempts to quantitatively answer the following question: given two sets of patterns with a given degree of dissimilarity, how easily will a network be able to "unify" their interpretation? This is quantified by the volume of the configurations of synaptic weights that classify the two sets in a similar manner. To show the applicability of our idea in a concrete setting, we compute this quantity for the perceptron, a simple binary classifier, using the classical statistical physics approach in the replica-symmetric ansatz. In this case, we show how an analytical expression measures the "distance-based capacity", the maximum load of patterns sustainable by the network, at fixed dissimilarity between patterns and fixed allowed number of errors. This curve indicates that generalization is possible at any distance, but with decreasing capacity. We propose that a distance-based definition of generalization may be useful in numerical experiments with real-world neural networks, and to explore computationally sub-dominant sets of synaptic solutions.


title: Frictional Damping in Biomimetic Scale Beam Oscillations

abstract: Stiff scales adorn the exterior surfaces of fishes, snakes, and many reptiles. They provide protection from external piercing attacks and control over global deformation behavior to aid locomotion, slithering, and swimming across a wide range of environmental condition. In this letter, we investigate the dynamic behavior of biomimetic scale substrates for further understanding the origins of the nonlinearity that involve various aspect of scales interaction, sliding kinematics, interfacial friction, and their combination. Particularly, we study the vibrational characteristics through an analytical model and numerical investigations for the case of a simply supported scale covered beam. Our results reveal for the first time that biomimetic scale beams exhibit viscous damping behavior even when only Coulomb friction is postulated for free vibrations. We anticipate and quantify the anisotropy in the damping behavior with respect to curvature. We also find that unlike static pure bending where friction increases bending stiffness, a corresponding increase in natural frequency for the dynamic case does not arise for simply supported beam. Since both scale geometry, distribution and interfacial properties can be easily tailored, our study indicates a biomimetic strategy to design exceptional synthetic materials with tailorable damping behavior.


title: Frobenius Pseudoprimes

abstract: The proliferation of probable prime tests in recent years has produced a plethora of definitions with the word ``pseudoprime'' in them. Examples include pseudoprimes, Euler pseudoprimes, strong pseudoprimes, Lucas pseudoprimes, strong Lucas pseudoprimes, extra strong Lucas pseudoprimes and Perrin pseudoprimes. Though these tests represent a wealth of ideas, they exist as a hodge-podge of definitions rather than as examples of a more general theory. It is the goal of this paper to present a way of viewing many of these tests as special cases of a general principle, as well as to re-formulate them in the context of finite fields. One aim of the reformulation is to enable the creations of stronger tests; another is to aid in proving results about large classes of pseudoprimes.


title: Availability of Hyperlinked Resources in Astrophysics Papers

abstract: Astrophysics papers often rely on software which may or may not be available, and URLs are often used as proxy citations for software and data. We extracted all URLs from two journals' 2015 research articles, removed those from certain long-term reliable domains, and tested the remainder to determine what percentage of these URLs were accessible in October 2018.


title: Fundamentals of Power Allocation Strategies for Downlink Multi-user NOMA with Target Rates

abstract: For downlink multi-user non-orthogonal multiple access (NOMA) systems with successive interference cancellation (SIC) receivers, and a base-station not possessing the instantaneous channel gains, the fundamental relationship between the target rates and power allocation is investigated. It is proven that the total interference from signals not removed by SIC has a fundamental upper-limit which is a function of the target rates, and the outage probability is one when exceeding this limit. The concept of well-behaved power allocation strategies is defined, and its properties are proven to be derived solely based on the target rates. The existence of power allocation strategies that enable NOMA to outperform OMA in per-user outage probability is proven, and are always well-behaved for the case when the outage probability performance of NOMA and OMA are equal for all users. The proposed SIC decoding order is then shown to the most energy efficient. The derivation of well-behaved power allocation strategies that have improved outage probability performance over OMA for each user is outlined. Simulations validate the theoretical results, demonstrating that NOMA systems can always outperform OMA systems in outage probability performance, without relying on the exact channel gains.


title: A Probable Prime Test With High Confidence

abstract: Monier and Rabin proved that an odd composite can pass the Strong Probable Prime Test for at most $\frac 14$ of the possible bases. In this paper, a probable prime test is developed using quadratic polynomials and the Frobenius automorphism. The test, along with a fixed number of trial divisions, ensures that a composite $n$ will pass for less than $\frac 1{7710}$ of the polynomials $x^2-bx-c$ with $\left(b^2+4c\over n\right)=-1$ and $\left(-c\over n\right)=1$. The running time of the test is asymptotically $3$ times that of the Strong Probable Prime Test.


title: LOFAR observations of 4C+19.44. On the discovery of low frequency spectral curvature in relativistic jet knots

abstract: We present the first LOFAR observations of the radio jet in the quasar 4C+19.44 (a.k.a. PKS 1354+19) obtained with the long baselines. The achieved resolution is very well matched to that of archival Jansky Very Large Array (JVLA) observations at higher radio frequencies as well as the archival X-ray images obtained with {\it Chandra}. We found that, for several knots along the jet, the radio flux densities measured at hundreds of MHz lie well below the values estimated by extrapolating the GHz spectra. This clearly indicates the presence of spectral curvature. Radio spectral curvature has been already observed in different source classes and/or extended radio structures and it has been often interpreted as due to intrinsic processes, as a curved particle energy distribution, rather than absorption mechanisms ({ Razin-Tsytovich} effect, free-free or synchrotron self absorption to name a few). Here we discuss our results according to the scenario where particles undergo stochastic acceleration mechanisms also in quasar jet knots.


title: There are Infinitely Many Perrin Pseudoprimes

abstract: This paper proves the existence of infinitely many Perrin pseudoprimes, as conjectured by Adams and Shanks in 1982. The theorem proven covers a general class of pseudoprimes based on recurrence sequences. The result uses ingredients of the proof of the infinitude of Carmichael numbers, along with zero-density estimates for Hecke L-functions.


title: A Universality Law For Sign Correlations of Eigenfunctions of Differential Operators

abstract: We establish a universality law for sequences of functions $\{w_n\}_{n \in \mathbb{N}}$ satisfying a form of WKB approximation on compact intervals. This includes eigenfunctions of generic Schr\"odinger operators, as well as Laguerre and Chebyshev polynomials. Given two distinct points $x, y \in \mathbb{R}$, we ask how often do $w_n(x)$ and $w_n(y)$ have the same sign. Asymptotically, one would expect this to be true half the time, but this turns out to not always be the case. Under certain natural assumptions, we prove that, for all $x \neq y$, $$ \frac{1}{3} \leq \lim_{N \to \infty} \frac{1}{N} \# \left\{0 \leq n < N: \mathrm{sgn}(w_n(x)) = \mathrm{sgn}(w_n(y)) \right\} \leq \frac{2}{3}, $$ and that these bounds are optimal, and can be attained. Our methods extend to other questions of similar flavor and we also discuss a number of open problems.


title: Does Homophily Make Socialbots More Influential? Exploring Infiltration Strategies

abstract: Socialbots are intelligent software controlling all the behavior of fake accounts in an online social network. They use artificial intelligence techniques to pass themselves off as human social media users. Socialbots exploit user trust to achieve their malicious goals, such as astroturfing, performing Sybil attacks, spamming, and harvesting private data. The first phase to countermeasure the malicious activities of the socialbots is studying their characteristics and revealing strategies they can employ to successfully infiltrate stealthily into target online social network. In this paper,we investigate the success of using different infiltration strategies in terms of infiltration performance and being stealthy. Every strategy is characterized by socialbots profile and behavioral characteristics. The findings from this study illustrate that assuming a specific taste for the tweets a socialbot retweets and/or likes make it more influential. Furthermore, the experimental results indicate that considering the presence of common characteristics and similarity increase the probability of being followed by other users. This is in complete agreement with homophily concept which is the tendency of individuals to associate and bond with similar others in social networks.


title: Data-driven Identification and Prediction of Power System Dynamics Using Linear Operators

abstract: In this paper, we propose linear operator theoretic framework involving Koopman operator for the data-driven identification of power system dynamics. We explicitly account for noise in the time series measurement data and propose robust approach for data-driven approximation of Koopman operator for the identification of nonlinear power system dynamics. The identified model is used for the prediction of state trajectories in the power system. The application of the framework is illustrated using an IEEE nine bus test system.


title: Shape of the quark gluon plasma droplet reflected in the high-p_T data

abstract: We show, through analytic arguments, numerical calculations, and comparison with experimental data, that the ratio of the high-p_T observables v_2/(1-R_AA) reaches a well-defined saturation value at high p_T, and that this ratio depends only on the spatial anisotropy of the quark gluon plasma (QGP) formed in ultrarelativistic heavy-ion collisions. With expected future reduction of experimental errors, the anisotropy extracted from experimental data will further constrain the calculations of initial particle production in heavy-ion collisions and thus test our understanding of QGP physics.


title: The Geometrical Trinity of Gravity

abstract: The geometrical nature of gravity emerges from the universality dictated by the equivalence principle. In the usual formulation of General Relativity, the geometrisation of the gravitational interaction is performed in terms of the spacetime curvature, which is now the standard interpretation of gravity. However, this is not the only possibility. In these notes we discuss two alternative, though equivalent, formulations of General Relativity in flat spacetimes, in which gravity is fully ascribed either to torsion or to non-metricity, thus putting forward the existence of three seemingly unrelated representations of the same underlying theory. Based on these three alternative formulations of General Relativity, we then discuss some extensions.


title: Simultaneous amplitude and phase modulation for wide-field nonlinear microscopy applications

abstract: In wide-field nonlinear microscopy, wavefront modulation by means of phase-only spatial light modulators (SLMs) allows achieving simultaneous two-photon excitation and fluorescence emission from specific region-of-interests (ROIs) of biological specimens. This is basically accomplished at the illumination path of the microscope by the reconstruction of computer generated holograms (CGHs) onto the sample plane. However, as two-photon absorption (TPA) is inherently an intensity-square dependent process and iterative Fourier transform algorithms (IFTAs) can only approximate the illumination of selected ROIs with the reconstructed CGHs, both signal acquisition and/or image formation can be largely affected by the spatial irregularities of the illumination patterns. In addition, the speckle associated with the superposition of coherent light at the selected ROIs prevents illumination strategies based on CGHs to be successfully used for large-area (more than 50x50 $\mu$m2) excitation tasks. To overcome these limitations, we propose an alternative complex illumination method (CIM) able to generate simultaneous nonlinear excitation of large-area ROIs with full control over the amplitude and phase of the optical wavefront. We experimentally demonstrate spatially uniform illumination, as well as structured illumination with user-defined intensity levels onto micrometric but large-area ROIs. Furthermore, a proof-of-concept experiment on wide-field second harmonic generation (SHG) is provided. We believe that the proposed CIM could find applications in wide-field nonlinear microscopy, particularly for speed up signal acquisition time or improve two-photon image formation.


title: Spinor field singular functions in QED with strong external backgrounds

abstract: We construct and study singular functions in strong-field $QED$ with two external electromagnetic fields that represent principally different types of external backgrounds, the first one belongs to the class of so-called $t$-potential electric steps (electric-like fields that are switched on and off at initial and final time instants), and the second one belongs to the class of so-called $x$-potential electric steps (time-independent electric-like fields of constant direction that are concentrated in a restricted spatial area). As the first background ($T$-constant electric field) is chosen an uniform electric field which acts during a finite time interval $T$ , whereas as the second background ($L$-constant electric field) is chosen a constant electric field confined between two capacitor plates separated by a large distance $L$. For the both cases we find \textrm{in}- and \textrm{out}-solutions of the Dirac equation in terms of light cone variables. With the help of these solutions, we construct Fock-Schwinger proper-time integral representations for all the singular functions that provide nonperturbative (with respect to the external backgrounds) calculations of any transition amplitudes and mean values of any physical quantities. Considering calculations in the $T$-constant field and in the $L$-constant field as different regularizations of the corresponding calculations in the constant uniform electric field, we have demonstrated their equivalence for sufficiently large $T$ and $L$.


title: Exploring A Cosmic-Ray Origin of the Multi-wavelength Emission in M31

abstract: A recent detection of spatially extended gamma-ray emission in the central region of the Andromeda galaxy (M31) has led to several possible explanations being put forth, including dark matter annihilation and millisecond pulsars. Another possibility is that the emission in M31 can be accounted for with a purely astrophysical cosmic-ray (CR) scenario. This scenario would lead to a rich multi-wavelength emission that can, in turn, be used to test it. Relativistic cosmic-ray electrons (CRe) in magnetic fields produce radio emission through synchrotron radiation, while X-rays and gamma rays are produced through inverse Compton scattering. Additionally, collisions of primary cosmic-ray protons (CRp) in the interstellar medium produce charged and neutral pions that then decay into secondary CRe (detectable through radiative processes) and gamma-rays. Here, we explore the viability of a CR origin for multi-wavelength emission in M31, taking into consideration three scenarios: a CR scenario dominated by primary CRe, one dominated by CRp and the resulting secondary CRe and gamma rays from neutral pion decay, and a final case in which both of these components exist simultaneously. We find that the multi-component model is the most promising, and is able to fit the multi-wavelength spectrum for a variety of astrophysical parameters consistent with previous studies of M31 and of cosmic-ray physics. However, the CR power injection implied by our models exceeds the estimated CR power injection from typical astrophysical cosmic-ray sources such as supernovae.


title: LSST Narrowband Filters

abstract: In this white paper, we present the scientific cases for adding narrowband optical filters to the Large Synoptic Survey Telescope (LSST). LSST is currently planning to observe the southern sky in 6 broadband optical filters. Three of the four LSST science themes would benefit from adding narrowband filter observations. We discuss the technical considerations of using narrowband filters with the LSST and lay out the scientific impact that would result on the study of AGB stars, emission line nebula (e.g., supernova remnants and planetary nebulae), photometric redshifts of galaxies, and the determination of stellar parameters.


title: Quantitative Ultrasound and B-mode Image Texture Features Correlate with Collagen and Myelin Content in Human Ulnar Nerve Fascicles

abstract: We investigate the usefulness of quantitative ultrasound (QUS) and B-mode texture features for characterization of ulnar nerve fascicles. Ultrasound data were acquired from cadaveric specimens using a nominal 30 MHz probe. Next, the nerves were extracted to prepare histology sections. 85 fascicles were matched between the B-mode images and the histology sections. For each fascicle image, we selected an intra-fascicular region of interest. We used histology sections to determine features related to the concentration of collagen and myelin, and ultrasound data to calculate backscatter coefficient (-24.89 dB $\pm$ 8.31), attenuation coefficient (0.92 db/cm-MHz $\pm$ 0.04), Nakagami parameter (1.01 $\pm$ 0.18) and entropy (6.92 $\pm$ 0.83), as well as B-mode texture features obtained via the gray level co-occurrence matrix algorithm. Significant Spearman's rank correlations between the combined collagen and myelin concentrations were obtained for the backscatter coefficient (R=-0.68), entropy (R=-0.51), and for several texture features. Our study demonstrates that QUS may potentially provide information on structural components of nerve fascicles.


title: Detecting GAN generated Fake Images using Co-occurrence Matrices

abstract: The advent of Generative Adversarial Networks (GANs) has brought about completely novel ways of transforming and manipulating pixels in digital images. GAN based techniques such as Image-to-Image translations, DeepFakes, and other automated methods have become increasingly popular in creating fake images. In this paper, we propose a novel approach to detect GAN generated fake images using a combination of co-occurrence matrices and deep learning. We extract co-occurrence matrices on three color channels in the pixel domain and train a model using a deep convolutional neural network (CNN) framework. Experimental results on two diverse and challenging GAN datasets comprising more than 56,000 images based on unpaired image-to-image translations (cycleGAN [1]) and facial attributes/expressions (StarGAN [2]) show that our approach is promising and achieves more than 99% classification accuracy in both datasets. Further, our approach also generalizes well and achieves good results when trained on one dataset and tested on the other.


title: Visual recognition in the wild by sampling deep similarity functions

abstract: Recognising relevant objects or object states in its environment is a basic capability for an autonomous robot. The dominant approach to object recognition in images and range images is classification by supervised machine learning, nowadays mostly with deep convolutional neural networks (CNNs). This works well for target classes whose variability can be completely covered with training examples. However, a robot moving in the wild, i.e., in an environment that is not known at the time the recognition system is trained, will often face \emph{domain shift}: the training data cannot be assumed to exhaustively cover all the within-class variability that will be encountered in the test data. In that situation, learning is in principle possible, since the training set does capture the defining properties, respectively dissimilarities, of the target classes. But directly training a CNN to predict class probabilities is prone to overfitting to irrelevant correlations between the class labels and the specific subset of the target class that is represented in the training set. We explore the idea to instead learn a Siamese CNN that acts as similarity function between pairs of training examples. Class predictions are then obtained by measuring the similarities between a new test instance and the training samples. We show that the CNN embedding correctly recovers the relative similarities to arbitrary class exemplars in the training set. And that therefore few, randomly picked training exemplars are sufficient to achieve good predictions, making the procedure efficient.


title: Ultrafast processing of pixel detector data with machine learning frameworks

abstract: Modern photon science performed at high repetition rate free-electron laser (FEL) facilities and beyond relies on 2D pixel detectors operating at increasing frequencies (towards 100 kHz at LCLS-II) and producing rapidly increasing amounts of data (towards TB/s). This data must be rapidly stored for offline analysis and summarized in real time. While at LCLS all raw data has been stored, at LCLS-II this would lead to a prohibitive cost; instead, enabling real time processing of pixel detector raw data allows reducing the size and cost of online processing, offline processing and storage by orders of magnitude while preserving full photon information, by taking advantage of the compressibility of sparse data typical for LCLS-II applications. We investigated if recent developments in machine learning are useful in data processing for high speed pixel detectors and found that typical deep learning models and autoencoder architectures failed to yield useful noise reduction while preserving full photon information, presumably because of the very different statistics and feature sets between computer vision and radiation imaging. However, we redesigned in Tensorflow mathematically equivalent versions of the state-of-the-art, "classical" algorithms used at LCLS. The novel Tensorflow models resulted in elegant, compact and hardware agnostic code, gaining 1 to 2 orders of magnitude faster processing on an inexpensive consumer GPU, reducing by 3 orders of magnitude the projected cost of online analysis at LCLS-II. Computer vision a decade ago was dominated by hand-crafted filters; their structure inspired the deep learning revolution resulting in modern deep convolutional networks; similarly, our novel Tensorflow filters provide inspiration for designing future deep learning architectures for ultrafast and efficient processing and classification of pixel detector images at FEL facilities.


title: Hyper-Resistive Model of Ultra High Energy Cosmic Ray Acceleration by Magnetically Collimated Jets Created by Active Galactic Nuclei

abstract: This is the fourth in a series of companion papers showing that, when an efficient dynamo can be maintained by accretion disks around supermassive black holes in Active Galactic Nuclei (AGNs), it will lead to the formation of a powerful, magnetically-collimated helix that could explain both the observed jet/radiolobe structures on very large scales and ultimately the enormous power inferred from the observed ultra high energy cosmic rays (UHECRs) with energies > 10^19 eV. Many timescales are involved in this process. Our hyper-resistive magnetohydrodynamic (MHD) model provides a bridge between General Relativistic MHD simulations of dynamo formation, on the short accretion timescale, and observational evidence of magnetic collimation of large-scale jets on astrophysical timescales. Given the final magnetic structure, we apply hyper-resistive kinetic theory to show how instability causes slowly-evolving magnetically-collimated jets to become the most powerful relativistic accelerators in the Universe. The model yields nine observables in reasonable agreement with observations: the jet length, radiolobe radius and apparent opening angle as observed by synchrotron radiation; the synchrotron total power, synchrotron wavelengths and maximum electron energy (TeVs); and the maximum UHECR energy, the cosmic ray energy spectrum and the cosmic ray intensity on Earth.


title: How Gas Accretion Feeds Galactic Disks

abstract: Numerous observations indicate that galaxies need a continuous gas supply to fuel star formation and explain the star formation history. However, direct observational evidence of gas accretion remains rare. Using the EAGLE cosmological hydrodynamic simulation suite, we study cold gas accretion onto galaxies and the observational signatures of the cold gas kinematics. For EAGLE galaxies at z=0.27, we find that cold gas accretes onto galaxies anisotropically with typical inflow speeds between 20 km s$^{-1}$ and 60 km s$^{-1}$. Most of these galaxies have comparable mass inflow rates and star formation rates, implying that the cold inflowing gas plausibly accounts for sustaining the star-forming activities of the galaxies. As motivation for future work to compare the cold gas kinematics with measurements from quasar sightline observations, we select an EAGLE galaxy with an extended cold gas disk, and we probe the cold gas using mock quasar sightlines. We demonstrate that by viewing the disk edge-on, sightlines at azimuthal angles below 10$^{\circ}$ and impact parameters out to 60 pkpc can detect cold gas that corotates with the galaxy disk. This example suggests cold gas disks that extend beyond the optical disks possibly explain the sightline observations that detect corotating cold gas near galaxy major axes.


title: Novel critical behavior of magnetization in URhSi:Similarities to uranium ferromagnetic superconductors UGe$_2$ and URhGe

abstract: We study the critical behavior of dc magnetization in the uranium ferromagnet URhSi around the paramagnetic to ferromagnetic phase transition at T_C~ 10 K with a modified Arrott plot, a Kouvel-Fisher plot, the critical isotherm analysis and the scaling analysis. URhSi is isostructural to uranium ferromagnetic superconductors URhGe and UCoGe. The critical exponent beta for the temperature dependence of the spontaneous magnetization below TC, gamma for the magnetic susceptibility, and delta for the magnetic isotherm at T_C in URhSi have been determined as beta = 0.300 +- 0.002, gamma = 1.00 +- 0.02, and delta = 4.38 +- 0.04 by the scaling analysis and the critical isotherm analysis. These critical exponents fulfill the Widom scaling law delta = 1+ gamma/beta. Magnetization has strong uniaxial magnetic anisotropy in the ferromagnetic state of URhSi. However, the universality class of the ferromagnetic transition does not belong to the 3D Ising system with short-range exchange interactions between magnetic moments (beta = 0.325, gamma = 1.241, and delta = 4.82). The obtained exponents in URhSi are similar to those in the uranium ferromagnetic superconductors UGe_2 and URhGe, and uranium ferromagnets UIr and U(Co_0.98Os_0.02)Al. We have previously reported the unconventional critical behavior of magnetization in the uranium ferromagnetic superconductors [N. Tateiwa et al. Phys. Rev. B 89, 064420 (2014)]. The universality class of the ferromagnetic transition in URhSi may belong to the same one in the uranium ferromagnetic superconductors and the uranium ferromagnets. The unconventional critical behavior of the magnetization in the uranium compounds cannot be understood with previous theoretical interpretations of critical phenomena. The absence of the superconductivity in URhSi is discussed from several viewpoints.


title: Formulas for Data-driven Control: Stabilization, Optimality and Robustness

abstract: In a paper by Willems and coauthors it was shown that persistently exciting data can be used to represent the input-output behavior of a linear system. Based on this fundamental result, we derive a parametrization of linear feedback systems that paves the way to solve important control problems using data-dependent Linear Matrix Inequalities only. The result is remarkable in that no explicit system's matrices identification is required. The examples of control problems we solve include the state and output feedback stabilization, and the linear quadratic regulation problem. We also discuss robustness to noise-corrupted measurements and show how the approach can be used to stabilize unstable equilibria of nonlinear systems.


title: Estimates for $n$-widths of sets of smooth functions on complex spheres

abstract: In this work we investigate $n$-widths of multiplier operators $\Lambda_*$ and $\Lambda$, defined for functions on the complex sphere $\Omega_d$ of $\mathbb{C}^d$, associated with sequences of multipliers of the type $\{\lambda_{m,n}^*\}_{m,n\in \mathbb{N}}$, $\lambda_{m,n}^*=\lambda(m+n)$ and $\{\lambda_{m,n}\}_{m,n\in \mathbb{N}}$, $\lambda_{m,n}=\lambda(\max\{m,n\})$, respectively, for a bounded function $\lambda$ defined on $[0,\infty)$. If the operators $\Lambda_{*}$ and $\Lambda$ are bounded from $L^p(\Omega_d)$ into $L^q(\Omega_d)$, $1\leq p,q\leq\infty$, and $U_p$ is the closed unit ball of $L^p(\Omega_d)$, we study lower and upper estimates for the $n$-widths of Kolmogorov, linear, of Gelfand and of Bernstein, of the sets $\Lambda_{*}U_p$ and $\Lambda U_p$ in $L^q(\Omega_d)$. As application we obtain, in particular, estimates for the Kolmogorov $n$-width of classes of Sobolev, of finitely differentiable, infinitely differentiable and analytic functions on the complex sphere, in $L^q(\Omega_d)$, which are order sharp in various important situations.


title: Breaking the limit: Super-Eddington accretion onto black holes and neutron stars

abstract: With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (Volonteri & Rees, 2005). Indeed, for any stellar remnant seed, the rate required would need to exceed the Eddington limit. This is the point at which the outward force produced by radiation pressure is equal to the gravitational attraction experienced by the in-falling matter. In principle, this implies that there is a maximum luminosity an object of mass $M$ can emit; assuming spherical accretion and that the opacity is dominated by Thompson scattering, this Eddington luminosity is $L_{\rm{E}} = 1.38 \times 10^{38} (M/M_{\odot})$ erg s$^{-1}$. In reality, it is known that this limit can be violated, due to non-spherical geometry or various kinds of instabilities. Nevertheless, the Eddington limit remains an important reference point, and many of the details of how accretion proceeds above this limit remain unclear. Understanding how this so-called super-Eddington accretion occurs is of clear cosmological importance, since it potentially governs the growth of the first supermassive black holes (SMBHs) and the impact this growth would have had on their host galaxies (`feedback') and the epoch of reionization, as well as improving our understanding of accretion physics more generally.


title: Hammerhead, an ultrahigh resolution ePix camera for wavelength-dispersive spectrometers

abstract: Wavelength-dispersive spectrometers (WDS) are often used in synchrotron and FEL applications where high energy resolution (in the order of eV) is important. Increasing WDS energy resolution requires increasing spatial resolution of the detectors in the dispersion direction. The common approaches with strip detectors or small pixel detectors are not ideal. We present a novel approach, with a sensor using rectangular pixels with a high aspect ratio (between strips and pixels, further called "strixels"), and strixel redistribution to match the square pixel arrays of typical ASICs while avoiding the considerable effort of redesigning ASICs. This results in a sensor area of 17.4 mm x 77 mm, with a fine pitch of 25 $\mu$m in the horizontal direction resulting in 3072 columns and 176 rows. The sensors use ePix100 readout ASICs, leveraging their low noise (43 e$^-$, or 180 eV rms). We present results obtained with a Hammerhead ePix100 camera, showing that the small pitch (25 $\mu$m) in the dispersion direction maximizes performance for both high and low photon occupancies, resulting in optimal WDS energy resolution. The low noise level at high photon occupancy allows precise photon counting, while at low occupancy, both the energy and the subpixel position can be reconstructed for every photon, allowing an ultrahigh resolution (in the order of 1 $\mu$m) in the dispersion direction and rejection of scattered beam and harmonics. Using strixel sensors with redistribution and flip-chip bonding to standard ePix readout ASICs results in ultrahigh position resolution ($\sim$1 $\mu$m) and low noise in WDS applications, leveraging the advantages of hybrid pixel detectors (high production yield, good availability, relatively inexpensive) while minimizing development complexity through sharing the ASIC, hardware, software and DAQ development with existing versions of ePix cameras.


title: Directional PointNet: 3D Environmental Classification for Wearable Robotics

abstract: Environmental information can provide reliable prior information about human motion intent, which can aid the subject with wearable robotics to walk in complex environments. Previous researchers have utilized 1D signal and 2D images to classify environments, but they may face the problems of self-occlusion. Comparatively, 3D point cloud can be more appropriate to depict environments, thus we propose a directional PointNet to classify 3D point cloud directly. By utilizing the orientation information of the point cloud, the directional PointNet can classify daily terrains, including level ground, up stairs, and down stairs, and the classification accuracy achieves 99% for testing set. Moreover, the directional PointNet is more efficient than the previous PointNet because the T-net, which is utilized to estimate the transformation of the point cloud, is removed in this research and the length of the global feature is optimized. The experimental results demonstrate that the directional PointNet can classify the environments robustly and efficiently.


title: Safe Coordination of Human-Robot Firefighting Teams

abstract: Wildfires are destructive and inflict massive, irreversible harm to victims' lives and natural resources. Researchers have proposed commissioning unmanned aerial vehicles (UAVs) to provide firefighters with real-time tracking information; yet, these UAVs are not able to reason about a fire's track, including current location, measurement, and uncertainty, as well as propagation. We propose a model-predictive, probabilistically safe distributed control algorithm for human-robot collaboration in wildfire fighting. The proposed algorithm overcomes the limitations of prior work by explicitly estimating the latent fire propagation dynamics to enable intelligent, time-extended coordination of the UAVs in support of on-the-ground human firefighters. We derive a novel, analytical bound that enables UAVs to distribute their resources and provides a probabilistic guarantee of the humans' safety while preserving the UAVs' ability to cover an entire fire.


title: Remarks on Enveloping Semigroups

abstract: The local structures of enveloping semigroups of simple groups are investigated. All J-coirreducible connected stabilizer submonoids are determined. The notion of a navel of a reductive monoid is introduced. The cross-section lattice of the enveloping monoid is shown to be atomic. In type A, the generating series for the number of $G\times G$-orbits is found.


title: The Betti Numbers of a Determinantal Variety

abstract: We determine the Poincar\'e polynomial of the determinantal variety $\{\det = 0\}$ in the projective space associated with the monoid of $n\times n$ matrices.


title: A Bottom-up Approach to Testing Hypotheses That Have a Branching Tree Dependence Structure, with False Discovery Rate Control

abstract: Modern statistical analyses often involve testing large numbers of hypotheses. In many situations, these hypotheses may have an underlying tree structure that not only helps determine the order that tests should be conducted but also imposes a dependency between tests that must be accounted for. Our motivating example comes from testing the association between a trait of interest and groups of microbes that have been organized into operational taxonomic units (OTUs) or amplicon sequence variants (ASVs). Given p-values from association tests for each individual OTU or ASV, we would like to know if we can declare that a certain species, genus, or higher taxonomic grouping can be considered to be associated with the trait. For this problem, a bottom-up testing algorithm that starts at the lowest level of the tree (OTUs or ASVs) and proceeds upward through successively higher taxonomic groupings (species, genus, family etc.) is required. We develop such a bottom-up testing algorithm that controls the error rate of decisions made at higher levels in the tree, conditional on findings at lower levels in the tree. We further show this algorithm controls the false discovery rate based on the global null hypothesis that no taxa are associated with the trait. By simulation, we also show that our approach is better at finding driver taxa, the highest level taxa below which there are dense association signals. We illustrate our approach using data from a study of the microbiome among patients with ulcerative colitis and healthy controls.


title: Massive Neutron Stars with a Color Superconducting Quark Matter Core

abstract: We construct the equation of state for high density neutron star matter at zero temperature using the two-flavor Nambu--Jona-Lasinio (NJL) model as an effective theory of QCD. We build nuclear matter, quark matter, and the mixed phases from the same NJL Lagrangian, which has been used to model free and in-medium hadrons as well as nuclear systems. A focus here is to determine if the same coupling constants in the scalar diquark and vector meson channels, which give a good description of nucleon structure and nuclear matter, can also be used for the color superconducting high density quark matter phase. We find that this is possible for the scalar diquark (pairing) interaction, but the vector meson interaction has to be reduced so that superconducting quark matter becomes the stable phase at high densities. We compare our equation of state with recent phenomenological parametrizations based on generic stability conditions for neutron stars. We find that the maximum mass of a neutron star, with a color superconducting quark matter core, exceeds $2.01 \pm 0.04\,M_\odot$ which is the value of the recently observed massive neutron star PSR J0348+0432. The mass-radius relation is also consistent with gravitational wave observations (GW170817).


title: Long-time asymptotic behavior for the discrete defocusing mKdV equation

abstract: In this article, we apply Deift-Zhou nonlinear steepest descent method to analyze the long-time asymptotic behavior of the solution for the discrete defocusing mKdV equation. This equation was proposed by Ablowitz and Ladik.


title: Reconstructing Extreme Space Weather from Planet Hosting Stars

abstract: The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.


title: Study to achieve environment adaptive software

abstract: Recently, heterogeneous hardware such as GPU and FPGA is used in many systems and also IoT devices are increased repidly. However, to utilize heterogeneous hardware, the hurdles are high because of much technical skills. In order to break down such a situation, we think it is required in the future that application programmers only need to write logics to be processed, then software will adapt to the environments with heterogeneous hardware, to make it easy to utilize heterogeneous hardware and IoT devices. Therefore, in this paper, we propose environment adaptive software to operate an once written application with high performance by automatically converting the code and configuring setting so that we can utilize GPU, FPGA and IoT devices in the location to be deployed. We explain a proceesing flow and elemental technologies to achieve environment adaptive software.


title: Learning Super-resolution 3D Segmentation of Plant Root MRI Images from Few Examples

abstract: Analyzing plant roots is crucial to understand plant performance in different soil environments. While magnetic resonance imaging (MRI) can be used to obtain 3D images of plant roots, extracting the root structural model is challenging due to highly noisy soil environments and low-resolution of MRI images. To improve both contrast and resolution, we adapt the state-of-the-art method RefineNet for 3D segmentation of the plant root MRI images in super-resolution. The networks are trained from few manual segmentations that are augmented by geometric transformations, realistic noise, and other variabilities. The resulting segmentations contain most root structures, including branches not extracted by the human annotator.


title: An Extremal Property of the Hexagonal Lattice

abstract: We describe an extremal property of the hexagonal lattice $\Lambda \subset \mathbb{R}^2$. Let $p$ denote the circumcenter of its fundamental triangle (a so-called deep hole) and let $A_r$ denote the set of lattice points that are at distance $r$ from $p$ \begin{equation} A_r = \left\{ \lambda \in \Lambda: \| \lambda - p \| = r\right\}. \end{equation} If $\Gamma$ is a small perturbation of $\Lambda$ in the space of lattices with fixed density and $C_r$ denotes the set of points in $A_r$ shifted to the new lattice, then \begin{equation} \sum_{\mu \in C_r}{ \| p - \mu\|} - \sum_{\lambda \in A_r}{ \| p - \lambda\|} \gtrsim r \, |A_r| \, d(\Lambda, \Gamma)^2, \end{equation} where $d(\Lambda, \Gamma)$ denotes the distance between the lattices: the hexagonal lattice has the property that `far away points are closer than they are for nearby lattices'. This has implications in the calculus of variations: assume \begin{equation} g_{\Gamma}(z) = \sum_{\gamma \in \Gamma} f( \|z - \gamma \|) \quad \mbox{ satisfies } \quad \min_{z \in \mathbb{R}^2} g_{\Lambda}(z) = g_{\Lambda}(p). \end{equation} For a certain class of compactly supported functions $f$, the hexagonal lattice $\Lambda$ is then a strict local maximizer of \begin{equation} \max_{\Gamma} \min_{z \in \mathbb{R}^2} \sum_{\gamma \in \Gamma}{f( \|z - \gamma\| )}, \end{equation} where the maximum runs over all lattices of fixed density.


title: Numerical radius inequalities of operator matrices with applications

abstract: We present upper and lower bounds for the numerical radius of $2 \times 2$ operator matrices which improves on the existing bound for the same. As an application of the results obtained we give a better estimation for the zeros of a polynomial.


title: Orthogonality and Numerical radius inequalities of operator matrices

abstract: We completely characterize Birkhoff-James orthogonality with respect to numerical radius norm in the space of bounded linear operators on a complex Hilbert space. As applications of the results obtained, we estimate lower bounds of numerical radius for $n\times n$ operator matrices, which improve on and generalize existing lower bounds. We also obtain a better lower bound of numerical radius for an upper triangular operator matrix.


title: The Challenges of Observing, Calibrating, and Modeling Stellar Spectral Energy Distributions

abstract: While optical and quantum efficiency are on the rise, and spectrographs becoming massively multiplexed, measuring spectral energy distributions of astronomical sources with accuracy remains a challenge. In addition to atmospheric refraction, extinction, variability, and limited apertures of instrument entrance slits and optical fibers, accurate calibration is hampered by issues such as a limited choice of reliable standard stars. Modeling stellar spectral energy distributions has seen good progress, but some weaknesses survive, especially for late-type stars. This article provides an overview of these matters and discusses observation, calibration, and modeling strategies for future large spectroscopic surveys.


title: On the classification and false alarm of invalid prefixes in RPKI based BGP route origin validation

abstract: BGP is the default inter-domain routing protocol in today's Internet, but has serious security vulnerabilities\cite{murphy2005bgp}. One of them is (sub)prefix hijacking. IETF standardizes RPKI to validate the AS origin but RPKI has a lot of problems\cite{heilman2014consent}\cite{cooper2013risk}\cite{gilad2017we}\cite{gilad2017maxlength}, among which is potential false alarm. Although some previous work\cite{gilad2017we}\cite{heilman2014consent} points it out explicitly or implicitly, further measurement and analysis remain to be done. Our work measures and analyzes the invalid prefixes systematically. We first classify the invalid prefixes into six different types and then analyze their stability. We show that a large proportion of the invalid prefixes very likely result from traffic engineering, IP address transfer and failing to aggregate rather than real hijackings.


title: Dirac series for $E_{6(-14)}$

abstract: Up to equivalence, this paper classifies all the irreducible unitary representations with non-zero Dirac cohomology for the simple Lie group $E_{6(-14)}$, which is of Hermitian symmetric type. Each FS-scattered Dirac series of $E_{6(-14)}$ is realized as a composition factor of certain $A_{\mathfrak{q}}(\lambda)$ module. Along the way, we have also obtained all the fully supported irreducible unitary representations of $E_{6(-14)}$ with integral infinitesimal characters.


title: A KAM Theorem for Higher Dimensional Reversible Nonlinear Schr\"{o}dinger Equations

abstract: In the paper, we prove an abstract KAM (Kolmogorov-Arnold-Moser) theorem for infinite dimensional reversible systems. Using this KAM theorem, we obtain the existence and linear stability of quasi-periodic solutions for a class of reversible (non-Hamiltonian) coupled nonlinear Schr\"{o}dinger systems on $d-$torus $\mathbb{T}^d$.


title: Biquandle Module Invariants of Oriented Surface-Links

abstract: We define invariants of oriented surface-links by enhancing the biquandle counting invariant using \textit{biquandle modules}, algebraic structures defined in terms of biquandle actions on commutative rings analogous to Alexander biquandles. We show that bead colorings of marked graph diagrams are preserved by Yoshikawa moves and hence define enhancements of the biquandle counting invariant for surface links. We provide examples illustrating the computation of the invariant and demonstrate that these invariants are not determined by the first and second Alexander elementary ideals and characteristic polynomials.


title: Domain Generalization by Solving Jigsaw Puzzles

abstract: Human adaptability relies crucially on the ability to learn and merge knowledge both from supervised and unsupervised learning: the parents point out few important concepts, but then the children fill in the gaps on their own. This is particularly effective, because supervised learning can never be exhaustive and thus learning autonomously allows to discover invariances and regularities that help to generalize. In this paper we propose to apply a similar approach to the task of object recognition across domains: our model learns the semantic labels in a supervised fashion, and broadens its understanding of the data by learning from self-supervised signals how to solve a jigsaw puzzle on the same images. This secondary task helps the network to learn the concepts of spatial correlation while acting as a regularizer for the classification task. Multiple experiments on the PACS, VLCS, Office-Home and digits datasets confirm our intuition and show that this simple method outperforms previous domain generalization and adaptation solutions. An ablation study further illustrates the inner workings of our approach.


title: Reproducibility in $G_0W_0$ Calculations for Solids

abstract: Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ and $W$, and then perturbatively corrected by the resultant $GW$ self-energy. In practice, there are several ways to construct the $GW$ self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO$_2$, reported $GW$ fundamental gaps can vary by more than 1 eV. In this work, we address the convergence and key approximations in contemporary $G_0W_0$ calculations, including frequency-integration schemes and the treatment of the Coulomb divergence in the exact-exchange term. We study several systems,and compare three different $GW$ codes: BerkeleyGW, Abinit and Yambo. We demonstrate, for the first time, that the same quasiparticle energies for systems in the condensed phase can be obtained with different codes, and we provide a comprehensive assessment of implementations of the $GW$ approximation.


title: Fermionic Phonons: Exact Analytic Results and Quantum Statistical Mechanics for a One Dimensional Harmonic Crystal

abstract: Analytic expressions for the energy eigenvalues and eigenfunctions of a one-dimensional harmonic crystal are obtained. The average energy and density profiles are obtained numerically as a function of temperature. A surprisingly large number of energy levels (eg.\ 5,000 for 4 particles) are required for reliable results at even moderate temperatures.


title: Astro2020 science white paper: The gravitational wave view of massive black holes

abstract: Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.


title: An extension of Rohrlich's Theorem to the $j$-function

abstract: In this paper, we extend Rohrlich's Theorem on the integral of logarithms of meromorphic functions to compute the inner product between such functions and polynomials in the $j$-function. We then show that the generating function for these inner products with $j_n$ is part of a weight two modular object.


title: Notions of Centralized and Decentralized Opacity in Linear Systems

abstract: We formulate notions of opacity for cyberphysical systems modeled as discrete-time linear time-invariant systems. A set of secret states is $k$-ISO with respect to a set of nonsecret states if, starting from these sets at time $0$, the outputs at time $k$ are indistinguishable to an adversarial observer. Necessary and sufficient conditions to ensure that a secret specification is $k$-ISO are established in terms of sets of reachable states. We also show how to adapt techniques for computing under-approximations and over-approximations of the set of reachable states of dynamical systems in order to soundly approximate k-ISO. Further, we provide a condition for output controllability, if $k$-ISO holds, and show that the converse holds under an additional assumption. We extend the theory of opacity for single-adversary systems to the case of multiple adversaries and develop several notions of decentralized opacity. We study the following scenarios: i) the presence or lack of a centralized coordinator, and ii) the presence or absence of collusion among adversaries. In the case of colluding adversaries, we derive a condition for nonopacity that depends on the structure of the directed graph representing the communication between adversaries. Finally, we relax the condition that the outputs be indistinguishable and define a notion of $\epsilon$-opacity, and also provide an extension to the case of nonlinear systems.


title: Large Deviations for the Single Server Queue and the Reneging Paradox

abstract: For the M/M/1+M model at the law-of-large-numbers scale, the long run reneging count per unit time does not depend on the individual (i.e., per customer) reneging rate. This paradoxical statement has a simple proof. Less obvious is a large deviations analogue of this fact, stated as follows: The decay rate of the probability that the long run reneging count per unit time is atypically large or atypically small does not depend on the individual reneging rate. In this paper, the sample path large deviations principle for the model is proved and the rate function is computed. Next, large time asymptotics for the reneging rate are studied for the case when the arrival rate exceeds the service rate. The key ingredient is a calculus of variations analysis of the variational problem associated with atypical reneging. A characterization of the aforementioned decay rate, given explicitly in terms of the arrival and service rate parameters of the model, is provided yielding a precise mathematical description of this paradoxical behavior.


title: A Provably Communication-Efficient Asynchronous Distributed Inference Method for Convex and Nonconvex Problems

abstract: This paper proposes and analyzes a communication-efficient distributed optimization framework for general nonconvex nonsmooth signal processing and machine learning problems under an asynchronous protocol. At each iteration, worker machines compute gradients of a known empirical loss function using their own local data, and a master machine solves a related minimization problem to update the current estimate. We prove that for nonconvex nonsmooth problems, the proposed algorithm converges with a sublinear rate over the number of communication rounds, coinciding with the best theoretical rate that can be achieved for this class of problems. Linear convergence is established without any statistical assumptions of the local data for problems characterized by composite loss functions whose smooth parts are strongly convex. Extensive numerical experiments verify that the performance of the proposed approach indeed improves -- sometimes significantly -- over other state-of-the-art algorithms in terms of total communication efficiency.


title: Refinement of the classification of weak Fano threefolds with sextic del Pezzo fibrations

abstract: We refine the classification of weak Fano threefolds with sextic del Pezzo fibrations by considering the Hodge numbers of them. By the refined classification result, such threefolds are classified into 17 cases. The main result of this paper is to show that there exists an example for each of 17 cases, except for 2 cases.


title: Secure Control under Partial Observability with Temporal Logic Constraints

abstract: This paper studies the synthesis of control policies for an agent that has to satisfy a temporal logic specification in a partially observable environment, in the presence of an adversary. The interaction of the agent (defender) with the adversary is modeled as a partially observable stochastic game. The search for policies is limited to over the space of finite state controllers, which leads to a tractable approach to determine policies. The goal is to generate a defender policy to maximize satisfaction of a given temporal logic specification under any adversary policy. We relate the satisfaction of the specification in terms of reaching (a subset of) recurrent states of a Markov chain. We then present a procedure to determine a set of defender and adversary finite state controllers of given sizes that will satisfy the temporal logic specification. We illustrate our approach with an example.


title: Fast Interactive Object Annotation with Curve-GCN

abstract: Manually labeling objects by tracing their boundaries is a laborious process. In Polygon-RNN++ the authors proposed Polygon-RNN that produces polygonal annotations in a recurrent manner using a CNN-RNN architecture, allowing interactive correction via humans-in-the-loop. We propose a new framework that alleviates the sequential nature of Polygon-RNN, by predicting all vertices simultaneously using a Graph Convolutional Network (GCN). Our model is trained end-to-end. It supports object annotation by either polygons or splines, facilitating labeling efficiency for both line-based and curved objects. We show that Curve-GCN outperforms all existing approaches in automatic mode, including the powerful PSP-DeepLab and is significantly more efficient in interactive mode than Polygon-RNN++. Our model runs at 29.3ms in automatic, and 2.6ms in interactive mode, making it 10x and 100x faster than Polygon-RNN++.


title: Energy and performance evaluation of reactive, proactive, and hybrid routing protocols in wireless mesh network

abstract: This paper evaluates the energy consumption of well-known routing protocols, along with other metrics such as throughput, packet delivery ratio (PDR), and delay in different scenarios. We consider two other metrics in order to capture the efficiency of the energy consumption: e-throughput which is the ratio between the consumed energy and the throughput; and the e-PDR which is the ratio between the consumed energy and the PDR. We compare four routing protocols: AODV, OLSR, and HWMP in Reactive and Proactive modes. The number of nodes is varying between 25 and 81 nodes, with different mobility models. Simulations are conducted using NS3 and the parameters of a real network interface card. From the results, AODV presents the lowest energy consumption and a better e-Throughput. OLSR provides a better e-PDR in mobile scenarios. With a smaller e-PDR and e-Throughput, the proactive mode of HWMP is more energy efficient than the reactive mode.


title: A block tangential Lanczos method for model reduction of large-scale first and second order dynamical systems

abstract: In this paper, we present a new approach for model reduction of large scale first and second order dynamical systems with multiple inputs and multiple outputs (MIMO). This approach is based on the projection of the initial problem onto tangential Krylov subspaces to produce a simpler reduced-order model that approximates well the behavior of the original model. We present an algorithm named: Adaptive Block Tangential Lanczos-type (ABTL) algorithm. We give some algebraic properties and present some numerical experiences to show the effectiveness of the proposed algorithms.


title: Spherical Principal Component Analysis

abstract: Principal Component Analysis (PCA) is one of the most important methods to handle high dimensional data. However, most of the studies on PCA aim to minimize the loss after projection, which usually measures the Euclidean distance, though in some fields, angle distance is known to be more important and critical for analysis. In this paper, we propose a method by adding constraints on factors to unify the Euclidean distance and angle distance. However, due to the nonconvexity of the objective and constraints, the optimized solution is not easy to obtain. We propose an alternating linearized minimization method to solve it with provable convergence rate and guarantee. Experiments on synthetic data and real-world datasets have validated the effectiveness of our method and demonstrated its advantages over state-of-art clustering methods.


title: Controlled enantioselective orientation of chiral molecules with an optical centrifuge

abstract: We report on the first experimental demonstration of enantioselective rotational control of chiral molecules with a laser field. In our experiments, two enantiomers of propylene oxide are brought to accelerated unidirectional rotation by means of an optical centrifuge. Using Coulomb explosion imaging, we show that the centrifuged molecules acquire preferential orientation perpendicular to the plane of rotation, and that the direction of this orientation depends on the relative handedness of the enantiomer and the rotating centrifuge field. The observed effect is in agreement with theoretical predictions and is reproduced in numerical simulations of the centrifuge excitation followed by Coulomb explosion of the centrifuged molecules. The demonstrated technique opens new avenues in optical enantioselective control of chiral molecules with a plethora of potential applications in differentiation, separation and purification of chiral mixtures.


title: Ontology Based Global and Collective Motion Patterns for Event Classification in Basketball Videos

abstract: In multi-person videos, especially team sport videos, a semantic event is usually represented as a confrontation between two teams of players, which can be represented as collective motion. In broadcast basketball videos, specific camera motions are used to present specific events. Therefore, a semantic event in broadcast basketball videos is closely related to both the global motion (camera motion) and the collective motion. A semantic event in basketball videos can be generally divided into three stages: pre-event, event occurrence (event-occ), and post-event. In this paper, we propose an ontology-based global and collective motion pattern (On_GCMP) algorithm for basketball event classification. First, a two-stage GCMP based event classification scheme is proposed. The GCMP is extracted using optical flow. The two-stage scheme progressively combines a five-class event classification algorithm on event-occs and a two-class event classification algorithm on pre-events. Both algorithms utilize sequential convolutional neural networks (CNNs) and long short-term memory (LSTM) networks to extract the spatial and temporal features of GCMP for event classification. Second, we utilize post-event segments to predict success/failure using deep features of images in the video frames (RGB_DF_VF) based algorithms. Finally the event classification results and success/failure classification results are integrated to obtain the final results. To evaluate the proposed scheme, we collected a new dataset called NCAA+, which is automatically obtained from the NCAA dataset by extending the fixed length of video clips forward and backward of the corresponding semantic events. The experimental results demonstrate that the proposed scheme achieves the mean average precision of 58.10% on NCAA+. It is higher by 6.50% than state-of-the-art on NCAA.


title: Can the dynamical Lamb effect be observed in a superconducting circuit?

abstract: The dynamical Lamb effect is predicted to arise in superconducting circuits when the coupling of a superconducting qubit with a resonator is periodically switched "on" and "off" nonadiabatically. We show that by using a superconducting circuit which allows to switch between longitudinal and transverse coupling of a qubit to a resonator, it is possible of to observe the dynamical Lamb effect. {The switching between longitudinal and transverse coupling can be achieved by modulating the magnetic flux through the circuit loops.} By solving the Schr\"{o}dinger equation for a qubit coupled to a resonator, we calculate the time evolution of the probability of excitation of the qubit and the creation of $n$ photons in the resonator due to the dynamical Lamb effect. The probability is maximum when the coupling is periodically switched between longitudinal and transverse using a square-wave or sinusoidal modulation of the magnetic flux with frequency equal to the sum of the average qubit and photon transition frequencies.


title: A self-consistent method to estimate the rate of compact binary coalescences with a Poisson mixture model

abstract: The recently published GWTC-1 - a journal article summarizing the search for gravitational waves (GWs) from coalescing compact binaries in data produced by the LIGO-Virgo network of ground-based detectors during their first and second observing runs - quoted estimates for the rates of binary neutron star, neutron star black hole binary, and binary black hole mergers, as well as assigned probabilities of astrophysical origin for various significant and marginal GW candidate events. In this paper, we delineate the formalism used to compute these rates and probabilities, which assumes that triggers above a low ranking statistic threshold, whether of terrestrial or astrophysical origin, occur as independent Poisson processes. In particular, we include an arbitrary number of astrophysical categories by redistributing, via mass-based template weighting, the foreground probabilities of candidate events, across source classes. We evaluate this formalism on synthetic GW data, and demonstrate that this method works well for the kind of GW signals observed during the first and second observing runs.


title: Classification of irreducible Harish-Chandra modules over gap-$p$ Virasoro algebras

abstract: We prove that any irreducible Harish-Chandra modules for a class of Lie algebras, which we call gap-$p$ Virasoro algebras, must be a highest weight module, a lowest weight module, or a module of intermediate series.These algebras are closely related to the Heisenberg-Virasoro algebra and the algebra of derivations over a quantum torus. They also contain subalgebras which are isomorphic to the Virasoro algebra $Vir$, but graded by $p\mathbb Z$(unlike $Vir$ by $\mathbb Z$).


title: Comments on $T\bar T$, $J\bar{T}$ and String Theory

abstract: These notes contain some aspects of the holographic duality between the perturbative superstring on current-current deformations of $AdS_3\times S^1\times{\cal N}$ and single-trace $T\bar T$ and $J\bar T$ deformed $CFT_2$.


title: Multiplicity moments of shifted Gompertz distribution in $e^{+}e^{-}$, $p\bar{p}$ and $pp$ collisions at high energies

abstract: In continuation of our earlier work, in which we analysed the charged particle multiplicities in leptonic and hadronic interactions at different center of mass energies in full phase space as well as in restricted phase space with the shifted Gompertz distribution, a detailed analysis of the normalized and factorial moments is reported here. A two-component model in which probability distribution function is obtained from the superposition of two shifted Gompertz distributions introduced in our earlier work has also been used for the analysis. This is the first analysis of the moments with the shifted Gompertz distribution. Analysis has also been done to predict the moments of multiplicity distribution for the electron-positron collisions at c.m. energy of 500 GeV at a future Collider.