title: A New Family of Algebraically Defined Graphs With Small Automorphism Group

abstract: Let $p$ be an odd prime, $q=p^e$, $e\ge 1$, and $\mathbb{F} = \mathbb{F_q}$ denote the finite field of $q$ elements. Let $f: \mathbb{F}^2\to \mathbb{F}$ and $g: \mathbb{F}^3\to \mathbb{F}$ be functions, and let $P$ and $L$ be two copies of the 3-dimensional vector space $\mathbb{F}^3$. Consider a bipartite graph $\Gamma _\mathbb{F} (f, g)$ with vertex partitions $P$ and $L$ and with edges defined as follows: for every $(p)=(p_1,p_2,p_3)\in P$ and every $[l]= [l_1,l_2,l_3]\in L$, $\{(p), [l]\} = (p)[l]$ is an edge in $\Gamma _\mathbb{F} (f, g)$ if $$p_2+l_2 =f(p_1,l_1) \;\;\;\text{and}\;\;\; p_3 + l_3 = g(p_1,p_2,l_1).$$ Given $\Gamma _\mathbb{F} (f, g)$, is it always possible to find a function $h:\mathbb{F}^2\to \mathbb{F}$ such that the graph $\Gamma _\mathbb{F} (f, h)$ with the same vertex set as $\Gamma _\mathbb{F} (f, g)$ and with edges $(p)[l]$ defined in a similar way by the system $$p_2+l_2 =f(p_1,l_1) \;\;\;\text{and}\;\;\; p_3 + l_3 = h(p_1,l_1),$$ is isomorphic to $\Gamma _\mathbb{F} (f, g)$ for infinitely many $q$? In this paper we show that the answer to the question is negative and the graphs $\Gamma_{\mathbb{F}_p}(p_1\ell_1, p_1\ell_1p_2(p_1 + p_2 + p_1p_2))$ provide such an example for $p \equiv 1 \pmod{3}$. Our argument is based on proving that the automorphism group of these graphs has order $p$, which is the smallest possible order of the automorphism group of graphs of the form $\Gamma_{\mathbb{F}}(f, g)$.


title: Near-infrared spectroscopic indices for unresolved stellar populations. II. Index measurements

abstract: We measured the equivalent width of a large set of near-infrared (NIR, 0.8--2.4$ \ \mu$m) line-strength indices in the XShooter medium-resolution spectra of the central regions of 14 galaxies. We found that two aluminum indices Al at 1.31 $\mu$m and Al1 at 1.67 $\mu$m and the two CO indices CO1 at 1.56 $\mu$m and CO4 at 1.64 $\mu$m are tightly correlated with the velocity dispersion. Moreover, the NIR Al and CO1 indices show strong correlations with the optical Mg2 and Mgb indices, which are usually adopted as $\alpha$/Fe-enhancement diagnostics. The molecular FeH1 index at 1.58 $\mu$m tightly correlates with the optical <Fe> and [MgFe]' indices, which are used as total metallicity diagnostics. The NIR Pa$\beta$ index at 1.28 $\mu$m has a behaviour similar to the optical H$\beta$ index, which is a diagnostic of mean age. We defined two new composite indices, <Al> and [AlFeH], as possible candidates to be used as NIR diagnostics of total metallicity and $\alpha$/Fe enhancement. The NIR <Al> index has a strong correlation with the optical Mg2 and Mgb indices, while the [AlFeH] index is tightly correlated with the optical <Fe> and [MgFe]' indices. The distribution of the data points in the NIR Pa$\beta$-<Al> and Pa$\beta$-[AlFeH] diagrams mimic that in the optical [MgFe]'-H$\beta$ and the Mgb-<Fe> diagrams, which are widely used to constraint the properties of the unresolved stellar populations. We concluded that some NIR line-strength indices could be useful in studying stellar populations as well as in fine-tuning stellar population models.


title: Robust Estimation of Effective Diffusions from Multiscale Data

abstract: We present a novel methodology based on filtered data and moving averages for estimating effective dynamics from observations of multiscale systems. We show in a semi-parametric framework of the Langevin type that our approach is asymptotically unbiased with respect to the theory of homogenization. Moreover, we demonstrate on a range of challenging numerical experiments that our method is accurate in extracting coarse-grained dynamics from multiscale data. In particular, the estimators we propose are more robust and require less knowledge of the full model than the standard technique of subsampling, which is widely employed in practice in this setting.


title: Time Evolution of ML-MCTDH Wavefunctions I: Gauge Conditions, Basis Functions, and Singularities

abstract: We derive a family of equations of motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than that used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOMs approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.


title: Time Evolution of ML-MCTDH Wavefunctions II: Application of the Projector Splitting Integrator

abstract: The multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) approach suffers from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of a regularization parameter. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to $10^6$ bath modes, and find that for these problems the PSI requires roughly 3-4 orders of magnitude fewer Hamiltonian evaluations and 2-3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH, and 2-3/1-2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to $\sim 1.3\times 10^9$ parameters to obtain accurate dynamics.


title: Tensor Monopoles in superconducting systems

abstract: Topology in general but also topological objects such as monopoles are a central concept in physics. They are prime examples for the intriguing physics of gauge theories and topological states of matter. Vector monopoles are already frequently discussed such as the well-established Dirac monopole in three dimensions. Less known are tensor monopoles giving rise to tensor gauge fields. Here we report that tensor monopoles can potentially be realized in superconducting multi-terminal systems using the phase differences between superconductors as synthetic dimensions. In a first proposal we suggest a circuit of superconducting islands featuring charge states to realize a tensor monopole. As a second example we propose a triple dot system coupled to multiple superconductors that also gives rise to such a topological structure. All proposals can be implemented with current experimental means and the monopole readily be detected by measuring the quantum geometry.


title: The Quantum Dynamics of Two-component Bose-Einstein Condensate: an $Sp(4,R)$ Symmetry Approach

abstract: The compact groups such as $SU(n)$ and $SO(n)$ groups have been heavily studied and applied in the study of quantum many body systems. However, the non-compact groups such as the real symplectic groups are less touched. In this paper, we will reveal that the quantum dynamics of two-component Bose-Einstein condensate can be described by a \emph{non-compact} real symplectic group $Sp(4,R)$. With this group, we can give a explicit form for the wavefunction in any time of the evolution, meanwhile, map this whole time evolution to a trajectory in a six-dimensional manifold. By introducing a polar coordinate, we can visualize this six-dimensional manifold in 2d unit disk and reveal the relation between the behavior of the trajectory in this manifold and the eigen-energies of the Hamiltonian. Furthermore, the time evolution of expectation value of a physical observable such as number operator is proven closely related to the behavior of the trajectory in this manifold.


title: NumGPT: Improving Numeracy Ability of Generative Pre-trained Models

abstract: Existing generative pre-trained language models (e.g., GPT) focus on modeling the language structure and semantics of general texts. However, those models do not consider the numerical properties of numbers and cannot perform robustly on numerical reasoning tasks (e.g., math word problems and measurement estimation). In this paper, we propose NumGPT, a generative pre-trained model that explicitly models the numerical properties of numbers in texts. Specifically, it leverages a prototype-based numeral embedding to encode the mantissa of the number and an individual embedding to encode the exponent of the number. A numeral-aware loss function is designed to integrate numerals into the pre-training objective of NumGPT. We conduct extensive experiments on four different datasets to evaluate the numeracy ability of NumGPT. The experiment results show that NumGPT outperforms baseline models (e.g., GPT and GPT with DICE) on a range of numerical reasoning tasks such as measurement estimation, number comparison, math word problems, and magnitude classification. Ablation studies are also conducted to evaluate the impact of pre-training and model hyperparameters on the performance.


title: Interests, Difficulties, Sentiments, and Tool Usages of Concurrency Developers: A Large-Scale Study on Stack Overflow

abstract: Context: Software developers are increasingly facing the challenges of writing code that is not only concurrent but also correct. Objective: To help these developers, it is necessary to understand concurrency topics they are interested in, their difficulty in finding answers for questions in these topics, their sentiment for these topics, and how they use concurrency tools and techniques to guarantee correctness. Method: We conduct a large-scale study on the entirety of Stack Overflow to understand interests, difficulties, sentiment, and tool usages of concurrency developers. We discuss the implications of our findings for the practice, research, and education of concurrent software development, and investigate the relation of our findings with the findings of the previous work. Results: A few findings of our study are: (1) questions that concurrency developers ask can be grouped into a hierarchy with 27 concurrency topics under 8 major categories, (2) thread safety is among the most popular concurrency topics and client-server concurrency is among the least popular, (3) irreproducible behavior is among the most difficult topics and memory consistency is among the least difficult, (4) data scraping is among the most positive concurrency topics and irreproducible behavior is among the most negative, (5) root cause identification has the most number of questions for usage of data race tools and alternative use has the least. Conclusion: The results of our study can not only help concurrency developers but also concurrency educators and researchers to better decide where to focus their efforts, by trading off one concurrency topic against another.


title: An Executable Structural Operational Formal Semantics for Python

abstract: Python is a popular high-level general-purpose programming language also heavily used by the scientific community. It supports a variety of different programming paradigms and is preferred by many for its ease of use. With the vision of harvesting static analysis techniques like abstract interpretation for Python, we develop a formal semantics for Python. A formal semantics is an important cornerstone for any sound static analysis technique. We base our efforts on the general framework of structural operational semantics yielding a small-step semantics in principle allowing for concurrency and interaction with an environment. The main contributions of this thesis are twofold: first, we develop a meta-theoretic framework for the formalization of structural operational semantics in tandem with the necessary tool support for the automated derivation of interpreters from such formal semantics, and, second, we validate the suitability of this approach for the formalization of modern programming languages developing a semantics for Python.


title: Measuring the Cosmic X-ray Background accurately

abstract: Measuring the Cosmic X-ray Background (CXB) is a key to understand the Active Galactic Nuclei population, their absorption distribution and their average spectra. However, hard X-ray instruments suffer from time-dependent backgrounds and cross-calibration issues. The uncertainty of the CXB normalization remain of the order of 20%. To obtain a more accurate measurement, the Monitor Vsego Neba (MVN) instrument was built in Russia but not yet launched to the ISS (arXiv:1410.3284). We follow the same ideas to develop a CXB detector made of four collimated spectrometers with a rotating obturator on top. The collimators block off-axis photons below 100 keV and the obturator modulates on-axis photons allowing to separate the CXB from the instrumental background. Our spectrometers are made of 20 mm thick CeBr$_{3}$ crystals on top of a SiPM array. One tube features a $\sim$20 cm$^2$ effective area and more energy coverage than MVN, leading to a CXB count rate improved by a factor of $\sim$10 and a statistical uncertainty $\sim$0.5% on the CXB flux. A prototype is being built and we are seeking for a launch opportunity.


title: Smart Traffic Monitoring System using Computer Vision and Edge Computing

abstract: Traffic management systems capture tremendous video data and leverage advances in video processing to detect and monitor traffic incidents. The collected data are traditionally forwarded to the traffic management center (TMC) for in-depth analysis and may thus exacerbate the network paths to the TMC. To alleviate such bottlenecks, we propose to utilize edge computing by equipping edge nodes that are close to cameras with computing resources (e.g. cloudlets). A cloudlet, with limited computing resources as compared to TMC, provides limited video processing capabilities. In this paper, we focus on two common traffic monitoring tasks, congestion detection, and speed detection, and propose a two-tier edge computing based model that takes into account of both the limited computing capability in cloudlets and the unstable network condition to the TMC. Our solution utilizes two algorithms for each task, one implemented at the edge and the other one at the TMC, which are designed with the consideration of different computing resources. While the TMC provides strong computation power, the video quality it receives depends on the underlying network conditions. On the other hand, the edge processes very high-quality video but with limited computing resources. Our model captures this trade-off. We evaluate the performance of the proposed two-tier model as well as the traffic monitoring algorithms via test-bed experiments under different weather as well as network conditions and show that our proposed hybrid edge-cloud solution outperforms both the cloud-only and edge-only solutions.


title: Polarization measurements of the Crab Pulsar with POLAR

abstract: POLAR is a dedicated Gamma-Ray Burst polarimeter making use of Compton-scattering which took data from the second Chinese spacelab, the Tiangong-2 from September 2016 to April 2017. It has a wide Field of View of $\sim6$ steradians and an effective area of $\sim400\ cm^2$ at 300 keV. These features make it one of the most sensitive instruments in its energy range (15-500 keV), and therefore capable of almost continuously monitoring persistent sources such as pulsars. Significant folded pulsation from both PSR B0531+21 (the Crab Pulsar) and PSR B1509-58 has been observed. Observations of the Crab Pulsar with POLAR have previously been used for phase-resolved spectroscopy of the Crab Pulsar to calibrate the instrumental responses of POLAR. In this work, we investigate a polarimetric joint-fitting method for observations of the Crab Pulsar with POLAR. Unlike a GRB observation with POLAR, the observations of the Crab Pulsar are complicated by multiple observational datasets during which the polarization plane rotates as well. So before fitting, we have to correct the modulation curves under different datasets, by taking into account the rotations of the Crab Pulsar's relative position in the detctor's local coordinate, and the changes of detector response in different datasets. Despite these difficulties and the low signal to background for such sources constraining, polarization measurements were possible with the POLAR data. We will present the methodology briefly, which could be applied to any wide FoV polarimeter, and polarization results of the Crab pulsar with POLAR. Finally, the inferred ability of pulsar detection with POLAR-2 (the successor of POLAR) will also be discussed.


title: Dynamics of a trapped ion in a quantum gas: effects of particle statistics

abstract: We study the quantum dynamics of an ion confined in a radiofrequency trap in interaction with either a Bose or spin-polarized Fermi gas. To this end, we derive quantum optical master equations in the limit of weak coupling and the Lamb-Dicke approximations. For the bosonic bath, we also include the so-called "Lamb-shift" correction to the ion trap due to the coupling to the quantum gas as well as the extended Fr\"ohlich interaction within the Bogolyubov approximation that have been not considered in previous studies. We calculate the ion kinetic energy for various atom-ion scattering lengths as well as gas temperatures by considering the intrinsic micromotion and we analyse the damping of the ion motion in the gas as a function of the gas temperature. We find that the ion's dynamics depends on the quantum statistics of the gas and that a fermionic bath enables to attain lower ionic energies.


title: PP-OCRv2: Bag of Tricks for Ultra Lightweight OCR System

abstract: Optical Character Recognition (OCR) systems have been widely used in various of application scenarios. Designing an OCR system is still a challenging task. In previous work, we proposed a practical ultra lightweight OCR system (PP-OCR) to balance the accuracy against the efficiency. In order to improve the accuracy of PP-OCR and keep high efficiency, in this paper, we propose a more robust OCR system, i.e. PP-OCRv2. We introduce bag of tricks to train a better text detector and a better text recognizer, which include Collaborative Mutual Learning (CML), CopyPaste, Lightweight CPUNetwork (LCNet), Unified-Deep Mutual Learning (U-DML) and Enhanced CTCLoss. Experiments on real data show that the precision of PP-OCRv2 is 7% higher than PP-OCR under the same inference cost. It is also comparable to the server models of the PP-OCR which uses ResNet series as backbones. All of the above mentioned models are open-sourced and the code is available in the GitHub repository PaddleOCR which is powered by PaddlePaddle.


title: Accurate sky signal reconstruction for ground-based spectroscopy with kinetic inductance detectors

abstract: Context. Wide-field spectrometers are needed to deal with current astrophysical challenges that require multiband observations at millimeter wavelengths. An example of these is the KIDs Interferometer Spectrum Survey (KISS), which uses two arrays of kinetic inductance detectors (KIDs) coupled to a Martin-Puplett interferometer (MPI). KISS has a wide instantaneous field of view (1 deg in diameter) and a spectral resolution up to 1.45 GHz in the 120-180 GHz electromagnetic band. The instrument is installed on the 2.25 m Q-U-I JOint TEnerife telescope in Teide Observatory (Tenerife, Canary Islands), at an altitude of 2395 m above sea level. Aims. This work presents an original readout modulation method developed to improve the sky signal reconstruction accuracy for types of instruments for which a fast sampling frequency is required both to remove atmospheric fluctuations and to perform full spectroscopic measurements on each sampled sky position. Methods. We first demonstrate the feasibility of this technique using simulations. Then, we apply such a scheme to on-sky calibration. Results. We show that the sky signal can be reconstructed to better than 0.5% for astrophysical sources, and to better than 2% for large background variations such as in "skydip", in an ideal noiseless scenario. The readout modulation method is validated by observations on-sky during the KISS commissioning campaign. Conclusions. We conclude that accurate photometry can be obtained for future KID-based MPI.


title: Toward Generating Sufficiently Valid Test Case Results: A Method for Systematically Assigning Test Cases to Test Bench Configurations in a Scenario-Based Test Approach for Automated Vehicles

abstract: To successfully launch automated vehicles into the consumer market, there must be credible proof that the vehicles will operate safely. However, finding a method to validate the vehicles' safe operation is a challenging problem. While scenario-based test approaches seem to be possible solutions, they require execution of a large number of test cases. Several test benches, ranging from actual test vehicles to partly or fully simulated environments, are available to execute these test cases. Each test bench provides different elements, which in turn, have different parameters and parameter ranges. The composition of elements with their specific parameter values at a specific test bench that is used to execute a test case is referred to as a test bench configuration. However, selecting the most suitable test bench configuration is difficult. The selected test bench configuration determines whether the execution of a specific test case provides sufficiently valid test case results with respect to the intended purpose, for example, validating a vehicle's safe operation. The effective and efficient execution of a large number of test cases requires a method for systematically assigning test cases to the most suitable test bench configuration. Based on a proposed method for classifying test bench configurations, we propose and illustrate a method for systematically assigning test cases to test bench configurations in a scenario-based test approach for automated vehicles. This assignment method allows for the effective and efficient execution of a large number of test cases while generating sufficiently valid test case results.


title: Interferometric sensing of a commercial geophone

abstract: We present a modified commercial L-4C geophone with interferometric readout that demonstrated a resolution 60 times lower than the included coil-magnet readout at low frequencies. The intended application for the modified sensor is in vibration isolation platforms that require improved performance at frequencies lower than 1 Hz. A controls and noise-model of an Advanced LIGO 'HAM-ISI' vibration isolation system was developed, and it shows that our sensor can reduce the residual vibration by a factor of 70 at 0.1 Hz


title: Congruency-Constrained TU Problems Beyond the Bimodular Case

abstract: A long-standing open question in Integer Programming is whether integer programs with constraint matrices with bounded subdeterminants are efficiently solvable. An important special case thereof are congruency-constrained integer programs $\min\{c^\top x\colon\ Tx\leq b,\ \gamma^\top x\equiv r\pmod{m},\ x\in\mathbb{Z}^n\}$ with a totally unimodular constraint matrix $T$. Such problems have been shown to be polynomial-time solvable for $m=2$, which led to an efficient algorithm for integer programs with bimodular constraint matrices, i.e., full-rank matrices whose $n\times n$ subdeterminants are bounded by two in absolute value. Whereas these advances heavily relied on existing results on well-known combinatorial problems with parity constraints, new approaches are needed beyond the bimodular case, i.e., for $m>2$. We make first progress in this direction through several new techniques. In particular, we show how to efficiently decide feasibility of congruency-constrained integer programs with a totally unimodular constraint matrix for $m=3$. Furthermore, for general $m$, our techniques also allow for identifying flat directions of infeasible problems, and deducing bounds on the proximity between solutions of the problem and its relaxation.


title: LAMOST Time-Domain Survey: First Results of four $K$2 plates

abstract: From Oct. 2019 to Apr. 2020, LAMOST performs a time-domain spectroscopic survey of four $K$2 plates with both low- and med-resolution observations. The low-resolution spectroscopic survey gains 282 exposures ($\approx$46.6 hours) over 25 nights, yielding a total of about 767,000 spectra, and the med-resolution survey takes 177 exposures ($\approx$49.1 hours) over 27 nights, collecting about 478,000 spectra. More than 70%/50% of low-resolution/med-resolution spectra have signal-to-noise ratio higher than 10. We determine stellar parameters (e.g., $T_{\rm eff}$, log$g$, [Fe/H]) and radial velocity (RV) with different methods, including LASP, DD-Payne, and SLAM. In general, these parameter estimations from different methods show good agreement, and the stellar parameter values are consistent with those of APOGEE. We use the $Gaia$ DR2 RV data to calculate a median RV zero point (RVZP) for each spectrograph exposure by exposure, and the RVZP-corrected RVs agree well with the APOGEE data. The stellar evolutionary and spectroscopic masses are estimated based on the stellar parameters, multi-band magnitudes, distances and extinction values. Finally, we construct a binary catalog including about 2700 candidates by analyzing their light curves, fitting the RV data, calculating the binarity parameters from med-resolution spectra, and cross-matching the spatially resolved binary catalog from $Gaia$ EDR3. The LAMOST TD survey is expected to get breakthrough in various scientific topics, such as binary system, stellar activity, and stellar pulsation, etc.


title: Recommendation Fairness: From Static to Dynamic

abstract: Driven by the need to capture users' evolving interests and optimize their long-term experiences, more and more recommender systems have started to model recommendation as a Markov decision process and employ reinforcement learning to address the problem. Shouldn't research on the fairness of recommender systems follow the same trend from static evaluation and one-shot intervention to dynamic monitoring and non-stop control? In this paper, we portray the recent developments in recommender systems first and then discuss how fairness could be baked into the reinforcement learning techniques for recommendation. Moreover, we argue that in order to make further progress in recommendation fairness, we may want to consider multi-agent (game-theoretic) optimization, multi-objective (Pareto) optimization, and simulation-based optimization, in the general framework of stochastic games.


title: Spectrally Extended Line Field Optical Coherence Tomography Angiography

abstract: Optical coherence tomography angiography (OCTA) has been established as a powerful tool for investigating vascular diseases and is expected to become a standard of care technology. However, its widespread clinical usage is hindered by technical gaps such as limited field of view (FOV), lack of quantitative flow information, and suboptimal motion correction. Here we report a new imaging platform, termed spectrally extended line field (SELF) OCTA that provides advanced solutions to the above-mentioned challenges. SELF-OCTA breaks the speed limitations and achieves two-fold gain in FOV without sacrificing signal strength through parallel image acquisition. Towards quantitative angiography, the 'frequency flow' imaging mechanism overcomes the imaging speed bottleneck by obviating the requirement for superfluous B-scans. In addition, the 'frequency flow' imaging mechanism facilitates OCTA-data based motion tracking with overlap between adjacent line fields. Since it can be implemented in any existing OCT device without significant hardware modification or affecting existing functions, we expect that SELF-OCTA will make non-invasive, wide field, quantitative, and low-cost angiographic imaging available to larger patient populations.


title: Acceleration of the order of convergence of a family of fractional fixed point methods and its implementation in the solution of a nonlinear algebraic system related to hybrid solar receivers

abstract: This paper presents a way to define, classify and accelerate the order of convergence of an uncountable family of fractional fixed point methods, which may be useful to continue expanding the applications of fractional operators. The proposed method to accelerate convergence is used in a fractional iterative method, and with the obtained method are solved simultaneously two nonlinear algebraic systems that depend on time-dependent parameters, and that allow obtaining the temperatures and efficiencies of a hybrid solar receiver. Finally, two uncountable families of fractional fixed point methods are presented, in which the proposed method to accelerate convergence can be implemented.

categories: math.NA cs.NA math.FA math.OA physics.app-ph physics.comp-ph


title: An eXtended Finite Element Method Implementation in COMSOL Multiphysics: Solid Mechanics

abstract: This paper presents the first time implementation of the eXtended Finite Element Method (XFEM) in the general purpose commercial software COMSOL Multiphysics. An enrichment strategy is proposed, consistent with the structure of the software. To this end, for each set of enrichment functions, an additional Solid Mechanics module is incorporated into the numerical framework, coupled with compatible modifications to the internal variables. The Linear Elastic Fracture Mechanics (LEFM) is exclusively adopted for the crack analysis. The model pre-processing, level set update, stress intensity factor calculation and crack propagation analysis are conducted by employing COMSOL's built-in features in conjunction with external MATLAB functions through COMSOL LiveLink. All implementational aspects and suggested remedies for the treatment of enriched elements, framework setup, evaluation of stress intensity factors, and numerical integration are described in detail. The accuracy and robustness of the proposed method are examined by several numerical examples for stationary and propagating crack problems in 2D and 3D settings. The results represent excellent agreement with available analytical, numerical and experimental observations in the literature. Keywords:XFEM; COMSOL Multiphysics; Crack analysis; Fracture propagation


title: PEEK: A Large Dataset of Learner Engagement with Educational Videos

abstract: Educational recommenders have received much less attention in comparison to e-commerce and entertainment-related recommenders, even though efficient intelligent tutors have great potential to improve learning gains. One of the main challenges in advancing this research direction is the scarcity of large, publicly available datasets. In this work, we release a large, novel dataset of learners engaging with educational videos in-the-wild. The dataset, named Personalised Educational Engagement with Knowledge Topics PEEK, is the first publicly available dataset of this nature. The video lectures have been associated with Wikipedia concepts related to the material of the lecture, thus providing a humanly intuitive taxonomy. We believe that granular learner engagement signals in unison with rich content representations will pave the way to building powerful personalization algorithms that will revolutionise educational and informational recommendation systems. Towards this goal, we 1) construct a novel dataset from a popular video lecture repository, 2) identify a set of benchmark algorithms to model engagement, and 3) run extensive experimentation on the PEEK dataset to demonstrate its value. Our experiments with the dataset show promise in building powerful informational recommender systems. The dataset and the support code is available publicly.


title: PAUSE: Positive and Annealed Unlabeled Sentence Embedding

abstract: Sentence embedding refers to a set of effective and versatile techniques for converting raw text into numerical vector representations that can be used in a wide range of natural language processing (NLP) applications. The majority of these techniques are either supervised or unsupervised. Compared to the unsupervised methods, the supervised ones make less assumptions about optimization objectives and usually achieve better results. However, the training requires a large amount of labeled sentence pairs, which is not available in many industrial scenarios. To that end, we propose a generic and end-to-end approach -- PAUSE (Positive and Annealed Unlabeled Sentence Embedding), capable of learning high-quality sentence embeddings from a partially labeled dataset. We experimentally show that PAUSE achieves, and sometimes surpasses, state-of-the-art results using only a small fraction of labeled sentence pairs on various benchmark tasks. When applied to a real industrial use case where labeled samples are scarce, PAUSE encourages us to extend our dataset without the liability of extensive manual annotation work.


title: Fair Comparison: Quantifying Variance in Resultsfor Fine-grained Visual Categorization

abstract: For the task of image classification, researchers work arduously to develop the next state-of-the-art (SOTA) model, each bench-marking their own performance against that of their predecessors and of their peers. Unfortunately, the metric used most frequently to describe a model's performance, average categorization accuracy, is often used in isolation. As the number of classes increases, such as in fine-grained visual categorization (FGVC), the amount of information conveyed by average accuracy alone dwindles. While its most glaring weakness is its failure to describe the model's performance on a class-by-class basis, average accuracy also fails to describe how performance may vary from one trained model of the same architecture, on the same dataset, to another (both averaged across all categories and at the per-class level). We first demonstrate the magnitude of these variations across models and across class distributions based on attributes of the data, comparing results on different visual domains and different per-class image distributions, including long-tailed distributions and few-shot subsets. We then analyze the impact various FGVC methods have on overall and per-class variance. From this analysis, we both highlight the importance of reporting and comparing methods based on information beyond overall accuracy, as well as point out techniques that mitigate variance in FGVC results.


title: Tight binding approximation for bilayer graphene and nanotube structures: from commensurability to incommensurability between the layers

abstract: One- and two-dimensional bilayer systems are examples of ultra-tunable quantum materials that are considered as the basis for the new generation of electronic and photonic devices. Here we develop a general theory of the electron band structure for such commensurate and incommensurate bilayer carbon structures within the tight binding approximation. To model the band structure of commensurate twisted bilayer graphene (TBLG), we apply the classic zone folding theory. The latter leads us to the construction of TBLG Hamiltonians in the basis of shifted Bloch wave functions (SBWF), which, in contrast to the usual Bloch functions, have the wave vector q shifted by a set of vectors Q_i. The dimension of the considered Hamiltonians is equal to 4T, where the factor T is a number of vertices Q_i of the folded reciprocal space falling into the original first Brillouin zone of any of the layers. We propose and discuss a method for choosing a reduced set of SBWFs to construct effective Hamiltonians that correctly describe the low-energy spectrum of commensurate TBLG. The flattening of low-energy bands with a decrease in twist angle is discussed. As we show, this spectrum results from interactions between the lowest-energy modes of the folded dispersion curves. The effective Hamiltonians for calculating the low-energy band structure of incommensurate TBLG and double-walled carbon nanotubes (DWCNTs) are constructed in a similar way. To test the developed theory, we calculate the energies of 105 intra-tube optical transitions in 29 DWCNTs and compare them with previously published and novel experimental data. We also apply the theory to calculate the energies of recently discovered inter-tube transitions. Geometrical conditions allowing this type of transitions are discussed. We show that these transitions occur in DWCNTs which layers have close chiral angles and the same handedness.


title: Idiosyncratic but not Arbitrary: Learning Idiolects in Online Registers Reveals Distinctive yet Consistent Individual Styles

abstract: An individual's variation in writing style is often a function of both social and personal attributes. While structured social variation has been extensively studied, e.g., gender based variation, far less is known about how to characterize individual styles due to their idiosyncratic nature. We introduce a new approach to studying idiolects through a massive cross-author comparison to identify and encode stylistic features. The neural model achieves strong performance at authorship identification on short texts and through an analogy-based probing task, showing that the learned representations exhibit surprising regularities that encode qualitative and quantitative shifts of idiolectal styles. Through text perturbation, we quantify the relative contributions of different linguistic elements to idiolectal variation. Furthermore, we provide a description of idiolects through measuring inter- and intra-author variation, showing that variation in idiolects is often distinctive yet consistent.


title: Analysis of Regularized Learning for Linear-functional Data in Banach Spaces

abstract: In this article, we study the whole theory of regularized learning for linear-functional data in Banach spaces including representer theorems, pseudo-approximation theorems, and convergence theorems. The input training data are composed of linear functionals in the predual space of the Banach space to represent the discrete local information of multimodel data and multiscale models. The training data and the multi-loss functions are used to compute the empirical risks to approximate the expected risks, and the regularized learning is to minimize the regularized empirical risks over the Banach spaces. The exact solutions of the original problems are approximated globally by the regularized learning even if the original problems are unknown or unformulated. In the convergence theorems, we show the convergence of the approximate solutions to the exact solutions by the weak* topology of the Banach space. Moreover, the theorems of the regularized learning are applied to solve many problems of machine learning such as support vector machines and artificial neural networks.


title: How much pretraining data do language models need to learn syntax?

abstract: Transformers-based pretrained language models achieve outstanding results in many well-known NLU benchmarks. However, while pretraining methods are very convenient, they are expensive in terms of time and resources. This calls for a study of the impact of pretraining data size on the knowledge of the models. We explore this impact on the syntactic capabilities of RoBERTa, using models trained on incremental sizes of raw text data. First, we use syntactic structural probes to determine whether models pretrained on more data encode a higher amount of syntactic information. Second, we perform a targeted syntactic evaluation to analyze the impact of pretraining data size on the syntactic generalization performance of the models. Third, we compare the performance of the different models on three downstream applications: part-of-speech tagging, dependency parsing and paraphrase identification. We complement our study with an analysis of the cost-benefit trade-off of training such models. Our experiments show that while models pretrained on more data encode more syntactic knowledge and perform better on downstream applications, they do not always offer a better performance across the different syntactic phenomena and come at a higher financial and environmental cost.


title: Stable Positron Acceleration in Thin, Warm, Hollow Plasma Channels

abstract: Hollow plasma channels are attractive for lepton acceleration because they provide intrinsic emittance preservation regimes. However, beam breakup instabilities dominate the dynamics. Here, we show that thin, warm hollow channels can sustain large-amplitude plasma waves ready for high-quality positron acceleration. We verify that the combination of warm electrons and thin hollow channel enables positron focusing structures. Such focusing wakefields unlock beam breakup damping mechanisms. We demonstrate that such channels emerge self-consistently during the long-term plasma dynamics in the blowout's regime aftermath, allowing for experimental demonstration.


title: The pyglaf argumentation reasoner (ICCMA2021)

abstract: The pyglaf reasoner takes advantage of circumscription to solve computational problems of abstract argumentation frameworks. In fact, many of these problems are reduced to circumscription by means of linear encodings, and a few others are solved by means of a sequence of calls to an oracle for circumscription. Within pyglaf, Python is used to build the encodings and to control the execution of the external circumscription solver, which extends the SAT solver glucose and implements algorithms taking advantage of unsatisfiable core analysis and incremental computation.


title: Concentration of the complexity of spherical pure $p$-spin models at arbitrary energies

abstract: We consider critical points of the spherical pure $p$-spin spin glass model with Hamiltonian $H_{N}\left(\boldsymbol{\sigma}\right)=\frac{1}{N^{\left(p-1\right)/2}}\sum_{i_{1},...,i_{p}=1}^{N}J_{i_{1},...,i_{p}}\sigma_{i_{1}}\cdots\sigma_{i_{p}}$, where $\boldsymbol{\sigma}=\left(\sigma_{1},...,\sigma_{N}\right)\in \mathbb{S}^{N-1}:=\left\{ \boldsymbol{\sigma}\in\mathbb{R}^{N}:\,\left\Vert \boldsymbol{\sigma}\right\Vert _{2}=\sqrt{N}\right\} $ and $J_{i_{1},...,i_{p}}$ are i.i.d standard normal variables. Using a second moment analysis, we prove that for $p\geq 32$ and any $E>-E_\infty$, where $E_\infty$ is the (normalized) ground state, the ratio of the number of critical points $\boldsymbol{\sigma}$ with $H_N(\boldsymbol{\sigma})\leq NE$ and its expectation asymptotically concentrates at $1$. This extends to arbitrary $E$ a similar conclusion of [Sub17a].


title: Instabilities of complex fluids with partially structured and partially random interactions

abstract: We develop a theory for thermodynamic instabilities of complex fluids composed of many interacting chemical species organised in families. This model includes partially structured and partially random interactions and can be solved exactly using tools from random matrix theory. The model exhibits three kinds of fluid instabilities: one in which the species form a condensate with a local density that depends on their family (family condensation); one in which species demix in two phases depending on their family (family demixing); and one in which species demix in a random manner irrespective of their family (random demixing). We determine the critical spinodal density of the three types of instabilities and find that the critical spinodal density is finite for both family condensation and family demixing, while for random demixing the critical spinodal density grows as the square root of the number of species. We use the developed framework to describe phase-separation instability of the cytoplasm induced by a change in pH.

categories: physics.bio-ph cond-mat.dis-nn cond-mat.soft cond-mat.stat-mech


title: A note on infinite number of exact Lagrangian fillings for spherical spuns

abstract: In this short note we discuss high-dimensional examples of Legendrian submanifolds of the standard contact Euclidean space with an infinite number of exact Lagrangian fillings up to Hamiltonian isotopy. They are obtained from the examples of Casals and Ng by applying to them the spherical spinning construction.


title: Aspartix-V21

abstract: In this solver description we present ASPARTIX-V, in its 2021 edition, which participates in the International Competition on Computational Models of Argumentation (ICCMA) 2021. ASPARTIX-V is capable of solving all classical (static) reasoning tasks part of ICCMA'21 and extends the ASPARTIX system suite by incorporation of recent ASP language constructs (e.g. conditional literals), domain heuristics within ASP, and multi-shot methods. In this light ASPARTIX-V deviates from the traditional focus of ASPARTIX on monolithic approaches (i.e., one-shot solving via a single ASP encoding) to further enhance performance.


title: Complete stationary spacelike surfaces in an $n$-dimensional Generalized Robertson-Walker spacetime

abstract: Several uniqueness results for non-compact complete stationary spacelike surfaces in an $n(\geq 3)$-dimensional Generalized Robertson Walker spacetime are obtained. In order to do that, we assume a natural inequality involving the Gauss curvature of the surface, the restrictions of the warping function and the sectional curvature of the fiber to the surface. This inequality gives the parabolicity of the surface. Using this property, a distinguished non-negative superharmonic function on the surface is shown to be constant, which implies that the stationary spacelike surface must be totally geodesic. Moreover, non-trivial examples of stationary spacelike surfaces in the four dimensional Lorentz-Minkowski spacetime are exposed to show that each of our assumptions is needed.


title: Locally Recoverable Streaming Codes for Packet-Erasure Recovery

abstract: Streaming codes are a class of packet-level erasure codes that are designed with the goal of ensuring recovery in low-latency fashion, of erased packets over a communication network. It is well-known in the streaming code literature, that diagonally embedding codewords of a $[\tau+1,\tau+1-a]$ Maximum Distance Separable (MDS) code within the packet stream, leads to rate-optimal streaming codes capable of recovering from $a$ arbitrary packet erasures, under a strict decoding delay constraint $\tau$. Thus MDS codes are geared towards the efficient handling of the worst-case scenario corresponding to the occurrence of $a$ erasures. In the present paper, we have an increased focus on the efficient handling of the most-frequent erasure patterns. We study streaming codes which in addition to recovering from $a>1$ arbitrary packet erasures under a decoding delay $\tau$, have the ability to handle the more common occurrence of a single-packet erasure, while incurring smaller delay $r<\tau$. We term these codes as $(a,\tau,r)$ locally recoverable streaming codes (LRSCs), since our single-erasure recovery requirement is similar to the requirement of locality in a coded distributed storage system. We characterize the maximum possible rate of an LRSC by presenting rate-optimal constructions for all possible parameters $\{a,\tau,r\}$. Although the rate-optimal LRSC construction provided in this paper requires large field size, the construction is explicit. It is also shown that our $(a,\tau=a(r+1)-1,r)$ LRSC construction provides the additional guarantee of recovery from the erasure of $h, 1 \leq h \leq a$, packets, with delay $h(r+1)-1$. The construction thus offers graceful degradation in decoding delay with increasing number of erasures.


title: Aspects of topological superconductivity in 2D systems: noncollinear magnetism, skyrmions, and higher-order topology

abstract: The review is aimed at highlighting the aspects of topological superconductivity in the absence of spin-orbit interaction in two-dimensional systems with long-range non-collinear spin ordering or magnetic skyrmions. Another purpose is to give a brief introduction to the new concept of topological superconductivity, i.e. higher-order topology in two-dimensional systems including spin-orbit coupled structures. The formation of Majorana modes due to magnetic textures is discussed. The role of effective triplet pairings and odd fermion parity of the ground state wave function in different systems is emphasized. We describe the peculiarities of the magnetic skyrmions, leading to the formation of the Majorana modes and defects on which the modes are localized. The problem of braiding in the two-dimensional systems, especially in higher-order topological superconductors, is considered.


title: Spectra of Bethe subalgebras of $Y(\mathfrak{gl}_n)$ in tame representations

abstract: We study the eigenproblem for Bethe subalgebras of the Yangian $Y(\mathfrak{gl}_n)$ in tame representations, i.e. in finite dimensional representations which admit Gelfand-Tsetlin bases. Namely, we prove that for any tensor product of skew modules $V=\otimes_{i=1}^k V_{\lambda_i \setminus \mu_i}(z_i)$ over the Yangian $Y(\mathfrak{gl}_n)$ with generic $z_i$'s, the family of Bethe subalgebras $B(X)$ with $X$ being a regular element of the maximal torus of $GL_n$ (or, more generally, with $X \in \overline{M_{0,n+2}}$) acts with a cyclic vector on $V$. Moreover, for $X$ in the real form of $\overline{M_{0,n+2}}$ which is the closure of regular unitary diagonal matrices we show, that the family of subalgebras $B(X)$ acts with simple spectrum on $\otimes_{i=1}^k V_{\lambda_i \setminus \mu_i}(z_i)$ for generic $z_i$'s where all $V_{\lambda_i \setminus \mu_i}(z_i)$ are Kirillov-Reshetikhin modules. In the subsequent paper we will use this to define a KR-crystal structure on the spectrum of a Bethe subalgebra on $V$.


title: Aspect-Controllable Opinion Summarization

abstract: Recent work on opinion summarization produces general summaries based on a set of input reviews and the popularity of opinions expressed in them. In this paper, we propose an approach that allows the generation of customized summaries based on aspect queries (e.g., describing the location and room of a hotel). Using a review corpus, we create a synthetic training dataset of (review, summary) pairs enriched with aspect controllers which are induced by a multi-instance learning model that predicts the aspects of a document at different levels of granularity. We fine-tune a pretrained model using our synthetic dataset and generate aspect-specific summaries by modifying the aspect controllers. Experiments on two benchmarks show that our model outperforms the previous state of the art and generates personalized summaries by controlling the number of aspects discussed in them.


title: Graphite superlubricity enabled by triboinduced nanocontacts

abstract: Colloidal probe Atomic Force Microscopy allows to explore sliding states of vanishing friction, i.e. superlubricity, in mesoscopic graphite contacts. Superlubricity is known to appear upon formation of a triboinduced transfer layer, originated by material transfer of graphene flakes from the graphitic substrate to the colloidal probe. Previous studies suggest that friction vanishes due to crystalline incommensurability at the newly formed interface. However this picture still lacks several details, such as the roles of the tribolayer roughness and of loading conditions. Hereafter we gain deeper insight into the tribological response of micrometric silica beads sliding on graphite under ambient conditions. We show that the tribotransferred flakes behave as lubricious nanoasperities with a twofold role. First, they decrease the silica-graphite true contact area, in fact causing a breakdown of adhesion and friction by one order of magnitude. Second, they govern mechanical dissipation through the specific energy landscape experienced by the topographically-highest triboinduced nanoasperity. Remarkably, such contact junctions can undergo a load-driven atomic-scale transition from continuous superlubric sliding to dissipative stick-slip, that agrees with the single-asperity Prandtl-Tomlinson model. Superlubricity in mesoscopic silica-graphite junctions may therefore arise from the load-controlled competition between interfacial crystalline incommensurability and contact pinning effects at one dominant nanoasperity.


title: Learning to Bid in Contextual First Price Auctions

abstract: In this paper, we investigate the problem about how to bid in repeated contextual first price auctions. We consider a single bidder (learner) who repeatedly bids in the first price auctions: at each time $t$, the learner observes a context $x_t\in \mathbb{R}^d$ and decides the bid based on historical information and $x_t$. We assume a structured linear model of the maximum bid of all the others $m_t = \alpha_0\cdot x_t + z_t$, where $\alpha_0\in \mathbb{R}^d$ is unknown to the learner and $z_t$ is randomly sampled from a noise distribution $\mathcal{F}$ with log-concave density function $f$. We consider both \emph{binary feedback} (the learner can only observe whether she wins or not) and \emph{full information feedback} (the learner can observe $m_t$) at the end of each time $t$. For binary feedback, when the noise distribution $\mathcal{F}$ is known, we propose a bidding algorithm, by using maximum likelihood estimation (MLE) method to achieve at most $\widetilde{O}(\sqrt{\log(d) T})$ regret. Moreover, we generalize this algorithm to the setting with binary feedback and the noise distribution is unknown but belongs to a parametrized family of distributions. For the full information feedback with \emph{unknown} noise distribution, we provide an algorithm that achieves regret at most $\widetilde{O}(\sqrt{dT})$. Our approach combines an estimator for log-concave density functions and then MLE method to learn the noise distribution $\mathcal{F}$ and linear weight $\alpha_0$ simultaneously. We also provide a lower bound result such that any bidding policy in a broad class must achieve regret at least $\Omega(\sqrt{T})$, even when the learner receives the full information feedback and $\mathcal{F}$ is known.


title: Black holes in $f(\mathbb Q)$ Gravity

abstract: We systematically study the field equations of $f(\mathbb Q)$ gravity for spherically symmetric and stationary metric-affine spacetimes. Such spacetimes are described by a metric as well as a flat and torsionless affine connection. In the Symmetric Teleparallel Equivalent of GR (STEGR), the connection is pure gauge and hence unphysical. However, in the non-linear extension $f(\Q)$, it is promoted to a dynamical field which changes the physics. Starting from a general metric-affine geometry, we construct the most general static and spherically symmetric forms of the metric and the affine connection. We then use these symmetry reduced geometric objects to prove that the field equations of $f(\Q)$ gravity admit GR solutions as well as beyond-GR solutions, contrary to what has been claimed in the literature. We formulate precise criteria, under which conditions it is possible to obtain GR solutions and under which conditions it is possible to obtain beyond-GR solutions. We subsequently construct several perturbative corrections to the Schwarzschild solution for different choices of $f(\Q)$, which in particular include a hair stemming from the now dynamical affine connection. We also present an exact beyond-GR vacuum solution. Lastly, we apply this method of constructing spherically symmetric and stationary solutions to $f(\T)$ gravity, which reproduces similar solutions but without a dynamical connection.


title: When differential privacy meets NLP: The devil is in the detail

abstract: Differential privacy provides a formal approach to privacy of individuals. Applications of differential privacy in various scenarios, such as protecting users' original utterances, must satisfy certain mathematical properties. Our contribution is a formal analysis of ADePT, a differentially private auto-encoder for text rewriting (Krishna et al, 2021). ADePT achieves promising results on downstream tasks while providing tight privacy guarantees. Our proof reveals that ADePT is not differentially private, thus rendering the experimental results unsubstantiated. We also quantify the impact of the error in its private mechanism, showing that the true sensitivity is higher by at least factor 6 in an optimistic case of a very small encoder's dimension and that the amount of utterances that are not privatized could easily reach 100% of the entire dataset. Our intention is neither to criticize the authors, nor the peer-reviewing process, but rather point out that if differential privacy applications in NLP rely on formal guarantees, these should be outlined in full and put under detailed scrutiny.


title: Thermodynamic geometry of static and rotating regular black holes in conformal massive gravity

abstract: A version of massive gravity, namely conformal massive gravity, is employed to study the behavior of thermodynamic geometry for both the static and the rotating regular black holes. Whereas in thermodynamic geometry singularity can be interpreted as the second phase transition, seeing such phenomena as heat capacity behavior for both cases is investigated in detail. In addition, other thermodynamic quantities like the entropy and the temperature are analyzed as well. Another important property of the black holes is the so-called stability, in which utilizing the first phase transition of the heat capacity is detailed, for both cases, say the static and the spinning black holes. It is also interestingly figured out that, the behavior of temperature of the uncharged black holes in the presence of hair parameter, $\lambda$, can mimic the charged ones. {The effects of scalar charge, $Q$, and hair parameters with both positive and negative signs and how it affects key parameters in the study of black holes are investigated as well.} To investigate the properties of the black holes both regular thermodynamics and thermodynamic geometry approaches are considered. Then one can observe that aiming to realize the physics of the black holes, many complicated phenomena can be analyzed dramatically easier by considering the latter approach.


title: Magnetic Fields and Accreting Giant Planets around PDS 70

abstract: The recent high spatial/spectral resolution observations have enabled constraining formation mechanisms of giant planets, especially at the final stages. The current interpretation of such observations is that these planets undergo magnetospheric accretion, suggesting the importance of planetary magnetic fields. We explore the properties of accreting, magnetized giant planets surrounded by their circumplanetary disks, using the physical parameters inferred for PDS 70 b/c. We compute the magnetic field strength and the resulting spin rate of giant planets, and find that these planets may possess dipole magnetic fields of either a few 10 G or a few 100 G; the former is the natural outcome of planetary growth and radius evolution, while the resulting spin rate cannot reproduce the observations. For the latter, a consistent picture can be drawn, where strong magnetic fields induced by hot planetary interiors lead both to magnetospheric accretion and to spin-down due to disk locking. We also compute the properties of circumplanetary disks in the vicinity of these planets, taking into account planetary magnetic fields. The resulting surface density becomes very low, compared with the canonical models, implying the importance of radial movement of satellite-forming materials. Our model predicts a positive gradient of the surface density, which invokes the traps for both satellite migration and radially drifting dust particles. This work thus concludes that the final formation stages of giant planets are similar to those of low-mass stars such as brown dwarfs, as suggested by recent studies.


title: Melting of three-sublattice order in triangular lattice Ising antiferromagnets: Power-law order, $Z_6$ parafermionic multicriticality, and weakly first order transitions

abstract: The nature of the thermal melting process by which triangular-lattice Ising antiferromagnets lose their low-temperature ferrimagnetic three-sublattice order depends on the range of the interactions: It changes character when second and third neighbour ferromagnetic interactions become comparable to the nearest-neighbour antiferromagnetic coupling. We present a detailed numerical characterization of the corresponding threshold at which two-step melting of three-sublattice order gives way to a direct first-order transition at which this order is lost. The multicritical behaviour at this threshold is argued to be in the universality class of the $Z_6$ parafermion conformal field theory with central charge $c=5/4$. The presence of this multicritical threshold influences the melting behaviour and long-wavelength properties over a fairly large range of parameters, and at temperatures that are of the same order as the exchange interactions. It is therefore of potential experimental relevance in the context of easy-axis triangular lattice antiferromagnets that display such low temperature ordering.


title: Contest Design with Threshold Objectives

abstract: We study contests where the designer's objective is an extension of the widely studied objective of maximizing the total output: The designer gets zero marginal utility from a player's output if the output of the player is very low or very high. We model this using two objective functions: binary threshold, where a player's contribution to the designer's utility is 1 if her output is above a certain threshold, and 0 otherwise; and linear threshold, where a player's contribution is linear if her output is between a lower and an upper threshold, and becomes constant below the lower and above the upper threshold. For both of these objectives, we study (1) rank-order allocation contests that use only the ranking of the players to assign prizes and (2) general contests that may use the numerical values of the players' outputs to assign prizes. We characterize the contests that maximize the designer's objective and indicate techniques to efficiently compute them. We also prove that for the linear threshold objective, a contest that distributes the prize equally among a fixed number of top-ranked players offers a factor-2 approximation to the optimal rank-order allocation contest.


title: First Responders Got Wings: UAVs to the Rescue of Localization Operations in Beyond 5G Systems

abstract: Natural and human-made disasters have dramatically increased during the last decades. Given the strong relationship between first responders localization time and the final number of deaths, the modernization of search-and-rescue operations has become imperative. In this context, Unmanned Aerial Vehicles (UAVs)-based solutions are the most promising candidates to take up on the localization challenge by leveraging on emerging technologies such as: Artificial Intelligence (AI), Reconfigurable Intelligent Surfaces (RIS) and Orthogonal Time Frequency Space (OTFS) modulations. In this paper, we capitalize on such recently available techniques by shedding light on the main challenges and future opportunities to boost the localization performance of state-of-the-art techniques to give birth to unprecedentedly effective missing victims localization solutions.


abstract: IEEE BigData 2021 Cup: Soft Sensing at Scale is a data mining competition organized by Seagate Technology, in association with the IEEE BigData 2021 conference. The scope of this challenge is to tackle the task of classifying soft sensing data with machine learning techniques. In this paper we go into the details of the challenge and describe the data set provided to participants. We define the metrics of interest, baseline models, and describe approaches we found meaningful which may be a good starting point for further analysis. We discuss the results obtained with our approaches and give insights on what potential challenges participants may run into. Students, researchers, and anyone interested in working on a major industrial problem are welcome to participate in the challenge!


title: A Dynamic Population Model of Strategic Interaction and Migration under Epidemic Risk

abstract: In this paper, we show how a dynamic population game can model the strategic interaction and migration decisions made by a large population of agents in response to epidemic prevalence. Specifically, we consider a modified susceptible-asymptomatic-infected-recovered (SAIR) epidemic model over multiple zones. Agents choose whether to activate (i.e., interact with others), how many other agents to interact with, and which zone to move to in a time-scale which is comparable with the epidemic evolution. We define and analyze the notion of equilibrium in this game, and investigate the transient behavior of the epidemic spread in a range of numerical case studies, providing insights on the effects of the agents' degree of future awareness, strategic migration decisions, as well as different levels of lockdown and other interventions. One of our key findings is that the strategic behavior of agents plays an important role in the progression of the epidemic and can be exploited in order to design suitable epidemic control measures.


title: Encoding and ordering X-cactuses

abstract: Phylogenetic networks are a generalization of evolutionary or phylogenetic trees that are commonly used to represent the evolution of species which cross with one another. A special type of phylogenetic network is an {\em $X$-cactus}, which is essentially a cactus graph in which all vertices with degree less than three are labelled by at least one element from a set $X$ of species. In this paper, we present a way to {\em encode} $X$-cactuses in terms of certain collections of partitions of $X$ that naturally arise from $X$-cactuses. Using this encoding, we also introduce a partial order on the set of $X$-cactuses (up to isomorphism), and derive some structural properties of the resulting partially ordered set. This includes an analysis of some properties of its least upper and greatest lower bounds. Our results not only extend some fundamental properties of phylogenetic trees to $X$-cactuses, but also provides a new approach to solving topical problems in phylogenetic network theory such as deriving consensus networks.


title: The Convergence of Blockchain, IoT and 6G: Potential, Opportunities, Challenges and Research Roadmap

abstract: The world is undergoing a profound transformation with the advent of intelligent information era. 6G networks envisioned being the game changer in next generation wireless communication systems that will address the challenges of limited information speed escalated with the augmentation of billions of data applications encountered by the current fifth generation (5G) networks. Some key radical technologies in 6G together with existing 5G candidate schemes will guarantee the expected quality of experience (QoE) to attain ubiquitous wireless connectivity for the Internet of Everything (IoE) ranging from the telecom industry to digital smart industries. Blockchain technology (BCT) has gained significant attention due to undertake the decentralization, transparency, spectrum resource scarcity, inherent privacy and security, poor interoperability, confidentiality, and emerging smart applications domains including Industrial IoT and Industry 4.0. The mismatch between the requirements of many data intensive disruptive IoT applications and 5G network capabilities steered the demand of decentralized BCT based 6G architecture. Inspired by these facts, this paper studies an extensive survey to draw a new direction of blockchain integration into 6G mobile networks, IoT technologies, and smart industries focusing the potential merits and challenges in terms of infrastructure sharing, computational loads, latency, bandwidth overhead, business model, sustainability goals, and edge intelligence. We highlighted the convergence of IoT in blockchain to enable intelligent distribution in future industrial IoT and the technical model of 6G networks to realize the successful deployment of BCT schemes. This paper pointed out the current intriguing challenges, canvassed the mitigation techniques, and plausible future research opportunities that may benefit the pursuit of this vision.


title: A Sub-GeV Low Mass Hidden Dark Sector of $SU(2)_H \times U(1)_X$

abstract: We present a detailed study of the non-abelian vector dark matter candidate $W^\prime$ with a MeV-GeV low mass range, accompanied by a dark photon $A^\prime$ and a dark $Z^\prime$ of similar masses, in the context of a gauged two-Higgs-doublet model with the hidden gauge group that has the same structure as the Standard Model electroweak gauge group. The stability of dark matter is protected by an accidental discrete $Z_2$ symmetry ($h$-parity) which was usually imposed ad hoc by hand. We examine the model by taking into account various experimental constraints including dark photon searches at NA48, NA64, E141, $\nu$-CAL, BaBar and LHCb experiments, electroweak precision data from LEP, relic density from Planck satellite, direct (indirect) detection of dark matter from CRESST-III, DarkSide-50, XENON1T (Fermi-LAT), and collider physics from the LHC. The theoretical requirements of bounded from below of the scalar potential and tree level perturbative unitarity of the scalar sector are also imposed. The viable parameter space of the model consistent with all the constraints is exhibited. While a dark $Z^\prime$ can be the dominant contribution in the relic density due to resonant annihilation of dark matter, a dark photon is crucial to dark matter direct detection. We also demonstrate that the parameter space can be further probed by various sub-GeV direct dark matter experimental searches at CDEX, NEWS-G and SuperCDMS in the near future.


title: On Relationalist Reconstructions of Quantum Theory

abstract: Why I'm not happy with how Relational Quantum Mechanics has addressed the reconstruction of quantum theory, and why you shouldn't be either.


title: Limits on Relativistic Quantum Measurement

abstract: Requiring causality on measurements in quantum field theory seems to impose strong conditions on a self-adjoint operator to be really measurable. This may seem limiting and artificial in the operator language of algebraic quantum field theory (AQFT), but is essential for a truly relativistic theory. Recent publications attempt to deal with this issue by including the apparatus into the formalism, connecting AQFT with measurement theory, but other options have been suggested. In this essay, I discuss the causality conditions on self-adjoint operators both in the language of AQFT and in the language of quantum information theory. I then present measurement theory in AQFT, modelling the apparatus as a quantum field with coupling to the measured system restricted to a region of spacetime. I highlight how this approach leads to a causally well behaved theory. Finally, I attempt to formulate the causality conditions on measurements in the Feynman path integral approach, using the concept of decoherent histories. I claim that the path integral approach has problems with causality similar to the operator based approaches and that even here causality is an a posteriori condition.


title: Optimizing Quantum Variational Circuits with Deep Reinforcement Learning

abstract: Quantum Machine Learning (QML) is considered to be one of the most promising applications of near term quantum devices. However, the optimization of quantum machine learning models presents numerous challenges arising from the imperfections of hardware and the fundamental obstacles in navigating an exponentially scaling Hilbert space. In this work, we evaluate the potential of contemporary methods in deep reinforcement learning to augment gradient based optimization routines in quantum variational circuits. We find that reinforcement learning augmented optimizers consistently outperform gradient descent in noisy environments. All code and pretrained weights are available to replicate the results or deploy the models at: https://github.com/lockwo/rl_qvc_opt.


title: Isogenies of certain K3 surfaces of rank 18

abstract: We construct geometric isogenies between three types of two-parameter families of K3 surfaces of Picard rank 18. One is the family of Kummer surfaces associated with Jacobians of genus-two curves admitting an elliptic involution, another is the family of Kummer surfaces associated with the product of two non-isogeneous elliptic curves, and the third is the twisted Legendre pencil. The isogenies imply the existence of algebraic correspondences between these K3 surfaces and prove that the associated four-dimensional Galois representations are isomorphic. We also apply our result to several subfamilies of Picard rank 19. The result generalizes work of van Geemen and Top.


title: Phase diagram of a model for topological superconducting wires

abstract: We calculate the phase diagram of a model for topological superconducting wires with local s-wave pairing, spin-orbit coupling $\vec{\lambda}$ and magnetic field $\vec{B}$ with arbitrary orientations. This model is a generalized lattice version of the one proposed by Lutchyn $\textit{et al.}$ [Phys. Rev. Lett. $\textbf{105}$ 077001 (2010)] and Oreg $\textit{et al.}$ [Phys. Rev. Lett. $\textbf{105}$ 177002 (2010)], who considered $\vec{\lambda}$ perpendicular to $\vec{B}$. The model has a topological gapped phase with Majorana zero modes localized at the ends of the wires. We determine analytically the boundary of this phase. When the directions of the spin-orbit coupling and magnetic field are not perpendicular, in addition to the topological phase and the gapped non topological phase, a gapless superconducting phase appears.


title: Precision measurement of gravitational frequency shift of radio signals using Rao-Cramer estimates

abstract: A method has been developed for precision measurement of the gravitational frequency shift of communication radio signals between the spacecraft and the ground tracking station based on the maximum likelihood principle, using the Rao-Cramer limit estimates for the kinematic parameters associated with orbital motion. Numerical illustrations of the efficiency of the method are presented using the example of data obtained in experiments with the Spectr-R satellite as part of the VLBI system in the "Radioastron" mission. A compensatory "on-line" technique for suppressing Doppler and atmospheric noise has been implemented due to the presence of two modes of communication in gravity sessions: unidirectional (1w) and looped (2w). Recipes for reducing the magnitude of systematic errors are discussed.


title: Configuration spaces over singular spaces -- I. Dirichlet-Form and Metric Measure Geometry

abstract: We construct a canonical differential structure on the configuration space $\Upsilon$ over a singular base space $X$ and with a general invariant measure $\mu$ on $\Upsilon$. We present an analytic structure on $\Upsilon$, constructing a strongly local Dirichlet form $\mathcal E$ on $L^2(\Upsilon, \mu)$ for $\mu$ in a large class of probability measures. We then investigate the geometric structure of the extended metric measure space $\Upsilon$ endowed with the $L^2$-transportation extended distance $\mathsf{d}_{\Upsilon}$ and with the measure $\mu$. By establishing Rademacher- and Sobolev-to-Lipschitz-type properties for $\mathcal E$, we finally provide a complete identification of the analytic and the geometric structure -- the canonical differential structure induced on $\Upsilon$ by $X$ and $\mu$ -- showing that $\mathcal E$ coincides with the Cheeger energy of $(\Upsilon,\mathsf{d}_{\Upsilon},\mu)$ and that the intrinsic distance of $\mathcal E$ coincides with $\mathsf{d}_{\Upsilon}$. The class of base spaces to which our results apply includes sub-Riemannian manifolds, RCD spaces, and path/loop spaces over Riemannian manifolds; as for $\mu$ our results include quasi-Gibbs measures, in particular: Poisson measures, canonical Gibbs measures, as well as some determinantal/permanental point processes (sine$_\beta$, Airy$_\beta$, Bessel$_{\alpha,\beta}$, Ginibre). A number of applications to interacting particle systems and infinite-dimensional metric measure geometry are also discussed. In particular, we prove the universality of the $L^2$-transportation distance $\mathsf{d}_{\Upsilon}$ for the Varadhan short-time asymptotics for diffusions on $\Upsilon$, regardless of the choice of $\mu$. Many of our results are new even in the case of configuration spaces over Euclidean spaces.


title: Localization, monoid sets and K-theory

abstract: We develop the K-theory of sets with an action of a pointed monoid (or monoid scheme), analogous to the $K$-theory of modules over a ring (or scheme). In order to form localization sequences, we construct the quotient category of a nice regular category by a Serre subcategory.


title: On the Convergence of Decentralized Adaptive Gradient Methods

abstract: Adaptive gradient methods including Adam, AdaGrad, and their variants have been very successful for training deep learning models, such as neural networks. Meanwhile, given the need for distributed computing, distributed optimization algorithms are rapidly becoming a focal point. With the growth of computing power and the need for using machine learning models on mobile devices, the communication cost of distributed training algorithms needs careful consideration. In this paper, we introduce novel convergent decentralized adaptive gradient methods and rigorously incorporate adaptive gradient methods into decentralized training procedures. Specifically, we propose a general algorithmic framework that can convert existing adaptive gradient methods to their decentralized counterparts. In addition, we thoroughly analyze the convergence behavior of the proposed algorithmic framework and show that if a given adaptive gradient method converges, under some specific conditions, then its decentralized counterpart is also convergent. We illustrate the benefit of our generic decentralized framework on a prototype method, i.e., AMSGrad, both theoretically and numerically.


title: Gravity-induced phase phenomena in plate-rod colloidal mixtures

abstract: Gravity can affect colloidal suspensions since for micrometer-sized particles gravitational and thermal energies can be comparable over vertical length scales of a few millimeters. In mixtures, each species possesses a different buoyant mass, which can make experimental results counter-intuitive and difficult to interpret. Here, we revisit from a theoretical perspective iconic sedimentation-diffusion-equilibrium experiments on colloidal plate-rod mixtures by van der Kooij and Lekkerkerker. We reproduce their findings, including the observation of five different mesophases in a single cuvette. Using sedimentation path theory, we incorporate gravity into a microscopic theory for the bulk of a plate-rod mixture. We also show how to disentangle the effects of gravity from sedimentation experiments to obtain the bulk behavior and make predictions that can be experimentally tested. These include changes in the sequence by altering the sample height. We demonstrate that both buoyant mass ratio and sample height form control parameters to study bulk phase behavior.


title: Pseudospin-one particles in the time-periodic dice lattice: A new approach to transport control

abstract: The controlling of the transmission in the pseudospin-one Dirac-Weyl systems offers a rich tool to study new concepts of massive Dirac electron tunneling by means of a time-dependent potential. The time-periodic potential is one of the experimental techniques to have more control over the tunneling effect. In this paper, we study the transmission coefficient for different sidebands to obtain total transmission. We show how the super Klein tunneling under special condition is independent of the incidence angle, oscillation amplitude, frequency, and barrier width. We consider a band gap opening with different locations of the flat band and modulate the resonances by tuning free parameters in our system.


title: Gauge-Free Duality in Pure Square Spin Ice: Topological Currents and Monopoles

abstract: We consider a pure square spin ice, that is a square ice where only nearest neighbors are coupled. A gauge-free duality between the perpendicular and collinear structure leads to a natural description in terms of topological currents and charges as the relevant degrees of freedom. That, in turn, can be expressed via a continuous field theory where the discrete spins are subsumed into entropic interactions among charges and currents. This approach produces structure factors, correlations, and susceptibilities for spins, monopoles, and currents. It also generalizes to non-zero temperature the height formalism of the disordered ground state. The framework can be applied to a zoology of recent experimental results, especially realizations on quantum annealers and can be expanded to include longer range interactions.


title: Core motifs predict dynamic attractors in combinatorial threshold-linear networks

abstract: Combinatorial threshold-linear networks (CTLNs) are a special class of inhibition-dominated TLNs defined from directed graphs. Like more general TLNs, they display a wide variety of nonlinear dynamics including multistability, limit cycles, quasiperiodic attractors, and chaos. In prior work, we have developed a detailed mathematical theory relating stable and unstable fixed points of CTLNs to graph-theoretic properties of the underlying network. Here we find that a special type of fixed points, corresponding to core motifs, are predictive of both static and dynamic attractors. Moreover, the attractors can be found by choosing initial conditions that are small perturbations of these fixed points. This motivates us to hypothesize that dynamic attractors of a network correspond to unstable fixed points supported on core motifs. We tested this hypothesis on a large family of directed graphs of size $n=5$, and found remarkable agreement. Furthermore, we discovered that core motifs with similar embeddings give rise to nearly identical attractors. This allowed us to classify attractors based on structurally-defined graph families. Our results suggest that graphical properties of the connectivity can be used to predict a network's complex repertoire of nonlinear dynamics.


title: Unsupervised Conversation Disentanglement through Co-Training

abstract: Conversation disentanglement aims to separate intermingled messages into detached sessions, which is a fundamental task in understanding multi-party conversations. Existing work on conversation disentanglement relies heavily upon human-annotated datasets, which are expensive to obtain in practice. In this work, we explore to train a conversation disentanglement model without referencing any human annotations. Our method is built upon a deep co-training algorithm, which consists of two neural networks: a message-pair classifier and a session classifier. The former is responsible for retrieving local relations between two messages while the latter categorizes a message to a session by capturing context-aware information. Both networks are initialized respectively with pseudo data built from an unannotated corpus. During the deep co-training process, we use the session classifier as a reinforcement learning component to learn a session assigning policy by maximizing the local rewards given by the message-pair classifier. For the message-pair classifier, we enrich its training data by retrieving message pairs with high confidence from the disentangled sessions predicted by the session classifier. Experimental results on the large Movie Dialogue Dataset demonstrate that our proposed approach achieves competitive performance compared to the previous supervised methods. Further experiments show that the predicted disentangled conversations can promote the performance on the downstream task of multi-party response selection.


title: ExCode-Mixed: Explainable Approaches towards Sentiment Analysis on Code-Mixed Data using BERT models

abstract: The increasing use of social media sites in countries like India has given rise to large volumes of code-mixed data. Sentiment analysis of this data can provide integral insights into people's perspectives and opinions. Developing robust explainability techniques which explain why models make their predictions becomes essential. In this paper, we propose an adequate methodology to integrate explainable approaches into code-mixed sentiment analysis.


title: nnFormer: Interleaved Transformer for Volumetric Segmentation

abstract: Transformer, the model of choice for natural language processing, has drawn scant attention from the medical imaging community. Given the ability to exploit long-term dependencies, transformers are promising to help atypical convolutional neural networks to overcome their inherent shortcomings of spatial inductive bias. However, most of recently proposed transformer-based segmentation approaches simply treated transformers as assisted modules to help encode global context into convolutional representations. To address this issue, we introduce nnFormer, a 3D transformer for volumetric medical image segmentation. nnFormer not only exploits the combination of interleaved convolution and self-attention operations, but also introduces local and global volume-based self-attention mechanism to learn volume representations. Moreover, nnFormer proposes to use skip attention to replace the traditional concatenation/summation operations in skip connections in U-Net like architecture. Experiments show that nnFormer significantly outperforms previous transformer-based counterparts by large margins on three public datasets. Compared to nnUNet, nnFormer produces significantly lower HD95 and comparable DSC results. Furthermore, we show that nnFormer and nnUNet are highly complementary to each other in model ensembling.


title: On the impact of MDP design for Reinforcement Learning agents in Resource Management

abstract: The recent progress in Reinforcement Learning applications to Resource Management presents MDPs without a deeper analysis of the impacts of design decisions on agent performance. In this paper, we compare and contrast four different MDP variations, discussing their computational requirements and impacts on agent performance by means of an empirical analysis. We conclude by showing that, in our experiments, when using Multi-Layer Perceptrons as approximation function, a compact state representation allows transfer of agents between environments, and that transferred agents have good performance and outperform specialized agents in 80\% of the tested scenarios, even without retraining.


title: Determination of $1p$ and $2p$ stripping excitation functions for $^{16}$O+$^{142}$Ce using a Recoil Mass Spectrometer

abstract: We report the first direct measurement of differential transfer cross sections using a Recoil Mass Spectrometer. Absolute differential $1p$ and $2p$-stripping cross sections at $\theta_\mathrm{c.m.}=180^\circ$ have been determined for the system $^{16}$O+$^{142}$Ce by detecting the heavier target-like ions at the focal plane of the Heavy Ion Reaction Analyzer. Focal plane spectra have been compared with the results of a semi-microscopic Monte-Carlo simulation to unambiguously identify the transfer channels. Transmission efficiency of the target-like ions through the spectrometer has also been estimated using the simulation. The methodology adopted in this work can be applied to other recoil separators. The measured excitation functions for the reactions $^{142}\mathrm{Ce(}^{16}\mathrm{O,}^{15}\mathrm{N)}^{143}\mathrm{Pr}$ and $^{142}\mathrm{Ce(}^{16}\mathrm{O,}^{14}\mathrm{C)}^{144}\mathrm{Nd}$ have been compared with coupled reaction channel calculations. An excellent matching between measurement and theory has been obtained. For $1p$-stripping, major contribution to the cross section has been found to be the transfer of a proton from $^{16}\mathrm{O}$ to the $2d_{\frac{5}{2}}$ excited state of $^{143}\mathrm{Pr}$, leaving behind $^{15}\mathrm{N}$ in the $1p_{\frac{1}{2}}$ ground state. Transfer of a cluster of two protons from $^{16}\mathrm{O}$ to the $2^{+}$ excited state of $^{144}\mathrm{Nd}$, resulting in $^{14}\mathrm{C}$ in the $0^{+}$ ground state, appears to be the most probable cause for $2p$-stripping. Measured transfer probabilities for $1p$ and $2p$ channels have been compared with Time-Dependent Hartree-Fock calculations. Proton stripping channels are found to be more favourable compared to neutron pick-up channels. However, the theory overpredicts measurement hinting at the need for extended approaches with explicit treatment of pairing correlations in the calculations.


title: Adaptive variational Bayes: Optimality, computation and applications

abstract: In this paper, we explore adaptive inference based on variational Bayes. Although several studies have been conducted to analyze the contraction properties of variational posteriors, there is still a lack of a general and computationally tractable variational Bayes method that performs adaptive inference. To fill this gap, we propose a novel adaptive variational Bayes framework, which can operate on a collection of models. The proposed framework first computes a variational posterior over each individual model separately and then combines them with certain weights to produce a variational posterior over the entire model. It turns out that this combined variational posterior is the closest member to the posterior over the entire model in a predefined family of approximating distributions. We show that the adaptive variational Bayes attains optimal contraction rates adaptively under very general conditions. We also provide a methodology to maintain the tractability and adaptive optimality of the adaptive variational Bayes even in the presence of an enormous number of individual models, such as sparse models. We apply the general results to several examples, including deep learning and sparse factor models, and derive new and adaptive inference results. In addition, we characterize an implicit regularization effect of variational Bayes and show that the adaptive variational posterior can utilize this.


title: Long-term multi-band photometric monitoring of Mrk 501

abstract: Radio-to-TeV observations of the bright nearby (z=0.034) blazar Markarian 501 (Mrk 501), performed from December 2012 to April 2018, are used to study the emission mechanisms in its relativistic jet. We examined the multi-wavelength variability and the correlations of the light curves obtained by eight different instruments, including the First G-APD Cherenkov Telescope (FACT), observing Mrk 501 in very high-energy (VHE) gamma-rays at TeV energies. We identified individual TeV and X-ray flares and found a sub-day lag between variability in these two bands. Simultaneous TeV and X-ray variations with almost zero lag are consistent with synchrotron self-Compton (SSC) emission, where TeV photons are produced through inverse Compton scattering. The characteristic time interval of 5-25 days between TeV flares is consistent with them being driven by Lense-Thirring precession.


title: Bivariate collocation for computing $R_{0}$ in epidemic models with two structures

abstract: Structured epidemic models can be formulated as first-order hyperbolic PDEs, where the "spatial" variables represent individual traits, called structures. For models with two structures, we propose a numerical technique to approximate $R_{0}$, which measures the transmissibility of an infectious disease and, rigorously, is defined as the dominant eigenvalue of a next-generation operator. Via bivariate collocation and cubature on tensor grids, the latter is approximated with a finite-dimensional matrix, so that its dominant eigenvalue can easily be computed with standard techniques. We use test examples to investigate experimentally the behavior of the approximation: the convergence order appears to be infinite when the corresponding eigenfunction is smooth, and finite for less regular eigenfunctions. To demonstrate the effectiveness of the technique for more realistic applications, we present a new epidemic model structured by demographic age and immunity, and study the approximation of $R_{0}$ in some particular cases of interest.


title: COCO Denoiser: Using Co-Coercivity for Variance Reduction in Stochastic Convex Optimization

abstract: First-order methods for stochastic optimization have undeniable relevance, in part due to their pivotal role in machine learning. Variance reduction for these algorithms has become an important research topic. In contrast to common approaches, which rarely leverage global models of the objective function, we exploit convexity and L-smoothness to improve the noisy estimates outputted by the stochastic gradient oracle. Our method, named COCO denoiser, is the joint maximum likelihood estimator of multiple function gradients from their noisy observations, subject to co-coercivity constraints between them. The resulting estimate is the solution of a convex Quadratically Constrained Quadratic Problem. Although this problem is expensive to solve by interior point methods, we exploit its structure to apply an accelerated first-order algorithm, the Fast Dual Proximal Gradient method. Besides analytically characterizing the proposed estimator, we show empirically that increasing the number and proximity of the queried points leads to better gradient estimates. We also apply COCO in stochastic settings by plugging it in existing algorithms, such as SGD, Adam or STRSAGA, outperforming their vanilla versions, even in scenarios where our modelling assumptions are mismatched.


title: Pareto Robust optimization on Euclidean vector spaces

abstract: Pareto efficiency for robust linear programs was introduced by Iancu and Trichakis in [9]. We generalize their approach and theoretical results to robust optimization problems in Euclidean spaces with affine uncertainty. Additionally, we demonstrate the value of this approach in an exemplary manner in the area of robust semidefinite programming (SDP). In particular, we prove that computing a Pareto robustly optimal solution for a robust SDP is tractable and illustrate the benefit of such solutions at the example of the maximal eigenvalue problem. Furthermore, we modify the famous algorithm of Goemans and Williamson [8] in order to compute cuts for the robust max-cut problem that yield an improved approximation guarantee in non-worst-case scenarios.


title: Route Reservation Architecture in Tandem Transportation Networks with the Asymmetric Cell Transmission Model

abstract: Traffic congestion in big cities has been proven to be a difficult problem with suboptimal effects in terms of driver delay and frustration, cost, and impact on the environment. In principle, many transportation networks lack a unified framework, which will coordinate the traffic in such a manner, in order to suppress congestion and at the same time improve the travel time of the users situated in it. The rapid advancements in information, communication, and computation technologies have given rise to more elaborate modeling frameworks, aiming to act as the coordination unit necessary to counter the issue of congestion in real-time conditions. Such actions might have an adverse effect on the efficiency of the network, prospectively leading to greater waiting time intervals for each individual driver. We propose a macroscopic model equipped with an underlying reservation feature, known as Route Reservation Architecture (RRA). Vehicles enter the network (mainstream-wise or from the on-ramps) as long as this inflow does not incur a density that exceeds the network's critical density. Those vehicles prospectively exceeding the critical density are stored as queues at the origin of the motorway stretch and within the on-ramps. Once there is sufficient space, for those vehicles to be accommodated by the respective cell of the network, they are discharged from their queueing instance at their respective origin, moving towards their assigned path freely. In previous works, this architecture has been only applied in microscopic simulations in the context of urban networks without the influence of source terms (on-ramps). When the critical density of the stretch is reached, the reservations are activated, instigating a waiting interval to the vehicles stored at queues, allowing vehicles only to enter the network at a later time instant, such that the critical density is not crossed.


title: Routing Emergency Vehicles in Arterial Road Networks using Real-time Mixed Criticality Systems*

abstract: Reducing the response time of Emergency Vehicles (EVs) has an undoubted advantage in saving life and property. Implementing pre-emption can aid in achieving it. EVs get unobstructed movement via pre-emption, usually by altering traffic signals and giving a green wave throughout the route. This approach of absolute pre-emption effects adversely on regular traffic by imposing unnecessary waiting. In this paper, we propose a novel emergency vehicle pre-emption (EVP) algorithm implemented in the Vehicular Ad-hoc Network (VANET) that can reduce the imposed undesirable waiting time, but still ascertains EVs meet target response time. We introduce mixed-criticality real-time system scheduling concept where different level of emergencies is mapped with different criticality levels and assign certain success assurance level to respective criticality. We implemented the EVP algorithm for an arterial traffic network and leveraged the use of valuable information that gets transmitted via VANET to make critical decisions. The proposed algorithm can significantly reduce the average waiting time of regular traffic. It also ascertains all EVs with different level of criticality meet target response time respective to their assurance level.