title: On the hybrid Davies like generator for quantum dissipation

abstract: We provide a class of quantum evolution beyond Markovian semigroup. This class is governed by a hybrid Davies like generator such that dissipation is controlled by a suitable memory kernel and decoherence by standard GKLS generator. These two processes commute and both of them commute with the unitary evolution controlled by the systems Hamiltonian. The corresponding memory kernel gives rise to semi-Markov evolution of the diagonal elements of the density matrix. However, the corresponding evolution needs not be completely positive. The role of decoherence generator is to restore complete positivity. Hence, to pose the dynamical problem one needs two processes generated by classical semi-Markov memory kernel and purely quantum decoherence generator. This scheme is illustrated for a qubit evolution.


title: Gravitational Redshift Tests with Atomic Clocks and Atom Interferometers

abstract: Atomic interference experiments can probe the gravitational redshift via the internal energy splitting of atoms and thus give direct access to test the universality of the coupling between matter-energy and gravity at different spacetime points. By including possible violations of the equivalence principle in a fully quantized treatment of all degrees of freedom, we characterize how the sensitivity to gravitational redshift violations arises in atomic clocks and atom interferometers, as well as their underlying limitations. Specifically, we show that: (i.) Contributions beyond linear order to trapping potentials lead to such a sensitivity of trapped atomic clocks. (ii.) While Bragg-type interferometers, even with a superposition of internal states, with state-independent, linear interaction potentials are at first insensitive to gravitational redshift tests, modified configurations, for example by relaunching the atoms, can mimic such tests tests under certain conditions. (iii.) Guided atom interferometers are comparable to atomic clocks. (iv.) Internal transitions lead to state-dependent interaction potentials through which light-pulse atom interferometers can become sensitive to gravitational redshift violations.


title: COSCO: Container Orchestration using Co-Simulation and Gradient Based Optimization for Fog Computing Environments

abstract: Intelligent task placement and management of tasks in large-scale fog platforms is challenging due to the highly volatile nature of modern workload applications and sensitive user requirements of low energy consumption and response time. Container orchestration platforms have emerged to alleviate this problem with prior art either using heuristics to quickly reach scheduling decisions or AI driven methods like reinforcement learning and evolutionary approaches to adapt to dynamic scenarios. The former often fail to quickly adapt in highly dynamic environments, whereas the latter have run-times that are slow enough to negatively impact response time. Therefore, there is a need for scheduling policies that are both reactive to work efficiently in volatile environments and have low scheduling overheads. To achieve this, we propose a Gradient Based Optimization Strategy using Back-propagation of gradients with respect to Input (GOBI). Further, we leverage the accuracy of predictive digital-twin models and simulation capabilities by developing a Coupled Simulation and Container Orchestration Framework (COSCO). Using this, we create a hybrid simulation driven decision approach, GOBI*, to optimize Quality of Service (QoS) parameters. Co-simulation and the back-propagation approaches allow these methods to adapt quickly in volatile environments. Experiments conducted using real-world data on fog applications using the GOBI and GOBI* methods, show a significant improvement in terms of energy consumption, response time, Service Level Objective and scheduling time by up to 15, 40, 4, and 82 percent respectively when compared to the state-of-the-art algorithms.


title: Enhanced weak-value amplification via photon recycling

abstract: In a quantum-noise limited system, weak-value amplification using post-selection normally does not produce more sensitive measurements than standard methods for ideal detectors: the increased weak value is compensated by the reduced power due to the small post-selection probability. Here we experimentally demonstrate recycled weak-value measurements using a pulsed light source and optical switch to enable nearly deterministic weak-value amplification of a mirror tilt. Using photon counting detectors, we demonstrate a signal improvement by a factor of $4.4 \pm 0.2$ and a signal-to-noise ratio improvement of $2.10 \pm 0.06$, compared to a single-pass weak-value experiment, and also compared to a conventional direct measurement of the tilt. The signal-to-noise ratio improvement could reach around 6 for the parameters of this experiment, assuming lower loss elements.


title: Selective mode conversion and rainbow trapping via graded elastic waveguides

abstract: We experimentally achieve wave mode conversion and rainbow trapping in an elastic waveguide loaded with an array of resonators. Rainbow trapping is a phenomenon that induces wave confinement as a result of a spatial variation of the wave velocity, here promoted by gently varying the length of consecutive resonators. By breaking the geometrical symmetry of the waveguide, we combine the wave speed reduction with a reflection mechanism that mode-converts flexural waves impinging on the array into torsional waves travelling along opposite directions. The framework presented herein may open opportunities in the context of wave manipulation through the realization of structural components with concurrent wave conversion and energy trapping capabilities.


title: Revising Johnson's table for the 21st century

abstract: What does it mean today to study a problem from a computational point of view? We focus on parameterized complexity and on Column 16 "Graph Restrictions and Their Effect" of D. S. Johnson's Ongoing guide, where several puzzles were proposed in a summary table with 30 graph classes as rows and 11 problems as columns. Several of the 330 entries remain unclassified into Polynomial or NP-complete after 35 years. We provide a full dichotomy for the Steiner Tree column by proving that the problem is NP-complete when restricted to Undirected Path graphs. We revise Johnson's summary table according to the granularity provided by the parameterized complexity for NP-complete problems.


title: Accurate outdoor ground truth based on total stations

abstract: In robotics, accurate ground-truth position fostered the development of mapping and localization algorithms through the creation of cornerstone datasets. In outdoor environments and over long distances, total stations are the most accurate and precise measurement instruments for this purpose. Most total station-based systems in the literature are limited to three Degrees Of Freedoms (DOFs), due to the use of a single-prism tracking approach. In this paper, we present preliminary work on measuring a full pose of a vehicle, bringing the referencing system to six DOFs. Three total stations are used to track in real time three prisms attached to a target platform. We describe the structure of the referencing system and the protocol for acquiring the ground truth with this system. We evaluated its precision in a variety of different outdoor environments, ranging from open-sky to forest trails, and compare this system with another popular source of reference position, the Real Time Kinematics (RTK) positioning solution. Results show that our approach is the most precise, reaching an average positional error of 10 mm and 0.6 deg. This difference in performance was particularly stark in environments where Global Navigation Satellite System (GNSS) signals can be weaker due to overreaching vegetation.


title: On the use of PN Ranging with High-rate Spectrally-efficient Modulations

abstract: In this paper, we study the feasibility of coupling the PN ranging with filtered high-order modulations, and investigate the simultaneous demodulation of a high-rate telemetry stream while tracking the PN ranging sequence. Accordingly, we design a receiver scheme that is able to perform a parallel cancellation, in closed-loop, of the ranging and the telemetry signal reciprocally. From our analysis, we find that the non-constant envelope property of the modulation causes an additional jitter on the PN ranging timing estimation that, on the other hand, can be limited by properly sizing the receiver loop bandwidth. Our study proves that the use of filtered high-order modulations combined with PN ranging outperforms the state-of-the-art in terms of spectral efficiency and achievable data rate, while having comparable ranging performance.


title: Implementing Reinforcement Learning Algorithms in Retail Supply Chains with OpenAI Gym Toolkit

abstract: From cutting costs to improving customer experience, forecasting is the crux of retail supply chain management (SCM) and the key to better supply chain performance. Several retailers are using AI/ML models to gather datasets and provide forecast guidance in applications such as Cognitive Demand Forecasting, Product End-of-Life, Forecasting, and Demand Integrated Product Flow. Early work in these areas looked at classical algorithms to improve on a gamut of challenges such as network flow and graphs. But the recent disruptions have made it critical for supply chains to have the resiliency to handle unexpected events. The biggest challenge lies in matching supply with demand. Reinforcement Learning (RL) with its ability to train systems to respond to unforeseen environments, is being increasingly adopted in SCM to improve forecast accuracy, solve supply chain optimization challenges, and train systems to respond to unforeseen circumstances. Companies like UPS and Amazon have developed RL algorithms to define winning AI strategies and keep up with rising consumer delivery expectations. While there are many ways to build RL algorithms for supply chain use cases, the OpenAI Gym toolkit is becoming the preferred choice because of the robust framework for event-driven simulations. This white paper explores the application of RL in supply chain forecasting and describes how to build suitable RL models and algorithms by using the OpenAI Gym toolkit.


title: A Study on Hand Proof for The Four-Color Theorem

abstract: For the four-color theorem that has been developed over one and half centuries, all people believe it right but without complete proof convincing all1-3. Former proofs are to find the basic four-colorable patterns on a planar graph to reduce a map coloring4-6, but the unavoidable set is almost limitless and required recoloring hardly implements by hand7-14. Another idea belongs to formal proof limited to logical operation15. However, recoloring or formal proof way may block people from discovering the inherent essence of a coloring graph. Defining creation and annihilation operations, we show that four colors are sufficient to color a map and how to color it. We find what trapped vertices and boundary-vertices are, and how they decide how many colors to be required in coloring arbitrary maps. We reveal that there is the fourth color for new adding vertex differing from any three coloring vertices in creation operation. To implement a coloring map, we also demonstrate how to color an arbitrary map by iteratively using creation and annihilation operations. We hope our hand proof is beneficial to understand the mechanisms of the four-color theorem.


title: Generation of an electromagnetic field nonminimally coupled to gravity during Higgs inflation

abstract: In the framework of Higgs inflation, we consider the electromagnetic field nonminimally coupled to gravity via the parity-preserving $\propto RF^{2}$ and parity-violating $\propto RF\tilde{F}$ terms. Using the perturbation theory to the leading order in these couplings, we study the generation of the electromagnetic field during the inflation stage. We derive the effective kinetic and axial coupling functions arising in the Einstein frame in the usual metric as well as Palatini formulations of gravity. For both formulations, we determine the power spectrum, energy density, and helicality of the generated electromagnetic fields for different values of the nonminimal coupling constants. Finally, we estimate the maximal present-day magnitude of the magnetic field as $10^{-14}-10^{-15}\,$G with the correlation length of order 10 pc.


title: Sample selection from a given dataset to validate machine learning models

abstract: The selection of a validation basis from a full dataset is often required in industrial use of supervised machine learning algorithm. This validation basis will serve to realize an independent evaluation of the machine learning model. To select this basis, we propose to adopt a "design of experiments" point of view, by using statistical criteria. We show that the "support points" concept, based on Maximum Mean Discrepancy criteria, is particularly relevant. An industrial test case from the company EDF illustrates the practical interest of the methodology.


title: Flicker Noise in Two-Dimensional Electron Gas

abstract: Using the method developed in a recent paper (Euro. Phys. J. B 92.8 (2019): 1-28) we consider $1/f$ noise in two-dimensional electron gas (2DEG). The electron coherence length of the system is considered as a basic parameter for discretizing the space, inside which the dynamics of electrons is described by quantum mechanics, while for length scales much larger than it the dynamics is semi-classical. For our model, which is based on the Thomas-Fermi-Dirac approximation, there are two control parameters: temperature $T$ and the disorder strength ($\Delta$). Our Monte Carlo studies show that the system exhibits $1/f$ noise related to the electronic avalanche size, which can serve as a model for describing the experimentally observed flicker noise in 2DEG. The power spectrum of our model scales with frequency with an exponent in the interval $0.3<\alpha_{PS}<0.6$. We numerically show that the electronic avalanches are scale-invariant with power-law behaviors in and out of the metal-insulator transition line.


title: Do Feature Attribution Methods Correctly Attribute Features?

abstract: Feature attribution methods are popular in interpretable machine learning. These methods compute the attribution of each input feature to represent its importance, but there is no consensus on the definition of "attribution", leading to many competing methods with little systematic evaluation, complicated in particular by the lack of ground truth attribution. To address this, we propose a dataset modification procedure to induce such ground truth. Using this procedure, we evaluate three common methods: saliency maps, rationales, and attentions. We identify several deficiencies and add new perspectives to the growing body of evidence questioning the correctness and reliability of these methods applied on datasets in the wild. We further discuss possible avenues for remedy and recommend new attribution methods to be tested against ground truth before deployment. The code is available at https://github.com/YilunZhou/feature-attribution-evaluation


title: A 0.502$\cdot$MaxCut Approximation using Quadratic Programming

abstract: We study the MaxCut problem for graphs $G=(V,E)$. The problem is NP-hard, there are two main approximation algorithms with theoretical guarantees: (1) the Goemans \& Williamson algorithm uses semi-definite programming to provide a 0.878MaxCut approximation (which, if the Unique Games Conjecture is true, is the best that can be done in polynomial time) and (2) Trevisan proposed an algorithm using spectral graph theory from which a 0.614MaxCut approximation can be obtained. We discuss a new approach using a specific quadratic program and prove that its solution can be used to obtain at least a 0.502MaxCut approximation. The algorithm seems to perform well in practice.


title: Certain generating functions for Cigler's polynomials

abstract: In this paper, we use the homogeneous $q$-operators [J. Difference Equ. Appl. {\bf20 } (2014), 837--851.] to derive Rogers formulas, extended Rogers formulas and Srivastava-Agarwal type bilinear generating functions for Cigler's polynomials [J. Difference Equ. Appl. {\bf 24} (2018), 479--502.]. Finally, we also derive two interesting transformation formulas between ${}_2\Phi_1, \, {}_2\Phi_2$ and ${}_3\Phi_2$.


title: Dynamical prediction of two meteorological factors using the deep neural network and the long short-term memory $(2)$

abstract: This paper presents the predictive accuracy using two-variate meteorological factors, average temperature and average humidity, in neural network algorithms. We analyze result in five learning architectures such as the traditional artificial neural network, deep neural network, and extreme learning machine, long short-term memory, and long-short-term memory with peephole connections, after manipulating the computer-simulation. Our neural network modes are trained on the daily time-series dataset during seven years (from 2014 to 2020). From the trained results for 2500, 5000, and 7500 epochs, we obtain the predicted accuracies of the meteorological factors produced from outputs in ten metropolitan cities (Seoul, Daejeon, Daegu, Busan, Incheon, Gwangju, Pohang, Mokpo, Tongyeong, and Jeonju). The error statistics is found from the result of outputs, and we compare these values to each other after the manipulation of five neural networks. As using the long-short-term memory model in testing 1 (the average temperature predicted from the input layer with six input nodes), Tonyeong has the lowest root mean squared error (RMSE) value of 0.866 $(%)$ in summer from the computer-simulation in order to predict the temperature. To predict the humidity, the RMSE is shown the lowest value of 5.732 $(%)$, when using the long short-term memory model in summer in Mokpo in testing 2 (the average humidity predicted from the input layer with six input nodes). Particularly, the long short-term memory model is is found to be more accurate in forecasting daily levels than other neural network models in temperature and humidity forecastings. Our result may provide a computer-simuation basis for the necessity of exploring and develping a novel neural network evaluation method in the future.


title: The aox--HeII EW Connection in Radio-Loud Quasars

abstract: Radio-loud quasars (RLQs) are known to produce excess X-ray emission, compared to radio-quiet quasars (RQQs) of the same luminosity, commonly attributed to jet-related emission. Recently, we found that the HeII EW and $\alpha_{\rm{ox}}$ in RQQs are strongly correlated, which suggests that their extreme-ultraviolet (EUV) and X-ray emission mechanisms are tightly related. Using 48 RLQs, we show that steep-spectrum radio quasars (SSRQs) and low radio-luminosity ($L_{\rm R}$) flat-spectrum radio quasars (FSRQs) follow the $\alpha_{\rm ox}$--HeII EW relation of RQQs. This suggests that the X-ray and EUV emission mechanisms in these types of RLQs is the same as in RQQs, and is not jet related. High-$L_{\rm R}$ FSRQs show excess X-ray emission given their HeII EW by a factor of $\approx$ 3.5, which suggests that only in this type of RLQ is the X-ray production likely jet related.


title: Guessing the buffer bound for k-synchronizability

abstract: A communicating system is $k$-synchronizable if all of the message sequence charts representing the executions can be divided into slices of $k$ sends followed by $k$ receptions. It was previously shown that, for a fixed given $k$, one could decide whether a communicating system is $k$-synchronizable. This result is interesting because the reachability problem can be solved for $k$-synchronizable systems. However, the decision procedure assumes that the bound $k$ is fixed. In this paper we improve this result and show that it is possible to decide if such a bound $k$ exists.


title: Bose-Fermi transmutation for one-dimensional harmonic trap

abstract: Using Density Matrix renormalization group (DMRG), we study the ground state properties of spin one-half fermions and scalar bosons in the soft-core limit, with weak s-wave inter and intra species interactions. We considered the system subject to one-dimensional (1D) optical lattice and a superimposed potential at zero temperature, in the framework of Bose-Fermi-Hubbard model. We found that for certain fillings and interaction parameters, a transmutation occurs between the ground states of bosons and fermions when the densities are exchanged. We too report that the density distributions of bosons and fermions overlap with each other in a bosonic and fermionic Mott plateau, when the interaction parameters fulfill the relationship $U_{BB}>U_{BF}<U_{BF}$. We also find that the fermions are repelled out of the central region of the trap for sufficiently strong $U_{FF}$ interaction, exhibiting phase separation of Bose and Fermi components.


title: Modeling liquid water by climbing up Jacob's ladder in density functional theory facilitated by using deep neural network potentials

abstract: Within the framework of Kohn-Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob's ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric-isothermal ensemble. In order to generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water towards the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional.


title: Transmission of electromagnetic waves through nonlinear optical layers

abstract: The excitation of soliton states in optical layers exhibiting Kerr nonlinearities is theoretically investigated. The optical transmission coefficient is obtained as a function of the nonlinearity power of an incident monochromatic electromagnetic wave polarized in the transverse electric configuration. In materials exhibiting a defocusing nonlinearity, soliton excitations are not observed if the nonlinearity power associated with the incident electromagnetic wave is below a well-defined value. At such a limiting value of the nonlinearity power, a soliton excitation with an electric-field amplitude independent of the position is found within the nonlinear layer, regardless of the incident-wave frequency value.


title: Nonparametric Test for Volatility in Clustered Multiple Time Series

abstract: Contagion arising from clustering of multiple time series like those in the stock market indicators can further complicate the nature of volatility, rendering a parametric test (relying on asymptotic distribution) to suffer from issues on size and power. We propose a test on volatility based on the bootstrap method for multiple time series, intended to account for possible presence of contagion effect. While the test is fairly robust to distributional assumptions, it depends on the nature of volatility. The test is correctly sized even in cases where the time series are almost nonstationary. The test is also powerful specially when the time series are stationary in mean and that volatility are contained only in fewer clusters. We illustrate the method in global stock prices data.


title: The Chern class for $K_3$ and the cyclic quantum dilogarithm

abstract: In this note we confirm the conjecture of Calegari, Garoufalidis and Zagier in their recent paper in Ann. Sci. \'{E}cole Normale Sup. that $R_\zeta=c_\zeta^2$, where $R_\zeta$ is their map on $K_3$ defined using the cyclic quantum dilogarithm and $c_\zeta$ is the Chern class map on $K_3$.


title: Regional poverty in Bulgaria in the period 2008-2019

abstract: Background: Poverty among the population of a country is one of the most disputable topics in social studies. Many researchers devote their work to identifying the factors that influence it most. Bulgaria is one of the EU member states with the highest poverty levels. Regional facets of social exclusion and risks of poverty among the population are a key priority of the National Development Strategy for the third decade of 21st century. In order to mitigate the regional poverty levels it is necessary for the social policy makers to pay more attention to the various factors expected to influence these levels. Results: Poverty reduction is observed in most areas of the country. The regions with obviously favorable developments are Sofia district, Pernik, Pleven, Lovech, Gabrovo, Veliko Tarnovo, Silistra, Shumen, Stara Zagora, Smolyan, Kyustendil and others. Increased levels of poverty are found for Razgrad and Montana districts. It was fond that the reduction in the risk of poverty is associated to the increase in employment, investment, and housing. Conclusion: The social policy making needs to be aware of the fact that the degree of exposition to risk of poverty and social exclusion significantly relates to the levels of regional employment, investment and housing.


title: Polytime reductions of AF-algebraic problems

abstract: We assess the computational complexity of several decision problems concerning (Murray-von Neumann) equivalence classes of projections of AF-algebras whose Elliott classifier is lattice-ordered. We construct polytime reductions among many of these problems.


title: Global Symmetries of the Kepler Problem

abstract: The global symmetry transformations generated by Runge-Lenz vector of twodimensional Kepler problem are explicitly described. They are given in terms of SU(2) left group multiplication with group elements being suitably parametrized by phase space points. The resulting non-linear action of SU(2) on the phase space is characterized in terms of the theory of nonlinear realizations developed in Phys. Rev. 177 (1969) 2239.


title: Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

abstract: We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.


title: Enhancing Safety of Students with Mobile Air Filtration during School Reopening from COVID-19

abstract: The paper discusses how robots enable occupant-safe continuous protection for students when schools reopen. Conventionally, fixed air filters are not used as a key pandemic prevention method for public indoor spaces because they are unable to trap the airborne pathogens in time in the entire room. However, by combining the mobility of a robot with air filtration, the efficacy of cleaning up the air around multiple people is largely increased. A disinfection co-robot prototype is thus developed to provide continuous and occupant-friendly protection to people gathering indoors, specifically for students in a classroom scenario. In a static classroom with students sitting in a grid pattern, the mobile robot is able to serve up to 14 students per cycle while reducing the worst-case pathogen dosage by 20%, and with higher robustness compared to a static filter. The extent of robot protection is optimized by tuning the passing distance and speed, such that a robot is able to serve more people given a threshold of worst-case dosage a person can receive.


title: Coulomb explosion imaging of carbon monoxide dimers

abstract: We report on experimental results obtained from collisions of slow highly charged Ar9+ ions with a carbon monoxide dimer (CO)2 target. A COLTRIMS setup and a Coulomb explosion imaging approach are used to reconstruct the structure of the CO dimers. The three dimensional structure is deduced from the 2-body and 3-body dissociation channels from which both the intermolecular bond length and the relative orientation of the two molecules are determined. For the 3-body channels, the experimental data are interpreted with the help of a classical model in which the trajectories of the three emitted fragments are numerically integrated. We measured the equilibrium intermolecular distance to be Re = 4.2 A. The orientation of both CO molecules with respect to the dimer axis is found to be quasi-isotropic due to the large vibrational temperature of the gas jet.


title: Genotype-Guided Radiomics Signatures for Recurrence Prediction of Non-Small-Cell Lung Cancer

abstract: Non-small cell lung cancer (NSCLC) is a serious disease and has a high recurrence rate after the surgery. Recently, many machine learning methods have been proposed for recurrence prediction. The methods using gene data have high prediction accuracy but require high cost. Although the radiomics signatures using only CT image are not expensive, its accuracy is relatively low. In this paper, we propose a genotype-guided radiomics method (GGR) for obtaining high prediction accuracy with low cost. We used a public radiogenomics dataset of NSCLC, which includes CT images and gene data. The proposed method is a two-step method, which consists of two models. The first model is a gene estimation model, which is used to estimate the gene expression from radiomics features and deep features extracted from computer tomography (CT) image. The second model is used to predict the recurrence using the estimated gene expression data. The proposed GGR method designed based on hybrid features which is combination of handcrafted-based and deep learning-based. The experiments demonstrated that the prediction accuracy can be improved significantly from 78.61% (existing radiomics method) and 79.14% (deep learning method) to 83.28% by the proposed GGR.


title: What Are Bayesian Neural Network Posteriors Really Like?

abstract: The posterior over Bayesian neural network (BNN) parameters is extremely high-dimensional and non-convex. For computational reasons, researchers approximate this posterior using inexpensive mini-batch methods such as mean-field variational inference or stochastic-gradient Markov chain Monte Carlo (SGMCMC). To investigate foundational questions in Bayesian deep learning, we instead use full-batch Hamiltonian Monte Carlo (HMC) on modern architectures. We show that (1) BNNs can achieve significant performance gains over standard training and deep ensembles; (2) a single long HMC chain can provide a comparable representation of the posterior to multiple shorter chains; (3) in contrast to recent studies, we find posterior tempering is not needed for near-optimal performance, with little evidence for a "cold posterior" effect, which we show is largely an artifact of data augmentation; (4) BMA performance is robust to the choice of prior scale, and relatively similar for diagonal Gaussian, mixture of Gaussian, and logistic priors; (5) Bayesian neural networks show surprisingly poor generalization under domain shift; (6) while cheaper alternatives such as deep ensembles and SGMCMC methods can provide good generalization, they provide distinct predictive distributions from HMC. Notably, deep ensemble predictive distributions are similarly close to HMC as standard SGLD, and closer than standard variational inference.


title: Integrating 6LoWPAN Security with RPL Using The Chained Secure Mode Framework

abstract: The IPv6 over Low-powered Wireless Personal Area Network (6LoWPAN) protocol was introduced to allow the transmission of Internet Protocol version 6 (IPv6) packets using the smaller-size frames of the IEEE 802.15.4 standard, which is used in many Internet of Things (IoT) networks. The primary duty of the 6LoWPAN protocol is packet fragmentation and reassembly. However, the protocol standard currently does not include any security measures, not even authenticating the fragments immediate sender. This lack of immediate-sender authentication opens the door for adversaries to launch several attacks on the fragmentation process, such as the buffer-reservation attacks that lead to a Denial of Service (DoS) attack and resource exhaustion of the victim nodes. This paper proposes a security integration between 6LoWPAN and the Routing Protocol for Low Power and Lossy Networks (RPL) through the Chained Secure Mode (CSM) framework as a possible solution. Since the CSM framework provides a mean of immediate-sender trust, through the use of Network Coding (NC), and an integration interface for the other protocols (or mechanisms) to use this trust to build security decisions, 6LoWPAN can use this integration to build a chain-of-trust along the fragments routing path. A proof-of-concept implementation was done in Contiki Operating System (OS), and its security and performance were evaluated against an external adversary launching a buffer-reservation attack. The results from the evaluation showed significant mitigation of the attack with almost no increase in power consumption, which presents the great potential for such integration to secure the forwarding process at the 6LoWPAN Adaptation Layer


title: The Raise Regression: Justification, properties and application

abstract: Multicollinearity produces an inflation in the variance of the Ordinary Least Squares estimators due to the correlation between two or more independent variables (including the constant term). A widely applied solution is to estimate with penalized estimators (such as the ridge estimator, the Liu estimator, etc.) which exchange the mean square error by the bias. Although the variance diminishes with these procedures, all seems to indicate that the inference is lost and also the goodness of fit. Alternatively, the raise regression (\cite{Garcia2011} and \cite{Salmeron2017}) allows the mitigation of the problems generated by multicollinearity but without losing the inference and keeping the coefficient of determination. This paper completely formalizes the raise estimator summarizing all the previous contributions: its mean square error, the variance inflation factor, the condition number, the adequate selection of the variable to be raised, the successive raising and the relation between the raise and the ridge estimator. As a novelty, it is also presented the estimation method, the relation between the raise and the residualization, it is analyzed the norm of the estimator and the behaviour of the individual and joint significance test and the behaviour of the mean square error and the coefficient of variation. The usefulness of the raise regression as alternative to mitigate the multicollinearity is illustrated with two empirical applications.


title: Parallel Projection---An Improved Return Mapping Algorithm for Finite Element Modeling of Shape Memory Alloys

abstract: We present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.


title: Surpassing the Energy Resolution Limit with ferromagnetic torque sensors

abstract: We discuss the fundamental noise limitations of a ferromagnetic torque sensor based on a levitated magnet in the tipping regime. We evaluate the optimal magnetic field resolution taking into account the thermomechanical noise and the mechanical detection noise at the standard quantum limit (SQL). We find that the Energy Resolution Limit (ERL), pointed out in recent literature as a relevant benchmark for most classes of magnetometers, can be surpassed by many orders of magnitude. Moreover, similarly to the case of a ferromagnetic gyroscope, it is also possible to surpass the standard quantum limit for magnetometry with independent spins, arising from spin-projection noise. Our finding indicates that magnetomechanical systems optimized for magnetometry can achieve a magnetic field resolution per unit volume several orders of magnitude better than any conventional magnetometer. We discuss possible implications, focusing on fundamental physics problems such as the search for exotic interactions beyond the standard model.


title: Predicate Invention by Learning From Failures

abstract: Discovering novel high-level concepts is one of the most important steps needed for human-level AI. In inductive logic programming (ILP), discovering novel high-level concepts is known as predicate invention (PI). Although seen as crucial since the founding of ILP, PI is notoriously difficult and most ILP systems do not support it. In this paper, we introduce POPPI, an ILP system that formulates the PI problem as an answer set programming problem. Our experiments show that (i) PI can drastically improve learning performance when useful, (ii) PI is not too costly when unnecessary, and (iii) POPPI can substantially outperform existing ILP systems.


title: Detect opinion-based groups and reveal polarisation in survey data

abstract: Networks, representing attitudinal survey data, expose the structure of opinion-based groups. We make use of these network projections to identify the groups reliably through community detection algorithms and to examine social-identity-based groups. Our goal is to present a method for revealing polarisation and opinion-based in attitudinal surveys. This method can be broken down into the following steps: data preparation, construction of similarity-based networks, algorithmic identification of opinion-based groups, and identification of important items for community structure. We assess the method's performance and possible scope for applying it to empirical data and to a broad range of synthetic data sets. The empirical data application points out possible conclusions (i.e. social-identity polarisation), whereas the synthetic data sets mark out the method's boundaries. Next to an application example on political attitude survey, our results suggest that the method works for various surveys but is also moderated by the efficacy of the community detection algorithms. Concerning the identification of opinion-based groups, we provide a solid method to rank the item's influence on group formation and as a group identifier. We discuss how this network approach for identifying polarisation can classify non-overlapping opinion-based groups even in the absence of extreme opinions.


title: Reply to Shvaika et al.: Presence of a boson peak in anharmonic phonon models with Akhiezer-type damping

abstract: We reply to the Comment by Svhaika, Ruocco, Schirmacher and collaborators. There were two accidental mistakes in our original paper (Phys. Rev. Lett. 112, 145501 (2019)), which have been now corrected. All the physical conclusions and results of the original paper, including the prediction of boson peak due to anharmonicity, remain valid in the corrected version.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.soft cond-mat.stat-mech


title: Photonic co-processors in HPC: using LightOn OPUs for Randomized Numerical Linear Algebra

abstract: Randomized Numerical Linear Algebra (RandNLA) is a powerful class of methods, widely used in High Performance Computing (HPC). RandNLA provides approximate solutions to linear algebra functions applied to large signals, at reduced computational costs. However, the randomization step for dimensionality reduction may itself become the computational bottleneck on traditional hardware. Leveraging near constant-time linear random projections delivered by LightOn Optical Processing Units we show that randomization can be significantly accelerated, at negligible precision loss, in a wide range of important RandNLA algorithms, such as RandSVD or trace estimators.


title: Discriminative-Generative Dual Memory Video Anomaly Detection

abstract: Recently, people tried to use a few anomalies for video anomaly detection (VAD) instead of only normal data during the training process. A side effect of data imbalance occurs when a few abnormal data face a vast number of normal data. The latest VAD works use triplet loss or data re-sampling strategy to lessen this problem. However, there is still no elaborately designed structure for discriminative VAD with a few anomalies. In this paper, we propose a DiscRiminative-gEnerative duAl Memory (DREAM) anomaly detection model to take advantage of a few anomalies and solve data imbalance. We use two shallow discriminators to tighten the normal feature distribution boundary along with a generator for the next frame prediction. Further, we propose a dual memory module to obtain a sparse feature representation in both normality and abnormality space. As a result, DREAM not only solves the data imbalance problem but also learn a reasonable feature space. Further theoretical analysis shows that our DREAM also works for the unknown anomalies. Comparing with the previous methods on UCSD Ped1, UCSD Ped2, CUHK Avenue, and ShanghaiTech, our model outperforms all the baselines with no extra parameters. The ablation study demonstrates the effectiveness of our dual memory module and discriminative-generative network.


title: Properties of the Support of the Capacity-Achieving Distribution of the Amplitude-Constrained Poisson Noise Channel

abstract: This work considers a Poisson noise channel with an amplitude constraint. It is well-known that the capacity-achieving input distribution for this channel is discrete with finitely many points. We sharpen this result by introducing upper and lower bounds on the number of mass points. Concretely, an upper bound of order $\mathsf{A} \log^2(\mathsf{A})$ and a lower bound of order $\sqrt{\mathsf{A}}$ are established where $\mathsf{A}$ is the constraint on the input amplitude. In addition, along the way, we show several other properties of the capacity and capacity-achieving distribution. For example, it is shown that the capacity is equal to $ - \log P_{Y^\star}(0)$ where $P_{Y^\star}$ is the optimal output distribution. Moreover, an upper bound on the values of the probability masses of the capacity-achieving distribution and a lower bound on the probability of the largest mass point are established. Furthermore, on the per-symbol basis, a nonvanishing lower bound on the probability of error for detecting the capacity-achieving distribution is established under the maximum a posteriori rule.


title: Anti-Poiseuille flow in neutral graphene

abstract: Hydrodynamic flow of charge carriers in graphene is an energy flow unlike the usual mass flow in conventional fluids. In neutral graphene, the energy flow is decoupled from the electric current, making it difficult to observe the hydrodynamic effects and measure the viscosity of the electronic fluid by means of electric current measurements. In particular, we show that the hallmark Poiseuille flow in a narrow channel cannot be driven by the electric field irrespective of boundary conditions at the channel edges. Nevertheless one can observe nonuniform current densities similarly to the case of the well-known ballistic-diffusive crossover. The standard diffusive behavior with the uniform current density across the channel is achieved under the assumptions of specular scattering on the channel boundaries. This flow can also be made nonuniform by applying weak magnetic fields. In this case, the curvature of the current density profile is determined by the quasiparticle recombination processes dominated by the disorder-assisted electron-phonon scattering -- the so-called supercollisions.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.stat-mech cond-mat.str-el physics.flu-dyn


title: Machine-Learning Assisted Optimization Strategies for Phase Change Materials Embedded within Electronic Packages

abstract: Leveraging the latent heat of phase change materials (PCMs) can reduce the peak temperatures and transient variations in temperature in electronic devices. But as the power levels increase, the thermal conduction pathway from the heat source to the heat sink limits the effectiveness of these systems. In this work, we evaluate embedding the PCM within the silicon device layer of an electronic device to minimize the thermal resistance between the source and the PCM to minimize this thermal resistance and enhance the thermal performance of the device. The geometry and material properties of the embedded PCM regions are optimized using a combination of parametric and machine learning algorithms. For a fixed geometry, considering commercially available materials, Solder 174 significantly outperforms other organic and metallic PCMs. Also with a fixed geometry, the optimal melting points to minimize the peak temperature is higher than the optimal melting point to minimize the amplitude of the transient temperature oscillation, and both optima increase with increasing heater power. Extending beyond conventional optimization strategies, genetic algorithms and particle swarm optimization with and without neural network surrogate models are used to enable optimization of many geometric and material properties. For the test case evaluated, the optimized geometries and properties are similar between all ML-assisted algorithms, but the computational time depends on the technique. Ultimately, the optimized design with embedded phase change materials reduces the maximum temperature rise by 19% and the fluctuations by up to 88% compared to devices without PCM.


title: Quantum f(R) gravity and AdS/CFT

abstract: We propose to study the entanglement entropy in braneworld modified gravity. We show that the $d+1$-dimensional action $I_{d+1}$ is the dual of $d-1$% -dimensional entanglement entropy $S_{d-1}^{EE}$. Moreover, we remark that the generalization of action-entropy shows us a new form of $f(R)$ gravity which is in good agreement with the most choices of $f(R)$ gravity in the literature. We show also that there are two copies of the AdS spaces in holographic entanglement $f(R)$ gravity. We have proposed that the time is a holographic projection of the hidden $d$-brane on the visible $d$-brane. To determine the geometry of this holographic projection, we have used the $% AdS_{2}$ geometry of the black hole. Namely, the past of events is memorized over an $\mathbb{S}^{1}$. We show that the time projection is a stereographic projection from a time sphere $\mathbb{S}^{d-2}$ on the hidden $d$-brane to the visible $d$-brane. We suggest that there is a difference between the time in classical gravity and the time in quantum gravity.


title: Customizable Reference Runtime Monitoring of Neural Networks using Resolution Boxes

abstract: Classification neural networks fail to detect inputs that do not fall inside the classes they have been trained for. Runtime monitoring techniques on the neuron activation pattern can be used to detect such inputs. We present an approach for monitoring classification systems via data abstraction. Data abstraction relies on the notion of box with a resolution. Box-based abstraction consists in representing a set of values by its minimal and maximal values in each dimension. We augment boxes with a notion of resolution and define their clustering coverage, which is intuitively a quantitative metric that indicates the abstraction quality. This allows studying the effect of different clustering parameters on the constructed boxes and estimating an interval of sub-optimal parameters. Moreover, we automatically construct monitors that leverage both the correct and incorrect behaviors of a system. This allows checking the size of the monitor abstractions and analyzing the separability of the network. Monitors are obtained by combining the sub-monitors of each class of the system placed at some selected layers. Our experiments demonstrate the effectiveness of our clustering coverage estimation and show how to assess the effectiveness and precision of monitors according to the selected clustering parameter and monitored layers.


title: Deterministic Algorithms for the Hidden Subgroup Problem

abstract: We present deterministic algorithms for the Hidden Subgroup Problem. The first algorithm, for abelian groups, achieves the same asymptotic worst-case query complexity as the optimal randomized algorithm, namely O($\sqrt{ n}\,$), where $n$ is the order of the group. The analogous algorithm for non-abelian groups comes within a $\sqrt{ \log n}$ factor of the optimal randomized query complexity. The best known randomized algorithm for the Hidden Subgroup Problem has expected query complexity that is sensitive to the input, namely O($\sqrt{ n/m}\,$), where $m$ is the order of the hidden subgroup. In the first version of this article (arXiv:2104.14436v1 [cs.DS]), we asked if there is a deterministic algorithm whose query complexity has a similar dependence on the order of the hidden subgroup. Prompted by this question, Ye and Li (arXiv:2110.00827v1 [cs.DS]) present deterministic algorithms for abelian groups which solve the problem with O($\sqrt{ n/m }\,$ ) queries, and find the hidden subgroup with O($\sqrt{ n (\log m) / m} + \log m$) queries. Moreover, they exhibit instances which show that in general, the deterministic query complexity of the problem may be o($\sqrt{ n/m } \,$), and that of finding the entire subgroup may also be o($\sqrt{ n/m } \,$) or even $\omega(\sqrt{ n/m } \,)$. We present a different deterministic algorithm for the Hidden Subgroup Problem that also has query complexity O($\sqrt{ n/m }\,$) for abelian groups. The algorithm is arguably simpler. Moreover, it works for non-abelian groups, and has query complexity O($\sqrt{ (n/m) \log (n/m) }\,$) for a large class of instances, such as those over supersolvable groups. We build on this to design deterministic algorithms to find the hidden subgroup for all abelian and some non-abelian instances, at the cost of a $\log m$ multiplicative factor increase in the query complexity.


title: Overlap Times in the Infinite Server Queue

abstract: Imagine, you enter a grocery store to buy food. How many peopledo you overlap with in this store? How much time do you overlap witheach person in the store? In this paper, we answer these questions bystudying the overlap times between customers in the infinite serverqueue. We compute in closed form the steady state distribution ofthe overlap time between a pair of customers and the distribution ofthe number of customers that an arriving customer will overlap with.Finally, we define a residual process that counts the number of over-lapping customers that overlap in the queue for at least{\delta}time unitsand compute its mean, variance, and distribution in the exponentialservice setting


title: Three-electron correlations in strong laser field ionization: Spin induced effects

abstract: Strong field processes in the non-relativistic regime are insensitive to the electron spin, i.e. the observables appear to be independent of this electron property. This does not have to be the case for several active electrons where Pauli principle may affect the their dynamics. We exemplify this statement studying model atoms with three active electrons interacting with strong pulsed radiation, using an ab-initio time-dependent Schr\"odinger equation on a grid. In our restricted dimensionality model we are able, for the first time, to analyse momenta correlations of the three outgoing electrons using Dalitz plots. We show that significant differences are obtained between model Neon and Nitrogen atoms. These differences are traced back to the different symmetries of the electronic wavefunctions, and directly related to the different initial state spin components.


title: High-Dimensional Neural Network Potentials for Magnetic Systems Using Spin-Dependent Atom-Centered Symmetry Functions

abstract: Machine learning potentials have emerged as a powerful tool to extend the time and length scales of first principles-quality simulations. Still, most machine learning potentials cannot distinguish different electronic spin orientations and thus are not applicable to materials in different magnetic states. Here, we propose spin-dependent atom-centered symmetry functions as a new type of descriptor taking the atomic spin degrees of freedom into account. When used as input for a high-dimensional neural network potential (HDNNP), accurate potential energy surfaces of multicomponent systems describing multiple magnetic states can be constructed. We demonstrate the performance of these magnetic HDNNPs for the case of manganese oxide, MnO. We show that the method predicts the magnetically distorted rhombohedral structure in excellent agreement with density functional theory and experiment. Its efficiency allows to determine the N\'{e}el temperature considering structural fluctuations, entropic effects, and defects. The method is general and is expected to be useful also for other types of systems like oligonuclear transition metal complexes.


title: Fourth- and fifth-order virial expansion of harmonically trapped fermions at unitarity

abstract: By generalizing our automated algebra approach from homogeneous space to harmonically trapped systems, we have calculated the fourth- and fifth-order virial coefficients of universal spin-1/2 fermions in the unitary limit, confined in an isotropic harmonic potential. We present results for said coefficients as a function of trapping frequency (or, equivalently, temperature), which compare favorably with previous Monte Carlo calculations (available only at fourth order) as well as with our previous estimates in the untrapped limit (high temperature, low frequency). We use our estimates of the virial expansion, together with resummation techniques, to calculate the compressibility and spin susceptibility.


title: A Careful Reassessment of Globular Cluster Multiple Population Radial Distributions with Sloan Digital Sky Survey and Johnson-Cousins Broadband Photometry

abstract: Inconsistencies regarding the nature of globular cluster multiple population radial distributions is a matter for concern given their role in testing or validating cluster dynamical evolution modeling. In this study, we present a re-analysis of eight globular cluster radial distributions using publicly available ground-based ugriz and UBVRI photometry; correcting for a systematic error identified in the literature. We detail the need for including and considering not only K-S probabilities but critical K-S statistic values as well when drawing conclusions from radial distributions, as well as the impact of sample incompleteness. Revised cumulative radial distributions are presented, and the literature of each cluster reviewed to provide a fuller picture of our results. We find that many multiple populations are not as segregated as once thought, and that there is a pressing need for better understanding of the spatial distributions of multiple populations in globular clusters.


title: Toric non-equalized flips associated to $\mathbb{C}^*$-actions

abstract: Starting from $\mathbb{C}^*$-actions on complex projective varieties, we construct and investigate birational maps among the corresponding extremal fixed point components. We study the case in which such birational maps are locally described by toric flips, either of Atiyah type or so called non-equalized. We relate this notion of toric flip with the property of the action being non-equalized. Moreover, we find explicit examples of rational homogeneous varieties admitting a $\mathbb{C}^*$-action whose weighted blow-up at the extremal fixed point components gives a birational map among two projective varieties that is locally a toric non-equalized flip.


title: A complete Effective Field Theory for Dark Matter

abstract: We present an effective field theory describing the relevant interactions of the Standard Model with an electrically neutral particle that can account for the dark matter in the Universe. The possible mediators of these interactions are assumed to be heavy. The dark matter candidates that we consider have spin 0, 1/2 or 1, belong to an electroweak multiplet with arbitrary isospin and hypercharge and their stability at cosmological scales is guaranteed by imposing a $\mathbb{Z}_2$ symmetry. We present the most general framework for describing the interaction of the dark matter with standard particles, and construct a general non-redundant basis of the gauge-invariant operators up to dimension six. The basis includes multiplets with non-vanishing hypercharge, which can also be viable DM candidates. We give two examples illustrating the phenomenological use of such a general effective framework. First, we consider the case of a scalar singlet, provide convenient semi-analytical expressions for the relevant dark matter observables, use present experimental data to set constraints on the Wilson coefficients of the operators, and show how the interplay of different operators can open new allowed windows in the parameter space of the model. Then we study the case of a lepton isodoublet, which involves co-annihilation processes, and we discuss the impact of the operators on the particle mass splitting and direct detection cross sections. These examples highlight the importance of the contribution of the various non-renormalizable operators, which can even dominate over the gauge interactions in certain cases.


title: Improving Simulations with Symmetry Control Neural Networks

abstract: The dynamics of physical systems is often constrained to lower dimensional sub-spaces due to the presence of conserved quantities. Here we propose a method to learn and exploit such symmetry constraints building upon Hamiltonian Neural Networks. By enforcing cyclic coordinates with appropriate loss functions, we find that we can achieve improved accuracy on simple classical dynamics tasks. By fitting analytic formulae to the latent variables in our network we recover that our networks are utilizing conserved quantities such as (angular) momentum.


title: Trakhtenbrot's Theorem in Coq: Finite Model Theory through the Constructive Lens

abstract: We study finite first-order satisfiability (FSAT) in the constructive setting of dependent type theory. Employing synthetic accounts of enumerability and decidability, we give a full classification of FSAT depending on the first-order signature of non-logical symbols. On the one hand, our development focuses on Trakhtenbrot's theorem, stating that FSAT is undecidable as soon as the signature contains an at least binary relation symbol. Our proof proceeds by a many-one reduction chain starting from the Post correspondence problem. On the other hand, we establish the decidability of FSAT for monadic first-order logic, i.e. where the signature only contains at most unary function and relation symbols, as well as the enumerability of FSAT for arbitrary enumerable signatures. To showcase an application of Trakhtenbrot's theorem, we continue our reduction chain with a many-one reduction from FSAT to separation logic. All our results are mechanised in the framework of a growing Coq library of synthetic undecidability proofs.


title: The effects of spanwise confinement on stratified shear instabilities

abstract: We consider the influence of transverse confinement on the instability properties of velocity and density distributions reminiscent of those pertaining to exchange flows in stratified inclined ducts, such as the recent experiment of Lefauve et al. (J. Fluid Mech. 848, 508-544, 2018). Using a normal mode streamwise and temporal expansion for flows in ducts with various aspect ratios $B$ and non-trivial transverse velocity profiles, we calculate two-dimensional (2D) dispersion relations with associated eigenfunctions varying in the 'crosswise' direction, in which the density varies, and the spanwise direction, both normal to the duct walls and to the flow direction. We also compare these 2D dispersion relations to the so-called one-dimensional (1D) dispersion relation obtained for spanwise invariant perturbations, for different aspect ratios $B$ and bulk Richardson numbers $Ri_b$. In this limited parameter space, the presence of lateral walls has a stabilizing effect. Furthermore, accounting for spanwise-varying perturbations results in a plethora of unstable modes, the number of which increases as the aspect ratio is increased. These modes present an odd-even regularity in their spatial structures, which is rationalized by comparison to the so-called one-dimensional oblique (1D-O) dispersion relation obtained for oblique waves. Finally, we show that in most cases, the most unstable 2D mode is the one that oscillates the least in the spanwise direction, as a consequence of viscous damping. However, in a limited region of the parameter space and in the absence of stratification, we show that a secondary mode with a more complex `twisted' structure dominated by crosswise vorticity becomes more unstable than the least oscillating Kelvin-Helmholtz mode associated with spanwise vorticity.


title: Parallel implementation of a compatible high-order meshless method for the Stokes' equations

abstract: A parallel implementation of a compatible discretization scheme for steady-state Stokes problems is presented in this work. The scheme uses generalized moving least squares to generate differential operators and apply boundary conditions. This meshless scheme allows a high-order convergence for both the velocity and pressure, while also incorporates finite-difference-like sparse discretization. Additionally, the method is inherently scalable: the stencil generation process requires local inversion of matrices amenable to GPU acceleration, and the divergence-free treatment of velocity replaces the traditional saddle point structure of the global system with elliptic diagonal blocks amenable to algebraic multigrid. The implementation in this work uses a variety of Trilinos packages to exploit this local and global parallelism, and benchmarks demonstrating high-order convergence and weak scalability are provided.


title: On smoothing of plurisubharmonic functions on unbounded domains

abstract: We prove that for every $n \ge 2$, there exists a pseudoconvex domain $\Omega \subset \mathbb{C}^n$ such that $\mathfrak{c}^0(\Omega) \subsetneq \mathfrak{c}^1(\Omega)$, where $\mathfrak{c}^k(\Omega)$ denotes the core of $\Omega$ with respect to $\mathcal{C}^k$-smooth plurisubharmonic functions on $\Omega$. Moreover, we show that there exists a bounded continuous plurisubharmonic function on $\Omega$ that is not the pointwise limit of a sequence of $\mathcal{C}^1$-smooth bounded plurisubharmonic functions on $\Omega$.


title: MUSE: Multi-faceted Attention for Signed Network Embedding

abstract: Signed network embedding is an approach to learn low-dimensional representations of nodes in signed networks with both positive and negative links, which facilitates downstream tasks such as link prediction with general data mining frameworks. Due to the distinct properties and significant added value of negative links, existing signed network embedding methods usually design dedicated methods based on social theories such as balance theory and status theory. However, existing signed network embedding methods ignore the characteristics of multiple facets of each node and mix them up in one single representation, which limits the ability to capture the fine-grained attentions between node pairs. In this paper, we propose MUSE, a MUlti-faceted attention-based Signed network Embedding framework to tackle this problem. Specifically, a joint intra- and inter-facet attention mechanism is introduced to aggregate fine-grained information from neighbor nodes. Moreover, balance theory is also utilized to guide information aggregation from multi-order balanced and unbalanced neighbors. Experimental results on four real-world signed network datasets demonstrate the effectiveness of our proposed framework.


title: Discontinuous Galerkin and $C^0$-IP finite element approximation of periodic Hamilton--Jacobi--Bellman--Isaacs problems with application to numerical homogenization

abstract: In the first part of the paper, we study the discontinuous Galerkin (DG) and $C^0$ interior penalty ($C^0$-IP) finite element approximation of the periodic strong solution to the fully nonlinear second-order Hamilton--Jacobi--Bellman--Isaacs (HJBI) equation with coefficients satisfying the Cordes condition. We prove well-posedness and perform abstract a posteriori and a priori analyses which apply to a wide family of numerical schemes. These periodic problems arise as the corrector problems in the homogenization of HJBI equations. The second part of the paper focuses on the numerical approximation to the effective Hamiltonian of ergodic HJBI operators via DG/$C^0$-IP finite element approximations to approximate corrector problems. Finally, we provide numerical experiments demonstrating the performance of the numerical schemes.


title: Like-charge attraction at short distances in a charge-asymmetric two-dimensional two-component plasma: Exact results

abstract: We determine exactly the short-distance effective potential between two "guest" charges immersed in a two-dimensional two-component charge-asymmetric plasma composed of positively ($q_1 = +1$) and negatively ($q_2 = -1/2$) charged point particles. The result is valid over the whole regime of stability, where the Coulombic coupling (dimensionless inverse temperature) $\beta <4$. At high Coulombic coupling $\beta>2$, this model features like-charge attraction. Also, there cannot be repulsion between opposite-charges at short-distances, at variance with large-distance interactions.


title: TGV-based restoration of Poissonian images with automatic estimation of the regularization parameter

abstract: The problem of restoring images corrupted by Poisson noise is common in many application fields and, because of its intrinsic ill posedness, it requires regularization techniques for its solution. The effectiveness of such techniques depends on the value of the regularization parameter balancing data fidelity and regularity of the solution. Here we consider the Total Generalized Variation regularization introduced in [SIAM J. Imag. Sci, 3(3), 492-526, 2010], which has demonstrated its ability of preserving sharp features as well as smooth transition variations, and introduce an automatic strategy for defining the value of the regularization parameter. We solve the corresponding optimization problem by using a 3-block version of ADMM. Preliminary numerical experiments support the proposed approach.


title: Late-time acceleration with a scalar field source: Observational constraints and statefinder diagnostics

abstract: This article discusses a dark energy cosmological model in the standard theory of gravity - general relativity with a broad scalar field as a source. Exact solutions of Einstein's field equations are derived by considering a particular form of deceleration parameter $q$, which shows a smooth transition from decelerated to accelerated phase in the evolution of the universe. The external datasets such as Hubble ($H(z)$) datasets, Supernovae (SN) datasets, and Baryonic Acoustic Oscillation (BAO) datasets are used for constraining the model par parameters appearing in the functional form of $q$. The transition redshift is obtained at $% z_{t}=0.67_{-0.36}^{+0.26}$ for the combined data set ($H(z)+SN+BAO$), where the model shows signature-flipping and is consistent with recent observations. Moreover, the present value of the deceleration parameter comes out to be $q_{0}=-0.50_{-0.11}^{+0.12}$ and the jerk parameter $% j_{0}=-0.98_{-0.02}^{+0.06}$ (close to 1) for the combined datasets, which is compatible as per Planck2018 results. The analysis also constrains the omega value i.e., $\Omega _{m_{0}}\leq 0.269$ for the smooth evolution of the scalar field EoS parameter. It is seen that energy density is higher for the effective energy density of the matter field than energy density in the presence of a scalar field. The evolution of the physical and geometrical parameters is discussed in some details with the model parameters' numerical constrained values. Moreover, we have performed the state-finder analysis to investigate the nature of dark energy.


title: Mechanisms of X-ray and gamma-ray emissions in lightning and spark discharges

abstract: X-ray and gamma-ray emissions observed in lightning and long sparks are usually connected with the bremsstrahlung of high-energy runaway electrons. Here, an alternative physical mechanism for producing X-ray and gamma-ray emissions caused by the polarization current and associated electromagnetic field moving with relativistic velocity along a curved discharge channel has been proposed. It is pointed out that lightning and spark discharges should also produce a coherent radio-frequency radiation. The influence of the conductivity and the radius of the lightning channel on the propagation velocity of electromagnetic waves, taking into account the absorption, have been investigated. The existence of fast electromagnetic surface waves propagating along the lightning discharge channel at a speed close to the speed of light in vacuum is shown. The possibility of the production of microwave, X-ray and gamma-ray emissions by a polarization current pulse moving along a curved path via synchrotron radiation mechanism during the lightning leader steps formation and the very beginning of the return stroke stage is pointed out. The existence of long tails in the power spectrum is shown, which explains observations of photon energies in the range of 10-100 MeV in the TGF, as well as measured power spectrum of laboratory spark discharge.


title: First principles studies of the surface and opto-electronic properties of ultra-thin t-Se

abstract: Selenium is a crucial earth-abundant and non-toxic semiconductor with a wide range of applications across the semiconductor industries. Selenium has drawn attention from scientific communities for its wide range of applicability: from photovoltaics to imaging devices. Its usage as a photosensitive material largely involves the synthesis of the amorphous phase (a-Se) via various experimental techniques. However, the ground state crystalline phase of this material, known as the trigonal selenium (\textit{t}-Se), is not extensively studied for its optimum electronic and optical properties. In this work, we present density functional theory (DFT) based systematic studies on the ultra-thin $(10\overline{1}0)$ surface slabs of \textit{t}-Se. We report the surface energies, work function, electronic and optical properties as a function of number of layers for $(10\overline{1}0)$ surface slabs to access its suitability for applications as a photosensitive material.


title: A game theoretic perspective on Bayesian multi-objective optimization

abstract: This chapter addresses the question of how to efficiently solve many-objective optimization problems in a computationally demanding black-box simulation context. We shall motivate the question by applications in machine learning and engineering, and discuss specific harsh challenges in using classical Pareto approaches when the number of objectives is four or more. Then, we review solutions combining approaches from Bayesian optimization, e.g., with Gaussian processes, and concepts from game theory like Nash equilibria, Kalai-Smorodinsky solutions and detail extensions like Nash-Kalai-Smorodinsky solutions. We finally introduce the corresponding algorithms and provide some illustrating results.


title: Planet-induced radio emission from the coronae of M dwarfs

abstract: There have recently been detections of radio emission from low-mass stars, some of which are indicative of star-planet interactions. Motivated by these exciting new results, here we present stellar wind models for the active planet-hosting M dwarf AU Mic. Our models incorporate the large-scale photospheric magnetic field map of the star, reconstructed using the Zeeman-Doppler Imaging method. We use our models to assess if planet-induced radio emission could be generated in the corona of AU Mic, through a mechanism analogous to the sub-Alfv\'enic Jupiter-Io interaction. In the case that AU Mic has a mass-loss rate of 27 times that of the Sun, we find that both planets b and c in the system can induce radio emission from 10 MHz to 3 GHz in the corona of the host star for the majority of their orbits, with peak flux densities of 10 mJy. Our predicted emission bears a striking similarity to that recently reported from GJ 1151 by Vedantham et al. (2020), which is indicative of being induced by a planet. Detection of such radio emission would allow us to place an upper limit on the mass-loss rate of the star.


title: Nonparametric Difference-in-Differences in Repeated Cross-Sections with Continuous Treatments

abstract: This paper studies the identification of causal effects of a continuous treatment using a new difference-in-difference strategy. Our approach allows for endogeneity of the treatment, and employs repeated cross-sections. It requires an exogenous change over time which affects the treatment in a heterogeneous way, stationarity of the distribution of unobservables and a rank invariance condition on the time trend. On the other hand, we do not impose any functional form restrictions or an additive time trend, and we are invariant to the scaling of the dependent variable. Under our conditions, the time trend can be identified using a control group, as in the binary difference-in-differences literature. In our scenario, however, this control group is defined by the data. We then identify average and quantile treatment effect parameters. We develop corresponding nonparametric estimators and study their asymptotic properties. Finally, we apply our results to the effect of disposable income on consumption.


title: Majorana bound states in semiconducting nanostructures

abstract: In this Tutorial, we give a pedagogical introduction to Majorana bound states (MBSs) arising in semiconducting nanostructures. We start by briefly reviewing the well-known Kitaev chain toy model in order to introduce some of the basic properties of MBSs before proceeding to describe more experimentally relevant platforms. Here, our focus lies on simple `minimal' models where the Majorana wave functions can be obtained explicitly by standard methods. In a first part, we review the paradigmatic model of a Rashba nanowire with strong spin-orbit interaction (SOI) placed in a magnetic field and proximitized by a conventional $s$-wave superconductor. We identify the topological phase transition separating the trivial phase from the topological phase and demonstrate how the explicit Majorana wave functions can be obtained in the limit of strong SOI. In a second part, we discuss MBSs engineered from proximitized edge states of two-dimensional (2D) topological insulators. We introduce the Jackiw-Rebbi mechanism leading to the emergence of bound states at mass domain walls and show how this mechanism can be exploited to construct MBSs. Due to their recent interest, we also include a discussion of Majorana corner states in 2D second-order topological superconductors. This Tutorial is mainly aimed at graduate students -- both theorists and experimentalists -- seeking to familiarize themselves with some of the basic concepts in the field.


title: Self-correcting longitudinal phase space in a multistage plasma accelerator

abstract: Plasma accelerators driven by intense laser or particle beams provide gigavolt-per-meter accelerating fields, promising to drastically shrink particle accelerators for high-energy physics and photon science. Applications such as linear colliders and free-electron lasers (FELs) require high energy and energy efficiency, but also high stability and beam quality. The latter includes low energy spread, which can be achieved by precise beam loading of the plasma wakefield using longitudinally shaped bunches, resulting in efficient and uniform acceleration. However, the plasma wavelength, which sets the scale for the region of very large accelerating fields to be 100 {\mu}m or smaller, requires bunches to be synchronized and shaped with extreme temporal precision, typically on the femtosecond scale. Here, a self-correction mechanism is introduced, greatly reducing the susceptibility to jitter. Using multiple accelerating stages, each with a small bunch compression between them, almost any initial bunch, regardless of current profile or injection phase, will self-correct into the current profile that flattens the wakefield, damping the relative energy spread and any energy offsets. As a consequence, staging can be used not only to reach high energies, but also to produce the exquisite beam quality and stability required for a variety of applications.


title: Twin Systems for DeepCBR: A Menagerie of Deep Learning and Case-Based Reasoning Pairings for Explanation and Data Augmentation

abstract: Recently, it has been proposed that fruitful synergies may exist between Deep Learning (DL) and Case Based Reasoning (CBR); that there are insights to be gained by applying CBR ideas to problems in DL (what could be called DeepCBR). In this paper, we report on a program of research that applies CBR solutions to the problem of Explainable AI (XAI) in the DL. We describe a series of twin-systems pairings of opaque DL models with transparent CBR models that allow the latter to explain the former using factual, counterfactual and semi-factual explanation strategies. This twinning shows that functional abstractions of DL (e.g., feature weights, feature importance and decision boundaries) can be used to drive these explanatory solutions. We also raise the prospect that this research also applies to the problem of Data Augmentation in DL, underscoring the fecundity of these DeepCBR ideas.


title: System inference via field inversion for the spatio-temporal progression of infectious diseases: Studies of COVID-19 in Michigan and Mexico

abstract: We present an approach to studying and predicting the spatio-temporal progression of infectious diseases. We treat the problem by adopting a partial differential equation (PDE) version of the Susceptible, Infected, Recovered, Deceased (SIRD) compartmental model of epidemiology, which is achieved by replacing compartmental populations by their densities. Building on our recent work (Computational Mechanics, 66, 1177, 2020), we replace our earlier use of global polynomial basis functions with those having local support, as epitomized in the finite element method, for the spatial representation of the SIRD parameters. The time dependence is treated by inferring constant parameters over time intervals that coincide with the time step in semi-discrete numerical implementations. In combination, this amounts to a scheme of field inversion of the SIRD parameters over each time step. Applied to data over ten months of 2020 for the pandemic in the US state of Michigan and to all of Mexico, our system inference via field inversion infers spatio-temporally varying PDE SIRD parameters that replicate the progression of the pandemic with high accuracy. It also produces accurate predictions, when compared against data, for a three week period into 2021. Of note is the insight that is suggested on the spatio-temporal variation of infection, recovery and death rates, as well as patterns of the population's mobility revealed by diffusivities of the compartments.


title: Analytic Spread of Filtrations and Symbolic Algebras

abstract: In this paper we define and explore the analytic spread $\ell(\mathcal I)$ of a filtration in a local ring. We show that, especially for divisorial and symbolic filtrations, some basic properties of the analytic spread of an ideal extend to filtrations, even when the filtration is non Noetherian. We also illustrate some significant differences between the analytic spread of a filtration and the analytic spread of an ideal with examples. In the case of an ideal $I$, we have the classical bounds $\mbox{ht}(I)\le\ell(I)\le \dim R$. The upper bound $\ell(\mathcal I)\le \dim R$ is true for filtrations $\mathcal I$, but the lower bound is not true for all filtrations. We show that for the filtration $\mathcal I$ of symbolic powers of a height two prime ideal $\mathfrak p$ in a regular local ring of dimension three (a space curve singularity), so that $\mbox{ht}(\mathcal I) =2$ and $\dim R=3$, we have that $0\le \ell(\mathcal I)\le 2$ and all values of 0,1 and 2 can occur. In the cases of analytic spread 0 and 1 the symbolic algebra is necessarily non-Noetherian. The symbolic algebra is non-Noetherian if and only if $\ell(\mathfrak p^{(n)})=3$ for all symbolic powers of $\mathfrak p$ and if and only if $\ell(\mathcal I_a)=3$ for all truncations $\mathcal I_a$ of $\mathcal I$.


title: Thickness Dependence of Magneto-transport Properties in Tungsten Ditelluride

abstract: We investigate the electronic structure of tungsten ditelluride (WTe$_2$) flakes with different thicknesses in magneto-transport studies. The temperature-dependent resistance and magnetoresistance (MR) measurements both confirm the breaking of carrier balance induced by thickness reduction, which suppresses the `turn-on' behavior and large positive MR. The Shubnikov-de-Haas oscillation studies further confirm the thickness-dependent change of electronic structure of WTe$_2$ and reveal a possible temperature-sensitive electronic structure change. Finally, we report the thickness-dependent anisotropy of Fermi surface, which reveals that multi-layer WTe$_2$ is an electronic 3D material and the anisotropy decreases as thickness decreases.


title: A Measure of Classicality

abstract: A striking feature of our fundamentally indeterministic quantum universe is its quasiclassical realm -- the wide range of time place and scale in which the deterministic laws of classical physics hold. Our quasiclassical realmis an emergent feature of the fundamental theories of our universe's quantum state and dynamics. There are many types of quasiclassical realms our Universe could exhibit characterized by different variables, different levels of coarse-graining, different locations in spacetime, different classical physics, and different levels of classicality.We propose a measure of classicality for quasiclassical realms, We speculate on the observable consequences of different levels of classicality especially for information gathering and utilizing systems (IGUSes) such ourselves as observers of the Universe.


title: 3D Human Action Representation Learning via Cross-View Consistency Pursuit

abstract: In this work, we propose a Cross-view Contrastive Learning framework for unsupervised 3D skeleton-based action Representation (CrosSCLR), by leveraging multi-view complementary supervision signal. CrosSCLR consists of both single-view contrastive learning (SkeletonCLR) and cross-view consistent knowledge mining (CVC-KM) modules, integrated in a collaborative learning manner. It is noted that CVC-KM works in such a way that high-confidence positive/negative samples and their distributions are exchanged among views according to their embedding similarity, ensuring cross-view consistency in terms of contrastive context, i.e., similar distributions. Extensive experiments show that CrosSCLR achieves remarkable action recognition results on NTU-60 and NTU-120 datasets under unsupervised settings, with observed higher-quality action representations. Our code is available at https://github.com/LinguoLi/CrosSCLR.


title: Towards a practical lip-to-speech conversion system using deep neural networks and mobile application frontend

abstract: Articulatory-to-acoustic (forward) mapping is a technique to predict speech using various articulatory acquisition techniques as input (e.g. ultrasound tongue imaging, MRI, lip video). The advantage of lip video is that it is easily available and affordable: most modern smartphones have a front camera. There are already a few solutions for lip-to-speech synthesis, but they mostly concentrate on offline training and inference. In this paper, we propose a system built from a backend for deep neural network training and inference and a fronted as a form of a mobile application. Our initial evaluation shows that the scenario is feasible: a top-5 classification accuracy of 74% is combined with feedback from the mobile application user, making sure that the speaking impaired might be able to communicate with this solution.


title: Star DGT: a Robust Gabor Transform for Speech Denoising

abstract: In this paper, we address the speech denoising problem, where Gaussian, pink and blue additive noises are to be removed from a given speech signal. Our approach is based on a redundant, analysis-sparse representation of the original speech signal. We pick an eigenvector of the Zauner unitary matrix and -- under certain assumptions on the ambient dimension -- we use it as window vector to generate a spark deficient Gabor frame. The analysis operator associated with such a frame, is a (highly) redundant Gabor transform, which we use as a sparsifying transform in denoising procedure. We conduct computational experiments on real-world speech data, using as baseline three Gabor transforms generated by state-of-the-art window vectors in time-frequency analysis and compare their performance to the proposed Gabor transform. The results show that our proposed redundant Gabor transform outperforms all others, consistently for all signals.


title: Search for hybrid morphology radio galaxies from FIRST survey at 1400 MHz

abstract: Hybrid Morphology Radio Sources (HyMoRS) are a very rare subclass of radio galaxies with apparent mixed FR morphology, i.e. these galaxies seem to have an FR-I structure on one side of the core and an FR-II structure on the other side of the core. We systematically searched for HyMoRS using Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-cm (FIRST) survey with 1400 MHz frequency and identified 33 candidate HyMoRS. Our finding significantly increased the known sample size of HyMoRS. HyMoRS may play an essential role in understanding the interaction of jets with the interstellar medium and the much-debated topic of the FR dichotomy. We identified optical counterparts for 29 sources in our catalogue. In our sample of sources, one source (J1106+1355) had quasar-like behaviour. Four sources were BRCLG (Brightest Cluster Galaxies) and six were LRG (Luminous Red Galaxies). We have estimated the spectral index and radio luminosity of HyMoRS in our catalogue, when possible. We found that J1136--0328 was the most luminous source in our sample ($\log L = 27.01$ W Hz$^{-1}$sr$^{-1}$). It was also the farthest HyMoRS (with a redshift $z$ = 0.82) in our sample. With the help of a large sample size of discovered sources, various statistical properties of detected galaxies were studied.


title: Impact of Encoding and Segmentation Strategies on End-to-End Simultaneous Speech Translation

abstract: Boosted by the simultaneous translation shared task at IWSLT 2020, promising end-to-end online speech translation approaches were recently proposed. They consist in incrementally encoding a speech input (in a source language) and decoding the corresponding text (in a target language) with the best possible trade-off between latency and translation quality. This paper investigates two key aspects of end-to-end simultaneous speech translation: (a) how to encode efficiently the continuous speech flow, and (b) how to segment the speech flow in order to alternate optimally between reading (R: encoding input) and writing (W: decoding output) operations. We extend our previously proposed end-to-end online decoding strategy and show that while replacing BLSTM by ULSTM encoding degrades performance in offline mode, it actually improves both efficiency and performance in online mode. We also measure the impact of different methods to segment the speech signal (using fixed interval boundaries, oracle word boundaries or randomly set boundaries) and show that our best end-to-end online decoding strategy is surprisingly the one that alternates R/W operations on fixed size blocks on our English-German speech translation setup.


title: Discrete Transformation Elasticity: An Approach to Design Lattice-based Polar Metamaterials

abstract: The transformation method is a powerful tool for providing the constitutive parameters of the transformed material in the new coordinates. In transformation elasticity, a general curvilinear change of coordinates transforms conventional Hooke's law into a different constitutive law in which the transformed material is not only anisotropic but also polar and chiral and no known elastic solid satisfies. However, this state-of-the-art description provides no insight as to what the underlying microstructure of this transformed material could be, the design of which is a major challenge in this field. The study aims to theoretically justify the fundamental need for the polar material by critically revisiting the discrete transformation method. The key idea is to let transformation gradient operate not only on the elastic properties but on the underlying architectures of the mechanical lattice. As an outstanding application, we leverage the proposed design paradigm to physically construct a polar lattice metamaterial for the observation of elastic carpet cloaking. Numerical simulations are then implemented to show excellent cloaking performance under different static and dynamic mechanical loads. The approach presented herein could promote and accelerate new designs of lattice topologies for transformation elasticity in particular and is able to be extended for realizing other emerging elastic properties and unlocking peculiar functions including statics and dynamics in general.


title: Generalized Abelian Gauge Field Theory under Rotor Model

abstract: Gauge field theory with rank-one field $T_{\mu}$ is a quantum field theory that describes the interaction of elementary spin-1 particles, of which being massless to preserve gauge symmetry. In this paper, we give a generalized, extended study of abelian gauge field theory under successive rotor model in general $D$-dimensional flat spacetime for spin-1 particles in the context of higher order derivatives. We establish a theorem that $n$ rotor contributes to the $\Box^n T^{\mu}$ fields in the integration-by-parts formalism of the action. This corresponds to the transformation of gauge field $T^{\mu} \rightarrow \Box^n T^{\mu}$ and gauge field strength $G_{\mu\nu}\rightarrow \Box^n G_{\mu\nu} $ in the action. The $n=0$ case restores back to the standard abelian gauge field theory. The equation of motion and Noether's conserved current of the theory are also studied.


title: Lusztig correspondence and the Gan-Gross-Prasad problem

abstract: The Gan-Gross-Prasad problem is to describe the restriction of representations of a classical group $G$ to smaller groups $H$ of the same kind. In this paper, we solved the Gan-Gross-Prasad problem over finite fields completely. In previous work \cite{LW1,LW2,LW3,Wang}, we study the Gan-Gross-Prasad problem for unipotent representations of finite classical groups. The main tools used are the Lusztig correspondence as well as a formula of Reeder \cite{R} for the pairings of Deligne-Lusztig characters. We give a reduction decomposition of Reeder's formula, and deduce the Gan-Gross-Prasad problem for arbitrary representations from the unipotent representations by Lusztig correspondence.


title: Learning Hamiltonian dynamics by reservoir computer

abstract: Reconstructing the KAM dynamics diagram of Hamiltonian system from the time series of a limited number of parameters is an outstanding question in nonlinear science, especially when the Hamiltonian governing the system dynamics are unknown. Here, we demonstrate that this question can be addressed by the machine learning approach knowing as reservoir computer (RC). Specifically, we show that without prior knowledge about the Hamilton's equations of motion, the trained RC is able to not only predict the short-term evolution of the system state, but also replicate the long-term ergodic properties of the system dynamics. Furthermore, by the architecture of parameter-aware RC, we also show that the RC trained by the time series acquired at a handful parameters is able to reconstruct the entire KAM dynamics diagram with a high precision by tuning a control parameter externally. The feasibility and efficiency of the learning techniques are demonstrated in two classical nonlinear Hamiltonian systems, namely the double-pendulum oscillator and the standard map. Our study indicates that, as a complex dynamical system, RC is able to learn from data the Hamiltonian.


title: Identifying modulation formats using integrated clustering algorithm

abstract: Modulation format identification (MFI) is crucial in next-generation optical networks such as cognitive optical networks. An integrated-clustering-algorithm-based MFI scheme in a coherent optical communication system is proposed herein. Numerical simulations are performed to test the performance of the scheme on the platform of a polarization domain multiplexing system at a symbol rate of 28 Gbaud. Simulations show that the MFI scheme can achieve an accuracy of 100% for five modulation formats considered in this study when the optical signal-to-noise ratio values are lower than the 7% forward error correction threshold. Tolerances to residual chromatic dispersion (CD) are discussed, and the results show that the proposed scheme can perform well as the residual CD changes over a wide range of values. Moreover, by selecting a number of key blocks, the proposed scheme is less complex than other clustering-based MFI schemes.


title: Reverse Shortest Path Problem for Unit-Disk Graphs

abstract: Given a set P of n points in the plane, the unit-disk graph G_{r}(P) with respect to a parameter r is an undirected graph whose vertex set is P such that an edge connects two points p, q \in P if the Euclidean distance between p and q is at most r (the weight of the edge is 1 in the unweighted case and is the distance between p and q in the weighted case). Given a value \lambda>0 and two points s and t of P, we consider the following reverse shortest path problem: computing the smallest r such that the shortest path length between s and t in G_r(P) is at most \lambda. In this paper, we present an algorithm of O(\lfloor \lambda \rfloor \cdot n \log n) time and another algorithm of O(n^{5/4} \log^{7/4} n) time for the unweighted case, as well as an O(n^{5/4} \log^{5/2} n) time algorithm for the weighted case. We also consider the L_1 version of the problem where the distance of two points is measured by the L_1 metric; we solve the problem in O(n \log^3 n) time for both the unweighted and weighted cases.


title: Analytical formulation for multidimensional continuous opinion models

abstract: Usually, opinion formation models assume that individuals have an opinion about a given topic which can change due to interactions with others. However, individuals can have different opinions in different topics and therefore n-dimensional models are best suited to deal with these cases. While there have been many efforts to develop analytical models for one dimensional opinion models, less attention has been paid to multidimensional ones. In this work, we develop an analytical approach for multidimensional models of continuous opinions where dimensions can be correlated or uncorrelated. We show that for any generic reciprocal interactions between agents, the mean value of initial opinion distribution is conserved. Moreover, for positive social influence interaction mechanisms, the variance of opinion distributions decreases with time and the system converges to a delta distributed function. In particular, we calculate the convergence time when agents get closer in a discrete quantity after interacting, showing a clear difference between correlated and uncorrelated cases.