title: Sequential construction and dimension reduction of Gaussian processes under constraints

abstract: Accounting for inequality constraints, such as boundedness, monotonicity or convexity, is challenging when modeling costly-to-evaluate black box functions. In this regard, finite-dimensional Gaussian process (GP) regression models bring a valuable solution, as they guarantee that the inequality constraints are satisfied everywhere. Nevertheless, these models are currently restricted to small dimensional situations (up to dimension 5). Addressing this issue, we introduce the MaxMod algorithm that sequentially inserts one-dimensional knots or adds active variables, thereby performing at the same time dimension reduction and efficient knot allocation. We prove the convergence of this algorithm. In intermediary steps of the proof, we propose the notion of multi-affine extension and study its properties. We also prove the convergence of finite-dimensional GPs, when the knots are not dense in the input space, extending the recent literature. With simulated and real data, we demonstrate that the MaxMod algorithm remains efficient in higher dimension (at least in dimension 20), and needs fewer knots than other constrained GP models from the state-of-the-art, to reach a given approximation error.


title: Analysis of a cross-diffusion model for rival gangs interaction in a city

abstract: We study a two-species cross-diffusion model that is inspired by a system of convection-diffusion equations derived from an agent-based model on a two-dimensional discrete lattice. The latter model has been proposed to simulate gang territorial development through the use of graffiti markings. We find two energy functionals for the system that allow us to prove a weak-stability result and identify equilibrium solutions. We show that under the natural definition of weak solutions, obtained from the weak-stability result, the system does not allow segregated solutions. Moreover, we present a result on the long-term behavior of solutions in the case when the product of the masses of the densities are smaller than a critical value. This result is complemented with numerical experiments.


title: Analysis of Hand-Crafted and Automatic-Learned Features for Glaucoma Detection Through Raw Circmpapillary OCT Images

abstract: Taking into account that glaucoma is the leading cause of blindness worldwide, we propose in this paper three different learning methodologies for glaucoma detection in order to elucidate that traditional machine-learning techniques could outperform deep-learning algorithms, especially when the image data set is small. The experiments were performed on a private database composed of 194 glaucomatous and 198 normal B-scans diagnosed by expert ophthalmologists. As a novelty, we only considered raw circumpapillary OCT images to build the predictive models, without using other expensive tests such as visual field and intraocular pressure measures. The results ratify that the proposed hand-driven learning model, based on novel descriptors, outperforms the automatic learning. Additionally, the hybrid approach consisting of a combination of both strategies reports the best performance, with an area under the ROC curve of 0.85 and an accuracy of 0.82 during the prediction stage.


title: Evidence of the four-quark nature of $f_0(980)$ and $f_0(500)$

abstract: There exists a great deal of concrete evidence in favor of the exotic four-quark nature of light scalars. At the same time, the further expansion of the area of the $q^2\bar q^2$ model validity for light scalars on ever new processes seems extremely interesting and important. We analyze the BESIII data on the decay $J/\psi\to \gamma\pi^0\pi^0$ and show that the results of this high-statistics experiment can be interpreted in favor of the four-quark nature of light scalar mesons $f_0(980)$ and $f_0(500)$.


title: Brotate and Tribike: Designing Smartphone Control for Cycling

abstract: The more people commute by bicycle, the higher is the number of cyclists using their smartphones while cycling and compromising traffic safety. We have designed, implemented and evaluated two prototypes for smartphone control devices that do not require the cyclists to remove their hands from the handlebars - the three-button device Tribike and the rotation-controlled Brotate. The devices were the result of a user-centred design process where we identified the key features needed for a on-bike smartphone control device. We evaluated the devices in a biking exercise with 19 participants, where users completed a series of common smartphone tasks. The study showed that Brotate allowed for significantly more lateral control of the bicycle and both devices reduced the cognitive load required to use the smartphone. Our work contributes insights into designing interfaces for cycling.


title: A Broad Absorption Line Outflow Associated with the Broad Emission Line Region in the Quasar SDSS J075133.35+134548.3

abstract: We report the discovery of unusual broad absorption lines (BALs) in the bright quasar SDSS J075133.35+134548.3 at z~1, and present a detailed study of BAL and emission lines using the archived and newly obtained optical and NIR spectroscopic data. The BALs are detected tentatively in Mg II, Al III and reliably in HeI*3889, HeI*10830, which show complexed velocity structures splitting into two components: a high-velocity component (HV), with a blueshifted velocity range of -9300 ~ -3500 km/s, can be reliably detected tentatively in Mg II, Al III and reliably in HeI* 10830, whereas it is undetectable in HeI* 3889; and low-Velocity component (LV), with the velocity of -3500 ~ -1800 km/s, is only detected in HeI*3889 and HeI* 10830. With the BALs from different ions, the HV outflowing gas can be constrained to have a density n_H~10^10.3-10^11.4 cm^-3, a column density N_H~10^21 cm^-2 and an ionization parameter U~10^-1.83 - 10^-1.72; and thus inferring a distance of R_HV~0.5 pc from the central continuum source with a monochromatic luminosity of 7.0*10^45 erg/s at 5100 A. This distance is remarkable similar to that of normal broad emission line region (BLR) estimated according to the luminosity-distance from the reverberation experiments, suggesting the association between the BLR and the HV BAL outflowing gas. Interestingly enough, a blueshifted component is also significantly detected in Mg II and Al III broad emission lines (BELs), and the Al III/Mg II of such a BEL component can be reproduced by the physical parameters inferred from the HV BAL gas. The LV BAL gas very likely have a larger column density and a higher ionization level than the HV BAL gas, and hence a smaller distance than that of the latter. Further spectroscopy with high S/N ratio and broader wavelength coverage is needed to confirm this and shed new light on the possible connection between BALs and BELs.


title: Spectral functions from the real-time functional renormalization group

abstract: We employ the functional renormalization group approach formulated on the Schwinger-Keldysh contour to calculate real-time correlation functions in scalar field theories. We provide a detailed description of the formalism, discuss suitable truncation schemes for real-time calculations as well as the numerical procedure to self-consistently solve the flow equations for the spectral function. Subsequently, we discuss the relations to other perturbative and non-perturbative approaches to calculate spectral functions, and present a detailed comparison and benchmark in $d=0+1$ dimensions.


title: Bifurcation analysis of the Hardy-Sobolev equation

abstract: In this paper, we prove existence of multiple non-radial solutions to the Hardy-Sobolev equation $$\begin{cases} -\Delta u-\displaystyle\frac \gamma{|x|^2}u=\displaystyle\frac{1}{|x|^s}|u|^{p_s-2}u & \text{ in } \mathbb{R}^N\setminus\{0\},\\ u\geq 0, & \end{cases}$$ where $N\geq 3$, $s\in[0,2)$, $p_s=\frac{2(N-s)}{N-2}$ and $\gamma\in (-\infty,\frac{(N-2)^2} 4)$. We extend results of E.N. Dancer, F. Gladiali, M. Grossi, Proc. Roy. Soc. Edinburgh Sect. A 147 (2017) where only the case $s=0$ is considered. Moreover, thanks to monotonicity properties of the solutions, we separate two branches of non-radial solutions.


title: Magnetically induced metal-insulator transition in Pb2CaOsO6

abstract: We report on the structural, magnetic, and electronic properties of two new double-perovskites synthesized under high pressure; Pb2CaOsO6 and Pb2ZnOsO6. Upon cooling below 80 K, Pb2CaOsO6 simultaneously undergoes a metal--insulator transition and develops antiferromagnetic order. Pb2ZnOsO6, on the other hand, remains a paramagnetic metal down to 2 K. The key difference between the two compounds lies in their crystal structure. The Os atoms in Pb2ZnOsO6 are arranged on an approximately face-centred cubic lattice with strong antiferromagnetic nearest-neighbor exchange couplings. The geometrical frustration inherent to this lattice prevents magnetic order from forming down to the lowest temperatures. In contrast, the unit cell of Pb2CaOsO6 is heavily distorted up to at least 500 K, including antiferroelectric-like displacements of the Pb and O atoms despite metallic conductivity above 80 K. This distortion relieves the magnetic frustration, facilitating magnetic order which in turn drives the metal--insulator transition. Our results suggest that the phase transition in Pb2CaOsO6 is spin-driven, and could be a rare example of a Slater transition.


title: QR-MIX: Distributional Value Function Factorisation for Cooperative Multi-Agent Reinforcement Learning

abstract: In Cooperative Multi-Agent Reinforcement Learning (MARL) and under the setting of Centralized Training with Decentralized Execution (CTDE), agents observe and interact with their environment locally and independently. With local observation and random sampling, the randomness in rewards and observations leads to randomness in long-term returns. Existing methods such as Value Decomposition Network (VDN) and QMIX estimate the value of long-term returns as a scalar that does not contain the information of randomness. Our proposed model QR-MIX introduces quantile regression, modeling joint state-action values as a distribution, combining QMIX with Implicit Quantile Network (IQN). However, the monotonicity in QMIX limits the expression of joint state-action value distribution and may lead to incorrect estimation results in non-monotonic cases. Therefore, we proposed a flexible loss function to approximate the monotonicity found in QMIX. Our model is not only more tolerant of the randomness of returns, but also more tolerant of the randomness of monotonic constraints. The experimental results demonstrate that QR-MIX outperforms the previous state-of-the-art method QMIX in the StarCraft Multi-Agent Challenge (SMAC) environment.


title: Mobile Apps for Foot Measurement: A Scoping Review

abstract: With the proliferation of smart phone, a major growth in the use of apps related to the health category, specifically those concerned with foot health can be observed. Although new, these apps are being used practically for scanning feet with an aim to providing accurate information about various properties of the human foot. With the availability of many 'foot scanning and measuring apps' in the app stores, the need for an evaluation system for such apps can be deemed necessary as little information regarding the evidence-based quality of these apps is available. To characterize the assessment of measurement techniques and essential software quality characteristics of mobile foot measuring apps, and determine their effectiveness for potential use as commercial professional tools for foot care health professionals such as pedorthists, podiatrists, orthotists and so on, to assist in measuring foot for custom shoes, and for individuals to enhance the awareness of foot health and hygiene and prevention of foot-related problems. An electronic search across Android and iOS app stores was conducted between July 2020 and August 2020 for apps related to foot measurement. Mobile apps with stated goals of foot measurement and general foot health were identified and selected by three independent raters and discrepancies regarding the selected apps were resolved via a fourth rater. Evaluation inferences found all apps failing to meet even half of the measurement-specific criteria required for the proper manufacturing of custom-made footwear. 23% (6/26) of apps were found to utilize either external scanners or advanced algorithms to reconstruct 3D models of user foot that can possibly be used for ordering custom-made footwear and medical casts for fitting irregular foot sizes and shapes. Overall, current apps for foot measurement do not follow any specific guidelines for measurement purposes.


title: Optimizing BLE-Like Neighbor Discovery

abstract: Neighbor discovery (ND) protocols are used for establishing a first contact between multiple wireless devices. The energy consumption and discovery latency of this procedure are determined by the parametrization of the protocol. In most existing protocols, reception and transmission are temporally coupled. Such schemes are referred to as \textit{slotted}, for which the problem of finding optimized parametrizations has been studied thoroughly in the literature. However, slotted approaches are not efficient in applications in which new devices join the network gradually and only the joining devices and a master node need to run the ND protocol simultaneously. For example, this is typically the case in IoT scenarios or Bluetooth Low Energy (BLE) piconets. Here, protocols in which packets are transmitted with periodic intervals (PI) can achieve significantly lower worst-case latencies than slotted ones. For this class of protocols, optimal parameter values remain unknown. To address this, we propose an optimization framework for PI-based BLE-like protocols, which translates any specified duty-cycle (and therefore energy budget) into a set of optimized parameter values. We show that the parametrizations resulting from one variant of our proposed scheme are optimal when one receiver discovers one transmitter, and no other parametrization or ND protocol - neither slotted nor slotless - can guarantee lower discovery latencies for a given duty-cycle in this scenario. Since the resulting protocol utilizes the channel more aggressively than other ND protocols, beacons will collide more frequently. Hence, due to collisions, the rate of successful discoveries gracefully decreases for larger numbers of devices discovering each other simultaneously. We also propose a scheme for configuring the BLE protocol (and not just BLE-\textit{like} protocols).


title: Rise of the Machines? Intraday High-Frequency Trading Patterns of Cryptocurrencies

abstract: This research analyses high-frequency data of the cryptocurrency market in regards to intraday trading patterns related to algorithmic trading and its impact on the European cryptocurrency market. We study trading quantitatives such as returns, traded volumes, volatility periodicity, and provide summary statistics of return correlations to CRIX (CRyptocurrency IndeX), as well as respective overall high-frequency based market statistics with respect to temporal aspects. Our results provide mandatory insight into a market, where the grand scale employment of automated trading algorithms and the extremely rapid execution of trades might seem to be a standard based on media reports. Our findings on intraday momentum of trading patterns lead to a new quantitative view on approaching the predictability of economic value in this new digital market.


title: Magnetic fields in star-forming systems (II): examining dust polarization, the Zeeman effect, and the Faraday rotation measure as magnetic field tracers

abstract: The degree to which the formation and evolution of clouds and filaments in the interstellar medium is regulated by magnetic fields remains an open question. Yet the fundamental properties of the fields (strength and 3D morphology) are not readily observable. We investigate the potential for recovering magnetic field information from dust polarization, the Zeeman effect, and the Faraday rotation measure ($RM$) in a SILCC-Zoom magnetohydrodynamic (MHD) filament simulation. The object is analyzed at the onset of star formation, and it is characterized by a line-mass of about M/L $\sim 63\ M_{\odot}\ pc^{-1}$ out to a radius of $1\,$pc and a kinked 3D magnetic field morphology. We generate synthetic observations via POLARIS radiative transfer (RT) post-processing, and compare with an analytical model of helical or kinked field morphology to help interpreting the inferred observational signatures. We show that the tracer signals originate close to the filament spine. We find regions along the filament where the angular-dependency with the line-of-sight (LOS) is the dominant factor and dust polarization may trace the underlying kinked magnetic field morphology. We also find that reversals in the recovered magnetic field direction are not unambiguously associated to any particular morphology. Other physical parameters, such as density or temperature, are relevant and sometimes dominant compared to the magnetic field structure in modulating the observed signal. We demonstrate that the Zeeman effect and the $RM$ recover the line-of-sight magnetic field strength to within a factor 2.1 - 3.4. We conclude that the magnetic field morphology may not be unambiguously determined in low-mass systems by observations of dust polarization, Zeeman effect, or $RM$, whereas the field strengths can be reliably recovered.


title: Impact of News on the Commodity Market: Dataset and Results

abstract: Over the last few years, machine learning based methods have been applied to extract information from news flow in the financial domain. However, this information has mostly been in the form of the financial sentiments contained in the news headlines, primarily for the stock prices. In our current work, we propose that various other dimensions of information can be extracted from news headlines, which will be of interest to investors, policy-makers and other practitioners. We propose a framework that extracts information such as past movements and expected directionality in prices, asset comparison and other general information that the news is referring to. We apply this framework to the commodity "Gold" and train the machine learning models using a dataset of 11,412 human-annotated news headlines (released with this study), collected from the period 2000-2019. We experiment to validate the causal effect of news flow on gold prices and observe that the information produced from our framework significantly impacts the future gold price.


title: Tactical Decision Making for Emergency Vehicles Based on A Combinational Learning Method

abstract: Increasing the response time of emergency vehicles(EVs) could lead to an immeasurable loss of property and life. On this account, tactical decision making for EVs' microscopic control remains an indispensable issue to be improved. In this paper, a rule-based avoiding strategy(AS) is devised, that CVs in the prioritized zone ahead of EV should accelerate or change their lane to avoid it. Besides, a novel DQN method with speed-adaptive compact state space (SC-DQN) is put forward to fit in EVs' high-speed feature and generalize in various road topologies. Afterward, the execution of AS feedback to the input of SC-DQN so that they joint organically as a combinational method. The following approach reveals that DRL could complement rule-based avoiding strategy in generalization, and on the contrary, the rule-based avoiding strategy could complement DRL in stability, and their combination could lead to less response time, lower collision rate and smoother trajectory.


title: A model-independent constraint on the Hubble constant with gravitational waves from the Einstein Telescope

abstract: In this paper, we investigate the expected constraints on the Hubble constant from the gravitational-wave standard sirens, in a cosmological-model-independent way. In the framework of the well-known Hubble law, the GW signal from each detected binary merger in the local universe ($z<0.10$) provides a measurement of luminosity distance $D_L$ and thus the Hubble constant $H_0$. Focusing on the simulated data of gravitational waves from the third-generation gravitational wave detector (the Einstein Telescope, ET), combined with the redshifts determined from electromagnetic counter parts and host galaxies, one can expect the Hubble constant to be constrained at the precision of $\sim 10^{-2}$ with 20 well-observed binary neutron star (BNS) mergers. Additional standard-siren measurements from other types of future gravitational-wave sources (NS-BH and BBH) will provide more precision constraints of this important cosmological parameter. Therefore, we obtain that future measurements of the luminosity distances of gravitational waves sources will be much more competitive than the current analysis, which makes it expectable more vigorous and convincing constraints on the Hubble constant in a cosmological-model-independent way.


title: Local probe of the interlayer coupling strength of few-layers SnSe by contact-resonance atomic force microscopy

abstract: The interlayer bonding in two dimensional materials is particularly important because it is not only related to their physical and chemical stability but also affects their mechanical, thermal, electronic, optical, and other properties. To address this issue, we report the direct characterization of the interlayer bonding in 2D SnSe using contact-resonance atomic force microscopy in this study. Site specific CR spectroscopy and CR force spectroscopy measurements are performed on both SnSe and its supporting SiO2 substrate comparatively. Based on the cantilever and contact mechanic models, the contact stiffness and vertical Young's modulus are evaluated in comparison with SiO2 as a reference material. The interlayer bonding of SnSe is further analyzed in combination with the semi-analytical model and density functional theory calculations. The direct characterization of interlayer interactions using this nondestructive methodology of CR AFM would facilitate a better understanding of the physical and chemical properties of 2D layered materials, specifically for interlayer intercalation and vertical heterostructures.


title: Anharmonicity in Raman-active phonon modes in atomically thin MoS$_2$

abstract: Phonon-phonon anharmonic effects have a strong influence on the phonon spectrum; most prominent manifestation of these effects are the softening (shift in frequency) and broadening (change in FWHM) of the phonon modes at finite temperature. Using Raman spectroscopy, we studied the temperature dependence of the FWHM and Raman shift of $\mathrm{E_{2g}^1}$ and $\mathrm{A_{1g}}$ modes for single-layer and natural bilayer MoS$_2$ over a broad range of temperatures ($8 < $T$ < 300$ K). Both the Raman shift and FWHM of these modes show linear temperature dependence for $T>100$ K, whereas they become independent of temperature for $T<100$ K. Using first-principles calculations, we show that three-phonon anharmonic effects intrinsic to the material can account for the observed temperature-dependence of the line-width of both the modes. It also plays an important role in determining the temperature-dependence of the frequency of the Raman modes. The observed evolution of the line-width of the A$_{1g}$ mode suggests that electron-phonon processes are additionally involved. From the analysis of the temperature-dependent Raman spectra of MoS$_2$ on two different substrates -- SiO$_2$ and hexagonal boron nitride, we disentangle the contributions of external stress and internal impurities to these phonon-related processes. We find that the renormalization of the phonon mode frequencies on different substrates is governed by strain and intrinsic doping. Our work establishes the role of intrinsic phonon anharmonic effects in deciding the Raman shift in MoS$_2$ irrespective of substrate and layer number.


title: A Security Architecture for Railway Signalling

abstract: We present the proposed security architecture Deutsche Bahn plans to deploy to protect its trackside safety-critical signalling system against cyber-attacks. We first present the existing reference interlocking system that is built using standard components. Next, we present a taxonomy to help model the attack vectors relevant for the railway environment. Building upon this, we present the proposed "compartmentalized" defence concept for securing the upcoming signalling systems.


title: Box dimensions of $(\times m, \times n)$-invariant sets

abstract: We study the box dimensions of sets invariant under the toral endomorphism $(x, y) \mapsto (m x \text{ mod } 1, \, n y \text{ mod } 1)$ for integers $n>m \geq 2$. The basic examples of such sets are Bedford-McMullen carpets and, more generally, invariant sets are modelled by subshifts on the associated symbolic space. When this subshift is topologically mixing and sofic the situation is well-understood by results of Kenyon and Peres. Moreover, other work of Kenyon and Peres shows that the Hausdorff dimension is generally given by a variational principle. Therefore, our work is focused on the box dimensions in the case where the underlying shift is not topologically mixing and sofic. We establish straightforward upper and lower bounds for the box dimensions in terms of entropy which hold for all subshifts and show that the upper bound is the correct value for coded subshifts whose entropy can be realised by words which can be freely concatenated, which includes many well-known families such as $\beta$-shifts, (generalised) $S$-gap shifts, and transitive sofic shifts. We also provide examples of transitive coded subshifts where the general upper bound fails and the box dimension is actually given by the general lower bound. In the non-transitive sofic setting, we provide a formula for the box dimensions which is often intermediate between the general lower and upper bounds.


title: 3D printed waveguides based on Photonic Crystal Fiber designs for complex fiber-end photonic devices

abstract: Optical waveguide segments based on geometrically unbound photonic crystal fibers (PCF) designs could be exploited as building blocks to realize miniaturized complex devices which implement advanced photonic operations. Here, we show how to fabricate optical waveguide segments with PCF designs by direct high-resolution 3D printing and how the combination of these segments can realise complex photonic devices. We demonstrate the unprecedented precision and flexibility of our method by fabricating the first-ever fiber polarising beam splitter based on PCFs. The device was directly printed in one step on the end-face of a standard single-mode fiber and was 210~$\mu$m-long, offering broadband operation in the optical telecommunications C-band. Our approach harnesses the potential of high-resolution 3D printing and of PCF designs paving the way for the development of novel miniaturised complex photonic systems, which will positively impact and advance optical telecommunications, sensor technology, and biomedical devices.


title: Strain energy function for isotropic non-linear elastic incompressible solids with linear finite strain response in shear and torsion

abstract: We find the strain energy function for isotropic incompressible solids exhibiting a linear relationship between shear stress and amount of shear, and between torque and amount of twist, when subject to large simple shear or torsion deformations. It is inclusive of the well-known neo-Hookean and the Mooney-Rivlin models, but also can accommodate other terms, as certain arbitrary functions of the principal strain invariants. Effectively, the extra terms can be used to account for several non-linear effects observed experimentally but not captured by the neo-Hookean and Mooney-Rivlin models, such as strain stiffening effects due to limiting chain extensibility.


title: Domain of Influence analysis: implications for Data Assimilation in space weather forecasting

abstract: Solar activity, ranging from the background solar wind to energetic coronal mass ejections (CMEs), is the main driver of the conditions in the interplanetary space and in the terrestrial space environment, known as space weather. A better understanding of the Sun-Earth connection carries enormous potential to mitigate negative space weather effects with economic and social benefits. Effective space weather forecasting relies on data and models. In this paper, we discuss some of the most used space weather models, and propose suitable locations for data gathering with space weather purposes. We report on the application of \textit{Representer analysis (RA)} and \textit{Domain of Influence (DOI) analysis} to three models simulating different stages of the Sun-Earth connection: the OpenGGCM and Tsyganenko models, focusing on solar wind - magnetosphere interaction, and the PLUTO model, used to simulate CME propagation in interplanetary space. Our analysis is promising for space weather purposes for several reasons. First, we obtain quantitative information about the most useful locations of observation points, such as solar wind monitors. For example, we find that the absolute values of the DOI are extremely low in the magnetospheric plasma sheet. Since knowledge of that particular sub-system is crucial for space weather, enhanced monitoring of the region would be most beneficial. Second, we are able to better characterize the models. Although the current analysis focuses on spatial rather than temporal correlations, we find that time-independent models are less useful for Data Assimilation activities than time-dependent models. Third, we take the first steps towards the ambitious goal of identifying the most relevant heliospheric parameters for modelling CME propagation in the heliosphere, their arrival time, and their geoeffectiveness at Earth.


title: Using the Navier-Cauchy equation for motion estimation in dynamic imaging

abstract: Tomographic image reconstruction is well understood if the specimen being studied is stationary during data acquisition. However, if this specimen changes during the measuring process, standard reconstruction techniques can lead to severe motion artefacts in the computed images. Solving a dynamic reconstruction problem therefore requires to model and incorporate suitable information on the dynamics in the reconstruction step to compensate for the motion. Many dynamic processes can be described by partial differential equations which thus could serve as additional information for the purpose of motion compensation. In this article, we consider the Navier-Cauchy equation which characterizes small elastic deformations and serves, for instance, as a model for respiratory motion. Our goal is to provide a proof-of-concept that by incorporating the deformation fields provided by this PDE, one can reduce the respective motion artefacts in the reconstructed image. To this end, we solve the Navier-Cauchy equation prior to the image reconstruction step using suitable initial and boundary data. Then, the thus computed deformation fields are incorporated into an analytic dynamic reconstruction method to compute an image of the unknown interior structure. The feasibility is illustrated with numerical examples from computerized tomography.


title: Analysis of the least sum-of-minimums estimator for switched systems

abstract: This paper considers a particular parameter estimator for switched systems and analyzes its properties. The estimator in question is defined as the map from the data set to the solution set of an optimization problem where the to-be-optimized cost function is a sum of pointwise infima over a finite set of sub-functions. This is a hard nonconvex problem. The paper studies some fundamental properties of this problem such as uniqueness of the solution or boundedness of the estimation error regardless of computational considerations. The interest of the analysis is to lay out the main influential properties of the data on the performance of this (ideal) estimator.


title: Tracing circumnuclear dense gas in H2O maser galaxies

abstract: A sample of 30 H2O extra-galactic maser galaxies with their published HCN(J=1-0) and HCO+(J=1-0) observations has been compiled to investigate the dense gas correlation with H2O maser emission. Our sample number exceeds the size of the previous HCN samples studied so far by a factor of three, and it is the first study on the possible relation with the dense gas tracer HCO+. We find a strong correlation between normalized H2O maser emission luminosity (LH2O\LCO) and normalized HCO+ line luminosity (LHCO+\LCO). Moreover, a weak correlation has been found between LH2O\LCO and normalized HCN line luminosity (LHCN\LCO). The sample is also studied after excluding Luminous and Ultraluminous infrared galaxy (U)LIRG sources, and the mentioned correlations are noticeably stronger. We show that 'Dense gas' fractions as obtained from HCN and HCO+ molecules tightly correlate with maser emission, especially for galaxies with normal IR luminosity(LIR< 10^11Lsun) and we show that HCO+ is a better 'dense gas' tracer than HCN. Further systematic studies of these dense gas tracers with higher transition level lines are vital to probe megamaser physical conditions and to accurately determining how maser emission interrelates with the dense gas.


title: Unmanned Aerial Vehicle Control Through Domain-based Automatic Speech Recognition

abstract: Currently, unmanned aerial vehicles, such as drones, are becoming a part of our lives and reaching out to many areas of society, including the industrialized world. A common alternative to control the movements and actions of the drone is through unwired tactile interfaces, for which different remote control devices can be found. However, control through such devices is not a natural, human-like communication interface, which sometimes is difficult to master for some users. In this work, we present a domain-based speech recognition architecture to effectively control an unmanned aerial vehicle such as a drone. The drone control is performed using a more natural, human-like way to communicate the instructions. Moreover, we implement an algorithm for command interpretation using both Spanish and English languages, as well as to control the movements of the drone in a simulated domestic environment. The conducted experiments involve participants giving voice commands to the drone in both languages in order to compare the effectiveness of each of them, considering the mother tongue of the participants in the experiment. Additionally, different levels of distortion have been applied to the voice commands in order to test the proposed approach when facing noisy input signals. The obtained results show that the unmanned aerial vehicle is capable of interpreting user voice instructions achieving an improvement in speech-to-action recognition for both languages when using phoneme matching in comparison to only using the cloud-based algorithm without domain-based instructions. Using raw audio inputs, the cloud-based approach achieves 74.81% and 97.04% accuracy for English and Spanish instructions respectively, whereas using our phoneme matching approach the results are improved achieving 93.33% and 100.00% accuracy for English and Spanish languages.


title: Beyond the brain: towards a mathematical modeling of emotions

abstract: Emotions are a central key for understanding human beings and of fundamental importance regarding their impact in human and animal behaviors. They have been for a long time a subject of study for various scholars including in particular philosophers and mystics. In modern science, the emotional phenomenon has attracted for a few decades an increasing number of studies, notably in the fields of Psychology, Psychiatry, Neuroscience and Biochemistry. However, since our perception of emotions is not, so far, directly detectable nor recordable by our measure instruments, Physics and Mathematics have not been so far used academically to provide a precise description of the phenomenon of feeling an emotion. Relying upon the works of O. Elahi and on the hypothesis that the human soul and its psyche may manifest in ourselves (in both conscious and unconscious manner) in an analog way as electromagnetic waves, we propose here a few mathematical descriptions consistent with the human personal experience, of the feeling and cognition of emotions. As far as we know, such a mathematical description has never been provided before. It allows a quantitative (intensity) and qualitative (nature of feelings/frequency) of the emotional phenomenon which provides a novel scientific approach of the nature of the mind, complementary to the on going research of physiological manifestation of emotions. We anticipate such an approach and the associated mathematical modeling to become an important tool to describe emotions and their subsequent behavior. In complement of the modeling of oscillations and brain dynamics, it provides a fruitful direction of research with potentially broad and deep impacts in both applied mathematics, physics, cognitive and behavioral sciences.


title: Towards Communication in a Curved Spacetime Geometry

abstract: The current race in quantum communication -- endeavouring to establish a global quantum network -- must account for special and general relativistic effects. The well-studied general relativistic effects include Shapiro time-delay, gravitational lensing, and frame dragging which all are due to how a mass distribution alters geodesics. Here, we report how the curvature of spacetime geometry affects the propagation of information carriers along an arbitrary geodesic. An explicit expression for the distortion onto the carrier wavefunction in terms of the Riemann curvature is obtained. Furthermore, we investigate this distortion for anti-de Sitter and Schwarzschild geometries. For instance, the spacetime curvature causes a 0.10~radian phase-shift for communication between Earth and the International Space Station on a monochromatic laser beam and quadrupole astigmatism can cause a 12.2 % cross-talk between structured modes traversing through the solar system. Our finding shows that this gravitational distortion is significant, and it needs to be either pre- or post-corrected at the sender or receiver to retrieve the information.


title: Phase Transition and Thermal Fluctuations of Quintessential Kerr-Newman-AdS Black Hole

abstract: This paper is devoted to analyzing the critical phenomenon and phase transition of quintessential Kerr-Newman-anti-de Sitter black hole in the framework of Maxwell equal-area law. For this purpose, we first derive thermodynamic quantities such as Hawking temperature, entropy and angular momentum in the context of extended phase space. These quantities satisfy Smarr-Gibbs-Dehum relation in the presence of quintessence matter. We then discuss the critical behavior of thermodynamic quantities through two approaches, i.e., van der Waals-like equation of state and Maxwell equal-area law. It is found that the latter approach is more effective to analyze the critical behavior of the complicated black holes. Using equal-area law, we also study phase diagram in $T-S$ plane and find an isobar which shows the coexistence region of two phases. We conclude that below the critical temperature, black holes show a similar phase transition as that of van der Waals fluid. Finally, we study the effects of thermal fluctuations on the stability of this black hole.


title: Labeled Mereological Set Theory

abstract: The methodology used here might provide a neat method of examining paradoxes and ways to circumvent them. Most of the known set theoretic paradoxes (Russell's, Cantor's, Burali-Forti's,..) can be paralleled here and examined. This account will shed the light on a particular application of this method that appears to elude paradoxes; an application that have shortcomings that will be illustrated here and suggestions to solutions are proposed. First I'll present the exposition of the theory, and then speak about its background and the aims behind it and how it can extend our knowledge of overcoming paradoxes.


title: Non-Bloch band theory and bulk-edge correspondence in non-Hermitian systems

abstract: In this paper, we review our non-Bloch band theory in one-dimensional non-Hermitian tight-binding systems. In our theory, it is shown that in non-Hermitian systems, the Brillouin zone is determined so as to reproduce continuum energy bands in a large open chain. By using simple models, we explain the concept of the non-Bloch band theory and the method to calculate the Brillouin zone. In particular, for the non-Hermitian Su-Schrieffer-Heeger model, the bulk-edge correspondence can be established between the topological invariant defined from our theory and existence of the topological edge states.


title: The effect of an exterior electric field on the instability of dielectric plates

abstract: We investigate the theoretical nonlinear response, Hessian stability, and possible wrinkling behaviour of a voltage-activated dielectric plate immersed in a tank filled with silicone oil. Fixed rigid electrodes are placed on the top and bottom of the tank, and an electric field is generated by a potential difference between the electrodes. We solve the associated incremental boundary value problem of superimposed, inhomogeneous small-amplitude wrinkles, signalling the onset of instability. We decouple the resulting bifurcation equation into symmetric and antisymmetric modes. For a neo-Hookean dielectric plate, we show that a potential difference between the electrodes can induce a thinning of the plate and thus an increase of its planar area, similar to the scenarios encountered when there is no silicone oil. However, we also find that, depending on the material and geometric parameters, an increasing applied voltage can also lead to a thickening of the plate, and thus a shrinking of its area. In that scenario, Hessian instability and wrinkling bifurcation may then occur spontaneously once some critical voltages are reached.


title: Integrated Photodynamic Raman Theranostics for Cancer Diagnosis, Treatment, and Post-Treatment Molecular Monitoring

abstract: Theranostics, the combination of diagnosis and therapy, has long held promise as a means to achieving personalised precision cancer treatments. However, despite its potential, theranostics has yet to realise significant clinical translation, largely due the complexity and overriding toxicity concerns of existing theranostic nanoparticle strategies. Here, we present an alternative nanoparticle-free theranostic approach based on simultaneous Raman spectroscopy and photodynamic therapy (PDT) in an integrated clinical platform for cancer theranostics. We detail the compatibility of Raman spectroscopy and PDT for cancer theranostics, whereby Raman spectroscopic diagnosis can be performed on PDT photosensitiser-positive cells and tissues without inadvertent photosensitiser activation/photobleaching or impaired diagnostic capacity. We further demonstrate that our theranostic platform enables in vivo tumour diagnosis, treatment, and post-treatment molecular monitoring in real-time. This system thus achieves effective theranostic performance, providing a promising new avenue towards the clinical realisation of theranostics.


title: Twisted differential K-characters and D-branes

abstract: We analyse in detail the language of partially non-abelian Deligne cohomology and of twisted differential K-theory, in order to describe the geometry of type II superstring backgrounds with D-branes. This description will also provide the opportunity to show some mathematical results of independent interest. In particular, we begin classifying the possible gauge theories on a D-brane or on a stack of D-branes using the intrinsic tool of long exact sequences. Afterwards, we recall how to construct two relevant models of differential twisted K-theory, paying particular attention to the dependence on the twisting cocycle within its cohomology class. In this way we will be able to define twisted K-homology and twisted Cheeger-Simons K-characters in the category of simply-connected manifolds, eliminating any unnatural dependence on the cocycle. The ambiguity left for non simply-connected manifolds will naturally correspond to the ambiguity in the gauge theory, following the previous classification. This picture will allow for a complete characterization of D-brane world-volumes, the Wess-Zumino action and topological D-brane charges within the K-theoretical framework, that can be compared step by step to the old cohomological classification. This has already been done for backgrounds with vanishing B-field; here we remove this hypothesis.

categories: hep-th math-ph math.AT math.DG math.KT math.MP


title: AoI Minimization in Status Update Control with Energy Harvesting Sensors

abstract: Information freshness is crucial for time-critical IoT applications, e.g., monitoring and control systems. We consider an IoT status update system with multiple users, multiple energy harvesting sensors, and a wireless edge node. The users receive time-sensitive information about physical quantities, each measured by a sensor. Users send requests to the edge node where a cache contains the most recently received measurements from each sensor. To serve a request, the edge node either commands the sensor to send a status update or retrieves the aged measurement from the cache. We aim at finding the best actions of the edge node to minimize the age of information of the served measurements. We model this problem as a Markov decision process and develop reinforcement learning (RL) algorithms: model-based value iteration and model-free Q-learning methods. We also propose a Q-learning method for the realistic case where the edge node is informed about the sensors' battery levels only via the status updates. The case under transmission limitations is also addressed. Furthermore, properties of an optimal policy are analytically characterized. Simulation results show that an optimal policy is a threshold-based policy and that the proposed RL methods significantly reduce the average cost compared to several baselines.


title: Electrostatically tunable axisymmetric vibrations of soft electro-active tubes

abstract: Due to their unique electromechanical coupling properties, soft electro-active (SEA) resonators are actively tunable, extremely suitable, and practically important for designing the next-generation acoustic and vibration treatment devices. In this paper, we investigate the electrostatically tunable axisymmetric vibrations of SEA tubes with different geometric sizes. We consider both axisymmetric torsional and longitudinal vibrations for an incompressible SEA cylindrical tube under inhomogeneous biasing fields induced by radial electric voltage and axial pre-stretch. We then employ the state-space method, which combines the state-space formalism in cylindrical coordinates with the approximate laminate technique, to derive the frequency equations for two separate classes of axisymmetric vibration of the tube subjected to appropriate boundary conditions. We perform numerical calculations to validate the convergence and accuracy of the state-space method and to illuminate that the axisymmetric vibration characteristics of SEA tubes may be tuned significantly by adjusting the electromechanical biasing fields as well as altering the tube geometry. The reported results provide a solid guidance for the proper design of tunable resonant devices composed of SEA tubes


title: Experimental Observation of Vortex Rings in a Bulk Magnet

abstract: Vortex rings are remarkably stable structures occurring in numerous systems: for example in turbulent gases, where they are at the origin of weather phenomena [1]; in fluids with implications for biology [2]; in electromagnetic discharges [3]; and in plasmas [4]. While vortex rings have also been predicted to exist in ferromagnets [5], they have not yet been observed. Using X-ray magnetic nanotomography [6], we imaged three-dimensional structures forming closed loops in a bulk micromagnet, each composed of a vortex-antivortex pair. Based on the magnetic vorticity, a quantity analogous to hydrodynamic vorticity, we identify these configurations as magnetic vortex rings. While such structures have been predicted to exist as transient states in exchange ferromagnets [5], the vortex rings we observe exist as stable, static configurations, whose stability we attribute to the dipolar interaction. In addition, we observe stable vortex loops intersected by magnetic singularities [7], at which the magnetisation within the vortex and antivortex cores reverses. We gain insight into the stability of these states through field and thermal equilibration protocols. These measurements pave the way for the observation of complex three-dimensional solitons in bulk magnets, as well as for the development of applications based on three-dimensional magnetic structures.


title: Anonymization of labeled TOF-MRA images for brain vessel segmentation using generative adversarial networks

abstract: Anonymization and data sharing are crucial for privacy protection and acquisition of large datasets for medical image analysis. This is a big challenge, especially for neuroimaging. Here, the brain's unique structure allows for re-identification and thus requires non-conventional anonymization. Generative adversarial networks (GANs) have the potential to provide anonymous images while preserving predictive properties. Analyzing brain vessel segmentation, we trained 3 GANs on time-of-flight (TOF) magnetic resonance angiography (MRA) patches for image-label generation: 1) Deep convolutional GAN, 2) Wasserstein-GAN with gradient penalty (WGAN-GP) and 3) WGAN-GP with spectral normalization (WGAN-GP-SN). The generated image-labels from each GAN were used to train a U-net for segmentation and tested on real data. Moreover, we applied our synthetic patches using transfer learning on a second dataset. For an increasing number of up to 15 patients we evaluated the model performance on real data with and without pre-training. The performance for all models was assessed by the Dice Similarity Coefficient (DSC) and the 95th percentile of the Hausdorff Distance (95HD). Comparing the 3 GANs, the U-net trained on synthetic data generated by the WGAN-GP-SN showed the highest performance to predict vessels (DSC/95HD 0.82/28.97) benchmarked by the U-net trained on real data (0.89/26.61). The transfer learning approach showed superior performance for the same GAN compared to no pre-training, especially for one patient only (0.91/25.68 vs. 0.85/27.36). In this work, synthetic image-label pairs retained generalizable information and showed good performance for vessel segmentation. Besides, we showed that synthetic patches can be used in a transfer learning approach with independent data. This paves the way to overcome the challenges of scarce data and anonymization in medical imaging.


title: Transfers of energy through fast diffusion channels in some resonant PDEs on the circle

abstract: In this paper we consider two classes of resonant Hamiltonian PDEs on the circle with non-convex (respect to actions) first order resonant Hamiltonian. We show that, for appropriate choices of the nonlinearities we can find time-independent linear potentials that enable the construction of solutions that undergo a prescribed growth in the Sobolev norms. The solutions that we provide follow closely the orbits of a nonlinear resonant model, which is a good approximation of the full equation. The non-convexity of the resonant Hamiltonian allows the existence of fast diffusion channels along which the orbits of the resonant model experience a large drift in the actions in the optimal time. This phenomenon induces a transfer of energy among the Fourier modes of the solutions which in turn is responsible for the growth of higher order Sobolev norms.


title: Densities on Dedekind domains, completions and Haar measure

abstract: Let $D$ be the ring of $S$-integers in a global field and $\hat{D}$ its profinite completion. We discuss the relation between density in $D$ and the Haar measure of $\hat{D}$: in particular, we ask when the density of a subset $X$ of $D$ is equal to the Haar measure of its closure in $\hat{D}$. In order to have a precise statement, we give a general definition of density which encompasses the most commonly used ones. Using it we provide a necessary and sufficient condition for the equality between density and measure which subsumes a criterion due to Poonen and Stoll. In another direction, we extend the Davenport-Erd\H{o}s theorem to every $D$ as above and offer a new interpretation of it as a "density=measure" result. Our point of view also provides a simple proof that in any $D$ the set of elements divisible by at most $k$ distinct primes has density 0 for any natural number $k$. Finally, we show that the closure of the set of prime elements of $D$ is the union of the group of units of $\hat{D}$ with a negligible part.


title: The Gelfand-Kazhdan criterion as a necessary and sufficient criterion

abstract: We show that under certain conditions the Gelfand-Kazhdan criterion for the Gelfand property is a necessary condition. We work in the generality of finite groups, however part of the argument carries over to p-adic and real groups.


title: Superlubric Schottky Generator in Microscale with High Current Density and Ultralong Life

abstract: Miniaturized or even microscale generators that could effectively and persistently converse weak and random mechanical energy from environments into electricity promise huge applications in the internet of things, sensor networks, big data, personal health systems, artificial intelligence, etc. However, such generators haven't appeared yet because either the current density, or persistence, or both of all reported attempts were too low to real applications. Here, we demonstrate a superlubric Schottky generator (SLSG) in microscale such that the sliding contact between a microsized graphite flake and an n-type silicon is in a structural superlubric state, namely a ultralow friction and wearless state. This SLSG generates a stable electrical current at a high density (~119 Am-2) for at least 5,000 cycles. Since no current decay and wear were observed during the entire experiment, we believe that the real persistence of the SLSG should be enduring or substantively unlimited. In addition, the observed results exclude the mechanism of friction excitation in our Schottky generator, and provide the first experimental support of the conjectured mechanism of depletion layer establishment and destruction (DLED). Furthermore, we demonstrate a physical process of the DLED mechanism by the use of a quasi-static semiconductor finite element simulation. Our work may guide and accelerate future SLSGs into real applications.


title: Using Spectral Submanifolds for Optimal Mode Selection in Model Reduction

abstract: Model reduction of large nonlinear systems often involves the projection of the governing equations onto linear subspaces spanned by carefully-selected modes. The criteria to select the modes relevant for reduction are usually problem-specific and heuristic. In this work, we propose a rigorous mode-selection criterion based on the recent theory of Spectral Submanifolds (SSM), which facilitates a reliable projection of the governing nonlinear equations onto modal subspaces. SSMs are exact invariant manifolds in the phase space that act as nonlinear continuations of linear normal modes. Our criterion identifies critical linear normal modes whose associated SSMs have locally the largest curvature. These modes should then be included in any projection-based model reduction as they are the most sensitive to nonlinearities. To make this mode selection automatic, we develop explicit formulas for the scalar curvature of an SSM and provide an open-source numerical implementation of our mode-selection procedure. We illustrate the power of this procedure by accurately reproducing the forced-response curves on three examples of varying complexity, including high-dimensional finite element models.


title: Time dependence of quantum oscillator excitation by electromagnetic pulses

abstract: The paper is devoted to the theoretical investigation of time dependences of quantum oscillator excitation by electromagnetic pulses for arbitrary values of field amplitude in the pulse. We consider the harmonic oscillator without relaxation and excitation between stationary states. The general formula for excitation of quantum states as function of time is derived in terms of instant energy of associated classical oscillator in the field of electromagnetic pulse. Thus, it is established that time dependence of quantum oscillator excitation is completely determined by the energy of the associated classical oscillator at given moment of time. Using derived expression time dependence of quantum oscillator excitation is investigated in details including total excitation from ground state, excitation from excited states, total excitation probability and corresponding excitation spectra.


title: Optimal Trajectory Planning for Cinematography with Multiple Unmanned Aerial Vehicles

abstract: This paper presents a method for planning optimal trajectories with a team of Unmanned Aerial Vehicles (UAVs) performing autonomous cinematography. The method is able to plan trajectories online and in a distributed manner, providing coordination between the UAVs. We propose a novel non-linear formulation for this challenging problem of computing multi-UAV optimal trajectories for cinematography; integrating UAVs dynamics and collision avoidance constraints, together with cinematographic aspects like smoothness, gimbal mechanical limits and mutual camera visibility. We integrate our method within a hardware and software architecture for UAV cinematography that was previously developed within the framework of the MultiDrone project; and demonstrate its use with different types of shots filming a moving target outdoors. We provide extensive experimental results both in simulation and field experiments. We analyze the performance of the method and prove that it is able to compute online smooth trajectories, reducing jerky movements and complying with cinematography constraints.


title: Closure of the entanglement gap at quantum criticality: The case of the Quantum Spherical Model

abstract: The study of entanglement spectra is a powerful tool to detect or elucidate universal behaviour in quantum many-body systems. We investigate the scaling of the entanglement (or Schmidt) gap $\delta\xi$, i.e., the lowest laying gap of the entanglement spectrum, at a two-dimensional quantum critical point. We focus on the paradigmatic quantum spherical model, which exhibits a second-order transition, and is mappable to free bosons with an additional external constraint. We analytically show that the Schmidt gap vanishes at the critical point, although only logarithmically. For a system on a torus and the half-system bipartition, the entanglement gap vanishes as $\pi^2/\ln(L)$, with $L$ the linear system size. The entanglement gap is nonzero in the paramagnetic phase and exhibits a faster decay in the ordered phase. The rescaled gap $\delta\xi\ln(L)$ exhibits a crossing for different system sizes at the transition, although logarithmic corrections prevent a precise verification of the finite-size scaling. Interestingly, the change of the entanglement gap across the phase diagram is reflected in the zero-mode eigenvector of the spin-spin correlator. At the transition quantum fluctuations give rise to a non-trivial structure of the eigenvector, whereas in the ordered phase it is flat. We also show that the vanishing of the entanglement gap at criticality can be qualitatively but not quantitatively captured by neglecting the structure of the zero-mode eigenvector.

categories: cond-mat.stat-mech cond-mat.str-el hep-th quant-ph


title: Physical properties of new MAX phase borides M2SB (M = Zr, Hf and Nb) in comparison with conventional MAX phase carbides M2SC (M = Zr, Hf and Nb): Comprehensive insights

abstract: In this article, a detailed study of the recently synthesized MAX phase borides M2SB (M = Zr, Hf and Nb) has been performed via first principles technique. Investigation of mechanical hardness, elastic anisotropy, optical properties, dynamical stability and thermal properties are considered for the first time. The estimated values of stiffness constants and elastic moduli are found in good agreement with available results. The micro and macro hardness (Hmicro and Hmacro) parameters are calculated. The Vickers hardness is also calculated using Mulliken population analysis. The electronic density of states and charge density mapping are used to explain the variation of stiffness constants, elastic moduli and hardness parameters among the studied ternary borides. The Nb2SB compound is found to show best combination of mechanical properties. Mixture of covalent and ionic bonding within these borides is explained using Mulliken population analysis. The direction dependent values of Youngs modulus, compressibility, shear modulus and Poissons ratio are visualized by 2D and 3D representations and different anisotropic factors are calculated. The important optical constants are calculated and analyzed. The metallic nature of the studied borides is confirmed from the DOS and optical properties. The reflectivity spectra reveal the potential use of Zr2SB as coating materials to diminish solar heating. The studied borides are dynamically stable as confirmed from the phonon dispersion curves. The characteristic thermodynamic properties have also been calculated and analyzed. The physical properties of corresponding 211 MAX phase carbides are also calculated for comparison with those of the titled ternary borides.


title: Hyperspectral Image Super-Resolution via Deep Prior Regularization with Parameter Estimation

abstract: Hyperspectral image (HSI) super-resolution is commonly used to overcome the hardware limitations of existing hyperspectral imaging systems on spatial resolution. It fuses a low-resolution (LR) HSI and a high-resolution (HR) conventional image of the same scene to obtain an HR HSI. In this work, we propose a method that integrates a physical model and deep prior information. Specifically, a novel, yet effective two-stream fusion network is designed to serve as a {regularizer} for the fusion problem. This fusion problem is formulated as an optimization problem whose solution can be obtained by solving a Sylvester equation. Furthermore, the regularization parameter is simultaneously estimated to automatically adjust contribution of the physical model and {the} learned prior to reconstruct the final HR HSI. Experimental results on {both simulated and real data} demonstrate the superiority of the proposed method over other state-of-the-art methods on both quantitative and qualitative comparisons.


title: DeepSun: Machine-Learning-as-a-Service for Solar Flare Prediction

abstract: Solar flare prediction plays an important role in understanding and forecasting space weather. The main goal of the Helioseismic and Magnetic Imager (HMI), one of the instruments on NASA's Solar Dynamics Observatory, is to study the origin of solar variability and characterize the Sun's magnetic activity. HMI provides continuous full-disk observations of the solar vector magnetic field with high cadence data that lead to reliable predictive capability; yet, solar flare prediction effort utilizing these data is still limited. In this paper, we present a machine-learning-as-a-service (MLaaS) framework, called DeepSun, for predicting solar flares on the Web based on HMI's data products. Specifically, we construct training data by utilizing the physical parameters provided by the Space-weather HMI Active Region Patches (SHARP) and categorize solar flares into four classes, namely B, C, M, X, according to the X-ray flare catalogs available at the National Centers for Environmental Information (NCEI). Thus, the solar flare prediction problem at hand is essentially a multi-class (i.e., four-class) classification problem. The DeepSun system employs several machine learning algorithms to tackle this multi-class prediction problem and provides an application programming interface (API) for remote programming users. To our knowledge, DeepSun is the first MLaaS tool capable of predicting solar flares through the Internet.


title: Probabilistic Gradients for Fast Calibration of Differential Equation Models

abstract: Calibration of large-scale differential equation models to observational or experimental data is a widespread challenge throughout applied sciences and engineering. A crucial bottleneck in state-of-the art calibration methods is the calculation of local sensitivities, i.e. derivatives of the loss function with respect to the estimated parameters, which often necessitates several numerical solves of the underlying system of partial or ordinary differential equations. In this paper we present a new probabilistic approach to computing local sensitivities. The proposed method has several advantages over classical methods. Firstly, it operates within a constrained computational budget and provides a probabilistic quantification of uncertainty incurred in the sensitivities from this constraint. Secondly, information from previous sensitivity estimates can be recycled in subsequent computations, reducing the overall computational effort for iterative gradient-based calibration methods. The methodology presented is applied to two challenging test problems and compared against classical methods.

categories: math.OC cs.LG cs.NA math.NA stat.CO stat.ME stat.ML


title: Geometry-aware neural solver for fast Bayesian calibration of brain tumor models

abstract: Modeling of brain tumor dynamics has the potential to advance therapeutic planning. Current modeling approaches resort to numerical solvers that simulate the tumor progression according to a given differential equation. Using highly-efficient numerical solvers, a single forward simulation takes up to a few minutes of compute. At the same time, clinical applications of tumor modeling often imply solving an inverse problem, requiring up to tens of thousands forward model evaluations when used for a Bayesian model personalization via sampling. This results in a total inference time prohibitively expensive for clinical translation. While recent data-driven approaches become capable of emulating physics simulation, they tend to fail in generalizing over the variability of the boundary conditions imposed by the patient-specific anatomy. In this paper, we propose a learnable surrogate for simulating tumor growth which maps the biophysical model parameters directly to simulation outputs, i.e. the local tumor cell densities, whilst respecting patient geometry. We test the neural solver on Bayesian tumor model personalization for a cohort of glioma patients. Bayesian inference using the proposed surrogate yields estimates analogous to those obtained by solving the forward model with a regular numerical solver. The near-real-time computation cost renders the proposed method suitable for clinical settings. The code is available at https://github.com/IvanEz/tumor-surrogate.


title: Bidirectional teleportation for underwater quantum communications

abstract: In this work, we evaluate the performance of a bidirectional teleportation protocol on an IBM-Q's quantum processor of six or more qubits. If the experiment is successful, we will implement this protocol between two submerged nuclear submarines on opposite sides of the ocean thanks to a satellite that generates and distributes entangled pairs, as well as transmits optical bits of disambiguation between both submarines.


title: The New Compact 341 Model: Higgs Decay Modes

abstract: New developments in the anomaly free compact 341 model are discussed and the higgs bosons decay modes are studied taking into account the contributions of new fermions, gauge bosons and scalar bosons predicted by the model. It is shown from signal strengths and the branching ratios of the various decay modes analysis and the LHC constraints that there is a room for this extended BSM model and it is viable.


title: $(k+1)$-potent Matrices in triangular matrix Groups and Incidence Algebras of Finite Posets

abstract: Let $\mathbb{K}$ be a field such that $char(\mathbb{K})\nmid k$ and $char(\mathbb{K})\nmid k+1$. We describe all $(k+1)$-potent matrices over the group of upper triangular matrix. In the case that $\mathbb{K}$ is a finite field we show how to compute the number of these elements in triangular matrix groups and use this formula to compute the number of $(k+1)$-potent elements in the Incidence Algebra $\mathcal{I}(X,\mathbb{K})$ where $X$ is a finite poset.


title: A Dynamic Histories Interpretation of Quantum Theory

abstract: The problem of how to interpret quantum mechanics has persisted for a century. The disconnect between the wavefunction state vector and what is observed in experimental apparati has had no shortage of explanations. But all explanations so far fall short of a compelling and complete interpretation. In this letter, I present a novel interpretation called dynamic histories. I show mathematically how quantum mechanics can be reinterpreted as deterministically evolving dynamical world lines in a 5D universe. Quantum probabilities can be then be reinterpreted as stemming from ignorance of the state of our own world line. Meanwhile, the lack of observed superposition in experimental apparati is explained in that we only live on a single history with a definite set of properties. Hence, superposition is not an actual state of a particle but a model of ignorance as in classical probability theory. This explains nonlocal effects without superluminal communication. I also discuss how this relates to 5D Kaluza-Klein theory.


title: Unifying Two Notions of Nonlocality in Quantum Theory

abstract: Ensembles containing orthogonal product states are found to be indistinguishable under local operations and classical communication (LOCC), thereby showing irreversibility in the preparation and distinguishing processes, which is commonly known as nonlocality without entanglement. On the other hand, correlations arising from incompatible measurements on entangled states lead to Bell-nonlocality. We unify these two concepts from the change in certain property incurred in the ensemble under a suitable global unitary transformation. Specifically, we prove that under controlled-NOT (CNOT) operation, a full product basis can create entangled states if and only if the full bases or any subspace of it become irreducible in the process of LOCC discrimination. The proposed criteria quantifies the amount of nonlocality associated with the sets of product states which are even incomplete. For a set having entangled states, we modify the quantity accordingly and show that it can provide an explanation for the phenomena of more nonlocality with less entanglement.


title: Dynamics of a modified Leslie-Gower predator-prey model with Allee effect on the prey and a generalist predator

abstract: A predator-prey model with functional response Holling type II, Allee effect in the prey and a generalist predator is considered. It is shown that the model with strong Allee effect has at most two positive equilibrium point in the first quadrant, one is always a saddle point and the other exhibits multi-stability phenomenon since the equilibrium point can be stable or unstable. While the model with weak Allee effect has at most three positive equilibrium point in the first quadrant, one is always a saddle point and the other two can be stable or unstable node. In addition, when the parameters vary in a small neighbourhood of system parameters the model undergoes to different bifurcations, such as saddle-node, Hopf and Bogadonov-Takens bifurcations. Moreover, numerical simulation is used to illustrate the impact in the stability of positive equilibrium point(s) by adding an Allee effect and an alternative food source for predators.


title: Binarized Neural Architecture Search for Efficient Object Recognition

abstract: Traditional neural architecture search (NAS) has a significant impact in computer vision by automatically designing network architectures for various tasks. In this paper, binarized neural architecture search (BNAS), with a search space of binarized convolutions, is introduced to produce extremely compressed models to reduce huge computational cost on embedded devices for edge computing. The BNAS calculation is more challenging than NAS due to the learning inefficiency caused by optimization requirements and the huge architecture space, and the performance loss when handling the wild data in various computing applications. To address these issues, we introduce operation space reduction and channel sampling into BNAS to significantly reduce the cost of searching. This is accomplished through a performance-based strategy that is robust to wild data, which is further used to abandon less potential operations. Furthermore, we introduce the Upper Confidence Bound (UCB) to solve 1-bit BNAS. Two optimization methods for binarized neural networks are used to validate the effectiveness of our BNAS. Extensive experiments demonstrate that the proposed BNAS achieves a comparable performance to NAS on both CIFAR and ImageNet databases. An accuracy of $96.53\%$ vs. $97.22\%$ is achieved on the CIFAR-10 dataset, but with a significantly compressed model, and a $40\%$ faster search than the state-of-the-art PC-DARTS. On the wild face recognition task, our binarized models achieve a performance similar to their corresponding full-precision models.


title: Performance Analysis of Model-Free Adaptive Control

abstract: We analyzed model-free adaptive control (MFAC) law through closed-loop function to widen its application range.


title: Analysis of the convergence of the degree distribution of contracting random networks towards a Poisson distribution using the relative entropy

abstract: We present analytical results for the structural evolution of random networks undergoing contraction processes via generic node deletion scenarios, namely, random deletion, preferential deletion and propagating deletion. Focusing on configuration model networks, which exhibit a given degree distribution $P_0(k)$ and no correlations, we show using a rigorous argument that upon contraction the degree distributions of these networks converge towards a Poisson distribution. To this end, we use the relative entropy $S_t=S[P_t(k) || \pi(k|\langle K \rangle_t)]$ of the degree distribution $P_t(k)$ of the contracting network at time $t$ with respect to the corresponding Poisson distribution $\pi(k|\langle K \rangle_t)$ with the same mean degree $\langle K \rangle_t$ as a distance measure between $P_t(k)$ and Poisson. The relative entropy is suitable as a distance measure since it satisfies $S_t \ge 0$ for any degree distribution $P_t(k)$, while equality is obtained only for $P_t(k) = \pi(k|\langle K \rangle_t)$. We derive an equation for the time derivative $dS_t/dt$ during network contraction and show that the relative entropy decreases monotonically to zero during the contraction process. We thus conclude that the degree distributions of contracting configuration model networks converge towards a Poisson distribution. Since the contracting networks remain uncorrelated, this means that their structures converge towards an Erd{\H o}s-R\'enyi (ER) graph structure, substantiating earlier results obtained using direct integration of the master equation and computer simulations [I. Tishby, O. Biham and E. Katzav, {\it Phys. Rev. E} {\bf 100}, 032314 (2019)]. We demonstrate the convergence for configuration model networks with degenerate degree distributions (random regular graphs), exponential degree distributions and power-law degree distributions (scale-free networks).


title: Sojourn probabilities in tubes and pathwise irreversibility for It\^o processes

abstract: The sojourn probability of an It\^o diffusion process, i.e. its probability to remain in the tubular neighborhood of a smooth path, is a central quantity in the study of path probabilities. For $N$-dimensional It\^o processes with state-dependent full-rank diffusion tensor, we derive a general expression for the sojourn probability in tubes whose radii are small but finite, and fixed by the metric of the ambient Euclidean space. The central quantity in our study is the exit rate at which trajectories leave the tube for the first time. This has an interpretation as a Lagrangian and can be measured directly in experiment, unlike previously defined sojourn probabilities which depend on prior knowledge of the state-dependent diffusivity. We find that while in the limit of vanishing tube radius the ratio of sojourn probabilities for a pair of distinct paths is in general divergent, the same for a path and its time-reversal is always convergent and finite. This provides a pathwise definition of irreversibility for It\^o processes that is agnostic to the state-dependence of the diffusivity. For one-dimensional systems we derive an explicit expression for our Lagrangian in terms of the drift and diffusivity, and find that our result differs from previously reported multiplicative-noise Lagrangians. We confirm our result by comparing to numerical simulations, and relate our theory to the Stratonovich Lagrangian for multiplicative noise. For one-dimensional systems, we discuss under which conditions the vanishing-radius limiting ratio of sojourn probabilities for a pair of forward and backward paths recovers the established pathwise entropy production. Finally, we demonstrate for our one-dimensional example system that the most probable tube for a barrier crossing depends sensitively on the tube radius, and hence on the tolerated amount of fluctuations around the smooth reference path.

categories: cond-mat.stat-mech cond-mat.mes-hall cond-mat.soft physics.chem-ph


title: Probabilistic Dipole Inversion for Adaptive Quantitative Susceptibility Mapping

abstract: A learning-based posterior distribution estimation method, Probabilistic Dipole Inversion (PDI), is proposed to solve the quantitative susceptibility mapping (QSM) inverse problem in MRI with uncertainty estimation. In PDI, a deep convolutional neural network (CNN) is used to represent the multivariate Gaussian distribution as the approximate posterior distribution of susceptibility given the input measured field. Such CNN is first trained on healthy subjects via posterior density estimation, where the training dataset contains samples from the true posterior distribution. Domain adaptations are then deployed on patient datasets with new pathologies not included in pre-training, where PDI updates the pre-trained CNN's weights in an unsupervised fashion by minimizing the Kullback-Leibler divergence between the approximate posterior distribution represented by CNN and the true posterior distribution from the likelihood distribution of a known physical model and pre-defined prior distribution. Based on our experiments, PDI provides additional uncertainty estimation compared to the conventional MAP approach, while addressing the potential issue of the pre-trained CNN when test data deviates from training.


title: Analysis of Deviance for Hypothesis Testing in Generalized Partially Linear Models

abstract: In this study, we develop nonparametric analysis of deviance tools for generalized partially linear models based on local polynomial fitting. Assuming a canonical link, we propose expressions for both local and global analysis of deviance, which admit an additivity property that reduces to analysis of variance decompositions in the Gaussian case. Chi-square tests based on integrated likelihood functions are proposed to formally test whether the nonparametric term is significant. Simulation results are shown to illustrate the proposed chi-square tests and to compare them with an existing procedure based on penalized splines. The methodology is applied to German Bundesbank Federal Reserve data.


title: Structured Equilibria for Dynamic Games with Asymmetric Information and Dependent Types

abstract: We consider a dynamic game with asymmetric information where each player observes privately a noisy version of a (hidden) state of the world V, resulting in dependent private observations. We study structured perfect Bayesian equilibria that use private beliefs in their strategies as sufficient statistics for summarizing their observation history. The main difficulty in finding the appropriate sufficient statistic (state) for the structured strategies arises from the fact that players need to construct (private) beliefs on other players' private beliefs on V, which in turn would imply that an infinite hierarchy of beliefs on beliefs needs to be constructed, rendering the problem unsolvable. We show that this is not the case: each player's belief on other players' beliefs on V can be characterized by her own belief on V and some appropriately defined public belief. We then specialize this setting to the case of a Linear Quadratic Gaussian (LQG) non-zero-sum game and we characterize linear structured PBE that can be found through a backward/forward algorithm akin to dynamic programming for the standard LQG control problem. Unlike the standard LQG problem, however, some of the required quantities for the Kalman filter are observation-dependent and thus cannot be evaluated off-line through a forward recursion.


title: Cooperative Formation of Autonomous Vehicles in Mixed Traffic Flow: Beyond Platooning

abstract: Cooperative formation and control of autonomous vehicles (AVs) promise increased efficiency and safety on public roads. In mixed traffic flow consisting of AVs and human-driven vehicles (HDVs), the prevailing platooning of multiple AVs is not the only choice for cooperative formation. In this paper, we investigate how different formations of AVs impact traffic performance from a set-function optimization perspective. We first reveal a stability invariance property and a diminishing improvement property of noncooperative formation when AVs adopt typical Adaptive Cruise Control (ACC) strategies. Then, we focus on the case of cooperative formation where the AV controllers are cooperatively designed %redesign the control strategies of AVs in different formations and investigate the optimal formation of multiple AVs using set-function optimization. Two predominant optimal formations, i.e., uniform distribution and platoon formation, emerge from extensive numerical experiments. Interestingly, platooning might have the least potential to improve traffic performance when HDVs have poor string stability behavior. These results suggest more opportunities for cooperative formation of AVs, beyond platooning, in practical mixed traffic flow.


title: The homotopy type of the space of algebraic loops on a toric variety

abstract: We investigate the homotopy type of the space of tuples of polynomials inducing base-point preserving algebraic maps from the circle S1 to a toric variety X{\Sigma}. In particular, we prove a homotopy stability result for this space by combining the Vassiliev spectral sequence and the scanning map.


title: The patch like model of galaxies formation: the virial paradox, core-cusp and missing satellite problems

abstract: The patch like model of the hierarchical galaxy formation in the $\Lambda$CDM cosmological model with small damping scale is considered. In this model galaxies and clusters of galaxies are identified with rare high density peaks, what suppresses the action of random factors in the vicinity of peaks and makes the process of halos formation more rapid and regular. High concentration of irregular subhalos surrounding the central peaks and their subsequent merging just after formation allows to consider this medium as a mixture of collisionless dispersed dark matter (DM) particles and collisional subhalos. Merging of these subhalos with the central dominating halo is accompanied by tidal destruction of the central cusp, what progressively shallows the density profile and promotes formation of super massive central black holes. In the framework of this model we can reproduce the observed correlation of mass and density of virialized galaxies and clusters of galaxies known as the virial paradox [4,5]. These correlations are closely linked with the composition of DM and the shape of the power spectrum what allows to restrict them using already available observations. These correlations put constraints on the HDM and WDM models and allow to test models of cosmological inflation. We confirm that the missing satellite problem is directly linked with the virial paradox and reheating of the Universe which increases temperature and entropy of the baryons, prevents formation of first stars and divides halos into two populations: the first one includes galaxies formed before reheating which are mainly concentrated in the vicinity of the massive ones while the second population -- numerous dark halos formed after reheating -- accumulates majority of DM but does not contain stars. Their spatial distribution is more homogeneous.


title: Wolff potential estimates for supersolutions of equations with generalized Orlicz growth

abstract: In this paper, we establish pointwise estimates for supersolutions of quasilinear elliptic equations with structural conditions involving a generalized Orlicz growth in terms of a Wolff type potential. As a consequence, under the extra assumption, we obtain that the supersolutions satisfy a Harnack inequality and local H\"older continuity.


title: Bandit Learning in Convex Non-Strictly Monotone Games

abstract: We address learning Nash equilibria in convex games under the payoff information setting. We consider the case in which the game pseudo-gradient is monotone but not necessarily strictly monotone. This relaxation of strict monotonicity enables application of learning algorithms to a larger class of games, such as, for example, a zero-sum game with a merely convex-concave cost function. We derive an algorithm whose iterates provably converge to the least-norm Nash equilibrium in this setting. {From the perspective of a single player using the proposed algorithm, we view the game as an instance of online optimization}. Through this lens, we quantify the regret rate of the algorithm and provide an approach to choose the algorithm's parameters to minimize the regret rate.


title: Completeness in Polylogarithmic Time and Space

abstract: Complexity theory can be viewed as the study of the relationship between computation and applications, understood the former as complexity classes and the latter as problems. Completeness results are clearly central to that view. Many natural algorithms resulting from current applications have polylogarithmic time (PolylogTime) or space complexity (PolylogSpace). The classical Karp notion of complete problem however does not plays well with these complexity classes. It is well known that PolylogSpace does not have complete problems under logarithmic space many-one reductions. In this paper we show similar results for deterministic and non-deterministic PolylogTime as well as for every other level of the polylogarithmic time hierarchy. We achieve that by following a different strategy based on proving the existence of proper hierarchies of problems inside each class. We then develop an alternative notion of completeness inspired by the concept of uniformity from circuit complexity and prove the existence of a (uniformly) complete problem for PolylogSpace under this new notion. As a consequence of this result we get that complete problems can still play an important role in the study of the interrelationship between polylogarithmic and other classical complexity classes.


title: Long-time asymptotics and stability for the sine-Gordon equation

abstract: In this paper, we study the long-time dynamics and stability properties of the sine-Gordon equation $$f_{tt}-f_{xx}+\sin f=0.$$ Firstly, we use the nonlinear steepest descent for Riemann-Hilbert problems to compute the long-time asymptotics of the solutions to the sine-Gordon equation whose initial condition belongs to some weighted Sobolev spaces. Secondly, we study the asymptotic stability of the sine-Gordon equation. It is known that the obstruction to the asymptotic stability of the sine-Gordon equation in the energy space is the existence of small breathers which is also closely related to the emergence of wobbling kinks. Combining the long-time asymptotics and a refined approximation argument, we analyze the asymptotic stability properties of the sine-Gordon equation in weighted energy spaces. Our stability analysis gives a criterion for the weight which is sharp up to the endpoint so that the asymptotic stability holds.


title: Study of the Galactic Rotation Based on Masers and Radio Stars with VLBI Measurements of Their Parallaxes

abstract: Based on published data, we have produced a sample of 256 radio sources whose trigonometric parallaxes and proper motions were measured by VLBI. This sample contains Galactic masers associated with massive protostars and stars in active star-forming regions. It also includes young low-mass stars from the Gould Belt region whose radio observations were performed in continuum. Based on this, most complete sample of sources to date, we have estimated the velocities $(U,V,W)_\odot$ and the parameters of the angular velocity of Galactic rotation $\Omega_0, \Omega^{(1)}_0, ... , \Omega^{(4)}_0$ and obtained a new estimate of the distance fromthe Sun to the Galactic center, $R_0=8.15^{+0.04}_{-0.20}$ kpc. The parameters of the Galactic spiral density wave have been found from the series of radial, $V_R$, and residual tangential, $\Delta V_{circ}$, velocities of stars. The amplitudes of the radial and tangential velocity perturbations are $f_R=7.0\pm0.9$ km s$^{-1}$ and $f_\theta=3.8\pm1.1$ km s$^{-1}$, the perturbation wavelengths are $\lambda_R=2.3\pm0.2$ kpc and $\lambda_\theta=2.0\pm0.4$ kpc, and the Sun's phases in the spiral density wave are $(\chi_\odot)_R=-163^\circ\pm9^\circ$ and $(\chi_\odot)_\theta=-137^\circ\pm10^\circ$ for the adopted four-armed spiral pattern.


title: Towards optimization techniques on diffeological spaces by generalizing Riemannian concepts

abstract: Diffeological spaces firstly introduced by J.M. Souriau in the 1980s are a natural generalization of smooth manifolds. However, optimization techniques are only known on manifolds so far. Generalizing these techniques to diffeological spaces is very challenging because of several reasons. One of the main reasons is that there are various definitions of tangent spaces which do not coincide. Additionally, one needs to deal with a generalization of a Riemannian space in order to define gradients which are indispensable for optimization methods. This paper is devoted to an optimization technique on diffeological spaces. Thus, one main aim of this paper is a suitable definition of a tangent space in view to optimization methods. Based on this definition, we present a diffeological Riemannian space and a diffeological gradient, which we need to formulate an optimization algorithm on diffeological spaces. Moreover, in order to be able to update the iterates in an optimization algorithm on diffeological spaces, we present a diffeological retraction and the Levi-Civita connection on diffeological spaces. We give examples for the novel objects and apply the presented diffeological algorithm to an optimization problem.


title: Real-World Snapshots vs. Theory: Questioning the t-Probing Security Model

abstract: Due to its sound theoretical basis and practical efficiency, masking has become the most prominent countermeasure to protect cryptographic implementations against physical side-channel attacks (SCAs). The core idea of masking is to randomly split every sensitive intermediate variable during computation into at least t+1 shares, where t denotes the maximum number of shares that are allowed to be observed by an adversary without learning any sensitive information. In other words, it is assumed that the adversary is bounded either by the possessed number of probes (e.g., microprobe needles) or by the order of statistical analyses while conducting higher-order SCA attacks (e.g., differential power analysis). Such bounded models are employed to prove the SCA security of the corresponding implementations. Consequently, it is believed that given a sufficiently large number of shares, the vast majority of known SCA attacks are mitigated. In this work, we present a novel laser-assisted SCA technique, called Laser Logic State Imaging (LLSI), which offers an unlimited number of contactless probes, and therefore, violates the probing security model assumption. This technique enables us to take snapshots of hardware implementations, i.e., extract the logical state of all registers at any arbitrary clock cycle with a single measurement. To validate this, we mount our attack on masked AES hardware implementations and practically demonstrate the extraction of the full-length key in two different scenarios. First, we assume that the location of the registers (key and/or state) is known, and hence, their content can be directly read by a single snapshot. Second, we consider an implementation with unknown register locations, where we make use of multiple snapshots and a SAT solver to reveal the secrets.


title: Unsupervised Part Discovery by Unsupervised Disentanglement

abstract: We address the problem of discovering part segmentations of articulated objects without supervision. In contrast to keypoints, part segmentations provide information about part localizations on the level of individual pixels. Capturing both locations and semantics, they are an attractive target for supervised learning approaches. However, large annotation costs limit the scalability of supervised algorithms to other object categories than humans. Unsupervised approaches potentially allow to use much more data at a lower cost. Most existing unsupervised approaches focus on learning abstract representations to be refined with supervision into the final representation. Our approach leverages a generative model consisting of two disentangled representations for an object's shape and appearance and a latent variable for the part segmentation. From a single image, the trained model infers a semantic part segmentation map. In experiments, we compare our approach to previous state-of-the-art approaches and observe significant gains in segmentation accuracy and shape consistency. Our work demonstrates the feasibility to discover semantic part segmentations without supervision.


title: The Dynamics of the Wide-Angle Tailed (WAT) galaxy cluster Abell 562

abstract: We present the first in depth dynamical analysis of the archetypal Wide-Angle Tailed (WAT) cluster Abell 562. We have combined Gemini observations with archival data from the literature to form a sample of 76 cluster members and derived a mean redshift of $0.1088 \pm 0.0004$ and a velocity dispersion of $919 \pm 116$ km s$^{-1}$. This relatively large velocity dispersion suggests either a very massive cluster ($M_{dyn} > 6.9 \times$ 10$^{14}$ $M_{sun}$) and/or a merger system. The merger model is supported by a non-Gaussian galaxy velocity distribution, an elongated spatial distribution of likely cluster members, and an elongated X-ray emitting gas. This scenario would generate the bulk flow motion of the intra-cluster medium that can exert enough ram pressure to bend the radio jets. Thus, our observations support the model in which a recent off-axis merger event produced the cluster wide conditions needed to shape the WAT in Abell 562.


title: The Unbalanced Gromov Wasserstein Distance: Conic Formulation and Relaxation

abstract: Comparing metric measure spaces (i.e. a metric space endowed with aprobability distribution) is at the heart of many machine learning problems. The most popular distance between such metric measure spaces is theGromov-Wasserstein (GW) distance, which is the solution of a quadratic assignment problem. The GW distance is however limited to the comparison of metric measure spaces endowed with a probability distribution. To alleviate this issue, we introduce two Unbalanced Gromov-Wasserstein formulations: a distance and a more tractable upper-bounding relaxation.They both allow the comparison of metric spaces equipped with arbitrary positive measures up to isometries. The first formulation is a positive and definite divergence based on a relaxation of the mass conservation constraint using a novel type of quadratically-homogeneous divergence. This divergence works hand in hand with the entropic regularization approach which is popular to solve large scale optimal transport problems. We show that the underlying non-convex optimization problem can be efficiently tackled using a highly parallelizable and GPU-friendly iterative scheme. The second formulation is a distance between mm-spaces up to isometries based on a conic lifting. Lastly, we provide numerical experiments onsynthetic examples and domain adaptation data with a Positive-Unlabeled learning task to highlight the salient features of the unbalanced divergence and its potential applications in ML.