title: Effects of centrality fluctuation and deuteron formation on proton number cumulant in Au+Au collisions at $\sqrt{s_\mathrm{NN}}$ = 3 GeV from JAM model

abstract: We studied the effects of centrality fluctuation and deuteron formation on the cumulants ($C_n$) and correlation functions ($\kappa_n$) of protons up to sixth order in most central ($b<3$ fm) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$ = 3 GeV from a microscopic transport model (JAM). The results are presented as a function of rapidity acceptance within transverse momentum $0.4<p_{T}<2 $ GeV/$c$. We compared the results obtained by centrality bin width correction (CBWC) using charged reference particle multiplicity with CBWC done using impact parameter bins. It was found that at low energies the centrality resolution for determining the collision centrality using charged particle multiplicities is not good enough to reduce the initial volume fluctuations effect for higher-order cumulant analysis. New methods need to be developed to classify events with high centrality resolution for heavy-ion collisions at low energies. Finally, we observed that the formation of deuteron will suppress the higher-order cumulants and correlation functions of protons and is found to be similar to the efficiency effect. This work can serve as a noncritical baseline for the QCD critical point search at the high baryon density region.


title: Machine Learning Towards Enabling Spectrum-as-a-Service Dynamic Sharing

abstract: The growth in wireless broadband users, devices, and novel applications has led to a significant increase in the demand for new radio frequency spectrum. This is expected to grow even further given the projection that the global traffic per year will reach 4.8 zettabytes by 2022. Moreover, it is projected that the number of Internet users will reach 4.8 billion and the number of connected devices will be close 28.5 billion devices. However, due to the spectrum being mostly allocated and divided, providing more spectrum to expand existing services or offer new ones has become more challenging. To address this, spectrum sharing has been proposed as a potential solution to improve spectrum utilization efficiency. Adopting effective and efficient spectrum sharing mechanisms is in itself a challenging task given the multitude of levels and techniques that can be integrated to enable it. To that end, this paper provides an overview of the different spectrum sharing levels and techniques that have been proposed in the literature. Moreover, it discusses the potential of adopting dynamic sharing mechanisms by offering Spectrum-as-a-Service architecture. Furthermore, it describes the potential role of machine learning models in facilitating the automated and efficient dynamic sharing of the spectrum and offering Spectrum-as-a-Service.


title: A Note On Inference for the Mixed Fractional Ornstein-Uhlenbeck Process with Drift

abstract: This paper is devoted to parameter estimation of the mixed fractional Ornstein-Uhlenbeck process with a drift. Large sample asymptotical properties of the Maximum Likelihood Estimator is deduced using the Laplace transform computations or the Cameron-Martin formula with extra part from \cite{CK19}


title: MOSGA: Modular Open-Source Genome Annotator

abstract: The generation of high-quality assemblies, even for large eukaryotic genomes, has become a routine task for many biologists thanks to recent advances in sequencing technologies. However, the annotation of these assemblies - a crucial step towards unlocking the biology of the organism of interest - has remained a complex challenge that often requires advanced bioinformatics expertise. Here we present MOSGA, a genome annotation framework for eukaryotic genomes with a user-friendly web-interface that generates and integrates annotations from various tools. The aggregated results can be analyzed with a fully integrated genome browser and are provided in a format ready for submission to NCBI. MOSGA is built on a portable, customizable, and easily extendible Snakemake backend, and thus, can be tailored to a wide range of users and projects. We provide MOSGA as a publicly free available web service at https://mosga.mathematik.uni-marburg.de and as a docker container at registry.gitlab.com/mosga/mosga:latest. Source code can be found at https://gitlab.com/mosga/mosga


title: An integrative smoothed particle hydrodynamics framework for modeling cardiac function

abstract: Mathematical modeling of cardiac function can provide augmented simulation-based diagnosis tool for complementing and extending human understanding of cardiac diseases which represent the most common cause of worldwide death. As the realistic starting-point for developing an unified meshless approach for total heart modeling, herein we propose an integrative smoothed particle hydrodynamics (SPH) framework for addressing the simulation of the principle aspects of cardiac function, including cardiac electrophysiology, passive mechanical response and electromechanical coupling. To that end, several algorithms, e.g., splitting reaction-by-reaction method combined with quasi-steady-state (QSS) solver , anisotropic SPH-diffusion discretization and total Lagrangian SPH formulation, are introduced and exploited for dealing with the fundamental challenges of developing integrative SPH framework for simulating cardiac function, namely, (i) the correct capturing of the stiff dynamics of the transmembrane potential and the gating variables , (ii) the stable predicting of the large deformations and the strongly anisotropic behavior of the myocardium, and (iii) the proper coupling of electrophysiology and tissue mechanics for electromechanical feedback. A set of numerical examples demonstrate the effectiveness and robustness of the present SPH framework, and render it a potential and powerful alternative that can augment current lines of total cardiac modeling and clinical applications.


title: Cohomology and conformal derivations of BiHom-Lie conformal superalgebras

abstract: In this paper, we introduce the notion of BiHom-Lie conformal superalgebras. We develop its representation theory and define the cohomology group with coefficients in a module. Finally, we introduce conformal derivations of BiHom-Lie conformal superalgebras and study some of their properties.


title: Electoral Accountability and Selection with Personalized Information Aggregation

abstract: We study a model of electoral accountability and selection whereby heterogeneous voters aggregate incumbent politician's performance data into personalized signals through paying limited attention. Extreme voters' signals exhibit an own-party bias, which hampers their ability to discern the good and bad performances of the incumbent. While this effect alone would undermine electoral accountability and selection, there is a countervailing effect stemming from partisan disagreement, which makes the centrist voter more likely to be pivotal. In case the latter's unbiased signal is very informative about the incumbent's performance, the combined effect on electoral accountability and selection can actually be a positive one. For this reason, factors that carry a negative connotation in every political discourse -- such as increasing mass polarization and shrinking attention span -- have ambiguous accountability and selection effects in general. Correlating voters' signals, if done appropriately, unambiguously improves electoral accountability and selection and, hence, voter welfare.


title: Development and comparison of spectral algorithms for numerical modeling of the quasi-static mechanical behavior of inhomogeneous materials

abstract: In the current work, a number of algorithms are developed and compared for the numerical solution of periodic (quasi-static) linear elastic mechanical boundary-value problems (BVPs) based on two different discretizations of Fourier series. The first is standard and based on the trapezoidal approximation of the Fourier mode integral, resulting in trapezoidal discretization (TD) of the truncated Fourier series. Less standard is the second discretization based on piecewise-constant approximation of the Fourier mode integrand, yielding a piecewise-constant discretization (PCD) of this series. Employing these, fixed-point algorithms are formulated via Green-function preconditioning (GFP) and finite-difference discretization (of differential operators; FDD). In particular, in the context of PCD, this includes an algorithm based on the so-called "discrete Green operator" (DGO) recently introduced by Eloh et al. (2019), which employs GFP, but not FDD. For computational comparisons, the (classic) benchmark case of a cubic inclusion embedded in a matrix (e.g., Suquet, 1997; Willot, 2015) is employed. Both discontinuous and smooth transitions in elastic stiffness at the matrix-inclusion (MI) interface are considered. In the context of both TD and PCD, a number of GFP- and FDD-based algorithms are developed. Among these, one based on so-called averaged-forward-backward-differencing (AFB) is shown to result in the greatest improvement in convergence rate. As it turns out, AFB is equivalent to the "rotated scheme" (R) of Willot (2015) in the context of TD. In the context of PCD, comparison of the performance and convergence behavior of AFB/R- and DGO-based algorithms shows that the former is more efficient than the latter for larger phase contrasts.


title: Polyethylene Identification in ocean water samples by means of 50 keV energy electron beam

abstract: The study presented hereafter shows a new methodology to reveal traces of polyethylene (the most common microplastic particles, known as a structure of C_2 H_4) in a sample of ocean water by the irradiation of a 50 keV, 1 uA electron beam. This is performed by analyzing the photon (produced by the electrons in water ) fluxes and spectra (i.e. fluxes as a function of photon energy) at different types of contaminated water with an adequate device and in particular looking at the peculiar interactions of electrons/photons with the potential abnormal atomic hydrogen (H), oxygen (O), carbon (C), phosphorus (P) compositions present in the water, as a function of living and not living organic organisms with a PO4 group RNA/DNA strands in a cluster configuration through a volumetric cells grid.


title: An equation of state for expanded metals

abstract: We present a model equation of states for expanded metals, which contains a pressure term due to a screened-Coulomb potential with a screening parameter reflecting the Mott-Anderson metal-to-nonmetal transition. As anticipated almost 80 years ago by Zel'dovich and Landau, this term gives rise to a second coexistence line in the phase diagram, indicating a phase separation between a metallic and a nonmetallic liquid.


title: Magnetization of the spin-1/2 Heisenberg antiferromagnet on the triangular lattice

abstract: After decades of debate, now there is a rough consensus that at zero temperature the spin-$1/2$ Heisenberg antiferromagnet on the triangular lattice is three-sublattice $120^\circ$ magnetically ordered, in contrast to a quantum spin liquid as originally proposed. However, there remains considerable discrepancy in the magnetization reported among various methods. To resolve this issue, in this work we revisit this model by the tensor-network state algorithm. The ground-state energy per bond $E_b$ and magnetization per spin $M_0$ in the thermodynamic limit are obtained with high precision. The former is estimated to be $E_b = -0.18334(10)$. This value agrees well with that from the series expansion. The three-sublattice magnetic order is firmly confirmed and the magnetization is determined as $M_0 = 0.161(5)$. It is about $32\%$ of its classical value and slightly below the lower bound from the series expansion. In comparison with the best estimated value by Monte Carlo and density-matrix renormalization group, our result is about $20\%$ smaller. This magnetic order is consistent with further analysis of the three-body correlation. Our work thus provides new benchmark results for this prototypical model.


title: A class of solutions of the two-dimensional Toda lattice equation

abstract: A method is proposed to systematically generate solutions of the two-dimensional Toda lattice equation in terms of previously known solutions $\phi\left(x,y\right)$ of the two-dimensional Laplace's equation. The two-dimensional solution of Nakamura's [J. Phys. Soc. Jpn. \textbf{52}, 380 (1983)] is shown to correspond to one particular choice of $\phi\left(x,y\right)$.

categories: nlin.SI math-ph math.MP physics.plasm-ph


title: Safe-by-Design Control for Euler-Lagrange Systems

abstract: Safety-critical control is characterized as ensuring constraint satisfaction for a given dynamical system. Recent developments in zeroing control barrier functions (ZCBFs) have provided a framework for ensuring safety of a superlevel set of a single constraint function. Euler-Lagrange systems represent many real-world systems including robots and vehicles, which must abide by safety-regulations, especially for use in human-occupied environments. These safety regulations include state constraints (position and velocity) and input constraints that must be respected at all times. ZCBFs are valuable for satisfying system constraints for general nonlinear systems, however their construction to satisfy state and input constraints is not straightforward. Furthermore, the existing barrier function methods do not address the multiple state constraints that are required for safety of Euler-Lagrange systems. In this paper, we propose a methodology to construct multiple, non-conflicting control barrier functions for Euler-Lagrange systems subject to input constraints to satisfy safety regulations, while concurrently taking into account robustness margins and sampling-time effects. The proposed approach consists of a sampled-data controller and an algorithm for barrier function construction to enforce safety (i.e satisfy position and velocity constraints). The proposed method is validated in simulation on a 2-DOF planar manipulator.


title: Federated Classification using Parsimonious Functions in Reproducing Kernel Hilbert Spaces

abstract: Federated learning forms a global model using data collected from a federation agent. This type of learning has two main challenges: the agents generally don't collect data over the same distribution, and the agents have limited capabilities of storing and transmitting data. Therefore, it is impractical for each agent to send the entire data over the network. Instead, each agent must form a local model and decide what information is fundamental to the learning problem, which will be sent to a central unit. The central unit can then form the global model using only the information received from the agents. We propose a method that tackles these challenges. First each agent forms a local model using a low complexity reproducing kernel Hilbert space representation. From the model the agents identify the fundamental samples which are sent to the central unit. The fundamental samples are obtained by solving the dual problem. The central unit then forms the global model. We show that the solution of the federated learner converges to that of the centralized learner asymptotically as the sample size increases. The performance of the proposed algorithm is evaluated using experiments with both simulated data and real data sets from an activity recognition task, for which the data is collected from a wearable device. The experimentation results show that the accuracy of our method converges to that of a centralized learner with increasing sample size.


title: Hybrid stars with sequential phase transitions: the emergence of the g$_2$ mode

abstract: Neutron stars are the densest objects in the Universe, with $M \sim 1.4 M_{\odot}$ and $R \sim 12$ km, and the equation of state associated to their internal composition is still unknown. The extreme conditions to which matter is subjected inside neutron stars could lead to a phase transition in their inner cores, giving rise to a hybrid compact object. The observation of $2M_{\odot}$ binary pulsars (PSR~J1614-2230, PSR~J0343$+$0432 and PSR~J0740$+$6620) strongly constraints theoretical models of the equation of state. Moreover, the detection of gravitational waves emitted during the binary neutron star merger, GW170817, and its electromagnetic counterpart, GRB170817A, impose additional constraints on the tidal deformability. In this work, we investigate hybrid stars with sequential phase transitions hadron-quark-quark in their cores. We assume that both phase transitions are sharp and analyse the rapid and slow phase conversion scenarios. For the outer core, we use modern hadronic equations of state. For the inner core we employ the constant speed of sound parametrization for quark matter. We analyze more than 3000 hybrid equations of state, taking into account the recent observational constraints from neutron stars. The effects of hadron-quark-quark phase transitions on the normal oscillation modes $f$ and $g$, are studied under the Cowling relativistic approximation. Our results show that, in the slow conversion regime, a second quark-quark phase transition gives rise to a new $g_2$~mode. We discuss the observational implications of our results associated to the gravitational waves detection and the possibility of detecting hints of sequential phase transitions and the associated $g_2$~mode.


title: Full compressible 3D magnetohydrodynamic simulation of solar wind

abstract: Identifying the heating mechanisms of the solar corona and the driving mechanisms of solar wind are key challenges in understanding solar physics. A full three-dimensional compressible magnetohydrodynamic (MHD) simulation was conducted to distinguish between the heating mechanisms in the fast solar wind above the open field region. Our simulation describes the evolution of the Alfv\'{e}nic waves, which includes the compressible effects from the photosphere to the heliospheric distance $s$ of 27 solar radii ($R_\odot$). The hot corona and fast solar wind were reproduced simultaneously due to the dissipation of the Alfv\'{e}n waves. The inclusion of the transition region and lower atmosphere enabled us to derive the solar mass loss rate for the first time by performing a full three-dimensional compressible MHD simulation. The Alfv\'{e}n turbulence was determined to be the dominant heating mechanism in the solar wind acceleration region ($s>1.3 R_\odot$), as suggested by previous solar wind models. In addition, shock formation and phase mixing are important below the lower transition region ($s<1.03R_\odot$) as well.


title: LACO: A Latency-Driven Network Slicing Orchestration in Beyond-5G Networks

abstract: Network Slicing is expected to become a game changer in the upcoming 5G networks and beyond, enlarging the telecom business ecosystem through still-unexplored vertical industry profits. This implies that heterogeneous service level agreements (SLAs) must be guaranteed per slice given the multitude of predefined requirements. In this paper, we pioneer a novel radio slicing orchestration solution that simultaneously provides-latency and throughput guarantees in a multi-tenancy environment. Leveraging on a solid mathematical framework, we exploit the exploration-vs-exploitation paradigm by means of a multi-armed-bandit-based(MAB) orchestrator, LACO, that makes adaptive resource slicing decisions with no prior knowledge on the traffic demand or channel quality statistics. As opposed to traditional MAB methods that are blind to the underlying system, LACO relies on system structure information to expedite decisions. After a preliminary simulations campaign empirically proving the validness of our solution, we provide a robust implementation of LACO using off-the-shelf equipment to fully emulate realistic network conditions:near-optimal results within affordable computational time are measured when LACO is in place.


title: Comparing anticyclotomic Selmer groups of positive coranks for congruent modular forms -- Part II

abstract: We study the Selmer group associated to a $p$-ordinary newform $f \in S_{2r}(\Gamma_0(N))$ over the anticyclotomic $\mathbb{Z}_p$-extension of an imaginary quadratic field $K/\mathbb{Q}$. Under certain assumptions, we prove that this Selmer group has no proper $\Lambda$-submodules of finite index. This generalizes work of Bertolini in the elliptic curve case. We also offer both a correction and an improvement to an earlier result on Iwasawa invariants of congruent modular forms by the present authors.


title: Nonexistence of variational minimizers related to a quasilinear singular problem in metric measure spaces

abstract: In this article we consider a variational problem related to a quasilinear singular problem and obtain a nonexistence result in a metric measure space with a doubling measure and a Poincar\'e inequality. Our method is purely variational and to the best of our knowledge, this is the first work concerning singular problems in a general metric setting.


title: Spin-induced scalarization of Kerr black holes with a massive scalar field

abstract: In the present paper we study the onset of the spin-induced scalarization of a Kerr black hole in scalar-Gauss-Bonnet gravity with a massive scalar field. Our approach is based on a (2+1) time evolution of the relevant linearized scalar field perturbation equation. We examine the region where the Kerr black hole becomes unstable giving rise to new scalarized rotating black holes with a massive scalar field. With increasing of the scalar field mass, the minimum value of the Gauss-Bonnet coupling parameter at which scalarization is possible, increases and thus the instability region shrinks. Interestingly, the introduction of scalar field mass does not change the critical minimal value of the black hole angular momentum $a_{\rm crit}/M$ where the instability of the Kerr black hole develops.


title: Accelerated Multi-Agent Optimization Method over Stochastic Networks

abstract: We propose a distributed method to solve a multi-agent optimization problem with strongly convex cost function and equality coupling constraints. The method is based on Nesterov's accelerated gradient approach and works over stochastically time-varying communication networks. We consider the standard assumptions of Nesterov's method and show that the sequence of the expected dual values converge toward the optimal value with the rate of $\mathcal{O}(1/k^2)$. Furthermore, we provide a simulation study of solving an optimal power flow problem with a well-known benchmark case.


title: Sequenced Route Query with Semantic Hierarchy

abstract: The trip planning query searches for preferred routes starting from a given point through multiple Point-of-Interests (PoI) that match user requirements. Although previous studies have investigated trip planning queries, they lack flexibility for finding routes because all of them output routes that strictly match user requirements. We study trip planning queries that output multiple routes in a flexible manner. We propose a new type of query called skyline sequenced route (SkySR) query, which searches for all preferred sequenced routes to users by extending the shortest route search with the semantic similarity of PoIs in the route. Flexibility is achieved by the {\it semantic hierarchy} of the PoI category. We propose an efficient algorithm for the SkySR query, bulk SkySR algorithm that simultaneously searches for sequenced routes and prunes unnecessary routes effectively. Experimental evaluations show that the proposed approach significantly outperforms the existing approaches in terms of response time (up to four orders of magnitude). Moreover, we develop a prototype service that uses the SkySR query, and conduct a user test to evaluate its usefulness.


title: Randomness Concerns When Deploying Differential Privacy

abstract: The U.S. Census Bureau is using differential privacy (DP) to protect confidential respondent data collected for the 2020 Decennial Census of Population & Housing. The Census Bureau's DP system is implemented in the Disclosure Avoidance System (DAS) and requires a source of random numbers. We estimate that the 2020 Census will require roughly 90TB of random bytes to protect the person and household tables. Although there are critical differences between cryptography and DP, they have similar requirements for randomness. We review the history of random number generation on deterministic computers, including von Neumann's "middle-square" method, Mersenne Twister (MT19937) (previously the default NumPy random number generator, which we conclude is unacceptable for use in production privacy-preserving systems), and the Linux /dev/urandom device. We also review hardware random number generator schemes, including the use of so-called "Lava Lamps" and the Intel Secure Key RDRAND instruction. We finally present our plan for generating random bits in the Amazon Web Services (AWS) environment using AES-CTR-DRBG seeded by mixing bits from /dev/urandom and the Intel Secure Key RDSEED instruction, a compromise of our desire to rely on a trusted hardware implementation, the unease of our external reviewers in trusting a hardware-only implementation, and the need to generate so many random bits.


title: Optical properties of a nonlinear magnetic charged rotating black hole surrounded by quintessence with a cosmological constant

abstract: In this paper,we discuss optical properties of the nonlinear magnetic charged black hole surrounded by quintessence with a non-zero cosmological constant $\Lambda$. Setting the state parameter $\omega=-3/2$ , we studied the horizon, the photon region and the shadow of this black hole. It turned out that for a fixed quintessential parameter $\gamma$, in a certain range, with the increase of the rotation parameter $a$ and magnetic charge $Q$, the inner horizon radius increases while the outer horizon radius decreases. And the cosmological horizon $r_{\Lambda}$ decrease when $\gamma$ or $\Lambda$ incease and increase slightly with increasing $a$ and $Q$. The shapes of photon region were then studied and depicted through graphical illustrations. Finally, we discussed the effects of the quintessential parameter $\gamma$ and the cosmological constant $\Lambda$ on the shadow cast by this balck hole with a fixed observer position.


title: Automatic Yara Rule Generation Using Biclustering

abstract: Yara rules are a ubiquitous tool among cybersecurity practitioners and analysts. Developing high-quality Yara rules to detect a malware family of interest can be labor- and time-intensive, even for expert users. Few tools exist and relatively little work has been done on how to automate the generation of Yara rules for specific families. In this paper, we leverage large n-grams ($n \geq 8$) combined with a new biclustering algorithm to construct simple Yara rules more effectively than currently available software. Our method, AutoYara, is fast, allowing for deployment on low-resource equipment for teams that deploy to remote networks. Our results demonstrate that AutoYara can help reduce analyst workload by producing rules with useful true-positive rates while maintaining low false-positive rates, sometimes matching or even outperforming human analysts. In addition, real-world testing by malware analysts indicates AutoYara could reduce analyst time spent constructing Yara rules by 44-86%, allowing them to spend their time on the more advanced malware that current tools can't handle. Code will be made available at https://github.com/NeuromorphicComputationResearchProgram .


title: Graphs and matrices: A translation of "Graphok \'{e}s matrixok" by D\'{e}nes K\H{o}nig (1931)

abstract: This paper, originally written in Hungarian by D\'{e}nes K\H{o}nig in 1931, proves that in a bipartite graph, the minimum vertex cover and the maximum matching have the same size. This statement is now known as K\H{o}nig's theorem. The paper also discusses the connection of graphs and matrices, then makes some observations about the combinatorial properties of the latter.


title: Products of locally cyclic groups

abstract: We consider groups of the form $G=AB$ with two locally cyclic subgroups $A$ and $B$. The structure of these groups is determined in the cases when $A$ and $B$ are both periodic or when one of them is periodic and the other is not. Together with a previous study of the case when $A$ and $B$ are torsion-free, this gives a complete classification of all groups that are the product of two locally cyclic groups. It is also shown that the product of a finite number of pairwise permutable periodic locally cyclic subgroups is a locally supersoluble group.


title: Analysis and Prediction of Deforming 3D Shapes using Oriented Bounding Boxes and LSTM Autoencoders

abstract: For sequences of complex 3D shapes in time we present a general approach to detect patterns for their analysis and to predict the deformation by making use of structural components of the complex shape. We incorporate long short-term memory (LSTM) layers into an autoencoder to create low dimensional representations that allow the detection of patterns in the data and additionally detect the temporal dynamics in the deformation behavior. This is achieved with two decoders, one for reconstruction and one for prediction of future time steps of the sequence. In a preprocessing step the components of the studied object are converted to oriented bounding boxes which capture the impact of plastic deformation and allow reducing the dimensionality of the data describing the structure. The architecture is tested on the results of 196 car crash simulations of a model with 133 different components, where material properties are varied. In the latent representation we can detect patterns in the plastic deformation for the different components. The predicted bounding boxes give an estimate of the final simulation result and their quality is improved in comparison to different baselines.


title: Payoff distribution in robust coalitional games on time-varying networks

abstract: In this paper, we consider a sequence of transferable utility (TU) coalitional games where the coalitional values are unknown but vary within certain bounds. As a solution to the resulting family of games, we formalise the notion of "robust CORE". Our main contribution is to design two distributed algorithms, namely, distributed payoff allocation and distributed bargaining, that converge to a consensual payoff distribution in the robust CORE. We adopt an operator-theoretic perspective to show convergence of both algorithms executed on time-varying communication networks. An energy storage optimization application motivates our framework for "robust coalitional games".


title: Improving Engagement of Animated Visualization with Visual Foreshadowing

abstract: Animated visualization is becoming increasingly popular as a compelling way to illustrate changes in time series data. However, maintaining the viewer's focus throughout the entire animation is difficult because of its time-consuming nature. Viewers are likely to become bored and distracted during the ever-changing animated visualization. Informed by the role of foreshadowing that builds the expectation in film and literature, we introduce visual foreshadowing to improve the engagement of animated visualizations. In specific, we propose designs of visual foreshadowing that engage the audience while watching the animation. To demonstrate our approach, we built a proof-of-concept animated visualization authoring tool that incorporates visual foreshadowing techniques with various styles. Our user study indicates the effectiveness of our foreshadowing techniques on improving engagement for animated visualization.


title: Casimir Energy for concentric $\delta$-$\delta'$ spheres

abstract: We study the vacuum interaction of a scalar field and two concentric spheres defined by a singular potential on their surfaces. The potential is a linear combination of the Dirac-$\delta$ and its derivative. The presence of the delta prime term in the potential causes that it behaves differently when it is seen from the inside or from the outside of the sphere. We study different cases for positive and negative values of the delta prime coupling, keeping positive the coupling of the delta. As a consequence, we find regions in the space of couplings, where the energy is positive, negative or zero. Moreover, the sign of the $\delta'$ couplings cause different behavior on the value of the Casimir energy for different values of the radii. This potential gives rise to general boundary conditions with limiting cases defining Dirichlet and Robin boundary conditions what allows us to simulate purely electric o purely magnetic spheres.


title: Comparison of the performance of high-order schemes based on the gas-kinetic and HLLC fluxes

abstract: In this paper, a comparison of the performance of two high-order finite volume methods based on the gas-kinetic scheme (GKS) and HLLC fluxes is carried out in structured rectangular mesh. For both schemes, the fifth-order WENO-AO reconstruction is adopted to achieve a high-order spatial accuracy. In terms of temporal discretization, a two-stage fourth-order (S2O4) time marching strategy is adopted for WENO5-AO-GKS scheme, and the fourth-order Runge-Kutta (RK4) method is employed for WENO5-AO-HLLC scheme. For the viscous flow computation, the GKS includes both inviscid and viscous fluxes in the evolution of a single cell interface gas distribution function. While for the WENO5-AO-HLLC scheme, the inviscid flux is provided by HLLC Riemann solver, and the viscous flux is discretized by a sixth-order central difference method. Based on the tests of forward Mach step and viscous shock tube, both schemes show outstanding shock capturing property. From the Titarev-Toro and double shear layer tests, WENO5-AO-GKS scheme seems to have a better resolution than WENO5-AO-HLLC scheme. Both schemes show excellent robustness in extreme cases, such as the Le Blanc problem. From the cases of the Noh problem and the compressible isotropic turbulence, WENO5-AO-GKS scheme shows favorite robustness. In the compressible isotropic turbulence and three-dimensional Taylor-Green vortex problems, WENO-AO-GKS can use a CFL number up to 0.5, instead of 0.3 for WENO5-AO-HLLC. In terms of computational efficiency, WENO5-AO-HLLC scheme is about 27% more expensive than WENO5-AO-GKS scheme in the two-dimensional viscous flow problems, but is about 15% faster in the three-dimensional case. Due to the multi-dimensionality, WENO5-AO-GKS scheme performs better than WENO5-AO-HLLC scheme in the laminar boundary layer and the double shear layer test.


title: Self-Supervised Scale Recovery for Monocular Depth and Egomotion Estimation

abstract: The self-supervised loss formulation for jointly training depth and egomotion neural networks with monocular images is well studied and has demonstrated state-of-the-art accuracy. One of the main limitations of this approach, however, is that the depth and egomotion estimates are only determined up to an unknown scale. In this paper, we present a novel scale recovery loss that enforces consistency between a known camera height and the estimated camera height, generating metric (scaled) depth and egomotion predictions. We show that our proposed method is competitive with other scale recovery techniques that require more information. Further, we demonstrate that our method facilitates network retraining within new environments, whereas other scale-resolving approaches are incapable of doing so. Notably, our egomotion network is able to produce more accurate estimates than a similar method which recovers scale at test time only.


title: Tight List-Sizes for Oblivious AVCs under Constraints

abstract: We study list-decoding over adversarial channels governed by oblivious adversaries (a.k.a. oblivious Arbitrarily Varying Channels (AVCs)). This type of adversaries aims to maliciously corrupt the communication without knowing the actual transmission from the sender. For any oblivious AVCs potentially with constraints on the sender's transmitted sequence and the adversary's noise sequence, we determine the exact value of the minimum list-size that can support a reliable communication at positive rate. This generalizes a classical result by Hughes (IEEE Transactions on Information Theory, 1997) and answers an open question posed by Sarwate and Gastpar (IEEE Transactions on Information Theory, 2012). A lower bound on the list-decoding capacity (whenever positive) is presented. Under a certain combinatorial conjecture, we also prove a matching upper bound. En route to a tight characterization of the list-decoding capacity, we propose a method for subcode construction towards the resolution of the combinatorial conjecture.


title: Thermally produced scalar singlet as leptonphilic dark matter

abstract: This work proposes a simple extension of the Standard Model where a scalar singlet and a leptonphilic scalar doublet exist. The scalar singlet could act as thermally produced dark matter. It annihilates into leptons through $s$-channel and benefit from the Breit--Winger resonance, so current thermally averaged cross-section are greatly enhanced. The constraints from the isotropic diffuse $\gamma$-ray background and correct dark matter relic abundance are analyzed. We find that the enhanced cross-section could fit the needed cross-section to explain the cosmic electron/positron excess.


title: On a connection used in deformation quantization

abstract: Using natural lifting operations, we give a coordinate-free proof of the fact that the connection used by Bordemann, Neumaier and Waldmann to construct the Fedosov standard ordered star product on the cotangent bundle of a Riemannian manifold is obtained by symplectification of the complete lift of the corresponding Levi-Civita connection, in the sense of Yano and Patterson. In terms of local coordinates, this has already been shown by Plebanski, Przanowski and Turrubiates.


title: Ergodic and non-ergodic dual-unitary quantum circuits with arbitrary local Hilbert space dimension

abstract: Dual-unitary quantum circuits can be used to construct 1+1 dimensional lattice models for which dynamical correlations of local observables can be explicitly calculated. We show how to analytically construct classes of dual-unitary circuits with any desired level of (non-)ergodicity for any dimension of the local Hilbert space, and present analytical results for thermalization to an infinite-temperature Gibbs state (ergodic) and a generalized Gibbs ensemble (non-ergodic). It is shown how a tunable ergodicity-inducing perturbation can be added to a non-ergodic circuit without breaking dual-unitarity, leading to the appearance of prethermalization plateaux for local observables.


title: ECG Beats Fast Classification Base on Sparse Dictionaries

abstract: Feature extraction plays an important role in Electrocardiogram (ECG) Beats classification system. Compared to other popular methods, VQ method performs well in feature extraction from ECG with advantages of dimensionality reduction. In VQ method, a set of dictionaries corresponding to segments of ECG beats is trained, and VQ codes are used to represent each heartbeat. However, in practice, VQ codes optimized by k-means or k-means++ exist large quantization errors, which results in VQ codes for two heartbeats of the same type being very different. So the essential differences between different types of heartbeats cannot be representative well. On the other hand, VQ uses too much data during codebook construction, which limits the speed of dictionary learning. In this paper, we propose a new method to improve the speed and accuracy of VQ method. To reduce the computation of codebook construction, a set of sparse dictionaries corresponding to wave segments of ECG beats is constructed. After initialized, sparse dictionaries are updated efficiently by Feature-sign and Lagrange dual algorithm. Based on those dictionaries, a set of codes can be computed to represent original ECG beats.Experimental results show that features extracted from ECG by our method are more efficient and separable. The accuracy of our method is higher than other methods with less time consumption of feature extraction


title: Linear Temporal Public Announcement Logic: a new perspective for reasoning about the knowledge of multi-classifiers

abstract: In this note, a formal transition system model called LTPAL to extract knowledge in a classification process is suggested. The model combines the Public Announcement Logic (PAL) and the Linear Temporal Logic (LTL). In the model, first, we consider classifiers, which capture single-framed data. Next, we took classifiers for data-stream data input into consideration. Finally, we formalize natural language properties in LTPAL with a video-stream object detection sample.


title: An isoperimetric problem with a competing nonlocal singular term

abstract: In this paper, we investigate the minimization of a functional in which the usual perimeter is competing with a nonlocal singular term comparable (but not necessarily equal to) a fractional perimeter. The motivation for this problem is a cell motility model introduced in some previous work by the first author. We establish several facts about global minimizers with a volume constraint. In particular we prove that minimizers exist and are radially symmetric for small mass, while minimizers cannot be radially symmetric for large mass. For large mass, we prove that the minimizing sequences either split into smaller sets that drift to infinity or develop fingers of a prescribed width. Finally, we connect these two alternatives to a related minimization problem for the optimal constant in a classical interpolation inequality (a Gagliardo-Nirenberg type inequality for fractional perimeter).


title: Spectro-temporal analysis of a sample of bursts from FRB 121102

abstract: FRB~121102 was the first Fast Radio Burst (FRB) that was shown to repeat. Since its discovery in 2012, more than two hundred bursts have been detected from the source. These bursts exhibit a diverse range of spectral and temporal characteristics and many questions about their origin and form remain unanswered. Here, we present a sample of radio bursts from FRB 121102 detected using the Lovell telescope at Jodrell Bank Observatory. We show four examples of bursts that show peculiar spectro-temporal characteristics and compare them with properties of bursts of FRB~121102 detected at other observatories. We report on a precursor burst that is separated by just 17~ms from the main burst, the shortest reported separation between two individual bursts to date. We also provide access to data for all the detections of FRB~121102 in this campaign.


title: DCTRGAN: Improving the Precision of Generative Models with Reweighting

abstract: Significant advances in deep learning have led to more widely used and precise neural network-based generative models such as Generative Adversarial Networks (GANs). We introduce a post-hoc correction to deep generative models to further improve their fidelity, based on the Deep neural networks using the Classification for Tuning and Reweighting (DCTR) protocol. The correction takes the form of a reweighting function that can be applied to generated examples when making predictions from the simulation. We illustrate this approach using GANs trained on standard multimodal probability densities as well as calorimeter simulations from high energy physics. We show that the weighted GAN examples significantly improve the accuracy of the generated samples without a large loss in statistical power. This approach could be applied to any generative model and is a promising refinement method for high energy physics applications and beyond.

categories: hep-ph hep-ex physics.data-an physics.ins-det stat.ML


title: Monotonicity of entropy for unimodal real quadratic rational maps

abstract: We show that the topological entropy is monotonic for unimodal interval maps which are obtained from the restriction of quadratic rational maps with real coefficients. This is done by ruling out the existence of certain post-critical curves in the moduli space of aforementioned maps, and confirms a conjecture made in [Fil19] based on experimental evidence.


title: The impact of COVID-19 on relative changes in aggregated mobility using mobile-phone data

abstract: Evaluating relative changes leads to additional insights which would remain hidden when only evaluating absolute changes. We analyze a dataset describing mobility of mobile phones in Austria before, during COVID-19 lock-down measures until recent. By applying compositional data analysis we show that formerly hidden information becomes available: we see that the elderly population groups increase relative mobility and that the younger groups especially on weekends also do not decrease their mobility as much as the others.


title: Nonrelativistic near-BPS corners of $\mathcal{N}=4$ super-Yang-Mills with $SU(1,1)$ symmetry

abstract: We consider limits of $\mathcal{N}=4$ super Yang-Mills (SYM) theory that approach BPS bounds and for which an $SU(1,1)$ structure is preserved. The resulting near-BPS theories become non-relativistic, with a $U(1)$ symmetry emerging in the limit that implies the conservation of particle number. They are obtained by reducing $\mathcal{N}=4$ SYM on a three-sphere and subsequently integrating out fields that become non-dynamical as the bounds are approached. Upon quantization, and taking into account normal-ordering, they are consistent with taking the appropriate limits of the dilatation operator directly, thereby corresponding to Spin Matrix theories, found previously in the literature. In the particular case of the $SU(1,1|1)$ near-BPS/Spin Matrix theory, we find a superfield formulation that applies to the full interacting theory. Moreover, for all the theories we find tantalizingly simple semi-local formulations as theories living on a circle. Finally, we find positive-definite expressions for the interactions in the classical limit for all the theories, which can be used to explore their strong coupling limits. This paper will have a companion paper in which we explore BPS bounds for which a $SU(2,1)$ structure is preserved.


title: AGN and Star-Formation Properties of Inside-out Assembled Galaxy Candidates at z<0.1

abstract: We study a sample of 48127 galaxies selected from the SDSS MPA-JHU catalogue, with $\log M_{\star}/M_{\odot} = 10.73 - 11.03$ and $z<0.1$. Local galaxies in this stellar mass range have been shown to have systematically shorter assembly times within their inner regions ($<0.5~R_{50}$) when compared to that of the galaxy as a whole, contrary to lower or higher mass galaxies which show consistent assembly times at all radii. Hence, we refer to these galaxies as Inside-Out Assembled Galaxy (IOAG) candidates. We find that the majority of IOAG candidates with well-detected emission lines are classified as either AGN (40%) or composite (40%) in the BPT diagram. We also find that the majority of our sources are located below the main sequence of star formation, and within the green valley or red sequence. Most BPT-classified star-forming IOAG candidates have spiral morphologies and are in the main sequence, whereas Seyfert 2 and composites have mostly spiral morphologies but quiescent star formation rates (SFRs). We argue that a high fraction of IOAG candidates seem to be in the process of quenching, moving from the blue cloud to the red sequence. Those classified as AGN have systematically lower SFRs than star-forming galaxies suggesting that AGN activity may be related to this quenching. However, the spiral morphology of these galaxies remains in place, suggesting that the central star-formation is suppressed before the morphological transformation occurs.


title: Modeling Curtailment in Germany: How Spatial Resolution Impacts Line Congestion

abstract: This paper investigates the effects of network constraints in energy system models at transmission level on renewable energy generation and curtailment as the network is being spatially aggregated. We seek to reproduce historically measured curtailment in Germany for the years 2013-2018 using an open model of the transmission system, PyPSA-Eur. Our simulations include spatial and temporal considerations, including congestion per line as well as curtailment per control zone and quarter. Results indicate that curtailment at high network resolution is significantly overestimated due to inaccurate allocation of electricity demand and renewable capacities to overloaded sites. However, high congestion rates of the transmission network decrease as the network is clustered to a smaller number of nodes, thus reducing curtailment. A measure to capture errors in the assignment of electricity demand and power plants is defined and hints towards a preferable spatial resolution. Thus, we are able to balance the effects of accurate node assignment and network congestion revealing that a reduced model can capture curtailment from recent historical data. This shows that it is possible to reduce the network to improve computation times and capture the most important effects of network constraints on variable renewable energy feed-in at the same time.


title: A sufficient condition for free-stream preserving in the nonlinear conservative finite difference schemes on curvilinear grids

abstract: In simulations of compressible flows, the conservative finite difference method (FDM) based on the nonlinear upwind schemes, e.g. WENO5, might violate free-stream preserving (FP), due to the loss of the geometric conservation law (GCL) identity when applied on the curvilinear grids. Although some techniques on FP have been proposed previously, no general rule is given for this issue. In this paper, by rearranging the upwind dissipation of the nonlinear schemes as a combination of sub-stencil reconstructions (taking WENO5 as an example), it can be proved that the upwind dissipation diminishes under the uniform flow condition if the metrics yield an identical value under the same schemes with these reconstructions, making the free-stream condition be preserved. According to this sufficient condition, the novel FP metrics are constructed for WENO5 and WENO7. By this means the original forms of these WENO schemes can be kept. In addition, the accuracy of these schemes can be retained as well with a simple accuracy compensation by replacing the central part fluxes with a high-order one. Various validations indicate that the present FP schemes retain the great capability to resolve the smooth regions accurately and capture the discontinuities robustly.


title: Modified estimator for the proportion of true null hypotheses under discrete setup with proven FDR control by the adaptive Benjamini-Hochberg procedure

abstract: Some crucial issues about a recently proposed estimator for the proportion of true null hypotheses ($\pi_0$) under discrete setup are discussed. An estimator for $\pi_0$ is introduced under the same setup. The estimator may be seen as a modification of a very popular estimator for $\pi_0$, originally proposed under the assumption of continuous test statistics. It is shown that adaptive Benjamini-Hochberg procedure remains conservative with the new estimator for $\pi_0$ being plugged in.


title: Nanomechanical damping via electron-assisted relaxation of two-level systems

abstract: We report on measurements of dissipation and frequency noise at millikelvin temperatures of nanomechanical devices covered with aluminum. A clear excess damping is observed after switching the metallic layer from superconducting to the normal state with a magnetic field. Beyond the standard model of internal tunneling systems coupled to the phonon bath, here we consider the relaxation to the conduction electrons together with the nature of the mechanical dispersion laws for stressed/unstressed devices. With these key ingredients, a model describing the relaxation of two-level systems inside the structure due to interactions with electrons and phonons with well separated timescales captures the data. In addition, we measure an excess 1/f-type frequency noise in the normal state, which further emphasizes the impact of conduction electrons.


title: Symmetry-Adapted Tight-Binding Electronic Structure Analysis of Carbon Nanotubes with Defects, Kinks, Twist, and Stretch

abstract: This paper applies a symmetry-adapted method to examine the influence of deformation and defects on the electronic structure and band structure in carbon nanotubes. First, the symmetry-adapted approach is used to develop the analog of Bloch waves. Building on this, the technique of perfectly-matched layers is applied to develop a method to truncate the computational domain of electronic structure calculations without spurious size effects. This provides an efficient and accurate numerical approach to compute the electronic structure and electromechanics of defects in nanotubes. The computational method is applied to study the effect of twist, stretch, and bending, with and without various types of defects, on the band structure of nanotubes. Specifically, the effect of stretch and twist on band structure in defect-free conducting and semiconducting nanotubes is examined, and the interaction with vacancy defects is elucidated. Next, the effect of localized bending or kinking on the electronic structure is studied. Finally, the paper examines the effect of 5-8-5 Stone-Wales defects. In all of these settings, the perfectly-matched layer method enables the calculation of localized non-propagating defect modes with energies in the bandgap of the defect-free nanotube.


title: A field theory in Randers-Finsler spacetime

abstract: Finsler geometry is a natural arena to investigate the physics of spacetimes with local Lorentz violating. The directional dependence of the Finsler metric provides a way to encode the Lorentz violating effects into the geometric structure of spacetime. Here, a classical field theory is proposed in a special Finsler geometry, the so-called Randers-Finsler spacetime, where the Lorentz violation is produced by a background vector field. By promoting the Randers-Finsler metric to a differential operator, a Finsler-invariant action for the scalar, gauge and fermions are proposed. The theory contains nonlocal terms, as in the Very Special Relativity based theories. By expanding the Lagrangian, minimal and nonminimal Standard Model Extension terms arises, revealing a perturbative Lorentz violation. For a CPT-even term, the Carrol-Field-Jackiw and derivative extensions are obtained.


title: Understanding Compositional Structures in Art Historical Images using Pose and Gaze Priors

abstract: Image compositions as a tool for analysis of artworks is of extreme significance for art historians. These compositions are useful in analyzing the interactions in an image to study artists and their artworks. Max Imdahl in his work called Ikonik, along with other prominent art historians of the 20th century, underlined the aesthetic and semantic importance of the structural composition of an image. Understanding underlying compositional structures within images is challenging and a time consuming task. Generating these structures automatically using computer vision techniques (1) can help art historians towards their sophisticated analysis by saving lot of time; providing an overview and access to huge image repositories and (2) also provide an important step towards an understanding of man made imagery by machines. In this work, we attempt to automate this process using the existing state of the art machine learning techniques, without involving any form of training. Our approach, inspired by Max Imdahl's pioneering work, focuses on two central themes of image composition: (a) detection of action regions and action lines of the artwork; and (b) pose-based segmentation of foreground and background. Currently, our approach works for artworks comprising of protagonists (persons) in an image. In order to validate our approach qualitatively and quantitatively, we conduct a user study involving experts and non-experts. The outcome of the study highly correlates with our approach and also demonstrates its domain-agnostic capability. We have open-sourced the code at https://github.com/image-compostion-canvas-group/image-compostion-canvas.


title: Preventing epidemics by wearing masks: An application to COVID-19

abstract: The goal of this work is to consider widespread use of face masks as a non-pharmaceutical control strategy for the Covid-19 pandemic. A SEIR model that divides the population into individuals that wear masks and those that do not is considered. After calculating the basic reproductive number by a next generation approach, a criterion for determining when an epidemic can be prevented by the use of masks only and the critical percentage of mask users for disease prevention in the population are derived. The results are then applied to real world data from the United States, Brazil and Italy.


title: Edge Degeneracy: Algorithmic and Structural Results

abstract: We consider a cops and robber game where the cops are blocking edges of a graph, while the robber occupies its vertices. At each round of the game, the cops choose some set of edges to block and right after the robber is obliged to move to another vertex traversing at most $s$ unblocked edges ($s$ can be seen as the speed of the robber). Both parts have complete knowledge of the opponent's moves and the cops win when they occupy all edges incident to the robbers position. We introduce the capture cost on $G$ against a robber of speed $s$. This defines a hierarchy of invariants, namely $\delta^{1}_{\rm e},\delta^{2}_{\rm e},\ldots,\delta^{\infty}_{\rm e}$, where $\delta^{\infty}_{\rm e}$ is an edge-analogue of the admissibility graph invariant, namely the {\em edge-admissibility} of a graph. We prove that the problem asking wether $\delta^{s}_{\rm e}(G)\leq k$, is polynomially solvable when $s\in \{1,2,\infty\}$ while, otherwise, it is NP-complete. Our main result is a structural theorem for graphs of bounded edge-admissibility. We prove that every graph of edge-admissibility at most $k$ can be constructed using $(\leq k)$-edge-sums, starting from graphs whose all vertices, except possibly from one, have degree at most $k$. Our structural result is approximately tight in the sense that graphs generated by this construction always have edge-admissibility at most $2k-1$. Our proofs are based on a precise structural characterization of the graphs that do not contain $\theta_{r}$ as an immersion, where $\theta_{r}$ is the graph on two vertices and $r$ parallel edges.


title: Recent advances in the 3D kinematic Babcock-Leighton solar dynamo modeling

abstract: In this review, we explain recent progress made in the Babcock-Leighton dynamo models for the Sun, which have been most successful to explain various properties of the solar cycle. In general, these models are 2D axisymmetric and the mean-field dynamo equations are solved in the meriodional plane of the Sun. Various physical processes (e.g., magnetic buoyancy and Babcock-Leighton mechanism) involved in these models are inherently 3D process and could not be modeled properly in a 2D framework. After pointing out limitations of 2D models (e.g., Mean-field Babcock-Leighton dynamo models and Surface Flux Transport models), we describe recently developed next-generation 3D dynamo models that implement more sophisticated flux emergence algorithm of buoyant flux tube rise through the convection zone and capture Babcock-Leighton process more realistically than previous 2D models. The detailed results from these 3D dynamo models including surface flux transport counterpart are presented. We explain the cycle irregularities that are reproduced in 3D dynamo models by introducing scattering around the tilt angle only. Some results by assimilating observed photospheric convective velocity fields into the 3D models are also discussed, pointing out the wide opportunity that these 3D models hold to deliver.


title: Frequency-domain calculation of Smith-Purcell radiation for metallic and dielectric gratings

abstract: The intensity of Smith-Purcell radiation from metallic and dielectric gratings (silicon, silica) is compared in a frequency-domain simulation. The numerical model is discussed and verified with the Frank-Tamm formula for Cherenkov radiation. For 30 keV electrons, rectangular dielectric gratings are less efficient than their metallic counterpart, by an order of magnitude for silicon, and two orders of magnitude for silica. For all gratings studied, radiation intensity oscillates with grating tooth height due to electromagnetic resonances in the grating. 3D and 2D numerical models are compared.


title: Route to logical strange nonchaotic attractors with single periodic force and noise

abstract: Strange nonchaotic attractors (SNAs) have been identified and studied in the literature exclusively in quasiperiodically driven nonlinear dynamical systems. It is an interesting question to ask whether they can be identified with other types of forcings as well, which still remains as an open problem. Here, we show that robust SNAs can be created by a small amount of noise in periodically driven nonlinear dynamical systems by a single force. The robustness of these attractors is tested by perturbing the system with logical signals leading to the emulation of different logical elements in the SNA regions.


title: Suppressing Frequency Fluctuations of Self-Sustained Vibrations in Underdamped Nonlinear Resonators

abstract: We consider frequency fluctuations in self-sustained oscillators based on nonlinear underdamped resonators. An important type of such resonators are nano- and micro-electro-mechanical systems. Various noise sources are considered, with the emphasis on the fundamentally unavoidable noise that comes along with dissipation from the coupling to a thermal reservoir. The formulation in terms of the action-angle variables of the resonator allows us to study a deeply nonlinear regime. In this regime the vibration frequency as a function of the action can have an extremum. We show that frequency fluctuations can be strongly reduced by choosing the operation point at this extremum. We suggest a practical implementation of a nanoresonator that has the appropriate property and show explicit results for the corresponding model.


title: Investigation of wakes generated by fractal plates in the compressible flow regime using large-eddy simulations

abstract: We investigate flows interacting with a square and a fractal shape multi-scale structures in the compressible regime for Mach numbers at subsonic and supersonic upstream conditions using large-eddy-simulations (LES). We also aim at identifying similarities and differences that these interactions have with corresponding interactions in the canonical incompressible flow problem. To account for the geometrical complexity associated with the fractal structures, we apply an immersed boundary method to model the no-slip boundary condition at the solid surfaces, with adequate mesh resolution in the vicinity of the small fractal features. We validate the numerical results through extensive comparisons with experimental wind tunnel measurements at a low Mach number. Similar to the incompressible flow case results, we find a break-up of the flow structures by the fractal plate and an increase in turbulent mixing in the downstream direction. As the Mach number increases, we observe noticeable wake meandering and higher spread rate of the wake in the lateral direction perpendicular to the streamwise-spanwise plane. Although not significant, we quantify the difference between the square and the fractal plates using two-point velocity correlations across the Mach number range. The wakes generated by the fractal plate in the compressible regime showed lower turbulent kinetic energy (TKE) and energy spectra levels compared to those of the square case. Moreover, results in terms of the near-field pressure spectra seem to indicate that the fractal plate has the potential to reduce the aerodynamic noise.


title: Further results on synergistic Lyapunov functions and hybrid feedback design through backstepping

abstract: We extend results on backstepping hybrid feedbacks by exploiting synergistic Lyapunov function and feedback (SLFF) pairs in a generalized form. Compared to existing results, we delineate SLFF pairs that are ready-made and do not require extra dynamic variables for backstepping. From an (weak) SLFF pair for an affine control system, we construct an SLFF pair for an extended system where the control input is produced through an integrator. The resulting hybrid feedback asymptotically stabilizes the extended system when the synergy gap for the original system is strictly positive. To highlight the versatility of SLFF pairs, we provide a result on the existence of a SLFF pair whenever a hybrid feedback stabilizer exists. The results are illustrated on the 3D pendulum.


title: PSO-PS: Parameter Synchronization with Particle Swarm Optimization for Distributed Training of Deep Neural Networks

abstract: Parameter updating is an important stage in parallelism-based distributed deep learning. Synchronous methods are widely used in distributed training the Deep Neural Networks (DNNs). To reduce the communication and synchronization overhead of synchronous methods, decreasing the synchronization frequency (e.g., every $n$ mini-batches) is a straightforward approach. However, it often suffers from poor convergence. In this paper, we propose a new algorithm of integrating Particle Swarm Optimization (PSO) into the distributed training process of DNNs to automatically compute new parameters. In the proposed algorithm, a computing work is encoded by a particle, the weights of DNNs and the training loss are modeled by the particle attributes. At each synchronization stage, the weights are updated by PSO from the sub weights gathered from all workers, instead of averaging the weights or the gradients. To verify the performance of the proposed algorithm, the experiments are performed on two commonly used image classification benchmarks: MNIST and CIFAR10, and compared with the peer competitors at multiple different synchronization configurations. The experimental results demonstrate the competitiveness of the proposed algorithm.


title: VisCode: Embedding Information in Visualization Images using Encoder-Decoder Network

abstract: We present an approach called VisCode for embedding information into visualization images. This technology can implicitly embed data information specified by the user into a visualization while ensuring that the encoded visualization image is not distorted. The VisCode framework is based on a deep neural network. We propose to use visualization images and QR codes data as training data and design a robust deep encoder-decoder network. The designed model considers the salient features of visualization images to reduce the explicit visual loss caused by encoding. To further support large-scale encoding and decoding, we consider the characteristics of information visualization and propose a saliency-based QR code layout algorithm. We present a variety of practical applications of VisCode in the context of information visualization and conduct a comprehensive evaluation of the perceptual quality of encoding, decoding success rate, anti-attack capability, time performance, etc. The evaluation results demonstrate the effectiveness of VisCode.


title: Search for intrinsic NALs in BAL/mini-BAL quasar spectra

abstract: Some fraction of narrow absorption lines are physically associated to the quasar/host-galaxy materials (i.e., intrinsic NALs) like those of BALs and mini-BALs. The relation between these three types of absorption lines has not been understood yet, however one interpretation is that these absorption features correspond to different inclination angles. In this study, we search for intrinsic NALs in 11 BAL/mini-BAL quasar spectra retrieved from VLT/UVES public archive, in order to test a possible relation of intrinsic NALs and BALs/mini-BALs in the geometry models. We use partial coverage analysis to separate intrinsic NALs from ones which are associated to cosmologically intervening materials like foreground galaxies and intergalactic medium (i.e., intervening NALs). We identify one reliable and two possible intrinsic NAL systems out of 36 NAL systems in 9 BAL/mini-BAL quasar spectra after removing two quasars without clear BAL features. In spite of a small sample size, we placed a lower limit on the fraction of BAL/mini-BAL quasars that have at least one intrinsic C IV NAL ($\sim33^{+33}_{-18}\%$). This can be interpreted that intrinsic NAL absorbers exist everywhere regardless of inclination angle. We found that one of the intrinsic NAL systems detected in SDSS J121549.80-003432.1 is located at a large radial distance of R > 130 kpc, using a method of photoionization model with ground/excited-state lines. Considering the wide range of intrinsic NAL absorber distribution in inclination angles and radial distances, it suggests that origins and geometry of them are more complicated than we expected.


title: CLF-Based Control for Hybrid Dynamical Systems

abstract: Pointwise minimum norm control laws for hybrid dynamical systems are proposed. Hybrid systems are given by differential equations capturing the continuous dynamics or flows, and by difference equations capturing the discrete dynamics or jumps. The proposed control laws are defined as the pointwise minimum norm selection from the set of inputs guaranteeing a decrease of a control Lyapunov function. The cases of individual and common inputs during flows and jumps, as well as when inputs enter through one of the system dynamics, are considered. Examples illustrate the results.


title: Understanding and Exploiting Dependent Variables with Deep Metric Learning

abstract: Deep Metric Learning (DML) approaches learn to represent inputs to a lower-dimensional latent space such that the distance between representations in this space corresponds with a predefined notion of similarity. This paper investigates how the mapping element of DML may be exploited in situations where the salient features in arbitrary classification problems vary over time or due to changing underlying variables. Examples of such variable features include seasonal and time-of-day variations in outdoor scenes in place recognition tasks for autonomous navigation and age/gender variations in human/animal subjects in classification tasks for medical/ethological studies. Through the use of visualisation tools for observing the distribution of DML representations per each query variable for which prior information is available, the influence of each variable on the classification task may be better understood. Based on these relationships, prior information on these salient background variables may be exploited at the inference stage of the DML approach by using a clustering algorithm to improve classification performance. This research proposes such a methodology establishing the saliency of query background variables and formulating clustering algorithms for better separating latent-space representations at run-time. The paper also discusses online management strategies to preserve the quality and diversity of data and the representation of each class in the gallery of embeddings in the DML approach. We also discuss latent works towards understanding the relevance of underlying/multiple variables with DML.


title: Cross-layer Band Selection and Routing Design for Diverse Band-aware DSA Networks

abstract: As several new spectrum bands are opening up for shared use, a new paradigm of \textit{Diverse Band-aware Dynamic Spectrum Access} (d-DSA) has emerged. d-DSA equips a secondary device with software defined radios (SDRs) and utilize whitespaces (or idle channels) in \textit{multiple bands}, including but not limited to TV, LTE, Citizen Broadband Radio Service (CBRS), unlicensed ISM. In this paper, we propose a decentralized, online multi-agent reinforcement learning based cross-layer BAnd selection and Routing Design (BARD) for such d-DSA networks. BARD not only harnesses whitespaces in multiple spectrum bands, but also accounts for unique electro-magnetic characteristics of those bands to maximize the desired quality of service (QoS) requirements of heterogeneous message packets; while also ensuring no harmful interference to the primary users in the utilized band. Our extensive experiments demonstrate that BARD outperforms the baseline dDSAaR algorithm in terms of message delivery ratio, however, at a relatively higher network latency, for varying number of primary and secondary users. Furthermore, BARD greatly outperforms its single-band DSA variants in terms of both the metrics in all considered scenarios.


title: A First Look at Zoombombing

abstract: Online meeting tools like Zoom and Google Meet have become central to our professional, educational, and personal lives. This has opened up new opportunities for large scale harassment. In particular, a phenomenon known as zoombombing has emerged, in which aggressors join online meetings with the goal of disrupting them and harassing their participants. In this paper, we conduct the first data-driven analysis of calls for zoombombing attacks on social media. We identify ten popular online meeting tools and extract posts containing meeting invitations to these platforms on a mainstream social network, Twitter, and on a fringe community known for organizing coordinated attacks against online users, 4chan. We then perform manual annotation to identify posts that are calling for zoombombing attacks, and apply thematic analysis to develop a codebook to better characterize the discussion surrounding calls for zoombombing. During the first seven months of 2020, we identify over 200 calls for zoombombing between Twitter and 4chan, and analyze these calls both quantitatively and qualitatively. Our findings indicate that the vast majority of calls for zoombombing are not made by attackers stumbling upon meeting invitations or bruteforcing their meeting ID, but rather by insiders who have legitimate access to these meetings, particularly students in high school and college classes. This has important security implications, because it makes common protections against zoombombing, such as password protection, ineffective. We also find instances of insiders instructing attackers to adopt the names of legitimate participants in the class to avoid detection, making countermeasures like setting up a waiting room and vetting participants less effective. Based on these observations, we argue that the only effective defense against zoombombing is creating unique join links for each participant.


title: QSAN: A Quantum-probability based Signed Attention Network for Explainable False Information Detection

abstract: False information detection on social media is challenging as it commonly requires tedious evidence-collecting but lacks available comparative information. Clues mined from user comments, as the wisdom of crowds, could be of considerable benefit to this task. However, it is non-trivial to capture the complex semantics from the contents and comments in consideration of their implicit correlations. Although deep neural networks have good expressive power, one major drawback is the lack of explainability. In this paper, we focus on how to learn from the post contents and related comments in social media to understand and detect the false information more effectively, with explainability. We thus propose a Quantum-probability based Signed Attention Network (QSAN) that integrates the quantum-driven text encoding and a novel signed attention mechanism in a unified framework. QSAN is not only able to distinguish important comments from the others, but also can exploit the conflicting social viewpoints in the comments to facilitate the detection. Moreover, QSAN is advantageous with its explainability in terms of transparency due to quantum physics meanings and the attention weights. Extensive experiments on real-world datasets show that our approach outperforms state-of-the-art baselines and can provide different kinds of user comments to explain why a piece of information is detected as false.


title: System Size Identification from Sinusoidal Probing in Diffusive Complex Networks

abstract: One of the most fundamental characteristic of a complex system is its size (or volume), which, in many modelling, is represented by the number of its individual components. Complex systems under investigation nowadays are typically large and/or time-varying, rendering their identification challenging. We propose here an accurate and efficient method to determine the size of (i.e., number of agents in) a complex, diffusively coupled dynamical system, that leverages the response of the system to an injected probing signal. For our method to be applicable, we rely on some assumptions on system's characteristics, namely, on the spectrum of the coupling graph and on the basin stability of its steady state(s). Even though such assumptions imply that our method cannot be applied to any instance of diffusively coupled group of dynamical agents, we argue that it covers relevant and interesting examples. Furthermore, the simplicity of the approach and its low computational complexity renders it very interesting for the systems to which it applies.


title: Optimal training of integer-valued neural networks with mixed integer programming

abstract: Recent work has shown potential in using Mixed Integer Programming (MIP) solvers to optimize certain aspects of neural networks (NNs). However the intriguing approach of training NNs with MIP solvers is under-explored. State-of-the-art-methods to train NNs are typically gradient-based and require significant data, computation on GPUs, and extensive hyper-parameter tuning. In contrast, training with MIP solvers does not require GPUs or heavy hyper-parameter tuning, but currently cannot handle anything but small amounts of data. This article builds on recent advances that train binarized NNs using MIP solvers. We go beyond current work by formulating new MIP models which improve training efficiency and which can train the important class of integer-valued neural networks (INNs). We provide two novel methods to further the potential significance of using MIP to train NNs. The first method optimizes the number of neurons in the NN while training. This reduces the need for deciding on network architecture before training. The second method addresses the amount of training data which MIP can feasibly handle: we provide a batch training method that dramatically increases the amount of data that MIP solvers can use to train. We thus provide a promising step towards using much more data than before when training NNs using MIP models. Experimental results on two real-world data-limited datasets demonstrate that our approach strongly outperforms the previous state of the art in training NN with MIP, in terms of accuracy, training time and amount of data. Our methodology is proficient at training NNs when minimal training data is available, and at training with minimal memory requirements -- which is potentially valuable for deploying to low-memory devices.


title: Quantum phenomenological gravitational dynamics: A general view from thermodynamics of spacetime

abstract: In this work we derive general quantum phenomenological equations of gravitational dynamics and analyse its features. The derivation uses the formalism developed in thermodynamics of spacetime and introduces low energy quantum gravity modifications to it. Quantum gravity effects are considered via modification of Bekenstein entropy by an extra logarithmic term in the area. This modification is predicted by several approaches to quantum gravity, including loop quantum gravity, string theory, AdS/CFT correspondence and generalised uncertainty principle phenomenology, giving our result a general character. The derived equations generalise classical equations of motion of unimodular gravity, instead of the ones of general relativity, and they contain at most second derivatives of the metric. We provide two independent derivations of the equations based on thermodynamics of local causal diamonds. First one uses Jacobson's maximal vacuum entanglement hypothesis, the second one Clausius entropy flux. Furthermore, we consider questions of diffeomorphism and local Lorentz invariance of the resulting dynamics and discuss its application to a simple cosmological model, finding a resolution of the classical singularity.


title: Maximal singular integral operators acting on noncommutative $L_p$-spaces

abstract: In this paper, we study the boundedness theory for maximal Calder\'on-Zygmund operators acting on noncommutative $L_p$-spaces. Our first result is a criterion for the weak type $(1,1)$ estimate of noncommutative maximal Calder\'on-Zygmund operators; as an application, we obtain the weak type $(1,1)$ estimates of operator-valued maximal singular integrals of convolution type under proper {regularity} conditions. These are the {\it first} noncommutative maximal inequalities for families of linear operators that can not be reduced to positive ones. For homogeneous singular integrals, the strong type $(p,p)$ ($1<p<\infty$) maximal estimates are shown to be true even for {rough} kernels. As a byproduct of the criterion, we obtain the noncommutative weak type $(1,1)$ estimate for Calder\'on-Zygmund operators with integral regularity condition that is slightly stronger than the H\"ormander condition; this evidences somewhat an affirmative answer to an open question in the noncommutative Calder\'on-Zygmund theory.


title: Completely Dark Photons from Gravitational Particle Production During Inflation

abstract: Starting with the de Broglie--Proca Lagrangian for a massive vector field, we calculate the number density of particles resulting from gravitational particle production (GPP) during inflation, with detailed consideration to the evolution of the number density through the reheating. We find plausible scenarios for the production of dark-photon dark matter of mass in a wide range, as low as a micro-electron volt to $10^{14}$ GeV. Gravitational particle production does not depend on any coupling of the dark photon to standard-model particles.


title: A Robotic Positive Psychology Coach to Improve College Students' Wellbeing

abstract: A significant number of college students suffer from mental health issues that impact their physical, social, and occupational outcomes. Various scalable technologies have been proposed in order to mitigate the negative impact of mental health disorders. However, the evaluation for these technologies, if done at all, often reports mixed results on improving users' mental health. We need to better understand the factors that align a user's attributes and needs with technology-based interventions for positive outcomes. In psychotherapy theory, therapeutic alliance and rapport between a therapist and a client is regarded as the basis for therapeutic success. In prior works, social robots have shown the potential to build rapport and a working alliance with users in various settings. In this work, we explore the use of a social robot coach to deliver positive psychology interventions to college students living in on-campus dormitories. We recruited 35 college students to participate in our study and deployed a social robot coach in their room. The robot delivered daily positive psychology sessions among other useful skills like delivering the weather forecast, scheduling reminders, etc. We found a statistically significant improvement in participants' psychological wellbeing, mood, and readiness to change behavior for improved wellbeing after they completed the study. Furthermore, students' personality traits were found to have a significant association with intervention efficacy. Analysis of the post-study interview revealed students' appreciation of the robot's companionship and their concerns for privacy.


title: Hybrid Thermal Machines: Generalized Thermodynamic Resources for Multitasking

abstract: Thermal machines perform useful tasks--such as producing work, cooling, or heating--by exchanging energy, and possibly additional conserved quantities such as particles, with reservoirs. Here we consider thermal machines that perform more than one useful task simultaneously, terming these "hybrid thermal machines". We outline their restrictions imposed by the laws of thermodynamics and we quantify their performance in terms of efficiencies. To illustrate their full potential, reservoirs that feature multiple conserved quantities, described by generalized Gibbs ensembles, are considered. A minimal model for a hybrid thermal machine is introduced, featuring three reservoirs and two conserved quantities, e.g., energy and particle number. This model can be readily implemented in a thermoelectric setup based on quantum dots, and hybrid regimes are accessible considering realistic parameters.


title: Refined approachability algorithms and application to regret minimization with global costs

abstract: Blackwell's approachability is a framework where two players, the Decision Maker and the Environment, play a repeated game with vector-valued payoffs. The goal of the Decision Maker is to make the average payoff converge to a given set called the target. When this is indeed possible, simple algorithms which guarantee the convergence are known. This abstract tool was successfully used for the construction of optimal strategies in various repeated games, but also found several applications in online learning. By extending an approach proposed by (Abernethy et al., 2011), we construct and analyze a class of Follow the Regularized Leader algorithms (FTRL) for Blackwell's approachability which are able to minimize not only the Euclidean distance to the target set (as it is often the case in the context of Blackwell's approachability) but a wide range of distance-like quantities. This flexibility enables us to apply these algorithms to closely minimize the quantity of interest in various online learning problems. In particular, for regret minimization with $\ell_p$ global costs, we obtain the first bounds with explicit dependence in $p$ and the dimension $d$.


title: Simplifying the design of multilevel thermal machines using virtual qubits

abstract: Quantum thermodynamics often deals with the dynamics of small quantum machines interfacing with a large and complex environment. Virtual qubits, collisional models and reset master equations have become highly useful tools for predicting the qualitative behaviour of two-dimensional target systems coupled to few-qubit machines and a thermal environment. While few successes in matching the simplified model parameters for all possible physical systems are known, the qualitative predictions still allow for a general design of quantum machines irrespective of the implementation. We generalise these tools by introducing multiple competing virtual qubits for modelling multi-dimensional systems coupled to larger and more complex machines. By simulating the full physical dynamics for targets with three dimensions, we uncover general properties of reset models that can be used as `dials' to correctly predict the qualitative features of physical changes in a realistic setup and thus design autonomous quantum machines beyond a few qubits. We then present a general analytic solution of the reset model for arbitrary-dimensional systems coupled to multi-qubit machines. Finally, we showcase an improved three-level laser as an exemplary application of our results.


title: Evidence for a soft-phonon-mode-driven Peierls-type distortion in Sc$_3$CoC$_4$

abstract: We provide experimental and theoretical evidence for the realization of the Peierls-type structurally distorted state in the quasi-one-dimensional superconductor Sc$_3$CoC$_4$ by a phonon-softening mechanism. The transition from the high- to the final low-temperature phase below 80K proceeds via an extended intermediate temperature regime between 80K and 150K characterized by phonon-driven atom displacements. In support of the low-dimensional character of the title compound we find a highly anisotropic correlation-length of these dynamic distortions.


title: Influence of disorder on a Bragg microcavity

abstract: Using the resonant-state expansion for leaky optical modes of a planar Bragg microcavity, we investigate the influence of disorder on its fundamental cavity mode. We model the disorder by randomly varying the thickness of the Bragg-pair slabs (composing the mirrors) and the cavity, and calculate the resonant energy and linewidth of each disordered microcavity exactly, comparing the results with the resonant-state expansion for a large basis set and within its first and second orders of perturbation theory. We show that random shifts of interfaces cause a growth of the inhomogeneous broadening of the fundamental mode that is proportional to the magnitude of disorder. Simultaneously, the quality factor of the microcavity decreases inversely proportional to the square of the magnitude of disorder. We also find that first-order perturbation theory works very accurately up to a reasonably large disorder magnitude, especially for calculating the resonance energy, which allows us to derive qualitatively the scaling of the microcavity properties with disorder strength.


title: Gauss maps of surfaces in 3-dimensional Heisenberg group

abstract: In this paper, we study the Gauss map of surfaces in 3-dimensional Heisenberg group using the Gans model of the hyperbolic plane. We establish a relationship between the tension field of the Gauss map and mean curvature of a surface in $\mathcal{H}_{3}$.


title: Graph neural networks-based Scheduler for Production planning problems using Reinforcement Learning

abstract: Reinforcement learning (RL) is increasingly adopted in job shop scheduling problems (JSSP). But RL for JSSP is usually done using a vectorized representation of machine features as the state space. It has three major problems: (1) the relationship between the machine units and the job sequence is not fully captured, (2) exponential increase in the size of the state space with increasing machines/jobs, and (3) the generalization of the agent to unseen scenarios. We present a novel framework - GraSP-RL, GRAph neural network-based Scheduler for Production planning problems using Reinforcement Learning. It represents JSSP as a graph and trains the RL agent using features extracted using a graph neural network (GNN). While the graph is itself in the non-euclidean space, the features extracted using the GNNs provide a rich encoding of the current production state in the euclidean space, which is then used by the RL agent to select the next job. Further, we cast the scheduling problem as a decentralized optimization problem in which the learning agent is assigned to all the production units and the agent learns asynchronously from the data collected on all the production units. The GraSP-RL is then applied to a complex injection molding production environment with 30 jobs and 4 machines. The task is to minimize the makespan of the production plan. The schedule planned by GraSP-RL is then compared and analyzed with a priority dispatch rule algorithm like first-in-first-out (FIFO) and metaheuristics like tabu search (TS) and genetic algorithm (GA). The proposed GraSP-RL outperforms the FIFO, TS, and GA for the trained task of planning 30 jobs in JSSP. We further test the generalization capability of the trained agent on two different problem classes: Open shop system (OSS) and Reactive JSSP (RJSSP) where our method produces results better than FIFO and comparable results to TS and GA.


title: Derived parabolic induction

abstract: The classical parabolic induction functor is a fundamental tool on the representation theoretic side of the Langlands program. In this article, we study its derived version. It was shown by the second author that the derived category of smooth $G$-representations over $k$, $G$ a $p$-adic reductive group and $k$ a field of characteristic $p$, is equivalent to the derived category of a certain differential graded $k$-algebra $H_G^\bullet$, whose zeroth cohomology is a classical Hecke algebra. This equivalence predicts the existence of a derived parabolic induction functor on the dg Hecke algebra side, which we construct in this paper. This relies on the theory of six-functor formalisms for differential graded categories developed by O.\ Schn\"urer. We also discuss the adjoint functors of derived parabolic induction.


title: New opportunities at the photon energy frontier

abstract: Ultra-peripheral collisions (UPCs) involving heavy ions and protons are the energy frontier for photon-mediated interactions. UPC photons can be used for many purposes, including probing low-$x$ gluons via photoproduction of dijets and vector mesons, probes of beyond-standard-model processes, such as those enabled by light-by-light scattering, and studies of two-photon production of the Higgs.


title: Microwave Superconductivity

abstract: We give a broad overview of the history of microwave superconductivity and explore the technological developments that have followed from the unique electrodynamic properties of superconductors. Their low loss properties enable resonators with high quality factors that can nevertheless handle extremely high current densities. This in turn enables superconducting particle accelerators, high-performance filters and analog electronics, including metamaterials, with extreme performance. The macroscopic quantum properties have enabled new generations of ultra-high-speed digital computing and extraordinarily sensitive detectors. The microscopic quantum properties have enabled large-scale quantum computers, which at their heart are essentially microwave-fueled quantum engines. We celebrate the rich history of microwave superconductivity and look to the promising future of this exciting branch of microwave technology.


title: Robust and Efficient Optimization Using a Marquardt-Levenberg Algorithm with R Package marqLevAlg

abstract: Implementations in R of classical general-purpose algorithms for local optimization generally have two major limitations which cause difficulties in applications to complex problems: too loose convergence criteria and too long calculation time. By relying on a Marquardt-Levenberg algorithm (MLA), a Newton-like method particularly robust for solving local optimization problems, we provide with marqLevAlg package an efficient and general-purpose local optimizer which (i) prevents convergence to saddle points by using a stringent convergence criterion based on the relative distance to minimum/maximum in addition to the stability of the parameters and of the objective function; and (ii) reduces the computation time in complex settings by allowing parallel calculations at each iteration. We demonstrate through a variety of cases from the literature that our implementation reliably and consistently reaches the optimum (even when other optimizers fail), and also largely reduces computational time in complex settings through the example of maximum likelihood estimation of different sophisticated statistical models.


title: Genealogies in bistable waves

abstract: We study a model of selection acting on a diploid population (one in which each individual carries two copies of each gene) living in one spatial dimension. We suppose a particular gene appears in two forms (alleles) $A$ and $a$, and that individuals carrying $AA$ have a higher fitness than $aa$ individuals, while $Aa$ individuals have a lower fitness than both $AA$ and $aa$ individuals. The proportion of advantageous $A$ alleles expands through the population approximately according to a travelling wave. We prove that on a suitable timescale, the genealogy of a sample of $A$ alleles taken from near the wavefront converges to a Kingman coalescent as the population density goes to infinity. This contrasts with the case of directional selection in which the corresponding limit is thought to be the Bolthausen-Sznitman coalescent. The proof uses 'tracer dynamics'.


title: Classical Yang-Mills observables from amplitudes

abstract: The double copy suggests that the basis of the dynamics of general relativity is Yang-Mills theory. Motivated by the importance of the relativistic two-body problem, we study the classical dynamics of colour-charged particle scattering from the perspective of amplitudes, rather than equations of motion. We explain how to compute the change of colour, and the radiation of colour, during a classical collision. We apply our formalism at next-to-leading order for the colour change and at leading order for colour radiation.


title: Structure and water attachment rates of ice in the atmosphere: role of nitrogen

abstract: In this work we perform computer simulations of the ice surface in order to elucidate the role of nitrogen in the crystal growth rates and crystal habits of snow in the atmosphere. In pure water vapor at temperatures typical of ice crystal formation in cirrus clouds, we find that basal and primary prismatic facets exhibit a layer of premelted ice, with thickness in the subnanometer range. For partial pressures of 1 bar, well above the expected values in the troposphere, we find that only small amounts of nitrogen are adsorbed. The adsorption takes place onto the premelted surface, and hardly any nitrogen dissolves within the premelting film. The premelting film thickness does not change either. We quantify the resulting change of the ice/vapor surface tension to be in the hundredth of mN/m and find that the structure of the pristine ice surface is not changed in a significant manner. We perform a trajectory analysis of colliding water molecules, and find that the attachment rates from direct ballistic collision are very close to unity irrespective of the nitrogen pressure. Nitrogen is however at sufficient density to deflect a fraction of trajectories with smaller distance than the mean free path. Our results show explicitly that the reported differences in growth rates measured in pure water vapor and a controlled nitrogen atmosphere are not related to a significant disruption of the ice surface due to nitrogen adsorption. On the contrary, we show clearly from our trajectory analysis that nitrogen slows down the crystal growth rates due to collisions between water molecules with bulk nitrogen gas. This clarifies the long standing controversy of the role of inert gases on crystal growth rates and demonstrates their influence is solely related to the diffusion limited flow of water vapor across the gas phase.

categories: physics.chem-ph cond-mat.mes-hall physics.ao-ph physics.comp-ph physics.geo-ph