title: Roots and Dynamics of Octonion Polynomials

abstract: This paper is devoted to several new results concerning (standard) octonion polynomials. The first is the determination of the roots of all right scalar multiples of octonion polynomials. The roots of left multiples are also discussed, especially over fields of characteristic not 2. We then turn to study the dynamics of monic quadratic real octonion polynomials, classifying the fixed points into attracting, repelling and ambivalent, and concluding with a discussion on the behavior of pseudo-periodic points.


title: Semi-modular forms from Fibonacci-Eisenstein series

abstract: In recent work, M. Just and the second author defined a class of "semi-modular forms" on $\mathbb C$, in analogy with classical modular forms, that are "half modular" in a particular sense; and constructed families of such functions as Eisenstein-like series using symmetries related to integer partitions. Looking for further natural examples of semi-modular behavior, here we construct a family of Eisenstein-like series to produce semi-modular forms, using symmetries related to Fibonacci numbers instead of partitions. We then consider other Lucas sequences that yield semi-modular forms.


title: Unconventional Weyl exceptional contours in non-Hermitian photonic continua

abstract: Unconventional Weyl points with topological charges higher than 1 can transform into various complex unconventional Weyl exceptional contours under non-Hermitian perturbations. However, theoretical studies of these exceptional contours have been limited to tight-binding models. Here, we propose to realize unconventional Weyl exceptional contours in photonic continua -- non-Hermitian anisotropic chiral plasma, based on ab initio calculation by Maxwell's equations. By perturbing in-plane permittivity, an unconventional Weyl point can transform into a quadratic Weyl exceptional circle, a Type-I Weyl exceptional chain with one chain point, a Type-II Weyl exceptional chain with two chain points, or other forms. Realistic metamaterials with effective constitutive parameters are proposed to implement these unconventional Weyl exceptional contours. Our work paves a way toward exploration of exotic physics of unconventional Weyl exceptional contours in non-Hermitian topological photonic continua.


title: Convergent evolution in a large cross-cultural database of musical scales

abstract: Scales, sets of discrete pitches that form the basis of melodies, are thought to be one of the most universal hallmarks of music. But we know relatively little about cross-cultural diversity of scales or how they evolved. To remedy this, we assemble a cross-cultural database (Database of Musical Scales: DaMuSc) of scale data, collected over the past century by various ethnomusicologists. Statistical analyses of the data highlight that certain intervals (e.g., the octave, fifth, second) are used frequently across cultures. Despite some diversity among scales, it is the similarities across societies which are most striking - most scales are found close to equidistant 5- and 7-note scales. We discuss the mechanisms of variation and selection in the evolution of scales, and how the assembled data may be used to examine the root causes of convergent evolution.


title: Favorable Interfacial Chemomechanics Enables Stable Cycling of High Li-Content Li-In/Sn Anodes in Sulfide Electrolyte Based Solid-State Batteries

abstract: Solid-state batteries (SSBs) can offer a paradigm shift in battery safety and energy density. Yet, the promise hinges on the ability to integrate high-performance electrodes with state-of-the-art solid electrolytes. For example, lithium (Li) metal, the most energy-dense anode candidate, suffers from severe interfacial chemomechanical issues that lead to cell failure. Li alloys of In/Sn are attractive alternatives, but their exploration has mostly been limited to the low capacity(low Li content)and In rich Li$_x$In (x$\leq$0.5). Here, the fundamental electro-chemo-mechanical behavior of Li-In and Li-Sn alloys of varied Li stoichiometries is unravelled in sulfide electrolyte based SSBs. The intermetallic electrodes developed through a controlled synthesis and fabrication technique display impressive (electro)chemical stability with Li$_6$PS$_5$Cl as the solid electrolyte and maintain nearly perfect interfacial contact during the electrochemical Li insertion/deinsertion under an optimal stack pressure. Their intriguing variation in the Li migration barrier with composition and its influence on the observed Li cycling overpotential is revealed through combined computational and electrochemical studies. Stable interfacial chemomechanics of the alloys allow long-term dendrite free Li cycling (>1000 h) at relatively high current densities (1 mA cm$^{-2}$) and capacities (1 mAh cm$^{-2}$), as demonstrated for Li$_{13}$In$_3$ and Li$_{17}$Sn$_4$, which are more desirable from a capacity and cost consideration compared to the low Li content analogues. The presented understanding can guide the development of high-capacity Li-In/Sn alloy anodes for SSBs.


title: Asymptotic expansions for the truncation error in Ramanujan-type series

abstract: Many of the fastest known algorithms to compute $\pi$ involve generalized hypergeometric series, such as the Ramanujan-Sato series. In this paper, we investigate the rates of convergence for several such series and we give asymptotic expansions for the error of finite approximation. For example, when using the first $n$ terms of the Chudnovskys' series, we obtain the finite approximation $\pi_n\approx \pi$. It is known that the truncation error satisfies $|\pi_n-\pi|\approx 53360^{-3n}.$ In this paper, we prove that the asymptotic expansion for the truncation error in the Chudnovskys' series is $$\left|\pi_n-\pi\right|=53360^{-3n}\cdot\frac{{106720}\sqrt{{10005}\pi}}{{1672209}\sqrt{n}}\cdot\exp\left(\frac{A_1}{n}+\frac{A_2}{n^2}+\frac{\delta_n}{n^3}\right),$$ with ${0.006907}<\delta_n<{0.008429}$ and the exact rational values of $A_1$ and $A_2$: $$A_1= -\frac{1781843197433}{7456754505816},$$ $$A_2= -\frac{1080096011925710088395}{3475199235000451148614116}.$$ Thus we demonstrate how to establish precise error bounds for the approximations for $\pi$ obtained through Ramanujan-like series for $1/\pi$. We also give asymptotic expansions for all known rational hypergeometric series for $1/\pi$ in the appendix.


title: Scoring a Goal optimally in a Soccer game under Liouville-like quantum gravity action

abstract: In this paper we present a new theoretical model to find out an optimal weight associated with a soccer player under the presence of a stochastic goal dynamics by using Feynman path integral method, where the action of a player is on $\sqrt{8/3}$-Liouville Quantum Gravity surface. Before determine the optimal weight we first establish an Infinitary logic which can deal with infinite variables on the strategy space then, a quantum formula of this logic has been developed and finally, based on this we are able show the existence of a Lefschetz-Hopf fixed point of this game. As in a competitive tournament, all possible standard strategies to score goals are known to the opposition team, a player's action is stochastic in nature and they would have some comparative advantage to score goals. Furthermore, conditions like uncertainties due to rain, dribbling and passing skill of a player, whether the match is a day match or a day-night match, advantages of having home crowd and asymmetric information of action profiles have been considered on the way to determine the optimal weight.


title: Massive feature extraction for explaining and foretelling hydroclimatic time series forecastability at the global scale

abstract: Statistical analyses and descriptive characterizations are sometimes assumed to be offering information on time series forecastability. Despite the scientific interest suggested by such assumptions, the relationships between descriptive time series features (e.g., temporal dependence, entropy, seasonality, trend and linearity features) and actual time series forecastability (quantified by issuing and assessing forecasts for the past) are scarcely studied and quantified in the literature. In this work, we aim to fill in this gap by investigating such relationships, and the way that they can be exploited for understanding hydroclimatic forecastability and its patterns. To this end, we follow a systematic framework bringing together a variety of -- mostly new for hydrology -- concepts and methods, including 57 descriptive features and nine seasonal time series forecasting methods (i.e., one simple, five exponential smoothing, two state space and one automated autoregressive fractionally integrated moving average methods). We apply this framework to three global datasets originating from the larger Global Historical Climatology Network (GHCN) and Global Streamflow Indices and Metadata (GSIM) archives. As these datasets comprise over 13 000 monthly temperature, precipitation and river flow time series from several continents and hydroclimatic regimes, they allow us to provide trustable characterizations and interpretations of 12-month ahead hydroclimatic forecastability at the global scale...


title: Large Deviations of Semi-supervised Learning in the Stochastic Block Model

abstract: In community detection on graphs, the semi-supervised learning problem entails inferring the ground-truth membership of each node in a graph, given the connectivity structure and a limited number of revealed node labels. Different subsets of revealed labels can in principle lead to higher or lower information gains and induce different reconstruction accuracies. In the framework of the dense stochastic block model, we employ statistical physics methods to derive a large deviation analysis for this problem, in the high-dimensional limit. This analysis allows the characterization of the fluctuations around the typical behaviour, capturing the effect of correlated label choices and yielding an estimate of their informativeness and their rareness among subsets of the same size. We find theoretical evidence of a non-monotonic relationship between reconstruction accuracy and the free energy associated to the posterior measure of the inference problem. We further discuss possible implications for active learning applications in community detection.


title: A generative model for age and income distribution

abstract: Each individual in society experiences an evolution of their income during their lifetime. Macroscopically, this dynamics creates a statistical relationship between age and income for each society. In this study, we investigate income distribution and its relationship with age and identify a stable joint distribution function for age and income within the United Kingdom and the United States. We demonstrate a flexible calibration methodology using panel and population surveys and capture the characteristic differences between the UK and the US populations. The model here presented can be utilised for forecasting income and planning pensions.


title: Use of analogies in science education, a systematic mapping study

abstract: This systematic mapping study consisted of tracking the scientific literature that addresses the issue of analogies as a didactic strategy in science teaching. An analogy can be understood as comparing an existing knowledge with a new knowledge to achieve a better understanding of the new knowledge as a result of the comparison of similarities; or in other words, use students' own concepts to introduce new concepts using comparisons between the two. The purpose of this study was to identify, analyze, synthesize and evaluate research works that touched on this topic, with this, to have knowledge about the models of uses of analogies, most used didactic strategies, research methodologies in this field and how to evaluate the learning effectiveness of working with analogies. The methodology that was used is the systematic mapping study; Five questions were posed that guided the information tracking process. Later, the electronic documents in English for the last twenty years were traced in five databases related to the educational field. Finally, it is concluded by responding to the purpose of the study where it is evident that, broadly speaking, the research methodologies in this field are quantitative as well as qualitative, to implement analogies, resources such as images, illustrations, textual indications and audiovisual aids are used, it is usually evaluated the effectiveness of using analogies with multiple choice tests, oral tests of creating analogies by students.


title: Testing the differentiated impact of the COVID-19 pandemic on air travel demand considering social inclusion

abstract: The economic downturn and the air travel crisis triggered by the recent coronavirus pandemic pose a substantial threat to the new consumer class of many emerging economies. In Brazil, considerable improvements in social inclusion have fostered the emergence of hundreds of thousands of first-time fliers over the past decades. We apply a two-step regression methodology in which the first step consists of identifying air transport markets characterized by greater social inclusion, using indicators of the local economies' income distribution, credit availability, and access to the Internet. In the second step, we inspect the drivers of the plunge in air travel demand since the pandemic began, differentiating markets by their predicted social inclusion intensity. After controlling for potential endogeneity stemming from the spread of COVID-19 through air travel, our results suggest that short and low-density routes are among the most impacted airline markets and that business-oriented routes are more impacted than leisure ones. Finally, we estimate that a market with 1 per cent higher social inclusion is associated with a 0.153 per cent to 0.166 per cent more pronounced decline in demand during the pandemic. Therefore, markets that have benefited from greater social inclusion in the country may be the most vulnerable to the current crisis.


title: Monte Carlo simulations as a route to compute probabilities

abstract: Monte Carlo simulations are based on the manipulation of random numbers to evaluate probable outcomes, with applicability in a variety of different fields. By assigning probabilities, which can be determined a priori, to various events, it is possible to track the evolution of the system over length and time scales which are not normally accessible to other simulation techniques. Monte Carlo simulations can provide insights, which can be used to develop more realistic models. In this work, these simulations are used to model a variety of probability problems normally encountered at the high school and undergraduate level. The simulations are used to introduce concepts related to system size, simulation runs (repeatability), and basic statistics. While many of the problems discussed here have analytical expressions, systems where easy analytical solutions are not available are also discussed. In this work, MATLAB is used to model these systems and the rand function in the software is used extensively for random number generation.


title: One-pot Liquid-Phase Synthesis of MoS$_2$-WS$_2$ van der Waals Heterostructures for Broadband Photodetection

abstract: Two dimensional (2D) van der Waals heterostructures (vdWHs) have their unique potential in facilitating the stacking of layers of different 2D materials for optoelectronic devices with superior characteristics at a reduced cost. However, the fabrication of large area all-2D heterostructures is still challenging towards realizing practical devices. In the present work, we have demonstrated a rapid yet simple, impurity free and highly efficient sonication-assisted chemical exfoliation approach to synthesize hybrid vdWHs based on 2D molybdenum disulphide (MoS$_2$) and tungsten disulphide (WS$_2$), with high yield. Microscopic and spectroscopic studies have confirmed the successful exfoliation of layered 2D materials and formation of their hybrid heterostructure. The co-existence of 2D MoS2 and WS2 in the vdW hybrid is established by optical absorption and Raman shift measurements along with their chemical stiochiometry determined by X-ray photoelectron spectroscopy. The spectral response of the vdWH/Si (2D/3D) heterojunction photodetector fabricated using the as-synthesized material is found to show superior broadband photoresponse compared to that shown by the individual 2D MoS$_2$ and WS$_2$ based devices. The peak responsivity is found to be ~2.15 A/W at a wavelength of ~560 nm for an applied bias of -5 V. The ease of fabrication and superior performance of the chemically synthesized vdWH-based devices have revealed their potential use for large area optoelectronic applications on Si compatible CMOS platforms.


title: Sea Ice Forecasting using Attention-based Ensemble LSTM

abstract: Accurately forecasting Arctic sea ice from subseasonal to seasonal scales has been a major scientific effort with fundamental challenges at play. In addition to physics-based earth system models, researchers have been applying multiple statistical and machine learning models for sea ice forecasting. Looking at the potential of data-driven sea ice forecasting, we propose an attention-based Long Short Term Memory (LSTM) ensemble method to predict monthly sea ice extent up to 1 month ahead. Using daily and monthly satellite retrieved sea ice data from NSIDC and atmospheric and oceanic variables from ERA5 reanalysis product for 39 years, we show that our multi-temporal ensemble method outperforms several baseline and recently proposed deep learning models. This will substantially improve our ability in predicting future Arctic sea ice changes, which is fundamental for forecasting transporting routes, resource development, coastal erosion, threats to Arctic coastal communities and wildlife.


title: An augmented Lagrangian algorithm for recovery of ice thickness in unidirectional flow using the Shallow Ice Approximation

abstract: A key parameter in ice flow modelling is basal slipping at the ice-bed interface as it can have a large effect on the resultant ice thickness. Unfortunately, its contribution to surface observations can be hard to distinguish from that of bed undulations. Therefore, inferring the ice thickness from surface measurements is an interesting and non-trivial inverse problem. This paper presents a method for recovering dually the ice thickness and the basal slip using only surface elevation and speed measurements. The unidirectional shallow ice approximation is first implemented to model steady state ice flow for given bedrock and basal slip profiles. This surface is then taken as synthetic observed data. An augmented Lagrangian algorithm is then used to find the diffusion coefficient which gives the best fit to observations. Combining this recovered diffusion with observed surface velocity, a simple Newton's method is used to recover both the ice thickness and basal slip. The method was successful in each test case and this implies that it should be possible to recover both of these parameters in two-dimensional cases also.

categories: physics.ao-ph math.OC physics.comp-ph physics.flu-dyn physics.geo-ph


title: An Efficient and Statistically Accurate Lagrangian Data Assimilation Algorithm with Applications to Discrete Element Sea Ice Models

abstract: Lagrangian data assimilation of complex nonlinear turbulent flows is an important but computationally challenging topic. In this article, an efficient data-driven statistically accurate reduced-order modeling algorithm is developed that significantly accelerates the computational efficiency of Lagrangian data assimilation. The algorithm starts with a Fourier transform of the high-dimensional flow field, which is followed by an effective model reduction that retains only a small subset of the Fourier coefficients corresponding to the energetic modes. Then a linear stochastic model is developed to approximate the nonlinear dynamics of each Fourier coefficient. Effective additive and multiplicative noise processes are incorporated to characterize the modes that exhibit Gaussian and non-Gaussian statistics, respectively. All the parameters in the reduced order system, including the multiplicative noise coefficients, are determined systematically via closed analytic formulae. These linear stochastic models succeed in forecasting the uncertainty and facilitate an extremely rapid data assimilation scheme. The new Lagrangian data assimilation is then applied to observations of sea ice floe trajectories that are driven by atmospheric winds and turbulent ocean currents. It is shown that observing only about $30$ non-interacting floes in a $200$km$\times200$km domain is sufficient to recover the key multi-scale features of the ocean currents. The additional observations of the floe angular displacements are found to be suitable supplements to the center-of-mass positions for improving the data assimilation skill. In addition, the observed large and small floes are more useful in recovering the large- and small-scale features of the ocean, respectively. The Fourier domain data assimilation also succeeds in recovering the ocean features in the areas where cloud cover obscures the observations.

categories: physics.ao-ph cs.NA math.NA physics.flu-dyn


title: A mechanism-based multi-trap phase field model for hydrogen assisted fracture

abstract: We present a new mechanistic, phase field-based formulation for predicting hydrogen embrittlement. The multi-physics model developed incorporates, for the first time, a Taylor-based dislocation model to resolve the mechanics of crack tip deformation. This enables capturing the role of dislocation hardening mechanisms in elevating the tensile stress, hydrogen concentration and dislocation trap density within tens of microns ahead of the crack tip. The constitutive strain gradient plasticity model employed is coupled to a phase field formulation, to simulate the fracture process, and to a multi-trap hydrogen transport model. The analysis of stationary and propagating cracks reveals that the modelling framework presented is capable of adequately capturing the sensitivity to the hydrogen concentration, the loading rate, the material strength and the plastic length scale. In addition, model predictions are compared to experimental data of notch tensile strength versus hydrogen content on a high-strength steel; a very good agreement is attained. We define and implement both atomistic-based and phenomenological hydrogen degradation laws and discuss similarities, differences and implications for the development of parameter-free hydrogen embrittlement models.


title: Recent advances in ADL, CutLang and adl2tnm

abstract: This paper presents an overview and features of an Analysis Description Language (ADL) designed for HEP data analysis. ADL is a domain specific, declarative language that describes the physics content of an analysis in a standard and unambiguous way, independent of any computing frameworks. It also describes infrastructures that render ADL executable, namely CutLang, a direct runtime interpreter (originally also a language), and adl2tnm, a transpiler converting ADL into C++ code. In ADL, analyses are described in human readable plain text files, clearly separating object, variable and event selection definitions in blocks, with a syntax that includes mathematical and logical operations, comparison and optimisation operators, reducers, four-vector algebra and commonly used functions. Recent studies demonstrate that adapting the ADL approach has numerous benefits for the experimental and phenomenological HEP communities. These include facilitating the abstraction, design, optimization, visualization, validation, combination, reproduction, interpretation and overall communication of the analysis contents and long term preservation of the analyses beyond the lifetimes of experiments. Here we also discuss some of the current ADL applications in physics studies and future prospects based on static analysis and differentiable programming.


title: Predictive and Prescriptive Performance of Bike-Sharing Demand Forecasts for Inventory Management

abstract: Bike-sharing systems are a rapidly developing mode of transportation and provide an efficient alternative to passive, motorized personal mobility. The asymmetric nature of bike demand causes the need for rebalancing bike stations, which is typically done during night time. To determine the optimal starting inventory level of a station for a given day, a User Dissatisfaction Function (UDF) models user pickups and returns as non-homogeneous Poisson processes with piece-wise linear rates. In this paper, we devise a deep generative model directly applicable in the UDF by introducing a variational Poisson recurrent neural network model (VP-RNN) to forecast future pickup and return rates. We empirically evaluate our approach against both traditional and learning-based forecasting methods on real trip travel data from the city of New York, USA, and show how our model outperforms benchmarks in terms of system efficiency and demand satisfaction. By explicitly focusing on the combination of decision-making algorithms with learning-based forecasting methods, we highlight a number of shortcomings in literature. Crucially, we show how more accurate predictions do not necessarily translate into better inventory decisions. By providing insights into the interplay between forecasts, model assumptions, and decisions, we point out that forecasts and decision models should be carefully evaluated and harmonized to optimally control shared mobility systems.


title: Spatio-temporal estimation of wind speed and wind power using machine learning: predictions, uncertainty and technical potential

abstract: The growth of wind generation capacities in the past decades has shown that wind energy can contribute to the energy transition in many parts of the world. Being highly variable and complex to model, the quantification of the spatio-temporal variation of wind power and the related uncertainty is highly relevant for energy planners. Machine Learning has become a popular tool to perform wind-speed and power predictions. However, the existing approaches have several limitations. These include (i) insufficient consideration of spatio-temporal correlations in wind-speed data, (ii) a lack of existing methodologies to quantify the uncertainty of wind speed prediction and its propagation to the wind-power estimation, and (iii) a focus on less than hourly frequencies. To overcome these limitations, we introduce a framework to reconstruct a spatio-temporal field on a regular grid from irregularly distributed wind-speed measurements. After decomposing data into temporally referenced basis functions and their corresponding spatially distributed coefficients, the latter are spatially modelled using Extreme Learning Machines. Estimates of both model and prediction uncertainties, and of their propagation after the transformation of wind speed into wind power, are then provided without any assumptions on distribution patterns of the data. The methodology is applied to the study of hourly wind power potential on a grid of 250 by 250 squared meters for turbines of 100 meters hub height in Switzerland, generating the first dataset of its type for the country. The potential wind power generation is combined with the available area for wind turbine installations to yield an estimate of the technical potential for wind power in Switzerland. The wind power estimate presented here represents an important input for planners to support the design of future energy systems with increased wind power generation.

categories: eess.SP cs.AI cs.LG physics.data-an stat.AP


title: U-GAT: Multimodal Graph Attention Network for COVID-19 Outcome Prediction

abstract: During the first wave of COVID-19, hospitals were overwhelmed with the high number of admitted patients. An accurate prediction of the most likely individual disease progression can improve the planning of limited resources and finding the optimal treatment for patients. However, when dealing with a newly emerging disease such as COVID-19, the impact of patient- and disease-specific factors (e.g. body weight or known co-morbidities) on the immediate course of disease is by and large unknown. In the case of COVID-19, the need for intensive care unit (ICU) admission of pneumonia patients is often determined only by acute indicators such as vital signs (e.g. breathing rate, blood oxygen levels), whereas statistical analysis and decision support systems that integrate all of the available data could enable an earlier prognosis. To this end, we propose a holistic graph-based approach combining both imaging and non-imaging information. Specifically, we introduce a multimodal similarity metric to build a population graph for clustering patients and an image-based end-to-end Graph Attention Network to process this graph and predict the COVID-19 patient outcomes: admission to ICU, need for ventilation and mortality. Additionally, the network segments chest CT images as an auxiliary task and extracts image features and radiomics for feature fusion with the available metadata. Results on a dataset collected in Klinikum rechts der Isar in Munich, Germany show that our approach outperforms single modality and non-graph baselines. Moreover, our clustering and graph attention allow for increased understanding of the patient relationships within the population graph and provide insight into the network's decision-making process.


title: Revisiting dissipative motion of a spinning heavy symmetric top and the rise of the top by friction

abstract: The dissipative motion and the rise of a heavy symmetrical top with a hemispherical peg are studied. A model taking the fixed point of the top as the center of the peg is considered when the top completely slips and the rolling motion is ignored. This is different from existing models like Jellet's one. Jellett's model and pure slipping are compared for different tops for the rise of the top.


title: On the dissipation at a shock wave in an elastic bar

abstract: This paper aims to quantitatively relate the energy dissipated at a shock wave in a nonlinearly elastic bar to the energy in the oscillations in two related dissipationless, dispersive systems. In contrast to a phase boundary, there is no kinetic relation associated with a shock wave. Three one-dimensional dynamic impact problems are studied: Problem 1 concerns a nonlinearly elastic bar, Problem 2 a discrete chain of particles, and Problem 3 a continuum with a strain gradient term in the constitutive relation. In the impact problem considered, the free boundary of each initially quiescent body is subjected to a sudden velocity that is then held constant for all subsequent time. There is energy dissipation at the shock in Problem 1 but Problems 2 and 3 are conservative. Problem 1 is solved analytically, Problem 2 numerically and an approximate solution to Problem 3 is constructed using modulation theory. The rate of increase of the oscillatory energy in Problems 2 and 3 are calculated and compared with the dissipation rate at the shock in Problem 1. The results indicate that the former is a good measure of the latter.


title: Nonlinear charge transport in highly polar semiconductors: GaN, AlN, InN and GaAs

abstract: In this paper, we present a collection of results focussing on the transport properties of doped direct-gap inverted-band highly polar III-nitride semiconductors (GaN, AlN, InN) and GaAs in the transient and steady state, calculated by using nonlinear quantum kinetic theory based on a non-equilibrium statistical ensemble formalism (NESEF). In the present paper, these results are compared with calculations usingMonteCarlo modelling simulations and experimental measurements. Both n-type and p-type materials, in the presence of intermediate to high electric fields, are considered for several temperatures and carrier concentrations.The agreement between the results obtained using nonlinear quantum kinetic theory, with those ofMonte Carlo calculations and experimental data is remarkably good, thus satisfactorily validating the NESEF.


title: Densely connected neural networks for nonlinear regression

abstract: Densely connected convolutional networks (DenseNet) behave well in image processing. However, for regression tasks, convolutional DenseNet may lose essential information from independent input features. To tackle this issue, we propose a novel DenseNet regression model where convolution and pooling layers are replaced by fully connected layers and the original concatenation shortcuts are maintained to reuse the feature. To investigate the effects of depth and input dimension of proposed model, careful validations are performed by extensive numerical simulation. The results give an optimal depth (19) and recommend a limited input dimension (under 200). Furthermore, compared with the baseline models including support vector regression, decision tree regression, and residual regression, our proposed model with the optimal depth performs best. Ultimately, DenseNet regression is applied to predict relative humidity, and the outcome shows a high correlation (0.91) with observations, which indicates that our model could advance environmental data analysis.

categories: cs.LG physics.ao-ph physics.data-an stat.ME


title: A SPA-based Manifold Learning Framework for Motor Imagery EEG Data Classification

abstract: The electroencephalography (EEG) signal is a non-stationary, stochastic, and highly non-linear bioelectric signal for which achieving high classification accuracy is challenging, especially when the number of subjects is limited. As frequently used solution, classifiers based on multilayer neural networks has to be implemented without large training data sets and careful tuning. This paper proposes a manifold learning framework to classify two types of EEG data from motor imagery (MI) tasks by discovering lower dimensional geometric structures. For feature extraction, it is implemented by Common Spatial Pattern (CSP) from the preprocessed EEG signals. In the neighborhoods of the features for classification, the local approximation to the support of the data is obtained, and then the features are assigned to the classes with the closest support. A spherical approximation (SPA) classifier is created using spherelets for local approximation, and the extracted features are classified with this manifold-based method. The SPA classifier achieves high accuracy in the 2008 BCI competition data, and the analysis shows that this method can significantly improve the decoding accuracy of MI tasks and exhibit strong robustness for small sample datasets. It would be simple and efficient to tune the two-parameters classifier for the online brain-computer interface(BCI)system.


title: Nonparametric posterior learning for emission tomography with multimodal data

abstract: We continue studies of the uncertainty quantification problem in emission tomographies such as PET or SPECT when additional multimodal data (e.g., anatomical MRI images) are available. To solve the aforementioned problem we adapt the recently proposed nonparametric posterior learning technique to the context of Poisson-type data in emission tomography. Using this approach we derive sampling algorithms which are trivially parallelizable, scalable and very easy to implement. In addition, we prove conditional consistency and tightness for the distribution of produced samples in the small noise limit (i.e., when the acquisition time tends to infinity) and derive new geometrical and necessary condition on how MRI images must be used. This condition arises naturally in the context of identifiability problem for misspecified generalized Poisson models. We also contrast our approach with Bayesian Markov Chain Monte Carlo sampling based on one data augmentation scheme which is very popular in the context of Expectation-Maximization algorithms for PET or SPECT. We show theoretically and also numerically that such data augmentation significantly increases mixing times for the Markov chain. In view of this, our algorithms seem to give a reasonable trade-off between design complexity, scalability, numerical load and assessment for the uncertainty.


title: Overview of the global semiconductor industry market

abstract: Mastering semiconductor technology is essential to insert any country into the trends of the future, such as smart cities, internet of things, space exploration, etc. In this paper we present the growing annual revenue of the semiconductor industry in the last 20 years and comment on the importance of mastering semiconductor production technology and its implications for the development of a nation.


title: On the practical application of negentropic polarization of thermodynamically equilibrium photon gas

abstract: In the course of the analysis of the axiomatic principles underlying statistical physics, the assumption about the limited correctness of the postulate of equiprobability of all available microstates of a closed system was tested. Earlier, the result of simulation modeling of the most probable stationary state of a stochastic isolated system, inside which isotropic monochrome radiation spontaneously acquired anisotropic polarization, was obtained (later the presence of such an effect was confirmed by a natural experiment). This article considers an already practically significant variant of suchlike a system. The conditions for the appearance of polarization anisotropy in a thermodynamically equilibrium medium initially filled with isotropic Planck radiation are predicted. Important consistent pattern of the applied mathematical model are noted.


title: The implications of gamma-ray photons from LHAASO on Lorentz symmetry

abstract: The Large High Altitude Air Shower Observatory (LHAASO) has reported the measurement of photons with high energy up to 1.42 PeV from twelve gamma-ray sources. We are concerned with the implications of LHAASO data on the fate of Lorenz symmetry at such high energy level, thus we consider the interaction of the gamma ray with those photons in cosmic microwave background (CMB), and compute the optical depth, the mean free path as well as the survival probability for photons from all these gamma-ray sources. Employing the threshold value predicted by the standard special relativity, it is found that the lowest survival probability for observed gamma ray photons is about 0.60, which is a fairly high value and implies that abundant photons with energy above the threshold value may reach the Earth without Lorentz symmetry violation. We conclude that it is still far to argue that the Lorentz symmetry would be violated due to the present observations from LHAASO.


title: Comparison of space weathering spectral changes induced by solar wind and micrometeoroid impacts using ion- and femtosecond-laser-irradiated olivine and pyroxene

abstract: Space weathering is a process that changes the surface of airless planetary bodies. Prime space weathering agents are solar wind irradiation and micrometeoroid bombardment. These processes alter planetary reflectance spectra and often modify their compositional diagnostic features. In this work we focused on simulating and comparing the spectral changes caused by solar wind irradiation and by micrometeoroid bombardment to gain a better understanding of these individual space weathering processes. We used olivine and pyroxene pellets as proxies for planetary materials. To simulate solar wind irradiation we used hydrogen, helium, and argon ions with energies from 5 to 40 keV and fluences of up to $10^{18}$ particles/cm$^2$. To simulate micrometeoroid bombardment we used individual femtosecond laser pulses. We analysed the corresponding evolution of different spectral parameters, which we determined by applying the Modified Gaussian Model, and we also conducted principal component analysis. The original mineralogy of the surface influences the spectral evolution more than the weathering agent, as seen from the diverse evolution of the spectral slope of olivine and pyroxene upon irradiation. The spectral slope changes seen in olivine are consistent with observations of A-type asteroids, while the moderate to no slope changes observed in pyroxene are consistent with asteroid (4) Vesta. We also observed some differences in the spectral effects induced by the two weathering agents. Ions simulating solar wind have a smaller influence on longer wavelengths of the spectra than laser irradiation simulating micrometeoroid impacts. This is most likely due to the different penetration depths of ions and laser pulses. Our results suggest that in some instances it might be possible to distinguish between the contributions of the two agents on a weathered surface.


title: Constrained Graphic Layout Generation via Latent Optimization

abstract: It is common in graphic design humans visually arrange various elements according to their design intent and semantics. For example, a title text almost always appears on top of other elements in a document. In this work, we generate graphic layouts that can flexibly incorporate such design semantics, either specified implicitly or explicitly by a user. We optimize using the latent space of an off-the-shelf layout generation model, allowing our approach to be complementary to and used with existing layout generation models. Our approach builds on a generative layout model based on a Transformer architecture, and formulates the layout generation as a constrained optimization problem where design constraints are used for element alignment, overlap avoidance, or any other user-specified relationship. We show in the experiments that our approach is capable of generating realistic layouts in both constrained and unconstrained generation tasks with a single model. The code is available at https://github.com/ktrk115/const_layout .


title: One-sided arc averaging geometries in time-distance local helioseismology

abstract: The study of solar oscillations (helioseismology) has been a very successful method of researching the Sun. Helioseismology teaches us about the structure and mean properties of the Sun. Together with mid-resolution data, the local properties were uncovered in quiet-Sun regions. However, magnetic fields affect the oscillations and prevent us from studying the properties of magnetically active regions with helioseismology. We aim to create a new methodology to suppress the negative effects of magnetic fields on solar oscillations and measure plasma properties close to active regions. The methodology consists of new averaging geometries, a non-linear approach to travel-time measurements, and a consistent inversion method that combines plasma flows and sound-speed perturbations. We constructed the one-sided arc averaging geometries and applied them to the non-linear approach of travel-time measurements. Using the one-sided arc travel times, we reconstructed the annulus travel times in a quiet-Sun region. We tested the methodology against the validated helioseismic inversion pipeline. We applied the new methodology for an inversion for surface horizontal flows in a region with a circular H-type sunspot. The inverted surface horizontal flows are comparable with the output of the coherent structure tracking, which is not strongly affected by the presence of the magnetic field. We show that the new methodology suppresses the negative effects of magnetic fields up to outer penumbra. We measure divergent flows with properties comparable to the moat flow. The new methodology can teach us about the depth structure of active regions and physical conditions that contribute to the evolution of the active regions.


title: Shallow Feature Matters for Weakly Supervised Object Localization

abstract: Weakly supervised object localization (WSOL) aims to localize objects by only utilizing image-level labels. Class activation maps (CAMs) are the commonly used features to achieve WSOL. However, previous CAM-based methods did not take full advantage of the shallow features, despite their importance for WSOL. Because shallow features are easily buried in background noise through conventional fusion. In this paper, we propose a simple but effective Shallow feature-aware Pseudo supervised Object Localization (SPOL) model for accurate WSOL, which makes the utmost of low-level features embedded in shallow layers. In practice, our SPOL model first generates the CAMs through a novel element-wise multiplication of shallow and deep feature maps, which filters the background noise and generates sharper boundaries robustly. Besides, we further propose a general class-agnostic segmentation model to achieve the accurate object mask, by only using the initial CAMs as the pseudo label without any extra annotation. Eventually, a bounding box extractor is applied to the object mask to locate the target. Experiments verify that our SPOL outperforms the state-of-the-art on both CUB-200 and ImageNet-1K benchmarks, achieving 93.44% and 67.15% (i.e., 3.93% and 2.13% improvement) Top-5 localization accuracy, respectively.


title: Few-Shot Domain Adaptation For End-to-End Communication

abstract: The problem of end-to-end learning of a communication system using an autoencoder -- consisting of an encoder, channel, and decoder modeled using neural networks -- has recently been shown to be an effective approach. A challenge faced in the practical adoption of this learning approach is that under changing channel conditions (e.g. a wireless link), it requires frequent retraining of the autoencoder in order to maintain a low decoding error rate. Since retraining is both time consuming and requires a large number of samples, it becomes impractical when the channel distribution is changing quickly. We propose to address this problem using a fast and sample-efficient (few-shot) domain adaptation method that does not change the encoder and decoder networks. Different from conventional training-time unsupervised or semi-supervised domain adaptation, here we have a trained autoencoder from a source distribution that we want to adapt (at test time) to a target distribution using only a small labeled dataset, and no unlabeled data. We focus on a generative channel model based on the Gaussian mixture density network (MDN), and propose a regularized, parameter-efficient adaptation of the MDN using a set of affine transformations. The learned affine transformations are then used to design an optimal transformation at the decoder input to compensate for the distribution shift, and effectively present to the decoder inputs close to the source distribution. Experiments on many simulated distribution changes common to the wireless setting, and a real mmWave FPGA testbed demonstrate the effectiveness of our method at adaptation using very few target domain samples. The code for our work can be found at: https://github.com/jayaram-r/domain-adaptation-autoencoder.


title: Intersection numbers on tropical Hassett spaces

abstract: We study the intersection of tropical psi-classes on tropical heavy/light Hassett spaces, generalising a result of Kerber--Markwig for tropical moduli spaces of rational stable curves with distinct marked points. Our computation reveals that the weight of a maximal cone in an intersection has a combinatorial intepretation in terms of the underlying tropical curve and it is always nonnegative. In particular, our result specialises to that, in top dimension, the tropical intersection product coincides with its classical counterpart.


title: The supercritical deformed Hermitian Yang--Mills equation on compact projective manifolds

abstract: In this paper, we extend a result of Gao Chen regarding the solvability of the twisted deformed Hermitian Yang-Mills equations on compact K\"ahler manifolds to allow for the twisting function to be non-constant and slightly negative in all dimensions. Using this result, we prove that the twisted supercritical dHYM equation on compact, projective manifolds can be solved provided certain numerical conditions are satisfied.


title: Optoelectronic and Electronic properties of cubic bulk and 001 surface of CsPbBr3

abstract: First Principles calculations are performed on the cubic bulk and its (001) surface of CsPbBr3 to reveal its electronic and optoelectronic properties. Various calculation schemes like GGA, GGA+SOC and the highly accurate modified Becke Johnson potential (mBJ) are applied. The electronic energy band gap calculations obtained from mBJ shows close approximity with the available experimental results. The cubic bulk phase is investigated for its mechanical stability and is found to bypass other structural phases. The optical absorption of both the bulk and the surface surpasses non-lead based cubic halides and zinc-based organic-inorganic hybrids.


title: Dust growth, fragmentation and self-induced dust traps in PHANTOM

abstract: We present the implementation of a dust growth and fragmentation module in the public Smoothed Particle Hydrodynamics (SPH) code PHANTOM. This module is made available for public use with this paper. The coagulation model considers locally monodisperse dust size distributions around single values that are carried by the SPH particles. Along with the presentation of the model, implementation and tests, we showcase growth and fragmentation in a few typical circumstellar disc simulations and revisit previous results. The module is also interfaced with the radiative transfer code MCFOST, which facilitates the comparison between simulations and ALMA observations by generating synthetic maps. Circumstellar disc simulations with growth and fragmentation reproduce the `self-induced dust trap' mechanism first proposed by Gonzalez et al., which supports its existence. Synthetic images of discs featuring this mechanism suggest it would be detectable by ALMA as a bright axisymmetric ring at several tens of au from the star. With this paper, our aim is to provide a public tool to be able to study and explore dust growth in a variety of applications related to planet formation.


title: Conveyor-mode single-electron shuttling in Si/SiGe for a scalable quantum computing architecture

abstract: Small spin-qubit registers defined by single electrons confined in Si/SiGe quantum dots operate successfully and connecting these would permit scalable quantum computation. Shuttling the qubit carrying electrons between registers is a natural choice for high-fidelity coherent links provided the overhead of control signals stays moderate. Our proof-of-principle demonstrates shuttling of a single electron by a propagating wave-potential in an electrostatically defined 420 nm long Si/SiGe quantum-channel. This conveyor-mode shuttling approach requires independent from its length only four sinusoidal control signals. We discuss the tuning of the signal parameters, detect the smoothness of the electron motion enabling the mapping of potential disorder and observe a high single-electron shuttling fidelity of $99.42\pm0.02\,\%$ including a reversal of direction. Our shuttling device can be readily embedded in industrial fabrication of Si/SiGe qubit chips and paves the way to solving the signal-fanout problem for a fully scalable semiconductor quantum-computing architecture.


title: Estimates for Interpolation Projectors and Related Problems in Computational Geometry

abstract: This paper contains a survey of results obtained by the authors mostly during the past few years and published by 2021. In particular, we present the best of known estimates of numerical characteristics related to the research theme. Sections: 1. Introduction. 2. The case when $n+1$ is an Hadamard number. 3. Estimates for the minimal absorption index of a cube by a simplex. 4. Estimates for the minimal norm of a projector in linear interpolation on a cube in ${\mathbb R}^n$. 5. Estimates of numbers $\xi_n^\prime$ and $\theta_n^\prime$. 6. Simplices satisfying the inclusions $S\subset Q_n\subset nS$. 7. Perfect simplices. 8. Equisecting simplices. 9. Properties of $(0,1)$-matrices of order $n$ having maximal determinant. 10. Problems for a simplex and a Euclidean ball. 11. Linear interpolation on a Euclidean ball. Bibliography: 56 titles. Keywords: simplex, cube, Euclidean ball, homothety, axial diameter, absorption index, Hadamard number, interpolation, projector, norm, estimate.


title: Small ball probability estimates for the H\"older semi-norm of the stochastic heat equation

abstract: We consider the stochastic heat equation on $[0,\,1]$ with periodic boundary conditions and driven by space-time white noise. Under various natural conditions, we study small ball probabilities for the H\"older semi-norms of the solutions, and provide near optimal bounds on these probabilities. As an application, we prove a support theorem in these H\"older semi-norms.


title: Shallow Attention Network for Polyp Segmentation

abstract: Accurate polyp segmentation is of great importance for colorectal cancer diagnosis. However, even with a powerful deep neural network, there still exists three big challenges that impede the development of polyp segmentation. (i) Samples collected under different conditions show inconsistent colors, causing the feature distribution gap and overfitting issue; (ii) Due to repeated feature downsampling, small polyps are easily degraded; (iii) Foreground and background pixels are imbalanced, leading to a biased training. To address the above issues, we propose the Shallow Attention Network (SANet) for polyp segmentation. Specifically, to eliminate the effects of color, we design the color exchange operation to decouple the image contents and colors, and force the model to focus more on the target shape and structure. Furthermore, to enhance the segmentation quality of small polyps, we propose the shallow attention module to filter out the background noise of shallow features. Thanks to the high resolution of shallow features, small polyps can be preserved correctly. In addition, to ease the severe pixel imbalance for small polyps, we propose a probability correction strategy (PCS) during the inference phase. Note that even though PCS is not involved in the training phase, it can still work well on a biased model and consistently improve the segmentation performance. Quantitative and qualitative experimental results on five challenging benchmarks confirm that our proposed SANet outperforms previous state-of-the-art methods by a large margin and achieves a speed about 72FPS.


title: Sequential Multivariate Change Detection with Calibrated and Memoryless False Detection Rates

abstract: Responding appropriately to the detections of a sequential change detector requires knowledge of the rate at which false positives occur in the absence of change. Setting detection thresholds to achieve a desired false positive rate is challenging. Existing works resort to setting time-invariant thresholds that focus on the expected runtime of the detector in the absence of change, either bounding it loosely from below or targeting it directly but with asymptotic arguments that we show cause significant miscalibration in practice. We present a simulation-based approach to setting time-varying thresholds that allows a desired expected runtime to be accurately targeted whilst additionally keeping the false positive rate constant across time steps. Whilst the approach to threshold setting is metric agnostic, we show how the cost of using the popular quadratic time MMD estimator can be reduced from $O(N^2B)$ to $O(N^2+NB)$ during configuration and from $O(N^2)$ to $O(N)$ during operation, where $N$ and $B$ are the numbers of reference and bootstrap samples respectively.


title: Lagrangian approach for the study of heat transfer in a nuclear reactor core using the SPH methodology

abstract: Numerical modeling simulations and the use of high-performance computing are fundamental for detailed safety analysis, control and operation of a nuclear reactor, allowing the study and analysis of problems related to thermal-hydraulics, neutronic, and the dynamic of fluids that are involved in these systems. In this work, we introduce the basis for the implementation of the smoothed particle hydrodynamics (SPH) approach to analyze heat transfer in a nuclear reactor core. Heat transfer by means of convection is of great importance in many engineering applications and especially in the analysis of heat transfer in nuclear reactors. As a first approach, the natural convection in the gap (space that exists between the fuel rod and the cladding) can be analyzed helping to reduce uncertainty in such calculations that usually relies on empirical correlations while using other numerical tools. The numerical method developed in this work was validated while comparing the results obtained in previous numerical simulations and experimental data reported in the literature showing that our implementation is suitable for the study of heat transfer in nuclear reactors. Numerical simulations were done with the DualSPHysics open source code that allows performing parallel calculations using different numbers of cores. The current implementation is a version written in CUDA (Compute Unified Device Architecture) that allows also the use of GPU processors (Graphics Processor Unit) to accelerate the calculations in parallel using a large number of cores contained in the GPU. This makes it possible to analyze large systems using a reasonable computer time. The obtained results verified and validated our method and allowed us to have a strong solver for future applications of heat transfer in nuclear reactors fuel inside the reactor cores.


title: Counterexample Classification

abstract: In model checking, when a given model fails to satisfy the desired specification, a typical model checker provides a counterexample that illustrates how the violation occurs. In general, there exist many diverse counterexamples that exhibit distinct violating behaviors, which the user may wish to examine before deciding how to repair the model. Unfortunately, obtaining this information is challenging in existing model checkers since (1) the number of counterexamples may be too large to enumerate one by one, and (2) many of these counterexamples are redundant, in that they describe the same type of violating behavior. In this paper, we propose a technique called counterexample classification. The goal of classification is to partition the space of all counterexamples into a finite set of counterexample classes, each of which describes a distinct type of violating behavior for the given specification. These classes are then presented as a summary of possible violating behaviors in the system, freeing the user from manually having to inspect or analyze numerous counterexamples to extract the same information. We have implemented a prototype of our technique on top of an existing formal modeling and verification tool, the Alloy Analyzer, and evaluated the effectiveness of the technique on case studies involving the well-known Needham-Schroeder protocol with promising results.


title: Cooperation dynamics under pandemic risks and heterogeneous economic interdependence

abstract: The spread of COVID-19 and ensuing containment measures have accentuated the profound interdependence among nations or regions. This has been particularly evident in tourism, one of the sectors most affected by uncoordinated mobility restrictions. The impact of this interdependence on the tendency to adopt less or more restrictive measures is hard to evaluate, more so if diversity in economic exposures to citizens' mobility are considered. Here, we address this problem by developing an analytical and computational game-theoretical model encompassing the conflicts arising from the need to control the economic effects of global risks, such as in the COVID-19 pandemic. The model includes the individual costs derived from severe restrictions imposed by governments, including the resulting economic interdependence among all the parties involved in the game. By using tourism-based data, the model is enriched with actual heterogeneous income losses, such that every player has a different economic cost when applying restrictions. We show that economic interdependence enhances cooperation because of the decline in the expected payoffs by free-riding parties (i.e., those neglecting the application of mobility restrictions). Furthermore, we show (analytically and through numerical simulations) that these cross-exposures can transform the nature of the cooperation dilemma each region or country faces, modifying the position of the fixed points and the size of the basins of attraction that characterize this class of games. Finally, our results suggest that heterogeneity among regions may be used to leverage the impact of intervention policies by ensuring an agreement among the most relevant initial set of cooperators.


title: A Machine learning approach for rapid disaster response based on multi-modal data. The case of housing & shelter needs

abstract: Along with climate change, more frequent extreme events, such as flooding and tropical cyclones, threaten the livelihoods and wellbeing of poor and vulnerable populations. One of the most immediate needs of people affected by a disaster is finding shelter. While the proliferation of data on disasters is already helping to save lives, identifying damages in buildings, assessing shelter needs, and finding appropriate places to establish emergency shelters or settlements require a wide range of data to be combined rapidly. To address this gap and make a headway in comprehensive assessments, this paper proposes a machine learning workflow that aims to fuse and rapidly analyse multimodal data. This workflow is built around open and online data to ensure scalability and broad accessibility. Based on a database of 19 characteristics for more than 200 disasters worldwide, a fusion approach at the decision level was used. This technique allows the collected multimodal data to share a common semantic space that facilitates the prediction of individual variables. Each fused numerical vector was fed into an unsupervised clustering algorithm called Self-Organizing-Maps (SOM). The trained SOM serves as a predictor for future cases, allowing predicting consequences such as total deaths, total people affected, and total damage, and provides specific recommendations for assessments in the shelter and housing sector. To achieve such prediction, a satellite image from before the disaster and the geographic and demographic conditions are shown to the trained model, which achieved a prediction accuracy of 62 %


title: Axion spectra and the associated x-ray spectra of low-mass stars

abstract: For the first time, we predict the axion spectra and their associated luminosities for several low-mass stars -- one and two solar masses stars in the main sequence and post-main sequence stages of evolution. Equally, we also compute the x-ray excess emission resulting from the conversion of axions back to photons in the presence of a strong magnetic field in the envelope of these stars. Hence, a given star will have a unique axion spectrum and $L_a$ axion luminosity. And if such star has a strong magnetic field in its stellar envelope, it will also show a characteristic x-ray spectrum and $L_{a\gamma}$ x-ray luminosity. Such radiation will add up to the x-ray electromagnetic spectrum and $L_X$ luminosity of the star. The present study focuses on axion models with an axion-photon coupling constant, $5\;10^{-11} {\rm GeV}^{-1}$, a value just below the most recent upper limit of $6.6\;10^{-11} {\rm GeV}^{-1}$ found by CAST and IAXO helioscopes. The range of axion parameters discussed here spans many axion models' parameter space, including the DFSZ and KSVZ models. We found that axions with a mass in the range $10^{-7}$ to $10^{-5}\; {\rm eV}$ and an axion-photon coupling constant of $5\;10^{-11} {\rm GeV}^{-1}$ produce an axion emission spectra with an averaged axion energy that varies from 1 to 5 KeV, and an $L_{a}$ ranging from $10^{-6} $ to $7\;10^{-3} \;L_\odot$. We also predict that $L_{a\gamma}$ varies from $5\;10^{-8}$ to $10^{-5} \;L_\odot$ for stars with an averaged magnetic field of $3 \;10^{4}\;{\rm G}$ in their atmospheres. Most of these $L_{a\gamma}$ predictions are larger than the $L_X$ observed in some stars. Therefore, such preliminary results show the potential of the next generation of stellar x-ray missions to constrain several classes of axion models.

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


title: Resilience of Well-structured Graph Transformation Systems

abstract: Resilience is a concept of rising interest in computer science and software engineering. For systems in which correctness w.r.t. a safety condition is unachievable, fast recovery is demanded. We investigate resilience problems of graph transformation systems. Our main contribution is the decidability of two resilience problems for well-structured graph transformation systems (with strong compatibility). We prove our results in the abstract framework of well-structured transition systems and apply them to graph transformation systems, incorporating also the concept of adverse conditions.


title: Make Scientific Conferences More Inclusive by Increasing Transparency in Speaker Selection

abstract: Scientific conferences play an important part in every scientists career. Here, I suggest an approach that can mitigate biases and help make scientific conferences more inclusive for everybody by increasing the transparency in speaker selection.


title: Rayleigh quotients of the level set manifolds related to the nonlinear PDE

abstract: The main topic of this note is a discussion of applicability conditions of the Nehari manifold method depending on the value of parameters of equations. As the main tool, we apply the nonlinear generalized Rayleigh quotient method.


title: Indexability and Rollout Policy for Multi-State Partially Observable Restless Bandits

abstract: Restless multi-armed bandits with partially observable states has applications in communication systems, age of information and recommendation systems. In this paper, we study multi-state partially observable restless bandit models. We consider three different models based on information observable to decision maker -- 1) no information is observable from actions of a bandit 2) perfect information from bandit is observable only for one action on bandit, there is a fixed restart state, i.e., transition occurs from all other states to that state 3) perfect state information is available to decision maker for both actions on a bandit and there are two restart state for two actions. We develop the structural properties. We also show a threshold type policy and indexability for model 2 and 3. We present Monte Carlo (MC) rollout policy. We use it for whittle index computation in case of model 2. We obtain the concentration bound on value function in terms of horizon length and number of trajectories for MC rollout policy. We derive explicit index formula for model 3. We finally describe Monte Carlo rollout policy for model 1 when it is difficult to show indexability. We demonstrate the numerical examples using myopic policy, Monte Carlo rollout policy and Whittle index policy. We observe that Monte Carlo rollout policy is good competitive policy to myopic.


title: Static analysis of ReLU neural networks with tropical polyhedra

abstract: This paper studies the problem of range analysis for feedforward neural networks, which is a basic primitive for applications such as robustness of neural networks, compliance to specifications and reachability analysis of neural-network feedback systems. Our approach focuses on ReLU (rectified linear unit) feedforward neural nets that present specific difficulties: approaches that exploit derivatives do not apply in general, the number of patterns of neuron activations can be quite large even for small networks, and convex approximations are generally too coarse. In this paper, we employ set-based methods and abstract interpretation that have been very successful in coping with similar difficulties in classical program verification. We present an approach that abstracts ReLU feedforward neural networks using tropical polyhedra. We show that tropical polyhedra can efficiently abstract ReLU activation function, while being able to control the loss of precision due to linear computations. We show how the connection between ReLU networks and tropical rational functions can provide approaches for range analysis of ReLU neural networks.


title: A Statistical Mechanics Approach to Describe Cell Re-orientation under Stretch

abstract: Experiments show that when a monolayer of cells cultured on an elastic substrate is subject to a cyclic stretch, cells tend to re-orient either perpendicularly or at an oblique angle with respect to the main direction of the stretch. Due to stochastic effects, however, the distribution of angles achieved by the cells is broader and, experimentally, histrograms over the interval [0, 90] are reported. Here we will determine the evolution and the stationary state of probability density functions describing the statistical distribution of the orientations of the cells using Fokker-Planck equations derived from microscopic rules for the evolution of the orientation of the cell. As a first attempt, we shall use a stochastic differential equation related to a very general elastic energy and we will show that the results of the time integration and of the stationary state of the related forward Fokker-Planck equation compare very well with experimental results obtained by different researchers. Then, in order to model more accurately the microscopic process of cell re-orientation, we consider discrete in time random processes that allow to recover Fokker- Planck equations through the well known technique of quasi-invariant limit. In particular, we shall introduce a non-local rule related to the evaluation of the state of stress experienced by the cell extending its protrusions, and a model of re-orientation as a result of an optimal control internally activated by the cell. Also in the latter case the results match very well with experiments.


title: Efficient Sparse Spherical k-Means for Document Clustering

abstract: Spherical k-Means is frequently used to cluster document collections because it performs reasonably well in many settings and is computationally efficient. However, the time complexity increases linearly with the number of clusters k, which limits the suitability of the algorithm for larger values of k depending on the size of the collection. Optimizations targeted at the Euclidean k-Means algorithm largely do not apply because the cosine distance is not a metric. We therefore propose an efficient indexing structure to improve the scalability of Spherical k-Means with respect to k. Our approach exploits the sparsity of the input vectors and the convergence behavior of k-Means to reduce the number of comparisons on each iteration significantly.


title: Gravitational influence of high power laser pulses

abstract: The study of the generation of metric perturbation in the laboratory presents an opportunity to observe and understand more easily the mechanisms at work in gravitation. The present study will focus on the metric perturbation generated by a light pulse, as it could be generated by a current ultra-high power laser. Although of very small magnitude, the potential thus generated has advantages over that generated by mass acceleration, such as the absence of noise due to non uniform acceleration or the ability to scale up the experiment. It is indeed easier to scale up an electromagnetic oscillation compared to a mechanical oscillator, which must either be made with a large accelerated mass or a lot of small masses, all in sync, which acceleration must furthermore be quadripolar. Generation of metric deformation by laser could therefore prove useful in the long-term establishment of a laboratory experiment for the generation and detection of gravitational waves.


title: PML and high-accuracy boundary integral equation solver for wave scattering by a locally defected periodic surface

abstract: This paper studies the PML method for wave scattering in a half space of homogeneous medium bounded by a two-dimensional, perfectly conducting, and locally defected periodic surface, and develops a high-accuracy boundary-integral-equation (BIE) solver. Along the vertical direction, we place a PML to truncate the unbounded domain onto a strip and prove that the PML solution converges linearly to the true solution in the physical subregion of the strip with the PML thickness. Laterally, we divide the unbounded strip into three regions: a region containing the defect and two semi-waveguide regions, separated by two vertical line segments. In both semi-waveguides, we prove the well-posedness of an associated scattering problem so as to well define a Neumann-to-Dirichlet (NtD) operator on the associated vertical segment. The two NtD operators, serving as exact lateral boundary conditions, reformulate the unbounded strip problem as a boundary value problem onto the defected region. Due to the periodicity of the semi-waveguides, both NtD operators turn out to be closely related to a Neumann-marching operator, governed by a nonlinear Riccati equation. It is proved that the Neumann-marching operators are contracting, so that the PML solution decays exponentially fast along both lateral directions. The consequences culminate in two opposite aspects. Negatively, the PML solution cannot exponentially converge to the true solution in the whole physical region of the strip. Positively, from a numerical perspective, the Riccati equations can now be efficiently solved by a recursive doubling procedure and a high-accuracy PML-based BIE method so that the boundary value problem on the defected region can be solved efficiently and accurately. Numerical experiments demonstrate that the PML solution converges exponentially fast to the true solution in any compact subdomain of the strip.

categories: math.NA cs.NA physics.comp-ph physics.optics


title: Entanglement in descendants

abstract: We study the single interval entanglement and relative entropies of conformal descendants in 2d CFT. Descendants contain non-trivial entanglement, though the entanglement entropy of the canonical primary in the free boson CFT contains no additional entanglement compared to the vacuum, we show that the entanglement entropy of the state created by its level one descendant is non-trivial and is identical to that of the $U(1)$ current in this theory. We determine the first sub-leading corrections to the short interval expansion of the entanglement entropy of descendants in a general CFT from their four point function on the n-sheeted plane. We show that these corrections are determined by multiplying squares of appropriate dressing factors to the corresponding corrections of the primary. Relative entropy between descendants of the same primary is proportional to the square of the difference of their dressing factors. We apply our results to a class of descendants of generalized free fields and descendants of the vacuum and show that their dressing factors are universal.


title: Adversarial Data Augmentation for Disordered Speech Recognition

abstract: Automatic recognition of disordered speech remains a highly challenging task to date. The underlying neuro-motor conditions, often compounded with co-occurring physical disabilities, lead to the difficulty in collecting large quantities of impaired speech required for ASR system development. To this end, data augmentation techniques play a vital role in current disordered speech recognition systems. In contrast to existing data augmentation techniques only modifying the speaking rate or overall shape of spectral contour, fine-grained spectro-temporal differences between disordered and normal speech are modelled using deep convolutional generative adversarial networks (DCGAN) during data augmentation to modify normal speech spectra into those closer to disordered speech. Experiments conducted on the UASpeech corpus suggest the proposed adversarial data augmentation approach consistently outperformed the baseline augmentation methods using tempo or speed perturbation on a state-of-the-art hybrid DNN system. An overall word error rate (WER) reduction up to 3.05\% (9.7\% relative) was obtained over the baseline system using no data augmentation. The final learning hidden unit contribution (LHUC) speaker adapted system using the best adversarial augmentation approach gives an overall WER of 25.89% on the UASpeech test set of 16 dysarthric speakers.


title: Groups as automorphisms of dessins d'enfants

abstract: It is known that every finite group can be represented as the full group of automorphisms of a suitable compact dessin d'enfant. In this paper, we give a constructive and easy proof that the same holds for any countable group by considering non-compact dessins. Moreover, we show that any tame action of a countable group is so realizable.


title: Multiphase flows simulation with the Smoothed Particle Hydrodynamics method

abstract: This work presents a new multiphase SPH model that includes the shifting algorithm and a variable smoothing length formalism to simulate multi-phase flows with accuracy and proper interphase management. The implementation was performed in the DualSPHysics code and validated for different canonical experiments, such as the single-phase and multiphase Poiseuille and Couette test cases. The method is accurate even for the multiphase case for which two phases are simulated. The shifting algorithm and the variable smoothing length formalism has been applied in the multiphase SPH model to improve the numerical results at the interphase even when it is highly deformed and non-linear effects become important. The obtained accuracy in the validation tests and the good interphase definition in the instability cases indicate an important improvement in the numerical results compared with single-phase and multiphase models where the shifting algorithm and the variable smoothing length formalism are not applied.


title: Chiral anomaly, induced current, and vacuum polarization tensor for a Dirac field in the presence of a defect

abstract: We evaluate the vacuum polarization tensor (VPT) for a massless Dirac field in 1+1 and 3+1 dimensions, in the presence of a particular kind of defect, which in a special limit imposes bag boundary conditions. We also show that the chiral anomaly in the presence of such a defect is the same as when no defects are present, both in 1+1 and 3+1 dimensions. This implies that the induced vacuum current in 1+1 dimensions due to the lowest order VPT is exact.


title: Extending Sticky-Datalog+/- via Finite-Position Selection Functions: Tractability, Algorithms, and Optimization

abstract: Weakly-Sticky(WS) Datalog+/- is an expressive member of the family of Datalog+/- program classes that is defined on the basis of the conditions of stickiness and weak-acyclicity. Conjunctive query answering (QA) over the WS programs has been investigated, and its tractability in data complexity has been established. However, the design and implementation of practical QA algorithms and their optimizations have been open. In order to fill this gap, we first study Sticky and WS programs from the point of view of the behavior of the chase procedure. We extend the stickiness property of the chase to that of generalized stickiness of the chase (GSCh) modulo an oracle that selects (and provides) the predicate positions where finitely values appear during the chase. Stickiness modulo a selection function S that provides only a subset of those positions defines sch(S), a semantic subclass of GSCh. Program classes with selection functions include Sticky and WS, and another syntactic class that we introduce and characterize, namely JWS, of jointly-weakly-sticky programs, which contains WS. The selection functions for these last three classes are computable, and no external, possibly non-computable oracle is needed. We propose a bottom-up QA algorithm for programs in the class sch(S), for a general selection function S. As a particular case, we obtain a polynomial-time QA algorithm for JWS and weakly-sticky programs. Unlike WS, JWS turns out to be closed under magic-sets query optimization. As a consequence, both the generic polynomial-time QA algorithm and its magic-set optimization can be particularized and applied to WS.


title: How Morphological Computation shapes Integrated Information in Embodied Agents

abstract: The Integrated Information Theory provides a quantitative approach to consciousness and can be applied to neural networks. An embodied agent controlled by such a network influences and is being influenced by its environment. This involves, on the one hand, morphological computation within goal directed action and, on the other hand, integrated information within the controller, the agent's brain. In this article, we combine different methods in order to examine the information flows among and within the body, the brain and the environment of an agent. This allows us to relate various information flows to each other. We test this framework in a simple experimental setup. There, we calculate the optimal policy for goal-directed behavior based on the "planning as inference" method, in which the information-geometric em-algorithm is used to optimize the likelihood of the goal. Morphological computation and integrated information are then calculated with respect to the optimal policies. Comparing the dynamics of these measures under changing morphological circumstances highlights the antagonistic relationship between these two concepts. The more morphological computation is involved, the less information integration within the brain is required. In order to determine the influence of the brain on the behavior of the agent it is necessary to additionally measure the information flow to and from the brain.


title: Spike timing regularity in vestibular afferent neurons: How ionic currents influence sensory encoding mechanisms

abstract: Primary vestibular neurons are categorized as either regularly or irregularly firing afferents that use rate and temporal sensory encoding strategies, respectively. While many factors influence firing in these neurons, recent work in mammalian vestibular afferents has demonstrated a rich diversity in ion channels that drive spiking regularity. Here, I review key ionic currents studied in vitro and demonstrate how they may enable sensory encoding strategies demonstrated in vivo.


title: Tree-Algorithms with Multi-Packet Reception and Successive Interference Cancellation

abstract: In this paper, we perform a thorough analysis of tree-algorithms with multi-packet reception (MPR) and successive interference cancellation (SIC), showing a number of novel results. We first derive the basic performance parameters, which are the expected length of the collision resolution interval and the normalized throughput, conditioned on the number of contending users. We then study their asymptotic behaviour, identifying an oscillatory component that amplifies with the increase in MPR. In the next step, we derive the throughput for the gated and windowed access, assuming Poisson arrivals. We show that for windowed access, the bound on maximum stable normalized throughput increases with the increase in MPR. his implies that investing in advanced physical capabilities, i.e., MPR and SIC pays off from the perspective of the medium access control algorithm.


title: The Substructures in Disks undergoing Vertical Shear Instability: II. Observational Predictions for the Dust Continuum

abstract: High-angular resolution observations at sub-millimeter/millimeter wavelengths of disks surrounding young stars have shown that their morphology is made of azimuthally-symmetric or point-symmetric substructures, in some cases with spiral arms, localized spur- or crescent-shaped features. The majority of theoretical studies with the aim of interpreting the observational results have focused on disk models with planets, under the assumption that the disk substructures are due to the disk-planet interaction. However, so far only in very few cases exoplanets have been detected in these systems. Furthermore, some substructures are expected to appear \textit{before} planets form, as they are necessary to drive the concentration of small solids which can lead to the formation of planetesimals. In this work we present observational predictions from high-resolution 3D radiative hydrodynamical models which follow the evolution of gas and solids in a protoplanetary disk. We focus on substructures in the distribution of millimeter-sized and smaller solid particles produced by the vertical shear instability. We show that their characteristics are compatible with some of the shallow gaps detected in recent observations at sub-mm/mm wavelengths, and present predictions for future observations with better sensitivity and angular resolution with ALMA and a Next Generation Very Large Array.


title: Disappearance of partonic collectivity in $\sqrt{s_{NN}}$ = 3 GeV Au+Au collisions at RHIC

abstract: We report on the measurements of directed flow $v_1$ and elliptic flow $v_2$ for hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{S}^0$, $p$, $\phi$, $\Lambda$ and $\Xi^{-}$) from Au+Au collisions at $\sqrt{s_{NN}}$ = 3\,GeV and $v_{2}$ for ($\pi^{\pm}$, $K^{\pm}$, $p$ and $\overline{p}$) at 27 and 54.4\,GeV with the STAR experiment. While at the two higher energy midcentral collisions the number-of-constituent-quark (NCQ) scaling holds, at 3\,GeV the $v_{2}$ at midrapidity is negative for all hadrons and the NCQ scaling is absent. In addition, the $v_1$ slopes at midrapidity for almost all observed hadrons are found to be positive, implying dominant repulsive baryonic interactions. The features of negative $v_2$ and positive $v_1$ slope at 3\,GeV can be reproduced with a baryonic mean-field in transport model calculations. These results imply that the medium in such collisions is likely characterized by baryonic interactions.


title: Spin Wave Computing using pre-recorded magnetization patterns

abstract: We propose a novel type of a spin wave computing device, based on a bilayer structure which includes a bias layer, made from a hard magnetic material and a propagation layer, made from a magnetic material with low damping, for example, Yttrium Garnet (YiG) or Permalloy. The bias layer maintains a stable pre-recorded magnetization pattern, generating a bias field with a desired spatial dependence, which in turn sets the equilibrium magnetization inside the propagation layer. When an external source applies an RF field or spinwave to the propagation layer, excited spin waves scatter on the magnetization's inhomogenuities, resulting in a complex interference behavior. One thus has the ability to adjust spin wave propagation properties simply by altering the magnetization in the bias layer. We demonstrate that the phenomenon can be utilized to perform a variety of computational operations, including Fourier Transform, Vector-Matrix multiplication and Grover search algorithm, with the operational parameters exceeding conventional designs by orders of magnitude.


title: Birth and decline of magma oceans. Part 1: erosion and deposition of crystal layers in evolving magmatic reservoirs

abstract: This paper is the first of a two companion papers presenting a theoretical and experimental study of the evolution of crystallizing magma oceans in planetesimals. We aim to understand the behavior of crystals formed in a convective magma and the implications of crystal segregation for the reservoir thermal and structural evolution. In particular, the goal is to constrain the possibility to form and preserve cumulates and/or flotation crusts by sedimentation/flotation of crystals. We first use lab-scale analog experiments to study the stability and the erosion of a floating lid composed of plastics beads over a convective viscous fluid volumetrically heated by microwave absorption. We propose an erosion law that depends only on two dimensionless numbers which govern these phenomena: (i) the Rayleigh-Roberts number, characterizing the strength of convection and (ii) the Shields number, that encompasses the physics of the flow-particle interaction. We further consider the formation of a cumulate at the base of the convective layer due to sedimentation of beads that are denser than the fluid. We find that particles deposition occurs at a velocity that scales with the Stokes velocity, a result consistent with previous experimental studies. The theoretical framework built on these experimental results is applied in a second paper on the evolution of magma oceans in planetesimals and the fate of particles in this convective environment.


title: Multi-phase Liver Tumor Segmentation with Spatial Aggregation and Uncertain Region Inpainting

abstract: Multi-phase computed tomography (CT) images provide crucial complementary information for accurate liver tumor segmentation (LiTS). State-of-the-art multi-phase LiTS methods usually fused cross-phase features through phase-weighted summation or channel-attention based concatenation. However, these methods ignored the spatial (pixel-wise) relationships between different phases, hence leading to insufficient feature integration. In addition, the performance of existing methods remains subject to the uncertainty in segmentation, which is particularly acute in tumor boundary regions. In this work, we propose a novel LiTS method to adequately aggregate multi-phase information and refine uncertain region segmentation. To this end, we introduce a spatial aggregation module (SAM), which encourages per-pixel interactions between different phases, to make full use of cross-phase information. Moreover, we devise an uncertain region inpainting module (URIM) to refine uncertain pixels using neighboring discriminative features. Experiments on an in-house multi-phase CT dataset of focal liver lesions (MPCT-FLLs) demonstrate that our method achieves promising liver tumor segmentation and outperforms state-of-the-arts.


title: Robust Acoustic Scene Classification in the Presence of Active Foreground Speech

abstract: We present an iVector based Acoustic Scene Classification (ASC) system suited for real life settings where active foreground speech can be present. In the proposed system, each recording is represented by a fixed-length iVector that models the recording's important properties. A regularized Gaussian backend classifier with class-specific covariance models is used to extract the relevant acoustic scene information from these iVectors. To alleviate the large performance degradation when a foreground speaker dominates the captured signal, we investigate the use of the iVector framework on Mel-Frequency Cepstral Coefficients (MFCCs) that are derived from an estimate of the noise power spectral density. This noise-floor can be extracted in a statistical manner for single channel recordings. We show that the use of noise-floor features is complementary to multi-condition training in which foreground speech is added to training signal to reduce the mismatch between training and testing conditions. Experimental results on the DCASE 2016 Task 1 dataset show that the noise-floor based features and multi-condition training realize significant classification accuracy gains of up to more than 25 percentage points (absolute) in the most adverse conditions. These promising results can further facilitate the integration of ASC in resource-constrained devices such as hearables.


title: I2V-GAN: Unpaired Infrared-to-Visible Video Translation

abstract: Human vision is often adversely affected by complex environmental factors, especially in night vision scenarios. Thus, infrared cameras are often leveraged to help enhance the visual effects via detecting infrared radiation in the surrounding environment, but the infrared videos are undesirable due to the lack of detailed semantic information. In such a case, an effective video-to-video translation method from the infrared domain to the visible light counterpart is strongly needed by overcoming the intrinsic huge gap between infrared and visible fields. To address this challenging problem, we propose an infrared-to-visible (I2V) video translation method I2V-GAN to generate fine-grained and spatial-temporal consistent visible light videos by given unpaired infrared videos. Technically, our model capitalizes on three types of constraints: 1)adversarial constraint to generate synthetic frames that are similar to the real ones, 2)cyclic consistency with the introduced perceptual loss for effective content conversion as well as style preservation, and 3)similarity constraints across and within domains to enhance the content and motion consistency in both spatial and temporal spaces at a fine-grained level. Furthermore, the current public available infrared and visible light datasets are mainly used for object detection or tracking, and some are composed of discontinuous images which are not suitable for video tasks. Thus, we provide a new dataset for I2V video translation, which is named IRVI. Specifically, it has 12 consecutive video clips of vehicle and monitoring scenes, and both infrared and visible light videos could be apart into 24352 frames. Comprehensive experiments validate that I2V-GAN is superior to the compared SOTA methods in the translation of I2V videos with higher fluency and finer semantic details. The code and IRVI dataset are available at https://github.com/BIT-DA/I2V-GAN.


title: Hardness and Approximation of Submodular Minimum Linear Ordering Problems

abstract: The minimum linear ordering problem (MLOP) generalizes well-known combinatorial optimization problems such as minimum linear arrangement and minimum sum set cover. MLOP seeks to minimize an aggregated cost $f(\cdot)$ due to an ordering $\sigma$ of the items (say $[n]$), i.e., $\min_{\sigma} \sum_{i\in [n]} f(E_{i,\sigma})$, where $E_{i,\sigma}$ is the set of items mapped by $\sigma$ to indices $[i]$. Despite an extensive literature on MLOP variants and approximations for these, it was unclear whether the graphic matroid MLOP was NP-hard. We settle this question through non-trivial reductions from mininimum latency vertex cover and minimum sum vertex cover problems. We further propose a new combinatorial algorithm for approximating monotone submodular MLOP, using the theory of principal partitions. This is in contrast to the rounding algorithm by Iwata, Tetali, and Tripathi [ITT2012], using Lov\'asz extension of submodular functions. We show a $(2-\frac{1+\ell_{f}}{1+|E|})$-approximation for monotone submodular MLOP where $\ell_{f}=\frac{f(E)}{\max_{x\in E}f(\{x\})}$ satisfies $1 \leq \ell_f \leq |E|$. Our theory provides new approximation bounds for special cases of the problem, in particular a $(2-\frac{1+r(E)}{1+|E|})$-approximation for the matroid MLOP, where $f = r$ is the rank function of a matroid. We further show that minimum latency vertex cover (MLVC) is $\frac{4}{3}$-approximable, by which we also lower bound the integrality gap of its natural LP relaxation, which might be of independent interest.


title: On sharp scattering threshold for the mass-energy double critical NLS via double track profile decomposition

abstract: The present paper is concerned with the large data scattering problem for the mass-energy double critical NLS \begin{align} i\partial_t u+\Delta u\pm |u|^{\frac{4}{d}}u\pm |u|^{\frac{4}{d-2}}u=0\tag{DCNLS} \end{align} in $H^1(\mathbb{R}^d)$ with $d\geq 3$. In the defocusing-defocusing regime, Tao, Visan and Zhang show that the unique solution of DCNLS is global and scattering in time for arbitrary initial data in $H^1(\mathbb{R}^d)$. This does not hold when at least one of the nonlinearities is focusing, due to the possible formation of blow-up and soliton solutions. However, precise thresholds for a solution of DCNLS being scattering were open in all the remaining regimes. Following the classical concentration compactness principle, we impose sharp scattering thresholds in terms of ground states for DCNLS in all the remaining regimes. The new challenge arises from the fact that the remainders of the standard $L^2$- or $\dot{H}^1$-profile decomposition fail to have asymptotically vanishing diagonal $L^2$- and $\dot{H}^1$-Strichartz norms simultaneously. To overcome this difficulty, we construct a double track profile decomposition which is capable to capture the low, medium and high frequency bubbles within a single profile decomposition and possesses remainders that are asymptotically small in both of the diagonal $L^2$- and $\dot{H}^1$-Strichartz spaces.


title: 2-D Directed Formation Control Based on Bipolar Coordinates

abstract: This work proposes a novel 2-D formation control scheme for acyclic triangulated directed graphs (a class of minimally acyclic persistent graphs) based on bipolar coordinates with (almost) global convergence to the desired shape. Prescribed performance control is employed to devise a decentralized control law that avoids singularities and introduces robustness against external disturbances while ensuring predefined transient and steady-state performance for the closed-loop system. Furthermore, it is shown that the proposed formation control scheme can handle formation maneuvering, scaling, and orientation specifications simultaneously. Additionally, the proposed control law is implementable in agents' arbitrarily oriented local coordinate frames using only low-cost onboard vision sensors, which are favorable for practical applications. Finally, a formation maneuvering simulation study verifies the proposed approach.