title: Inference of the chromospheric magnetic field configuration of solar plage using the Ca II 8542 {\AA} line

abstract: It has so far proven impossible to reproduce all aspects of the solar plage chromosphere in quasi-realistic numerical models. The magnetic field configuration in the lower atmosphere is one of the few free parameters in such simulations. The literature only offers proxy-based estimates of the field strength, as it is difficult to obtain observational constraints in this region. Sufficiently sensitive spectro-polarimetric measurements require a high signal-to-noise ratio, spectral resolution, and cadence, which are at the limit of current capabilities. We use critically sampled spectro-polarimetric observations of the \cair line obtained with the CRISP instrument of the Swedish 1-m Solar Telescope to study the strength and inclination of the chromospheric magnetic field of a plage region. This will provide direct physics-based estimates of these values, which could aid modelers to put constraints on plage models. We increased the signal-to-noise ratio of the data by applying several methods including deep learning and PCA. We estimated the noise level to be $1\cdot10^{-3} I_c$. We then used STiC, a non-local thermodynamic equilibrium (NLTE) inversion code to infer the atmospheric structure and magnetic field pixel by pixel. We are able to infer the magnetic field strength and inclination for a plage region and for fibrils in the surrounding canopy. In the plage we report an absolute field strength of $|B| =440 \pm 90$ G, with an inclination of $10^\circ \pm 16^\circ$ with respect to the local vertical. This value for $|B|$ is roughly double of what was reported previously, while the inclination matches previous studies done in the photosphere. In the fibrillar region we found $|B| = 300 \pm 50$ G, with an inclination of $50^\circ \pm 13^\circ$.


title: Stratification of canopy magnetic fields in a plage region. Constraints from a spatially-regularized weak-field approximation method

abstract: The role of magnetic fields in the chromospheric heating problem remains greatly unconstrained. Most theoretical predictions from numerical models rely on a magnetic configuration, field strength and connectivity whose details have not been well established with observational studies. High-resolution studies of chromospheric magnetic fields in plage are very scarce or non-existent in general. Our aim is to study the stratification of the magnetic field vector in plage regions. We use high-spatial resolution full-Stokes observations acquired with CRISP instrument at the Swedish 1-m Solar Telescope in the Mg I $\lambda$5173, Na I $\lambda$5896 and Ca II $\lambda$8542 lines. We have developed a spatially-regularized weak-field approximation (WFA) method based on the idea of spatial regularization. This method allows for a fast computation of magnetic field maps for an extended field of view. The fidelity of this new technique has been assessed using a snapshot from a realistic 3D magnetohydrodynamics simulation. We have derived the depth-stratification of the line-of-sight component of the magnetic field from the photosphere to the chromosphere in a plage region. The magnetic fields are concentrated in the intergranular lanes in the photosphere and expand horizontally toward the chromosphere, filling all the space and forming a canopy. Our results suggest that the lower boundary of this canopy must be located around 400-600 km from the photosphere. The mean canopy total magnetic field strength in the lower chromosphere ($z\approx760$ km) is 658 G. At $z=1160$ km we estimate $<B_\parallel>\approx 417$ G. We propose a modification to the WFA that improves its applicability to data with worse signal-to-noise ratio. These methods provide a quick and reliable way of studying multi-layer magnetic field observations without the many difficulties inherent to other inversion methods.


title: First-Order Model-Checking in Random Graphs and Complex Networks

abstract: Complex networks are everywhere. They appear for example in the form of biological networks, social networks, or computer networks and have been studied extensively. Efficient algorithms to solve problems on complex networks play a central role in today's society. Algorithmic meta-theorems show that many problems can be solved efficiently. Since logic is a powerful tool to model problems, it has been used to obtain very general meta-theorems. In this work, we consider all problems definable in first-order logic and analyze which properties of complex networks allow them to be solved efficiently. The mathematical tool to describe complex networks are random graph models. We define a property of random graph models called $\alpha$-power-law-boundedness. Roughly speaking, a random graph is $\alpha$-power-law-bounded if it does not admit strong clustering and its degree sequence is bounded by a power-law distribution with exponent at least $\alpha$ (i.e. the fraction of vertices with degree $k$ is roughly $O(k^{-\alpha})$). We solve the first-order model-checking problem (parameterized by the length of the formula) in almost linear FPT time on random graph models satisfying this property with $\alpha \ge 3$. This means in particular that one can solve every problem expressible in first-order logic in almost linear expected time on these random graph models. This includes for example preferential attachment graphs, Chung-Lu graphs, configuration graphs, and sparse Erd\H{o}s-R\'{e}nyi graphs. Our results match known hardness results and generalize previous tractability results on this topic.


title: Rank-metric codes over arbitrary Galois extensions and rank analogues of Reed-Muller codes

abstract: This paper extends the study of rank-metric codes in extension fields $\mathbb{L}$ equipped with an arbitrary Galois group $G = \mathrm{Gal}(\mathbb{L}/\mathbb{K})$. We propose a framework for studying these codes as subspaces of the group algebra $\mathbb{L}[G]$, and we relate this point of view with usual notions of rank-metric codes in $\mathbb{L}^N$ or in $\mathbb{K}^{N\times N}$, where $N = [\mathbb{L} : \mathbb{K}]$. We then adapt the notion of error-correcting pairs to this context, in order to provide a non-trivial decoding algorithm for these codes. We then focus on the case where $G$ is abelian, which leads us to see codewords as elements of a multivariate skew polynomial ring. We prove that we can bound the dimension of the vector space of zeroes of these polynomials, depending of their degree. This result can be seen as an analogue of Alon-F\"uredi theorem -- and by means, of Schwartz-Zippel lemma -- in the rank metric. Finally, we construct the counterparts of Reed-Muller codes in the rank metric, and we give their parameters. We also show the connection between these codes and classical Reed-Muller codes in the case where $\mathbb{L}$ is a Kummer extension.


title: Case study: Mapping potential informal settlements areas in Tegucigalpa with machine learning to plan ground survey

abstract: Data collection through censuses is conducted every 10 years on average in Latin America, making it difficult to monitor the growth and support needed by communities living in these settlements. Conducting a field survey requires logistical resources to be able to do it exhaustively. The increasing availability of open data, high-resolution satellite images, and free software to process them allow us to be able to do so in a scalable way based on the analysis of these sources of information. This case study shows the collaboration between Dymaxion Labs and the NGO Techo to employ machine learning techniques to create the first informal settlements census of Tegucigalpa, Honduras.


title: Iwasawa theory for quadratic Hilbert modular forms

abstract: We study the Iwasawa main conjecture for quadratic Hilbert modular forms over the p-cyclotomic tower. Using an Euler system in the cohomology of Siegel modular varieties, we prove the "Kato divisibility" of the Iwasawa main conjecture under certain technical hypotheses. By comparing this result with the opposite divisibility due to Wan, we obtain the full Main Conjecture over the cyclotomic Zp-extension. As a consequence, we prove new cases of the Bloch--Kato conjecture for quadratic Hilbert modular forms, and of the equivariant Birch--Swinnerton-Dyer conjecture in analytic rank 0 for elliptic curves over real quadratic fields twisted by Dirichlet characters.


title: TweetsCOV19 -- A Knowledge Base of Semantically Annotated Tweets about the COVID-19 Pandemic

abstract: Publicly available social media archives facilitate research in the social sciences and provide corpora for training and testing a wide range of machine learning and natural language processing methods. With respect to the recent outbreak of the Coronavirus disease 2019 (COVID-19), online discourse on Twitter reflects public opinion and perception related to the pandemic itself as well as mitigating measures and their societal impact. Understanding such discourse, its evolution, and interdependencies with real-world events or (mis)information can foster valuable insights. On the other hand, such corpora are crucial facilitators for computational methods addressing tasks such as sentiment analysis, event detection, or entity recognition. However, obtaining, archiving, and semantically annotating large amounts of tweets is costly. In this paper, we describe TweetsCOV19, a publicly available knowledge base of currently more than 8 million tweets, spanning October 2019 - April 2020. Metadata about the tweets as well as extracted entities, hashtags, user mentions, sentiments, and URLs are exposed using established RDF/S vocabularies, providing an unprecedented knowledge base for a range of knowledge discovery tasks. Next to a description of the dataset and its extraction and annotation process, we present an initial analysis and use cases of the corpus.


title: Statistical inference of assortative community structures

abstract: We develop a principled methodology to infer assortative communities in networks based on a nonparametric Bayesian formulation of the planted partition model. We show that this approach succeeds in finding statistically significant assortative modules in networks, unlike alternatives such as modularity maximization, which systematically overfits both in artificial as well as in empirical examples. In addition, we show that our method is not subject to a resolution limit, and can uncover an arbitrarily large number of communities, as long as there is statistical evidence for them. Our formulation is amenable to model selection procedures, which allow us to compare it to more general approaches based on the stochastic block model, and in this way reveal whether assortativity is in fact the dominating large-scale mixing pattern. We perform this comparison with several empirical networks, and identify numerous cases where the network's assortativity is exaggerated by traditional community detection methods, and we show how a more faithful degree of assortativity can be identified.

categories: physics.soc-ph cond-mat.dis-nn cs.SI stat.ML


title: Experimental validation of coronagraphic focal-plane wavefront sensing for future segmented space telescopes

abstract: Direct imaging of Earth-like planets from space requires dedicated observatories, combining large segmented apertures with instruments and techniques such as coronagraphs, wavefront sensors, and wavefront control in order to reach the high contrast of 10^10 that is required. The complexity of these systems would be increased by the segmentation of the primary mirror, which allows for the larger diameters necessary to image Earth-like planets but also introduces specific patterns in the image due to the pupil shape and segmentation and making high-contrast imaging more challenging. Among these defects, the phasing errors of the primary mirror are a strong limitation to the performance. In this paper, we focus on the wavefront sensing of segment phasing errors for a high-contrast system, using the COronagraphic Focal plane wave-Front Estimation for Exoplanet detection (COFFEE) technique. We implemented and tested COFFEE on the High-contrast imaging for Complex Aperture Telescopes (HiCAT) testbed, in a configuration without any coronagraph and with a classical Lyot coronagraph, to reconstruct errors applied on a 37 segment mirror. We analysed the quality and limitations of the reconstructions. We demonstrate that COFFEE is able to estimate correctly the phasing errors of a segmented telescope for piston, tip, and tilt aberrations of typically 100nm RMS. We also identified the limitations of COFFEE for the reconstruction of low-order wavefront modes, which are highly filtered by the coronagraph. This is illustrated using two focal plane mask sizes on HiCAT. We discuss possible solutions, both in the hardware system and in the COFFEE optimizer, to mitigate these issues.


title: Gray tensor products and lax functors of $(\infty,2)$-categories

abstract: We give a definition of the Gray tensor product in the setting of scaled simplicial sets which is associative and forms a left Quillen bifunctor with respect to the bicategorical model category of Lurie. We then introduce a notion of oplax functor in this setting, and use it in order to characterize the Gray tensor product by means of a universal property. A similar characterization was used by Gaitsgory and Rozenblyum in their definition of the Gray product, thus giving a promising lead for comparing the two settings.


title: Percolation Thresholds for Robust Network Connectivity

abstract: Communication networks, power grids, and transportation networks are all examples of networks whose performance depends on reliable connectivity of their underlying network components even in the presence of usual network dynamics due to mobility, node or edge failures, and varying traffic loads. Percolation theory quantifies the threshold value of a local control parameter such as a node occupation (resp., deletion) probability or an edge activation (resp., removal) probability above (resp., below) which there exists a giant connected component (GCC), a connected component comprising of a number of occupied nodes and active edges whose size is proportional to the size of the network itself. Any pair of occupied nodes in the GCC is connected via at least one path comprised of active edges and occupied nodes. The mere existence of the GCC itself does not guarantee that the long-range connectivity would be robust, e.g., to random link or node failures due to network dynamics. In this paper, we explore new percolation thresholds that guarantee not only spanning network connectivity, but also robustness. We define and analyze four measures of robust network connectivity, explore their interrelationships, and numerically evaluate the respective robust percolation thresholds for the 2D square lattice.


title: Quantumized Microwave Detection Based on $\Lambda$-Type Three-level Superconducting System: HMM Modeling and Performance Prediction

abstract: We adopt artificial $\Lambda$-type three-level system with superconducting devices for microwave signal detection, where the signal intensity reaches the level of discrete photons instead of continuous waveform. Based on the state transition principles of the three-level system, we propose a statistical model for microwave signal detection. Moreover, we investigate the achievable transmission rate and signal detection based on the statistical model. It is predicted that the proposed detection can achieve significantly higher sensitivity compared with the currently deployed 4G/5G communication system. We further characterize the received signal considering the saturation phonomenon, which reveals negligible performance degradation caused by saturation under weak received power regime.


title: A Data-driven Market Simulator for Small Data Environments

abstract: Neural network based data-driven market simulation unveils a new and flexible way of modelling financial time series without imposing assumptions on the underlying stochastic dynamics. Though in this sense generative market simulation is model-free, the concrete modelling choices are nevertheless decisive for the features of the simulated paths. We give a brief overview of currently used generative modelling approaches and performance evaluation metrics for financial time series, and address some of the challenges to achieve good results in the latter. We also contrast some classical approaches of market simulation with simulation based on generative modelling and highlight some advantages and pitfalls of the new approach. While most generative models tend to rely on large amounts of training data, we present here a generative model that works reliably in environments where the amount of available training data is notoriously small. Furthermore, we show how a rough paths perspective combined with a parsimonious Variational Autoencoder framework provides a powerful way for encoding and evaluating financial time series in such environments where available training data is scarce. Finally, we also propose a suitable performance evaluation metric for financial time series and discuss some connections of our Market Generator to deep hedging.

categories: q-fin.ST cs.LG q-fin.CP q-fin.MF stat.ML


title: Examining the Effect of COVID-19 on Foreign Exchange Rate and Stock Market -- An Applied Insight into the Variable Effects of Lockdown on Indian Economy

abstract: Since March 25, 2020, India had been under a nation-wide lockdown announced as a response to the spread of SARS-CoV-2 and COVID-19 and has resorted to a process of 'unlocking' the lockdown over the past couple of months. This work attempts to examine the effect of novel coronavirus 2019 (COVID-19) and its resulting disease, the COVID-19, on the foreign exchange rates and stock market performances of India using secondary data over a span of 112 days spanning between March 11 and June 30, 2020. The study explores whether the causal relationships and directions among the growth rate of confirmed cases (GROWTHC), exchange rate (GEX) and SENSEX value (GSENSEX) are remaining the same across different pre and post-lockdown phases, attempting to capture any potential changes over time via the vector autoregressive (VAR) models. A positive correlation is found between the growth rate of confirmed cases and the growth rate of exchange rate, and a negative correlation between the growth rate of confirmed cases and the growth rate of SENSEX value. However, on applying a vector autoregressive (VAR) model, it is observed that an increase in the confirmed COVID-19 cases causes no significant change in the values of the exchange rate and SENSEX index. The result varies if the analysis is split across different time periods - before lockdown, the four phases of lockdown, and the first phase of unlock. Nuanced and sensible interpretations of the numeric results indicate significant variability across time in terms of the relation between the variables of interest. The detailed knowledge about the varying patterns of dependence could potentially help the policy makers and investors of India in order to develop their strategies to cope up with the situation.


title: Active Janus particles in a complex plasma

abstract: Active Janus particles suspended in a plasma were studied experimentally. The Janus particles were micron-size plastic microspheres, one half of which was coated with a thin layer of platinum. They were suspended in the plasma sheath of a radio-frequency discharge in argon at low pressure. The Janus particles moved in characteristic looped trajectories suggesting a combination of spinning and circling motion; their interactions led to the emergence of rich dynamics characterized by non-Maxwellian velocity distribution. The particle propulsion mechanism is discussed, the main force driving the particle motion is identified as photophoretic force.


title: Neural Networks Based Beam Codebooks: Learning mmWave Massive MIMO Beams that Adapt to Deployment and Hardware

abstract: Millimeter wave (mmWave) and massive MIMO systems are intrinsic components of 5G and beyond. These systems rely on using beamforming codebooks for both initial access and data transmission. Current beam codebooks, however, generally consist of a large number of narrow beams that scan all possible directions, even if these directions are never used. This leads to very large training overhead. Further, these codebooks do not normally account for the hardware impairments or the possible non-uniform array geometries, and their calibration is an expensive process. To overcome these limitations, this paper develops an efficient online machine learning framework that learns how to adapt the codebook beam patterns to the specific deployment, surrounding environment, user distribution, and hardware characteristics. This is done by designing a novel complex-valued neural network architecture in which the neuron weights directly model the beamforming weights of the analog phase shifters, accounting for the key hardware constraints such as the constant-modulus and quantized-angles. This model learns the codebook beams through online and self-supervised training avoiding the need for explicit channel state information. This respects the practical situations where the channel is either unavailable, imperfect, or hard to obtain, especially in the presence of hardware impairments. Simulation results highlight the capability of the proposed solution in learning environment and hardware aware beam codebooks, which can significantly reduce the training overhead, enhance the achievable data rates, and improve the robustness against possible hardware impairments.


title: Liouville theorems for a stationary and non-stationary coupled system of liquid crystal flows in local Morrey spaces

abstract: We consider here the simplified Ericksen-Leslie system on the whole three-dimensional space. This system deals with the incompressible Navier-Stokes equations strongly coupled with a harmonic map flow which models the dynamical behavior for nematic liquid crystals. For both the stationary (time independing) case and the non-stationary (time depending) case, using the fairly general framework of a kind of local Morrey spaces, we obtain some a priori conditions on the unknowns of this coupled system to prove that they vanish identically. This results are known as Liouville-type theorems. As a bi-product, our theorems also improve some well-known results on Liouville-type theorems for the particular case of classical Navier-Stokes equations


title: The radio/X-ray burst from SGR 1935+2154: radiation mechanisms and the possible QPOs

abstract: Recently, a fast radio burst (FRB)-like event is found to be associated with a Galactic magnetar, SGR 1935+2154, accompanied by an X-ray burst. We find this radio burst challenges the typical emission mechanisms involving magnetars, which includes coherent curvature radiation from charged bunches, fast magnetosonic (FMS) wave, synchrotron maser from shocks, and the pulsar-like mechanism for low-twist magnetars. More specifically, we find that (1) the X-rays are most-likely to be produced inside the magnetosphere. (2) For the coherent curvature radiation from the decay of Alfv\'en wave, it will generally predict a duration ($\lesssim0.1$ ms) smaller than observations, because of the strong twists of magnetic field lines and the internal damping of Alfv\'en waves. (3) The FMS wave model predicts a very low emission frequency $\nu_{\rm p}\sim0.03$ MHz $\ll$ GHz, unless it is produced inside the magnetosphere. But the absorption effect of the magnetospheric FMS wave model remains to be studied. (4) The synchrotron maser model is challenged, because observations show that the peaks in both X-ray and radio light curves are with the same temporal separation $\Delta t_{\rm FRB}=\Delta t_\gamma\approx0.03$ s, while it would predict $\Delta t_{\rm FRB}\ll\Delta t_\gamma$. (5) It seems to be difficult to directly apply the low-twist pulsar-like mechanism to flaring magnetars, as magnetar activity can significantly deform the magnetosphere. (6) We suggested four possibilities to study the general properties of FRBs for future observations, especially the possibility of identifying quasi-periodic oscillations with period $\sim1-10$ ms in double/multiple-peaked FRBs.


title: Growth of finitely generated simple Lie algebras

abstract: We realize any submultiplicative increasing function which is equivalent to a polynomial proportion of itself as the growth function of a finitely generated simple Lie algebra. As an application, we resolve two open problems posed by Petrogradsky on Lie algebras with intermediate growth.


title: Source Code Comments: Overlooked in the Realm of Code Clone Detection

abstract: Reusing code can produce duplicate or near-duplicate code clones in code repositories. Current code clone detection techniques, like Program Dependence Graphs, rely on code structure and their dependencies to detect clones. These techniques are expensive, using large amounts of processing power, time, and memory. In practice, programmers often utilize code comments to comprehend and reuse code, as comments carry important domain knowledge. But current code detection techniques ignore code comments, mainly due to the ambiguity of the English language. Recent advances in information retrieval techniques may have the potential to utilize code comments for clone detection. We investigated this by empirically comparing the accuracy of detecting clones with solely comments versus solely source code (without comments) on the JHotDraw package, which contains 315 classes and 27K lines of code. To detect clones at the file level, we used a topic modeling technique, Latent Dirichlet Allocation, to analyze code comments and GRAPLE -- utilizing Program Dependency Graph -- to analyze code. Our results show 94.86 recall and 84.21 precision with Latent Dirichlet Allocation and 28.7 recall and 55.39 precision using GRAPLE. We found Latent Dirichlet Allocation generated false positives in cases where programs lacked quality comments. But this limitation can be addressed by using a hybrid approach: utilizing code comments at the file level to reduce the clone set and then using Program Dependency Graph-based techniques at the method level to detect precise clones. Our further analysis across Java and Python packages, Java Swing and PyGUI, found a recall of 74.86\% and a precision of 84.21\%. Our findings call for reexamining the assumptions regarding the use of code comments in current clone detection techniques.


title: Near Mean Motion Resonance of Terrestrial Planet Pair Induced by Giant Planet: Application to Kepler-68 System

abstract: In this work, we investigate configuration formation of two inner terrestrial planets near mean motion resonance (MMRs) induced by the perturbation of a distant gas-giant for the Kepler-68 system, by conducting thousands of numerical simulations. The results show that the formation of terrestrial planets is relevant to the speed of Type I migration, the mass of planets, and the existence of giant planet. The mass and eccentricity of the giant planet may play a crucial role in shaping the final configuration of the system. The inner planet pair can be trapped in 5:3 or 7:4 MMRs if the giant planet revolves the central star with an eccentric orbit, which is similar to the observed configuration of Kepler-68. Moreover, we find that the eccentricity of the middle planet can be excited to roughly 0.2 if the giant planet is more massive than 5 $M_J$,otherwise the terrestrial planets are inclined to remain near-circular orbits. Our study may provide a likely formation scenario for the planetary systems that harbor several terrestrial planets near MMRs inside and one gas-giant exterior to them.


title: Existence and structure of symmetric Beltrami flows on compact $3$-manifolds

abstract: We show that for almost every given symmetry transformation of a Riemannian manifold there exists an eigenvector field of the curl operator, corresponding to a non-zero eigenvalue, which obeys the symmetry. More precisely, given a smooth, compact, oriented Riemannian $3$-manifold $(\bar{M},g)$ with (possibly empty) boundary and a smooth flow of isometries $\phi_t:\bar{M}\rightarrow \bar{M}$ we show that, if $\bar{M}$ has non-empty boundary or if the infinitesimal generator is not purely harmonic, there is a smooth vector field $X$, tangent to the boundary, which is an eigenfield of curl and satisfies $(\phi_t)_{*}X=X$, i.e. is invariant under the pushforward of the symmetry transformation. We then proceed to show that if the quantities involved are real analytic and $(\bar{M},g)$ has non-empty boundary, then Arnold's structure theorem applies to all eigenfields of curl, which obey a symmetry and appropriate boundary conditions. More generally we show that the structure theorem applies to all real analytic vector fields of non-vanishing helicity which obey some nontrivial symmetry. A byproduct of our proof is a characterisation of the flows of real analytic Killing fields on compact, connected, orientable $3$-manifolds with and without boundary.


title: Interference Cancellation Based Channel Estimation for Massive MIMO Systems with Time Shifted Pilots

abstract: In massive multiple-input multiple-output (MIMO) systems with time shifted pilot (TSP) schemes, the inter-group interference caused by the pilot contamination can be eliminated when the number of base station (BS) antennas M approaches infinity. However, M is finite in practice and the effectiveness of the TSP is limited by channel estimation errors. In this paper, it is analytically shown that the mean square channel estimation error (MSCEE) of the TSP is dominated by the intergroup data interference. To reduce the MSCEE in the finite antenna massive MIMO systems, an interference cancellation based channel estimation for the TSP (IC-TSP) is proposed, where the dominant inter-group data interference is canceled based on BS cooperation. To show the advantage of the IC-TSP, the additional overhead of IC-TSP is evaluated by considering different M and the coherence time of BS-BS channels. Furthermore, the impact of sectorization and compressed sensing based BS-BS channel estimation are also discussed. We show that when 128 < M < 2048, with the inter-group data interference from the nearest two cell layers being canceled, the IC-TSP achieves a spectral efficiency gain of more than 1.2 bps/Hz over the TSP.


title: Using rational homology circles to construct rational homology balls

abstract: Motivated by Akbulut-Larson's construction of Brieskorn spheres bounding rational homology 4-balls, we explore plumbed 3-manifolds that bound rational homology circles and use them to construct infinite families of rational homology 3-spheres that bound rational homology 4-balls. Some of these rational homology 3-spheres are new examples of integer homology 3-spheres that bound rational homology 4-balls, but do not bound integer homology 4-balls. In particular, we find infinite families of torus bundles over the circle that bound rational homology circles, provide a simple method for constructing more general plumbed 3-manifolds that bound rational homology circles, and show that, for example, $-1$-surgery along any unknotting number one knot $K$ with a positive crossing that can be switched to unknot $K$ bounds a rational homology 4-ball.


title: Quantum Computing for Finance: State of the Art and Future Prospects

abstract: This article outlines our point of view regarding the applicability, state-of-the-art, and potential of quantum computing for problems in finance. We provide an introduction to quantum computing as well as a survey on problem classes in finance that are computationally challenging classically and for which quantum computing algorithms are promising. In the main part, we describe in detail quantum algorithms for specific applications arising in financial services, such as those involving simulation, optimization, and machine learning problems. In addition, we include demonstrations of quantum algorithms on IBM Quantum back-ends and discuss the potential benefits of quantum algorithms for problems in financial services. We conclude with a summary of technical challenges and future prospects.


title: Model-Independent Perspectives on Coupled Dark Energy and the Swampland

abstract: We present a general model-independent approach to study the coupled dark energy and the string Swampland criteria. We show how the dark sector interaction is degenerated with the equation of state of dark energy in the context of the expansion of the Universe. With priors for either of them, the dynamics of dark energy and the dark sector interactions can be reconstructed together with the bounds of the Swampland criteria. Combining cosmic chronometers, baryon acoustic oscillation (BAO), and Type Ia supernovae our results suggest a mild $1 \sigma$ significance of dark sector interactions at low redshift for the coupled quintessence. The Lyman-$\alpha$ BAO at $z=2.34$ leads a $2 \sigma$ signal of nonzero interactions at high redshift. The implications for coupled quintessence are discussed.


title: Uncovering the Connections Between Adversarial Transferability and Knowledge Transferability

abstract: Knowledge transferability, or transfer learning, has been widely adopted to allow a pre-trained model in the source domain to be effectively adapted to downstream tasks in the target domain. It is thus important to explore and understand the factors affecting knowledge transferability. In this paper, as the first work, we analyze and demonstrate the connections between knowledge transferability and another important phenomenon--adversarial transferability, \emph{i.e.}, adversarial examples generated against one model can be transferred to attack other models. Our theoretical studies show that adversarial transferability indicates knowledge transferability and vice versa. Moreover, based on the theoretical insights, we propose two practical adversarial transferability metrics to characterize this process, serving as bidirectional indicators between adversarial and knowledge transferability. We conduct extensive experiments for different scenarios on diverse datasets, showing a positive correlation between adversarial transferability and knowledge transferability. Our findings will shed light on future research about effective knowledge transfer learning and adversarial transferability analyses.


title: Blockchain-Aided Flow Insertion and Verification in Software Defined Networks

abstract: The Internet of Things (IoT) connected by Software Defined Networking (SDN) promises to bring great benefits to cyber-physical systems. However, the increased attack surface offered by the growing number of connected vulnerable devices and complex nature of SDN control plane applications could overturn the huge benefits of such a system. This paper addresses the vulnerability of some unspecified security flaw in the SDN control plane application (such as a zero-day software vulnerability) which can be exploited to insert malicious flow rules in the switch that do not match network policies. Specifically, we propose a blockchain-as-a-service (BaaS) based framework that supports switch flow verification and insertion; and additionally provides straightforward deployment of blockchain technology within an existing SDN infrastructure. While use of an external BaaS brings straightforward deployment, it obscures knowledge of the blockchain agents who are responsible for flow conformance testing through a smart blockchain contract, leading to potential exploitation. Thus, we design a strategy to prevent the blockchain agents from acting arbitrarily, as this would result in what is termed a "moral hazard". We achieve this by developing a novel mathematical model of the fair reward scheme based on game theory. To understand the performance of our system, we evaluate our model using a Matlab based simulation framework. The simulation results demonstrate that the proposed algorithm balances the needs of the blockchain agents to maximise the overall social welfare, i.e. the sum of profits across all parties.


title: Taming neural networks with TUSLA: Non-convex learning via adaptive stochastic gradient Langevin algorithms

abstract: Artificial neural networks (ANNs) are typically highly nonlinear systems which are finely tuned via the optimization of their associated, non-convex loss functions. In many cases, the gradient of any such loss function has superlinear growth, making the use of the widely-accepted (stochastic) gradient descent methods, which are based on Euler numerical schemes, problematic. We offer a new learning algorithm based on an appropriately constructed variant of the popular stochastic gradient Langevin dynamics (SGLD), which is called tamed unadjusted stochastic Langevin algorithm (TUSLA). We also provide a nonasymptotic analysis of the new algorithm's convergence properties in the context of non-convex learning problems with the use of ANNs. Thus, we provide finite-time guarantees for TUSLA to find approximate minimizers of both empirical and population risks. The roots of the TUSLA algorithm are based on the taming technology for diffusion processes with superlinear coefficients as developed in \citet{tamed-euler, SabanisAoAP} and for MCMC algorithms in \citet{tula}. Numerical experiments are presented which confirm the theoretical findings and illustrate the need for the use of the new algorithm in comparison to vanilla SGLD within the framework of ANNs.


title: Surrogate Models for Rainfall Nowcasting

abstract: Nowcasting (or short-term weather forecasting) is particularly important in the case of extreme events as it helps prevent human losses. Many of our activities, however, also depend on the weather. Therefore, nowcasting has shown to be useful in many different domains. Currently, immediate rainfall forecasts in France are calculated using the Arome-NWC model developed by M\'et\'eo-France, which is a complex physical model. Arome-NWC forecasts are stored with a 15 minute time interval. A higher time resolution is, however, desirable for other meteorological applications. Complex model calculations, such as Arome-NWC, can be very expensive and time consuming. A surrogate model aims at producing results which are very close to the ones obtained using a complex model, but with largely reduced calculation times. Building a surrogate model requires only a few calculations with the real model. Once the surrogate model is built, further calculations can be quickly realized. In this study, we propose to build surrogate models for immediate rainfall forecasts with two different approaches: combining Proper Orthogonal Decomposition (POD) and Kriging, or combining POD and Random Forest (RF). We show that results obtained with our surrogate models are not only close to the ones obtained by Arome-NWC, but they also have a higher time resolution (1 minute) with a reduced calculation time.


title: Plumbene on a magnetic substrate: a combined STM and DFT study

abstract: As heavy analog of graphene, plumbene is a two-dimensional material with strong spin-orbit coupling effects. Using scanning tunneling microscopy (STM), we observe that Pb forms a flat honeycomb lattice on an Fe monolayer on Ir(111). In contrast, without the Fe layer, a c(2x4) structure of Pb on Ir(111) is found. We use density functional theory (DFT) calculations to rationalize these findings and analyze the impact of the hybridization on the plumbene band structure. In the unoccupied states the splitting of the Dirac cone by spin-orbit interaction is clearly observed while in the occupied states of the freestanding plumbene we find a band inversion that leads to the formation of a topologically non-trivial gap. Exchange splitting as mediated by the strong hybridization with the Fe layer drives a quantum spin Hall to quantum anomalous Hall state transition.


title: Generalized Maslov indices for non-Hamiltonian systems

abstract: We extend the definition of the Maslov index to a broad class of non-Hamiltonian dynamical systems. To do this, we introduce a family of topological spaces--which we call Maslov-Arnold spaces--that share key topological features with the Lagrangian Grassmannian, and hence admit a similar index theory. This family contains the Lagrangian Grassmannian, and much more. We construct a family of examples, called hyperplane Maslov-Arnold spaces, that are dense in the Grassmannian, and hence are much larger than the Lagrangian Grassmannian (which is a submanifold of positive codimension). The resulting index is then used to study eigenvalue problems for non-symmetric reaction-diffusion systems. A highlight of our analysis is a topological interpretation of the Turing instability: the bifurcation that occurs as one increases the ratio of diffusion coefficients corresponds to a change in the generalized Maslov index.


title: Mobility operator service capacity sharing contract design to risk-pool against network disruptions

abstract: We propose a new mechanism to design risk-pooling contracts between operators to facilitate horizontal cooperation to mitigate those costs and improve service resilience during disruptions. We formulate a novel two-stage stochastic multicommodity flow model to determine the cost savings of a coalition under different disruption scenarios and solve it using L-shaped method along with sample average approximation. Computational tests of the L-shaped method against deterministic equivalent method with sample average approximation are conducted for network instances with up to 64 nodes, 10 OD pairs, and 1024 scenarios. The results demonstrate that the solution algorithm only becomes computationally effective for larger size instances (above 128 nodes) and that SAA maintains a close approximation. The proposed model is applied to a regional multi-operator network in the Randstad area of the Netherlands, for four operators, 40 origin-destination pairs, and over 1400 links where disruption data is available. Using the proposed method, we identify stable cost allocations among four operating agencies that could yield a 66% improvement in overall network performance over not having any risk-pooling contract in place. Furthermore, the model allows policymakers to evaluate the sensitivity of any one operator's bargaining power to different network structures and disruption scenario distributions, as we illustrate for the HTM operator in Randstad.


title: Smart and Secure Wireless Communications via Reflecting Intelligent Surfaces: A Short Survey

abstract: With the emergence of the internet of things (IoT) technology, wireless connectivity should be more ubiquitous than ever. In fact, the availability of wireless connection everywhere comes with security threats that, unfortunately, cannot be handled by conventional cryptographic solutions alone, especially in heterogeneous and decentralized future wireless networks. In general, physical layer security (PLS) helps in bridging this gap by taking advantage of the fading propagation channel. Moreover, the adoption of reconfigurable intelligent surfaces (RIS) in wireless networks makes the PLS techniques more efficient by involving the channel into the design loop. In this paper, we conduct a comprehensive literature review on the RIS-assisted PLS for future wireless communications. We start by introducing the basic concepts of RISs and their different applications in wireless communication networks and the most common PLS performance metrics. Then, we focus on the review and classification of RIS-assisted PLS applications, exhibiting multiple scenarios, system models, objectives, and methodologies. In fact, most of the works in this field formulate an optimization problem to maximize the secrecy rate (SR) or secrecy capacity (SC) at a legitimate user by jointly optimizing the beamformer at the transmitter and the RIS's coefficients, while the differences are in the adopted methodology to optimally/sub-optimally approach the solution. We finalize this survey by presenting some insightful recommendations and suggesting open problems for future research extensions.


title: On Hodge polynomials of Singular Character Varieties

abstract: Let $\mathcal{X}_{\Gamma}G:=\mathrm{Hom}(\Gamma,G)/\!/G$ be the $G$-character variety of $\Gamma$, where $G$ is a complex reductive group and $\Gamma$ a finitely presented group. We introduce new techniques for computing Hodge-Deligne and Serre polynomials of $\mathcal{X}_{\Gamma}G$, and present some applications, focusing on the cases when $\Gamma$ is a free or free abelian group. Detailed constructions and proofs of the main results will appear elsewhere.


title: Exploring New Physics with O(keV) Electron Recoils in Direct Detection Experiments

abstract: Motivated by the recent XENON1T results, we explore various new physics models that can be discovered through searches for electron recoils in O(keV)-threshold direct-detection experiments. First, we consider the absorption of light bosons, either as dark matter relics or being produced directly in the Sun. In the latter case, we find that keV mass bosons produced in the Sun provide an adequate fit to the data but are excluded by stellar cooling constraints. We address this tension by introducing a novel Chameleon-like axion model, which can explain the excess while evading the stellar bounds. We find that absorption of bosonic dark matter provides a viable explanation for the excess only if the dark matter is a dark photon or an axion. In the latter case, photophobic axion couplings are necessary to avoid X-ray constraints. Second, we analyze models of dark matter-electron scattering to determine which models might explain the excess. Standard scattering of dark matter with electrons is generically in conflict with data from lower-threshold experiments. Momentum-dependent interactions with a heavy mediator can fit the data with dark matter mass heavier than a GeV but are generically in tension with collider constraints. Next, we consider dark matter consisting of two (or more) states that have a small mass splitting. The exothermic (down)scattering of the heavier state to the lighter state can fit the data for keV mass splittings. Finally, we consider a subcomponent of dark matter that is accelerated by scattering off cosmic rays, finding that dark matter interacting though an O(100 keV)-mass mediator can fit the data. The cross sections required in this scenario are, however, typically challenged by complementary probes of the light mediator. Throughout our study, we implement an unbinned Monte Carlo analysis and use an improved energy reconstruction of the XENON1T events.


title: Product formulas and convolutions for two-dimensional Laplace-Beltrami operators: beyond the trivial case

abstract: We introduce the notion of a family of convolution operators associated with a given elliptic partial differential operator. Such a convolution structure is shown to exist for a general class of Laplace-Beltrami operators on two-dimensional manifolds endowed with cone-like metrics. This structure gives rise to a convolution semigroup representation for the Markovian semigroup generated by the Laplace-Beltrami operator. In the particular case of the operator $\mathcal{L} = \partial_x^2 + {1 \over 2x} \partial_x + {1 \over x} \partial_\theta^2$ on $\mathbb{R}^+ \times \mathbb{T}$, we deduce the existence of a convolution structure for a two-dimensional integral transform whose kernel and inversion formula can be written in closed form in terms of confluent hypergeometric functions. The results of this paper can be interpreted as a natural extension of the theory of one-dimensional generalized convolutions to the framework of multiparameter eigenvalue problems.


title: Reverse Quantum Speed Limit: How Slow Quantum Battery can Discharge?

abstract: We introduce the notion of reverse quantum speed limit for arbitrary quantum evolution which answers a fundamental question: ``how slow a quantum system can evolve in time?" Using the geometrical approach to quantum mechanics, the reverse speed limit follows from the fact that the gauge invariant length of the reference section is always greater than the Fubini-Study distance on the projective Hilbert space of the quantum system. We illustrate the reverse speed limit for two-level quantum systems with an external driving Hamiltonian and show that our results hold well. We find several examples where our bound is tight. We also find one practical application of the reverse speed limit in discharging process of quantum batteries which answers the question: ``how slow quantum batteries can discharge?" Our result provides a lower bound on the average discharging power of quantum batteries.


title: Particle Shape Influences Settling and Sorting Behavior in Microfluidic Domains

abstract: We present a new numerical model to simulate settling trajectories of discretized individual or a mixture of particles of different geometrical shapes in a quiescent fluid and their flow trajectories in a flowing fluid. Simulations unveiled diverse particle settling trajectories as a function of their geometrical shape and density. The effects of the surface concavity of a boomerang particle and aspect ratio of a rectangular particle on the periodicity and amplitude of oscillations in their settling trajectories were numerically captured. Use of surrogate circular particles for settling or flowing of a mixture of non-circular particles were shown to miscalculate particle velocities by a factor of 0.9-2.2 and inaccurately determine the particles' trajectories. In a microfluidic chamber with particles of different shapes and sizes, simulations showed that steady vortices do not necessarily always control particle entrapments, nor do larger particles get selectively and consistently entrapped in steady vortices. Strikingly, a change in the shape of large particles from circular to elliptical resulted in stronger entrapments of smaller circular particles, but enhanced outflows of larger particles, which could be an alternative microfluidics-based method for sorting and separation of particles of different sizes and shapes.


title: Conjugation Curvature in Solvable Baumslag-Solitar Groups

abstract: For an element in $BS(1,n) = \langle t,a | tat^{-1} = a^n \rangle$ written in the normal form $t^{-u}a^vt^w$ with $u,w \geq 0$ and $v \in \mathbb{Z}$, we exhibit a geodesic word representing the element and give a formula for its word length with respect to the generating set $\{t,a\}$. Using this word length formula, we prove that there are sets of elements of positive density of positive, negative and zero conjugation curvature, as defined by Bar Natan, Duchin and Kropholler.


title: A mechanism to promote social behaviour in household load balancing

abstract: Reducing the peak energy consumption of households is essential for the effective use of renewable energy sources, in order to ensure that as much household demand as possible can be met by renewable sources. This entails spreading out the use of high-powered appliances such as dishwashers and washing machines throughout the day. Traditional approaches to this problem have relied on differential pricing set by a centralised utility company. But this mechanism has not been effective in promoting widespread shifting of appliance usage. Here we consider an alternative decentralised mechanism, where agents receive an initial allocation of time-slots to use their appliances and can then exchange these with other agents. If agents are willing to be more flexible in the exchanges they accept, then overall satisfaction, in terms of the percentage of agents time-slot preferences that are satisfied, will increase. This requires a mechanism that can incentivise agents to be more flexible. Building on previous work, we show that a mechanism incorporating social capital - the tracking of favours given and received - can incentivise agents to act flexibly and give favours by accepting exchanges that do not immediately benefit them. We demonstrate that a mechanism that tracks favours increases the overall satisfaction of agents, and crucially allows social agents that give favours to outcompete selfish agents that do not under payoff-biased social learning. Thus, even completely self-interested agents are expected to learn to produce socially beneficial outcomes.


title: Critical 3-hypergraphs (detailed version)

abstract: Given a 3-hypergraph $H$, a subset $M$ of $V(H)$ is a module of $H$ if for each $e\in E(H)$ such that $e\cap M\neq\emptyset$ and $e\setminus M\neq\emptyset$, there exists $m\in M$ such that $e\cap M=\{m\}$ and for every $n\in M$, we have $(e\setminus\{m\})\cup\{n\}\in E(H)$. For example, $\emptyset$, $V(H)$ and $\{v\}$, where $v\in V(H)$, are modules of $H$, called trivial. A 3-hypergraph is prime if all its modules are trivial. Furthermore, a prime 3-hypergraph is critical if all its induced subhypergraphs, obtained by removing one vertex, are not prime. We characterize the critical 3-hypergraphs.


title: Blockchain-Based Applications in Higher Education: A Systematic Mapping Study

abstract: The utilisation of blockchain has moved beyond digital currency to other fields such as health, the Internet of Things, and education. In this paper, we present a systematic mapping study to collect and analyse relevant research on blockchain technology related to the higher education field. The paper concentrates on two main themes. First, it examines state of the art in blockchain-based applications that have been developed for educational purposes. Second, it summarises the challenges and research gaps that need to be addressed in future studies.


title: From oscillation dip to oscillation valley in atmospheric neutrino experiments

abstract: Atmospheric neutrino experiments can show the "oscillation dip" feature in data, due to their sensitivity over a large $L/E$ range. In experiments that can distinguish between neutrinos and antineutrinos, like INO, oscillation dips can be observed in both these channels separately. We present the dip-identification algorithm employing a data-driven approach -- one that uses the asymmetry in the upward-going and downward-going events, binned in the reconstructed $L/E$ of muons -- to demonstrate the dip, which would confirm the oscillation hypothesis. We further propose, for the first time, the identification of an "oscillation valley" in the reconstructed ($E_\mu$,$\,\cos\theta_\mu$) plane, feasible for detectors like ICAL having excellent muon energy and direction resolutions. We illustrate how this two-dimensional valley would offer a clear visual representation and test of the $L/E$ dependence, the alignment of the valley quantifying the atmospheric mass-squared difference. Owing to the charge identification capability of the ICAL detector at INO, we always present our results using $\mu^{-}$ and $\mu^{+}$ events separately. Taking into account the statistical fluctuations and systematic errors, and varying oscillation parameters over their currently allowed ranges, we estimate the precision to which atmospheric neutrino oscillation parameters would be determined with the 10-year simulated data at ICAL using our procedure.


title: Bayesian Analysis of the $^{70}$Zn$(d, ^3\!\text{He}) ^{69}$Cu Transfer Reaction

abstract: Transfer reactions provide information about the single-particle nature of nuclear levels. In particular, the differential cross sections from these measurements are sensitive to the angular momentum of the transferred particle and the spectroscopic factor of the populated level. However, the process of extracting these properties is subject to uncertainties, both from experimental and theoretical sources. By integrating the distorted wave Born approximation into a Bayesian model, we propagate these uncertainties through to the spectroscopic factors and orbital angular momentum values. We use previously reported data of the proton pickup reaction $^{70}$Zn$(d, ^3\!\text{He}) ^{69}$Cu as an example. By accounting for uncertainties in the experimental data, optical model parameters, and reaction mechanism, we find that the extracted spectroscopic factors for low lying states of $^{69}$Cu are subject to large, asymmetric uncertainties ranging from $35 \%$ to $108 \%$. Additionally, Bayesian model comparison is employed to assign probabilities to each of the allowed angular momentum transfers. This method confirms the assignments for many states, but suggests that the data for a state lying at $3.70$ MeV is better characterized by an $\ell = 3$ transfer, rather than the previously reported $\ell = 2$.


title: On-demand generation of higher-order Fock states in quantum-dot--cavity systems

abstract: The on-demand preparation of higher-order Fock states is of fundamental importance in quantum information sciences. We propose and compare different protocols to generate higher-order Fock states in solid state quantum-dot--cavity systems. The protocols make use of a series of laser pulses to excite the quantum dot exciton and off-resonant pulses to control the detuning between dot and cavity. Our theoretical studies include dot and cavity loss processes as well as the pure-dephasing type coupling to longitudinal acoustic phonons in a numerically complete fashion. By going beyond the two-level approximation for quantum dots, we study the impact of a finite exchange splitting, the impact of a higher energetic exciton state, and an excitation with linearly polarized laser pulses leading to detrimental occupations of the biexciton state. We predict that under realistic conditions, a protocol which keeps the cavity at resonance with the quantum dot until the desired target state is reached is able to deliver fidelities to the Fock state $| 5\rangle$ well above $40\,\%$.


title: Critical behavior for point monopole and dipole electric impurities in uniformily and uniaxially strained graphene

abstract: We revisit the problem of bound states in graphene under the influence of point electric monopole and dipole impurity potentials extended to the case in which the membrane of this material is uniformly and uniaxially strained, which leads to a redefinition of the charge and dipole moment, respectively. By considering an anisotropic Fermi velocity, we analytically solve the resulting Dirac equation for each potential. We observe that the effect of the anisotropy is to promote or inhibit the critical behavior known to occur for each kind of impurity, depending on the direction along which strain is applied: both the atomic collapse, for the monopole impurity, and the emergence of cascades of infinitely many bound states with a universal Efimov-like scaling, for the dipole impurity, are phenomena that occur under less or more restrictive conditions due to strain.


title: On the monotonicity of the critical time in the Constrained-degree percolation model

abstract: The Constrained-degree percolation model was introduced in [B.N.B. de Lima, R. Sanchis, D.C. dos Santos, V. Sidoravicius, and R. Teodoro, Stoch. Process. Appl. (2020)], where it was proven that this model has a non-trivial phase transition on a square lattice. We study the Constrained-degree percolation model on the $d$-dimensional hypercubic lattice ($\mathbb{Z}^d$) and, via numerical simulations, found evidence that the critical time $t_{c}^{d}(k)$ is monotonous not increasing in the constrained $k$ if $d=3,4$, like it is when $d=2$. We verify that the lowest constrained value $k$ such that the system exhibits a phase transition is $k=3$ and that the correlation critical exponent $\nu$ for the Constrained-degree percolation model and ordinary Bernoulli percolation are the same.


title: A new proof of the growth rate of the solvable Baumslag-Solitar groups

abstract: We exhibit a regular language of geodesics for a large set of elements of $BS(1,n)$ and show that the growth rate of this language is the growth rate of the group. This provides a straightforward calculation of the growth rate of $BS(1,n)$, which was initially computed by Collins, Edjvet and Gill in [5]. Our methods are based on those we develop in [8] to show that $BS(1,n)$ has a positive density of elements of positive, negative and zero conjugation curvature, as introduced by Bar-Natan, Duchin and Kropholler in [1].


title: Perturbed Sachdev-Ye-Kitaev model: a polaron in the hyperbolic plane

abstract: We study the SYK$_4$ model with a weak SYK$_2$ term of magnitude $\Gamma$ beyond the simplest perturbative limit considered previously. For intermediate values of the perturbation strength, $J/N \ll \Gamma \ll J/\sqrt{N}$, fluctuations of the Schwarzian mode are suppressed, and the SYK$_4$ mean-field solution remains valid beyond the timescale $t_0 \sim N/J$ up to $t_* \sim J/\Gamma^2$. Out-of-time-order correlation function displays at short time intervals exponential growth with maximal Lyapunov exponent $2\pi T$, but its prefactor scales as $T$ at low temperatures $T \leq \Gamma$.


title: Smooth Adversarial Training

abstract: It is commonly believed that networks cannot be both accurate and robust, that gaining robustness means losing accuracy. It is also generally believed that, unless making networks larger, network architectural elements would otherwise matter little in improving adversarial robustness. Here we present evidence to challenge these common beliefs by a careful study about adversarial training. Our key observation is that the widely-used ReLU activation function significantly weakens adversarial training due to its non-smooth nature. Hence we propose smooth adversarial training (SAT), in which we replace ReLU with its smooth approximations to strengthen adversarial training. The purpose of smooth activation functions in SAT is to allow it to find harder adversarial examples and compute better gradient updates during adversarial training. Compared to standard adversarial training, SAT improves adversarial robustness for "free", i.e., no drop in accuracy and no increase in computational cost. For example, without introducing additional computations, SAT significantly enhances ResNet-50's robustness from 33.0% to 42.3%, while also improving accuracy by 0.9% on ImageNet. SAT also works well with larger networks: it helps EfficientNet-L1 to achieve 82.2% accuracy and 58.6% robustness on ImageNet, outperforming the previous state-of-the-art defense by 9.5% for accuracy and 11.6% for robustness. Models are available at https://github.com/cihangxie/SmoothAdversarialTraining.


title: Auditory streaming emerges from fast excitation and slow delayed inhibition

abstract: In the auditory streaming paradigm alternating sequences of pure tones can be perceived as a single galloping rhythm (integration) or as two sequences with separated low and high tones (segregation). Although studied for decades, the neural mechanisms underlining this perceptual grouping of sound remains a mystery. With the aim of identifying a plausible minimal neural circuit that captures this phenomenon, we propose a firing rate model with two periodically forced neural populations coupled by fast direct excitation and slow delayed inhibition. By analyzing the model in a non-smooth, slow-fast regime we analytically prove the existence of a rich repertoire of dynamical states and of their parameter dependent transitions. We impose plausible parameter restrictions and link all states with perceptual interpretations. Regions of stimulus parameters occupied by states linked with each percept matches those found in behavioral experiments. Our model suggests that slow inhibition masks the perception of subsequent tones during segregation (forward masking), while fast excitation enables integration for large pitch differences between the two tones.


title: Between-Domain Instance Transition Via the Process of Gibbs Sampling in RBM

abstract: In this paper, we present a new idea for Transfer Learning (TL) based on Gibbs Sampling. Gibbs sampling is an algorithm in which instances are likely to transfer to a new state with a higher possibility with respect to a probability distribution. We find that such an algorithm can be employed to transfer instances between domains. Restricted Boltzmann Machine (RBM) is an energy based model that is very feasible for being trained to represent a data distribution and also for performing Gibbs sampling. We used RBM to capture data distribution of the source domain and use it in order to cast target instances into new data with a distribution similar to the distribution of source data. Using datasets that are commonly used for evaluation of TL methods, we show that our method can successfully enhance target classification by a considerable ratio. Additionally, the proposed method has the advantage over common DA methods that it needs no target data during the process of training of models.


title: Anderson Acceleration for Nonconvex ADMM Based on Douglas-Rachford Splitting

abstract: The alternating direction multiplier method (ADMM) is widely used in computer graphics for solving optimization problems that can be nonsmooth and nonconvex. It converges quickly to an approximate solution, but can take a long time to converge to a solution of high-accuracy. Previously, Anderson acceleration has been applied to ADMM, by treating it as a fixed-point iteration for the concatenation of the dual variables and a subset of the primal variables. In this paper, we note that the equivalence between ADMM and Douglas-Rachford splitting reveals that ADMM is in fact a fixed-point iteration in a lower-dimensional space. By applying Anderson acceleration to such lower-dimensional fixed-point iteration, we obtain a more effective approach for accelerating ADMM. We analyze the convergence of the proposed acceleration method on nonconvex problems, and verify its effectiveness on a variety of computer graphics problems including geometry processing and physical simulation.


title: Graph Convolutional Neural Networks for analysis of EEG signals, BCI application

abstract: Decoding brain signals has gained many attention and has found much applications in recent years such as Brain Computer Interfaces, communicating with controlling external devices using the user's intentions, occupies an emerging field with the potential of changing the world, with diverse applications from rehabilitation to human augmentation. This being said brain signal analysis, EEG brain signal analysis in particular, is a challenging task. With the advances and achievements in the field of deep learning in problem solving with using only raw data, few attempts has been carried in recent years, to apply deep learning to tackle EEG among other types of brain signals. In this study, we propose a novel loss function, called DeepCSP to extend the classical Common Spatial Patterns to a non linear, differentiable module to serve as the loss function to enforce linearly separable latent representations of EEG signals belonging to different classes in an end to end manner on raw signals without the need to perform extensive feature engineering. With recent generalizations of deep learning methods to work on arbitrarily structured graphs and the introduced loss we have proposed two light weight models to decode EEG signals and carried experiments to show their performance.


title: Relatively free algebras of finite rank

abstract: Let $\mathbb{K}$ be a field of characteristic zero and $B=B_0+B_1$ a finite dimensional associative superalgebra. In this paper we investigate the polynomial identities of the relatively free algebras of finite rank of the variety $\mathfrak V$ defined by the Grassmann envelope of $B$. We also consider the $k$-th Grassmann Envelope of $B$, $G^{(k)}(B)$, constructed with the $k$-generated Grassmann algebra, instead of the infinite dimensional Grassmann algebra. We specialize our studies for the algebra $UT_2(G)$ and $UT_2(G^{(k)})$, which can be seen as the Grassmann envelope and $k$-th Grassmann envelope, respectively, of the superalgebra $UT_2(\mathbb{K}[u])$, where $u^2=1$.


title: Moser's Method and Conservative Extensions of Diffeomorphisms

abstract: This paper shall be concerned with three main results. After a brief recollection of basic symplectic geometry, we prove using Moser's homotopy method a special case of the Strong Darboux Theorem found, for instance, in Theorem 21.1.6 of [Hor]. Next, we'll prove two conservative extension results for a diffeomorphism on a circle. One uses Moser's homotopy method but loses a degree of regularity. The other uses the method of generating functions as found in [BCW] and [BGV]. Finally, we'll prove a conservative extension result for a "diffeomorphism" defined on the boundary of $(0, 1)^2$ and use the techniques developed there and by [M] to prove an ambient Dacarogna-Moser Theorem.


title: Decomposable Pauli diagonal maps and Tensor Squares of Qubit Maps

abstract: It is a well-known result due to E. St{\o}rmer that every positive qubit map is decomposable into a sum of a completely positive map and a completely copositive map. Here, we generalize this result to tensor squares of qubit maps. Specifically, we show that any positive tensor product of a qubit map with itself is decomposable. This solves a recent conjecture by S. Fillipov and K. Magadov. We contrast this result with examples of non-decomposable positive maps arising as the tensor product of two distinct qubit maps or as the tensor square of a decomposable map from a qubit to a ququart. To show our main result, we reduce the problem to Pauli diagonal maps. We then characterize the cone of decomposable ququart Pauli diagonal maps by determining all 252 extremal rays of ququart Pauli diagonal maps that are both completely positive and completely copositive. These extremal rays split into three disjoint orbits under a natural symmetry group, and two of these orbits contain only entanglement breaking maps. Finally, we develop a general combinatorial method to determine the extremal rays of Pauli diagonal maps that are both completely positive and completely copositive between multi-qubit systems using the ordered spectra of their Choi matrices. Classifying these extremal rays beyond ququarts is left as an open problem.

categories: quant-ph math-ph math.FA math.MP math.OA


title: Machine Learning-Assisted UAV Operations with UTM: Requirements, Challenges, and Solutions

abstract: Unmanned aerial vehicles (UAVs) are emerging in commercial spaces and will support many applications and services, such as smart agriculture, dynamic network deployment, and network coverage extension, surveillance and security. The unmanned aircraft system (UAS) traffic management (UTM) provides a framework for safe UAV operation integrating UAV controllers and central data bases via a communications network. This paper discusses the challenges and opportunities for machine learning (ML) for effectively providing critical UTM services. We introduce the four pillars of UTM---operation planning, situational awareness, status and advisors and security---and discuss the main services, specific opportunities for ML and the ongoing research. We conclude that the multi-faceted operating environment and operational parameters will benefit from collected data and data-driven algorithms, as well as online learning to face new UAV operation situations.


title: The thermodynamics of large-N QCD and the nature of metastable phases

abstract: In the limit of a large number of colors (N), both Yang-Mills and quantum chromodynamics are expected to have a first-order phase transition separating a confined hadronic phase and a deconfined plasma phase. One aspect of this separation is that at large N, one can unambiguously identify a plasma regime that is strongly coupled. The existence of a first-order transition suggests that the hadronic phase can be superheated and the plasma phase supercooled. The supercooled deconfined plasma present at large N, if it exists, has the remarkable property that it has negative absolute pressure -- i.e. a pressure below that of the vacuum. For energy densities of order unity in a 1/N expansion but beyond the endpoint of the hadronic superheated phase, a description of homogeneous matter composed of ordinary hadrons with masses of order unity in a 1/N expansion can exist, and acts as though it has a temperature of $T_H$ in order unity. However, the connection between the canonical and microcanonical descriptions breaks down and the system cannot fully equilibrate as $N \rightarrow \infty$. Rather, in a hadronic description, energy is pushed to hadrons with masses that are arbitrarily large. The thermodynamic limit of large volumes becomes subtle for such systems: the energy density is no longer intensive. These conclusions follow provided that standard large N scaling rules hold, the system at large N undergoes a generic first-order phase transition between the hadronic and plasma phases and that the mesons and glueballs follow a Hagedorn-type spectrum.


title: TOI-1728b: The Habitable-zone Planet Finder confirms a warm super Neptune orbiting an M dwarf host

abstract: We confirm the planetary nature of TOI-1728b using a combination of ground-based photometry, near-infrared Doppler velocimetry and spectroscopy with the Habitable-zone Planet Finder.TOI-1728 is an old, inactive M0 star with \teff{} $= 3980^{+31}_{-32}$ K, which hosts a transiting super Neptune at an orbital period of $\sim$ 3.49 days. Joint fitting of the radial velocities and TESS and ground-based transits yields a planetary radius of $5.05_{-0.17}^{+0.16}$ R$_{\oplus}$, mass $26.78_{-5.13}^{+5.43}$ M$_{\oplus}$ and eccentricity $0.057_{-0.039}^{+0.054}$. We estimate the stellar properties, and perform a search for He 10830 \AA absorption during the transit of this planet and claim a null detection with an upper limit of 1.1$\%$ with 90\% confidence. A deeper level of He 10830 \AA ~ absorption has been detected in the planet atmosphere of GJ 3470b, a comparable gaseous planet. TOI-1728b is the largest super Neptune -- the intermediate subclass of planets between Neptune and the more massive gas-giant planets -- discovered around an M dwarf. With its relatively large mass and radius, TOI-1728 represents a valuable datapoint in the M-dwarf exoplanet mass-radius diagram, bridging the gap between the lighter Neptune-sized planets and the heavier Jovian planets known to orbit M-dwarfs. With a low bulk density of $1.14_{-0.24}^{+0.26}$ g/cm$^3$, and orbiting a bright host star (J $\sim 9.6$, V $\sim 12.4$), TOI-1728b is also a promising candidate for transmission spectroscopy both from the ground and from space, which can be used to constrain planet formation and evolutionary models.


title: Estimating Displaced Populations from Overhead

abstract: We introduce a deep learning approach to perform fine-grained population estimation for displacement camps using high-resolution overhead imagery. We train and evaluate our approach on drone imagery cross-referenced with population data for refugee camps in Cox's Bazar, Bangladesh in 2018 and 2019. Our proposed approach achieves 7.02% mean absolute percent error on sequestered camp imagery. We believe our experiments with real-world displacement camp data constitute an important step towards the development of tools that enable the humanitarian community to effectively and rapidly respond to the global displacement crisis.


title: Tensor Programs II: Neural Tangent Kernel for Any Architecture

abstract: We prove that a randomly initialized neural network of *any architecture* has its Tangent Kernel (NTK) converge to a deterministic limit, as the network widths tend to infinity. We demonstrate how to calculate this limit. In prior literature, the heuristic study of neural network gradients often assumes every weight matrix used in forward propagation is independent from its transpose used in backpropagation (Schoenholz et al. 2017). This is known as the *gradient independence assumption (GIA)*. We identify a commonly satisfied condition, which we call *Simple GIA Check*, such that the NTK limit calculation based on GIA is correct. Conversely, when Simple GIA Check fails, we show GIA can result in wrong answers. Our material here presents the NTK results of Yang (2019a) in a friendly manner and showcases the *tensor programs* technique for understanding wide neural networks. We provide reference implementations of infinite-width NTKs of recurrent neural network, transformer, and batch normalization at https://github.com/thegregyang/NTK4A.


title: The existence and unicity of numerical solution of initial value problems by Walsh polynomials approach

abstract: Chen and Hsiao gave the numerical solution of initial value problems of systems of linear differential equations with constant coefficients by Walsh polynomials approach. This result was improved by G\'at and Toledo for initial value problems of differential equations with variable coefficients on the interval $[0,1[$ and initial value $\xi=0$. In the present paper we discuss the general case while $\xi$ can take any arbitrary value in the interval $[0,1[$. We show the existence and uniform convergence of the numerical solution, as well.


title: Lifted Disjoint Paths with Application in Multiple Object Tracking

abstract: We present an extension to the disjoint paths problem in which additional \emph{lifted} edges are introduced to provide path connectivity priors. We call the resulting optimization problem the lifted disjoint paths problem. We show that this problem is NP-hard by reduction from integer multicommodity flow and 3-SAT. To enable practical global optimization, we propose several classes of linear inequalities that produce a high-quality LP-relaxation. Additionally, we propose efficient cutting plane algorithms for separating the proposed linear inequalities. The lifted disjoint path problem is a natural model for multiple object tracking and allows an elegant mathematical formulation for long range temporal interactions. Lifted edges help to prevent id switches and to re-identify persons. Our lifted disjoint paths tracker achieves nearly optimal assignments with respect to input detections. As a consequence, it leads on all three main benchmarks of the MOT challenge, improving significantly over state-of-the-art.


title: Prediction with Approximated Gaussian Process Dynamical Models

abstract: The modeling and simulation of dynamical systems is a necessary step for many control approaches. Using classical, parameter-based techniques for modeling of modern systems, e.g., soft robotics or human-robot interaction, is often challenging or even infeasible due to the complexity of the system dynamics. In contrast, data-driven approaches need only a minimum of prior knowledge and scale with the complexity of the system. In particular, Gaussian process dynamical models (GPDMs) provide very promising results for the modeling of complex dynamics. However, the control properties of these GP models are just sparsely researched, which leads to a "blackbox" treatment in modeling and control scenarios. In addition, the sampling of GPDMs for prediction purpose respecting their non-parametric nature results in non-Markovian dynamics making the theoretical analysis challenging. In this article, we present approximated GPDMs which are Markov and analyze their control theoretical properties. Among others, the approximated error is analyzed and conditions for boundedness of the trajectories are provided. The outcomes are illustrated with numerical examples that show the power of the approximated models while the the computational time is significantly reduced.


title: Practical Trade-Offs for the Prefix-Sum Problem

abstract: Given an integer array A, the prefix-sum problem is to answer sum(i) queries that return the sum of the elements in A[0..i], knowing that the integers in A can be changed. It is a classic problem in data structure design with a wide range of applications in computing from coding to databases. In this work, we propose and compare several and practical solutions to this problem, showing that new trade-offs between the performance of queries and updates can be achieved on modern hardware.


title: Did You Remember to Test Your Tokens?

abstract: Authentication is a critical security feature for confirming the identity of a system's users, typically implemented with help from frameworks like Spring Security. It is a complex feature which should be robustly tested at all stages of development. Unit testing is an effective technique for fine-grained verification of feature behaviors that is not widely-used to test authentication. Part of the problem is that resources to help developers unit test security features are limited. Most security testing guides recommend test cases in a "black box" or penetration testing perspective. These resources are not easily applicable to developers writing new unit tests, or who want a security-focused perspective on coverage. In this paper, we address these issues by applying a grounded theory-based approach to identify common (unit) test cases for token authentication through analysis of 481 JUnit tests exercising Spring Security-based authentication implementations from 53 open source Java projects. The outcome of this study is a developer-friendly unit testing guide organized as a catalog of 53 test cases for token authentication, representing unique combinations of 17 scenarios, 40 conditions, and 30 expected outcomes learned from the data set in our analysis. We supplement the test guide with common test smells to avoid. To verify the accuracy and usefulness of our testing guide, we sought feedback from selected developers, some of whom authored unit tests in our dataset.


title: STORM: Foundations of End-to-End Empirical Risk Minimization on the Edge

abstract: Empirical risk minimization is perhaps the most influential idea in statistical learning, with applications to nearly all scientific and technical domains in the form of regression and classification models. To analyze massive streaming datasets in distributed computing environments, practitioners increasingly prefer to deploy regression models on edge rather than in the cloud. By keeping data on edge devices, we minimize the energy, communication, and data security risk associated with the model. Although it is equally advantageous to train models at the edge, a common assumption is that the model was originally trained in the cloud, since training typically requires substantial computation and memory. To this end, we propose STORM, an online sketch for empirical risk minimization. STORM compresses a data stream into a tiny array of integer counters. This sketch is sufficient to estimate a variety of surrogate losses over the original dataset. We provide rigorous theoretical analysis and show that STORM can estimate a carefully chosen surrogate loss for the least-squares objective. In an exhaustive experimental comparison for linear regression models on real-world datasets, we find that STORM allows accurate regression models to be trained.


title: Strong Spin-Orbit Torque effect on magnetic defects due to topological surface state electrons in Bi$_{2}$Te$_{3}$

abstract: We investigate the spin-orbit torque exerted on the magnetic moments of the transition-metal impurities Cr, Mn, Fe and Co, embedded in the surface of the topological insulator Bi$_{2}$Te$ _{3} $, in response to an electric field and a consequent electrical current flow in the surface. The multiple scattering problem of electrons off impurity atoms is solved by first-principles calculations within the full-potential relativistic Korringa-Kohn-Rostoker (KKR) Green function method, while the spin-orbit torque calculations are carried out by combining the KKR method with the semiclassical Boltzmann transport equation. We analyze the correlation of the spin-orbit torque to the spin accumulation and spin flux in the defects. We compare the torque on different magnetic impurities and unveil the effect of resonant scattering. In addition, we calculate the resistivity and the Joule heat as a function of the torque in these systems. We predict that the Mn/Bi$_{2}$Te$_{3}$ is optimal among the studied systems.


title: Anomaly Detection using Deep Reconstruction and Forecasting for Autonomous Systems

abstract: We propose self-supervised deep algorithms to detect anomalies in heterogeneous autonomous systems using frontal camera video and IMU readings. Given that the video and IMU data are not synchronized, each of them are analyzed separately. The vision-based system, which utilizes a conditional GAN, analyzes immediate-past three frames and attempts to predict the next frame. The frame is classified as either an anomalous case or a normal case based on the degree of difference estimated using the prediction error and a threshold. The IMU-based system utilizes two approaches to classify the timestamps; the first being an LSTM autoencoder which reconstructs three consecutive IMU vectors and the second being an LSTM forecaster which is utilized to predict the next vector using the previous three IMU vectors. Based on the reconstruction error, the prediction error, and a threshold, the timestamp is classified as either an anomalous case or a normal case. The composition of algorithms won runners up at the IEEE Signal Processing Cup anomaly detection challenge 2020. In the competition dataset of camera frames consisting of both normal and anomalous cases, we achieve a test accuracy of 94% and an F1-score of 0.95. Furthermore, we achieve an accuracy of 100% on a test set containing normal IMU data, and an F1-score of 0.98 on the test set of abnormal IMU data.


title: "On-the-fly" calculation of the Vibrational Sum-frequency Generation Spectrum at the Air-water Interface

abstract: In the present work, we provide an electronic structure based method for the "on-the-fly" determination of vibrational sum frequency generation (v-SFG) spectra. The predictive power of this scheme is demonstrated at the air-water interface. While the instantaneous fluctuations in dipole moment are obtained using the maximally localized Wannier functions, the fluctuations in polarizability are approximated to be proportional to the second moment of Wannier functions. The spectrum henceforth obtained captures the signatures of hydrogen bond stretching, bending, as well as low-frequency librational modes.


title: On quantum toroidal algebra of type $A_1$

abstract: In this paper we introduce a new quantum algebra which specializes to the $2$-toroidal Lie algebra of type $A_1$. We prove that this quantum toroidal algebra has a natural triangular decomposition, a (topological) Hopf algebra structure and a vertex operator realization.


title: Ising model and s-embeddings of planar graphs

abstract: We discuss the notion of s-embeddings $\mathcal{S}=\mathcal{S}_\mathcal{X}$ of planar graphs carrying a nearest-neighbor Ising model. The construction of $\mathcal{S}_\mathcal{X}$ is based upon a choice of a global complex-valued solution $\mathcal{X}$ of the propagation equation for Kadanoff-Ceva fermions. Each choice of $\mathcal{X}$ provides an interpretation of all other fermionic observables as s-holomorphic functions on $\mathcal{S}_\mathcal{X}$. We set up a general framework for the analysis of such functions on s-embeddings $\mathcal{S}^\delta$ with $\delta\to 0$. Throughout this analysis, a key role is played by the functions $\mathcal{Q}^\delta$ associated with $\mathcal{S}^\delta$, the so-called origami maps in the bipartite dimer model terminology. In particular, we give an interpretation of the mean curvature of the limit of discrete surfaces $(\mathcal{S}^\delta;\mathcal{Q}^\delta)$ viewed in the Minkowski space $\mathbb R^{2,1}$ as the mass in the Dirac equation describing the continuous limit of the model. We then focus on the simplest situation when $\mathcal{S}^\delta$ have uniformly bounded lengths/angles and $\mathcal{Q}^\delta=O(\delta)$; as a particular case this includes all critical Ising models on doubly periodic graphs via their canonical s-embeddings. In this setup we prove RSW-type crossing estimates for the random cluster representation of the model and the convergence of basic fermionic observables. The proof relies upon a new strategy as compared to the already existing literature, it also provides a quantitative estimate on the speed of convergence.


title: Learning compositional functions via multiplicative weight updates

abstract: Compositionality is a basic structural feature of both biological and artificial neural networks. Learning compositional functions via gradient descent incurs well known problems like vanishing and exploding gradients, making careful learning rate tuning essential for real-world applications. This paper proves that multiplicative weight updates satisfy a descent lemma tailored to compositional functions. Based on this lemma, we derive Madam -- a multiplicative version of the Adam optimiser -- and show that it can train state of the art neural network architectures without learning rate tuning. We further show that Madam is easily adapted to train natively compressed neural networks by representing their weights in a logarithmic number system. We conclude by drawing connections between multiplicative weight updates and recent findings about synapses in biology.


title: An Uncertainty Management Framework for Integrated Gas-Electric Energy Systems

abstract: In many parts of the world, electric power systems have seen a significant shift towards generation from renewable energy and natural gas. Because of their ability to flexibly adjust power generation in real time, gas-fired power plants are frequently seen as the perfect partner for variable renewable generation. However, this reliance on gas generation increases interdependence and propagates uncertainty between power grids and gas pipelines, and brings coordination and uncertainty management challenges. To address these issues, we propose an uncertainty management framework for uncertain, but bounded gas consumption by gas-fired power plants. The admissible ranges are computed based on a joint optimization problem for the combined gas and electricity networks, which involves chance-constrained scheduling for the electric grid and a novel robust optimization formulation for the natural gas network. This formulation ensures feasibility of the integrated system with a high probability, while providing a tractable numerical formulation. A key advance with respect to existing methods is that our method is based on a physically accurate, validated model for transient gas pipeline flows. Our case study benchmarks our proposed formulation against methods that ignore how reserve activation impacts the fuel use of gas power plants, and only consider predetermined gas consumption. The results demonstrate the importance of considering uncertainty to avoid operating constraint violations and curtailment of gas to the generators.


title: A new class of minimal asymptotic bases

abstract: A set $A$ of nonnegative integers is an asymptotic basis of order $h$ if every sufficiently large integer can be represented as the sum of $h$ not necessarily distinct elements of $A$. The asymptotic basis $A$ is minimal if removing any element of $A$ destroys every representation of infinitely many integers, and so $A\setminus \{a\}$ is not an asymptotic basis of order $h$ for all $a\in A$. In this paper, a new class of minimal asymptotic bases is constructed.


title: Dynamically Mitigating Data Discrepancy with Balanced Focal Loss for Replay Attack Detection

abstract: It becomes urgent to design effective anti-spoofing algorithms for vulnerable automatic speaker verification systems due to the advancement of high-quality playback devices. Current studies mainly treat anti-spoofing as a binary classification problem between bonafide and spoofed utterances, while lack of indistinguishable samples makes it difficult to train a robust spoofing detector. In this paper, we argue that for anti-spoofing, it needs more attention for indistinguishable samples over easily-classified ones in the modeling process, to make correct discrimination a top priority. Therefore, to mitigate the data discrepancy between training and inference, we propose D3M, to leverage a balanced focal loss function as the training objective to dynamically scale the loss based on the traits of the sample itself. Besides, in the experiments, we select three kinds of features that contain both magnitude-based and phase-based information to form complementary and informative features. Experimental results on the ASVspoof2019 dataset demonstrate the superiority of the proposed methods by comparison between our systems and top-performing ones. Systems trained with the balanced focal loss perform significantly better than conventional cross-entropy loss. With complementary features, our fusion system with only three kinds of features outperforms other systems containing five or more complex single models by 22.5% for min-tDCF and 7% for EER, achieving a min-tDCF and an EER of 0.0124 and 0.55% respectively. Furthermore, we present and discuss the evaluation results on real replay data apart from the simulated ASVspoof2019 data, indicating that research for anti-spoofing still has a long way to go. Source code, analysis data, and other details are publicly available at https://github.com/asvspoof/D3M.


title: ALMA observations of transient heating in a solar active region

abstract: We aim to investigate the temperature enhancements and formation heights of impulsive heating phenomena in solar active-regions such as Ellerman bombs (EBs), ultraviolet bursts (UVBs), and flaring active-region fibrils (FAFs) using interferometric observations in the millimeter (mm) continuum provided by the Atacama Large Millimeter/submillimeter Array (ALMA). We examined 3 mm signatures of heating events identified in Solar Dynamics Observatory (SDO) observations of an active region and compared the results with synthetic spectra from a 3D radiative magnetohydrodynamic simulation. We estimated the contribution from the corona to the mm brightness using differential emission measure analysis. We report the null detection of EBs in the 3 mm continuum at $\sim1.2$" spatial resolution, which is evidence that they are sub-canopy events that do not significantly contribute to heating the upper chromosphere. In contrast, we find the active region to be populated with multiple compact, bright, flickering mm bursts -- reminiscent of UVBs. The high brightness temperatures of up to $\sim14200$ K in some events have a significant contribution (up to $\sim$7%) from the corona. We also detect FAF-like events in the 3 mm continuum that show rapid motions of $>10000\,$K plasma launched with high plane-of-sky velocities ($37-340\rm\,km\,s^{-1}$) from bright kernels. The mm FAFs are the brightest class of warm canopy fibrils that connect magnetic regions of opposite polarities. The simulation confirms that ALMA should be able to detect the mm counterparts of UVBs and small flares and thus provide a complementary diagnostic for localized heating in the solar chromosphere.


title: Habit learning supported by efficiently controlled network dynamics in naive macaque monkeys

abstract: Primates display a marked ability to learn habits in uncertain and dynamic environments. The associated perceptions and actions of such habits engage distributed neural circuits. Yet, precisely how such circuits support the computations necessary for habit learning remain far from understood. Here we construct a formal theory of network energetics to account for how changes in brain state produce changes in sequential behavior. We exercise the theory in the context of multi-unit recordings spanning the caudate nucleus, prefrontal cortex, and frontal eyefields of female macaque monkeys engaged in 60-180 sessions of a free scan task that induces motor habits. The theory relies on the determination of effective connectivity between recording channels, and on the stipulation that a brain state is taken to be the trial-specific firing rate across those channels. The theory then predicts how much energy will be required to transition from one state into another, given the constraint that activity can spread solely through effective connections. Consistent with the theory's predictions, we observed smaller energy requirements for transitions between more similar and more complex trial saccade patterns, and for sessions characterized by less entropic selection of saccade patterns. Using a virtual lesioning approach, we demonstrate the resilience of the observed relationships between minimum control energy and behavior to significant disruptions in the inferred effective connectivity. Our theoretically principled approach to the study of habit learning paves the way for future efforts examining how behavior arises from changing patterns of activity in distributed neural circuitry.


title: Multimarginal Wasserstein Barycenter for Stain Normalization and Augmentation

abstract: Variations in hematoxylin and eosin (H&E) stained images (due to clinical lab protocols, scanners, etc) directly impact the quality and accuracy of clinical diagnosis, and hence it is important to control for these variations for a reliable diagnosis. In this work, we present a new approach based on the multimarginal Wasserstein barycenter to normalize and augment H&E stained images given one or more references. Specifically, we provide a mathematically robust way of naturally incorporating additional images as intermediate references to drive stain normalization and augmentation simultaneously. The presented approach showed superior results quantitatively and qualitatively as compared to state-of-the-art methods for stain normalization. We further validated our stain normalization and augmentations in the nuclei segmentation task on a publicly available dataset, achieving state-of-the-art results against competing approaches.