title: A concept for the extraction of the most refractory elements at CERN-ISOLDE as carbonyl complex ions

abstract: We introduce a novel thick-target concept tailored to the extraction of refractory 4d and 5d transition metal radionuclides of molybdenum, technetium, ruthenium and tungsten for radioactive ion beam production. Despite the more than 60-year old history of thick-target ISOL mass-separation facilities like ISOLDE, the extraction of the most refractory elements as radioactive ion beam has so far not been successful. In ordinary thick ISOL targets, their radioisotopes produced in the target are stopped within the condensed target material and have to diffuse through a solid material. Here, we present a concept which overcomes limitations associated with this method. We exploit the recoil momentum of nuclear reaction products for their release from the solid target material. They are thermalized in a carbon monoxide-containing atmosphere, in which volatile carbonyl complexes form readily at ambient temperature and pressure. This compound serves as volatile carrier for transport to the ion source. Excess carbon monoxide is removed by cryogenic gas separation to enable low pressures in the source region, in which the species are ionized and hence made available for radioactive ion beam formation. The setup is operated in batch mode, with the aim to extract isotopes having half-lives of at least several seconds. We report parameter studies of the key processes of the method, which validate this concept and which define the parameters for the setup. This would allow for the first time the extraction of radioactive molybdenum, tungsten and several other transition metals at thick-target ISOL facilities.


title: Stochastic evolution equations driven by cylindrical stable noise

abstract: We prove existence and uniqueness of a mild solution of a stochastic evolution equation driven by a standard $\alpha$-stable cylindrical L\'evy process defined on a Hilbert space for $\alpha \in (1,2)$. The coefficients are assumed to map between certain domains of fractional powers of the generator present in the equation. The solution is constructed as a weak limit of the Picard iteration using tightness arguments. Existence of strong solution is obtained by a general version of the Yamada--Watanabe theorem.


title: Non-Hermitian quantum gases: a platform for imaginary time crystals

abstract: One of the most important applications of quantum mechanics is the thermodynamic description of quantum gases. Despite the fundamental importance of this topic, a comprehensive description of the thermodynamic properties of non-Hermitian quantum gases is still lacking. Here, we investigate the properties of bosonic and fermionic non-Hermitian systems at finite temperatures. We show that non-Hermitian systems exihibit oscillations both in temperature and imaginary time. As such, they can be a possible platform to realize an imaginary time crystal (iTC) phase. The Hatano-Nelson model is identified as a simple lattice model to reveal this effect. In addition, we show that the conditions for the iTC to be manifest are the same as the conditions for the presence of disorder points, where the correlation functions show oscillating behavior. This analysis makes clear that our realization of iTC is effectively a way to filter one specific Matsubara mode. In this realization, the Matsubara frequency, that enters as a mathematical tool to compute correlation functions for finite temperatures, can be measured experimentally.


title: Bayesian I-optimal designs for choice experiments with mixtures

abstract: Discrete choice experiments are frequently used to quantify consumer preferences by having respondents choose between different alternatives. Choice experiments involving mixtures of ingredients have been largely overlooked in the literature, even though many products and services can be described as mixtures of ingredients. As a consequence, little research has been done on the optimal design of choice experiments involving mixtures. The only existing research has focused on D-optimal designs, which means that an estimation-based approach was adopted. However, in experiments with mixtures, it is crucial to obtain models that yield precise predictions for any combination of ingredient proportions. This is because the goal of mixture experiments generally is to find the mixture that optimizes the respondents' utility. As a result, the I-optimality criterion is more suitable for designing choice experiments with mixtures than the D-optimality criterion because the I-optimality criterion focuses on getting precise predictions with the estimated statistical model. In this paper, we study Bayesian I-optimal designs, compare them with their Bayesian D-optimal counterparts, and show that the former designs perform substantially better than the latter in terms of the variance of the predicted utility.


title: Follow-up of 27 radio-quiet gamma-ray pulsars at 110-190 MHz using the international LOFAR station FR606

abstract: The Fermi Large Area Telescope has detected over 260 gamma-ray pulsars. About one quarter of these are labeled as radio-quiet. In the population of nonrecycled gamma-ray pulsars, the fraction of radio-quiet pulsars is higher, about one half. Most radio observations of gamma-ray pulsars have been performed at frequencies between 300 MHz and 2 GHz. However, pulsar radio fluxes increase rapidly with decreasing frequency, and their radio beams often broaden at low frequencies. As a consequence, some of these pulsars might be detectable at low radio frequencies even when no radio flux is detected above 300 MHz. Our aim is to test this hypothesis with low-frequency radio observations. We have observed 27 Fermi-discovered gamma-ray pulsars with the international LOw Frequency ARray (LOFAR) station FR606 in single-station mode. We used the LOFAR high band antenna (HBA) band (110-190 MHz). On average, we use 9 h of observation per target after the removal of affected datasets, resulting in a sensitivity for pulse-averaged flux on the order of 1-10 mJy. We do not detect radio pulsations from any of the 27 sources, and we establish stringent upper limits on their low-frequency radio fluxes. These nondetections are compatible with the upper limits derived from radio observations at other frequencies. We also determine the pulsars' geometry from the gamma-ray profiles to see for which pulsars the low-frequency radio beam is expected to cross Earth. This set of observations provides the most constraining upper limits on the flux density at 150 MHz for 27 radio-quiet gamma-ray pulsars. In spite of the beam-widening expected at low radio frequencies, most of our nondetections can be explained by an unfavorable viewing geometry; for the remaining observations, especially those of pulsars detected at higher frequencies, the nondetection is compatible with insufficient sensitivity.


title: Ellipsotopes: Combining Ellipsoids and Zonotopes for Reachability Analysis and Fault Detection

abstract: Ellipsoids are a common representation for reachability analysis, because they can be transformed efficiently under affine maps, and allow conservative approximation of Minkowski sums, which let one incorporate uncertainty and linearization error in a dynamical system by expanding the size of the reachable set. Zonotopes, a type of symmetric, convex polytope, are similarly frequently used due to efficient numerical implementation of affine maps and exact Minkowski sums. Both of these representations also enable efficient, convex collision detection for fault detection or formal verification tasks, wherein one checks if the reachable set of a system collides (i.e., intersects) with an unsafe set. However, both representations often result in conservative representations for reachable sets of arbitrary systems, and neither is closed under intersection. Recently, representations such as constrained zonotopes and constrained polynomial zonotopes have been shown to overcome some of these conservativeness challenges, and are closed under intersection. However, constrained zonotopes can not represent shapes with smooth boundaries such as ellipsoids, and constrained polynomial zonotopes can require solving a non-convex program for collision checking or fault detection. This paper introduces ellipsotopes, a set representation that is closed under affine maps, Minkowski sums, and intersections. Ellipsotopes combine the advantages of ellipsoids and zonotopes while ensuring convex collision checking. The utility of this representation is demonstrated on several examples.


title: Local Fourier Analysis of P-Multigrid for High-Order Finite Element Operators

abstract: Multigrid methods are popular for solving linear systems derived from discretizing PDEs. Local Fourier Analysis (LFA) is a technique for investigating and tuning multigrid methods. P-multigrid is popular for high-order or spectral finite element methods, especially on unstructured meshes. In this paper, we introduce LFAToolkit.jl, a new Julia package for LFA of high-order finite element methods. LFAToolkit.jl analyzes preconditioning techniques for arbitrary systems of second order PDEs and supports mixed finite element methods. Specifically, we develop LFA of p-multigrid with arbitrary second-order PDEs using high-order finite element discretizations and examine the performance of Jacobi and Chebyshev smoothing for two-grid schemes with aggressive p-coarsening. A natural extension of this LFA framework is the analysis of h-multigrid for finite element discretizations or finite difference discretizations that can be represented in the language of finite elements. With this extension, we can replicate previous work on the LFA of h-multigrid for arbitrary order discretizations using a convenient and extensible abstraction. Examples in one, two, and three dimensions are presented to validate our LFA of p-multigrid for the Laplacian and linear elasticity.


title: Cosmological Strangeness Abundance

abstract: We investigate the strange particle composition of the early Universe in the hadron epoch $T_h\approx 150\ge T\ge 10$\,MeV. We study strangeness yield in thermal and chemical equilibrium constrained by prescribed entropy per baryon in a charge neutral and strangeness neutral $\langle s-\bar s\rangle$ Universe. Turning to kinetic processes in a Hubble expanding Universe, we determine conditions at which individual strangeness producing reactions fall out of detailed balance between decay and back-reaction strangeness production rates in presence of decreasing temperature $T$; we allow for weak, electromagnetic, and strong interaction processes. The weak interaction $\mu^\pm+\nu_{\mu}\rightarrow K^\pm$ freezeout is at $T_f^{K^\pm}=33.8\,\mathrm{MeV}$; the electromagnetic process $l^-+l^+\rightarrow\phi$ freezeout is at $T_f^\phi=23\sim25\,\mathrm{MeV}$; and the hadronic reaction $\pi+\pi\rightarrow K$ freezeout is at $T_f^K=19.8\,\mathrm{MeV}$.


title: Optimally Hiding Object Sizes with Constrained Padding

abstract: Among the most challenging traffic-analysis attacks to confound are those leveraging the sizes of objects downloaded over the network. In this paper we systematically analyze this problem under realistic constraints regarding the padding overhead that the object store is willing to incur. We give algorithms to compute privacy-optimal padding schemes -- specifically that minimize the network observer's information gain from a downloaded object's padded size -- in several scenarios of interest: per-object padding, in which the object store responds to each request for an object with the same padded copy; per-request padding, in which the object store pads an object anew each time it serves that object; and a scenario unlike the previous ones in that the object store is unable to leverage a known distribution over the object queries. We provide constructions for privacy-optimal padding in each case, compare them to recent contenders in the research literature, and evaluate their performance on practical datasets.


title: Sub-nK thermometry of an interacting $d$-dimensional homogeneous Bose gas

abstract: We propose experimentally feasible means for non-destructive thermometry of homogeneous Bose Einstein condensates in different spatial dimensions ($d\in\{1,2,3\}$). Our impurity based protocol suggests that the fundamental error bound on thermometry at the sub nano Kelvin domain depends highly on the dimension, in that the higher the dimension the better the precision. Furthermore, sub-optimal thermometry of the condensates by using measurements that are experimentally feasible is explored. We specifically focus on measuring position and momentum of the impurity that belong to the family of Gaussian measurements. We show that, generally, experimentally feasible measurements are far from optimal, except in 1D, where position measurements are indeed optimal. This makes realistic experiments perform very well at few nano Kelvin temperatures for all dimensions, and at sub nano Kelvin temperatures in the one dimensional scenario. These results take a significant step towards experimental realisation of probe-based quantum thermometry of Bose Einstein condensates, as it deals with them in one, two and three dimensions and uses feasible measurements applicable in current experimental setups.


title: Finite Horizon Privacy of Stochastic Dynamical Systems: A Synthesis Framework for Dependent Gaussian Mechanisms

abstract: We address the problem of synthesizing distorting mechanisms that maximize privacy of stochastic dynamical systems. Information about the system state is obtained through sensor measurements. This data is transmitted to a remote station through an unsecured/public communication network. We aim to keep part of the system state private (a private output); however, because the network is unsecured, adversaries might access sensor data and input signals, which can be used to estimate private outputs. To prevent an accurate estimation, we pass sensor data and input signals through a distorting (privacy-preserving) mechanism before transmission, and send the distorted data to the trusted user. These mechanisms consist of a coordinate transformation and additive dependent Gaussian vectors. We formulate the synthesis of the distorting mechanisms as a convex program, where we minimize the mutual information (our privacy metric) between an arbitrarily large sequence of private outputs and the disclosed distorted data for desired distortion levels -- how different actual and distorted data are allowed to be.


title: Localisable Monads

abstract: Monads govern computational side-effects in programming semantics. They can be combined in a ''bottom-up'' way to handle several instances of such effects. Indexed monads and graded monads do this in a modular way. Here, instead, we equip monads with fine-grained structure in a ''top-down'' way, using techniques from tensor topology. This provides an intrinsic theory of local computational effects without needing to know how constituent effects interact beforehand. Specifically, any monoidal category decomposes as a sheaf of local categories over a base space. We identify a notion of localisable monads which characterises when a monad decomposes as a sheaf of monads. Equivalently, localisable monads are formal monads in an appropriate presheaf 2-category, whose algebras we characterise. Three extended examples demonstrate how localisable monads can interpret the base space as locations in a computer memory, as sites in a network of interacting agents acting concurrently, and as time in stochastic processes.


title: $t$-channel singularities in cosmology and particle physics

abstract: The $t$-channel singularity of a cross-section for a binary $2 \to 2$ scattering process occurs when a particle exchanged in the $t$ channel is kinematically allowed to be on its mass-shell, that is when the process can be viewed as a sequence of two physical subprocesses: a two-body decay $1\to 2$ and an inverse decay $2 \to 1$. We derive conditions for the singularity to be present in a binary process. A class of divergent cross sections for Standard Model processes has been determined and illustrated by the weak analog of the Compton scattering $Ze^{-}\to Ze^{-}$. After critically reviewing regularization methods proposed in literature, we discuss singular processes which occur in a medium composed of particles in thermal equilibrium. The medium is shown to regulate the singularities naturally as particles acquire a finite width. We demonstrate a possible cosmological application by calculating thermally averaged cross section for an elastic scattering within a simple scalar model. The cross section, which is divergent in vacuum, becomes finite when the process occurs in a thermal bath due to the imaginary self-energy of the $t$-channel mediator computed within the Keldysh-Schwinger formalism.


title: Factor Representation and Decision Making in Stock Markets Using Deep Reinforcement Learning

abstract: Deep Reinforcement learning is a branch of unsupervised learning in which an agent learns to act based on environment state in order to maximize its total reward. Deep reinforcement learning provides good opportunity to model the complexity of portfolio choice in high-dimensional and data-driven environment by leveraging the powerful representation of deep neural networks. In this paper, we build a portfolio management system using direct deep reinforcement learning to make optimal portfolio choice periodically among S\&P500 underlying stocks by learning a good factor representation (as input). The result shows that an effective learning of market conditions and optimal portfolio allocations can significantly outperform the average market.


title: Biphoton phase-space correlations from Gouy-phase measurements using double slits

abstract: Quantum correlations encoded in photonic Laguerre-Gaussian modes were shown to be related to the Gouy phase shifts (D. Kawase et al., Phys. Rev. Lett. 101, 050501 (2008)) allowing for a non-destructive manipulation of photonic quantum states. In this work we exploit the relation between phase space correlations of biphotons produced by spontaneously parametric down conversion (SPDC) as encoded in the logarithmic negativity (LN) and the Gouy phase as they are diffracted through an asymmetrical double slit setup. Using an analytical approach based on a double-gaussian approximation for type-I SPDC biphotons, we show that measurements of Gouy phase differences provide information on their phase space entanglement variation, governed by the physical parameters of the experiment and expressed by the LN via covariance matrix elements.


title: Calibrating the Nelson-Siegel-Svensson Model by Genetic Algorithm

abstract: Accurately fitting the term structure of interest rates is critical to central banks and other market participants. The Nelson-Siegel and Nelson-Siegel-Svensson models are probably the best-known models for this purpose due to their intuitive appeal and simple representation. However, this simplicity comes at a price. The difficulty in calibrating these models is twofold. Firstly, the objective function being minimized during the calibration procedure is nonlinear and has multiple local optima. Secondly, there is strong co-dependence among the model parameters. As a result, their estimated values behave erratically over time. To avoid these problems, we apply a heuristic optimization method, specifically the Genetic Algorithm approach, and show that it is able to construct reliable interest rate curves and stable model parameters over time, regardless of the shape of the curves.


title: The solution space structure of planted constraint satisfaction problems with growing domains

abstract: Planting a solution into the random RB model, which is a prototype of random constraint satisfaction problem (CSP) with growing domains, can generate very hard satisfiable CSP benchmarks. We study the solution space structure of the planted RB model. With constraint density growing, we find that this model goes through four phase transitions. In the replica symmetric phase, what we call the independent phase transition occurs, after which the planted cluster (cluster containing the planted solution) is separated from the giant cluster. Then the solutions except that in the planted cluster go through the same clustering phase transition and the same satisfiability phase transition as the random RB model. The planted cluster goes through the isolated phase transition, after which the planted cluster contains only one solution. This phase diagram provides strong evidence that this model can generate very hard satisfiable CSP benchmarks. For over constraint instances (where the constraint density is very large), we find that the configuration space has only a single energy valley, which makes the instances tractable. Experiments using Belief Propagation confirm the locations of the clustering, satisfiability (by configurations outside the planted cluster), and isolated phase transition points.


title: Spectral properties of substitutions on compact alphabets

abstract: We consider substitutions on compact alphabets and provide sufficient conditions for the diffraction to be pure point, absolutely continuous and singular continuous. This allows one to construct examples for which the Koopman operator on the associated function space has specific spectral components. For abelian bijective substitutions, we provide a dichotomy result regarding the spectral type of the diffraction. We also provide the first example of a substitution that has countably infinite Lebesgue spectral components and countably infinite singular continuous components. Lastly, we give a non-constant length substitution on a countably infinite alphabet that gives rise to substitutive Delone sets of infinite type. This extends the spectral theory of substitutions on finite alphabets and Delone sets of finite type with inflation symmetry.


title: HTTP2vec: Embedding of HTTP Requests for Detection of Anomalous Traffic

abstract: Hypertext transfer protocol (HTTP) is one of the most widely used protocols on the Internet. As a consequence, most attacks (i.e., SQL injection, XSS) use HTTP as the transport mechanism. Therefore, it is crucial to develop an intelligent solution that would allow to effectively detect and filter out anomalies in HTTP traffic. Currently, most of the anomaly detection systems are either rule-based or trained using manually selected features. We propose utilizing modern unsupervised language representation model for embedding HTTP requests and then using it to classify anomalies in the traffic. The solution is motivated by methods used in Natural Language Processing (NLP) such as Doc2Vec which could potentially capture the true understanding of HTTP messages, and therefore improve the efficiency of Intrusion Detection System. In our work, we not only aim at generating a suitable embedding space, but also at the interpretability of the proposed model. We decided to use the current state-of-the-art RoBERTa, which, as far as we know, has never been used in a similar problem. To verify how the solution would work in real word conditions, we train the model using only legitimate traffic. We also try to explain the results based on clusters that occur in the vectorized requests space and a simple logistic regression classifier. We compared our approach with the similar, previously proposed methods. We evaluate the feasibility of our method on three different datasets: CSIC2010, CSE-CIC-IDS2018 and one that we prepared ourselves. The results we show are comparable to others or better, and most importantly - interpretable.


title: Q-Pain: A Question Answering Dataset to Measure Social Bias in Pain Management

abstract: Recent advances in Natural Language Processing (NLP), and specifically automated Question Answering (QA) systems, have demonstrated both impressive linguistic fluency and a pernicious tendency to reflect social biases. In this study, we introduce Q-Pain, a dataset for assessing bias in medical QA in the context of pain management, one of the most challenging forms of clinical decision-making. Along with the dataset, we propose a new, rigorous framework, including a sample experimental design, to measure the potential biases present when making treatment decisions. We demonstrate its use by assessing two reference Question-Answering systems, GPT-2 and GPT-3, and find statistically significant differences in treatment between intersectional race-gender subgroups, thus reaffirming the risks posed by AI in medical settings, and the need for datasets like ours to ensure safety before medical AI applications are deployed.


title: Identifiability of species network topologies from genomic sequences using the logDet distance

abstract: Inference of network-like evolutionary relationships between species from genomic data must address the interwoven signals from both gene flow and incomplete lineage sorting. The heavy computational demands of standard approaches to this problem severely limit the size of datasets that may be analyzed, in both the number of species and the number of genetic loci. Here we provide a theoretical pointer to more efficient methods, by showing that logDet distances computed from genomic-scale sequences retain sufficient information to recover network relationships in the level-1 ultrametric case. This result is obtained under the Network Multispecies Coalescent model combined with a mixture of General Time-Reversible sequence evolution models across individual gene trees, but does not depend on partitioning sequences by genes. Thus under standard stochastic models statistically justifiable inference of network relationships from sequences can be accomplished without consideration of individual genes or gene trees.


title: On Adaptive Grad-Div Parameter Selection

abstract: We propose, analyze and test a new adaptive penalty scheme that picks the penalty parameter $\epsilon$ element by element small where $\nabla\cdot u^h$ is large. We start by analyzing and testing the new scheme on the most simple but interesting setting, the Stokes problem. Finally, we extend and test the algorithm on the incompressible Navier Stokes equation on complex flow problems. Tests indicate that the new adaptive-$\epsilon$ penalty method algorithm predicts flow behavior accurately. The scheme is developed in the penalty method but also can be used to pick a grad-div stabilization parameter.


title: Degree bounds for invariant skew polynomials

abstract: When we consider the action of a finite group on a polynomial ring, a polynomial unchanged by the action is called an invariant polynomial. A famous result of Noether states that in characteristic zero the maximal degree of a minimal invariant polynomial is bounded above by the order of the group. Our work establishes that the same bound holds for invariant skew polynomials in the exterior algebra. Our approach to the problem relies on a theorem of Derksen that connects invariant theory to the study of ideals of subspace arrangements. We adapt his proof over the polynomial ring to the exterior algebra, reducing the question to establishing a bound on the Castelnuovo-Mumford regularity of intersections of linear ideals in the exterior algebra. We prove the required regularity bound using tools from representation theory. In particular, the proof relies on the existence of a functor on the category of polynomial functors that translates resolutions of ideals of subspace arrangements over the polynomial ring to resolutions of ideals of subspace arrangements over the exterior algebra.


title: Normalized Augmented Inverse Probability Weighting with Neural Network Predictions

abstract: The estimation of Average Treatment Effect (ATE) as a causal parameter is carried out in two steps, where in the first step, the treatment and outcome are modeled to incorporate the potential confounders, and in the second step, the predictions are inserted into the ATE estimators such as the Augmented Inverse Probability Weighting (AIPW) estimator. Due to the concerns regarding the nonlinear or unknown relationships between confounders and the treatment and outcome, there has been an interest in applying non-parametric methods such as Machine Learning (ML) algorithms instead. Some literature proposes to use two separate Neural Networks (NNs) where there's no regularization on the network's parameters except the Stochastic Gradient Descent (SGD) in the NN's optimization. Our simulations indicate that the AIPW estimator suffers extensively if no regularization is utilized. We propose the normalization of AIPW (referred to as nAIPW) which can be helpful in some scenarios. nAIPW, provably, has the same properties as AIPW, that is, the double-robustness and orthogonality properties. Further, if the first step algorithms converge fast enough, under regulatory conditions, nAIPW will be asymptotically normal. We also compare the performance of AIPW and nAIPW in terms of the bias and variance when small to moderate L1 regularization is imposed on the NNs.


title: An Empirical Evaluation of End-to-End Polyphonic Optical Music Recognition

abstract: Previous work has shown that neural architectures are able to perform optical music recognition (OMR) on monophonic and homophonic music with high accuracy. However, piano and orchestral scores frequently exhibit polyphonic passages, which add a second dimension to the task. Monophonic and homophonic music can be described as homorhythmic, or having a single musical rhythm. Polyphonic music, on the other hand, can be seen as having multiple rhythmic sequences, or voices, concurrently. We first introduce a workflow for creating large-scale polyphonic datasets suitable for end-to-end recognition from sheet music publicly available on the MuseScore forum. We then propose two novel formulations for end-to-end polyphonic OMR -- one treating the problem as a type of multi-task binary classification, and the other treating it as multi-sequence detection. Building upon the encoder-decoder architecture and an image encoder proposed in past work on end-to-end OMR, we propose two novel decoder models -- FlagDecoder and RNNDecoder -- that correspond to the two formulations. Finally, we compare the empirical performance of these end-to-end approaches to polyphonic OMR and observe a new state-of-the-art performance with our multi-sequence detection decoder, RNNDecoder.


title: Atomic scale control of spin current transmission at interfaces

abstract: Spin transmission at ferromagnet/heavy metal interfaces is of vital importance for many spintronic devices. Usually the spin current transmission is limited by the spin mixing conductance and loss mechanisms such as spin memory loss. In order to understand these effects, we study the interface transmission when an insulating interlayer is inserted between the ferromagnet and the heavy metal. For this we measure the inverse spin Hall voltage generated from optically injected spin current pulses as well as the magnitude of the spin pumping using ferromagnetic resonance. From our results we conclude that significant spin memory loss only occurs for 5d metals with less than half filled d-shell.


title: On Exponential Utility and Conditional Value-at-Risk as Risk-Averse Performance Criteria

abstract: The standard approach to risk-averse control is to use the Exponential Utility (EU) functional, which has been studied for several decades. Like other risk-averse utility functionals, EU encodes risk aversion through an increasing convex mapping $\varphi$ of objective costs to subjective costs. An objective cost is a realization $y$ of a random variable $Y$. In contrast, a subjective cost is a realization $\varphi(y)$ of a random variable $\varphi(Y)$ that has been transformed to measure preferences about the outcomes. For EU, the transformation is $\varphi(y) = \exp(\frac{-\theta}{2}y)$, and under certain conditions, the quantity $\varphi^{-1}(E(\varphi(Y)))$ can be approximated by a linear combination of the mean and variance of $Y$. More recently, there has been growing interest in risk-averse control using the Conditional Value-at-Risk (CVaR) functional. In contrast to the EU functional, the CVaR of a random variable $Y$ concerns a fraction of its possible realizations. If $Y$ is a continuous random variable with finite $E(|Y|)$, then the CVaR of $Y$ at level $\alpha$ is the expectation of $Y$ in the $\alpha \cdot 100 \%$ worst cases. Here, we study the applications of risk-averse functionals to controller synthesis and safety analysis through the development of numerical examples, with emphasis on EU and CVaR. Our contribution is to examine the decision-theoretic, mathematical, and computational trade-offs that arise when using EU and CVaR for optimal control and safety analysis. We are hopeful that this work will advance the interpretability and elucidate the potential benefits of risk-averse control technology.


title: Nonconvex Factorization and Manifold Formulations are Almost Equivalent in Low-rank Matrix Optimization

abstract: In this paper, we consider the geometric landscape connection of the widely studied manifold and factorization formulations in low-rank positive semidefinite (PSD) and general matrix optimization. We establish a sandwich relation on the spectrum of Riemannian and Euclidean Hessians at first-order stationary points (FOSPs). As a result of that, we obtain an equivalence on the set of FOSPs, second-order stationary points (SOSPs) and strict saddles between the manifold and the factorization formulations. In addition, we show the sandwich relation can be used to transfer more quantitative geometric properties from one formulation to another. Similarities and differences in the landscape connection under the PSD case and the general case are discussed. To the best of our knowledge, this is the first geometric landscape connection between the manifold and the factorization formulations for handling rank constraints, and it provides a geometric explanation for the similar empirical performance of factorization and manifold approaches in low-rank matrix optimization observed in the literature. In the general low-rank matrix optimization, the landscape connection of two factorization formulations (unregularized and regularized ones) is also provided. By applying these geometric landscape connections, in particular, the sandwich relation, we are able to solve unanswered questions in literature and establish stronger results in the applications on geometric analysis of phase retrieval, well-conditioned low-rank matrix optimization, and the role of regularization in factorization arising from machine learning and signal processing.

categories: math.OC cs.IT cs.LG eess.SP math.IT stat.ML


title: Long-range transport of 2D excitons with acoustic waves

abstract: Excitons are elementary optical excitation in semiconductors. The ability to manipulate and transport these quasiparticles would enable excitonic circuits and devices for quantum photonic technologies. Recently, interlayer excitons in 2D semiconductors have emerged as a promising candidate for engineering excitonic devices due to their long lifetime, large exciton binding energy, and gate tunability. However, the charge-neutral nature of the excitons leads to weak response to the in-plane electric field and thus inhibits transport beyond the diffusion length. Here, we demonstrate the directional transport of interlayer excitons in bilayer WSe2 driven by the propagating potential traps induced by surface acoustic waves (SAW). We show that at 100 K, the SAW-driven excitonic transport is activated above a threshold acoustic power and reaches 20 mm, a distance at least ten times longer than the diffusion length and only limited by the device size. Temperature-dependent measurement reveals the transition from the diffusion-limited regime at low temperature to the acoustic field-driven regime at elevated temperature. Our work shows that acoustic waves are an effective, contact-free means to control exciton dynamics and transport, promising for realizing 2D materials-based excitonic devices such as exciton transistors, switches, and transducers up to room temperature.


title: Black hole spectroscopy horizons for current and future gravitational wave detectors

abstract: Black hole spectroscopy is the proposal to observe multiple quasinormal modes in the ringdown of a binary black hole merger. In addition to the fundamental quadrupolar mode, overtones and higher harmonics may be present and detectable in the gravitational wave signal, allowing for tests of the no-hair theorem. We analyze in detail the strengths and weaknesses of the standard Rayleigh criterion supplied with a Fisher matrix error estimation, and we find that the criterion is useful, but too restrictive. Therefore we motivate the use of a conservative high Bayes factor threshold to obtain the black hole spectroscopy horizons of current and future detectors, i.e., the distance (averaged in sky location and binary inclination) up to which one or more additional modes can be detected and confidently distinguished from each other. We set up all of our searches for additional modes starting at $t = 10(M_1+M_2)$ after the peak amplitude in simulated signals of circular nonspinning binaries. An agnostic multimode analysis allows us to rank the subdominant modes: for nearly equal mass binaries we find $(\ell, m, n) = (2,2,1)$ and $(3,3,0)$ and, for very asymmetric binaries, $(3,3,0)$ and $(4,4,0)$, for the secondary and tertiary modes, respectively. At the current estimated rates for heavy stellar mass binary black hole mergers, with primary masses between 45 and 100 solar masses, we expect an event rate of mergers within our conservative estimate for the $(2,2,1)$ spectroscopy horizon of $0.03 - 0.10\ {\rm yr}^{-1}$ for LIGO at design sensitivity and $(0.6 - 2.4) \times 10^3\ {\rm yr}^{-1}$ for the future third generation ground-based detector Cosmic Explorer.


title: Solo-learn: A Library of Self-supervised Methods for Visual Representation Learning

abstract: This paper presents solo-learn, a library of self-supervised methods for visual representation learning. Implemented in Python, using Pytorch and Pytorch lightning, the library fits both research and industry needs by featuring distributed training pipelines with mixed-precision, faster data loading via Nvidia DALI, online linear evaluation for better prototyping, and many additional training tricks. Our goal is to provide an easy-to-use library comprising a large amount of Self-supervised Learning (SSL) methods, that can be easily extended and fine-tuned by the community. solo-learn opens up avenues for exploiting large-budget SSL solutions on inexpensive smaller infrastructures and seeks to democratize SSL by making it accessible to all. The source code is available at https://github.com/vturrisi/solo-learn.


title: How Can Datacenters Join the Smart Grid to Address the Climate Crisis? Using simulation to explore power and cost effects of direct participation in the energy market

abstract: Amidst the climate crisis, the massive introduction of renewable energy sources has brought tremendous challenges to both the power grid and its surrounding markets. As datacenters have become ever-larger and more powerful, they play an increasingly significant role in the energy arena. With their unique characteristics, datacenters have been proved to be well-suited for regulating the power grid yet currently provide little, if any, such active response. This problem is due to issues such as unsuitability of the market design, high complexity of the currently proposed solutions, as well as the potential risks thereof. This work aims to provide individual datacenters with insights on the feasibility and profitability of directly participating in the energy market. By modelling the power system of datacenters, and by conducting simulations on real-world datacenter traces, we demonstrate the substantial financial incentive for individual datacenters to directly participate in both the day-ahead and the balancing markets. In turn, we suggest a new short-term, direct scheme of market participation for individual datacenters in place of the current long-term, inactive participation. Furthermore, we develop a novel proactive DVFS scheduling algorithm that can both reduce energy consumption and save energy costs during the market participation of datacenters. Also, in developing this scheduler, we propose an innovative combination of machine learning methods and the DVFS technology that can provide the power grid with indirect demand response (DR). Our experimental results strongly support that individual datacenters can and should directly participate in the energy market both to save their energy costs and to curb their energy consumption, whilst providing the power grid with indirect DR.


title: Fluid-mediated sources of granular temperature at finite Reynolds numbers

abstract: We derive analytical solutions for hydrodynamic sources and sinks to granular temperature in moderately dense suspensions of elastic particles at finite Reynolds numbers. Modeling the neighbor-induced drag disturbances with a Langevin equation allows an exact solution for the joint fluctuating acceleration-velocity distribution function $P\left(v^{\prime},a^{\prime};t\right)$. Quadrant-conditioned covariance integrals of $P\left(v^{\prime},a^{\prime};t\right)$ yield the hydrodynamic source and sink that dictate the evolution of granular temperature. Analytical predictions are in agreement with benchmark data obtained from particle-resolved direct numerical simulations and show promise as a general theory from gas--solid to bubbly flows.


title: Armour: Generalizable Compact Self-Attention for Vision Transformers

abstract: Attention-based transformer networks have demonstrated promising potential as their applications extend from natural language processing to vision. However, despite the recent improvements, such as sub-quadratic attention approximation and various training enhancements, the compact vision transformers to date using the regular attention still fall short in comparison with its convnet counterparts, in terms of \textit{accuracy,} \textit{model size}, \textit{and} \textit{throughput}. This paper introduces a compact self-attention mechanism that is fundamental and highly generalizable. The proposed method reduces redundancy and improves efficiency on top of the existing attention optimizations. We show its drop-in applicability for both the regular attention mechanism and some most recent variants in vision transformers. As a result, we produced smaller and faster models with the same or better accuracies.


title: Approximate Q-conditional symmetries of partial differential equations

abstract: Following a recently introduced approach to approximate Lie symmetries of differential equations which is consistent with the principles of perturbative analysis of differential equations containing small terms, we analyze the case of approximate $Q$--conditional symmetries. An application of the method to a hyperbolic variant of a reaction--diffusion--convection equation is presented.


title: On qualitative analysis of a discrete time SIR epidemical model

abstract: The main purpose of this paper is to study the local dynamics and bifurcations of a discrete-time SIR epidemiological model. The existence and stability of disease-free and endemic fixed points are investigated along with a fairly complete classification of the systems bifurcations, in particular, a complete analysis on local stability and codimension 1 bifurcations in the parameter space. Sufficient conditions for positive trajectories are given. The existence of a 3-cycle is shown, which implies the existence of cycles of arbitrary length by the celebrated Sharkovskii's theorem. Generacity of some bifurcations is examined both analytically and through numerical computations. Bifurcation diagrams along with numerical simulations are presented. The system turns out to have both rich and interesting dynamics.


title: Lifted contact dynamics for efficient optimal control of rigid body systems with contacts

abstract: We propose a novel and efficient lifting approach for the optimal control of rigid-body systems with contacts to improve the convergence properties of Newton-type methods. To relax the high nonlinearity, we consider the state, acceleration, contact forces, and control input torques, as optimization variables and the inverse dynamics and acceleration constraints on the contact frames as equality constraints. We eliminate the update of the acceleration, contact forces, and their dual variables from the linear equation to be solved in each Newton-type iteration in an efficient manner. As a result, the computational cost per Newton-type iteration is almost identical to that of the conventional non-lifted Newton-type iteration that embeds contact dynamics in the state equation. We conducted numerical experiments on the whole-body optimal control of various quadrupedal gaits subject to the friction cone constraints considered in interior-point methods and demonstrated that the proposed method can significantly increase the convergence speed to more than twice that of the conventional non-lifted approach.


title: Revisiting geodesic observers in cosmology

abstract: Geodesic observers in cosmology are revisited. The coordinates based on freely falling observers introduced by Gautreau in de Sitter and Einstein-de Sitter spaces (and, previously, by Gautreau and Hoffmann in Schwarzschild space) are extended to general FLRW universes. We identify situations in which the relation between geodesic and comoving coordinates can be expressed explicitly in terms of elementary functions. In general, geodesic coordinates in cosmology turn out to be rather cumbersome and limited to the region below the apparent horizon.


title: Re-examination of the Time Structure of the SN1987A Neutrino Burst Data in Kamiokande-II

abstract: The seven seconds' gap in the Kamiokande-II SN1987A neutrino data is reexamined.


title: Featured Team Automata

abstract: We propose featured team automata to support variability in the development and analysis of teams, which are systems of reactive components that communicate according to specified synchronisation types. A featured team automaton concisely describes a family of concrete product models for specific configurations determined by feature selection. We focus on the analysis of communication-safety properties, but doing so product-wise quickly becomes impractical. Therefore, we investigate how to lift notions of receptiveness (no message loss) to the level of family models. We show that featured (weak) receptiveness of featured team automata characterises (weak) receptiveness for all product instantiations. A prototypical tool supports the developed theory.


title: Weakly Supervised Foreground Learning for Weakly Supervised Localization and Detection

abstract: Modern deep learning models require large amounts of accurately annotated data, which is often difficult to satisfy. Hence, weakly supervised tasks, including weakly supervised object localization~(WSOL) and detection~(WSOD), have recently received attention in the computer vision community. In this paper, we motivate and propose the weakly supervised foreground learning (WSFL) task by showing that both WSOL and WSOD can be greatly improved if groundtruth foreground masks are available. More importantly, we propose a complete WSFL pipeline with low computational cost, which generates pseudo boxes, learns foreground masks, and does not need any localization annotations. With the help of foreground masks predicted by our WSFL model, we achieve 72.97% correct localization accuracy on CUB for WSOL, and 55.7% mean average precision on VOC07 for WSOD, thereby establish new state-of-the-art for both tasks. Our WSFL model also shows excellent transfer ability.


title: Numerical Algorithms for Water Waves with Background Flow over Obstacles and Topography

abstract: We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply-connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

categories: physics.flu-dyn cs.NA math-ph math.AP math.MP math.NA


title: The IGRINS YSO Survey I. Stellar parameters of pre-main sequence stars in Taurus-Auriga

abstract: We present fundamental parameters for 110 canonical K- & M-type (1.3$-$0.13$M_\odot$) Taurus-Auriga young stellar objects (YSOs). The analysis produces a simultaneous determination of effective temperature ($T_{\rm eff}$), surface gravity ($\log$ g), magnetic field strength (B), and projected rotational velocity ($v \sin i$). Our method employed synthetic spectra and high-resolution (R$\sim$45,000) near-infrared spectra taken with the Immersion GRating INfrared Spectrometer (IGRINS) to fit specific K-band spectral regions most sensitive to those parameters. The use of these high-resolution spectra reduces the influence of distance uncertainties, reddening, and non-photospheric continuum emission on the parameter determinations. The median total (fit + systematic) uncertainties were 170 K, 0.28 dex, 0.60 kG, 2.5 km s$^{-1}$ for $T_{\rm eff}$, $\log$ g, B, and $v \sin i$, respectively. We determined B for 41 Taurus YSOs (upper limits for the remainder) and find systematic offsets (lower $T_{\rm eff}$, higher $\log$ g and $v \sin i$) in parameters when B is measurable but not considered in the fit. The average $\log$ g for the Class II and Class III objects differs by 0.23$\pm$0.05dex, which is consistent with Class III objects being the more evolved members of the star-forming region. However, the dispersion in $\log$ g is greater than the uncertainties, which highlights how the YSO classification correlates with age ($\log$ g), yet there are exceptionally young (lower $\log$ g) Class III YSOs and relatively old (higher $\log$ g) Class II YSOs with unexplained evolutionary histories. The spectra from this work are provided in an online repository along with TW Hydrae Association (TWA) comparison objects and the model grid used in our analysis.


title: Interaction of Electromagnetic field with a Gravitational wave in Minkowski and de-Sitter space-time

abstract: It is important to re-examine some of the interactions of gravitational waves with matter due to the recent discovery of the waves in LIGO. In a previous paper with Morales, interaction of gravitational waves was studied with scalar and neutrino fields in Minkowski space-time. In this paper, we find that electromagnetic waves have a similar interaction with gravitational wave in a flat background, and we comment on the implications of this for a gravitational wave detector. Primordial gravitational waves were generated in the early universe and thus it is important to study interaction of the gravitational waves in cosmological backgrounds. In this paper we discuss the interaction of electromagnetic waves and gravitational waves in a de-Sitter background. We find that the electromagnetic wave is perturbed and modulated. This result is new and will be useful in studying inhomogeneities in Cosmic Microwave Background.


title: Monte Carlo Tree Search for high precision manufacturing

abstract: Monte Carlo Tree Search (MCTS) has shown its strength for a lot of deterministic and stochastic examples, but literature lacks reports of applications to real world industrial processes. Common reasons for this are that there is no efficient simulator of the process available or there exist problems in applying MCTS to the complex rules of the process. In this paper, we apply MCTS for optimizing a high-precision manufacturing process that has stochastic and partially observable outcomes. We make use of an expert-knowledge-based simulator and adapt the MCTS default policy to deal with the manufacturing process.


title: A Comparative Study of Atmospheric Chemistry with VULCAN

abstract: We present an update of the open-source photochemical kinetics code VULCAN (Tsai et al. 2017; https://github.com/exoclime/VULCAN) to include C-H-N-O-S networks and photochemistry. Additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross-sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an intercomparison of HD 189733b models between Moses et al. (2011), Venot et al. (2012), and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep troposphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and ammonia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP-33b, HD 189733b, GJ 436b, and 51 Eridani b. We look into the effects of the C/O ratio and chemistry of titanium/vanadium species for WASP-33b; we revisit HD 189733b for the effects of sulfur and carbon condensation; the effects of internal heating and vertical mixing ($K_{\textrm{zz}}$) are explored for GJ 436b; we test updated planetary properties for 51 Eridani b with S$_8$ condensates. We find sulfur can couple to carbon or nitrogen and impact other species such as hydrogen, methane, and ammonia. The observable features of the synthetic spectra and trends in the photochemical haze precursors are discussed for each case.


title: Impact of heavy-quark production measurements in the CT18 global QCD analysis of PDFs

abstract: We discuss the impact of heavy-flavor production measurements in semi-inclusive deep inelastic scattering at HERA on the CTEQ-TEA PDFs. In particular, we study the impact of the latest charm and bottom production measurements from the H1 and ZEUS collaborations on the gluon, and the interplay of these measurements with the data ensemble of the recent CT18 global QCD analysis.


title: Revisiting the Okubo-Marshak argument

abstract: Modular localization and the theory of string-localized fields have revolutionized several key aspects of quantum field theory. They reinforce the contention that local symmetry emerges directly from quantum theory, but global gauge invariance remains in general an unwarranted assumption, to be examined case by case. Armed with those modern tools, we reconsider here the classical Okubo-Marshak argument on the non-existence of a "strong CP problem" in quantum chromodynamics.


title: Composing Recurrent Spiking Neural Networks using Locally-Recurrent Motifs and Risk-Mitigating Architectural Optimization

abstract: In neural circuits, recurrent connectivity plays a crucial role in network function and stability. However, existing recurrent spiking neural networks (RSNNs) are often constructed by random connections without optimization. While RSNNs can produce rich dynamics that are critical for memory formation and learning, systemic architectural optimization of RSNNs is still an open challenge. We aim to enable systematic design of large RSNNs via a new scalable RSNN architecture and automated architectural optimization. We compose RSNNs based on a layer architecture called Sparsely-Connected Recurrent Motif Layer (SC-ML) that consists of multiple small recurrent motifs wired together by sparse lateral connections. The small size of the motifs and sparse inter-motif connectivity leads to an RSNN architecture scalable to large network sizes. We further propose a method called Hybrid Risk-Mitigating Architectural Search (HRMAS) to systematically optimize the topology of the proposed recurrent motifs and SC-ML layer architecture. HRMAS is an alternating two-step optimization process by which we mitigate the risk of network instability and performance degradation caused by architectural change by introducing a novel biologically-inspired "self-repairing" mechanism through intrinsic plasticity. The intrinsic plasticity is introduced to the second step of each HRMAS iteration and acts as unsupervised fast self-adaptation to structural and synaptic weight modifications introduced by the first step during the RSNN architectural "evolution". To the best of the authors' knowledge, this is the first work that performs systematic architectural optimization of RSNNs. Using one speech and three neuromorphic datasets, we demonstrate the significant performance improvement brought by the proposed automated architecture optimization over existing manually-designed RSNNs.


title: Explicit RIP matrices: an update

abstract: Leveraging recent advances in additive combinatorics, we exhibit explicit matrices satisfying the Restricted Isometry Property with better parameters. Namely, for $\varepsilon=3.26\cdot 10^{-7}$, large $k$ and $k^{2-\varepsilon} \le N\le k^{2+\varepsilon}$, we construct $n \times N$ RIP matrices of order $k$ with $k = \Omega( n^{1/2+\varepsilon/4} )$.


title: Pits, Uplifts and Small Chaos Features on Europa: Morphologic and Morphometric Evidence for Intrusive Upwelling and Lower Limits to Ice Shell Thickness

abstract: One of the clearest but unresolved questions for Europa is the thickness of its icy shell. Europa's surface is resplendent with geological features that bear on this question, and ultimately on its interior, geological history, and astrobiological potential. We characterize the size and topographic expression of circular and subcircular features created by endogenic thermal and tectonic disturbances on Europa: pits, uplifts, and small, subcircular chaos. We utilize the medium-resolution Galileo regional maps (RegMaps), as well as high-resolution regions, digital elevation models derived from albedo-controlled photoclinometry, and in some cases stereo-controlled photoclinometry. While limited in extent, the high-resolution images are extremely valuable for detecting smaller features and for overall geomorphological analysis. A peak in the size-distribution for all features is found at ~5-6 km in diameter and no pits smaller than 3.3 km in diameter were found in high resolution images. Additionally, there is a trend for larger pits to be deeper, and larger uplifts to be higher. Our data support a diapiric or intracrustal sill interpretation (as opposed to purely non-intrusive, melt-through models) and place a lower limit on ice shell thickness at the time of feature formation of 3-to-8 km, assuming isostasy and depending on the composition of the ice and underlying ocean.


title: Counting essential minimal surfaces in closed negatively curved n-manifolds

abstract: For closed odd-dimensional manifolds with sectional curvature less or equal than -1, we define the minimal surface entropy that counts the number of surface subgroups. It attains the minimum if and only if the metric is hyperbolic. Moreover, we compute the entropy associated with other locally symmetric spaces and cusped hyperbolic 3-manifolds.


title: Rapid Convex Optimization of Centroidal Dynamics using Block Coordinate Descent

abstract: In this paper we explore the use of block coordinate descent (BCD) to optimize the centroidal momentum dynamics for dynamically consistent multi-contact behaviors. The centroidal dynamics have recently received a large amount of attention in order to create physically realizable motions for robots with hands and feet while being computationally more tractable than full rigid body dynamics models. Our contribution lies in exploiting the structure of the dynamics in order to simplify the original non-convex problem into two convex subproblems. We iterate between these two subproblems for a set number of iterations or until a consensus is reached. We explore the properties of the proposed optimization method for the centroidal dynamics and verify in simulation that motions generated by our approach can be tracked by the quadruped Solo12. In addition, we compare our method to a recently proposed convexification using a sequence of convex relaxations as well as a more standard interior point method used in the off- the-shelf solver IPOPT to show that our approach finds similar, if not better, trajectories (in terms of cost), and is more than four times faster than both approaches. Finally, compared to previous approaches, we note its practicality due to the convex nature of each subproblem which allows our method to be used with any off-the-shelf quadratic programming solver.


title: Extremal polynomials on a Jordan arc

abstract: Let $\Gamma$ be a $C^{2+}$ Jordan arc and let $\Gamma_0$ be the open arc which consists of interior points of $\Gamma$. We find concrete upper and lower bounds for the limit of Widom factors for $L_2(\mu)$ extremal polynomials on $\Gamma$ which was given in [18]. In addition, we show that the upper bound for the limit supremum of Widom factors for the weighted Chebyshev polynomials which was obtained in [18] can be improved once two normal derivatives of the Green function do not agree at one point $z\in \Gamma_0$. We also show that if $\Gamma_0$ is not analytic then we have improved upper bounds.


title: Goldilocks: Consistent Crowdsourced Scalar Annotations with Relative Uncertainty

abstract: Human ratings have become a crucial resource for training and evaluating machine learning systems. However, traditional elicitation methods for absolute and comparative rating suffer from issues with consistency and often do not distinguish between uncertainty due to disagreement between annotators and ambiguity inherent to the item being rated. In this work, we present Goldilocks, a novel crowd rating elicitation technique for collecting calibrated scalar annotations that also distinguishes inherent ambiguity from inter-annotator disagreement. We introduce two main ideas: grounding absolute rating scales with examples and using a two-step bounding process to establish a range for an item's placement. We test our designs in three domains: judging toxicity of online comments, estimating satiety of food depicted in images, and estimating age based on portraits. We show that (1) Goldilocks can improve consistency in domains where interpretation of the scale is not universal, and that (2) representing items with ranges lets us simultaneously capture different sources of uncertainty leading to better estimates of pairwise relationship distributions.


title: On optimality of barrier dividend control under endogenous regime switching with application to Chapter 11 bankruptcy

abstract: Motivated by recent developments in risk management based on the U.S. bankruptcy code, we revisit the De Finetti's optimal dividend problem by incorporating the reorganization process and regulator's intervention documented in Chapter 11 bankruptcy. The resulting surplus process, bearing financial stress towards the more subtle concept of bankruptcy, corresponds to a non-standard spectrally negative Levy process with endogenous regime switching. Some explicit expressions of the expected present values under a barrier strategy, new to the literature, are established in terms of scale functions. With the help of these expressions, when the tail of the Levy measure is log-convex, the optimal dividend control is shown to be of the barrier type and the associated optimal barrier can be identified using scale functions of spectrally negative Levy processes. Some financial implications are also discussed in an illustrative example.


title: Discovering Supernovae at Epoch of Reionization with Nancy Grace Roman Space Telescope

abstract: Massive stars play critical roles for the reionization of the Universe. Individual massive stars at the reionization epoch (z > 6) are too faint to observe and quantify their contributions to reionization. Some massive stars, however, explode as superluminous supernovae (SLSNe) or pair-instability supernovae (PISNe) that are luminous enough to observe even at z > 6 and allow for the direct characterization of massive star properties at the reionization epoch. In addition, hypothetical long-sought-after PISNe are expected to be present preferentially at high redshifts, and their discovery will have a tremendous impact on our understanding of massive star evolution and the formation of stellar mass black holes. The near-infrared Wide Field Instrument on Nancy Grace Roman Space Telescope will excel at discovering such rare high-redshift supernovae. In this work, we investigate the best survey strategy to discover and identify SLSNe and PISNe at z > 6 with Roman. We show that the combination of the F158 and F213 filters can clearly separate both SLSNe and PISNe at z > 6 from nearby supernovae through their colors and magnitudes. The limiting magnitudes are required to be 27.0 mag and 26.5 mag in the F158 and F213 filters, respectively, to identify supernovae at z > 6. If we conduct a 10 deg2 transient survey with these limiting magnitudes for 5 years with a cadence of one year, we expect to discover 22.5 +- 2.8 PISNe and 3.1 +- 0.3 SLSNe at z > 6, depending on the cosmic star-formation history. The same survey is estimated to discover 76.1 +- 8.2 PISNe and 9.1 +- 0.9 SLSNe at 5 < z < 6. Such a supernova survey requires the total observational time of approximately 525 hours in 5 years. The legacy data acquired with the survey will also be beneficial for many different science cases including the study of high-redshift galaxies.

categories: astro-ph.HE astro-ph.CO astro-ph.GA astro-ph.SR


title: Deformation Recovery Control and Post-Impact Trajectory Replanning for Collision-Resilient Mobile Robots

abstract: The paper focuses on collision-inclusive motion planning for impact-resilient mobile robots. We propose a new deformation recovery and replanning strategy to handle collisions that may occur at run-time. Contrary to collision avoidance methods that generate trajectories only in conservative local space or require collision checking that has high computational cost, our method directly generates (local) trajectories with imposing only waypoint constraints. If a collision occurs, our method then estimates the post-impact state and computes from there an intermediate waypoint to recover from the collision. To achieve so, we develop two novel components: 1) a deformation recovery controller that optimizes the robot's states during post-impact recovery phase, and 2) a post-impact trajectory replanner that adjusts the next waypoint with the information from the collision for the robot to pass through and generates a polynomial-based minimum effort trajectory. The proposed strategy is evaluated experimentally with an omni-directional impact-resilient wheeled robot. The robot is designed in house, and it can perceive collisions with the aid of Hall effect sensors embodied between the robot's main chassis and a surrounding deflection ring-like structure.


title: Twist-3 Double Spin Asymmetries in Drell-Yan Processes

abstract: We study double spin asymmetries in Drell-Yan processes in which one initial hadron is transversely polarized and another one is longitudinally polarized. The complete part of the hadronic tensor relevant to asymmetries is derived. This part consists of twist-2 and twist-3 parton distributions and is gauge invariant. We construct some observables which can be used to extract these parton distributions from experimental measurements.


title: Online Training of Spiking Recurrent Neural Networks with Phase-Change Memory Synapses

abstract: Spiking recurrent neural networks (RNNs) are a promising tool for solving a wide variety of complex cognitive and motor tasks, due to their rich temporal dynamics and sparse processing. However training spiking RNNs on dedicated neuromorphic hardware is still an open challenge. This is due mainly to the lack of local, hardware-friendly learning mechanisms that can solve the temporal credit assignment problem and ensure stable network dynamics, even when the weight resolution is limited. These challenges are further accentuated, if one resorts to using memristive devices for in-memory computing to resolve the von-Neumann bottleneck problem, at the expense of a substantial increase in variability in both the computation and the working memory of the spiking RNNs. To address these challenges and enable online learning in memristive neuromorphic RNNs, we present a simulation framework of differential-architecture crossbar arrays based on an accurate and comprehensive Phase-Change Memory (PCM) device model. We train a spiking RNN whose weights are emulated in the presented simulation framework, using a recently proposed e-prop learning rule. Although e-prop locally approximates the ideal synaptic updates, it is difficult to implement the updates on the memristive substrate due to substantial PCM non-idealities. We compare several widely adapted weight update schemes that primarily aim to cope with these device non-idealities and demonstrate that accumulating gradients can enable online and efficient training of spiking RNN on memristive substrates.


title: Searching for a Galactic component in the IceCube track-like neutrino events

abstract: Searches for spatial associations between high-energy neutrinos observed at the IceCube Neutrino Observatory and known astronomical objects may hold the key to establishing the neutrinos' origins and the origins of hadronic cosmic rays. While extragalactic sources like the blazar TXS 0506+056 merit significant attention, Galactic sources may also represent part of the puzzle. Here, we explore whether open clusters and supernova remnants in the Milky Way contribute measurably to the IceCube track-like neutrino events above 200 TeV. By searching for positional coincidences with catalogs of known astronomical objects, we can identify and investigate neutrino events whose origins are potentially Galactic. We use Monte Carlo randomization together with models of the Galactic plane in order to determine whether these coincidences are more likely to be causal associations or random chance. In all analyses presented, the number of coincidences detected was found to be consistent with the null hypothesis of chance coincidence. Our results imply that the combined contribution of Galactic open clusters and supernova remnants to the track-like neutrino events detected at IceCube is well under 30%. This upper limit is compatible with the results presented in other Galactic neutrino studies.


title: Neural Scene Decoration from a Single Photograph

abstract: Furnishing and rendering indoor scenes has been a long-standing task for interior design, where artists create a conceptual design for the space, build a 3D model of the space, decorate, and then perform rendering. Although the task is important, it is tedious and requires tremendous effort. In this paper, we introduce a new problem of domain-specific indoor scene image synthesis, namely neural scene decoration. Given a photograph of an empty indoor space and a list of decorations with layout determined by user, we aim to synthesize a new image of the same space with desired furnishing and decorations. Neural scene decoration can be applied to create conceptual interior designs in a simple yet effective manner. Our attempt to this research problem is a novel scene generation architecture that transforms an empty scene and an object layout into a realistic furnished scene photograph. We demonstrate the performance of our proposed method by comparing it with conditional image synthesis baselines built upon prevailing image translation approaches both qualitatively and quantitatively. We conduct extensive experiments to further validate the plausibility and aesthetics of our generated scenes. Our implementation is available at \url{https://github.com/hkust-vgd/neural_scene_decoration}.


title: On the Robustness of Domain Adaption to Adversarial Attacks

abstract: State-of-the-art deep neural networks (DNNs) have been proved to have excellent performance on unsupervised domain adaption (UDA). However, recent work shows that DNNs perform poorly when being attacked by adversarial samples, where these attacks are implemented by simply adding small disturbances to the original images. Although plenty of work has focused on this, as far as we know, there is no systematic research on the robustness of unsupervised domain adaption model. Hence, we discuss the robustness of unsupervised domain adaption against adversarial attacking for the first time. We benchmark various settings of adversarial attack and defense in domain adaption, and propose a cross domain attack method based on pseudo label. Most importantly, we analyze the impact of different datasets, models, attack methods and defense methods. Directly, our work proves the limited robustness of unsupervised domain adaptation model, and we hope our work may facilitate the community to pay more attention to improve the robustness of the model against attacking.


title: Leaf Recognition Using Convolutional Neural Networks Based Features

abstract: There is a warning light for the loss of plant habitats worldwide that entails concerted efforts to conserve plant biodiversity. Thus, plant species classification is of crucial importance to address this environmental challenge. In recent years, there is a considerable increase in the number of studies related to plant taxonomy. While some researchers try to improve their recognition performance using novel approaches, others concentrate on computational optimization for their framework. In addition, a few studies are diving into feature extraction to gain significantly in terms of accuracy. In this paper, we propose an effective method for the leaf recognition problem. In our proposed approach, a leaf goes through some pre-processing to extract its refined color image, vein image, xy-projection histogram, handcrafted shape, texture features, and Fourier descriptors. These attributes are then transformed into a better representation by neural network-based encoders before a support vector machine (SVM) model is utilized to classify different leaves. Overall, our approach performs a state-of-the-art result on the Flavia leaf dataset, achieving the accuracy of 99.58\% on test sets under random 10-fold cross-validation and bypassing the previous methods. We also release our codes (Scripts are available at https://github.com/dinhvietcuong1996/LeafRecognition) for contributing to the research community in the leaf classification problem.


title: What's Wrong with the Bottom-up Methods in Arbitrary-shape Scene Text Detection

abstract: The latest trend in the bottom-up perspective for arbitrary-shape scene text detection is to reason the links between text segments using Graph Convolutional Network (GCN). Notwithstanding, the performance of the best performing bottom-up method is still inferior to that of the best performing top-down method even with the help of GCN. We argue that this is not mainly caused by the limited feature capturing ability of the text proposal backbone or GCN, but by their failure to make a full use of visual-relational features for suppressing false detection, as well as the sub-optimal route-finding mechanism used for grouping text segments. In this paper, we revitalize the classic text detection frameworks by aggregating the visual-relational features of text with two effective false positive/negative suppression mechanisms. First, dense overlapping text segments depicting the `characterness' and `streamline' of text are generated for further relational reasoning and weakly supervised segment classification. Here, relational graph features are used for suppressing false positives/negatives. Then, to fuse the relational features with visual features, a Location-Aware Transfer (LAT) module is designed to transfer text's relational features into visual compatible features with a Fuse Decoding (FD) module to enhance the representation of text regions for the second step suppression. Finally, a novel multiple-text-map-aware contour-approximation strategy is developed, instead of the widely-used route-finding process. Experiments conducted on five benchmark datasets, i.e., CTW1500, Total-Text, ICDAR2015, MSRA-TD500, and MLT2017 demonstrate that our method outperforms the state-of-the-art performance when being embedded in a classic text detection framework, which revitalises the superb strength of the bottom-up methods.


title: Deep Learning Chromatic and Clique Numbers of Graphs

abstract: Deep neural networks have been applied to a wide range of problems across different application domains with great success. Recently, research into combinatorial optimization problems in particular has generated much interest in the machine learning community. In this work, we develop deep learning models to predict the chromatic number and maximum clique size of graphs, both of which represent classical NP-complete combinatorial optimization problems encountered in graph theory. The neural networks are trained using the most basic representation of the graph, the adjacency matrix, as opposed to undergoing complex domain-specific feature engineering. The experimental results show that deep neural networks, and in particular convolutional neural networks, obtain strong performance on this problem.


title: Infinite growth in vorticity gradient of compactly supported planar vorticity near Lamb dipole

abstract: We prove linear in time filamentation for perturbations of the Lamb dipole, which is a traveling wave solution to the incompressible Euler equations in $\mathbb{R}^2$. The main ingredient is a recent nonlinear orbital stability result by Abe-Choi. As a consequence, we obtain linear in time growth for the vorticity gradient for all times, for certain smooth and compactly supported initial vorticity in $\mathbb{R}^2$. The construction carries over to some generalized SQG equations.


title: Improving Distinction between ASR Errors and Speech Disfluencies with Feature Space Interpolation

abstract: Fine-tuning pretrained language models (LMs) is a popular approach to automatic speech recognition (ASR) error detection during post-processing. While error detection systems often take advantage of statistical language archetypes captured by LMs, at times the pretrained knowledge can hinder error detection performance. For instance, presence of speech disfluencies might confuse the post-processing system into tagging disfluent but accurate transcriptions as ASR errors. Such confusion occurs because both error detection and disfluency detection tasks attempt to identify tokens at statistically unlikely positions. This paper proposes a scheme to improve existing LM-based ASR error detection systems, both in terms of detection scores and resilience to such distracting auxiliary tasks. Our approach adopts the popular mixup method in text feature space and can be utilized with any black-box ASR output. To demonstrate the effectiveness of our method, we conduct post-processing experiments with both traditional and end-to-end ASR systems (both for English and Korean languages) with 5 different speech corpora. We find that our method improves both ASR error detection F 1 scores and reduces the number of correctly transcribed disfluencies wrongly detected as ASR errors. Finally, we suggest methods to utilize resulting LMs directly in semi-supervised ASR training.


title: Gaussian Process Regression and Classification using International Classification of Disease Codes as Covariates

abstract: International Classification of Disease (ICD) codes are widely used for encoding diagnoses in electronic health records (EHR). Automated methods have been developed over the years for predicting biomedical responses using EHR that borrow information among diagnostically similar patients. Relatively less attention has been paid to developing patient similarity measures that model the structure of ICD codes and the presence of multiple chronic conditions, where a chronic condition is defined as a set of ICD codes. Motivated by this problem, we first develop a type of string kernel function for measuring similarity between a pair of subsets of ICD codes, which uses the definition of chronic conditions. Second, we extend this similarity measure to define a family of covariance functions on subsets of ICD codes. Using this family, we develop Gaussian process (GP) priors for Bayesian nonparametric regression and classification using diagnoses and other demographic information as covariates. Markov chain Monte Carlo (MCMC) algorithms are used for posterior inference and predictions. The proposed methods are free of any tuning parameters and are well-suited for automated prediction of continuous and categorical biomedical responses that depend on chronic conditions. We evaluate the practical performance of our method on EHR data collected from 1660 patients at the University of Iowa Hospitals and Clinics (UIHC) with six different primary cancer sites. Our method has better sensitivity and specificity than its competitors in classifying different primary cancer sites and estimates the marginal associations between chronic conditions and primary cancer sites.


title: Particle resolved simulation of sediment transport by a hybrid parallel approach

abstract: Sediment transport over an erodible sediment bed is studied by particle resolved simulations with a hybrid parallel approach. To overcome the challenges of load imbalance in the traditional domain decomposition method when encountering highly uneven distributions of particles in sediment transport, the parallel approach of Darmana et al.(2006) originally developed for point particle simulations is modified and implemented into particle resolved simulations. A novel memory optimization technique is proposed to reduce the memory requirement of the hybrid approach for spherical particles with equal size. The present hybrid parallel approach shows good scalability and high parallel efficiency in a challenging sediment transport test case with more than a million spherical particles. Our code is validated by several benchmark cases, and the results show good agreement with experimental and computational data in the literature. Furthermore, a turbulent flow over an erodible sediment bed is simulated. An extraction method is proposed to distinguish the saltating and rolling particles and extract impact and rebound information of the particle-mobile bed interaction. The probability distribution functions (PDF) of several saltation parameters such as velocity, angle, and spanwise angular velocity of impact and rebound events are presented. Splash functions are established for the particle-mobile bed interaction in the turbulent flow, which was rarely investigated in the experiments and is helpful to model the complex particle-bed interactions in turbulent flow.


title: Topology optimization of the support structure for heat dissipation in additive manufacturing

abstract: A support structure is required to successfully create structural parts in the powder bed fusion process for additive manufacturing. In this study, we present the topology optimization of a support structure that improves the heat dissipation in the building process. First, we construct a numerical method that obtains the temperature field in the building process, represented by the transient heat conduction phenomenon with the volume heat flux. Next, we formulate an optimization problem for maximizing heat dissipation and develop an optimization algorithm that incorporates a level-set-based topology optimization. A sensitivity of the objective function is derived using the adjoint variable method. Finally, several numerical examples are provided to demonstrate the effectiveness and validity of the proposed method.


title: Pseudo-Riemannian Lie groups of modified H-type

abstract: We define a class of Riemannian and pseudo-Riemannian 2-step nilpotent Lie groups with nondegenerate centers that generalize the H-type groups of Kaplan. Examples are given and geometric properties are investigated.


title: Video Similarity and Alignment Learning on Partial Video Copy Detection

abstract: Existing video copy detection methods generally measure video similarity based on spatial similarities between key frames, neglecting the latent similarity in temporal dimension, so that the video similarity is biased towards spatial information. There are methods modeling unified video similarity in an end-to-end way, but losing detailed partial alignment information, which causes the incapability of copy segments localization. To address the above issues, we propose the Video Similarity and Alignment Learning (VSAL) approach, which jointly models spatial similarity, temporal similarity and partial alignment. To mitigate the spatial similarity bias, we model the temporal similarity as the mask map predicted from frame-level spatial similarity, where each element indicates the probability of frame pair lying right on the partial alignments. To further localize partial copies, the step map is learned from the spatial similarity where the elements indicate extending directions of the current partial alignments on the spatial-temporal similarity map. Obtained from the mask map, the start points extend out into partial optimal alignments following instructions of the step map. With the similarity and alignment learning strategy, VSAL achieves the state-of-the-art F1-score on VCDB core dataset. Furthermore, we construct a new benchmark of partial video copy detection and localization by adding new segment-level annotations for FIVR-200k dataset, where VSAL also achieves the best performance, verifying its effectiveness in more challenging situations. Our project is publicly available at https://pvcd-vsal.github.io/vsal/.


title: Quasisymmetric magnetic fields in asymmetric toroidal domains

abstract: We explore the existence of quasisymmetric magnetic fields in asymmetric toroidal domains. These vector fields can be identified with a class of magnetohydrodynamic equilibria in the presence of pressure anisotropy. First, using Clebsch potentials, we derive a system of two coupled nonlinear first order partial differential equations expressing a family of quasisymmetric magnetic fields in bounded domains. In regions where flux surfaces and surfaces of constant field strength are not tangential, this system can be further reduced to a single degenerate nonlinear second order partial differential equation with externally assigned initial data. Then, we exhibit regular quasisymmetric vector fields which correspond to local solutions of anisotropic magnetohydrodynamics in asymmetric toroidal domains such that tangential boundary conditions are fulfilled on a portion of the bounding surface. The problems of boundary shape and locality are also discussed. We find that symmetric magnetic fields can be fitted into asymmetric domains, and that the mathematical difficulty encountered in the derivation of global quasisymmetric magnetic fields lies in the topological obstruction toward global extension affecting local solutions of the governing nonlinear first order partial differential equations.


title: Transfer Learning for Pose Estimation of Illustrated Characters

abstract: Human pose information is a critical component in many downstream image processing tasks, such as activity recognition and motion tracking. Likewise, a pose estimator for the illustrated character domain would provide a valuable prior for assistive content creation tasks, such as reference pose retrieval and automatic character animation. But while modern data-driven techniques have substantially improved pose estimation performance on natural images, little work has been done for illustrations. In our work, we bridge this domain gap by efficiently transfer-learning from both domain-specific and task-specific source models. Additionally, we upgrade and expand an existing illustrated pose estimation dataset, and introduce two new datasets for classification and segmentation subtasks. We then apply the resultant state-of-the-art character pose estimator to solve the novel task of pose-guided illustration retrieval. All data, models, and code will be made publicly available.


title: Gravitational wave echoes search with combs

abstract: Gravitational wave echoes may provide a smoking gun signal for new physics in the immediate vicinity of black holes. As a quasiperiodic signal in time, echoes are characterized by the nearly constant time delay, and its precise measurement can help reveal a Planck-scale deviation outside of the would-be horizon. Different search methods have been developed for this quasiperiodic signal, while the searches suffer from large theoretical uncertainties of the echo waveform associated with the near-horizon physics. On the other hand, a coherent combination of a large number of pulses gives rise to a generic narrow resonance structure for the echo amplitude in frequency. The quasiperiodic resonance structure sets a complementary search target for echoes, and the time delay is inversely related to the average resonance spacing. A uniform comb has been proposed to look for the resonance structure in a rather model-independent way. In this paper, we develop a Bayesian algorithm to search for the resonance structure based on combs, where a phase-marginalized likelihood plays an essential role. The algorithm is validated with signal injections in detector noise from Advanced LIGO. With special treatments of the non-Gaussian artifacts, the noise outliers of the log Bayes factor distribution are properly removed. An echo signal not significantly below noise is detectable, and the time delay can be determined to very high precision. We perform the proposed search on real gravitational wave strain data of the first observing run of Advanced LIGO. We find no clear evidence of a comblike structure for GW150914 and GW151012.


title: Unsupervised Domain Adaptation for Retinal Vessel Segmentation with Adversarial Learning and Transfer Normalization

abstract: Retinal vessel segmentation plays a key role in computer-aided screening, diagnosis, and treatment of various cardiovascular and ophthalmic diseases. Recently, deep learning-based retinal vessel segmentation algorithms have achieved remarkable performance. However, due to the domain shift problem, the performance of these algorithms often degrades when they are applied to new data that is different from the training data. Manually labeling new data for each test domain is often a time-consuming and laborious task. In this work, we explore unsupervised domain adaptation in retinal vessel segmentation by using entropy-based adversarial learning and transfer normalization layer to train a segmentation network, which generalizes well across domains and requires no annotation of the target domain. Specifically, first, an entropy-based adversarial learning strategy is developed to reduce the distribution discrepancy between the source and target domains while also achieving the objective of entropy minimization on the target domain. In addition, a new transfer normalization layer is proposed to further boost the transferability of the deep network. It normalizes the features of each domain separately to compensate for the domain distribution gap. Besides, it also adaptively selects those feature channels that are more transferable between domains, thus further enhancing the generalization performance of the network. We conducted extensive experiments on three regular fundus image datasets and an ultra-widefield fundus image dataset, and the results show that our approach yields significant performance gains compared to other state-of-the-art methods.


title: M\"obius functions of higher rank and Dirichlet series

abstract: We introduce M\"obius functions of higher rank, a new class of arithmetic functions so that the classical M\"obius function is of rank 2. With this idea, we evaluate Dirichlet series on the sum of the reciprocal square of all $r$-free numbers. For the proof, the Riemann zeta function and cyclotomic polynomials play a key role.


title: Combining Attention with Flow for Person Image Synthesis

abstract: Pose-guided person image synthesis aims to synthesize person images by transforming reference images into target poses. In this paper, we observe that the commonly used spatial transformation blocks have complementary advantages. We propose a novel model by combining the attention operation with the flow-based operation. Our model not only takes the advantage of the attention operation to generate accurate target structures but also uses the flow-based operation to sample realistic source textures. Both objective and subjective experiments demonstrate the superiority of our model. Meanwhile, comprehensive ablation studies verify our hypotheses and show the efficacy of the proposed modules. Besides, additional experiments on the portrait image editing task demonstrate the versatility of the proposed combination.


title: Asymptotic stability of the superposition of viscous contact wave with rarefaction waves for the compressible Navier-Stokes-Maxwell equations

abstract: We study the large-time asymptotic behavior of solutions toward the combination of a viscous contact wave with two rarefaction waves for the compressible non-isentropic Navier-Stokes equations coupling with the Maxwell equations through the Lorentz force (called the Navier-Stokes-Maxwell equations). It includes the electrodynamic effects into the dissipative structure of the hyperbolic-parabolic system and turns out to be more complicated than that in the simpler compressible Navier-Stokes equations. Based on a new observation of the specific structure of the Maxwell equations in the Lagrangian coordinates, we prove that this typical composite wave pattern is time-asymptotically stable for the Navier-Stokes-Maxwell equations under some smallness conditions on the initial perturbations and wave strength, and also under the assumption that the dielectric constant is bounded. The main result is proved by using elementary energy methods. This is the first result about the nonlinear stability of the combination of two different wave patterns for the compressible Navier-Stokes-Maxwell equations.


title: Regret theory under fear of the unknown

abstract: It is common to encounter the situation with uncertainty for decision makers (DMs) in dealing with a complex decision making problem. The existing evidence shows that people usually fear the extreme uncertainty named as the unknown. This paper reports the modified version of the typical regret theory by considering the fear experienced by DMs for the unknown. Based on the responses of undergraduate students to the hypothetical choice problems with an unknown outcome, some experimental evidences are observed and analyzed. The framework of the modified regret theory is established by considering the effects of an unknown outcome. A fear function is equipped and some implications are proved. The behavioral foundation of the modified regret theory is further developed by modifying the axiomatic properties of the existing one as those based on the utility function; and it is recalled as the utility-based behavioral foundation for convenience. The application to the medical decision making with an unknown risk is studied and the effects of the fear function are investigated. The observations reveal that the existence of an unknown outcome could enhance, impede or reverse the preference relation of people in a choice problem, which can be predicted by the developed regret theory.


title: Effects of nonlocal dispersal strategies and heterogeneous environment on total population

abstract: In this paper, we consider the following single species model with nonlocal dispersal strategy $$ d\mathcal {L} [\theta] (x,t) + \theta(x,t) [m(x)- \theta(x,t)]=0 , $$ where $\mathcal {L}$ denotes the nonlocal diffusion operator, and investigate how the dispersal rate of the species and the distribution of resources affect the total population. First, we show that the upper bound for the ratio between total population and total resource is $C\sqrt{d}$. Moreover, examples are constructed to indicate that this upper bound is optimal. Secondly, for a type of simplified nonlocal diffusion operator, we prove that if $ \frac{\sup m}{\inf m}<\frac{\sqrt{5}+1}{2}$, the total population as a function of dispersal rate $d$ admits exactly one local maximum point in $\displaystyle (\inf m, \sup m)$. These results reveal essential discrepancies between local and nonlocal dispersal strategies.