title: Unsupervised Domain Adaptation for Cross-Subject Few-Shot Neurological Symptom Detection

abstract: Modern machine learning tools have shown promise in detecting symptoms of neurological disorders. However, current approaches typically train a unique classifier for each subject. This subject-specific training scheme requires long labeled recordings from each patient, thus failing to detect symptoms in new patients with limited recordings. This paper introduces an unsupervised domain adaptation approach based on adversarial networks to enable few-shot, cross-subject epileptic seizure detection. Using adversarial learning, features from multiple patients were encoded into a subject-invariant space and a discriminative model was trained on subject-invariant features to make predictions. We evaluated this approach on the intracranial EEG (iEEG) recordings from 9 patients with epilepsy. Our approach enabled cross-subject seizure detection with a 9.4\% improvement in 1-shot classification accuracy compared to the conventional subject-specific scheme.


title: Graphs with constant adjacency dimension

abstract: For a set W of vertices and a vertex v in a graph G, the k-vector r2(v|W) = (aG(v,w1),...,aG(v,wk)) is the adjacency representation of v with respect to W, where W = {w1,...,wk} and aG(x,y) is the minimum of 2 and the distance between the vertices x and y. The set W is an adjacency resolving set for G if distinct vertices of G have distinct adjacency representations with respect to W. The minimum cardinality of an adjacency resolving set for G is its adjacency dimension. It is clear that the adjacency dimension of an n-vertex graph G is between 1 and n-1. The graphs with adjacency dimension 1 and n-1 are known. All graphs with adjacency dimension 2, and all n-vertex graphs with adjacency dimension n-2 are studied in this paper. In terms of the diameter and order of G, a sharp upper bound is found for adjacency dimension of G. Also, a sharp lower bound for adjacency dimension of G is obtained in terms of order of G. Using these two bounds, all graphs with adjacency dimension 2, and all n-vertex graphs with adjacency dimension n-2 are characterized.


title: Active learning based generative design for the discovery of wide bandgap materials

abstract: Active learning has been increasingly applied to screening functional materials from existing materials databases with desired properties. However, the number of known materials deposited in the popular materials databases such as ICSD and Materials Project is extremely limited and consists of just a tiny portion of the vast chemical design space. Herein we present an active generative inverse design method that combines active learning with a deep variational autoencoder neural network and a generative adversarial deep neural network model to discover new materials with a target property in the whole chemical design space. The application of this method has allowed us to discover new thermodynamically stable materials with high band gap (SrYF$_5$) and semiconductors with specified band gap ranges (SrClF$_3$, CaClF$_5$, YCl$_3$, SrC$_2$F$_3$, AlSCl, As$_2$O$_3$), all of which are verified by the first principle DFT calculations. Our experiments show that while active learning itself may sample chemically infeasible candidates, these samples help to train effective screening models for filtering out materials with desired properties from the hypothetical materials created by the generative model. The experiments show the effectiveness of our active generative inverse design approach.


title: Multifractality Breaking from Bounded Random Measures

abstract: Multifractal systems usually have singularity spectra defined on bounded sets of H\"older exponents. As a consequence, their associated multifractal scaling exponents are expected to depend linearly upon statistical moment orders at high enough orders -- a phenomenon referred to as the {\it{linearization effect}}. Motivated by general ideas taken from models of turbulent intermittency and focusing on the case of two-dimensional systems, we investigate the issue within the framework of Gaussian multiplicative chaos. As verified by means of Monte Carlo simulations, it turns out that the linearization effect can be accounted for by Liouville-like random measures defined in terms of upper-bounded scalar fields. The coarse-grained statistical properties of Gaussian multiplicative chaos are furthermore found to be preserved in the linear regime of the scaling exponents. As a related application, we look at the problem of turbulent circulation statistics, and obtain a remarkably accurate evaluation of circulation statistical moments, recently determined with the help of massive numerical simulations.


title: Geometric effects on the electronic structure of curved nanotubes and curved graphene: the case of the helix, catenary, helicoid, and catenoid

abstract: Since electrons in a ballistic regime perceive a carbon nanotube or a graphene layer structure as a continuous medium, we can use the study of the quantum dynamics of one electron constrained to a curve or surface to obtain a qualitative description of the conduction electrons' behavior. The confinement process of a quantum particle to a curve or surface leads us, in the so-called "confining potential formalism" (CPF), to a geometry-induced potential (GIP) in the effective Schr\"odinger equation. With these considerations, this work aims to study in detail the consequences of constraining a quantum particle to a helix, catenary, helicoid, or catenoid, exploring the relations between these curves and surfaces using differential geometry. Initially, we use the variational method to estimate the energy of the particle in its ground state, and thus, we obtain better approximations with the use of the confluent Heun function through numerical calculations. Thus, we conclude that a quantum particle constrained to an infinite helix has its angular momentum quantized due to the geometry of the curve, while in the cases of the catenary, helicoid, and catenoid the particle can be found either in a single bound state or in excited states which constitute a continuous energy band. Additionally, we propose measurements of physical observables capable of discriminating the topologies of the studied surfaces, in the context of topological metrology.


title: Rosenbluth separation of $ \boldsymbol {d\sigma_{L}/dt}$ and $\boldsymbol {d\sigma_{T}/dt}$ in $\boldsymbol \pi^{\bf 0}$ deeply virtual electroproduction from the proton

abstract: We report on a Rosenbluth separation using previously published data by the CLAS collaboration in Hall B, Jefferson Lab for exclusive $\pi^{0}$ deeply virtual electroproduction (DVEP) from the proton at a mean $Q^{2}$ of $\approx$ 2 (GeV/c)$^{2}$. The central question we address is the applicability of factorization in $\pi^0$ DVEP at these kinematics. The results of our Rosenbluth separation clearly demonstrate the dominance of the longitudinal contribution to the cross section. The extracted longitudinal and transverse contributions are in agreement with previous data from Hall A at Jefferson Lab, but over a much wider $-t$ range (0.12 - 1.8 (GeV/c)$^{2}$). The measured dominance of the longitudinal contribution at $Q^{2} \approx$ 2 (GeV/c)$^{2}$ is consistent with the expectation of the handbag factorization theorem. We find that $\sigma_L(t) \sim 1/(-t)$ for $-t >$ 0.5 (GeV/c)$^2$. Determination of both longitudinal and transverse contributions to the deeply virtual $\pi^{0}$ electroproduction cross section allows extraction of additional GPDs.


title: Tall Towers on the Moon

abstract: The lunar South pole likely contains significant amounts of water in the permanently shadowed craters there. Extracting this water for life support at a lunar base or to make rocket fuel would take large amounts of power, of order Gigawatts. A natural place to obtain this power are the "Peaks of Eternal Light", that lie a few kilometers away on the crater rims and ridges above the permanently shadowed craters. The amount of solar power that could be captured depends on how tall a tower can be built to support the photovoltaic panels. The low gravity, lack of atmosphere, and quiet seismic environment of the Moon suggests that towers could be built much taller than on Earth. Here we look at the limits to building tall concrete towers on the Moon. We choose concrete as the capital cost of transporting large masses of iron or carbon fiber to the Moon is presently so expensive that profitable operation of a power plant is unlikely. Concrete instead can be manufactured in situ from the lunar regolith. We find that, with minimum wall thicknesses (20 cm), towers up to several kilometers tall are stable. The mass of concrete needed, however, grows rapidly with height, from $\sim$ 760 mt at 1 km to $\sim$ 4,100 mt at 2 km to $\sim 10^5$ mt at 7 km and $\sim 10^6$ mt at 17 km.


title: A novel technique for real-time ion identification and energy measurement for in situ space instrumentation

abstract: The AGILE (Advanced enerGetic Ion eLectron tElescope) project focuses on the development of a compact low-cost space-based instrument to measure the intensities of charged particles and ions in space. Using multiple layers of fast silicon sensors and custom front-end electronics, the instrument is designed for real-time particle identification of a large variety of elements from H to Fe and spanning energies from 1 to 100 MeV per nucleon. The robust method proposed in this work uses key defining features of electronic signals generated by charged particles (ions) traveling through silicon layers to reliably identify and characterize particles in situ. AGILE will use this real-time pulse shape discrimination technique for the first time in space based instrumentation.


title: Tuning the Reactivity of Nanoenergetic Gas Generators Based on Bismuth and Iodine oxidizers

abstract: There is a growing interest on novel energetic materials called Nanoenergetic Gas-Generators (NGGs) which are potential alternatives to traditional energetic materials including pyrotechnics, propellants, primers and solid rocket fuels. NGGs are formulations that utilize metal powders as a fuel and oxides or hydroxides as oxidizers that can rapidly release large amount of heat and gaseous products to generate shock waves. The heat and pressure discharge, impact sensitivity, long term stability and other critical properties depend on the particle size and shape, as well as assembling procedure and intermixing degree between the components. The extremely high energy density and the ability to tune the dynamic properties of the energetic system makes NGGs ideal candidates to dilute or replace traditional energetic materials for emerging applications. In terms of energy density, performance and controllability of dynamic properties, the energetic materials based on bismuth and iodine compounds are exceptional among the NGGs. The thermodynamic calculations and experimental study confirm that NGGs based on iodine and bismuth compounds mixed with aluminum nanoparticles are the most powerful formulations to date and can be used potentially in microthrusters technology with high thrust-to-weight ratio with controlled combustion and exhaust velocity for space applications. The resulting nano thermites generated significant value of pressure discharge up to 14.8 kPa m3/g. They can also be integrated with carbon nanotubes to form laminar composite yarns with high power actuation of up to 4700 W/kg, or be used in other emerging applications such as biocidal agents to effectively destroy harmful bacteria in seconds, with 22 mg/m2 minimal content over infected area.


title: The impact of kinetic inductance on the critical current oscillations of nanobridge SQUIDs

abstract: In this work, we study the current phase relation ($C\Phi R$) of lithographically fabricated molybdenum germanium (MoGe) nanobridges, which is intimately linked to the nanobridge kinetic inductance. We do this by imbedding the nanobridges in a SQUID. We observe that for temperatures far below $T_c$, the $C\Phi R$ is linear as long as the condensate is not weakened by the presence of supercurrent. We demonstrate lithographic control over the nanobridge kinetic inductance, which scales with the nanobridge aspect ratio. This allows to tune the SQUID $I_c(B)$ characteristic. The SQUID properties that can be controlled in this way include the SQUID sensitivity and the positions of the critical current maxima. These observations can be of use for the design and operation of future superconducting devices such as magnetic memories or flux qubits.


title: Learning Human-like Hand Reaching for Human-Robot Handshaking

abstract: One of the first and foremost non-verbal interactions that humans perform is a handshake. It has an impact on first impressions as touch can convey complex emotions. This makes handshaking an important skill for the repertoire of a social robot. In this paper, we present a novel framework for learning reaching behaviours for human-robot handshaking behaviours for humanoid robots solely using third-person human-human interaction data. This is especially useful for non-backdrivable robots that cannot be taught by demonstrations via kinesthetic teaching. Our approach can be easily executed on different humanoid robots. This removes the need for re-training, which is especially tedious when training with human-interaction partners. We show this by applying the learnt behaviours on two different humanoid robots with similar degrees of freedom but different shapes and control limits.


title: Positive scalar curvature on $\mathbf{Pin}^\pm$- and $\mathbf{Spin}^c$-manifolds

abstract: It is well-known that spin structures and Dirac operators play a crucial role in the study of positive scalar curvature metrics (psc-metrics) on compact manifolds. Here we consider a class of non-spin manifolds with "almost spin" structure, namely those with spin$^c$ or pin$^\pm$-structures. It turns out that in those cases (under natural assumptions on such a manifold $M$), the index of a relevant Dirac operator completely controls existence of a psc-metric which is $S^1$- or $C_2$-invariant near a "special submanifold" $B$ of $M$. This submanifold $B\subset M$ is dual to the complex (respectively, real) line bundle $L$ which determines the spin$^c$ or pin$^\pm$ structure on $M$. We also show that these manifold pairs $(M,B)$ can be interpreted as "manifolds with fibered singularities" equipped with "well-adapted psc-metrics". This survey is based on our recent work as well as on our joint work with Paolo Piazza.


title: 2-switch transition on unicyclic graphs and pseudoforest

abstract: In the present work we prove that given any two unicycle graphs (pseudoforests) that share the same degree sequence there is a finite sequence of 2-switches transforming one into the other such that all the graphs in the sequence are also unicyclic graphs (pseudoforests).


title: Monte Carlo-based estimation of patient absorbed dose in 99mTc-DMSA, -MAG3, and -DTPA SPECT imaging using the University of Florida (UF) phantoms

abstract: SPECT imaging is one of the common modalities used to examine kidney disease. The radiopharmaceuticals such as 99mTc-MAG3, 99mTc-DTPA, and 99mTc-DMSA are commonly employed for pediatric patients. It is generally accepted that children organs are more sensitive to radiation than adults, due to their growth rate. Therefore, evaluating the absorbed dose in children to avoid irrecoverable damage is highly crucial. In this work, absorbed dose by different organs of children within the SPECT imaging for the above-mentioned radiopharmaceuticals were estimated through the use of Monte Carlo simulation for patients at the ages of 4, 8, 11, and 14-years old. The patients absorbed dose was estimated based on the specific radiopharmaceutical uptake patterns and the age-related human body models. To this end, we used the University of Florida (UF) voxel-wise baby phantoms and the Monte Carlo simulation for the pediatric patients at 4, 8, 11, and 14 years of age. The accuracy of the dosimetric outcomes was evaluated versus the ICRP 128 data. The highest and lowest doses absorbed by the kidneys were 0.48 mGy-MBq and 0.0042 mGy-MBq obtained for 4-year-old and 14-year-old children using 99mTc-DMSA and 99mTc-MAG3 radiopharmaceuticals, respectively. Moreover, in these cases, the relative errors of simulation results compared to the reference ICRP data were 12.1% and 0.95%, respectively. The simulation results were in good agreement with the ICRP 128 data. The results showed that the highest absorbed dose was by kidneys and when 99mTc-DMSA was used. The highest and lowest absorbed dose in the organs occurs when 99mTc-DMSA and 99mTc-MAG3 (in normal renal function) are used, respectively. On the other hand, the results indicated that the dose received by the organs decreases with age.


title: The Mechanism of Scale-Invariance

abstract: A system is invariant with respect to an input transformation if we can transform any dynamic input by this function and obtain the same output dynamics after adjusting the initial conditions appropriately. Often, the set of all such input transformations forms a Lie group, the most prominent examples being scale-invariant ($u\mapsto e^pu$, $p\in\mathbb{R}$) and translational-invariant ($u\mapsto pu$) systems, the latter comprising linear systems with transfer function zeros at the origin. Here, we derive a necessary and sufficient normal form for invariant systems and, by analyzing this normal form, provide a complete characterization of the mechanism by which invariance can be achieved. In this normal form, all invariant systems (i) estimate the applied input transformation by means of an integral feedback, and (ii) then apply the inverse of this estimate to the input before processing it in any other way. We demonstrate our results based on three examples: a scale-invariant "feed-forward loop", a bistable switch, and a system resembling the core of the mammalian circadian network.


title: eQuilibrator 3.0 -- a platform for the estimation of thermodynamic constants

abstract: eQuilibrator (equilibrator.weizmann.ac.il) is a calculator for biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts ${\sim}1000$ distinct monthly users, its design could not accommodate larger compound databases and it lacked an application programming interface (API) for integration in other tools developed by the systems biology community. Here, we report a new python-based package for eQuilibrator, that comes with many new features such as a 50-fold larger compound database, the ability to add novel compound structures, improvements in speed and memory use, and correction for Mg2+ ion concentrations. Moreover, it adds the ability to compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modeling. We foresee that these improvements will make thermodynamic modeling more accessible and facilitate the integration of eQuilibrator into other software platforms.


title: On Surrogate Learning for Linear Stability Assessment of Navier-Stokes Equations with Stochastic Viscosity

abstract: We study linear stability of solutions to the Navier\textendash Stokes equations with stochastic viscosity. Specifically, we assume that the viscosity is given in the form of a~stochastic expansion. Stability analysis requires a solution of the steady-state Navier-Stokes equation and then leads to a generalized eigenvalue problem, from which we wish to characterize the real part of the rightmost eigenvalue. While this can be achieved by Monte Carlo simulation, due to its computational cost we study three surrogates based on generalized polynomial chaos, Gaussian process regression and a shallow neural network. The results of linear stability analysis assessment obtained by the surrogates are compared to that of Monte Carlo simulation using a set of numerical experiments.


title: Scaling up Mean Field Games with Online Mirror Descent

abstract: We address scaling up equilibrium computation in Mean Field Games (MFGs) using Online Mirror Descent (OMD). We show that continuous-time OMD provably converges to a Nash equilibrium under a natural and well-motivated set of monotonicity assumptions. This theoretical result nicely extends to multi-population games and to settings involving common noise. A thorough experimental investigation on various single and multi-population MFGs shows that OMD outperforms traditional algorithms such as Fictitious Play (FP). We empirically show that OMD scales up and converges significantly faster than FP by solving, for the first time to our knowledge, examples of MFGs with hundreds of billions states. This study establishes the state-of-the-art for learning in large-scale multi-agent and multi-population games.


title: The Phantom Alignment Strength Conjecture: Practical use of graph matching alignment strength to indicate a meaningful graph match

abstract: The alignment strength of a graph matching is a quantity that gives the practitioner a measure of the correlation of the two graphs, and it can also give the practitioner a sense for whether the graph matching algorithm found the true matching. Unfortunately, when a graph matching algorithm fails to find the truth because of weak signal, there may be "phantom alignment strength" from meaningless matchings that, by random noise, have fewer disagreements than average (sometimes substantially fewer); this alignment strength may give the misleading appearance of significance. A practitioner needs to know what level of alignment strength may be phantom alignment strength and what level indicates that the graph matching algorithm obtained the true matching and is a meaningful measure of the graph correlation. The {\it Phantom Alignment Strength Conjecture} introduced here provides a principled and practical means to approach this issue. We provide empirical evidence for the conjecture, and explore its consequences.


title: Rates of multivariate normal approximation for statistics in geometric probability

abstract: We employ stabilization methods and second order Poincar\'e inequalities to establish rates of multivariate normal convergence for a large class of vectors $(H_s^{(1)},...,H_s^{(m)})$, $s \geq 1$, of statistics of marked Poisson processes on $\mathbb{R}^d$, $d \geq 2$, as the intensity parameter $s$ tends to infinity. Our results are applicable whenever the constituent functionals $H_s^{(i)}$, $i\in\{1,...,m\}$, are expressible as sums of exponentially stabilizing score functions satisfying a moment condition. The rates are for the $d_2$-, $d_3$-, and $d_{convex}$-distances. When we compare with a centered Gaussian random vector, whose covariance matrix is given by the asymptotic covariances, the rates are in general unimprovable and are governed by the rate of convergence of $s^{-1} {\rm Cov}( H_s^{(i)}, H_s^{(j)})$, $i,j\in\{1,...,m\}$, to the limiting covariance, shown to be of order $s^{-1/d}$. We use the general results to deduce rates of multivariate normal convergence for statistics arising in random graphs and topological data analysis as well as for multivariate statistics used to test equality of distributions. Some of our results hold for stabilizing functionals of Poisson input on suitable metric spaces.


title: Analysis of critical parameters for nonrelativistic models in symmetric nuclear matter

abstract: In this work we have analyzed several features of symmetric nuclear matter (SNM) at finite temperature described by different zero- and finite-range nonrelativistic families of models, namely, Skyrme, Gogny, Momentum-dependent interaction (MDI), Michigan three-range Yukawa (M3Y) and Simple Effective Interaction (SEI). We have calculated the critical parameters (CP) associated to the liquid-gas phase coexistence for nuclear matter from these parametrizations and show that they are in agreement with their experimental and theoretical values obtained in the literature. Our study also points out to a strong evidence of universality presented by the hadronic models, namely, model independence in the gaseous phase and distinguishability among different interactions in the liquid phase. We have performed a correlation study among different CP and SNM properties. Such studies involving different finite range interactions are scarce in literature. The analyzed models show an overall increasing trend of the critical temperature as a function of critical pressure.


title: Multi Split Conformal Prediction

abstract: Split conformal prediction is a computationally efficient method for performing distribution-free predictive inference in regression. It involves, however, a one-time random split of the data, and the result depends on the particular split. To address this problem, we propose multi split conformal prediction, a simple method based on Markov's inequality to aggregate single split conformal prediction intervals across multiple splits.


title: The influence of the physical coefficients of a Bresse system with one singular local viscous damping in the longitudinal displacement on its stabilization

abstract: In this paper, we investigate the stabilization of a linear Bresse system with one singular local frictional damping acting in the longitudinal displacement, under fully Dirichlet boundary conditions. First, we prove the strong stability of our system. Next, using a frequency domain approach combined with the multiplier method, we establish the exponential stability of the solution if and only if the three waves have the same speed of propagation. On the contrary, we prove that the energy of our system decays polynomially with rates $t^{-1}$ or $t^{-\frac{1}{2}}$.


title: A Hierarchical Spike-and-Slab Model for Pan-Cancer Survival Using Pan-Omic Data

abstract: Pan-omics, pan-cancer analysis has advanced our understanding of the molecular heterogeneity of cancer, expanding what was known from single-cancer or single-omics studies. However, pan-cancer, pan-omics analyses have been limited in their ability to use information from multiple sources of data (e.g., omics platforms) and multiple sample sets (e.g., cancer types) to predict important clinical outcomes, like overall survival. We address the issue of prediction across multiple high-dimensional sources of data and multiple sample sets by using exploratory results from BIDIFAC+, a method for integrative dimension reduction of bidimensionally-linked matrices, in a predictive model. We apply a Bayesian hierarchical model that performs variable selection using spike-and-slab priors which are modified to allow for the borrowing of information across clustered data. This method is used to predict overall patient survival from the Cancer Genome Atlas (TCGA) using data from 29 cancer types and 4 omics sources. Our model selected patterns of variation identified by BIDIFAC+ that differentiate clinical tumor subtypes with markedly different survival outcomes. We also use simulations to evaluate the performance of the modified spike-and-slab prior in terms of its variable selection accuracy and prediction accuracy under different underlying data-generating frameworks. Software and code used for our analysis can be found at https://github.com/sarahsamorodnitsky/HierarchicalSS_PanCanPanOmics/ .


title: Physical-Layer Security via Distributed Beamforming in the Presence of Adversaries with Unknown Locations

abstract: We study the problem of securely communicating a sequence of information bits with a client in the presence of multiple adversaries at unknown locations in the environment. We assume that the client and the adversaries are located in the far-field region, and all possible directions for each adversary can be expressed as a continuous interval of directions. In such a setting, we develop a periodic transmission strategy, i.e., a sequence of joint beamforming gain and artificial noise pairs, that prevents the adversaries from decreasing their uncertainty on the information sequence by eavesdropping on the transmission. We formulate a series of nonconvex semi-infinite optimization problems to synthesize the transmission strategy. We show that the semi-definite program (SDP) relaxations of these nonconvex problems are exact under an efficiently verifiable sufficient condition. We approximate the SDP relaxations, which are subject to infinitely many constraints, by randomly sampling a finite subset of the constraints and establish the probability with which optimal solutions to the obtained finite SDPs and the semi-infinite SDPs coincide. We demonstrate with numerical simulations that the proposed periodic strategy can ensure the security of communication in scenarios in which all stationary strategies fail to guarantee security.


title: On the Subbagging Estimation for Massive Data

abstract: This article introduces subbagging (subsample aggregating) estimation approaches for big data analysis with memory constraints of computers. Specifically, for the whole dataset with size $N$, $m_N$ subsamples are randomly drawn, and each subsample with a subsample size $k_N\ll N$ to meet the memory constraint is sampled uniformly without replacement. Aggregating the estimators of $m_N$ subsamples can lead to subbagging estimation. To analyze the theoretical properties of the subbagging estimator, we adapt the incomplete $U$-statistics theory with an infinite order kernel to allow overlapping drawn subsamples in the sampling procedure. Utilizing this novel theoretical framework, we demonstrate that via a proper hyperparameter selection of $k_N$ and $m_N$, the subbagging estimator can achieve $\sqrt{N}$-consistency and asymptotic normality under the condition $(k_Nm_N)/N\to \alpha \in (0,\infty]$. Compared to the full sample estimator, we theoretically show that the $\sqrt{N}$-consistent subbagging estimator has an inflation rate of $1/\alpha$ in its asymptotic variance. Simulation experiments are presented to demonstrate the finite sample performances. An American airline dataset is analyzed to illustrate that the subbagging estimate is numerically close to the full sample estimate, and can be computationally fast under the memory constraint.


title: A weighted POD-reduction approach for parametrized PDE-constrained Optimal Control Problems with random inputs and applications to environmental sciences

abstract: Reduced basis approximations of Optimal Control Problems (OCPs) governed by steady partial differential equations (PDEs) with random parametric inputs are analyzed and constructed. Such approximations are based on a Reduced Order Model, which in this work is constructed using the method of weighted Proper Orthogonal Decomposition. This Reduced Order Model then is used to efficiently compute the reduced basis approximation for any outcome of the random parameter. We demonstrate that such OCPs are well-posed by applying the adjoint approach, which also works in the presence of admissibility constraints and in the case of non linear-quadratic OCPs, and thus is more general than the conventional Lagrangian approach. We also show that a step in the construction of these Reduced Order Models, known as the aggregation step, is not fundamental and can in principle be skipped for noncoercive problems, leading to a cheaper online phase. Numerical applications in three scenarios from environmental science are considered, in which the governing PDE is steady and the control is distributed. Various parameter distributions are taken, and several implementations of the weighted Proper Orthogonal Decomposition are compared by choosing different quadrature rules.


title: Setup of high resolution thermal expansion measurements in closed cycle cryostats using capacitive dilatometers

abstract: We present high resolution thermal expansion measurement data obtained with high relative sensitivity of $\Delta$L$/$L = 10$^{-9}$ and accuracy of $\pm$2$\%$ using closed cycle refrigerators employing two different dilatometers. The data quality is in excellent agreement with those obtained using wet liquid helium based systems, demonstrating great technological possibilities for future thermal expansion measurements in view of the depleting resource of liquid helium. The cryogenic environment was achieved using two different cryostats that use pulse tube and Gifford-Mcmahon coolers as the cryocoolers. Both the dilatometers employ a spring movement for achieving the parallel movement of the capacitor plates. $Dilatometer \#1$ was built in-house based on a published design while $dilatometer \#2$ was obtained commercially. Cell calibration for $dilatometer \#1$ was done using copper and minimal deviation of the cell effect from the published values were found. Linear thermal expansion coefficient $\alpha$ obtained using both dilatometers was evaluated using two different techniques, namely, numerical differentiation and derivative of a polynomial fit. The resultant $\alpha$ obtained for metals silver and aluminium showed excellent match with published values obtained on systems using wet cryostats. Finite element method simulations were performed for understanding the spring movement in each dilatometer using which the effect of different forces$/$pressures on the displacement of the spring was studied. Finally, we report thermal expansion measurements done on single crystals of two high temperature superconductors YBa$_2$Cu$_{3-x}$Al$_x$O$_{6+\delta}$ and Bi$_2$Sr$_2$CaCu$_2$O$_{8+x}$ along the c-axis and found very good match with published data obtained using wet liquid helium based cryostats.


title: TransCT: Dual-path Transformer for Low Dose Computed Tomography

abstract: Low dose computed tomography (LDCT) has attracted more and more attention in routine clinical diagnosis assessment, therapy planning, etc., which can reduce the dose of X-ray radiation to patients. However, the noise caused by low X-ray exposure degrades the CT image quality and then affects clinical diagnosis accuracy. In this paper, we train a transformer-based neural network to enhance the final CT image quality. To be specific, we first decompose the noisy LDCT image into two parts: high-frequency (HF) and low-frequency (LF) compositions. Then, we extract content features (X_{L_c}) and latent texture features (X_{L_t}) from the LF part, as well as HF embeddings (X_{H_f}) from the HF part. Further, we feed X_{L_t} and X_{H_f} into a modified transformer with three encoders and decoders to obtain well-refined HF texture features. After that, we combine these well-refined HF texture features with the pre-extracted X_{L_c} to encourage the restoration of high-quality LDCT images with the assistance of piecewise reconstruction. Extensive experiments on Mayo LDCT dataset show that our method produces superior results and outperforms other methods.


title: The impact of homophily on digital proximity tracing

abstract: We study how homophily of human physical interactions affects the efficacy of digital proximity tracing. Analytical results show a non monotonous dependence of the reproduction number with respect to the mixing rate between individuals that adopt the contact tracing app and the ones that do not. Furthermore, we find regimes in which the attack rate has local optima, minima or monotonously varies with the mixing rate. We corroborate our findings with Monte Carlo simulations on a primary-school network. This study provides a mathematical basis to better understand how homophily in health behavior shapes the dynamics of epidemics.


title: Counting Representations of Quivers Respecting Nilpotent Relations over Finite Fields

abstract: This paper presents analogous results of Hua [7][8] on numbers of representations of quivers over finite fields which respect nilpotent relations under certain assumptions. A closed formula which counts isomorphism classes of absolutely indecomposable representations with given dimension vectors is given and a $q$-deformation of Weyl-Kac denominator identity is established. In principle, if the numbers of representations are known, then the numbers of isomorphism classes of absolutely indecomposable representations are known.


title: Detection of Malicious Android Applications: Classical Machine Learning vs. Deep Neural Network Integrated with Clustering

abstract: Today anti-malware community is facing challenges due to the ever-increasing sophistication and volume of malware attacks developed by adversaries. Traditional malware detection mechanisms are not able to cope-up with next-generation malware attacks. Therefore in this paper, we propose effective and efficient Android malware detection models based on machine learning and deep learning integrated with clustering. We performed a comprehensive study of different feature reduction, classification and clustering algorithms over various performance metrics to construct the Android malware detection models. Our experimental results show that malware detection models developed using Random Forest eclipsed deep neural network and other classifiers on the majority of performance metrics. The baseline Random Forest model without any feature reduction achieved the highest AUC of 99.4%. Also, the segregating of vector space using clustering integrated with Random Forest further boosted the AUC to 99.6% in one cluster and direct detection of Android malware in another cluster, thus reducing the curse of dimensionality. Additionally, we found that feature reduction in detection models does improve the model efficiency (training and testing time) many folds without much penalty on the effectiveness of the detection model.


abstract: In this paper, we prove the Khavinson conjecture for hyperbolic harmonic functions on the unit ball. This conjecture was partially solved in \cite{JKM2020}.


title: Haidinger-Michelson rings in white light

abstract: We use coherence theory to explain why it is necessary to modify the conventional setup of a Michelson interferometer to obtain Haidinger rings with an extended source of white light. The modification consists of introducing a glass slide into one of the two arms of the interferometer. This insertion circumvents the drastic restriction imposed by the low temporal coherence of white light, which prevents the observation of interference rings with the conventional setup. In order to understand this restriction, we developed and implemented a criterion for observing interference fringes. The modified setup also makes it possible to perform a spectral interferometry experiment to analyze the output of the interferometer and determine the refractive index of the glass slide over the whole visible spectrum. The fit of measured data using Sellmeier's law gives the extrapolated value of the refractive index to the IR $n_\mathrm{IR}=1.5307\pm0.0003$ and the value of the characteristic wavelength $\lambda_0=(164.5\pm0.4)\;\mathrm{nm}$ of the oscillator responsible for the dispersion.


title: Update on the multimodal pathophysiological dataset of gradual cerebral ischemia in a cohort of juvenile pigs: auditory, sensory and high-frequency sensory evoked potentials

abstract: We expand from a spontaneous to an evoked potentials (EP) data set of brain electrical activities as electrocorticogram (ECoG) and electrothalamogram (EThG) in juvenile pig under various sedation, ischemia and recovery states. This EP data set includes three stimulation paradigms: auditory (AEP, 40 and 2000 Hz), sensory (SEP, left and right maxillary nerve) and high-frequency oscillations (HFO) SEP. This permits derivation of electroencephalogram (EEG) biomarkers of corticothalamic communication under these conditions. The data set is presented in full band sampled at 2000 Hz. We provide technical validation of the evoked responses for the states of sedation, ischemia and recovery. This extended data set now permits mutual inferences between spontaneous and evoked activities across the recorded modalities. Future studies on the dataset may contribute to the development of new brain monitoring technologies, which will facilitate the prevention of neurological injuries.


title: Order of torsion for reduction of linearly independent points for a family of Drinfeld modules

abstract: Let $q$ be a power of the prime number $p$, let $K={\mathbb F}_q(t)$, and let $r\ge 2$ be an integer. For points ${\mathbf a}, {\mathbf b}\in K$ which are $\mathbb{F}_q$-linearly independent, we show that there exist positive constants $N_0$ and $c_0$ such that for each integer $\ell\ge N_0$ and for each generator $\tau$ of ${\mathbb F}_{q^\ell}/{\mathbb F}_q$, we have that for all except $N_0$ values $\lambda\in{\overline{\mathbb{F}_q}}$, the corresponding specializations ${\mathbf a}, {\mathbf b}(\tau)$ and ${\mathbf b}(\tau)$ cannot have orders of degrees less than $c_0\log\log\ell$ as torsion points for the Drinfeld module $\Phi^{(\tau,\lambda)}:\mathbb{F}_q[T] {\longrightarrow} {\mathrm{End}}_{\overline{\mathbb{F}_q}}({\mathbb G}_a)$ (where ${\mathbb G}_a$ is the additive group scheme), given by $\Phi^{(\tau,\lambda)}_T(x)=\tau x+\lambda x^q + x^{q^r}$.


title: Four-body decays $B_{(s)} \rightarrow (K\pi)_{S/P} (K\pi)_{S/P}$ in the perturbative QCD approach

abstract: We investigate the four-body decays $B_{(s)} \rightarrow (K\pi)_{S/P} (K\pi)_{S/P}$ in the $K\pi$-pair invariant mass region around the $K^*(892)$ resonance in the perturbative QCD (PQCD) approach, where $(K\pi)_{S/P}$ denotes a $S$- or $P$-wave $K\pi$ configuration. The contributions from the $P$-wave resonance $K^*(892)$ and from possible $S$-wave components are included into the parametrization of the $K\pi$ two-meson distribution amplitudes, which are the nonperturbative inputs in this formalism. We calculate the branching ratio for each resonance and observe sizable $S$-wave contributions to the $B_s$ modes, well consistent with recently available LHCb data. We also deduce the polarization fractions for the $K^*\bar{K}^*$ final states, and compare our predictions to those in previous PQCD analyses of the corresponding two-body $B\to K^*\bar{K}^*$ decays. The polarization puzzle associated with the two $U$-spin related channels $B^0\rightarrow K^{*0}\bar{K}^{*0}$ and $B^0_s\rightarrow K^{*0}\bar{K}^{*0}$ still exists. In addition to direct $CP$ asymmetries, triple-product asymmetries and $S$-wave induced direct $CP$ asymmetries originating from the interference among various helicity amplitudes, are presented. It is shown that true triple-product asymmetries are rather small, while direct $CP$ asymmetries and $S$-wave-induced direct $CP$ asymmetries are significant in some decays. Our results will be subject to stringent tests with precise data from $B$ factories in the near future.


title: Identification of Significant Permissions for Efficient Android Malware Detection

abstract: Since Google unveiled Android OS for smartphones, malware are thriving with 3Vs, i.e. volume, velocity, and variety. A recent report indicates that one out of every five business/industry mobile application leaks sensitive personal data. Traditional signature/heuristic-based malware detection systems are unable to cope up with current malware challenges and thus threaten the Android ecosystem. Therefore recently researchers have started exploring machine learning and deep learning based malware detection systems. In this paper, we performed a comprehensive feature analysis to identify the significant Android permissions and propose an efficient Android malware detection system using machine learning and deep neural network. We constructed a set of $16$ permissions ($8\%$ of the total set) derived from variance threshold, auto-encoders, and principal component analysis to build a malware detection engine that consumes less train and test time without significant compromise on the model accuracy. Our experimental results show that the Android malware detection model based on the random forest classifier is most balanced and achieves the highest area under curve score of $97.7\%$, which is better than the current state-of-art systems. We also observed that deep neural networks attain comparable accuracy to the baseline results but with a massive computational penalty.


title: DMPC: A Data-and Model-Driven Approach to Predictive Control

abstract: This work presents DMPC (Data-and Model-Driven Predictive Control) to solve control problems in which some of the constraints or parts of the objective function are known, while others are entirely unknown to the controller. It is assumed that there is an exogenous ``black box'' system, e.g. a machine learning technique, that predicts the value of the unknown functions for a given trajectory. DMPC (1) provides an approach to merge both the model-based and black-box systems; (2) can cope with very little data and is sample efficient, building its solutions based on recently generated trajectories; and (3) improves its cost in each iteration until converging to an optimal trajectory, typically needing only a few trials even for nonlinear dynamics and objectives. Theoretical analysis of the algorithm is presented, proving that the quality of the trajectory does not worsen with each new iteration, as well as providing bounds on the complexity. We apply the DMPC algorithm to the motion planning of an autonomous vehicle with nonlinear dynamics.


title: Erdos-Renyi laws for exponentially and polynomially mixing dynamical systems

abstract: Erdos-Renyi limit laws give the length scale of a time-window over which time-averages in Birkhoff sums have a non-trivial almost-sure limit. We establish Erdos-Renyi type limit laws for Holder observables on dynamical systems modeled by Young Towers with exponential and polynomial tails. This extends earlier results on Erdos-Renyi limit laws to a broad class of dynamical systems with some degree of hyperbolicity.


title: Partitioned difference families: the storm has not yet passed

abstract: Two years ago, we alarmed the scientific community about the large number of bad papers in the literature on {\it zero difference balanced functions}, where direct proofs of seemingly new results are presented in an unnecessarily lengthy and convoluted way. Indeed, these results had been proved long before and very easily in terms of difference families. In spite of our report, papers of the same kind continue to proliferate. Regrettably, a further attempt to put the topic in order seems unavoidable. While some authors now follow our recommendation of using the terminology of {\it partitioned difference families}, their methods are still the same and their results are often trivial or even wrong. In this note, we show how a very recent paper of this type can be easily dealt with.


title: Spectra of variants of distance matrices of graphs and digraphs: a survey

abstract: Distance matrices of graphs were introduced by Graham and Pollack in 1971 to study a problem in communications. Since then, there has been extensive research on the distance matrices of graphs -- a 2014 survey by Aouchiche and Hansen on spectra of distance matrices of graphs lists more than 150 references. In the last ten years, variants such as the distance Laplacian, the distance signless Laplacian, and the normalized distance Laplacian matrix of a graph have been studied. After a brief description of the early history of the distance matrix and its motivating problem, this survey focuses on comparing and contrasting techniques and results for the four types of distance matrices. Digraphs are treated separately after the discussion of graphs, including discussion of similarities and differences between graphs and digraphs. New results are presented that complement existing results, including results for some the matrices on unimodality of characteristic polynomials for graphs, preservation of parameters by cospectrality for graphs, and bounds on spectral radii for digraphs.


title: Maximum Approximate Bernstein Likelihood Estimation of Densities in a Two-sample Semiparametric Model

abstract: Maximum likelihood estimators are proposed for the parameters and the densities in a semiparametric density ratio model in which the nonparametric baseline density is approximated by the Bernstein polynomial model. The EM algorithm is used to obtain the maximum approximate Bernstein likelihood estimates. Simulation study shows that the performance of the proposed method is much better than the existing ones. The proposed method is illustrated by real data examples. Some asymptotic results are also presented and proved.


title: The index of $\mathbb{T}^{\text{an}}$ in $\mathbb{T}$

abstract: For a prime level $N$, we compute the index of the anemic Hecke algebra of weight 2 level $\Gamma_0(N)$ forms inside the full Hecke algebra. We prove that the index is fully described by the presence of weight 1 Katz forms of the same level.


title: Rejection of Smooth GPS Time Synchronization Attacks via Sparse Techniques

abstract: This paper presents a novel time synchronization attack (TSA) model for the Global Positioning System (GPS) based on clock data behavior changes in a higher-order derivative domain. Further, the time synchronization attack rejection and mitigation based on sparse domain (TSARM-S) is presented. TSAs affect stationary GPS receivers in applications where precise timing is required, such as cellular communications, financial transactions, and monitoring of the electric power grid. In the present work, the clock bias and clock drift are monitored at higher-order clock data derivatives where the TSA is seen as a sparse spike-like event. The smoothness of the attack relates to the derivative order where the sparsity is observed. The proposed method jointly estimates a dynamic solution for GPS timing and rejects behavior changes based on such sparse events. An evaluation procedure is presented for two testbeds, namely a commercial receiver and a software-defined radio. Further, the proposed method is evaluated against distinct real-dataset Texas Spoofing Test Battery (TEXBAT) scenarios. Combined synthetic and real-data results show an average RMS clock bias error of 12.08 m for the SDR platform, and 45.74 m for the commercial device. Further, the technique is evaluated against state-of-the-art mitigation techniques and in a spoofing-plus-multipath scenario for robustness. Finally, TSARM-S can be potentially optimized and implemented in commercial devices via a firmware upgrade.


title: On Concentration Inequalities for Vector-Valued Lipschitz Functions

abstract: We derive two upper bounds for the probability of deviation of a vector-valued Lipschitz function of a collection of random variables from its expected value. The resulting upper bounds can be tighter than bounds obtained by a direct application of a classical theorem due to Bobkov and G\"{o}tze.


title: OpenICS: Open Image Compressive Sensing Toolbox and Benchmark

abstract: We present OpenICS, an image compressive sensing toolbox that includes multiple image compressive sensing and reconstruction algorithms proposed in the past decade. Due to the lack of standardization in the implementation and evaluation of the proposed algorithms, the application of image compressive sensing in the real-world is limited. We believe this toolbox is the first framework that provides a unified and standardized implementation of multiple image compressive sensing algorithms. In addition, we also conduct a benchmarking study on the methods included in this framework from two aspects: reconstruction accuracy and reconstruction efficiency. We wish this toolbox and benchmark can serve the growing research community of compressive sensing and the industry applying image compressive sensing to new problems as well as developing new methods more efficiently. Code and models are available at https://github.com/PSCLab-ASU/OpenICS. The project is still under maintenance, and we will keep this document updated.


title: Spin and charge interconversion in Dirac semimetal thin films

abstract: We report spin-to-charge and charge-to-spin conversion at room temperature in heterostructure devices that interface an archetypal Dirac semimetal, Cd3As2, with a metallic ferromagnet, Ni0.80Fe0.20 (permalloy). The spin-charge interconversion is detected by both spin torque ferromagnetic resonance and ferromagnetic resonance driven spin pumping. Analysis of the symmetric and anti-symmetric components of the mixing voltage in spin torque ferromagnetic resonance and the frequency and power dependence of the spin pumping signal show that the behavior of these processes is consistent with previously reported spin-charge interconversion mechanisms in heavy metals, topological insulators, and Weyl semimetals. We find that the efficiency of spin-charge interconversion in Cd3As2/permalloy bilayers can be comparable to that in heavy metals. We discuss the underlying mechanisms by comparing our results with first principles calculations.


title: Feedback Coding for Active Learning

abstract: The iterative selection of examples for labeling in active machine learning is conceptually similar to feedback channel coding in information theory: in both tasks, the objective is to seek a minimal sequence of actions to encode information in the presence of noise. While this high-level overlap has been previously noted, there remain open questions on how to best formulate active learning as a communications system to leverage existing analysis and algorithms in feedback coding. In this work, we formally identify and leverage the structural commonalities between the two problems, including the characterization of encoder and noisy channel components, to design a new algorithm. Specifically, we develop an optimal transport-based feedback coding scheme called Approximate Posterior Matching (APM) for the task of active example selection and explore its application to Bayesian logistic regression, a popular model in active learning. We evaluate APM on a variety of datasets and demonstrate learning performance comparable to existing active learning methods, at a reduced computational cost. These results demonstrate the potential of directly deploying concepts from feedback channel coding to design efficient active learning strategies.


title: Simultaneous Tactile Exploration and Grasp Refinement for Unknown Objects

abstract: This paper addresses the problem of simultaneously exploring an unknown object to model its shape, using tactile sensors on robotic fingers, while also improving finger placement to optimise grasp stability. In many situations, a robot will have only a partial camera view of the near side of an observed object, for which the far side remains occluded. We show how an initial grasp attempt, based on an initial guess of the overall object shape, yields tactile glances of the far side of the object which enable the shape estimate and consequently the successive grasps to be improved. We propose a grasp exploration approach using a probabilistic representation of shape, based on Gaussian Process Implicit Surfaces. This representation enables initial partial vision data to be augmented with additional data from successive tactile glances. This is combined with a probabilistic estimate of grasp quality to refine grasp configurations. When choosing the next set of finger placements, a bi-objective optimisation method is used to mutually maximise grasp quality and improve shape representation during successive grasp attempts. Experimental results show that the proposed approach yields stable grasp configurations more efficiently than a baseline method, while also yielding improved shape estimate of the grasped object.


title: Nonlinearity and wideband parametric amplification in an NbTiN microstrip transmission line

abstract: The nonlinear response associated with the current dependence of the superconducting kinetic inductance was studied in capacitively shunted NbTiN microstrip transmission lines. It was found that the inductance per unit length of one microstrip line could be changed by up to 20% by applying a DC current, corresponding to a single pass time delay of 0.7 ns. To investigate nonlinear dissipation, Bragg reflectors were placed on either end of a section of this type of transmission line, creating resonances over a range of frequencies. From the change in the resonance linewidth and amplitude with DC current, the ratio of the reactive to the dissipative response of the line was found to be 788. The low dissipation makes these transmission lines suitable for a number of applications that are microwave and millimeter-wave band analogues of nonlinear optical processes. As an example, by applying a millimeter-wave pump tone, very wide band parametric amplification was observed between about 3 and 34 GHz. Use as a current variable delay line for an on-chip millimeter-wave Fourier transform spectrometer is also considered.


title: Achieving Competitive Play Through Bottom-Up Approach in Semantic Segmentation

abstract: With the renaissance of neural networks, object detection has slowly shifted from a bottom-up recognition problem to a top-down approach. Best in class algorithms enumerate a near-complete list of objects and classify each into object/not object. In this paper, we show that strong performance can still be achieved using a bottom-up approach for vision-based object recognition tasks and achieve competitive video game play. We propose PuckNet, which is used to detect four extreme points (top, left, bottom, and right-most points) and one center point of objects using a fully convolutional neural network. Object detection is then a purely keypoint-based appearance estimation problem, without implicit feature learning or region classification. The method proposed herein performs on-par with the best in class region-based detection methods, with a bounding box AP of 36.4% on COCO test-dev. In addition, the extreme points estimated directly resolve into a rectangular object mask, with a COCO Mask AP of 17.6%, outperforming the Mask AP of vanilla bounding boxes. Guided segmentation of extreme points further improves this to 32.1% Mask AP. We applied the PuckNet vision system to the SuperTuxKart video game to test it's capacity to achieve competitive play in dynamic and co-operative multiplayer environments.


title: Emotion recognition techniques with rule based and machine learning approaches

abstract: Emotion recognition using digital image processing is a multifarious task because facial emotions depend on warped facial features as well as on gender, age, and culture. Furthermore, there are several factors such as varied illumination and intricate settings that increase complexity in facial emotion recognition. In this paper, we used four salient facial features, Eyebrows, Mouth opening, Mouth corners, and Forehead wrinkles to identifying emotions from normal, occluded and partially-occluded images. We have employed rule-based approach and developed new methods to extract aforementioned facial features similar to local bit patterns using novel techniques. We propose new methods to detect eye location, eyebrow contraction, and mouth corners. For eye detection, the proposed methods are Enhancement of Cr Red (ECrR) and Suppression of Cr Blue (SCrB) which results in 98% accuracy. Additionally, for eyebrow contraction detection, we propose two techniques (1) Morphological Gradient Image Intensity (MGII) and (2) Degree of Curvature Line (DCL). Additionally, we present a new method for mouth corners detection. For classification purpose, we use an individual classifier, majority voting (MV) and weighted majority voting (WMV) methods which mimic Human Emotions Sensitivity (HES). These methods are straightforward to implement, improve the accuracy of results, and work best for emotion recognition using partially occluded images. It is ascertained from the results that our method outperforms previous approaches. Overall accuracy rates are around 94%. The processing time on one image using processor core i5 is ~0.12 sec.


title: Software Training in HEP

abstract: Long term sustainability of the high energy physics (HEP) research software ecosystem is essential for the field. With upgrades and new facilities coming online throughout the 2020s this will only become increasingly relevant throughout this decade. Meeting this sustainability challenge requires a workforce with a combination of HEP domain knowledge and advanced software skills. The required software skills fall into three broad groups. The first is fundamental and generic software engineering (e.g. Unix, version control,C++, continuous integration). The second is knowledge of domain specific HEP packages and practices (e.g., the ROOT data format and analysis framework). The third is more advanced knowledge involving more specialized techniques. These include parallel programming, machine learning and data science tools, and techniques to preserve software projects at all scales. This paper dis-cusses the collective software training program in HEP and its activities led by the HEP Software Foundation (HSF) and the Institute for Research and Innovation in Software in HEP (IRIS-HEP). The program equips participants with an array of software skills that serve as ingredients from which solutions to the computing challenges of HEP can be formed. Beyond serving the community by ensuring that members are able to pursue research goals, this program serves individuals by providing intellectual capital and transferable skills that are becoming increasingly important to careers in the realm of software and computing, whether inside or outside HEP


title: Distribution Grid Modeling Using Smart Meter Data

abstract: The knowledge of distribution grid models, including topologies and line impedances, is essential to grid monitoring, control and protection. However, this information is often unavailable, incomplete or outdated. The increasing deployment of smart meters (SMs) provides a unique opportunity to address this issue. This paper proposes a two-stage data-driven framework for distribution grid modeling using SM data. In the first stage, we propose to identify the topology via reconstructing a weighted Laplacian matrix of distribution networks, which is mathematically proven to be robust against moderately heterogeneous R/X profiles. In the second stage, we develop nonlinear least absolute deviations (LAD) and least squares (LS) regression models to estimate line impedances of single branches based on a nonlinear inverse power flow, which is then embedded within a bottom-up sweep algorithm to achieve the identification across the network in a branch-wise manner. Because the estimation models are inherently non-convex programs and NP-hard, we specially address their tractable convex relaxations and verify the exactness. In addition, we design a conductor library to significantly narrow down the solution space. Numerical results on the modified IEEE 13-bus, 37-bus and 69-bus test feeders validate the effectiveness of the proposed methods.


title: Antipredator behavior in the rock-paper-scissors model

abstract: When faced with an imminent risk of predation, many animals react to escape consumption. Antipredator strategies are performed by individuals acting as a group to intimidate predators and minimize the damage when attacked. We study the antipredator prey response in spatial tritrophic systems with cyclic species dominance using the rock-paper-scissors game. The impact of the antipredator behavior is local, with the predation probability reducing exponentially with the number of preys in the predator's neighborhood. In contrast to the standard Lotka-Volterra implementation of the rock-paper-scissors model, where no spiral waves appear, our outcomes show that the antipredator behavior leads to spiral patterns from random initial conditions. The results show that the predation risk decreases exponentially with the level of antipredator strength. Finally, we investigate the coexistence probability and verify that antipredator behavior may jeopardize biodiversity for high mobility. Our findings may help biologists to understand ecosystems formed by species whose individuals behave strategically to resist predation.


title: A simple method for improving the accuracy of Chung-Lu random graph generation

abstract: Random graph models play a central role in network analysis. The Chung-Lu model, which connects nodes based on their expected degrees is of particular interest. It is widely used to generate null-graph models with expected degree sequences as well as implicitly define network measures such as modularity. Despite its popularity, practical methods for generating instances of Chung-Lu model-based graphs do relatively poor jobs in terms of accurately realizing many degree sequences. We introduce a simple method for improving the accuracy of Chung-Lu graph generation. Our method uses a Poisson approximation to define a linear system describing the expected degree sequence to be output from the model using standard generation techniques. We then use the inverse of this system to determine an appropriate input corresponding to the desired output. We give a closed form expression for this inverse and show that it may be used to drastically reduce error for many degree distributions.


title: Model-Agnostic Defense for Lane Detection against Adversarial Attack

abstract: Susceptibility of neural networks to adversarial attack prompts serious safety concerns for lane detection efforts, a domain where such models have been widely applied. Recent work on adversarial road patches have successfully induced perception of lane lines with arbitrary form, presenting an avenue for rogue control of vehicle behavior. In this paper, we propose a modular lane verification system that can catch such threats before the autonomous driving system is misled while remaining agnostic to the particular lane detection model. Our experiments show that implementing the system with a simple convolutional neural network (CNN) can defend against a wide gamut of attacks on lane detection models. With a 10% impact to inference time, we can detect 96% of bounded non-adaptive attacks, 90% of bounded adaptive attacks, and 98% of patch attacks while preserving accurate identification at least 95% of true lanes, indicating that our proposed verification system is effective at mitigating lane detection security risks with minimal overhead.


title: Generalized cyclotomic mappings: Switching between polynomial, cyclotomic, and wreath product form

abstract: This paper is concerned with so-called index $d$ generalized cyclotomic mappings of a finite field $\mathbb{F}_q$, which are functions $\mathbb{F}_q\rightarrow\mathbb{F}_q$ that agree with a suitable monomial function $x\mapsto ax^r$ on each coset of the index $d$ subgroup of $\mathbb{F}_q^{\ast}$. We discuss two important rewriting procedures in the context of generalized cyclotomic mappings and present applications thereof that concern index $d$ generalized cyclotomic permutations of $\mathbb{F}_q$ and pertain to cycle structures, the classification of $(q-1)$-cycles and involutions, as well as inversion.


title: Graded covering of a supermanifold I. The case of a Lie supergroup

abstract: We generalize the Donagi and Witten construction of a first obstruction class for splitting of a supermanifold via differential operators using the theory of $n$-fold vector bundles and graded manifolds. Applying the generalized Donagi--Witten construction we obtain a family of embeddings of the category of supermanifolds into the category of $n$-fold vector bundles and into the category of graded manifolds. This leads to a realization of any non-split supermanifold in terms of a collection of vector bundles and some morphism between them. Further we study the images of these embeddings into the category of graded manifolds in the case of a Lie supergroup and a Lie superalgebra. We show that these images satisfy universal property of a graded covering or a graded semicovering.


title: Estimating dust attenuation from galactic spectra. II. Stellar and gas attenuation in star-forming and diffuse ionized gas regions in MaNGA

abstract: We investigate the dust attenuation in both stellar populations and ionized gas in kpc-scale regions in nearby galaxies, using integral field spectroscopy data from MaNGA MPL-9. We identify star-forming (HII) and diffuse ionized gas (DIG) regions from MaNGA datacubes. From the stacked spectrum of each region, we measure the stellar attenuation, $E(B-V)_{\rm star}$, using the technique developed by Li et al.(2020), as well as the gas attenuation, $E(B-V)_{\rm gas}$, from the Balmer decrement. We then examine the correlation of $E(B-V)_{\rm star}$, $E(B-V)_{\rm gas}$, $E(B-V)_{\rm gas}-E(B-V)_{\rm star}$ and $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$ with 16 regional/global properties, and for regions with different $\rm H{\alpha}$ surface brightnesses ($\Sigma_{\rm H\alpha}$). We find a stronger correlation between $E(B-V)_{\rm star}$ and $E(B-V)_{\rm gas}$ in regions of higher $\Sigma_{\rm H\alpha}$. Luminosity-weighted age ($t_L$) is found to be the property that is the most strongly correlated with $E(B-V)_{\rm star}$, and consequently with $E(B-V)_{\rm gas}-E(B-V)_{\rm star}$ and $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$. At fixed $\Sigma_{\rm H\alpha}$, $\log_{10}t_L$ is linearly and negatively correlated with $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$ at all ages. Gas-phase metallicity and ionization level are important for the attenuation in the gas. Our results indicate that the ionizing source for DIG regions is likely distributed in the outer-skirt of galaxies, while for HII regions our results can be well explained by the two-component dust model of Charlot & Fall (2000).


title: Smooth Convex Optimization using Sub-Zeroth-Order Oracles

abstract: We consider the problem of minimizing a smooth, Lipschitz, convex function over a compact, convex set using sub-zeroth-order oracles: an oracle that outputs the sign of the directional derivative for a given point and a given direction, an oracle that compares the function values for a given pair of points, and an oracle that outputs a noisy function value for a given point. We show that the sample complexity of optimization using these oracles is polynomial in the relevant parameters. The optimization algorithm that we provide for the comparator oracle is the first algorithm with a known rate of convergence that is polynomial in the number of dimensions. We also give an algorithm for the noisy-value oracle that incurs a regret of $\tilde{\mathcal{O}}(n^{3.75} T^{0.75})$ (ignoring the other factors and logarithmic dependencies) where $n$ is the number of dimensions and $T$ is the number of queries.


title: Learning Proposals for Probabilistic Programs with Inference Combinators

abstract: We develop operators for construction of proposals in probabilistic programs, which we refer to as inference combinators. Inference combinators define a grammar over importance samplers that compose primitive operations such as application of a transition kernel and importance resampling. Proposals in these samplers can be parameterized using neural networks, which in turn can be trained by optimizing variational objectives. The result is a framework for user-programmable variational methods that are correct by construction and can be tailored to specific models. We demonstrate the flexibility of this framework by implementing advanced variational methods based on amortized Gibbs sampling and annealing.


title: Route Lengths in Invariant Spatial Tree Networks

abstract: Is there a constant $r_0$ such that, in any invariant tree network linking rate-$1$ Poisson points in the plane, the mean within-network distance between points at Euclidean distance $r$ is infinite for $r > r_0$? We prove a slightly weaker result. This is a continuum analog of a result of Benjamini et al (2001) on invariant spanning trees of the integer lattice.


title: Extreme space weather events caused by super active regions during solar cycles 21-24

abstract: Extreme space weather events including $\ge$X5.0 flares, ground level enhancement (GLE) events and super geomagnetic storms (Dst $\le$ -250 nT) caused by super active regions (SARs) during solar cycles 21-24 were studied. The total number of $\ge$X5.0 solar flares was 62, 41 of them were X5.0-X9.9 flares and 21 of them were $\ge$X10.0 flares. We found that 83.9\% of the $\ge$X5.0 flares were produced by SARs. 78.05\% of the X5.0-X9.9 and 95.24\% of the $\ge$X10.0 solar flares were produced by SARs. 46 GLEs registered during solar cycles 21-24, and 25 GLEs were caused by SARs, indicating that 54.3\% of the GLEs were caused by SARs. 24 super geomagnetic storms were recorded during solar cycles 21-24, and 12 of them were caused by SARs, namely 50\% of the super geomagnetic storms are caused by SARs. It is found that only 29 SARs can produce $\ge$X5.0 flares, 15 SARs can produce GLEs and 10 SARs can produce super geomagnetic storms. Of the 51 SARs, only 33 SARs can produce at least one extreme space weather event, while none of the rest 18 SARs can produce an extreme space weather event. There were only 4 SARs, each of them can produce not only a $\ge$X5.0 flare, but also a GLE event and a super geomagnetic storm. Most of the extreme space weather events caused by the SARs appeared during solar cycles 22 and 23, especially for GLE events and super geomagnetic storms. The longitudinal distributions of source locations for the extreme space weather events caused by SARs were also studied.


title: Robust learning under clean-label attack

abstract: We study the problem of robust learning under clean-label data-poisoning attacks, where the attacker injects (an arbitrary set of) correctly-labeled examples to the training set to fool the algorithm into making mistakes on specific test instances at test time. The learning goal is to minimize the attackable rate (the probability mass of attackable test instances), which is more difficult than optimal PAC learning. As we show, any robust algorithm with diminishing attackable rate can achieve the optimal dependence on $\epsilon$ in its PAC sample complexity, i.e., $O(1/\epsilon)$. On the other hand, the attackable rate might be large even for some optimal PAC learners, e.g., SVM for linear classifiers. Furthermore, we show that the class of linear hypotheses is not robustly learnable when the data distribution has zero margin and is robustly learnable in the case of positive margin but requires sample complexity exponential in the dimension. For a general hypothesis class with bounded VC dimension, if the attacker is limited to add at most $t>0$ poison examples, the optimal robust learning sample complexity grows almost linearly with $t$.


title: Secondary Homological Stability for Unordered Configuration Spaces

abstract: Secondary homological stability is a recently discovered stability pattern for the homology of a sequence of spaces exhibiting homological stability in a range where homological stability does not hold. We prove secondary homological stability for the homology of the unordered configuration spaces of a connected manifold. The main difficulty is the case that the manifold is closed because there are no obvious maps inducing stability and the homology eventually is periodic instead of stable. We resolve this issue by constructing a new chain-level stabilization map for configuration spaces.


title: Convolutional Normalization: Improving Deep Convolutional Network Robustness and Training

abstract: Normalization techniques have become a basic component in modern convolutional neural networks (ConvNets). In particular, many recent works demonstrate that promoting the orthogonality of the weights helps train deep models and improve robustness. For ConvNets, most existing methods are based on penalizing or normalizing weight matrices derived from concatenating or flattening the convolutional kernels. These methods often destroy or ignore the benign convolutional structure of the kernels; therefore, they are often expensive or impractical for deep ConvNets. In contrast, we introduce a simple and efficient "Convolutional Normalization" (ConvNorm) method that can fully exploit the convolutional structure in the Fourier domain and serve as a simple plug-and-play module to be conveniently incorporated into any ConvNets. Our method is inspired by recent work on preconditioning methods for convolutional sparse coding and can effectively promote each layer's channel-wise isometry. Furthermore, we show that our ConvNorm can reduce the layerwise spectral norm of the weight matrices and hence improve the Lipschitzness of the network, leading to easier training and improved robustness for deep ConvNets. Applied to classification under noise corruptions and generative adversarial network (GAN), we show that the ConvNorm improves the robustness of common ConvNets such as ResNet and the performance of GAN. We verify our findings via numerical experiments on CIFAR and ImageNet.


title: BEAUTY Powered BEAST

abstract: We study distribution-free goodness-of-fit tests with the proposed Binary Expansion Approximation of UniformiTY (BEAUTY) approach. This method generalizes the renowned Euler's formula, and approximates the characteristic function of any copula through a linear combination of expectations of binary interactions from marginal binary expansions. This novel theory enables a unification of many important tests of independence via approximations from specific quadratic forms of symmetry statistics, where the deterministic weight matrix characterizes the power properties of each test. To achieve a robust power, we examine test statistics with data-adaptive weights, referred to as the Binary Expansion Adaptive Symmetry Test (BEAST). Using properties of the binary expansion filtration, we demonstrate that the Neyman-Pearson test of uniformity can be approximated by an oracle weighted sum of symmetry statistics. The BEAST with this oracle provides a useful benchmark of feasible power. To approach this oracle power, we devise the BEAST through a regularized resampling approximation of the oracle test. The BEAST improves the empirical power of many existing tests against a wide spectrum of common alternatives and delivers a clear interpretation of dependency forms when significant.

categories: stat.ME cs.LG math.ST stat.AP stat.ML stat.TH


title: Bayesian filtering for nonlinear stochastic systems using holonomic gradient method with integral transform

abstract: This paper proposes a symbolic-numeric Bayesian filtering method for a class of discrete-time nonlinear stochastic systems to achieve high accuracy with a relatively small online computational cost. The proposed method is based on the holonomic gradient method (HGM), which is a symbolic-numeric method to evaluate integrals efficiently depending on several parameters. By approximating the posterior probability density function (PDF) of the state as a Gaussian PDF, the update process of its mean and variance can be formulated as evaluations of several integrals that exactly take into account the nonlinearity of the system dynamics. An integral transform is used to evaluate these integrals more efficiently using the HGM than our previous method. Further, a numerical example is provided to demonstrate the efficiency of the proposed method compared to other existing methods.

categories: math.NA cs.NA cs.SC cs.SY eess.SP eess.SY


title: Token-Modification Adversarial Attacks for Natural Language Processing: A Survey

abstract: Many adversarial attacks target natural language processing systems, most of which succeed through modifying the individual tokens of a document. Despite the apparent uniqueness of each of these attacks, fundamentally they are simply a distinct configuration of four components: a goal function, allowable transformations, a search method, and constraints. In this survey, we systematically present the different components used throughout the literature, using an attack-independent framework which allows for easy comparison and categorisation of components. Our work aims to serve as a comprehensive guide for newcomers to the field and to spark targeted research into refining the individual attack components.


title: Modeling and prediction of COVID-19 in the United States considering population behavior and vaccination

abstract: COVID-19 has devastated the entire global community. Vaccines present an opportunity to mitigate the pandemic; however, the effect of vaccination coupled with the behavioral response of the population is not well understood. We propose a model that incorporates two important dynamically varying population behaviors: level of caution and sense of safety. Level of caution increases with the number of infectious cases, while an increasing sense of safety with increased vaccination lowers precautionary behaviors. To the best of our knowledge, this is the first model that can effectively reproduce the complete time history of COVID-19 infections for various regions of the United States and provides relatable measures of dynamic changes in the population behavior and disease transmission rates. We propose a parameter d_I as a direct measure of a population's caution against an infectious disease, that can be obtained from the ongoing new infectious cases. The model provides a method for quantitative measure of critical infectious disease attributes for a population including highest disease transmission rate, effective disease transmission rate, and disease related precautionary behavior. We predict future COVID-19 pandemic trends in the United States accounting for vaccine rollout and behavioral response. Although a high rate of vaccination is critical to quickly end the pandemic, we find that a return towards pre-pandemic social behavior due to increased sense of safety during vaccine deployment, can cause an alarming surge in infections. Our results indicate that at the current rate of vaccination, the new infection cases for COVID-19 in the United States will approach zero by the end of August 2021. The model can be used for predicting future epidemic and pandemic dynamics before and during vaccination.


title: Role of unitary correlation operator on high-momentum antisymmetrized molecular dynamics using bare NN interaction for 3H and 4He

abstract: We extend the high-momentum antisymmetrized molecular dynamics (HMAMD) by incorporating the short-range part of the unitary correlation operator method (UCOM) as the variational method of finite nuclei. In this HMAMD+UCOM calculation of light nuclei, the HMAMD is mainly in charge of the tensor correlation with up to the four-body correlation, while the short-range correlation is further improved by using the UCOM. The binding energies of the 3H and 4He nuclei are calculated with this HMAMD+UCOM using the AV8' bare nucleon-nucleon (NN) interaction. The different roles of the short-range and tensor correlations from the HMAMD and UCOM are analyzed in the numerical results. Compared with the previous calculations based on the different variational methods, this newly extended HMAMD+UCOM method can almost provide the consistent results with the ab initio results.


title: Some thoughts on pseudoprimes

abstract: We consider several problems about pseudoprimes. First, we look at the issue of their distribution in residue classes. There is a literature on this topic in the case that the residue class is coprime to the modulus. Here we provide some robust statistics in both these cases and the general case. In particular we tabulate all even pseudoprimes to $10^{16}$. Second, we prove a recent conjecture of Ordowski: the set of integers $n$ which are a pseudoprime to some base which is a proper divisor of $n$ has an asymptotic density.


title: Consequential LCA for territorial and multimodal transportation policies: method and application to the free-floating e-scooter disruption in Paris

abstract: The indirect environmental impacts of transport disruptions in urban mobility are frequently overlooked due to a lack of appropriate assessment methods. Consequential Life Cycle Assessment (CLCA) is a method to capture the environmental consequences of the entire cause and effect chain of these disruptions but has never been adapted to transportat disruption at the city scale. This paper proposes a mathematical formalization of CLCA applied to a territorial mobility change. The method is applied to quantify the impact on climate change of the breakthrough of free-floating e-scooters (FFES) in Paris. A FFES user survey is conducted to estimate the modal shifts due to FFES. Trip substitutions from all the Parisian modes concerned are considered - personal or shared bicycles and motor scooters, private car, taxi and ride-hailing, bus, streetcar, metro and RER (the Paris metropolitan area mass rapid transit system). All these Parisian modes are assessed for the first time using LCA. Final results estimate that over one year, the FFES generated an extra thirteen thousand tons of CO2eq under an assumption of one million users, mainly due to major shifts coming from lower-emitting modes (60% from the metro and the RER, 22% from active modes). Recommendations are given to enhance their carbon footprint. A scenario analysis shows that increasing the lifetime mileage is insufficient to get a positive balance: reducing drastically servicing emissions is also required. A sensitivity analysis switching the French electricity mix for eleven other country mixes suggests a better climate change effect of the FFES in similar metropolitan areas with higher electricity carbon intensity, such as in Germany and China. Finally, the novelty and the limits of the method are discussed, as well as the results and the role of e-scooters, micromobility, and shared vehicles towards a sustainable mobility.


title: On the global small solution of 2-D Prandtl system with initial data in the optimal Gevrey class

abstract: Motivated by \cite{DG19}, we prove the global existence and large time behavior of small solutions to 2-D Prandtl system for data with Gevrey 2 regularity in the $x$ variable and Sobolev regularity in the $y$ variable. In particular, we extend the global well-posedness result in \cite{PZ5} for 2-D Prandtl system with analytic data to data with optimal Gevery regularity in the sense of \cite{Ger1}.


title: Enhancing hierarchical surrogate-assisted evolutionary algorithm for high-dimensional expensive optimization via random projection

abstract: By remarkably reducing real fitness evaluations, surrogate-assisted evolutionary algorithms (SAEAs), especially hierarchical SAEAs, have been shown to be effective in solving computationally expensive optimization problems. The success of hierarchical SAEAs mainly profits from the potential benefit of their global surrogate models known as "blessing of uncertainty" and the high accuracy of local models. However, their performance leaves room for improvement on highdimensional problems since now it is still challenging to build accurate enough local models due to the huge solution space. Directing against this issue, this study proposes a new hierarchical SAEA by training local surrogate models with the help of the random projection technique. Instead of executing training in the original high-dimensional solution space, the new algorithm first randomly projects training samples onto a set of low-dimensional subspaces, then trains a surrogate model in each subspace, and finally achieves evaluations of candidate solutions by averaging the resulting models. Experimental results on six benchmark functions of 100 and 200 dimensions demonstrate that random projection can significantly improve the accuracy of local surrogate models and the new proposed hierarchical SAEA possesses an obvious edge over state-of-the-art SAEAs


title: Decision Making in Monopoly using a Hybrid Deep Reinforcement Learning Approach

abstract: Learning to adapt and make real-time informed decisions in a dynamic and complex environment is a challenging problem. Monopoly is a popular strategic board game that requires players to make multiple decisions during the game. Decision-making in Monopoly involves many real-world elements such as strategizing, luck, and modeling of opponent's policies. In this paper, we present novel representations for the state and action space for the full version of Monopoly and define an improved reward function. Using these, we show that our deep reinforcement learning agent can learn winning strategies for Monopoly against different fixed-policy agents. In Monopoly, players can take multiple actions even if it is not their turn to roll the dice. Some of these actions occur more frequently than others, resulting in a skewed distribution that adversely affects the performance of the learning agent. To tackle the non-uniform distribution of actions, we propose a hybrid approach that combines deep reinforcement learning (for frequent but complex decisions) with a fixed policy approach (for infrequent but straightforward decisions). Experimental results show that our hybrid agent outperforms a standard deep reinforcement learning agent by 30% in the number of games won against fixed-policy agents.


title: Meta-learning One-class Classifiers with Eigenvalue Solvers for Supervised Anomaly Detection

abstract: Neural network-based anomaly detection methods have shown to achieve high performance. However, they require a large amount of training data for each task. We propose a neural network-based meta-learning method for supervised anomaly detection. The proposed method improves the anomaly detection performance on unseen tasks, which contains a few labeled normal and anomalous instances, by meta-training with various datasets. With a meta-learning framework, quick adaptation to each task and its effective backpropagation are important since the model is trained by the adaptation for each epoch. Our model enables them by formulating adaptation as a generalized eigenvalue problem with one-class classification; its global optimum solution is obtained, and the solver is differentiable. We experimentally demonstrate that the proposed method achieves better performance than existing anomaly detection and few-shot learning methods on various datasets.


title: No dependence of radio properties of brightest group galaxies on the luminosity gap

abstract: We study the radio and optical properties of the brightest group galaxies (BGGs) in a sample of galaxy groups from the SDSS DR7. The luminosity difference between the BGG and the second ranked galaxy in the group (known as the luminosity, or magnitude, gap) has been used as a probe for the level of galaxy interaction for the BGG within the group. We study the properties of BGGs with magnitude gaps in the range 0-2.7 magnitudes, in order to investigate any relation between luminosity gap and the radio properties of the BGG. In order to eliminate selection biases, we ensure that all variations in stellar mass are accounted for. We then confirm that, at fixed stellar mass, there are no significant variations in the optical properties of the BGGs over the full range of luminosity gaps studied. We compare these optical results with the EAGLE hydrodynamical simulations and find broad consistency with the observational data. Using EAGLE we also confirm that no trends begin to arise in the simulated data at luminosity gaps beyond our observational limits. Finally, we find that, at fixed stellar mass, the fraction of BGGs that are radio-loud also shows no trends as a function of luminosity gap. We examine how the BGG offset from the center of group may affect the radio results and find no significant trend for the fraction of radio-loud BGGs with magnitude gap in either the BGG samples with greater or less than 100kpc offset from the center of group.


title: Accelerating Recommendation System Training by Leveraging Popular Choices

abstract: Recommender models are commonly used to suggest relevant items to a user for e-commerce and online advertisement-based applications. These models use massive embedding tables to store numerical representation of items' and users' categorical variables (memory intensive) and employ neural networks (compute intensive) to generate final recommendations. Training these large-scale recommendation models is evolving to require increasing data and compute resources. The highly parallel neural networks portion of these models can benefit from GPU acceleration however, large embedding tables often cannot fit in the limited-capacity GPU device memory. Hence, this paper deep dives into the semantics of training data and obtains insights about the feature access, transfer, and usage patterns of these models. We observe that, due to the popularity of certain inputs, the accesses to the embeddings are highly skewed with a few embedding entries being accessed up to 10000x more. This paper leverages this asymmetrical access pattern to offer a framework, called FAE, and proposes a hot-embedding aware data layout for training recommender models. This layout utilizes the scarce GPU memory for storing the highly accessed embeddings, thus reduces the data transfers from CPU to GPU. At the same time, FAE engages the GPU to accelerate the executions of these hot embedding entries. Experiments on production-scale recommendation models with real datasets show that FAE reduces the overall training time by 2.3x and 1.52x in comparison to XDL CPU-only and XDL CPU-GPU execution while maintaining baseline accuracy


title: Enhanced Superconductivity in the Se-substituted 1T-PdTe$_2$

abstract: Two-dimensional transition metal dichalcogenide PdTe$_2$ recently attracts much attention due to its phase coexistence of type-II Dirac semimetal and type-I superconductivity. Here we report a 67 % enhancement of superconducting transition temperature in the 1T-PdSeTe in comparison to that of PdTe2 through partial substitution of Te atoms by Se. The superconductivity has been unambiguously confirmed by the magnetization, resistivity and specific heat measurements. 1T-PdSeTe shows type-II superconductivity with large anisotropy and non-bulk superconductivity nature with volume fraction ~ 20 % estimated from magnetic and heat capacity measurements. 1T-PdSeTe expands the family of superconducting transition metal dichalcogenides and thus provides additional insights for understanding superconductivity and topological physics in the 1T-PdTe$_2$ system


title: Nambu dynamics and its noncanonical Hamiltonian representation in many degrees of freedom systems

abstract: Nambu dynamics is a generalized Hamiltonian dynamics of more than two variables, whose time evolutions are given by the Nambu bracket, a generalization of the canonical Poisson bracket. Nambu dynamics can always be represented in the form of noncanonical Hamiltonian dynamics by defining the noncanonical Poisson bracket by means of the Nambu bracket. For the time evolution to be consistent, the Nambu bracket must satisfy the fundamental identity, while the noncanonical Poisson bracket must satisfy the Jacobi identity. However, in many degrees of freedom systems, it is well known that the fundamental identity does not hold. In the present paper we show that, even if the fundamental identity is violated, the Jacobi identity for the corresponding noncanonical Hamiltonian dynamics could hold. As an example, we evaluate these identities for a semiclassical system of two coupled oscillators.


title: Formation of Organic Color Centers in Air-Suspended Carbon Nanotubes Using Vapor-Phase Reaction

abstract: Organic color centers in single-walled carbon nanotubes have demonstrated exceptional ability to generate single photons at room temperature in the telecom range. Combining the color centers with pristine air-suspended tubes would be desirable for improved performance, but all current synthetic methods occur in solution which makes them incompatible. Here we demonstrate formation of color centers in air-suspended nanotubes using vapor-phase reaction. Functionalization is directly verified on the same nanotubes by photoluminescence spectroscopy, with unambiguous statistics from more than a few thousand individual nanotubes. The color centers show a strong diameter-dependent emission intensity, which can be explained with a theoretical model for chemical reactivity taking into account strain along the tube curvature. We are also able to estimate the defect density by comparing the experiments with simulations based on a one-dimensional diffusion equation, whereas the analysis of diameter dependent peak energies gives insight to the nature of the dopant states. Time-resolved measurements show a longer lifetime for color center emission compared to E$_{11}$ exciton states. Our results highlight the influence of the tube structure on vapor-phase reactivity and emission properties, providing guidelines for development of high-performance near-infrared quantum light sources.


title: An informed thought experiment exploring the potential for a paradigm shift in aquatic food production

abstract: The Neolithic Revolution began c. 10000 years ago and is characterised by the ultimate, near complete transition from hunting and gathering to agricultural food production on land. The Neolithic Revolution is thought to have been catalysed by a combination of local population pressure, cultural diffusion, property rights and climate change. We undertake a thought experiment that examines trends in these key hypothesised catalysts and patters of today to explore whether society could be on a path towards another paradigm shift in food production: away from hunting of wild fish towards a transition to mostly fish farming. We find similar environmental and cultural pressures have driven the rapid rise of aquaculture, during a period that has now been coined the Blue Revolution, providing impetus for such a transition in coming decades to centuries. We also highlight the interacting and often mutually reinforcing impacts of 1)technological and scientific advancement, 2)environmental awareness and collective action and 3)globalisation and trade influencing the trajectory and momentum of the Blue Revolution. We present two qualitative narratives that broadly fall within two future trajectories: 1)a ubiquitous aquaculture transition and 20commercial aquaculture and fisheries coexistence. This scenarios approach aims to encourage logical, forward thinking, and innovative solutions to complex systems dynamics. Scenario-based thought experiments are useful to explore large scale questions, increase the accessibility to a wider readership and ideally catalyse discussion around proactive governance mechanisms. We argue the future is not fixed and society now has greater foresight and capacity to choose the workable balance between fisheries sand aquaculture that supports economic, environmental, cultural and social objectives through combined planning, policies and management.