title: Silicon crystals for steering of high-intensity particle beams at ultra-high energy accelerators

abstract: Experimental results and simulation models show that crystals might play a relevant role for the development of new generations of high-energy and high-intensity particle accelerators and might disclose innovative possibilities at existing ones. In this paper we describe the most advanced manufacturing techniques of crystals suitable for operations at ultra-high energy and ultra-high intensity particle accelerators, reporting as an example of potential applications the collimation of the particle beams circulating in the Large Hadron Collider at CERN, which will be upgraded through the addition of bent crystals in the frame of the High Luminosity Large Hadron Collider project.


title: Modelling Excess Mortality in Covid-19-like Epidemics

abstract: We develop an agent-based model to assess the cumulative number of deaths during hypothetical Covid-19-like epidemics for various non-pharmaceutical intervention strategies. We consider local and non-local modes of disease transmission. The first simulates transmission through social contacts in the vicinity of the place of residence while the second through social contacts in public places: schools, hospitals, airports, etc., where many people meet, who live in remote geographic locations. Epidemic spreading is modeled as a discrete-time stochastic process on random geometric networks. We use the Monte-Carlo method in the simulations. The~following assumptions are made. The basic reproduction number is 2.5 and the infectious period lasts approximately ten days. Infections lead to SARS in about one percent of cases, which are likely to lead to respiratory default and death, unless the patient receives an appropriate medical treatment. The~healthcare system capacity is simulated by the availability of respiratory ventilators or intensive care beds. Some parameters of the model, like mortality rates or the number of respiratory ventilators per 100000 inhabitants, are chosen to simulate the real values for the USA and Poland. In the simulations we compare `do-nothing' strategy with mitigation strategies based on social distancing and reducing social mixing. We study epidemics in the pre-vaccine era, where immunity is obtained only by infection. The model applies only to epidemics for which reinfections are rare and can be neglected. The results of the simulations show that strategies that slow the development of an epidemic too much in the early stages do not significantly reduce the overall number of deaths in the long term, but increase the duration of the epidemic. In particular, a~hybrid strategy where lockdown is held for some time and is then completely released, is inefficient.


title: A Polynomial Kernel for Line Graph Deletion

abstract: The line graph of a graph $G$ is the graph $L(G)$ whose vertex set is the edge set of $G$ and there is an edge between $e,f\in E(G)$ if $e$ and $f$ share an endpoint in $G$. A graph is called line graph if it is a line graph of some graph. We study the Line-Graph-Edge Deletion problem, which asks whether we can delete at most $k$ edges from the input graph $G$ such that the resulting graph is a line graph. More precisely, we give a polynomial kernel for Line-Graph-Edge Deletion with $\mathcal{O}(k^{5})$ vertices. This answers an open question posed by Falk H\"{u}ffner at Workshop on Kernels (WorKer) in 2013.


title: Self-Attention Networks for Intent Detection

abstract: Self-attention networks (SAN) have shown promising performance in various Natural Language Processing (NLP) scenarios, especially in machine translation. One of the main points of SANs is the strength of capturing long-range and multi-scale dependencies from the data. In this paper, we present a novel intent detection system which is based on a self-attention network and a Bi-LSTM. Our approach shows improvement by using a transformer model and deep averaging network-based universal sentence encoder compared to previous solutions. We evaluate the system on Snips, Smart Speaker, Smart Lights, and ATIS datasets by different evaluation metrics. The performance of the proposed model is compared with LSTM with the same datasets.


title: When the regularized Lovelock tensors are kinetically coupled to scalar field

abstract: The recently proposed regularized Lovelock tensors are kinetically coupled to the scalar field. The resulting equation of motion is second order. In particular, it is found that when the $p=3$ regularized Lovelock tensor is kinetically coupled to the scalar field, the scalar field is the potential candidate of cosmic dark energy.


title: Reproducibility of `COST Reference Microplasma Jets'

abstract: Atmospheric pressure plasmas have been ground-breaking for plasma science and technologies, due to their significant application potential in many fields, including medicinal, biological, and environmental applications. This is predominantly due to their efficient production and delivery of chemically reactive species under ambient conditions. One of the challenges in progressing the field is comparing plasma sources and results across the community and the literature. To address this a reference plasma source was established during the `Biomedical Applications of Atmospheric Pressure Plasmas' EU COST Action MP1101. It is crucial that reference sources are reproducible. Here, we present the reproducibility and variance across multiple sources through examining various characteristics, including: absolute atomic oxygen densities, absolute ozone densities, electrical characteristics, optical emission spectroscopy, temperature measurements, and bactericidal activity. The measurements demonstrate that the tested COST jets are reproducible, within the standard deviation, for each measurement.


title: A lateral semicircular canal segmentation based geometric calibration for human temporal bone CT Image

abstract: Computed Tomography (CT) of the temporal bone has become an important method for diagnosing ear diseases. Due to the different posture of the subject and the settings of CT scanners, the CT image of the human temporal bone should be geometrically calibrated to ensure the symmetry of the bilateral anatomical structure. Manual calibration is a time-consuming task for radiologists and an important pre-processing step for further computer-aided CT analysis. We propose an automatic calibration algorithm for temporal bone CT images. The lateral semicircular canals (LSCs) are segmented as anchors at first. Then, we define a standard 3D coordinate system. The key step is the LSC segmentation. We design a novel 3D LSC segmentation encoder-decoder network, which introduces a 3D dilated convolution and a multi-pooling scheme for feature fusion in the encoding stage. The experimental results show that our LSC segmentation network achieved a higher segmentation accuracy. Our proposed method can help to perform calibration of temporal bone CT images efficiently.


title: Surface-acoustic-wave-induced unconventional superconducting pairing

abstract: Unconventional superconductivity is usually associated with symmetry breaking in the system. Here we consider a simple setup consisting of a solid state material with conduction electrons and an applied surface acoustic wave (SAW), that can break time and spatial translation symmetries. We study the symmetries of the possible SAW-induced order parameters, showing that even-frequency spin-triplet odd-parity order parameter can occur. We suggest different methods of how to engineer the symmetries of the order parameters using SAWs and the applications of such setups.


title: VPNet: Variable Projection Networks

abstract: We introduce VPNet, a novel model-driven neural network architecture based on variable projection (VP). Applying VP operators to neural networks results in learnable features, interpretable parameters, and compact network structures. This paper discusses the motivation and mathematical background of VPNet and presents experiments. The VPNet approach was evaluated in the context of signal processing, where we classified a synthetic dataset and real electrocardiogram (ECG) signals. Compared to fully connected and one-dimensional convolutional networks, VPNet offers fast learning ability and good accuracy at a low computational cost of both training and inference. Based on these advantages and the promising results obtained, we anticipate a profound impact on the broader field of signal processing, in particular on classification, regression and clustering problems.


title: Stochastic Geometry-based Analysis of LEO Satellite Communication Systems

abstract: This letter studies the performance of a low-earth orbit (LEO) satellite communication system where the locations of the LEO satellites are modeled as a binomial point process (BPP) on a spherical surface. In particular, we study the user coverage probability for a scenario where satellite gateways (GWs) are deployed on the ground to act as a relay between the users and the LEO satellites. We use tools from stochastic geometry to derive the coverage probability for the described setup assuming that LEO satellites are placed at n different altitudes, given that the number of satellites at each altitude ak is Nk for all k. To resemble practical scenarios where satellite communication can play an important role in coverage enhancement, we compare the performance of the considered setup with a scenario where the users are solely covered by a fiber-connected base station (referred to as anchored base station or ABS in the rest of the paper) at a relatively far distance, which is a common challenge in rural and remote areas. Using numerical results, we show the performance gain, in terms of coverage probability, at rural and remote areas when LEO satellite communication systems are adopted. Finally, we draw multiple system-level insights regarding the density of GWs required to outperform the ABS, as well as the number of LEO satellites and their altitudes.


title: Saffman-Taylor Fingers at Intermediate Noise

abstract: We study Saffman-Taylor flow in the presence of intermediate noise numerically by using both a boundary-integral approach as well as the Kadanoff-Liang modified Diffusion-Limited Aggregation model that incorporates surface tension and reduced noise. For little to no noise, both models result reproduce the well-known Saffman-Taylor finger. We compare both models in the region of intermediate noise where we get occasional tip-splitting events, focusing on the ensemble-average. We show that as the noise in the system is increased, the mean behavior in both models approaches the $\cos^2(\pi y/W)$ transverse density profile far behind the leading front. We also investigate how the noise scales and affects both models.


title: Exact Solutions of the DKP Oscillator in 3D Spaces with Extended Uncertainty Principle

abstract: We present the exact solution of the three-dimensional Duffin--Kemmer--Petiau oscillator for both spin 0 and spin 1 cases, with the presence of minimal uncertainty in momentum in anti--de Sitter model. We use the representation of vector spherical harmonics and the Nikiforov--Uvarov method to determine exactly the energy eigenvalues and the eigenfunctions in all cases. Our study of the energy spectrum allows us to define a new interpretation of natural and unnatural parity states of the vector particle and we show the crucial role played by the spin--orbit coupling in this differentiation between the parities.


title: EdgeKV: Decentralized, scalable, and consistent storage for the edge

abstract: Edge computing moves the computation closer to the data and the data closer to the user to overcome the high latency communication of cloud computing. Storage at the edge allows data access with high speeds that enable latency-sensitive applications in areas such as autonomous driving and smart grid. However, several distributed services are typically designed for the cloud and building an efficient edge-enabled storage system is challenging because of the distributed and heterogeneous nature of the edge and its limited resources. In this paper, we propose EdgeKV, a decentralized storage system designed for the network edge. EdgeKV offers fast and reliable storage, utilizing data replication with strong consistency guarantees. With a location-transparent and interface-based design, EdgeKV can scale with a heterogeneous system of edge nodes. We implement a prototype of the EdgeKV modules in Golang and evaluate it in both the edge and cloud settings on the Grid'5000 testbed. We utilize the Yahoo! Cloud Serving Benchmark (YCSB) to analyze the system's performance under realistic workloads. Our evaluation results show that EdgeKV outperforms the cloud storage setting with both local and global data access with an average write response time and throughput improvements of 26% and 19% respectively under the same settings. Our evaluations also show that EdgeKV can scale with the number of clients, without sacrificing performance. Finally, we discuss the energy efficiency improvement when utilizing edge resources with EdgeKV instead of a centralized cloud.


title: Rethinking Positional Encoding in Language Pre-training

abstract: In this work, we investigate the positional encoding methods used in language pre-training (e.g., BERT) and identify several problems in the existing formulations. First, we show that in the absolute positional encoding, the addition operation applied on positional embeddings and word embeddings brings mixed correlations between the two heterogeneous information resources. It may bring unnecessary randomness in the attention and further limit the expressiveness of the model. Second, we question whether treating the position of the symbol \texttt{[CLS]} the same as other words is a reasonable design, considering its special role (the representation of the entire sentence) in the downstream tasks. Motivated from above analysis, we propose a new positional encoding method called \textbf{T}ransformer with \textbf{U}ntied \textbf{P}ositional \textbf{E}ncoding (TUPE). In the self-attention module, TUPE computes the word contextual correlation and positional correlation separately with different parameterizations and then adds them together. This design removes the mixed and noisy correlations over heterogeneous embeddings and offers more expressiveness by using different projection matrices. Furthermore, TUPE unties the \texttt{[CLS]} symbol from other positions, making it easier to capture information from all positions. Extensive experiments and ablation studies on GLUE benchmark demonstrate the effectiveness of the proposed method. Codes and models are released at https://github.com/guolinke/TUPE.


title: Dark side of Weyl gravity

abstract: With the help of a functional renormalization group, we study the dynamical breakdown of scale invariance in quantum Weyl gravity by starting from the UV fixed point that we assume to be Gaussian. To this end, we resort to two classes of Bach-flat backgrounds, namely maximally symmetric spacetimes and Ricci-flat backgrounds in the improved one-loop scheme. We show that apart from a genuine IR fixed point that is reached at a zero value of the running scale, the renormalization group flow also exhibits bouncing behavior. We demonstrate that the IR fixed point found is IR-stable in the space of the considered couplings. As a next step, we analyze physics in the broken phase. In particular, we show that in the low-energy sector of the broken phase, the theory looks like Starobinsky $f(R)$ gravity with a gravi-cosmological constant that has a negative sign in comparison to the usual matter-induced cosmological constant. We discuss implications for cosmic inflation and highlight a non-trivial relation between Starobinsky's parameter and the gravi-cosmological constant. Salient issues, including the scheme independence of the IR fixed point and the role of trace anomaly, are also discussed.


title: Exact Solutions for a Quantum Ring with a Dipolar Impurity

abstract: We study analytically a system made up of a quantum ring with a dipolar impurity and under the effect of an Aharonov-Bohm field. We calculate the exact values of the energies and we also get the exact expressions of the wave functions.


title: A Survey Mobility Management in 5G Networks

abstract: For the next wireless communication systems, the major challenges are to provide ubiquitous wireless access abilities and maintain the quality of service and seamless mobility management for mobile communication devices in heterogeneous networks. Due to the rapid growth of mobile users, this demand becomes more challenging where the users always require seamless connectivity while they move to other places at any time. As the number of users increases, the network load also increases, the handover process needs to be performed in an efficient way. However, in many situations, the handover blocking, and unnecessary handover frequently happens, then affect the network and reduces its performance. The problem arises with the movement of mobile users between base stations while the link connectivity becomes weaker and the mobile node tries to switch to another base station to have a better link quality during a call for higher QoS (Quality of Service). In this survey, it is aimed to gather together the studies that help to improve mobility management processes in 5G (Fifth Generation) cellular networks.


title: Answer Ranking for Product-Related Questions via Multiple Semantic Relations Modeling

abstract: Many E-commerce sites now offer product-specific question answering platforms for users to communicate with each other by posting and answering questions during online shopping. However, the multiple answers provided by ordinary users usually vary diversely in their qualities and thus need to be appropriately ranked for each question to improve user satisfaction. It can be observed that product reviews usually provide useful information for a given question, and thus can assist the ranking process. In this paper, we investigate the answer ranking problem for product-related questions, with the relevant reviews treated as auxiliary information that can be exploited for facilitating the ranking. We propose an answer ranking model named MUSE which carefully models multiple semantic relations among the question, answers, and relevant reviews. Specifically, MUSE constructs a multi-semantic relation graph with the question, each answer, and each review snippet as nodes. Then a customized graph convolutional neural network is designed for explicitly modeling the semantic relevance between the question and answers, the content consistency among answers, and the textual entailment between answers and reviews. Extensive experiments on real-world E-commerce datasets across three product categories show that our proposed model achieves superior performance on the concerned answer ranking task.


title: A Benson-Type Algorithm for Bounded Convex Vector Optimization Problems with Vertex Selection

abstract: We present an algorithm for approximately solving bounded convex vector optimization problems. The algorithm provides both an outer and an inner polyhedral approximation of the upper image. It is a modification of the primal algorithm presented by L\"ohne, Rudloff, and Ulus in 2014. There, vertices of an already known outer approximation are successively cut off to improve the approximation error. We propose a new and efficient selection rule for deciding which vertex to cut off. Numerical examples are provided which illustrate that this method may solve fewer scalar problems overall and therefore may be faster while achieving the same approximation quality.


title: Exact Solutions of D-dimensional Klein-Gordon Oscillator with Snyder-de Sitter Algebra

abstract: We study the effects of Snyder-de Sitter commutation relations on relativistic bosons by solving analytically in the momentum space representation the Klein-Gordon oscillator in arbitrary dimensions. The exact bound states spectrum and the corresponding momentum space wave functions are obtained using Gegenbauer polynomials in one dimension space and Jacobi polynomials in D dimensions case. Finally, we study the thermodynamic properties of the system in the high temperature regime where we found that the corrections increase the free energy but decrease the energy, the entropy and the specific heat which is no longer constant. This work extends the part concerning the Klein-Gordon oscillator for the Snyder-de Sitter case studied in two-dimensional space in J. Math. Phys. 60, 013505 (2019).


title: A Multilevel Approach to Training

abstract: We propose a novel training method based on nonlinear multilevel minimization techniques, commonly used for solving discretized large scale partial differential equations. Our multilevel training method constructs a multilevel hierarchy by reducing the number of samples. The training of the original model is then enhanced by internally training surrogate models constructed with fewer samples. We construct the surrogate models using first-order consistency approach. This gives rise to surrogate models, whose gradients are stochastic estimators of the full gradient, but with reduced variance compared to standard stochastic gradient estimators. We illustrate the convergence behavior of the proposed multilevel method to machine learning applications based on logistic regression. A comparison with subsampled Newton's and variance reduction methods demonstrate the efficiency of our multilevel method.


title: Exploiting Diffuse Multipath in 5G SLAM

abstract: 5G millimeter wave (mmWave) signals can be used to jointly localize the receiver and map the propagation environment in vehicular networks, which is a typical simultaneous localization and mapping (SLAM) problem. Mapping the environment is challenging, due to measurements comprising both specular and diffuse multipath components, and diffuse multipath is usually considered as a perturbation. We here propose a novel method to utilize all available multipath signals from each landmark for mapping and incorporate this into a Poisson multi-Bernoulli mixture for the 5G SLAM problem. Simulation results demonstrate the efficacy of the proposed scheme.


title: Layer Sparsity in Neural Networks

abstract: Sparsity has become popular in machine learning, because it can save computational resources, facilitate interpretations, and prevent overfitting. In this paper, we discuss sparsity in the framework of neural networks. In particular, we formulate a new notion of sparsity that concerns the networks' layers and, therefore, aligns particularly well with the current trend toward deep networks. We call this notion layer sparsity. We then introduce corresponding regularization and refitting schemes that can complement standard deep-learning pipelines to generate more compact and accurate networks.


title: RISE Controller Tuning and System Identification Through Machine Learning for Human Lower Limb Rehabilitation via Neuromuscular Electrical Stimulation

abstract: Neuromuscular electrical stimulation (NMES) has been effectively applied in many rehabilitation treatments of individuals with spinal cord injury (SCI). In this context, we introduce a novel, robust, and intelligent control-based methodology to closed-loop NMES systems. Our approach utilizes a robust control law to guarantee system stability and machine learning tools to optimize both the controller parameters and system identification. Regarding the latter, we introduce the use of past rehabilitation data to build more realistic data-driven identified models. Furthermore, we apply the proposed methodology for the rehabilitation of lower limbs using a control technique named the robust integral of the sign of the error (RISE), an offline improved genetic algorithm optimizer, and neural network models. Although in the literature, the RISE controller presented good results on healthy subjects, without any fine-tuning method, a trial and error approach would quickly lead to muscle fatigue for individuals with SCI. In this paper, for the first time, the RISE controller is evaluated with two paraplegic subjects in one stimulation session and with seven healthy individuals in at least two and at most five sessions. The results showed that the proposed approach provided a better control performance than empirical tuning, which can avoid premature fatigue on NMES-based clinical procedures.


title: Five-dimensional para-CR manifolds and contact projective geometry in dimension three

abstract: We study invariant properties of $5$-dimensional para-CR structures whose Levi form is degenerate in precisely one direction and which are $2$-nondegenerate. We realize that two, out of three, primary (basic) para-CR invariants of such structures are the classical differential invariants known to Monge (1810) and to Wuenschmann (1905) \[ M(G) := 40G_{ppp}^3-45G_{pp}G_{ppp}G_{pppp}+9G_{pp}^2G_{ppppp}, \quad W(H) := 9D^2H_r-27DH_p-18H_rDH_r+18H_pH_r+4H_r^3+54H_z. \] The vanishing $M(G) \equiv 0$ provides a local necessary and sufficient condition for the graph of a function in the $(p,G)$-plane to be contained in a conic, while the vanishing $W(H) \equiv 0$ gives an if-and-only-if condition for a 3rd order ODE to define a natural Lorentzian geometry on the space of its solutions. Mainly, we give a geometric interpretation of the third basic invariant of our class of para-CR structures, the simplest one, of lowest order, and of mixed nature $N(G,H):=2G_{ppp}+G_{pp}H_{rr}$. We establish that the vanishing $N(G,H) \equiv 0$ gives an if-and-only-if condition for the two $3$-dimensional quotients of the para-CR manifold by its two canonical integrable rank-$2$ distributions, to be equipped with contact projective geometries. A curious transformation between the Wuenschmann invariant and the Monge invariant, first noted by us in arXiv:2003.08166, is also discussed, and its mysteries are further revealed.


title: Localization Uncertainty Estimation for Anchor-Free Object Detection

abstract: Since many safety-critical systems, such as surgical robots and autonomous driving cars operate in unstable environments with sensor noise and incomplete data, it is desirable for object detectors to take the localization uncertainty into account. However, there are several limitations of the existing uncertainty estimation methods for anchor-based object detection. 1) They model the uncertainty of the heterogeneous object properties with different characteristics and scales, such as location (center point) and scale (width, height), which could be difficult to estimate. 2) They model box offsets as Gaussian distributions, which is not compatible with the ground truth bounding boxes that follow the Dirac delta distribution. 3) Since anchor-based methods are sensitive to anchor hyper-parameters, their localization uncertainty could also be highly sensitive to the choice of hyper-parameters. To tackle these limitations, we propose a new localization uncertainty estimation method called UAD for anchor-free object detection. Our method captures the uncertainty in four directions of box offsets (left, right, top, bottom) that are homogeneous, so that it can tell which direction is uncertain, and provide a quantitative value of uncertainty in [0, 1]. To enable such uncertainty estimation, we design a new uncertainty loss, negative power log-likelihood loss, to measure the localization uncertainty by weighting the likelihood loss by its IoU, which alleviates the model misspecification problem. Furthermore, we propose an uncertainty-aware focal loss for reflecting the estimated uncertainty to the classification score. Experimental results on COCO datasets demonstrate that our method significantly improves FCOS, by up to 1.8 points, without sacrificing computational efficiency.


title: Combination of Tensor Network States and Green's function Monte Carlo

abstract: We propose an approach to study the ground state of quantum many-body systems in which Tensor Network States (TNS), specifically Projected Entangled Pair States (PEPS), and Green's function Monte Carlo (GFMC) are combined. PEPS, by design, encode the area law which governs the scaling of entanglement entropy in quantum systems with short range interactions but are hindered by the high computational complexity scaling with bond dimension (D). GFMC is a highly efficient method, but it usually suffers from the infamous negative sign problem which can be avoided by the fixed node approximation in which a guiding wave function is utilized to modify the sampling process. The trade-off for the absence of negative sign problem is the introduction of systematic error by guiding wave function. In this work, we combine these two methods, PEPS and GFMC, to take advantage of both of them. PEPS are very accurate variational wave functions, while at the same time, only contractions of single-layer tensor network are necessary in GFMC, which reduces the cost substantially. Moreover, energy obtained in GFMC is guaranteed to be variational and lower than the variational energy of the guiding PEPS wave function. Benchmark results of $J_1$-$J_2$ Heisenberg model on square lattice are provided.


title: Root finding algorithms and persistence of Jordan centrality in growing random trees

abstract: We consider models of growing random trees $\{\mathcal{T}_f(n):n\geq 1\}$ with model dynamics driven by an attachment function $f:\mathbb{Z}_+\to \mathbb{R}_+$. At each stage a new vertex enters the system and connects to a vertex $v$ in the current tree with probability proportional to $f(\text{degree}(v))$. The main goal of this study is to understand the performance of root finding algorithms. A large body of work (e.g. the work of Bubeck, Devroye and Lugosi or Jog and Loh) has emerged in the last few years in using techniques based on the Jordan centrality measure and its variants to develop root finding algorithms. Given an unlabelled unrooted tree, one computes the Jordan centrality for each vertex in the tree and for a fixed budget $K$ outputs the optimal $K$ vertices (as measured by Jordan centrality). Under general conditions on the attachment function $f$, we derive necessary and sufficient bounds on the budget $K(\epsilon)$ in order to recover the root with probability at least $1-\epsilon$. For canonical examples such as linear preferential attachment and uniform attachment, these general results give matching upper and lower bounds for the budget. We also prove persistence of the optimal $K$ Jordan centers for any $K$, i.e. the existence of an almost surely finite random time $n^*$ such that for $n \geq n^*$ the identity of the $K$-optimal Jordan centers in $\{\mathcal{T}_f(n):n\geq n^*\}$ does not change, thus describing robustness properties of this measure. Key technical ingredients in the proofs of independent interest include sufficient conditions for the existence of exponential moments for limits of (appropriately normalized) continuous time branching processes within which the models $\{\mathcal{T}_f(n):n\geq n^*\}$ can be embedded, as well as rates of convergence results to these limits.


title: Lovelock black holes surrounded by dark fluid in power-Yang-Mills massive gravity

abstract: We consider a model where massive static spherically symmetric black hole, in the presence of power-Yang-Mills source, is surrounded by a dark fluid with non-linear equation of state. In this set up we construct a new class of magnetized Lovelock black hole solutions of the gravitational field equations. In particular, we work out the metric functions in both D-dimensional massive Einstein and massive Gauss-Bonnet gravities. We study thermodynamics of these black holes also and show that the mass and associated thermodynamic quantities like Hawking temperature and heat capacity depend on parameters of the dark fluid and the power-Yang-Mills magnetic source. We note that the entropy does not satisfy the area law in Gauss-Bonnet and higher order Lovelock black holes. Furthermore, phase transitions of black holes in each case are also discussed.


title: Phase space structure and escape time dynamics in a Van der Waals model for exothermic reactions

abstract: We study the phase space structures that control the transport in a classical Hamiltonian model for a chemical reaction. This model has been proposed to study the yield of products in an ultracold exothermic reaction. In the considered model, two elements determine the evolution of the system: a Van der Waals force and short-range force associated with the many-body interactions. In the previous work has been used small random periodic changes in the direction of the momentum to simulate the short-range many-body interactions. In the present work, random Gaussian bumps have been added to the Van der Waals potential energy simulate the short-range effects between the particles in the system. We compare both variants of the model and explain their differences similarities and differences from a phase space perspective. In order to visualize the structures that direct the dynamics in the phase space, we construct a natural Lagrangian descriptor for Hamiltonian systems based on the Maupertuis action $S_0$.


title: A symplectic look at the Fargues-Fontaine curve

abstract: This paper discusses homological mirror symmetry for the Fargues-Fontaine curve of equal characteristic.


title: Algebra over generalized rings

abstract: For a commutative ring $A$, we have the category of (bounded-below) chain complexes of $A$-modules $Ch_{+}(A\mymod)$, a closed symmetric monoidal category with a compatible stable Quillen model structure. The associated homotopy category is the derived category $\mathbbm{D}(A\mymod)$, where one inverts all the quasi-isomorphisms, and it has the good description as the chain complexes made up of projective $A$-module in each dimension, and chain maps taken up to chain homotopy. We give here the analogous theory for a (commutative) generalized ring in the sense of \cite{MR3605614}. We refer to the new concept as ``$\aset$''. For an ordinary commutative ring $A$, an $A$-set is just an $A$-module in the usual meaning, and our construction will be equivalent to $\mathbbm{D}(A\mymod)$. For the initial object of the category of generalized rings $\mathbb{F}$ ``the field with one element'', we obtain the category of symmetric spectra, and the associated stable homotopy category with its smash product (an $\mathbb{F}$-set is just a pointed set, i.e. a set $X$ with a distinguish element $O_X\in X$). Thus the analogous theories of stable homotopy and of chain complexes of modules over a commutative ring appear as two sides of the same coin, and moreover, they appear in a context where they interact (via the forgetful functor and its left adjoint - the base change functor). For the ``real integers'' $A=\Z_{\R}$, the $\Z_{\R}$-sets include the symmetric convex subsets of $\R$-vector spaces. We also give the global theory of the derived category of $\bigo_X$-sets, for a generalized scheme $X$, in a way that is based on the local projective model structure.


title: Multiple list colouring of $3$-choice critical graphs

abstract: A graph $G$ is called $3$-choice critical if $G$ is not $2$-choosable but any proper subgraph is $2$-choosable. A characterization of $3$-choice critical graphs was given by Voigt in [On list Colourings and Choosability of Graphs, Habilitationsschrift, Tu Ilmenau(1998)]. Voigt conjectured that if $G$ is a bipartite $3$-choice critical graph, then $G$ is $(4m, 2m)$-choosable for every integer $m$. This conjecture was disproved by Meng, Puleo and Zhu in [On (4, 2)-Choosable Graphs, Journal of Graph Theory 85(2):412-428(2017)]. They showed that if $G=\Theta_{r,s,t}$ where $r,s,t$ have the same parity and $\min\{r,s,t\} \ge 3$, or $G=\Theta_{2,2,2,2p}$ with $p \ge 2$, then $G$ is bipartite $3$-choice critical, but not $(4,2)$-choosable. On the other hand, all the other bipartite 3-choice critical graphs are $(4,2)$-choosable. This paper strengthens the result of Meng, Puleo and Zhu and shows that all the other bipartite $3$-choice critical graphs are $(4m,2m)$-choosable for every integer $m$.


title: A Real-Time Dispatching Strategy for Shared Automated Electric Vehicles with Performance Guarantees

abstract: Real-time vehicle dispatching operations in traditional car-sharing systems is an already computationally challenging scheduling problem. Electrification only exacerbates the computational difficulties as charge level constraints come into play. To overcome this complexity, we employ an online minimum drift plus penalty (MDPP) approach for SAEV systems that (i) does not require a priori knowledge of customer arrival rates to the different parts of the system (i.e. it is practical from a real-world deployment perspective), (ii) ensures the stability of customer waiting times, (iii) ensures that the deviation of dispatch costs from a desirable dispatch cost can be controlled, and (iv) has a computational time-complexity that allows for real-time implementation. Using an agent-based simulator developed for SAEV systems, we test the MDPP approach under two scenarios with real-world calibrated demand and charger distributions: 1) a low-demand scenario with long trips, and 2) a high-demand scenario with short trips. The comparisons with other algorithms under both scenarios show that the proposed online MDPP outperforms all other algorithms in terms of both reduced customer waiting times and vehicle dispatching costs.


title: Extensions of classical hypergeometric identities of Bailey and Whipple

abstract: We obtain extensions of classical hypergeometric identities of Bailey and Whipple that transform nearly-poised and very-well-poised series to Saalsch\"utzian series, Saalsch\"utzian series to Saalsch\"utzian series, and very-well-poised and nearly-poised series to very-well-poised series. We employ a method in which summations and transformations of lower-order series are used to obtain transformations of higher-order series. By taking limits, we also obtain extensions of two classical quadratic transformations of Whipple and Bailey. Furthermore, we show how a number of other well-known results regarding hypergeometric series follow as special cases of our results.


title: Shadow Removal by a Lightness-Guided Network with Training on Unpaired Data

abstract: Shadow removal can significantly improve the image visual quality and has many applications in computer vision. Deep learning methods based on CNNs have become the most effective approach for shadow removal by training on either paired data, where both the shadow and underlying shadow-free versions of an image are known, or unpaired data, where shadow and shadow-free training images are totally different with no correspondence. In practice, CNN training on unpaired data is more preferred given the easiness of training data collection. In this paper, we present a new Lightness-Guided Shadow Removal Network (LG-ShadowNet) for shadow removal by training on unpaired data. In this method, we first train a CNN module to compensate for the lightness and then train a second CNN module with the guidance of lightness information from the first CNN module for final shadow removal. We also introduce a loss function to further utilise the colour prior of existing data. Extensive experiments on widely used ISTD, adjusted ISTD and USR datasets demonstrate that the proposed method outperforms the state-of-the-art methods with training on unpaired data.


title: Analogical Image Translation for Fog Generation

abstract: Image-to-image translation is to map images from a given \emph{style} to another given \emph{style}. While exceptionally successful, current methods assume the availability of training images in both source and target domains, which does not always hold in practice. Inspired by humans' reasoning capability of analogy, we propose analogical image translation (AIT). Given images of two styles in the source domain: $\mathcal{A}$ and $\mathcal{A}^\prime$, along with images $\mathcal{B}$ of the first style in the target domain, learn a model to translate $\mathcal{B}$ to $\mathcal{B}^\prime$ in the target domain, such that $\mathcal{A}:\mathcal{A}^\prime ::\mathcal{B}:\mathcal{B}^\prime$. AIT is especially useful for translation scenarios in which training data of one style is hard to obtain but training data of the same two styles in another domain is available. For instance, in the case from normal conditions to extreme, rare conditions, obtaining real training images for the latter case is challenging but obtaining synthetic data for both cases is relatively easy. In this work, we are interested in adding adverse weather effects, more specifically fog effects, to images taken in clear weather. To circumvent the challenge of collecting real foggy images, AIT learns with synthetic clear-weather images, synthetic foggy images and real clear-weather images to add fog effects onto real clear-weather images without seeing any real foggy images during training. AIT achieves this zero-shot image translation capability by coupling a supervised training scheme in the synthetic domain, a cycle consistency strategy in the real domain, an adversarial training scheme between the two domains, and a novel network design. Experiments show the effectiveness of our method for zero-short image translation and its benefit for downstream tasks such as semantic foggy scene understanding.


title: Offline Handwritten Chinese Text Recognition with Convolutional Neural Networks

abstract: Deep learning based methods have been dominating the text recognition tasks in different and multilingual scenarios. The offline handwritten Chinese text recognition (HCTR) is one of the most challenging tasks because it involves thousands of characters, variant writing styles and complex data collection process. Recently, the recurrent-free architectures for text recognition appears to be competitive as its highly parallelism and comparable results. In this paper, we build the models using only the convolutional neural networks and use CTC as the loss function. To reduce the overfitting, we apply dropout after each max-pooling layer and with extreme high rate on the last one before the linear layer. The CASIA-HWDB database is selected to tune and evaluate the proposed models. With the existing text samples as templates, we randomly choose isolated character samples to synthesis more text samples for training. We finally achieve 6.81% character error rate (CER) on the ICDAR 2013 competition set, which is the best published result without language model correction.


title: Scattering from controlled defects in woodpile photonic crystals

abstract: Photonic crystals with a sufficiently high refractive index contrast display partial or full band gaps. However, imperfections in the metamaterial cause light scattering and extinction of the interfering propagating waves. Positive as well as negative defect volumes may contribute to this kind of optical perturbation. In this study, we fabricate and characterize three-dimensional woodpile photonic crystals, with a pseudo-bandgap for near-infrared optical wavelengths. By direct laser writing, we intentionally introduce random defects in the periodic structure. We show that we can model random defect scattering by considering the difference between the disordered and the regular structure. Our findings pave the way towards better control and understanding of the role of defects in photonic materials that will be crucial for their usability in potential applications.


title: Site percolation thresholds on triangular lattice with complex neighborhoods

abstract: We determine thresholds $p_c$ for random site percolation on a triangular lattice for neighbourhoods containing nearest (NN), next-nearest (2NN), next-next-nearest (3NN), next-next-next-nearest (4NN) and next-next-next-next-nearest (5NN) neighbours, and their combinations forming regular hexagons (3NN+2NN+NN, 5NN+4NN+NN, 5NN+4NN+3NN+2NN, 5NN+4NN+3NN+2NN+NN). We use a fast Monte Carlo algorithm, by Newman and Ziff [M. E. J. Newman and R. M. Ziff, Physical Review E 64, 016706 (2001)], for obtaining the dependence of the largest cluster size on occupation probability. The method is combined with a method, by Bastas et al. [N. Bastas, K. Kosmidis, P. Giazitzidis, and M. Maragakis, Physical Review E 90, 062101 (2014)], of estimating thresholds from low statistics data. The estimated values of percolation thresholds are $p_c(\text{4NN})=0.192410(43)$, $p_c(\text{3NN+2NN})=0.232008(38)$, $p_c(\text{5NN+4NN})=0.140286(5)$, $p_c(\text{3NN+2NN+NN})=0.215484(19)$, $p_c(\text{5NN+4NN+NN})=0.131792(58)$, $p_c(\text{5NN+4NN+3NN+2NN})=0.117579(41)$, $p_c(\text{5NN+4NN+3NN+2NN+NN})=0.115847(21)$. The method is tested on the standard case of site percolation on triangular lattice, where $p_c(\text{NN})=p_c(\text{2NN})=p_c(\text{3NN})=p_c(\text{5NN})=\frac{1}{2}$ is recovered with five digits accuracy $p_c(\text{NN})=0.500029(46)$ by averaging over one thousand lattice realisations only.


title: Massive particle pair production and oscillation in Friedman Universe: reheating energy and entropy, and cold dark matter

abstract: Suppose that the early Universe starts with a cosmological $\Lambda$-term originating from quantum spacetime at the Planck scale. Dark energy drives inflation and reheating by reducing its value for massive particle-antiparticle pairs production and oscillation, resulting in a holographic and massive pair plasma state. The back-and-forth reaction of dark energy and massive pairs slows inflation to its end and starts reheating by rapidly producing stable and unstable pairs. We introduce the Boltzmann-type rate equation describing the back-and-forth reaction. It forms a close set with Friedman equations and reheating equations for unstable pairs decay to relativistic particles. The numerical solutions show preheating, massive pairs dominated and genuine reheating episodes. We obtain the reheating temperature and entropy in terms of the tensor-to-scalar ratio $0 < r < 0.047$ consistently with observations. Stable massive pairs represent cold dark matter particles and weakly interact with dark energy. The resultant cold dark matter abundance $\Omega_c\sim 10^{-1}$ is about a constant in time.


title: Quantum Engineering of Superdark Excited States in Arrays of Atoms

abstract: We suggest a regular method of achieving an extremely long lifetime of a collective singly excited state in a generic small-size ensemble of N identical atoms. The decay rate \Gamma_N of such a `superdark' state can be as small as \Gamma_N \propto \Gamma(r/\lambda)^{2(N-1)} (\Gamma is the radiative decay rate of an individual atom, r and \lambda are the system size and the wavelength of the radiation, respectively), i.e., considerably smaller than in any of the systems suggested up to now. The method is based on a special fine tuning of the atomic Hamiltonian: namely, on a proper position-dependent adjustment of atomic transition frequencies. So chosen set of the control parameters is sufficient to ensure the minimum of the spontaneous decay rate of the engineered state in a generic ensemble of atoms (`qubits').


title: Application of Statistical Methods in Software Engineering: Theory and Practice

abstract: The experimental evaluation of the methods and concepts covered in software engineering has been increasingly valued. This value indicates the constant search for new forms of assessment and validation of the results obtained in Software Engineering research. Results are validated in studies through evaluations, which in turn become increasingly stringent. As an alternative to aid in the verification of the results, that is, whether they are positive or negative, we suggest the use of statistical methods. This article presents some of the main statistical techniques available, as well as their use in carrying out the implementation of data analysis in experimental studies in Software Engineering. This paper presents a practical approach proving statistical techniques through a decision tree, which was created in order to facilitate the understanding of the appropriate statistical method for each data analysis situation. Actual data from the software projects were employed to demonstrate the use of these statistical methods. Although it is not the aim of this work, basic experimentation and statistics concepts will be presented, as well as a concrete indication of the applicability of these techniques.


title: Thermal Regelation of a Single Particle

abstract: In this paper we study the process termed thermal regelation at single particle level. We experimentally extract the physical parameters that determine the interactions between ice and the silica particle for the regelation process. As the premelted layer around the particle plays a key role here, we use the extracted parameters to calculate the thickness of this thin premelted layer, at different given particle sizes.


title: Masking the general population might attenuate COVID-19 outbreaks

abstract: The effect of masking the general population on a COVID-19 epidemic is estimated by computer simulation using two separate state-of-the-art web-based softwares, one of them calibrated for the SARS-CoV-2 virus. The questions addressed are these: 1. Can mask use by the general population limit the spread of SARS-CoV-2 in a country? 2. What types of masks exist, and how elaborate must a mask be to be effective against COVID-19? 3. Does the mask have to be applied early in an epidemic? 4. A brief general discussion of masks and some possible future research questions regarding masks and SARS-CoV-2. Results are as follows: (1) The results indicate that any type of mask, even simple home-made ones, may be effective. Masks use seems to have an effect in lowering new patients even the protective effect of each mask (here dubbed "one-mask protection") is low. Strict adherence to mask use does not appear to be critical. However, increasing the one-mask protection to > 50% was found to be advantageous. Masks seemed able to reduce overflow of capacity, e.g. of intensive care. As the default parameters of the software included another intervention, it seems possible to combine mask and other interventions. (2) Masks do seem to reduce the number of new cases even if introduced at a late stage in an epidemic. However, early implementation helps reduce the cumulative and total number of cases. (3) The simulations suggest that it might be possible to eliminate a COVID-19 outbreak by widespread mask use during a limited period. The results from these simulations are encouraging, but do not necessarily represent the real-life situation, so it is suggested that clinical trials of masks are now carried out while continuously monitoring effects and side-effects.


title: Medium effects on $\Xi^-$ production in the nuclear ($K^-$, $K^+$) reaction

abstract: We study theoretically medium effects on $\Xi^-$ production in the ($K^-$, $K^+$) reaction, using the optimal Fermi-averaging procedure which describes the Fermi motion of a nucleon on the on-energy-shell $K^-p\to K^+\Xi^-$ reaction condition in nuclei. The result shows the strong energy and angular dependence of the in-medium $K^-p\to K^+\Xi^-$ cross section,which affects significantly the shape and magnitude of the production spectrum for $\Xi^-$ hypernuclear states in the ($K^-$, $K^+$) reaction on a nuclear target.The application to the $\Xi^-$ quasi-free production via the ($K^-$, $K^+$) reaction on a $^{12}$C target is also discussed in a Fermi gas model.


title: Deep Orthogonal Decompositions for Convective Nowcasting

abstract: Near-term prediction of the structured spatio-temporal processes driving our climate is of profound importance to the safety and well-being of millions, but the prounced nonlinear convection of these processes make a complete mechanistic description even of the short-term dynamics challenging. However, convective transport provides not only a principled physical description of the problem, but is also indicative of the transport in time of informative features which has lead to the recent successful development of ``physics free'' approaches to the now-casting problem. In this work we demonstrate that their remains an important role to be played by physically informed models, which can successfully leverage deep learning (DL) to project the process onto a lower dimensional space on which a minimal dynamical description holds. Our approach synthesises the feature extraction capabilities of DL with physically motivated dynamics to outperform existing model free approaches, as well as state of the art hybrid approaches, on complex real world datasets including sea surface temperature and precipitation.


title: Constraints of the formation and abundances of methyl carbamate, a glycine isomer, in hot corinos

abstract: Methyl carbamate CH$_3$OC(O)NH$_2$ is an isomer of glycine. Quantum chemical analyses show that methyl carbamate is more stable isomer than glycine. Because of this, there could be a higher chance for methyl carbamte to exist in the interstellar medium as compared to glycine. Despite immense searches, till now glycine has not been detected in the ISM, therefore it is worthwhile to search its isomer methyl carbamate. In this paper, we present the constraints of methyl carbamate formation under the interstellar conditions. Large complex organic molecules are favorably produced in hot-corino environments of low mass protostars. We for the first time carried out astrochemical modeling focusing on the formation of methyl carbamate in physical conditions similar to hot-corino objects. Consequently, we examined ALMA archival data for existing spectral line observations toward hot corinos NGC1333 IRAS 4A2 and IRAS 16293B. Within the common spectral range towards these sources, we found three features are possibly related to the spectral transitions of methyl carbamate and consequently estimate the upper limit of column densities. Results of chemical modeling are consistent with the observational upper limit of estimated column density/abundance toward the sources. This may hint the validation of the proposed formation mechanism. Future observations using telescope like ngVLA may confirm the presence of MC toward the hot corinos.


title: Thorn forking, weak normality, and theories with selectors

abstract: We discuss the role of weakly normal formulas in the theory of thorn forking, as part of a commentary on the paper "Thorn forking and stable forking" by Ealy and Onshuus (Rev. acad. colomb. cienc. exact. fis. nat. vol.40 no.157 Bogot\'a Oct./Dec. 2016). We also give a counterexample to Corollary 4.2 from that paper, and in the process discuss "theories with selectors".


title: Improving VQA and its Explanations \\ by Comparing Competing Explanations

abstract: Most recent state-of-the-art Visual Question Answering (VQA) systems are opaque black boxes that are only trained to fit the answer distribution given the question and visual content. As a result, these systems frequently take shortcuts, focusing on simple visual concepts or question priors. This phenomenon becomes more problematic as the questions become complex that requires more reasoning and commonsense knowledge. To address this issue, we present a novel framework that uses explanations for competing answers to help VQA systems select the correct answer. By training on human textual explanations, our framework builds better representations for the questions and visual content, and then reweights confidences in the answer candidates using either generated or retrieved explanations from the training set. We evaluate our framework on the VQA-X dataset, which has more difficult questions with human explanations, achieving new state-of-the-art results on both VQA and its explanations.


title: FDA3 : Federated Defense Against Adversarial Attacks for Cloud-Based IIoT Applications

abstract: Along with the proliferation of Artificial Intelligence (AI) and Internet of Things (IoT) techniques, various kinds of adversarial attacks are increasingly emerging to fool Deep Neural Networks (DNNs) used by Industrial IoT (IIoT) applications. Due to biased training data or vulnerable underlying models, imperceptible modifications on inputs made by adversarial attacks may result in devastating consequences. Although existing methods are promising in defending such malicious attacks, most of them can only deal with limited existing attack types, which makes the deployment of large-scale IIoT devices a great challenge. To address this problem, we present an effective federated defense approach named FDA3 that can aggregate defense knowledge against adversarial examples from different sources. Inspired by federated learning, our proposed cloud-based architecture enables the sharing of defense capabilities against different attacks among IIoT devices. Comprehensive experimental results show that the generated DNNs by our approach can not only resist more malicious attacks than existing attack-specific adversarial training methods, but also can prevent IIoT applications from new attacks.


title: Higher derivative scalar-tensor theory from the spatially covariant gravity: a linear algebraic analysis

abstract: We investigate the ghostfree scalar-tensor theory with a timelike scalar field, with derivatives of the scalar field up to the third order and with the Riemann tensor up to the quadratic order. We build two types of linear spaces. One is the set of linearly independent generally covariant scalar-tensor monomials, the other is the set of linearly independent spatially covariant gravity monomials. We argue that these two types of linear space are isomorphic to each other in the sense of gauge fixing/recovering procedures. We then identify the subspaces in the spatially covariant gravity, which are spanned by linearly independent monomials built of the extrinsic and intrinsic curvature, the lapse function as well as their spatial derivatives, up to the fourth order in the total number of derivatives. The vectors in these subspaces, i.e., spatially covariant polynomials, automatically propagate at most three degrees of freedom. As a result, their images under the gauge recovering mappings are automatically the subspaces of scalar-tensor theory that propagate up to three degrees of freedom as long as the scalar field is timelike. The mappings from the spaces of spatially covariant gravity to the spaces of scalar-tensor theory are encoded in the projection matrices, of which we also derived the expressions explicitly. Our formalism and results can be useful in deriving the generally covariant higher derivative scalar-tensor theory without ghost(s).


title: Non-Stationary Multi-layered Gaussian Priors for Bayesian Inversion

abstract: In this article, we study Bayesian inverse problems with multi-layered Gaussian priors. We first describe the conditionally Gaussian layers in terms of a system of stochastic partial differential equations. We build the computational inference method using a finite-dimensional Galerkin method. We show that the proposed approximation has a convergence-in-probability property to the solution of the original multi-layered model. We then carry out Bayesian inference using the preconditioned Crank--Nicolson algorithm which is modified to work with multi-layered Gaussian fields. We show via numerical experiments in signal deconvolution and computerized X-ray tomography problems that the proposed method can offer both smoothing and edge preservation at the same time.


title: Comparison of ablators for the polar direct drive exploding pusher platform

abstract: We examine the performance of pure boron, boron carbide, high density carbon, and boron nitride ablators in the polar direct drive exploding pusher (PDXP) platform. The platform uses the polar direct drive configuration at the National Ignition Facility to drive high ion temperatures in a room temperature capsule and has potential applications for plasma physics studies and as a neutron source. The higher tensile strength of these materials compared to plastic enables a thinner ablator to support higher gas pressures, which could help optimize its performance for plasma physics experiments, while ablators containing boron enable the possiblity of collecting addtional data to constrain models of the platform. Applying recently developed and experimentally validated equation of state models for the boron materials, we examine the performance of these materials as ablators in 2D simulations, with particular focus on changes to the ablator and gas areal density, as well as the predicted symmetry of the inherently 2D implosion.


title: Radial Lattice Quantization of 3D $\phi^4$ Field Theory

abstract: The quantum extension of classical finite elements, referred to as quantum finite elements ({\bf QFE})~\cite{Brower:2018szu,Brower:2016vsl}, is applied to the radial quantization of 3d $\phi^4$ theory on a simplicial lattice for the $\mathbb R \times \mathbb S^2$ manifold. Explicit counter terms to cancel the one- and two-loop ultraviolet defects are implemented to reach the quantum continuum theory. Using the Brower-Tamayo~\cite{Brower:1989mt} cluster Monte Carlo algorithm, numerical results support the QFE ansatz that the critical conformal field theory (CFT) is reached in the continuum with the full isometries of $\mathbb R \times \mathbb S^2$ restored. The Ricci curvature term, while technically irrelevant in the quantum theory, is shown to dramatically improve the convergence opening, the way for high precision Monte Carlo simulation to determine the CFT data: operator dimensions, trilinear OPE couplings and the central charge.


title: Deep Bayesian Quadrature Policy Optimization

abstract: We study the problem of obtaining accurate policy gradient estimates using a finite number of samples. Monte-Carlo methods have been the default choice for policy gradient estimation, despite suffering from high variance in the gradient estimates. On the other hand, more sample efficient alternatives like Bayesian quadrature methods have received little attention due to their high computational complexity. In this work, we propose deep Bayesian quadrature policy gradient (DBQPG), a computationally efficient high-dimensional generalization of Bayesian quadrature, for policy gradient estimation. We show that DBQPG can substitute Monte-Carlo estimation in policy gradient methods, and demonstrate its effectiveness on a set of continuous control benchmarks. In comparison to Monte-Carlo estimation, DBQPG provides (i) more accurate gradient estimates with a significantly lower variance, (ii) a consistent improvement in the sample complexity and average return for several deep policy gradient algorithms, and, (iii) the uncertainty in gradient estimation that can be incorporated to further improve the performance.


title: Improved Small Domain Estimation via Compromise Regression Weights

abstract: Shrinkage estimates of small domain parameters typically utilize a combination of a noisy "direct" estimate that only uses data from a specific small domain and a more stable regression estimate. When the regression model is misspecified, estimation performance for the noisier domains can suffer due to substantial shrinkage towards a poorly estimated regression surface. In this paper, we introduce a new class of robust, empirically-driven regression weights that target estimation of the small domain means under potential misspecification of the global regression model. Our regression weights are a convex combination of the model-based weights associated with the best linear unbiased predictor (BLUP) and those associated with the observed best predictor (OBP). The compromise parameter in this convex combination is found by minimizing a novel, unbiased estimate of the mean-squared prediction error for the small domain means, and we label the associated small domain estimates the "compromise best predictor", or CBP. Using a data-adaptive mixture for the regression weights enables the CBP to possess the robustness of the OBP while retaining the main advantages of the EBLUP whenever the regression model is correct. We demonstrate the use of the CBP in an application estimating gait speed in older adults.


title: Auger Recombination Coefficients in Type-I Mid-Infrared InGaAsSb Quantum Well Lasers

abstract: From a systematic study of the threshold current density as a function of temperature and hydrostatic pressure, in conjunction with theoretical analysis of the gain and threshold carrier density, we have determined the wavelength dependence of the Auger recombination coefficients in InGaAsSb/GaSb quantum well lasers emitting in the 1.7-3.2 $\mu$m wavelength range. From hydrostatic pressure measurements, the non-radiative component of threshold currents for individual lasers was determined continuously as a function of wavelength. The results are analysed to determine the Auger coefficients quantitatively. This procedure involves calculating the threshold carrier density based on device properties, optical losses, and estimated Auger contribution to the total threshold current density. We observe a minimum in the Auger rate around 2.1 $\mu$m. A strong increase with decreasing mid-infrared wavelength (< 2 $\mu$m) indicates the prominent role of inter-valence Auger transitions to the split-off hole band (CHSH process). Above 2 $\mu$m, the increase with wavelength is approximately exponential due to CHCC or CHLH Auger recombination, limiting long wavelength operation. The observed dependence is consistent with that derived by analysing literature values of lasing thresholds for type-I InGaAsSb quantum well diodes. Over the wavelength range considered, the Auger coefficient varies from a minimum of $\leq$ 1x10$ ^{16}$cm$^{4}$s$^{-1}$ at 2.1 $\mu$m to ~8x10$^{16}$cm$^{4}$s$^{-1}$ at 3.2 $\mu$m.

categories: physics.app-ph cond-mat.mes-hall cond-mat.mtrl-sci physics.optics


title: Valid model-free spatial prediction

abstract: Predicting the response at an unobserved location is a fundamental problem in spatial statistics. Given the difficulty in modeling spatial dependence, especially in non-stationary cases, model-based prediction intervals are at risk of misspecification bias that can negatively affect their validity. Here we present a new approach for model-free nonparametric spatial prediction based on the conformal prediction machinery. Our key observation is that spatial data can be treated as exactly or approximately exchangeable in a wide range of settings. In particular, under an infill asymptotic regime, we prove that the response values are, in a certain sense, locally approximately exchangeable for a broad class of spatial processes, and we develop a local spatial conformal prediction algorithm that yields valid prediction intervals without strong model assumptions like stationarity. Numerical examples with both real and simulated data confirm that the proposed conformal prediction intervals are valid and generally more efficient than existing model-based procedures for large datasets across a range of non-stationary and non-Gaussian settings.


title: Variational Autoencoding of PDE Inverse Problems

abstract: Specifying a governing physical model in the presence of missing physics and recovering its parameters are two intertwined and fundamental problems in science. Modern machine learning allows one to circumvent these, via emulators and surrogates, but in doing so disregards prior knowledge and physical laws that are especially important for small data regimes, interpretability, and decision making. In this work we fold the mechanistic model into a flexible data-driven surrogate to arrive at a physically structured decoder network. This provides accelerated inference for the Bayesian inverse problem, and can act as a drop-in regulariser that encodes a-priori physical information. We employ the variational form of the PDE problem and introduce stochastic local approximations as a form of model based data augmentation. We demonstrate both the accuracy and increased computational efficiency of the framework on real world settings and structured spatial processes.


title: The Wold-type decomposition for $m$-isometries

abstract: The aim of this paper is to study the Wold-type decomposition in the class of $m$-isometries. One of our main results establishes an equivalent condition for an analytic $m$-isometry to admit the Wold-type decomposition for $m\ge2$. In particular, we introduce the $k$-kernel condition which we use to characterize analytic $m$-isometric operators which are unitarily equivalent to unilateral operator valued weighted shifts for $m\ge2$. As a result, we also show that $m$-isometric composition operators on directed graphs with one circuit containing only one element are not unitarily equivalent to unilateral weighted shifts. We also provide a characterization of $m$-isometric unilateral operator valued weighted shifts with positive and commuting weights.


title: Investigating and Mitigating Degree-Related Biases in Graph Convolutional Networks

abstract: Graph Convolutional Networks (GCNs) show promising results for semi-supervised learning tasks on graphs, thus become favorable comparing with other approaches. Despite the remarkable success of GCNs, it is difficult to train GCNs with insufficient supervision. When labeled data are limited, the performance of GCNs becomes unsatisfying for low-degree nodes. While some prior work analyze successes and failures of GCNs on the entire model level, profiling GCNs on individual node level is still underexplored. In this paper, we analyze GCNs in regard to the node degree distribution. From empirical observation to theoretical proof, we confirm that GCNs are biased towards nodes with larger degrees with higher accuracy on them, even if high-degree nodes are underrepresented in most graphs. We further develop a novel Self-Supervised-Learning Degree-Specific GCN (SL-DSGC) that mitigate the degree-related biases of GCNs from model and data aspects. Firstly, we propose a degree-specific GCN layer that captures both discrepancies and similarities of nodes with different degrees, which reduces the inner model-aspect biases of GCNs caused by sharing the same parameters with all nodes. Secondly, we design a self-supervised-learning algorithm that creates pseudo labels with uncertainty scores on unlabeled nodes with a Bayesian neural network. Pseudo labels increase the chance of connecting to labeled neighbors for low-degree nodes, thus reducing the biases of GCNs from the data perspective. Uncertainty scores are further exploited to weight pseudo labels dynamically in the stochastic gradient descent for SL-DSGC. Experiments on three benchmark datasets show SL-DSGC not only outperforms state-of-the-art self-training/self-supervised-learning GCN methods, but also improves GCN accuracy dramatically for low-degree nodes.


title: Application of variational policy gradient to atomic-scale materials synthesis

abstract: Atomic-scale materials synthesis via layer deposition techniques present a unique opportunity to control material structures and yield systems that display unique functional properties that cannot be stabilized using traditional bulk synthetic routes. However, the deposition process itself presents a large, multidimensional space that is traditionally optimized via intuition and trial and error, slowing down progress. Here, we present an application of deep reinforcement learning to a simulated materials synthesis problem, utilizing the Stein variational policy gradient (SVPG) approach to train multiple agents to optimize a stochastic policy to yield desired functional properties. Our contributions are (1) A fully open source simulation environment for layered materials synthesis problems, utilizing a kinetic Monte-Carlo engine and implemented in the OpenAI Gym framework, (2) Extension of the Stein variational policy gradient approach to deal with both image and tabular input, and (3) Developing a parallel (synchronous) implementation of SVPG using Horovod, distributing multiple agents across GPUs and individual simulation environments on CPUs. We demonstrate the utility of this approach in optimizing for a material surface characteristic, surface roughness, and explore the strategies used by the agents as compared with a traditional actor-critic (A2C) baseline. Further, we find that SVPG stabilizes the training process over traditional A2C. Such trained agents can be useful to a variety of atomic-scale deposition techniques, including pulsed laser deposition and molecular beam epitaxy, if the implementation challenges are addressed.


title: Design of a $\beta$-Ga$_2$O$_3$ Schottky Barrier Diode With p-type III-Nitride Guard Ring for Enhanced Breakdown

abstract: This work presents the electrostatic analysis of a novel Ga$_2$O$_3$ vertical Schottky diode with three different guard ring configurations to reduce the peak electric field at the metal edges. Highly doped p-type GaN, p-type nonpolar AlGaN and polarization doped graded p-AlGaN are simulated and analyzed as the guard ring material, which forms a heterojunction with the Ga$_2$O$_3$ drift layer. Guard ring with non-polar graded p-AlGaN with a bandgap larger than Ga$_2$O$_3$ is found to show the best performance in terms of screening the electric field at the metal edges. The proposed guard ring configuration is also compared with a reported Ga$_2$O$_3$ Schottky diode with no guard ring and a structure with a high resistive Nitrogen-doped guard ring. The optimized design is predicted to have breakdown voltage as high as 5.3 kV and a specific on-resistance of 3.55 m$\Omega$-cm$^2$ which leads to an excellent power figure of merit of 7.91 GW/cm$^2$.


title: Expressive Power of Invariant and Equivariant Graph Neural Networks

abstract: Various classes of Graph Neural Networks (GNN) have been proposed and shown to be successful in a wide range of applications with graph structured data. In this paper, we propose a theoretical framework able to compare the expressive power of these GNN architectures. The current universality theorems only apply to intractable classes of GNNs. Here, we prove the first approximation guarantees for practical GNNs, paving the way for a better understanding of their generalization. Our theoretical results are proved for invariant GNNs computing a graph embedding (permutation of the nodes of the input graph does not affect the output) and equivariant GNNs computing an embedding of the nodes (permutation of the input permutes the output). We show that Folklore Graph Neural Networks (FGNN), which are tensor based GNNs augmented with matrix multiplication are the most expressive architectures proposed so far for a given tensor order. We illustrate our results on the Quadratic Assignment Problem (a NP-Hard combinatorial problem) by showing that FGNNs are able to learn how to solve the problem, leading to much better average performances than existing algorithms (based on spectral, SDP or other GNNs architectures). On a practical side, we also implement masked tensors to handle batches of graphs of varying sizes.


title: I can attend a meeting too! Towards a human-like telepresence avatar robot to attend meeting on your behalf

abstract: Telepresence robots are used in various forms in various use-cases that helps to avoid physical human presence at the scene of action. In this work, we focus on a telepresence robot that can be used to attend a meeting remotely with a group of people. Unlike a one-to-one meeting, participants in a group meeting can be located at a different part of the room, especially in an informal setup. As a result, all of them may not be at the viewing angle of the robot, a.k.a. the remote participant. In such a case, to provide a better meeting experience, the robot should localize the speaker and bring the speaker at the center of the viewing angle. Though sound source localization can easily be done using a microphone-array, bringing the speaker or set of speakers at the viewing angle is not a trivial task. First of all, the robot should react only to a human voice, but not to the random noises. Secondly, if there are multiple speakers, to whom the robot should face or should it rotate continuously with every new speaker? Lastly, most robotic platforms are resource-constrained and to achieve a real-time response, i.e., avoiding network delay, all the algorithms should be implemented within the robot itself. This article presents a study and implementation of an attention shifting scheme in a telepresence meeting scenario which best suits the needs and expectations of the collocated and remote attendees. We define a policy to decide when a robot should rotate and how much based on real-time speaker localization. Using user satisfaction study, we show the efficacy and usability of our system in the meeting scenario. Moreover, our system can be easily adapted to other scenarios where multiple people are located.


title: Selling sex: what determines rates and popularity? An analysis of 11,500 online profiles

abstract: Sex work, or the exchange of sexual services for money or goods, is ubiquitous across eras and cultures. However, the practice of selling sex is often hidden due to stigma and the varying legal status of sex work. Online platforms that sex workers use to advertise services have become an increasingly important means of studying a market that is largely hidden. Although prior literature has primarily shed light on sex work from a public health or policy perspective (focusing largely on female sex workers), there are few studies that empirically research patterns of service provision in online sex work. This study investigated the determinants of pricing and popularity in the market for commercial sexual services online by using data from the largest UK network of online sexual services, a platform that is the industry-standard for sex workers. While the size of these influences varies across genders, nationality, age and the services provided are shown to be primary drivers of rates and popularity in sex work.


title: Two Views of the Radius Gap and the Role of Light Curve Fitting

abstract: Recently, several groups have resolved a gap that bifurcates planets between the size of Earth and Neptune into two populations. The location and depth of this feature is an important signature of the physical processes that form and sculpt planets. In particular, planets residing in the radius gap are valuable probes of these processes as they may be undergoing the final stages of envelope loss. Here, we discuss two views of the radius gap by Fulton & Petigura (2018; F18) and Van Eylen et al. (2018; V18). In V18, the gap is wider and more devoid of planets. This is due, in part, to V18's more precise measurements of planet radius $R_p$. Thanks to Gaia, uncertainties in stellar radii Rstar are no longer the limiting uncertainties in determining $R_p$ for the majority of Kepler planets; instead, errors in $R_p/R_\star$ dominate. V18's analysis incorporated short-cadence photometry along with constraints on mean stelar density that enabled more accurate determinations of $R_p/R_\star$. In the F18 analysis, less accurate $R_p/R_\star$ blurs the boundary the radius gap. The differences in $R_p/R_\star$ are largest at high impact parameter ($b \gtrsim 0.8$) and often exceed 10%. This motivates excluding high-$b$ planets from demographic studies, but identifying such planets from long-cadence photometry alone is challenging. We show that transit duration can serve as an effective proxy, and we leverage this information to enhance the contrast between the super-Earth and sub-Neptune populations.


title: Social Distancing is Good for Points too!

abstract: The nearest-neighbor rule is a well-known classification technique that, given a training set P of labeled points, classifies any unlabeled query point with the label of its closest point in P. The nearest-neighbor condensation problem aims to reduce the training set without harming the accuracy of the nearest-neighbor rule. FCNN is the most popular algorithm for condensation. It is heuristic in nature, and theoretical results for it are scarce. In this paper, we settle the question of whether reasonable upper-bounds can be proven for the size of the subset selected by FCNN. First, we show that the algorithm can behave poorly when points are too close to each other, forcing it to select many more points than necessary. We then successfully modify the algorithm to avoid such cases, thus imposing that selected points should "keep some distance". This modification is sufficient to prove useful upper-bounds, along with approximation guarantees for the algorithm.


title: The maximum regularity property of the steady Stokes problem associated with a flow through a profile cascade

abstract: We deal with a steady Stokes-type problem, associated with a flow of a Newtonian incompressible fluid through a spatially periodic profile cascade. The used mathematical model is based on the reduction to one spatial period, represented by a bounded 2D domain Omega. The corresponding Stokes-type problem is formulated by means of the Stokes equation, equation of continuity and three types of boundary conditions: the conditions of periodicity on the curves Gamma0 and Gamma1, the Dirichlet boundary conditions Gamma-in and Gamma-p and an artificial "do nothing"-type boundary condition on Gamma-out. (See Fig. 1.) We explain on the level of weak solutions the sense in which the last condition is satisfied. We show that, although domain Omega is not smooth and different types of boundary conditions meet in the corners of Omega, the considered problem has a strong solution with the so called maximum regularity property.


title: Irregular conformal blocks, Painlev\'e III and the blow-up equations

abstract: We study the relation of irregular conformal blocks with the Painlev\'e III$_3$ equation. The functional representation for the quasiclassical irregular block is shown to be consistent with the BPZ equations of conformal field theory and the Hamilton-Jacobi approach to Painlev\'e III$_3$. It leads immediately to a limiting case of the blow-up equations for dual Nekrasov partition function of 4d pure supersymmetric gauge theory, which can be even treated as a defining system of equations for both $c=1$ and $c\to\infty$ conformal blocks. We extend this analysis to the domain of strong-coupling regime where original definition of conformal blocks and Nekrasov functions is not known and apply the results to spectral problem of the Matheiu equations. Finally, we propose a construction of irregular conformal blocks in the strong coupling region by quantization of Painlev\'e III$_3$ equation, and obtain in this way a general expression, reproducing $c=1$ and quasiclassical $c\to\infty$ results as its particular cases. We have also found explicit integral representations for $c=1$ and $c=-2$ irregular blocks at infinity for some special points.


title: Spin- and Voltage-dependent emission from Intra- and Intermolecular TADF OLEDs

abstract: Organic light emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) utilize molecular systems with a small energy splitting between singlet and triplet states. This can either be realized in intramolecular charge transfer states of molecules with near-orthogonal donor and acceptor moieties or in intermolecular exciplex states formed between a suitable combination of individual donor and acceptor materials. Here, we investigate 4,4'-(9H,9'H-[3,3'-bicarbazole]-9,9'-diyl)bis(3-(trifluoromethyl) benzonitrile) (pCNBCzoCF3), which shows intramolecular TADF but can also form exciplex states in combination with 4,4',4''-tris[phenyl(m-tolyl)amino]triphenylamine (m-MTDATA). Orange emitting exciplex-based OLEDs additionally generate a sky-blue emission from the intramolecular emitter with an intensity that can be voltage-controlled. We apply electroluminescence detected magnetic resonance (ELDMR) to study the thermally activated spin-dependent triplet to singlet up-conversion in operating devices. Thereby, we can investigate intermediate excited states involved in OLED operation and derive the corresponding activation energy for both, intra- and intermolecular based TADF. Furthermore, we give a lower estimate for the extent of the triplet wavefunction to be >1.2 nm. Photoluminescence detected magnetic resonance (PLDMR) reveals the population of molecular triplets in optically excited thin films. Overall, our findings allow us to draw a comprehensive picture of the spin-dependent emission from intra- and intermolecular TADF OLEDs.


title: Double Degenerate Bose-Fermi Mixture of Strontium and Lithium

abstract: We report on the attainment of a degenerate Fermi gas of $\rm^{6}Li$ in contact with a Bose-Einstein condensate (BEC) of $^{84}$Sr. A degeneracy of $T/T_F=0.33(3)$ is observed with $1.6\times10^5$ $^{6}$Li atoms in the two lowest energy hyperfine states together with an almost pure BEC of $3.1\times10^5$ $^{84}$Sr atoms. The elastic s-wave scattering length between $^6$Li and $^{84}$Sr is estimated to be $|a_{\rm^{6}Li-\rm^{84}Sr}|=(7.1_{-1.7}^{+2.6})a_0$ ($a_0$ being the Bohr radius) from measured interspecies thermalization rates in an optical dipole trap.


title: Physics-aware registration based auto-encoder for convection dominated PDEs

abstract: We design a physics-aware auto-encoder to specifically reduce the dimensionality of solutions arising from convection-dominated nonlinear physical systems. Although existing nonlinear manifold learning methods seem to be compelling tools to reduce the dimensionality of data characterized by a large Kolmogorov n-width, they typically lack a straightforward mapping from the latent space to the high-dimensional physical space. Moreover, the realized latent variables are often hard to interpret. Therefore, many of these methods are often dismissed in the reduced order modeling of dynamical systems governed by the partial differential equations (PDEs). Accordingly, we propose an auto-encoder type nonlinear dimensionality reduction algorithm. The unsupervised learning problem trains a diffeomorphic spatio-temporal grid, that registers the output sequence of the PDEs on a non-uniform parameter/time-varying grid, such that the Kolmogorov n-width of the mapped data on the learned grid is minimized. We demonstrate the efficacy and interpretability of our approach to separate convection/advection from diffusion/scaling on various manufactured and physical systems.

categories: math.DS cs.LG cs.NA math.NA physics.flu-dyn


title: Hydrodynamics of quantum corrections to the Coulomb interaction via the third rank tensor evolution equation: Application to the Langmuir waves and the spin-electron-acoustic waves

abstract: If we study the quantum effects in plasmas in terms of traditional hydrodynamics via the continuity and Euler equations we find the quantum Bohm potential and the force of spin-spin interaction. However, if we extend the set hydrodynamic equations beyond the 13-moments approximation, and include the third rank tensor evolution equation along with the pressure evolution equation, we obtain the quantum corrections to the Coulomb interaction. It is found in contrast with the fact that hydrodynamic equations for the higher rank tensors do not contain interaction in the classic plasmas studied in the selfconsistent (meanfield) approximation. Therefore, we present the quantum hydrodynamic model, where the quantum effects are studied beyond the quantum Bohm potential. Developed model is considered in two regimes: all electrons in plasmas are considered as the single fluid, and the separate spin evolution regime, where electrons with different spin projections are considered as two different fluids. To illustrate the fundamental meaning of found quantum effects we demonstrate their contribution in the spectrum of the Langmuir waves and the spin-electron-acoustic waves. It is worth to mention that the application of the pressure evolution equation ensures that the contribution of pressure in the Langmuir wave spectrum is proportional to $(3/5)v_{Fe}^{2}$, unlike $(1/3)v_{Fe}^{2}$ appearing from hydrodynamics based on the continuity and Euler equations, where $v_{Fe}$ is the Fermi velocity. Same correction corresponds on other plasmas phenomena like the speed of sound for spin-electron-acoustic waves. Moreover, it is found that novel quantum effects provide the novel wave solutions.


title: Learning Goals from Failure

abstract: We introduce a framework that predicts the goals behind observable human action in video. Motivated by evidence in developmental psychology, we leverage video of unintentional action to learn video representations of goals without direct supervision. Our approach models videos as contextual trajectories that represent both low-level motion and high-level action features. Experiments and visualizations show our trained model is able to predict the underlying goals in video of unintentional action. We also propose a method to "automatically correct" unintentional action by leveraging gradient signals of our model to adjust latent trajectories. Although the model is trained with minimal supervision, it is competitive with or outperforms baselines trained on large (supervised) datasets of successfully executed goals, showing that observing unintentional action is crucial to learning about goals in video. Project page: https://aha.cs.columbia.edu/


title: From the open Heisenberg model to the Landau-Lifshitz equation

abstract: Magnetic systems can be described by the classical Landau-Lifshitz (LL) equation or the fully quantum open Heisenberg model. Using the Lindblad master equation and the mean-field approximation, we demonstrate that the open Heisenberg model is reduced to a generalized LL equation. The open dynamic is modeled using spin-boson interactions with a common bosonic reservoir at thermal equilibrium. By tracing out the bosonic degrees of freedom, we obtain two different decoherence mechanisms: on-site dissipation and an effective spin-spin interaction mediated by bosons. Using our approach, we perform hysteresis calculations, closely connected with the Stoner-Wohlfarth theory. We compare the exact numerical master equation and the mean-field model, revealing the role of correlations originated by non-local interactions. Our work opens new horizons for the study of the LL dynamics from an open quantum formalism.


title: Mechanical Evaluation Of Selective Laser Melted Ni-Rich Niti: Compression, Tension, And Torsion

abstract: Two unique behaviors of superelasticity and shape memory effect have made shape memory alloy such as NiTi, an interesting alloy for different applications. Recently, additive manufacturing (AM) as a powerful tool for fabricating NiTi has become of interest to make complex geometries. Selective laser melting (SLM) is an AM method that not only provides flexibility to make complex 3D shapes but also it can be possible to tailor the thermomechanical properties of the parts just by changing the process parameters. The non-homogeneous microstructure of as-fabricated parts as well as asymmetric mechanical behavior of NiTi, make it important to study the properties of SLM NiTi parts under different loading condition. In this study, Ni50.8Ti (at. %) powder was utilized to fabricate cube, dog-bone, and tube by SLM technique. The transformation temperatures (TTs) of samples were measured by the differential scanning calorimetry (DSC) method and the variation of TTs was discussed. Three coupons were tested mechanically under compression, tension, and torsion. In-situ digital image correlation (DIC) was employed to measure and monitor the strain field of samples during the mechanical tests. The strain distribution showed localized strain for all three samples. The equivalent stress/strain was calculated to compare the result of compressive, tensile, and torsional responses and the significant asymmetric behavior was shown and discussed.


title: Imaging the expanding knotty structure in the close environment of the LBV star $\eta$ Carinae

abstract: $\eta$~Car is one of the most massive stars in the Galaxy. It underwent a massive eruption in the 19th century, which produced the impressive bipolar Homunculus nebula now surrounding it. The central star is an eccentric binary with a period of 5.54\,years. Although the companion has not been detected directly, it causes time-variable ionization and colliding-wind X-ray emission. By characterizing the complex structure and kinematics of the ejecta close to the star, we aim to constrain past and present mass loss of $\eta$~Car. $\eta$~Car is observed with the extreme adaptive optics instrument SPHERE at the Very Large Telescope, using its polarimetric mode in the optical with the ZIMPOL camera. A spatial resolution of 20\,mas was achieved, i.e. very close to the presumed 13 mas apastron separation of the companion star. We detect new structures within the inner arcsecond to the star (2\,300\,au at a 2.3\,kpc distance). We can relate these structures to the eruption near 1890 by tracking their proper motions derived from our new images and historical images over a 30\,years time span. Besides, we find a fan-shaped structure in the inner 200~au to the star in the H$\alpha$ line, that could potentially be associated with the wind collision zone of the two stars.


title: Non-vanishing for cubic $L$--functions

abstract: We prove that there is a positive proportion of $L$-functions associated to cubic characters over $\mathbb{F}_q[T]$ that do not vanish at the critical point $s=1/2$. This is achieved by computing the first mollified moment using techniques previously developed by the authors in their work on the first moment of cubic $L$-functions, and by obtaining a sharp upper bound for the second mollified moment, building on work of Lester-Radziwill, which in turn develops further ideas from the work of Soundararajan, Harper, and Radziwill-Soundararajan. We work in the non-Kummer setting when $q \equiv 2\pmod{3}$, but our results could be translated into the Kummer setting when $q\equiv 1\pmod{3}$ as well as into the number field case (assuming the Generalized Riemann Hypothesis). Our positive proportion of non-vanishing is explicit, but extremely small, due to the fact that the implied constant in the upper bound for the mollified second moment is very large.


title: A Study of One-Parameter Regularization Methods for Mathematical Programs with Vanishing Constraints

abstract: Mathematical programs with vanishing constraints (MPVCs) are a class of nonlinear optimization problems with applications to various engineering problems such as truss topology design and robot motion planning. MPVCs are difficult problems from both a theoretical and numerical perspective: the combinatorial nature of the vanishing constraints often prevents standard constraint qualifications and optimality conditions from being attained; moreover, the feasible set is inherently nonconvex, and often has no interior around points of interest. In this paper, we therefore study and compare four regularization methods for the numerical solution of MPVCS. Each method depends on a single regularization parameter, which is used to embed the original MPVC into a sequence of standard nonlinear programs. Convergence results for these methods based on both exact and approximate stationary of the subproblems are established under weak assumptions. The improved regularity of the subproblems is studied by providing sufficient conditions for the existence of KKT multipliers. Numerical experiments, based on applications in truss topology design and an optimal control problem from aerothermodynamics, complement the theoretical analysis and comparison of the regularization methods. The computational results highlight the benefit of using regularization over applying a standard solver directly, and they allow us to identify two promising regularization schemes.


title: Inertial manifolds via spatial averaging revisited

abstract: The paper gives a comprehensive study of inertial manifolds for semilinear parabolic equations and their smoothness using the spatial averaging method suggested by G. Sell and J. Mallet-Paret. We present a universal approach which covers the most part of known results obtained via this method as well as gives a number of new ones. Among our applications are reaction-diffusion equations, various types of generalized Cahn-Hilliard equations, including fractional and 6th order Cahn-Hilliard equations and several classes of modified Navier-Stokes equations including the Leray-$\alpha$ regularization, hyperviscous regularization and their combinations. All of the results are obtained in 3D case with periodic boundary conditions.


title: Minimizing The Maximum Distance Traveled To Form Patterns With Systems of Mobile Robots

abstract: In the pattern formation problem, robots in a system must self-coordinate to form a given pattern, regardless of translation, rotation, uniform-scaling, and/or reflection. In other words, a valid final configuration of the system is a formation that is \textit{similar} to the desired pattern. While there has been no shortage of research in the pattern formation problem under a variety of assumptions, models, and contexts, we consider the additional constraint that the maximum distance traveled among all robots in the system is minimum. Existing work in pattern formation and closely related problems are typically application-specific or not concerned with optimality (but rather feasibility). We show the necessary conditions any optimal solution must satisfy and present a solution for systems of three robots. Our work also led to an interesting result that has applications beyond pattern formation. Namely, a metric for comparing two triangles where a distance of $0$ indicates the triangles are similar, and $1$ indicates they are \emph{fully dissimilar}.


title: Physical pendulum model: Fractional differential equation and memory effects

abstract: A detailed analysis of three pendular motion models is presented. Inertial effects, self-oscillation, and memory, together with non-constant moment of inertia, hysteresis and negative damping are shown to be required for the comprehensive description of the free pendulum oscillatory regime. The effects of very high initial amplitudes, friction in the roller bearing axle, drag, and pendulum geometry are also analysed and discussed. The model that consists of a fractional differential equation provides both the best explanation of, and the best fits to, experimental high resolution and long-time data gathered from standard action-camera videos.