title: Singular anisotropic elliptic equations with gradient-dependent lower order terms

abstract: We prove the existence of a solution to a singular anisotropic elliptic equation in a bounded open subset $\Omega$ of $\mathbb R^N$ with $N\ge 2$, subject to a homogeneous boundary condition: \begin{equation} \label{eq0} \left\{ \begin{array}{ll} \mathcal A u+ \Phi(u,\nabla u)=\Psi(u,\nabla u)+ \mathfrak{B} u \quad& \mbox{in } \Omega,\\ u=0 & \mbox{on } \partial\Omega. \end{array} \right. \end{equation} Here $ \mathcal A u=-\sum_{j=1}^N |\partial_j u|^{p_j-2}\partial_j u$ is the anisotropic $\overrightarrow{p}$-Laplace operator, while $\mathfrak B$ is an operator from $W_0^{1,\overrightarrow{p}}(\Omega)$ into $W^{-1,\overrightarrow{p}'}(\Omega)$ satisfying suitable, but general, structural assumptions. $\Phi$ and $\Psi$ are gradient-dependent nonlinearities whose models are the following: \begin{equation*} \label{phi}\Phi(u,\nabla u):=\left(\sum_{j=1}^N \mathfrak{a}_j |\partial_j u|^{p_j}+1\right)|u|^{m-2}u, \quad \Psi(u,\nabla u):=\frac{1}{u}\sum_{j=1}^N |u|^{\theta_j} |\partial_j u|^{q_j}. \end{equation*} We suppose throughout that, for every $1\leq j\leq N$, \begin{equation*}\label{ass} \mathfrak{a}_j\geq 0, \quad \theta_j>0, \quad 0\leq q_j<p_j, \quad 1<p_j,m\quad \mbox{and}\quad p<N, \end{equation*} and we distinguish two cases: 1) for every $1\leq j\leq N$, we have $\theta_j\geq 1$; 2) there exists $1\leq j\leq N$ such that $\theta_j<1$. In this last situation, we look for non-negative solutions of \eqref{eq0}.


title: QCD phase diagram at finite isospin chemical potential and temperature in an IR-improved soft-wall AdS/QCD model

abstract: We study the phase transition between pion condensed phase and normal phase, as well as chiral phase transition in a two flavor($\mathcal{N}_f=2$) IR- improved soft-wall AdS/QCD model at finite isospin chemical potential $\mu_I$ and temperature $T$. By self-consistently solving the equations of motion, we obtain the phase diagram in the plane of $\mu_I$ and $T$. The pion condensation appears together with a massless Nambu-Goldstone boson $m_{\pi_1}(T_c, \mu_I^c)=0$, which is very likely to be a second-order phase transition with mean-field critical exponents in small $\mu_I$ region. When $T=0$, the critical isospin chemical potential approximates to vacuum pion mass $\mu_I^c \approx m_0$. The pion condensed phase exists in an arched area and the boundary of the chiral crossover intersects the pion condensed phase at a tri-critical point. Qualitatively, the results are in good agreement with previous studies from Lattice simulations and model calculations.


title: Probabilistic Reasoning across the Causal Hierarchy

abstract: We propose a formalization of the three-tier causal hierarchy of association, intervention, and counterfactuals as a series of probabilistic logical languages. Our languages are of strictly increasing expressivity, the first capable of expressing quantitative probabilistic reasoning -- including conditional independence and Bayesian inference -- the second encoding do-calculus reasoning for causal effects, and the third capturing a fully expressive do-calculus for arbitrary counterfactual queries. We give a corresponding series of finitary axiomatizations complete over both structural causal models and probabilistic programs, and show that satisfiability and validity for each language are decidable in polynomial space.


title: Deep Plastic Surgery: Robust and Controllable Image Editing with Human-Drawn Sketches

abstract: Sketch-based image editing aims to synthesize and modify photos based on the structural information provided by the human-drawn sketches. Since sketches are difficult to collect, previous methods mainly use edge maps instead of sketches to train models (referred to as edge-based models). However, sketches display great structural discrepancy with edge maps, thus failing edge-based models. Moreover, sketches often demonstrate huge variety among different users, demanding even higher generalizability and robustness for the editing model to work. In this paper, we propose Deep Plastic Surgery, a novel, robust and controllable image editing framework that allows users to interactively edit images using hand-drawn sketch inputs. We present a sketch refinement strategy, as inspired by the coarse-to-fine drawing process of the artists, which we show can help our model well adapt to casual and varied sketches without the need for real sketch training data. Our model further provides a refinement level control parameter that enables users to flexibly define how "reliable" the input sketch should be considered for the final output, balancing between sketch faithfulness and output verisimilitude (as the two goals might contradict if the input sketch is drawn poorly). To achieve the multi-level refinement, we introduce a style-based module for level conditioning, which allows adaptive feature representations for different levels in a singe network. Extensive experimental results demonstrate the superiority of our approach in improving the visual quality and user controllablity of image editing over the state-of-the-art methods.


title: Theoretical Analysis of Double Differential Cross Section of Proton, Deuteron and Triton for $p+^7$Li Reaction at 14 MeV

abstract: Based on the statistical theory of light nucleus reactions (STLN), the description of the complicated emission processes of proton and light composite charged particles are further improved through considering the effects of Coulomb barriers both in incident and different outgoing reaction channels. And the analysis of the reaction channels including the sequential and simultaneous emission processes for $p + ^7$Li reaction is performed in detail. So the partial spectra of all of outgoing particles are also obtained for different reaction processes. The calculated double differential cross sections of total outgoing proton, deuteron and triton at $E_p = 14$ MeV agree well with the available experimental data for different outgoing angles. The ENDF-6 formatted data, which includes all of the reaction cross sections, elastic angular distributions, double differential cross sections of nucleon and light composite charged particles for $p + ^7$Li reaction, are also obtained by PUNF code.


title: A Generalized Probabilistic Learning Approach for Multi-Fidelity Uncertainty Propagation in Complex Physical Simulations

abstract: Two of the most significant challenges in uncertainty quantification pertain to the high computational cost for simulating complex physical models and the high dimension of the random inputs. In applications of practical interest, both of these problems are encountered, and standard methods either fail or are not feasible. To overcome the current limitations, we present a generalized formulation of a Bayesian multi-fidelity Monte-Carlo (BMFMC) framework that can exploit lower-fidelity model versions in a small data regime. The goal of our analysis is an efficient and accurate estimation of the complete probabilistic response for high-fidelity models. BMFMC circumvents the curse of dimensionality by learning the relationship between the outputs of a reference high-fidelity model and potentially several lower-fidelity models. While the continuous formulation is mathematically exact and independent of the low-fidelity model's accuracy, we address challenges associated with the small data regime (i.e., only a small number of 50 to 300 high-fidelity model runs can be performed). Specifically, we complement the formulation with a set of informative input features at no extra cost. Despite the inaccurate and noisy information that some low-fidelity models provide, we demonstrate that accurate and certifiable estimates for the quantities of interest can be obtained for uncertainty quantification problems in high stochastic dimensions, with significantly fewer high-fidelity model runs than state-of-the-art methods for uncertainty quantification. We illustrate our approach by applying it to challenging numerical examples such as Navier-Stokes flow simulations and fluid-structure interaction problems.


title: Whitney Approximation for Smooth CW Complex

abstract: We show a Whitney Approximation Theorem for a continuous map from a manifold to a smooth CW complex. This enables us to show that a topological CW complex is homotopy equivalent to a smooth CW complex in a category of topological spaces. It is also shown that, for any open covering of a smooth CW complex, there exists a partition of unity subordinate to the open covering.


title: Supervised Hyperalignment for multi-subject fMRI data alignment

abstract: Hyperalignment has been widely employed in Multivariate Pattern (MVP) analysis to discover the cognitive states in the human brains based on multi-subject functional Magnetic Resonance Imaging (fMRI) datasets. Most of the existing HA methods utilized unsupervised approaches, where they only maximized the correlation between the voxels with the same position in the time series. However, these unsupervised solutions may not be optimum for handling the functional alignment in the supervised MVP problems. This paper proposes a Supervised Hyperalignment (SHA) method to ensure better functional alignment for MVP analysis, where the proposed method provides a supervised shared space that can maximize the correlation among the stimuli belonging to the same category and minimize the correlation between distinct categories of stimuli. Further, SHA employs a generalized optimization solution, which generates the shared space and calculates the mapped features in a single iteration, hence with optimum time and space complexities for large datasets. Experiments on multi-subject datasets demonstrate that SHA method achieves up to 19% better performance for multi-class problems over the state-of-the-art HA algorithms.


title: Hybridization between the ligand $p$ band and Fe-3$d$ orbitals in the p-type ferromagnetic semiconductor (Ga,Fe)Sb

abstract: (Ga,Fe)Sb is a promising ferromagnetic semiconductor for practical spintronic device applications because its Curie temperature ($T_{\rm C}$) is above room temperature. However, the origin of ferromagnetism with high $T_{\rm C}$ remains to be elucidated. Here, we use soft x-ray angle-resolved photoemission spectroscopy (SX-ARPES) to investigate the valence-band (VB) structure of (Ga$_{0.95}$,Fe$_{0.05}$)Sb including the Fe-3$d$ impurity band (IB), to unveil the mechanism of ferromagnetism in (Ga,Fe)Sb. We find that the VB dispersion in (Ga$_{0.95}$,Fe$_{0.05}$)Sb observed by SX-ARPES is similar to that of GaSb, indicating that the doped Fe atoms hardly affect the band dispersion. The Fe-3$d$ resonant ARPES spectra demonstrate that the Fe-3$d$ IB crosses the Fermi level ($E_{\rm F}$) and hybridizes with the VB of GaSb. These observations indicate that the VB structure of (Ga$_{0.95}$,Fe$_{0.05}$)Sb is consistent with that of the IB model which is based on double-exchange interaction between the localized 3$d$ electrons of the magnetic impurities. The results indicate that the ferromagnetism in (Ga,Fe)Sb is formed by the hybridization of the Fe-3$d$ IB with the ligand $p$ band of GaSb.


title: Probing the concept of line tension down to the nanoscale

abstract: A novel mechanical approach is developed to explore by means of atom-scale simulation the concept of line tension at a solid-liquid-vapor contact line as well as its dependence on temperature, confinement, and solid/fluid interactions. More precisely, by estimating the stresses exerted along and normal to a straight contact line formed within a partially wet pore, the line tension can be estimated while avoiding the pitfalls inherent to the geometrical scaling methodology based on hemispherical drops. The line tension for Lennard-Jones fluids is found to follow a generic behavior with temperature and chemical potential effects that are all included in a simple contact angle parameterization. Former discrepancies between theoretical modeling and molecular simulation are resolved, and the line tension concept is shown to be robust down to molecular confinements. The same qualitative behavior is observed for water but the line tension at the wetting transition diverges or converges towards a finite value depending on the range of the solid/fluid interactions at play.


title: Second law of black hole thermodynamics

abstract: If simple entropy in the Bekenstein-Hawking area law for a Schwarzschild black hole is replaced with 'negative' quantum conditional entropy, which quantifies quantum entanglement, of positive-energy particles of the black hole relative to its outside, a paradox with the original pair-creation picture of Hawking radiation, the first law for black hole mechanics and quantum mechanics is resolved. However, there was no way to judge experimentally which area law is indeed adopted by black holes. Here, with the no-hair conjecture, we derive the perfect picture of a second law of black hole thermodynamics for any black hole from the modified area law, rather than Bekenstein's generalized one from the original area law. The second law is testable with an event horizon telescope, in contrast to Bekenstein's. If this is confirmed, the modified area law could be exalted to the first example of fundamental equations in physics which cannot be described without the concept of quantum information.

categories: hep-th astro-ph.HE cond-mat.stat-mech gr-qc quant-ph


title: Drop Dynamics of Viscoelastic Filament

abstract: The stretching of viscoelastic polymer solutions close to break-up can create attached drops on a filament, whose properties and dynamics are little understood. The stretching of capillary bridges and the consecutive filament, until its breakup, can be quantified using diameter-space-time diagrams, which demonstrate hierarchy, as well as, asymmetry of satellite drops around a big central drop. All drops experience migration, oscillation and merging. In addition, the position of the minimum diameter on the filament is determined, along with the number of drops, their positions, the diameters of drops and the filament breakup time. The maximum number of drops on the filament can be predicted using the Deborah number. The diagrams also quantify the large Hencky strains in the filaments before pinch-off. The obtained minimum diameter is used to measure the extensional viscosity, which indicates the effect of polymer concentration and direction of filament thinning.


title: Self-Supervised Fast Adaptation for Denoising via Meta-Learning

abstract: Under certain statistical assumptions of noise, recent self-supervised approaches for denoising have been introduced to learn network parameters without true clean images, and these methods can restore an image by exploiting information available from the given input (i.e., internal statistics) at test time. However, self-supervised methods are not yet combined with conventional supervised denoising methods which train the denoising networks with a large number of external training samples. Thus, we propose a new denoising approach that can greatly outperform the state-of-the-art supervised denoising methods by adapting their network parameters to the given input through selfsupervision without changing the networks architectures. Moreover, we propose a meta-learning algorithm to enable quick adaptation of parameters to the specific input at test time. We demonstrate that the proposed method can be easily employed with state-of-the-art denoising networks without additional parameters, and achieve state-of-the-art performance on numerous benchmark datasets.


title: 60 years of attosecond physics at ICPEAC: from collisions to ultrashort pulses

abstract: The field of attosecond physics has seen an almost explosive growth since the early 2000's and represents by now an increasing fraction of contributions to the bi-annual series of International Conferences of Photonic, Electronic, and Atomic Collisions (ICPEAC). The latter is anything but a coincidence as many of the underlying concepts of electronic and photonic dynamics are closely intertwined with atomic-scale collision processes. We illustrate this fruitful connection and its implications with the help of a few prototypical examples of current topical interest.


title: Phase resolved joint spectra tomography of a ring resonator photon pair source using a silicon photonic chip

abstract: The exponential growth of photonic quantum technologies is driving the demand of tools for measuring the quality of their information carriers. One of the most prominent is Stimulated Emission Tomography (SET), which uses classical coherent fields to measure the Joint Spectral Amplitude (JSA) of photon pairs with high speed and resolution. While the modulus of the JSA can be directly addressed from a single intensity measurement, the retrieval of the Joint Spectral Phase (JSP) is far more challenging and received minor attentions. However, a wide class of spontaneous sources of technological relevance, as chip integrated micro-resonators, have a JSP with a rich structure, that carries correlations hidden in the intensity domain. Here, using a compact and reconfigurable silicon photonic chip, it is measured for the first time the complex JSA of a micro-ring resonator photon pair source. The photonic circuit coherently excites the ring and a reference waveguide, and the interferogram formed by their stimulated fields is used to map the ring JSP through a novel phase reconstruction technique. This tool complements the traditionally bulky and sophisticated methods implemented so far, simultaneously minimizing the set of required resources.


title: A spectral geometry for the Standard Model without the fermion doubling

abstract: We propose a simple model of noncommutative geometry to describe the structure of the Standard Model, which satisfies spin${}_c$ condition, has no fermion doubling, does not lead to the possibility of color symmetry breaking and explains the CP-violation as the failure of the reality condition for the Dirac operator.


abstract: Thermodynamic properties of graphene bilayers are studied by path-integral molecular dynamics (PIMD) simulations, considering quantization of vibrational modes and anharmonic effects. Bilayer graphene has been studied at temperatures between 12 and 1500~K for zero external stress, using the LCBOPII effective potential. We concentrate on the thermal expansion, in-plane and out-of-plane compressibility, and specific heat. Additional insight into the meaning of our results for bilayer graphene is obtained from a comparison with data obtained from PIMD simulations for monolayer graphene and graphite. They are also analyzed in view of experimental data for graphite. Zero-point and thermal effects on the in-plane and "real" area of bilayer graphene are studied. The thermal expansion coefficient $\alpha_{xy}$ of the in-plane area is negative at low temperatures and positive for $T \gtrsim$ 800~K. The minimum $\alpha_{xy}$ is $-6.6 \times 10^{-6}$ K$^{-1}$ at $T \approx 220$~K. Both in-plane ($\chi_{xy}$) and out-of-plane ($\chi_z$) compressibilities of graphene bilayers are found to increase for rising temperature, and turn out to be lower than that corresponding to monolayer graphene and higher than those found for graphite. At 300 K, we find for the bilayer $\chi_{xy} = 9.5 \times 10^{-2}$ \AA$^2$/eV and $\chi_z = 2.97 \times 10^{-2}$ GPa$^{-1}$. Results for the specific heat obtained from the simulations are compared with those given by a harmonic approximation for the vibrational modes. This approach is noticeably accurate at temperatures lower than 200~K.


title: A Hybrid Volume-Surface-Wire Integral Equation for the Anisotropic Forward Problem in Electroencephalography

abstract: Solving the electroencephalography (EEG) forward problem is a fundamental step in a wide range of applications including biomedical imaging techniques based on inverse source localization. State-of-the-art electromagnetic solvers resort to a computationally expensive volumetric discretization of the full head to account for its complex and heterogeneous electric profile. The more efficient, popular in biomedical imaging circles, but unfortunately oversimplifying Boundary Element Method (BEM) relies instead on a piecewise-uniform approximation that severely curbs its application in high resolution EEGs. This contribution lifts the standard BEM contraints by treating the local anisotropies with adequate wire and thin volume integral equations that are tailored to specific structures of the fibrous white matter and the inhomogeneous skull. The proposed hybrid integral equation formulation thereby avoids the full volumetric discretization of the head medium and allows for a realistic and efficient BEM-like solution of the anisotropic EEG forward problem. The accuracy and flexibility of the proposed formulation is demonstrated through numerical experiments involving both canonical and realistic MRI-based head models.


title: Fast Adaptation to Super-Resolution Networks via Meta-Learning

abstract: Conventional supervised super-resolution (SR) approaches are trained with massive external SR datasets but fail to exploit desirable properties of the given test image. On the other hand, self-supervised SR approaches utilize the internal information within a test image but suffer from computational complexity in run-time. In this work, we observe the opportunity for further improvement of the performance of SISR without changing the architecture of conventional SR networks by practically exploiting additional information given from the input image. In the training stage, we train the network via meta-learning; thus, the network can quickly adapt to any input image at test time. Then, in the test stage, parameters of this meta-learned network are rapidly fine-tuned with only a few iterations by only using the given low-resolution image. The adaptation at the test time takes full advantage of patch-recurrence property observed in natural images. Our method effectively handles unknown SR kernels and can be applied to any existing model. We demonstrate that the proposed model-agnostic approach consistently improves the performance of conventional SR networks on various benchmark SR datasets.


title: Multirobot Coverage of Linear Modular Environments

abstract: Multirobot systems for covering environments are increasingly used in applications like cleaning, industrial inspection, patrolling, and precision agriculture. The problem of covering a given environment using multiple robots can be naturally formulated and studied as a multi-Traveling Salesperson Problem (mTSP). In a mTSP, the environment is represented as a graph and the goal is to find tours (starting and ending at the same depot) for the robots in order to visit all the vertices with minimum global cost, namely the length of the longest tour. The mTSP is an NP-hard problem for which several approximation algorithms have been proposed. These algorithms usually assume generic environments, but tighter approximation bounds can be reached focusing on specific environments. In this paper, we address the case of environments composed of sub-parts, called modules, that can be reached from each other only through some linking structures. Examples are multi-floor buildings, in which the modules are the floors and the linking structures are the staircases or the elevators, and floors of large hotels or hospitals, in which the modules are the rooms and the linking structures are the corridors. We focus on linear modular environments, with the modules organized sequentially, presenting an efficient (with polynomial worst-case time complexity) algorithm that finds a solution for the mTSP whose cost is within a bounded distance from the cost of the optimal solution. The main idea of our algorithm is to allocate disjoint "blocks" of adjacent modules to the robots, in such a way that each module is covered by only one robot. We experimentally compare our algorithm against some state-of-the-art algorithms for solving mTSPs in generic environments and show that it is able to provide solutions with lower makespan and spending a computing time several orders of magnitude shorter.


title: Population-Guided Parallel Policy Search for Reinforcement Learning

abstract: In this paper, a new population-guided parallel learning scheme is proposed to enhance the performance of off-policy reinforcement learning (RL). In the proposed scheme, multiple identical learners with their own value-functions and policies share a common experience replay buffer, and search a good policy in collaboration with the guidance of the best policy information. The key point is that the information of the best policy is fused in a soft manner by constructing an augmented loss function for policy update to enlarge the overall search region by the multiple learners. The guidance by the previous best policy and the enlarged range enable faster and better policy search. Monotone improvement of the expected cumulative return by the proposed scheme is proved theoretically. Working algorithms are constructed by applying the proposed scheme to the twin delayed deep deterministic (TD3) policy gradient algorithm. Numerical results show that the constructed algorithm outperforms most of the current state-of-the-art RL algorithms, and the gain is significant in the case of sparse reward environment.


title: Spatial-Temporal Transformer Networks for Traffic Flow Forecasting

abstract: Traffic forecasting has emerged as a core component of intelligent transportation systems. However, timely accurate traffic forecasting, especially long-term forecasting, still remains an open challenge due to the highly nonlinear and dynamic spatial-temporal dependencies of traffic flows. In this paper, we propose a novel paradigm of Spatial-Temporal Transformer Networks (STTNs) that leverages dynamical directed spatial dependencies and long-range temporal dependencies to improve the accuracy of long-term traffic forecasting. Specifically, we present a new variant of graph neural networks, named spatial transformer, by dynamically modeling directed spatial dependencies with self-attention mechanism to capture realtime traffic conditions as well as the directionality of traffic flows. Furthermore, different spatial dependency patterns can be jointly modeled with multi-heads attention mechanism to consider diverse relationships related to different factors (e.g. similarity, connectivity and covariance). On the other hand, the temporal transformer is utilized to model long-range bidirectional temporal dependencies across multiple time steps. Finally, they are composed as a block to jointly model the spatial-temporal dependencies for accurate traffic prediction. Compared to existing works, the proposed model enables fast and scalable training over a long range spatial-temporal dependencies. Experiment results demonstrate that the proposed model achieves competitive results compared with the state-of-the-arts, especially forecasting long-term traffic flows on real-world PeMS-Bay and PeMSD7(M) datasets.


title: On Optimal Locally Repairable Codes and Generalized Sector-Disk Codes

abstract: Optimal locally repairable codes with information locality are considered. Optimal codes are constructed, whose length is also order-optimal with respect to a new bound on the code length derived in this paper. The length of the constructed codes is super-linear in the alphabet size, which improves upon the well known pyramid codes, whose length is only linear in the alphabet size. The recoverable erasure patterns are also analyzed for the new codes. Based on the recoverable erasure patterns, we construct generalized sector-disk (GSD) codes, which can recover from disk erasures mixed with sector erasures in a more general setting than known sector-disk (SD) codes. Additionally, the number of sectors in the constructed GSD codes is super-linear in the alphabet size, compared with known SD codes, whose number of sectors is only linear in the alphabet size.


title: Bismuth iron garnet: ab initio study of electronic properties

abstract: Bismuth iron garnet (BIG), i.e. Bi3Fe5O12, is a strong ferrimagnet that also possess outstanding magneto-optical properties such as the largest known Faraday rotation. These properties are related with the distribution of magnetic moments on octahedral and tetrahedral sites, the presence of spin gaps in the density of state and a strong spin-orbit coupling. In this work, first-principles ab initio calculations are performed to study the structural, electronic and magnetic properties of BIG using Density Functional Theory with Hubbard+U (DFT+U) correction including spin-orbit coupling and HSE06 hybrid functional. We found that the presence of spin gaps in the electronic structure results from the interplay between exchange and correlation effects and the crystal field strengths for tetrahedral and octahedral iron sublattices. The DFT+U treatment tends to close the spin-gaps for larger U due to over-localization effects, notably in the octahedral site. On the other hand, the hybrid functional confirms the occurrences of three spin gaps in the iron states of the conduction band as expected from optical measurements. A strong exchange splitting at the top of the valence bands associated with a lone-pair type mixture of O p and Bi s,p states is also obtained. Similar exchange splitting was not previously observed for other iron based garnets, such as for yttrium iron garnet. It follows that hole doping, as obtained by Ca substitution at Bi sites, results in a full spin polarized density at the Fermi energy. This work helps to shed more light on the theoretical comprehension of the properties of BIG and opens the route towards the use of advanced Many Body calculations to predict the magneto-optical coupling effects in BIG in a direct comparison with the experimental measurements.


title: Cluster-based network model

abstract: We propose an automatable data-driven methodology for robust nonlinear reduced-order modelling from time-resolved snapshot data. In the kinematical coarse-graining, the snapshots are clustered into few centroids representable for the whole ensemble. The dynamics is conceptualized as a directed network, where the centroids represent nodes and the directed edges denote possible finite-time transitions. The transition probabilities and times are inferred from the snapshot data. The resulting cluster-based network model constitutes a deterministic-stochastic grey-box model resolving the coherent-structure evolution. This model is motivated by limit-cycle dynamics, illustrated for the chaotic Lorenz attractor and successfully demonstrated for the laminar two-dimensional mixing layer featuring Kelvin-Helmholtz vortices and vortex pairing, and for an actuated turbulent boundary layer with complex dynamics. Cluster-based network modelling opens a promising new avenue with unique advantages over other model-order reductions based on clustering or proper orthogonal decomposition.


title: Conversational Search for Learning Technologies

abstract: Conversational search is based on a user-system cooperation with the objective to solve an information-seeking task. In this report, we discuss the implication of such cooperation with the learning perspective from both user and system side. We also focus on the stimulation of learning through a key component of conversational search, namely the multimodality of communication way, and discuss the implication in terms of information retrieval. We end with a research road map describing promising research directions and perspectives.


title: Inverse flexoelectret effect: bending dielectrics by a uniform electric field

abstract: It is highly desirable to discover an electromechanical coupling that allows a dielectric material to generate curvature in response to a uniform electric field, which would add a new degree of freedom for designing actuators. Flexoelectricity, a two-way coupling between polarization and strain gradient, is a good candidate. But its applications are usually limited to the nanoscale due to its inherent size dependence. Here, an inverse flexoelectret effect in silicone elastomers is introduced to overcome this limitation. Based on this idea, a flexing actuator which can generate large curvature at the millimeter length scale is fabricated and shown to have excellent actuation performance comparable with current nanoscale flexoelectric actuators. Theoretical analysis indicates that the new phenomenon originates from the interplay of electrets and Maxwell stress. This work opens an avenue for applying macroscopic flexoelectricity in actuators and flexible electronics.


title: First Resolved Dust Continuum Measurements of Individual Giant Molecular Clouds in the Andromeda Galaxy

abstract: In our local Galactic neighborhood, molecular clouds are best studied using a combination of dust measurements, to determine robust masses, sizes and internal structures of the clouds, and molecular-line observations to determine cloud kinematics and chemistry. We present here the first results of a program designed to extend such studies to nearby galaxies beyond the Magellanic Clouds. Utilizing the wideband upgrade of the Submillimeter Array (SMA) at 230 GHz we have obtained the first continuum detections of the thermal dust emission on sub-GMC scales ($\sim$ 15 pc) within the Andromeda galaxy (M31). These include the first resolved continuum detections of dust emission from individual GMCs beyond the Magellanic Clouds. Utilizing a powerful capability of the SMA, we simultaneously recorded CO(2-1) emission with identical $(u,\,v)$ coverage, astrometry and calibration, enabling the first measurements of the CO conversion factor, $\alpha_{\rm\,CO(2-1)}$, toward individual GMCs across an external galaxy. Our direct measurement yields an average CO--to--dust mass conversion factor of $\alpha^\prime_{\rm CO-dust} = 0.042\pm0.018$ $M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the $J= 2-1$ transition. This value does not appear to vary with galactocentric radius. Assuming a constant gas-to-dust ratio of 136, the resulting $\alpha_{\rm CO}$ $=$ 5.7 $\pm$ 2.4 $M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the 2-1 transition is in excellent agreement with that of Milky Way GMCs, given the uncertainties. Finally, using the same analysis techniques, we compare our results with observations of the local Orion molecular clouds, placed at the distance of M31 and simulated to appear as they would if observed by the SMA.


title: Towards Coding for Human and Machine Vision: A Scalable Image Coding Approach

abstract: The past decades have witnessed the rapid development of image and video coding techniques in the era of big data. However, the signal fidelity-driven coding pipeline design limits the capability of the existing image/video coding frameworks to fulfill the needs of both machine and human vision. In this paper, we come up with a novel image coding framework by leveraging both the compressive and the generative models, to support machine vision and human perception tasks jointly. Given an input image, the feature analysis is first applied, and then the generative model is employed to perform image reconstruction with features and additional reference pixels, in which compact edge maps are extracted in this work to connect both kinds of vision in a scalable way. The compact edge map serves as the basic layer for machine vision tasks, and the reference pixels act as a sort of enhanced layer to guarantee signal fidelity for human vision. By introducing advanced generative models, we train a flexible network to reconstruct images from compact feature representations and the reference pixels. Experimental results demonstrate the superiority of our framework in both human visual quality and facial landmark detection, which provide useful evidence on the emerging standardization efforts on MPEG VCM (Video Coding for Machine).


title: On the Metallicity Gradient in the Galactic Disk

abstract: The problem of the chemical composition gradient in the galactic disk is studied based on a sample of metallicity estimates of open star clusters, using Gaia DR2-improved distance estimates. A clearly non-monotonic variation was observed in the average metallicity of clusters with increasing galactocentric distance. One can clearly see the metallicity jump of 0.22 in [Fe/H] at a Galactocentric distance of about 9.5 kpc, which appears to be linked to the outer boundary of the thinnest and youngest component of the galactic disk. The absence of a significant metallicity gradient in the internal (R<9 kpc) and external (R>10 kpc) regions of the disk demonstrates the absence of noticeable metal enrichment at times of the order of the ages corresponding to those of the disk regions under consideration. Observational data show that the disk experiences noticeable metal enrichment only during the starburst epochs. No significant dependence was found between the average metallicity and the age of the clusters.


title: A blazar as the likely counterpart to 4FGL J0647.7-4418 instead of a gamma-ray binary

abstract: The persistent gamma-ray source 4FGL J0647.7-4418 is tentatively associated in the latest Fermi catalogue with the sub-dwarf O-type X-ray binary HD 49798. However, an AGN candidate is also mentioned as an alternative identification in updated versions of the catalogue accompanying paper. If the first association were correct, this would add HD 49798 to the handful of currently known gamma-ray binaries and, therefore, represent a significant breakthrough not only because of a new member addition, but also because of the apparent white dwarf companion in this system. Despite these perspectives, here we show that the stellar association is likely wrong and that the proposed AGN object, well inside the Fermi 95% confidence ellipse, is a more conceivable counterpart candidate to the Fermi source due to its strong blazar similarities.


title: The Effect of a positive cosmological constant on the bounce of Loop Quantum Cosmology

abstract: We provide an analytical solution to the quantum dynamics of a flat Friedmann-Lema\^itre- Robertson-Walker model with a massless scalar field in the presence of a small and positive cosmological constant, in the context of Loop Quantum Cosmology. We use a perturbative treatment with respect to the model without a cosmological constant, which is exactly solvable. Our solution is approximate, but it is precisely valid at the high curvature regime where quantum gravity corrections are important. We compute explicitly the evolution of the expectation value of the volume. For semiclassical states characterized by a Gaussian spectral profile, the introduction of a positive cosmological constant displaces the bounce of the solvable model to lower volumes and to higher values of the scalar field. These displacements are state dependent, and in particular, they depend on the peak of the Gaussian profile, which measures the momentum of the scalar field. Moreover, for those semiclassical states, the bounce remains symmetric, as in the vanishing cosmological constant case. However, we show that the behavior of the volume is more intricate for generic states, leading in general to a non-symmetric bounce.


title: Influence of vacancy defects in solid high-order harmonic generation

abstract: The present work discusses the impact of vacancy defects in solid high-order harmonic generation. The total energy cutoff of the high-order harmonic spectrum increases as a function of concentration of vacancy defects, and the total spectrum gradually turns into a single slanted spectrum without having an abrupt transition between primary and secondary plateaus. The spectral intensity of the below-band-gap harmonics in a solid with vacancy defects is enhanced significantly in comparison to the harmonics in a pristine solid. The changes in the harmonic spectra are understood in terms of their effective band structures. The presence of vacancy defects breaks the translational symmetry of the unit cell locally. As a consequence of this, new defect states appear, which open additional paths for the electron dynamics. The ill-resolved electron trajectories in the Gabor profile confirm the interference of additional paths. Moreover, the single slanted high-order harmonic spectrum carries a unique signature of vacancy defects in comparison to the high-order harmonic spectrum corresponding to solids with defects such as underdoping or overdoping.

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


title: Online Memorization of Random Firing Sequences by a Recurrent Neural Network

abstract: This paper studies the capability of a recurrent neural network model to memorize random dynamical firing patterns by a simple local learning rule. Two modes of learning/memorization are considered: The first mode is strictly online, with a single pass through the data, while the second mode uses multiple passes through the data. In both modes, the learning is strictly local (quasi-Hebbian): At any given time step, only the weights between the neurons firing (or supposed to be firing) at the previous time step and those firing (or supposed to be firing) at the present time step are modified. The main result of the paper is an upper bound on the probability that the single-pass memorization is not perfect. It follows that the memorization capacity in this mode asymptotically scales like that of the classical Hopfield model (which, in contrast, memorizes static patterns). However, multiple-rounds memorization is shown to achieve a higher capacity (with a nonvanishing number of bits per connection/synapse). These mathematical findings may be helpful for understanding the functions of short-term memory and long-term memory in neuroscience.


title: GRIDS: Interactive Layout Design with Integer Programming

abstract: Grid layouts are used by designers to spatially organise user interfaces when sketching and wireframing. However, their design is largely time consuming manual work. This is challenging due to combinatorial explosion and complex objectives, such as alignment, balance, and expectations regarding positions. This paper proposes a novel optimisation approach for the generation of diverse grid-based layouts. Our mixed integer linear programming (MILP) model offers a rigorous yet efficient method for grid generation that ensures packing, alignment, grouping, and preferential positioning of elements. Further, we present techniques for interactive diversification, enhancement, and completion of grid layouts (Figure 1). These capabilities are demonstrated using GRIDS1, a wireframing tool that provides designers with real-time layout suggestions. We report findings from a ratings study (N = 13) and a design study (N = 16), lending evidence for the benefit of computational grid generation during early stages of design.


title: First observation of electric-quadrupole infrared transitions in water vapour

abstract: Molecular absorption of infrared radiation is generally due to ro-vibrational electric-dipole transitions. Electric-quadrupole transitions may still occur, but they are typically a million times weaker than electric-dipole transitions, rendering their observation extremely challenging. In polyatomic or polar diatomic molecules, ro-vibrational quadrupole transitions have never been observed. Here, we report the first direct detection of quadrupole transitions in water vapor. The detected quadrupole lines have intensity largely above the standard dipole intensity cut-off of spectroscopic databases and thus are important for accurate atmospheric and astronomical remote sensing.


title: New Derivation for Gaussian Mixture Model Parameter Estimation: MM Based Approach

abstract: In this letter, we revisit the problem of maximum likelihood estimation (MLE) of parameters of Gaussian Mixture Model (GMM) and show a new derivation for its parameters. The new derivation, unlike the classical approach employing the technique of expectation-maximization (EM), is straightforward and doesn't invoke any hidden or latent variables and calculation of the conditional density function. The new derivation is based on the approach of minorization-maximization and involves finding a tighter lower bound of the log-likelihood criterion. The update steps of the parameters, obtained via the new derivation, are same as the update steps obtained via the classical EM algorithm.


title: Strong linkage for function fields of surfaces

abstract: Over a global field any finite number of central simple algebras of exponent dividing $m$ is split by a common cyclic field extension of degree $m$. We show that the same property holds for function fields of two-dimensional excellent schemes over a henselian local domain of dimension one or two with algebraically closed residue field.


title: Architecture and Security of SCADA Systems: A Review

abstract: Pipeline bursting, production lines shut down, frenzy traffic, trains confrontation, nuclear reactor shut down, disrupted electric supply, interrupted oxygen supply in ICU - these catastrophic events could result because of an erroneous SCADA system/ Industrial Control System(ICS). SCADA systems have become an essential part of automated control and monitoring of many of the Critical Infrastructures (CI). Modern SCADA systems have evolved from standalone systems into sophisticated complex, open systems, connected to the Internet. This geographically distributed modern SCADA system is vulnerable to threats and cyber attacks. In this paper, we first review the SCADA system architectures that have been proposed/implemented followed by attacks on such systems to understand and highlight the evolving security needs for SCADA systems. A short investigation of the current state of intrusion detection techniques in SCADA systems is done , followed by a brief study of testbeds for SCADA systems. The cloud and Internet of things (IoT) based SCADA systems are studied by analysing the architecture of modern SCADA systems. This review paper ends by highlighting the critical research problems that need to be resolved to close the gaps in the security of SCADA systems.


title: Block-Lanczos density-matrix renormalization-group approach to spin transport in Heisenberg chains coupled to leads

abstract: We adapt the block-Lanczos density-matrix renormalization-group technique to study the spin transport in a spin chain coupled to two non-interacting fermionic leads. As an example, we consider leads described by two-dimensional tight-binding models on a square lattice. Although the simulations are carried out using a chain representation of the leads, observables in the original two-dimensional lattice can be calculated by reversing the block-Lanczos transformation. This is demonstrated for leads with Rashba spin-orbit coupling.


title: Knotted Portals in Virtual Reality

abstract: This article describes a software named "KnotPortal" for visualization of branched covers of knots based on an idea by Bill Thurston to view a knot as a portal to other universes. Our implementation allows users to explore these knotted portals either on a screen or in virtual reality using a head-mounted device with room-tracking. This allows users not only to see these glued different worlds, but experience them by being able to walk through the portals. This software can be used to enable students to learn about knots, gluing, (branched) covers, or just to have a fun look at portals.


title: Performance of the ALICE electromagnetic calorimeters in LHC Runs 1 and 2 and upgrade projects

abstract: ALICE (A Large Ion Collider Experiment) at the LHC is designed for studies of nuclear matter at extreme temperatures and energy densities, so called Quark-Gluon Plasma (QGP). Two detectors for measurements of electromagnetic signals, high-granularity photon spectrometer PHOS and large acceptance electromagnetic calorimeter EMCal/DCal, are incorporated into the ALICE apparatus. In order to enhance the data sample of the high-energy part of the spectra, a dedicated triggers on high energetic photons or jets are used. High precision measurements of neutral mesons require fine energy and timing calibrations. Photon signal purity is improved using dedicated photon PID criteria based on cluster shape and anti-track matching with the ALICE central tracking system, as well as a special detector located in front of PHOS called CPV. PHOS and EMCal participated in LHC Run 1 (2009-2013) and Run 2 (2015-2018) during data taking periods of pp, p-Pb and Pb-Pb collisions. We give an overview of their performance in Runs 1 and 2 in low and high multiplicity environments as well as the upgrade plans for future LHC runs.


title: Observed $\Omega_b$ spectrum and meson-baryon molecular states

abstract: We observe that four peaks seen in the high energy part of the $\Omega_b$ spectrum of the recent LHCb experiment are in remarkable agreement with predictions made for molecular $\Omega_b$ states stemming from the meson-baryon interaction, with an approach that applied to the $\Omega_c$ states gave rise to three states in good agreement with experiment in masses and widths. While the statistical significance of the peaks prevents us from claims of states at the present time, the agreement found should be an incentive to look at this experiment with increased statistics to give an answer to this suggestive idea.


title: Immediate extensions of valuation rings and ultrapowers

abstract: We describe the immediate extensions of a one dimensional valuation ring $V$ which could be embedded in some separation of a ultrapower of $V$ with respect to a certain ultrafilter. For such extensions holds a kind of Artin's approximation.


title: Reflections and powers of multisorted minions

abstract: Classes of multisorted minions closed under extensions, reflections, and direct powers are considered from a relational point of view. As a generalization of a result of Barto, Opr\v{s}al, and Pinsker, the closure of a multisorted minion is characterized in terms of constructions on multisorted relation pairs which are invariant for minions.


title: Privacy-Preserving Deep Learning Computation for Geo-Distributed Medical Big-Data Platforms

abstract: This paper proposes a distributed deep learning framework for privacy-preserving medical data training. In order to avoid patients' data leakage in medical platforms, the hidden layers in the deep learning framework are separated and where the first layer is kept in platform and others layers are kept in a centralized server. Whereas keeping the original patients' data in local platforms maintain their privacy, utilizing the server for subsequent layers improves learning performance by using all data from each platform during training.


title: Defiducialization: Providing Experimental Measurements for Accurate Fixed-Order Predictions

abstract: An experimental procedure is proposed to perform measurements of differential cross sections which can be compared to fixed-order QCD predictions with improved accuracy. The procedure can be applied to the Drell-Yan cross-section measurements which are differential in the boson transverse momentum. An example analysis is performed using the ATLAS measurement of the $Z$-boson production cross section at center-of-mass energy of $8$ TeV. The resulting full phase space measurement of the cross section differential in the boson rapidity is compared to theoretical predictions computed with next-to-next-to leading-order accuracy in QCD.


title: Billiards in ellipses revisited

abstract: We prove some recent experimental observations of D. Reznik concerning periodic billiard orbits in ellipses. For example, the sum of cosines of the angles of a periodic billiard polygon remains constant in the one-parameter family of such polygons (that exist due to the Poncelet porism). In our proofs, we use geometric and complex analytic methods.


title: Multipass SAR Interferometry Based on Total Variation Regularized Robust Low Rank Tensor Decomposition

abstract: Multipass SAR interferometry (InSAR) techniques based on meter-resolution spaceborne SAR satellites, such as TerraSAR-X or COSMO-Skymed, provide 3D reconstruction and the measurement of ground displacement over large urban areas. Conventional method such as Persistent Scatterer Interferometry (PSI) usually requires a fairly large SAR image stack (usually in the order of tens), in order to achieve reliable estimates of these parameters. Recently, low rank property in multipass InSAR data stack was explored and investigated in our previous work. By exploiting this low rank prior, more accurate estimation of the geophysical parameters can be achieved, which in turn can effectively reduce the number of interferograms required for a reliable estimation. Based on that, this paper proposes a novel tensor decomposition method in complex domain, which jointly exploits low rank and variational prior of the interferometric phase in InSAR data stacks. Specifically, a total variation (TV) regularized robust low rank tensor decomposition method is exploited for recovering outlier-free InSAR stacks. We demonstrate that the filtered InSAR data stacks can greatly improve the accuracy of geophysical parameters estimated from real data. Moreover, this paper demonstrates for the first time in the community that tensor-decomposition-based methods can be beneficial for large-scale urban mapping problems using multipass InSAR. Two TerraSAR-X data stacks with large spatial areas demonstrate the promising performance of the proposed method.


title: Generation of rotating solutions in Einstein-scalar gravity

abstract: We revisit the problem of the rotating generalization of the Fisher-Janis-Newman-Winicour solution of the minimal Einstein-scalar theory proving that the two previously proposed solutions do not satisfy the equations of motion. We also derive several new spinning solutions with the scalar charge, which are endowed with oblate deformation.


title: Breakdown of the tensor component in the Skyrme energy density functional

abstract: The tensor force, as an important component of strong nuclear force, generates a variety of intriguing effects ranging from few-body systems to neutron stars. It is responsible for the nucleon-nucleon correlation beyond mean-field approximation, and is accordingly proved to play no role in the standard Skyrme energy density functionals in the present work. Therefore, the Skyrme's original tensor interaction that is extensively-employed presently is invalid. As an alternative strategy, we introduced a central interaction, i.e., the $\bm{\sigma }_{1}\cdot \bm{\sigma }_{2}$ term, to improve the description of experimental single-particle structure, and to address its effect, we established two Skyrme interactions IMP1 and IMP2 complemented by the calibrated charge-violating interactions. The central $\bm{\sigma }_{1}\cdot \bm{\sigma }_{2}$ interaction turns out to substantially improve the description of shell evolution in Sn isotopes and $N=82$ isotones.


title: Computing Curved Area Labels in Near-Real Time

abstract: In the Area Labeling Problem one is after placing the label of a geographic area. Given the outer boundary of the area and an optional set of holes. The goal is to find a label position such that the label spans the area and is conform to its shape. The most recent research in this field from Barrault in 2001 proposes an algorithm to compute label placements based on curved support lines. His solution has some drawbacks as he is evaluating many very similar solutions. Furthermore he needs to restrict the search space due to performance issues and therefore might miss interesting solutions. We propose a solution that evaluates the search space more broadly and much more efficient. To achieve this we compute a skeleton of the polygon. The skeleton is pruned such that edges close to the boundary polygon are removed. In the so pruned skeleton we choose a set of candidate paths to be longest distinct subpaths of the graph. Based on these candidates the label support lines are computed and the label positions evaluated. Keywords: Area lettering \and Automated label placement \and Digital cartography \and Geographic information sciences \and Geometric Optimization.


title: Sunspot Drawings at Kodaikanal Observatory: A Representative Results on Hemispheric Sunspot Numbers and Area Measurements

abstract: The importance of the periodicity in sunspot appearance was well recognized by the mid of 19th century. Several observatories around the globe have made the record of sunspots in the form of drawings and preserved them safely for posterity. At the Kodaikanal Observatory (KO), the sunspot observations have begun in 1905. In those times observations were recorded on photographic plates and after the development of those plates in the laboratory, the drawings of the same were made on the Stonyhurst grids. In these drawings, called {\em sun charts}, different features on the sun's disk, e.g., sunspots, plages, filaments, prominences, etc. were clearly identified and visually marked with different colors. We have collected 111 years of sunspot drawing spanning over 10 solar cycles. These sunspot drawings were carefully stitched to make bound volumes, each for every 6-months. The drawings are kept at the Kodaikanal library for scientific use and analysis. In this article, we describe briefly the process of drawing, methods of counting sunspot numbers and measurement of sunspot area using square grids. We have collected the data for the Northern and Southern hemispheres separately. From the collected data, we compute the sunspot number and area and compare it with Royal Greenwich Observatory (RGO) and Sunspot Index and Long-term Solar Observations (SILSO) data. The results show that the measurement of the sunspot number is underestimated by about 40\%. The KO monthly averaged sunspot number data of both hemispheres is normalized with the RGO monthly averaged total sunspot number data.


title: Convergence of the Hesse-Koszul flow on compact Hessian manifolds

abstract: We study the long time behavior of the Hesse-Koszul flow on compact Hessian manifolds. When the first affine Chern class is negative, we prove that the flow converges to the unique Hesse-Einstein metric. We also derive a convergence result for a twisted Hesse-Koszul flow on any compact Hessian manifold. These results give alternative proofs for the existence of the unique Hesse-Einstein metric by Cheng-Yau and Caffarelli-Viaclovsky as well as the real Calabi theorem by Cheng-Yau, Delano\"e and Caffarelli-Viaclovsky.


title: Killing Stubborn Mutants with Symbolic Execution

abstract: We introduce SeMu, a Dynamic Symbolic Execution technique that generates test inputs capable of killing stubborn mutants (killable mutants that remain undetected after a reasonable amount of testing). SeMu aims at mutant propagation (triggering erroneous states to the program output) by incrementally searching for divergent program behaviours between the original and the mutant versions. We model the mutant killing problem as a symbolic execution search within a specific area in the programs' symbolic tree. In this framework, the search area is defined and controlled by parameters that allow scalable and cost-effective mutant killing. We integrate SeMu in KLEE and experimented with Coreutils (a benchmark frequently used in symbolic execution studies). Our results show that our modelling plays an important role in mutant killing. Perhaps more importantly, our results also show that, within a two-hour time limit, SeMu kills 37% of the stubborn mutants, where KLEE kills none and where the mutant infection strategy (strategy suggested by previous research) kills 17%.


title: Neural Network Tomography

abstract: Network tomography, a classic research problem in the realm of network monitoring, refers to the methodology of inferring unmeasured network attributes using selected end-to-end path measurements. In the research community, network tomography is generally investigated under the assumptions of known network topology, correlated path measurements, bounded number of faulty nodes/links, or even special network protocol support. The applicability of network tomography is considerably constrained by these strong assumptions, which therefore frequently position it in the theoretical world. In this regard, we revisit network tomography from the practical perspective by establishing a generic framework that does not rely on any of these assumptions or the types of performance metrics. Given only the end-to-end path performance metrics of sampled node pairs, the proposed framework, NeuTomography, utilizes deep neural network and data augmentation to predict the unmeasured performance metrics via learning non-linear relationships between node pairs and underlying unknown topological/routing properties. In addition, NeuTomography can be employed to reconstruct the original network topology, which is critical to most network planning tasks. Extensive experiments using real network data show that comparing to baseline solutions, NeuTomography can predict network characteristics and reconstruct network topologies with significantly higher accuracy and robustness using only limited measurement data.


title: Binary and Multitask Classification Model for Dutch Anaphora Resolution: Die/Dat Prediction

abstract: The correct use of Dutch pronouns 'die' and 'dat' is a stumbling block for both native and non-native speakers of Dutch due to the multiplicity of syntactic functions and the dependency on the antecedent's gender and number. Drawing on previous research conducted on neural context-dependent dt-mistake correction models (Heyman et al. 2018), this study constructs the first neural network model for Dutch demonstrative and relative pronoun resolution that specifically focuses on the correction and part-of-speech prediction of these two pronouns. Two separate datasets are built with sentences obtained from, respectively, the Dutch Europarl corpus (Koehn 2015) - which contains the proceedings of the European Parliament from 1996 to the present - and the SoNaR corpus (Oostdijk et al. 2013) - which contains Dutch texts from a variety of domains such as newspapers, blogs and legal texts. Firstly, a binary classification model solely predicts the correct 'die' or 'dat'. The classifier with a bidirectional long short-term memory architecture achieves 84.56% accuracy. Secondly, a multitask classification model simultaneously predicts the correct 'die' or 'dat' and its part-of-speech tag. The model containing a combination of a sentence and context encoder with both a bidirectional long short-term memory architecture results in 88.63% accuracy for die/dat prediction and 87.73% accuracy for part-of-speech prediction. More evenly-balanced data, larger word embeddings, an extra bidirectional long short-term memory layer and integrated part-of-speech knowledge positively affects die/dat prediction performance, while a context encoder architecture raises part-of-speech prediction performance. This study shows promising results and can serve as a starting point for future research on machine learning models for Dutch anaphora resolution.


title: New hierarchy of multiple soliton solutions for the (2+1)-dimensional Sawada-Kotera equation

abstract: A new transformation $u=4 ({\rm ln}f)_x$ that can formulate a quintic linear equation and a pair of Hirota's bilinear equations for the (2+1)-dimensional Sawada-Kotera (2DSK) equation is reported firstly, which enables one to obtain a new hierarchy of multiple soliton solutions of the 2DSK equation. It tells a crucial fact that a nonlinear partial differential equation could possess two hierarchies of multiple soliton solutions and the 2DSK equation is the first and only one found in this paper. The quintic linear equation is solved by a pair of Hirota's bilinear equations, of which one is the (2+1)-dimensional bilinear SK equation obtained by $u=2 ({\rm ln}f)_{x}$, and the other is the bilinear KdV equation. The (1+1)-dimensional SK equation does not possess this property. As another example, a (3+1)-dimensional nonlinear partial differential equation possessing a pair of Hirota's bilinear equations, however only bearing one hierarchy of multiple soliton solutions is studied.


title: Regular $3$-polytopes of order $2^np$

abstract: In [Problems on polytopes, their groups, and realizations, Periodica Math. Hungarica 53 (2006) 231-255] Schulte and Weiss proposed the following problem: {\em Characterize regular polytopes of orders $2^np$ for $n$ a positive integer and $p$ an odd prime}. In this paper, we first prove that if a $3$-polytope of order $2^np$ has Schl\"afli type $\{k_1, k_2\}$, then $p \mid k_1$ or $p \mid k_2$. This leads to two classes, up to duality, for the Schl\"afli type, namely Type (1) where $k_1=2^sp$ and $k_2=2^t$ and Type (2) where $k_1=2^sp$ and $k_2=2^tp$. We then show that there exists a regular $3$-polytope of order $2^np$ with Type (1) when $s\geq 2$, $t\geq 2$ and $n\geq s+t+1$ coming from a general construction of regular $3$-polytopes of order $2^n\ell_1\ell_2$ with Schl\"afli type $\{2^s\ell_1,2^t\ell_2\}$ when $s, t \geq 2$, $n\geq s+t+1$, and both $\ell_1$ and $\ell_2$ are odd. For Type (2), we prove that there exists a regular $3$-polytope of order $2^n\cdot 3$ with Schl\"afli type $\{6, 6\}$ when $n \ge 5$ coming from a general construction of regular $3$-polytopes of Schl\"afli type $\{6,6\}$ with orders $192m^3$, $384m^3$ or $768m^3$, for any positive integer $m$.


title: Independent Domination in Subcubic Graphs

abstract: A set $S$ of vertices in a graph $G$ is a dominating set if every vertex not in $S$ is adjacent to a vertex in $S$. If, in addition, $S$ is an independent set, then $S$ is an independent dominating set. The independent domination number $i(G)$ of $G$ is the minimum cardinality of an independent dominating set in $G$. In 2013 Goddard and Henning [Discrete Math 313 (2013), 839--854] conjectured that if $G$ is a connected cubic graph of order $n$, then $i(G) \le \frac{3}{8}n$, except if $G$ is the complete bipartite graph $K_{3,3}$ or the $5$-prism $C_5 \, \Box \, K_2$. Further they construct two infinite families of connected cubic graphs with independent domination three-eighths their order. They remark that perhaps it is even true that for $n > 10$ these two families are only families for which equality holds. In this paper, we provide a new family of connected cubic graphs $G$ of order $n$ such that $i(G) = \frac{3}{8}n$. We also show that if $G$ is a subcubic graph of order $n$ with no isolated vertex, then $i(G) \le \frac{1}{2}n$, and we characterize the graphs achieving equality in this bound.


title: Type Ia Supernova Sub-classes and Progenitor Origin

abstract: This paper presents a short review on the current state of SN Ia progenitor origin. Type Ia supernova explosions are observed to be widely diverse in peak luminosity, lightcurve width and shape, spectral features, and host stellar population environment. In the last decade alone, theoretical simulations and observational data have come together to seriously challenge the long-standing paradigm that all SNe Ia arise from explosions of Chandrasekhar mass white dwarfs. In this review I highlight some of the major developments (and changing views) of our understanding of the nature of SN Ia progenitor systems. I give a brief overview of binary star configurations and their plausible explosion mechanisms, and infer links between some of the various (observationally-categorized) SN Ia sub-classes and their progenitor origins from a theoretical standpoint.


title: Counting King Permutations on the Cylinder

abstract: We call a permutation $\sigma=[\sigma_1,\dots,\sigma_n] \in S_n$ a {\em cylindrical king permutation} if $ |\sigma_i-\sigma_{i+1}|>1$ for each $1\leq i \leq n-1$ and $|\sigma_1-\sigma_n|>1$. We present some results regarding the distribution of the cylindrical king permutations, including some interesting recursions. We also calculate their asymptotic proportion in the set of the 'king permutations', i.e. the ones which satisfy only the first of the two conditions above. With this aim we define a new parameter on permutations, namely, the number of {\em cyclic bonds} which is a modification of the number of bonds. In addition, we present some results regarding the distribution of this parameter.


title: Bond-based peridynamics does not converge to hyperelasticity as the horizon goes to zero

abstract: Bond-based peridynamics is a nonlocal continuum model in Solid Mechanics in which the energy of a deformation is calculated through a double integral involving pairs of points in the reference and deformed configurations. It is known how to calculate the {\Gamma}-limit of this model when the horizon (maximum interaction distance between the particles) tends to zero, and the limit turns out to be a (local) vector variational problem defined in a Sobolev space, of the type appearing in (classical) hyperelasticity. In this paper, we impose frame-indifference and isotropy in the model and find that very few hyperelastic functionals are {\Gamma}-limits of the bond-based peridynamics model. In particular, Mooney-Rivlin materials are not recoverable through this limit procedure.


title: Generative Pseudo-label Refinement for Unsupervised Domain Adaptation

abstract: We investigate and characterize the inherent resilience of conditional Generative Adversarial Networks (cGANs) against noise in their conditioning labels, and exploit this fact in the context of Unsupervised Domain Adaptation (UDA). In UDA, a classifier trained on the labelled source set can be used to infer pseudo-labels on the unlabelled target set. However, this will result in a significant amount of misclassified examples (due to the well-known domain shift issue), which can be interpreted as noise injection in the ground-truth labels for the target set. We show that cGANs are, to some extent, robust against such "shift noise". Indeed, cGANs trained with noisy pseudo-labels, are able to filter such noise and generate cleaner target samples. We exploit this finding in an iterative procedure where a generative model and a classifier are jointly trained: in turn, the generator allows to sample cleaner data from the target distribution, and the classifier allows to associate better labels to target samples, progressively refining target pseudo-labels. Results on common benchmarks show that our method performs better or comparably with the unsupervised domain adaptation state of the art.


title: Understanding video streaming algorithms in the wild

abstract: While video streaming algorithms are a hot research area, with interesting new approaches proposed every few months, little is known about the behavior of the streaming algorithms deployed across large online streaming platforms that account for a substantial fraction of Internet traffic. We thus study adaptive bitrate streaming algorithms in use at 10 such video platforms with diverse target audiences. We collect traces of each video player's response to controlled variations in network bandwidth, and examine the algorithmic behavior: how risk averse is an algorithm in terms of target buffer; how long does it takes to reach a stable state after startup; how reactive is it in attempting to match bandwidth versus operating stably; how efficiently does it use the available network bandwidth; etc. We find that deployed algorithms exhibit a wide spectrum of behaviors across these axes, indicating the lack of a consensus one-size-fits-all solution. We also find evidence that most deployed algorithms are tuned towards stable behavior rather than fast adaptation to bandwidth variations, some are tuned towards a visual perception metric rather than a bitrate-based metric, and many leave a surprisingly large amount of the available bandwidth unused.


title: Dynamics of the Taylor shift on Bergman spaces

abstract: The Taylor (backward) shift on Bergman spaces $A^p(\om)$ for general open sets $\om$ in the extended complex plane shows rich variety concerning its dynamical behaviour. Different aspects are worked out, where in the case $p<2$ a recent result of Bayart and Matheron plays a central role.


title: Comparison of conformal gravitation theories by solutions analogous to Reissner-Nordstr\"om ones

abstract: An exact solution was produced for one of the versions of conformal -invariant gravitation theories (conformal geometrodynamics - CG) for the body with a mass and an electric charge. The solution is analogous to the Reissner-Nordstr\"{o}m solution in the general relativity theory. Its similar solution was earlier produced by Mannheit and Kazanas for the different version of the conformal-invariant gravitation theory -- Weyl conformal gravitation theory (Weyl CGt). Similar solutions do not comply. Issues of conformity with regard to the similar solutions are discussed. A conclusion is made that only the similar solution produced from conformal geometrodynamics (CG) equations corresponds to electrodynamic equations.


title: Tetrahedral Symmetry A$_{\bf 4}$ in Anti-SU(5) GUT

abstract: We construct a flavor model in an anti-SU(5) GUT with a tetrahedral symmetry $A_4$. We choose a basis where $Q_{text{em}}=-\frac13$ quarks and charged leptons are already mass eigenstates. This choice is possible from the $A_4$ symmetry. Then, matter representation $\overline{10}_{-1}^{\rm\, matter}$ contains both a quark doublet and a heavy neutrino $N$, which enables us to use the $A_4$ symmetry to both $Q_{text{em}}=+\frac23$ quark masses and neutrino masses (through the see-saw via $N$). This is made possible because the anti-SU(5) breaking is achieved by the Higgs fields transforming as anti-symmetric representations of SU(5), $\overline{10}_{-1}^H\oplus 10_{+1}^H$, reducing the rank-5 anti-SU(5) group down to the rank-4 standard model group \smg. For possible mass matrices, the $A_4$ symmetry predictions on mass matrices at field theory level are derived. Finally, an illustration from string compactification is presented.


title: Optical spectroscopy of type-2 LINERs

abstract: Type-2 Low-ionization Narrow Emission-line Regions (LINERs) have been optically classified with the Palomar data as not presenting a broad component in the Balmer emission lines associated to the Broad Line Region (BLR) of the Active Galactic Nuclei (AGN). We unveil the presence of different kinematic components of emission lines in the nuclear region of a sample of local (z$\leq$0.022) type-2 LINERs. We focus on the analysis of their true nature by means of the detection (or not) of a broad component originated in the BLR. Additionally, we search for the possible presence of non-rotational motions such as outflows. We have applied a decomposition of the nuclear emission lines using a spectroscopic analysis of the optical spectra of 9 type-2 LINERs of intermediate-resolution data retrieved from the Hubble Space Telescope (HST) archive. The study is completed with archival spectra from the Double Spectrograph from the Palomar Observatory. The emission line fitting reveals the presence of a broad component associated to the BLR in 6 out of the 9 galaxies for the space-based data, and for 2 out of the 8 from the ground-based spectra. The velocity dispersion of two galaxies (NGC 4486 and NGC 4594) measured in HST/STIS data suggest the presence of outflows. The results indicate that the spatial resolution plays a major role in the detection of the BLR, as it appears diluted in the ground-based data (even after removing stellar contribution). This is also true for the emission line diagnostics, as the contaminant light contributes to lower emission line ratios towards the star-forming area of standard BPTs. We propose to reclassify NGC4594 as a type-1 LINER, since a BLR component is seen in both space- and ground-based spectra. We find ambiguous results for the BLR component of NGC 4486. The modest outflow detection in our sample may indicate that they are not as frequent as for type-1 LINERs.


title: On Decoding Using Codewords of the Dual Code

abstract: We present novel decoding schemes for hard and soft decision decoding of block codes using the minimal weight codewords of the dual code. The decoding schemes will be described for cyclic codes where polynomials can be used, however, the modification for non-cyclic codes is possible and straight forward. The hard decision decoding calculates syndrome polynomials which are the product of the received polynomial with dual codewords. Proper cyclic shifts of these syndrome polynomials are obtained and the non-zero positions are counted componentwise for these shifts. The values of this counting are a reliability measure and can be used for locating the error and also the non-error positions. This reliability measure is the basis for various variants of hard decision decoding algorithms. Decoding schemes with iterative error reduction are possible as well as information set decoding using the inherent reliability information of the measure even if there is no reliability information from the channel. Further, we will show how reliability information from the channel can be included in order to obtain soft decision decoding schemes. We derive the relation between bit flipping, believe propagation, and majority logic decoding to the novel schemes. As examples to illustrate the functioning we use BCH and Reed-Muller codes as examples for binary codes, and RS codes for non-binary codes. Besides the Plotkin construction we recall a known result that Reed-Muller codes punctured by one position are cyclic and thus, are equivalent to special cases of BCH codes. Simulation results for hard and soft decision decoding will be given for several examples and compared with results from literature. Finally, we analyze the soft decision decoding of the Plotkin construction and derive that one of the two codes uses a $3$ dB better channel (also known as channel polarization).


title: Expected Improvement versus Predicted Value in Surrogate-Based Optimization

abstract: Surrogate-based optimization relies on so-called infill criteria (acquisition functions) to decide which point to evaluate next. When Kriging is used as the surrogate model of choice (also called Bayesian optimization), one of the most frequently chosen criteria is expected improvement. We argue that the popularity of expected improvement largely relies on its theoretical properties rather than empirically validated performance. Few results from the literature show evidence, that under certain conditions, expected improvement may perform worse than something as simple as the predicted value of the surrogate model. We benchmark both infill criteria in an extensive empirical study on the `BBOB' function set. This investigation includes a detailed study of the impact of problem dimensionality on algorithm performance. The results support the hypothesis that exploration loses importance with increasing problem dimensionality. A statistical analysis reveals that the purely exploitative search with the predicted value criterion performs better on most problems of five or higher dimensions. Possible reasons for these results are discussed. In addition, we give an in-depth guide for choosing the infill criteria based on prior knowledge about the problem at hand, its dimensionality, and the available budget.


title: Shape optimization of a weighted two-phase Dirichlet eigenvalue

abstract: Let $m$ be a bounded function and $\alpha$ a nonnegative parameter. This article is concerned with the first eigenvalue $\lambda\_\alpha(m)$ of the drifted Laplacian type operator $\mathcal L\_m$ given by $\mathcal L\_m(u)= -\operatorname{div} \left((1+\alpha m)\nabla u\right)-mu$ on a smooth bounded domain, with Dirichlet boundary conditions. Assuming uniform pointwise and integral bounds on $m$, we investigate the issue of minimizing $\lambda\_\alpha(m)$ with respect to $m$. Such a problem is related to the so-called "two phase extremal eigenvalue problem" and arises naturally, for instance in population dynamics where it is related to the survival ability of a species in a domain. We prove that unless the domain is a ball, this problem has no "regular" solution. We then provide a careful analysis in the case of a ball by: (1) characterizing the solution among all radially symmetric resources distributions, with the help of a new method involving a homogenized version of the problem; (2) proving in a more general setting, a stability result for the centered distribution of resources with the help of a monotonicity principle for second order shape derivatives which significantly simplifies the analysis.


title: Segregation with Social Linkages: Evaluating Schelling's Model with Networked Individuals

abstract: This paper generalizes the original Schelling (1969, 1971a,b, 2006) model of racial and residential segregation to a context of variable externalities due to social linkages. In a setting in which individuals' utility function is a convex combination of a heuristic function a la Schelling, of the distance to friends, and of the cost of moving, the prediction of the original model gets attenuated: the segregation equilibria are not the unique solutions. While the cost of distance has a monotonic pro-status-quo effect, equivalent to that of models of migration and gravity models, if friends and neighbours are formed following independent processes the location of friends in space generates an externality that reinforces the initial configuration if the distance to friends is minimal, and if the degree of each agent is high. The effect on segregation equilibria crucially depends on the role played by network externalities.


title: Computing Persistent Homology with Various Coefficient Fields in a Single Pass

abstract: This article introduces an algorithm to compute the persistent homology of a filtered complex with various coefficient fields in a single matrix reduction. The algorithm is output-sensitive in the total number of distinct persistent homological features in the diagrams for the different coefficient fields. This computation allows us to infer the prime divisors of the torsion coefficients of the integral homology groups of the topological space at any scale, hence furnishing a more informative description of topology than persistence in a single coefficient field. We provide theoretical complexity analysis as well as detailed experimental results. The code is part of the Gudhi software library.


title: Analysis of the $\Omega_b(6316)$, $\Omega_b(6330)$, $\Omega_b(6340)$ and $\Omega_b(6350)$ with QCD sum rules

abstract: In this article, we introduce an explicit P-wave to construct three-quark currents to study the P-wave $\Omega_b$ states with the full QCD sum rules. The predicted masses have a hierarchy if the same parameters are chosen and favor assigning the $\Omega_b(6316)$, $\Omega_b(6330)$, $\Omega_b(6340)$ and $\Omega_b(6350)$ to be the P-wave $\Omega_b$ states with the $J^P={\frac{3}{2}}^-$, ${\frac{1}{2}}^-$, ${\frac{5}{2}}^-$ and ${\frac{3}{2}}^-$, respectively.


title: LibreSocial: A Peer-to-Peer Framework for Online Social Networks

abstract: Distributed online social networks (DOSNs) were first proposed to solve the problem of privacy, security and scalability. A significant amount of research was undertaken to offer viable DOSN solutions that were capable of competing with the existing centralized OSN applications such as Facebook, LinkedIn and Instagram. This research led to the emergence of the use of peer-to-peer (P2P) networks as a possible solution, upon which several OSNs such as LifeSocial.KOM, Safebook, PeerSoN among others were based. In this paper, we define the basic requirements for a P2P OSN. We then revisit one of the first P2P-based OSNs, LifeSocial.KOM, that is now called LibreSocial, which evolved in the past years to address the challenges of running a completely decentralized social network. Over the course of time, several essential new technologies have been incorporated within LibreSocial for better functionalities. We describe the architecture and each individual component of LibreSocial and point out how LibreSocial meets the basic requirements for a fully functional distributed OSN.


title: Time-Reversal Symmetry and the Structure of Superconducting Order Parameter of Nearly Ferromagnetic Spin-Triplet Superconductor UTe2

abstract: The space-group approach to the wavefunction of a Cooper pair is extended to Shubnikov groups of type III and applied to newly discovered nearly ferromagnetic superconductor UTe2. A nonunitary triplet SOP (superconducting order parameter), manifested experimentally, was obtained for Shubnikov group m'm'm . The nodal structure of triplet nonunitary SOP is investigated and it was shown that in the case of Shubnikov III symmetry nodal structure of ESP (equal spin pairing) state is in some cases unstable, i.e. changes when a complex phase winding is added. The irreducible corepresentations of a Cooper pair wave function corresponding to m'm'm symmetry and to all experimental data for UTe2 are identified.


title: Performance of NavIC for studying the ionosphere at an EIA region in India

abstract: This paper emphasizes on NavIC's performance in ionospheric studies over the Indian subcontinent region. The study is performed using data of one year (2017-18) at IIT Indore, a location near the northern crest of Equatorial Ionization Anomaly (EIA). It has been observed that even without the individual error corrections, the results are within $\pm20\%$ of NavIC VTEC estimates observed over the 1\ensuremath{^{\circ}} x 1\ensuremath{^{\circ}} grid of IPP surrounding the GPS VTEC estimates for most of the time. Additionally, ionospheric response during two distinct geomagnetic storms (September 08 and 28, 2017) at the same location and other IGS stations covering the Indian subcontinent using both GPS and NavIC has also been presented. This analysis revealed similar variations in TEC during the geomagnetic storms of September 2017, indicating the suitability of NavIC to study space weather events along with the ionospheric studies over the Indian subcontinent.

categories: physics.space-ph astro-ph.EP astro-ph.SR physics.geo-ph physics.plasm-ph


title: The Multi-Blade Boron-10-based neutron detector performance using a focusing reflectometer

abstract: The Multi-Blade is a Boron-10-based neutron detector designed for neutron reflectometers and developed for the two instruments (Estia and FREIA) planned for the European Spallation Source in Sweden. A reflectometry demonstrator has been installed at the AMOR reflectometer at the Paul Scherrer Institut (PSI - Switzerland). The setup exploits the Selene guide concept and it can be considered a scaled-down demonstrator of Estia. The results of these tests are discussed. It will be shown how the characteristics of the Multi-Blade detector are features that allow the focusing reflectometry operation mode. Additionally the performance of the Multi-Blade, in terms of rate capability, exceeds current state-of-the-art technology. The improvements with respect to the previous prototypes are also highlighted; from background considerations to the linear and angular uniformity response of the detector.


title: Clustering Approaches for Global Minimum Variance Portfolio

abstract: The only input to attain the portfolio weights of global minimum variance portfolio (GMVP) is the covariance matrix of returns of assets being considered for investment. Since the population covariance matrix is not known, investors use historical data to estimate it. Even though sample covariance matrix is an unbiased estimator of the population covariance matrix, it includes a great amount of estimation error especially when the number of observed data is not much bigger than number of assets. As it is difficult to estimate the covariance matrix with high dimensionality all at once, clustering stocks is proposed to come up with covariance matrix in two steps: firstly, within a cluster and secondly, between clusters. It decreases the estimation error by reducing the number of features in the data matrix. The motivation of this dissertation is that the estimation error can still remain high even after clustering, if a large amount of stocks is clustered together in a single group. This research proposes to utilize a bounded clustering method in order to limit the maximum cluster size. The result of experiments shows that not only the gap between in-sample volatility and out-of-sample volatility decreases, but also the out-of-sample volatility gets reduced. It implies that we need a bounded clustering algorithm so that maximum clustering size can be precisely controlled to find the best portfolio performance.


title: Lie algebra modules which are locally finite over the semi-simple part

abstract: For a finite-dimensional Lie algebra $\mathfrak{L}$ over $\mathbb{C}$ with a fixed Levi decomposition $\mathfrak{L} = \mathfrak{g} \oplus \mathfrak{r}$ where $\mathfrak{g}$ is semi-simple, we investigate $\mathfrak{L}$-modules which decompose, as $\mathfrak{g}$-modules, into a direct sum of simple finite-dimensional $\mathfrak{g}$-modules with finite multiplicities. We call such modules $\mathfrak{g}$-Harish-Chandra modules. We give a complete classification of simple $\mathfrak{g}$-Harish-Chandra modules for the Takiff Lie algebra associated to $\mathfrak{g} = \mathfrak{sl}_2$, and for the Schr\"{o}dinger Lie algebra, and obtain some partial results in other cases. An adapted version of Enright's and Arkhipov's completion functors plays a crucial role in our arguments. Moreover, we calculate the first extension groups of infinite-dimensional simple $\mathfrak{g}$-Harish-Chandra modules and their annihilators in the universal enveloping algebra, for the Takiff $\mathfrak{sl}_2$ and the Schr\"{o}dinger Lie algebra. In the general case, we give a sufficient condition for the existence of infinite-dimensional simple $\mathfrak{g}$-Harish-Chandra modules.


title: How to trap a gradient flow

abstract: We consider the problem of finding an $\varepsilon$-approximate stationary point of a smooth function on a compact domain of $\mathbb{R}^d$. In contrast with dimension-free approaches such as gradient descent, we focus here on the case where $d$ is finite, and potentially small. This viewpoint was explored in 1993 by Vavasis, who proposed an algorithm which, for any fixed finite dimension $d$, improves upon the $O(1/\varepsilon^2)$ oracle complexity of gradient descent. For example for $d=2$, Vavasis' approach obtains the complexity $O(1/\varepsilon)$. Moreover for $d=2$ he also proved a lower bound of $\Omega(1/\sqrt{\varepsilon})$ for deterministic algorithms (we extend this result to randomized algorithms). Our main contribution is an algorithm, which we call gradient flow trapping (GFT), and the analysis of its oracle complexity. In dimension $d=2$, GFT closes the gap with Vavasis' lower bound (up to a logarithmic factor), as we show that it has complexity $O\left(\sqrt{\frac{\log(1/\varepsilon)}{\varepsilon}}\right)$. In dimension $d=3$, we show a complexity of $O\left(\frac{\log(1/\varepsilon)}{\varepsilon}\right)$, improving upon Vavasis' $O\left(1 / \varepsilon^{1.2} \right)$. In higher dimensions, GFT has the remarkable property of being a logarithmic parallel depth strategy, in stark contrast with the polynomial depth of gradient descent or Vavasis' algorithm. In this higher dimensional regime, the total work of GFT improves quadratically upon the only other known polylogarithmic depth strategy for this problem, namely naive grid search. We augment this result with another algorithm, named \emph{cut and flow} (CF), which improves upon Vavasis' algorithm in any fixed dimension.


title: Optimal Control of Perfect Plasticity Part I: Stress Tracking

abstract: The paper is concerned with an optimal control problem governed by the rate-independent system of quasi-static perfect elasto-plasticity. The objective is to optimize the stress field by controlling the displacement at prescribed parts of the boundary. The control thus enters the system in the Dirichlet boundary conditions. Therefore, the safe load condition is automatically fulfilled so that the system admits a solution, whose stress field is unique. This gives rise to a well defined control-to-state operator, which is continuous but not Gateaux-differentiable. The control-to-state map is therefore regularized, first by means of the Yosida regularization and then by a second smoothing in order to obtain a smooth problem. The approximation of global minimizers of the original non-smooth optimal control problem is shown and optimality conditions for the regularized problem are established. A numerical example illustrates the feasibility of the smoothing approach.


title: Closed-loop deep learning: generating forward models with back-propagation

abstract: A reflex is a simple closed loop control approach which tries to minimise an error but fails to do so because it will always react too late. An adaptive algorithm can use this error to learn a forward model with the help of predictive cues. For example a driver learns to improve their steering by looking ahead to avoid steering in the last minute. In order to process complex cues such as the road ahead deep learning is a natural choice. However, this is usually only achieved indirectly by employing deep reinforcement learning having a discrete state space. Here, we show how this can be directly achieved by embedding deep learning into a closed loop system and preserving its continuous processing. We show specifically how error back-propagation can be achieved in z-space and in general how gradient based approaches can be analysed in such closed loop scenarios. The performance of this learning paradigm is demonstrated using a line-follower both in simulation and on a real robot that show very fast and continuous learning.


title: The Onset of 3D Magnetic Reconnection and Heating in the Solar Corona

abstract: Magnetic reconnection, a fundamentally important process in many aspects of astrophysics, is believed to be initiated by the tearing instability of an electric current sheet, a region where magnetic field abruptly changes direction and electric currents build up. Recent studies have suggested that the amount of magnetic shear in these structures is a critical parameter for the switch-on nature of magnetic reconnection in the solar atmosphere, at fluid spatial scales much larger than kinetic scales. We present results of simulations of reconnection in 3D current sheets with conditions appropriate to the solar corona. Using high-fidelity simulations, we follow the evolution of the linear and non-linear 3D tearing instability, leading to reconnection. We find that, depending on the parameter space, magnetic shear can play a vital role in the onset of significant energy release and heating via non-linear tearing. Two regimes in our study exist, dependent on whether the current sheet is longer or shorter than the wavelength of the fastest growing parallel mode (in the corresponding infinite system), thus determining whether sub-harmonics are present in the actual system. In one regime, where the fastest growing parallel mode has sub-harmonics, the non-linear interaction of these sub-harmonics and the coalescence of 3D plasmoids dominates the non-linear evolution, with magnetic shear playing only a weak role in the amount of energy released. In the second regime, where the fastest growing parallel mode has no-sub-harmonics, then only strongly sheared current sheets, where oblique mode are strong enough to compete with the dominant parallel mode, show any significant energy release. We expect both regimes to exist on the Sun, and so our results have important consequences for the the question of reconnection onset in different solar physics applications.


title: Flow-driven collapse of lubricant-infused surfaces

abstract: Lubricant-infused surfaces in an outer liquid flow generally reduce viscous drag. However, owing to the meniscus deformation, the infused state could collapse. Here we discuss the transition between infused and collapsed states of transverse shallow grooves, considering the capillary number, liquid/lubricant viscosity ratio and aspect ratio of the groove as parameters for inducing this transition. It is found that depending on the depth of the grooves, two different scenarios occur. A collapse of lubricant-infused surfaces could happen due to a depinning of the meniscus from the front groove edge. However, for very shallow textures, the meniscus contacts the bottom wall before such a depinning could occur. Our interpretation could help avoiding this generally detrimental effect in various applications.


title: Examining the secondary product origin of cosmic ray positrons with the latest AMS-02 data

abstract: Measurements of cosmic-ray (CR) positron fraction by PAMELA and other experiments have found an excess above 10 GeV relative to the standard predictions for secondary production in the interstellar medium (ISM). Although the excess has been mostly suggested to arise from some primary sources of positrons (such as pulsars and or annihilating dark matter particles), the almost constant flux ratio of $e^{+}/ \bar{p}$ argues for an alternative possibility that the excess positrons and antiprotons up to the highest energies are secondary products generated in hadronic interactions. Recently, Yang \& Aharonian (2019) revisit this possibility by assuming the presence of an additional population of CR nuclei sources. Here we examine this secondary product scenario using the \texttt{DRAGON} code, where the radiative loss of positrons is taken into account consistently. We confirm that the CR proton spectrum and the antiproton data can be explained by assuming the presence of an additional population of CR sources. However, the corresponding positron spectrum deviates from the measured data significantly above 100 GeV due to the strong radiative cooling. This suggests that, although hadronic interactions can explain the antiproton data, the corresponding secondary positron flux is still not enough to account for the AMS data. Hence contribution from some primary positron sources, such as pulsars or dark matter, is non-negligible.