title: The Non-Iterative Transformation Method

abstract: The Blasius flow is the idealized flow of a viscous fluid past an infinitesimally thick, semi-infinite flat plate. The definition of a non-iterative transformation method for the celebrated Blasius problem is due to T{\"o}pfer and dates more than a century ago. Here we define a non-iterative transformation method for Blasius equation with a moving wall, a slip flow condition or a surface gasification. The defined method allows us to deal with classes of problems in boundary layer theory that, depending on a parameter, admit multiple or no solutions. This approach is particularly convenient when the main interest is on the behaviour of the considered models with respect to the involved parameter. The obtained numerical results are found to be in good agreement with those available in literature.


title: Cooperative Problem Solving: an experience of high-school teaching updating

abstract: We present the results of an experience of teaching updating dispensed to Italian high-school physics teachers to promote the application of the Cooperative Problem Solving method as an useful strategy to improve physics learning at high-school level and to foster the development of problem solving skills. Beside analysing the method and discussing the ways to propose and apply it in a high-school context, the teachers experienced the method acting both as learners and as tutors of student group learners. Students and teachers evaluated as positive the experience, mainly focusing on cooperation within the group by information exchange and the application of a solution scheme. The ex-post analysis of the students' performance in applying the method to solve some rich context text showed the need of improving critical sense with respect to achieved results to fully exploit the strategy and develop their problem solving skills. Finally, an analysis on gender differences and scholar distribution of students is presented.


title: Total mass and limits of quasi-local mass at future null infinity

abstract: An idealized observer of an astronomical event is situated at future null infinity, where light rays emitted from the source approach. Mathematically, null infinity corresponds to the portion of the spacetime boundary defined by equivalence classes of null geodesics. But what can we observe at future null infinity? In the note, we start by reviewing the description of null infinity in terms of the Bondi-Sachs coordinate system. In particular, we reformulate and extend the invariance of Bondi mass and the mass loss formula in terms of a modified mass aspect 2-form. At the end, we discuss new results about the limit of quasilocal mass of unit spheres at null infinity in a Bondi-Sachs coordinate system. This is based on joint work with Po-Ning Chen, Jordan Keller, Ye-Kai Wang, and Shing-Tung Yau.


title: Learning Meta Face Recognition in Unseen Domains

abstract: Face recognition systems are usually faced with unseen domains in real-world applications and show unsatisfactory performance due to their poor generalization. For example, a well-trained model on webface data cannot deal with the ID vs. Spot task in surveillance scenario. In this paper, we aim to learn a generalized model that can directly handle new unseen domains without any model updating. To this end, we propose a novel face recognition method via meta-learning named Meta Face Recognition (MFR). MFR synthesizes the source/target domain shift with a meta-optimization objective, which requires the model to learn effective representations not only on synthesized source domains but also on synthesized target domains. Specifically, we build domain-shift batches through a domain-level sampling strategy and get back-propagated gradients/meta-gradients on synthesized source/target domains by optimizing multi-domain distributions. The gradients and meta-gradients are further combined to update the model to improve generalization. Besides, we propose two benchmarks for generalized face recognition evaluation. Experiments on our benchmarks validate the generalization of our method compared to several baselines and other state-of-the-arts. The proposed benchmarks will be available at https://github.com/cleardusk/MFR.


title: A Novel Online Action Detection Framework from Untrimmed Video Streams

abstract: Online temporal action localization from an untrimmed video stream is a challenging problem in computer vision. It is challenging because of i) in an untrimmed video stream, more than one action instance may appear, including background scenes, and ii) in online settings, only past and current information is available. Therefore, temporal priors, such as the average action duration of training data, which have been exploited by previous action detection methods, are not suitable for this task because of the high intra-class variation in human actions. We propose a novel online action detection framework that considers actions as a set of temporally ordered subclasses and leverages a future frame generation network to cope with the limited information issue associated with the problem outlined above. Additionally, we augment our data by varying the lengths of videos to allow the proposed method to learn about the high intra-class variation in human actions. We evaluate our method using two benchmark datasets, THUMOS'14 and ActivityNet, for an online temporal action localization scenario and demonstrate that the performance is comparable to state-of-the-art methods that have been proposed for offline settings.


title: Periodic solutions of a system of nonlinear difference equations with periodic coefficients

abstract: In this paper it is dealt with the following system of difference equations x_{n+1}=((a_{n})/(x_{n}))+((b_{n})/(y_{n})), y_{n+1}=((c_{n})/(x_{n}))+((d_{n})/(y_{n})), n in N_0, where the initial values x_0,y_0 are positive real numbers and the coefficients (a_{n})_{n>=0}, (b_{n})_{n>=0}, (c_{n})_{n>=0}, (d_{n})_{n>=0} are two-periodic sequences with positive terms. The system is an extention of a system that every positive solution is two-periodic or converges to its a two-periodic solution. Here, the long-term behavior of posistive solutions of the system is examined by using a new method to solve the system.


title: Anomalous supply shortages from dynamic pricing in on-demand mobility

abstract: Dynamic pricing schemes are increasingly employed across industries to maintain a self-organized balance of demand and supply. However, throughout complex dynamical systems, unintended collective states exist that may compromise their function. Here we reveal how dynamic pricing may induce demand-supply imbalances instead of preventing them. Combining game theory and time series analysis of dynamic pricing data from on-demand ride-hailing services, we explain this apparent contradiction. We derive a phase diagram demonstrating how and under which conditions dynamic pricing incentivizes collective action of ride-hailing drivers to induce anomalous supply shortages. By disentangling different timescales in price time series of ride-hailing services at 137 locations across the globe, we identify characteristic patterns in the price dynamics reflecting these anomalous supply shortages. Our results provide systemic insights for the regulation of dynamic pricing, in particular in publicly accessible mobility systems, by unraveling under which conditions dynamic pricing schemes promote anomalous supply shortages.


title: On H-sober spaces and H-sobrifications of $T_0$ spaces

abstract: In this paper, we provide a uniform approach to $d$-spaces, sober spaces and well-filtered spaces, and develop a general framework for dealing with all these spaces. For a subset system H, the theory of H-sober spaces and super H-sober spaces is established, and a direct construction of the H-sobrifcations and super H-sobrifications of $T_0$ spaces is given. Therefore, the category of all H-sober spaces is reflective in $\mathbf{Top}_0$, and the category of all super H-sober spaces is also reflective in $\mathbf{Top}_0$ if H has a natural property (called property M). It is shown that the H-sobrification preserves finite products of $T_0$ spaces, and the super H-sobrification preserves finite products of $T_0$ spaces if H has property M.


title: Longitudinal Dynamics Model Identification of an Electric Car Based on Real Response Approximation

abstract: Obtaining a realistic and accurate model of the longitudinal dynamics is key for a good speed control of a self-driving car. It is also useful to simulate the longitudinal behavior of the vehicle with high fidelity. In this paper, a straightforward and generic method for obtaining the friction, braking and propulsion forces as a function of speed, throttle input and brake input is proposed. Experimental data is recorded during tests over the full speed range to estimate the forces, to which the corresponding curves are adjusted. A simple and direct balance of forces in the direction tangent to the ground is used to obtain an estimation of the real forces involved. Then a model composed of approximate spline curves that fit the results is proposed. Using splines to model the dynamic response has the advantage of being quick and accurate, avoiding the complexity of parameter identification and tuning of non-linear responses embedding the internal functionalities of the car, like ABS or regenerative brake. This methodology has been applied to LS2N's electric Renault Zoe but can be applied to any other electric car. As shown in the experimental section, a comparison between the estimated acceleration of the car using the model and the real one over a wide range of speeds along a trip of about $10km/h$ reveals only $0.35m/s^2$ of error standard deviation in a range of $\pm{2}m/s^2$ which is very encouraging.


title: Formal Scenario-Based Testing of Autonomous Vehicles: From Simulation to the Real World

abstract: We present a new approach to automated scenario-based testing of the safety of autonomous vehicles, especially those using advanced artificial intelligence-based components, spanning both simulation-based evaluation as well as testing in the real world. Our approach is based on formal methods, combining formal specification of scenarios and safety properties, algorithmic test case generation using formal simulation, test case selection for track testing, executing test cases on the track, and analyzing the resulting data. Experiments with a real autonomous vehicle at an industrial testing facility support our hypotheses that (i) formal simulation can be effective at identifying test cases to run on the track, and (ii) the gap between simulated and real worlds can be systematically evaluated and bridged.


title: Geometric Approaches to Increase the Expressivity of Deep Neural Networks for MR Reconstruction

abstract: Recently, deep learning approaches have been extensively investigated to reconstruct images from accelerated magnetic resonance image (MRI) acquisition. Although these approaches provide significant performance gain compared to compressed sensing MRI (CS-MRI), it is not clear how to choose a suitable network architecture to balance the trade-off between network complexity and performance. Recently, it was shown that an encoder-decoder convolutional neural network (CNN) can be interpreted as a piecewise linear basis-like representation, whose specific representation is determined by the ReLU activation patterns for a given input image. Thus, the expressivity or the representation power is determined by the number of piecewise linear regions. As an extension of this geometric understanding, this paper proposes a systematic geometric approach using bootstrapping and subnetwork aggregation using an attention module to increase the expressivity of the underlying neural network. Our method can be implemented in both k-space domain and image domain that can be trained in an end-to-end manner. Experimental results show that the proposed schemes significantly improve reconstruction performance with negligible complexity increases.


title: Methodology and Performance of the Two-Year Galactic Plane Scanning Survey of Insight-HXMT

abstract: The Galactic plane scanning survey is one of the main scientific objectives of the Hard X-ray Modulation Telescope (known as Insight-HXMT). During the two-year operation of Insight-HXMT, more than 1000 scanning observations have been performed and the whole Galactic plane ($\rm 0^{\circ}<l<360^{\circ}$, $\rm -10^{\circ}<b<10^{\circ}$) has been covered completely. We summarize the Galactic plane scanning survey of Insight-HXMT for two years, including the characteristics of the scanning data, the data analysis process and the preliminary results of the Low-Energy telescope, the Medium-Energy telescope and the High-Energy telescope. With the light curve PSF fitting method, the fluxes of the known sources in the scanned area as well as the flux errors are obtained for each scanning observation. From the relationships of SNRs and fluxes, the $5\sigma$ sensitivities of three telescopes of Insight-HXMT are estimated as $\rm \sim7.6\times10^{-11}~erg cm^{-2}~s^{-1}$ ($\rm 3 mCrab,~1-6 keV$), $\rm \sim4.0\times10^{-10}~erg~cm^{-2}~s^{-1}$ ($\rm 20~mCrab,~7-40 keV$) and $\rm \sim2.6\times10^{-10}~erg cm^{-2}~s^{-1}$ ($\rm 18 mCrab,~25-100 keV$) for an individual scanning observation of $2-3$ hours, respectively. Up to September 2019, more than 800 X-ray sources with various types are monitored by the three telescopes and their long-term light curves with three energy bands are obtained to make further scientific analyses.


title: On a 'philosophical' question about Banach envelopes

abstract: We show how to construct nonlocally convex quasi-Banach spaces $X$ whose dual separates the points of a dense subspace of $X$ but does not separate the points of $X$. Our examples connect with a question raised by Pietsch [About the Banach envelope of $l_{1,\infty}$, Rev. Mat. Complut. 22 (2009), 209-226] and shed light into the unexplored class of quasi-Banach spaces with nontrivial dual which do not have sufficiently many functionals to separate the points of the space.


title: A Benchmarking Study of Embedding-based Entity Alignment for Knowledge Graphs

abstract: Entity alignment seeks to find entities in different knowledge graphs (KGs) that refer to the same real-world object. Recent advancement in KG embedding impels the advent of embedding-based entity alignment, which encodes entities in a continuous embedding space and measures entity similarities based on the learned embeddings. In this paper, we conduct a comprehensive experimental study of this emerging field. We survey 23 recent embedding-based entity alignment approaches and categorize them based on their techniques and characteristics. We also propose a new KG sampling algorithm, with which we generate a set of dedicated benchmark datasets with various heterogeneity and distributions for a realistic evaluation. We develop an open-source library including 12 representative embedding-based entity alignment approaches, and extensively evaluate these approaches, to understand their strengths and limitations. Additionally, for several directions that have not been explored in current approaches, we perform exploratory experiments and report our preliminary findings for future studies. The benchmark datasets, open-source library and experimental results are all accessible online and will be duly maintained.


title: Homotopy properties of the space $OS_{f}(X)$

abstract: For a given compact Hausdorff space $X$, we construct the space $OS_{f}(X)$ of normed, order-preserving, weakly additive, positively homogeneous and semi-additive functionals (for brevity, semi-additive functionals) and it is proved that the hyperspace $\exp\, X$ of the space $X$ is a deformation retract of the constructed space. Further we show that the shapes of the spaces $OS_{f}(X)$ and $\exp\, X$ coincide. We establish that if the space $\exp\, X$ is contractible, then the space $OS_{f}(X)$ is also contractible.


title: Learning to Optimize Autonomy in Competence-Aware Systems

abstract: Interest in semi-autonomous systems (SAS) is growing rapidly as a paradigm to deploy autonomous systems in domains that require occasional reliance on humans. This paradigm allows service robots or autonomous vehicles to operate at varying levels of autonomy and offer safety in situations that require human judgment. We propose an introspective model of autonomy that is learned and updated online through experience and dictates the extent to which the agent can act autonomously in any given situation. We define a competence-aware system (CAS) that explicitly models its own proficiency at different levels of autonomy and the available human feedback. A CAS learns to adjust its level of autonomy based on experience to maximize overall efficiency, factoring in the cost of human assistance. We analyze the convergence properties of CAS and provide experimental results for robot delivery and autonomous driving domains that demonstrate the benefits of the approach.


title: NP-Completeness Results for Graph Burning on Geometric Graphs

abstract: Graph burning runs on discrete time steps. The aim is to burn all the vertices in a given graph in the least number of time steps. This number is known to be the burning number of the graph. The spread of social influence, an alarm, or a social contagion can be modeled using graph burning. The less the burning number, the faster the spread. Optimal burning of general graphs is NP-Hard. There is a 3-approximation algorithm to burn general graphs where as better approximation factors are there for many sub classes. Here we study burning of grids; provide a lower bound for burning arbitrary grids and a 2-approximation algorithm for burning square grids. On the other hand, burning path forests, spider graphs, and trees with maximum degree three is already known to be NP-Complete. In this article we show burning problem to be NP-Complete on connected interval graphs, permutation graphs and several other geometric graph classes as corollaries.


title: A laser-plasma interaction experiment for solar burst studies

abstract: A new experimental platform based on laser-plasma interaction is proposed to explore the fundamental processes of wave coupling at the origin of interplanetary radio emissions. It is applied to the study of electromagnetic (EM) emission at twice the plasma frequency ($2\omega_p$) observed during solar bursts and thought to result from the coalescence of two Langmuir waves (LWs). In the interplanetary medium, the first LW is excited by electron beams, while the second is generated by electrostatic decay of Langmuir waves. In the present experiment, instead of an electron beam, an energetic laser propagating through a plasma excites the primary LW, with characteristics close to those at near-Earth orbit. The EM radiation at $2\omega_p$ is observed at different angles. Its intensity, spectral evolution and polarization confirm the LW-coalescence scenario.


title: NSBchain: A Secure Blockchain Framework for Network Slicing Brokerage

abstract: With the advent of revolutionary technologies, such as virtualization and softwarization, a novel concept for 5G networks and beyond has been unveiled: Network Slicing. Initially driven by the research community, standardization bodies as 3GPP have embraced it as a promising solution to revolutionize the traditional mobile telecommunication market by enabling new business models opportunities. Network Slicing is envisioned to open up the telecom market to new players such as Industry Verticals, e.g. automotive, smart factories, e-health, etc. Given the large number of potential new business players, dubbed as network tenants, novel solutions are required to accommodate their needs in a cost-efficient and secure manner. In this paper, we propose NSBchain, a novel network slicing brokering (NSB) solution, which leverages on the widely adopted Blockchain technology to address the new business models needs beyond traditional network sharing agreements. NSBchain defines a new entity, the Intermediate Broker (IB), which enables Infrastructure Providers (InPs) to allocate network resources to IBs through smart contracts and IBs to assign and re-distribute their resources among tenants in a secure, automated and scalable manner. We conducted an extensive performance evaluation by means of an open-source blockchain platform that proves the feasibility of our proposed framework considering a large number of tenants and two different consensus algorithms.


title: Uniform rectifiability implies Varopoulos extensions

abstract: We construct extensions of Varopolous type for functions $f \in \text{BMO}(E)$, for any uniformly rectifiable set $E$ of codimension one. More precisely, let $\Omega \subset \mathbb{R}^{n+1}$ be an open set satisfying the corkscrew condition, with an $n$-dimensional uniformly rectifiable boundary $\partial \Omega$, and let $\sigma := \mathcal{H}^n\lfloor_{\partial \Omega}$ denote the surface measure on $\partial \Omega$. We show that if $f \in \text{BMO}(\partial \Omega,d\sigma)$ with compact support on $\partial \Omega$, then there exists a smooth function $V$ in $\Omega$ such that $|\nabla V(Y)| \, dY$ is a Carleson measure with Carleson norm controlled by the BMO norm of $f$, and such that $V$ converges in some non-tangential sense to $f$ almost everywhere with respect to $\sigma$. Our results should be compared to recent geometric characterizations of $L^p$-solvability and of BMO-solvability of the Dirichlet problem, by Azzam, the first author, Martell, Mourgoglou and Tolsa and by the first author and Le, respectively. In combination, this latter pair of results shows that one can construct, for all $f \in C_c(\partial \Omega)$, a harmonic extension $u$, with $|\nabla u(Y)|^2 \text{dist}(Y,\partial \Omega) \, dY $ a Carleson measure controlled by the BMO norm of $f$, only in the presence of an appropriate quantitative connectivity condition.


title: Electrical generation and detection of terahertz signal based on spin-wave emission from ferrimagnets

abstract: Terahertz (THz) signals, mainly generated by photonic or electronic approaches, are being sought for various applications, whereas the development of magnetic source might be a necessary step to harness the magnetic nature of electromagnetic radiation. We show that the relativistic effect on the current-driven domain-wall motion induces THz spin-wave emission in ferrimagnets. The required current density increases dramatically in materials with strong exchange interaction and rapidly exceeds 1012 A m-2, leading to the device breakdown and thus the lack of experimental evidence. By translating the collective magnetization oscillations into voltage signals, we propose a three-terminal device for the electrical detection of THz spin wave. Through material engineering, wide frequency range from 264 GHz to 1.1 THz and uniform continuous signals with improved output power can be obtained. As a reverse effect, the spin wave generated in this system is able to move ferrimagnetic domain wall. Our work provides guidelines for the experimental verification of THz spin wave, and could stimulate the design of THz spintronic oscillators for wideband applications as well as the all-magnon spintronic devices.


title: Around a theorem of F. Dyson and A. Lenard: Energy Equilibria for Point Charge Distributions in Classical Electrostatics

abstract: We discuss several results in electrostatics: Onsager's inequality, an extension of Earnshaw's theorem, and a result stemming from the celebrated conjecture of Maxwell on the number of points of electrostatic equilibrium. Whenever possible, we try to provide a brief historical context and references.

categories: math-ph math.CV math.MP physics.class-ph


title: Automated calculation and convergence of defect transport tensors

abstract: Defect transport is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewise sampling completeness, to build atomistic kinetic Monte Carlo models on a state space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for defect transport tensors and provide a convergence metric by calculating the Kullback-Leiber divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

categories: physics.comp-ph cond-mat.mtrl-sci cond-mat.stat-mech physics.data-an


title: Adversarial Reverse Mapping of Equilibrated Condensed-Phase Molecular Structures

abstract: A tight and consistent link between resolutions is crucial to further expand the impact of multiscale modeling for complex materials. We herein tackle the generation of condensed molecular structures as a refinement -- backmapping -- of a coarse-grained structure. Traditional schemes start from a rough coarse-to-fine mapping and perform further energy minimization and molecular dynamics simulations to equilibrate the system. In this study we introduce DeepBackmap: A deep neural network based approach to directly predict equilibrated molecular structures for condensed-phase systems. We use generative adversarial networks to learn the Boltzmann distribution from training data and realize reverse mapping by using the coarse-grained structure as a conditional input. We apply our method to a challenging condensed-phase polymeric system. We observe that the model trained in a melt has remarkable transferability to the crystalline phase. The combination of data-driven and physics-based aspects of our architecture help reach temperature transferability with only limited training data.


title: Robust Task and Motion Planning for Long-Horizon Architectural Construction Planning

abstract: Integrating robotic systems in architectural and construction processes is of core interest to increase the efficiency of the building industry. Automated planning for such systems enables design analysis tools and facilitates faster design iteration cycles for designers and engineers. However, generic task-and-motion planning (TAMP) for long-horizon construction processes is beyond the capabilities of current approaches. In this paper, we develop a multi-agent TAMP framework for long horizon problems such as constructing a full-scale building. To this end we extend the Logic-Geometric Programming framework by sampling-based motion planning,a limited horizon approach, and a task-specific structural stability optimization that allow an effective decomposition of the task. We show that our framework is capable of constructing a large pavilion built from several hundred geometrically unique building elements from start to end autonomously.


title: Tug-of-War in a Double-Nanopore System

abstract: We simulate a tug-of-war (TOW) scenario for a model double-stranded DNA threading through a double nanopore (DNP) system. The DNA, simultaneously captured at both pores is subject to two equal and opposite forces $-\vec{f}_L= \vec{f}_R$ (TOW), where $\vec{f}_L$ and $\vec{f}_R$ are the forces applied to the left and the right pore respectively. Even though the net force on the DNA polymer $\Delta \vec{f}_{LR}=\vec{f}_L+ \vec{f}_R=0$, the mean first passage time (MFPT) $\langle \tau \rangle$ depends on the magnitude of the TOW forces $ \left | f_L \right | = \left |f_R \right | = f_{LR}$. We qualitatively explain this dependence of $\langle \tau \rangle$ on $f_{LR}$ from the known results for the single-pore translocation of a triblock copolymer. We demonstrate that the time of flight (TOF) of a monomer with index $m$ ($\langle \tau_{LR}(m) \rangle$) from one pore to the other exhibits quasi-periodic structure commensurate with the distance between the pores $d_{LR}$. Finally, we study the case $\Delta \vec{f}_{LR}=\vec{f}_L+ \vec{f}_R \ne 0$, and qualitatively reproduce the experimental result of the dependence of the MFPT on $\Delta\vec{f}_{LR}$. For a moderate bias, the MFPT for the DNP system for a chain length $N$ follows the same scaling ansatz as that of for the single nanopore, $\langle \tau \rangle = \left( AN^{1+\nu} + \eta_{pore}N \right) \left(\Delta f_{LR}\right)^{-1}$, where $\eta_{pore}$ is the pore friction, which enables us to estimate $\langle \tau \rangle $ for a long chain. Our Brownian dynamics simulation studies provide fundamental insights and valuable information about the details of the translocation speed obtained from $\langle \tau_{LR}(m) \rangle$, and accuracy of the translation of the data obtained in the time-domain to units of genomic distances.


title: Characterizing and Avoiding Problematic Global Optima of Variational Autoencoders

abstract: Variational Auto-encoders (VAEs) are deep generative latent variable models consisting of two components: a generative model that captures a data distribution p(x) by transforming a distribution p(z) over latent space, and an inference model that infers likely latent codes for each data point (Kingma and Welling, 2013). Recent work shows that traditional training methods tend to yield solutions that violate modeling desiderata: (1) the learned generative model captures the observed data distribution but does so while ignoring the latent codes, resulting in codes that do not represent the data (e.g. van den Oord et al. (2017); Kim et al. (2018)); (2) the aggregate of the learned latent codes does not match the prior p(z). This mismatch means that the learned generative model will be unable to generate realistic data with samples from p(z)(e.g. Makhzani et al. (2015); Tomczak and Welling (2017)). In this paper, we demonstrate that both issues stem from the fact that the global optima of the VAE training objective often correspond to undesirable solutions. Our analysis builds on two observations: (1) the generative model is unidentifiable - there exist many generative models that explain the data equally well, each with different (and potentially unwanted) properties and (2) bias in the VAE objective - the VAE objective may prefer generative models that explain the data poorly but have posteriors that are easy to approximate. We present a novel inference method, LiBI, mitigating the problems identified in our analysis. On synthetic datasets, we show that LiBI can learn generative models that capture the data distribution and inference models that better satisfy modeling assumptions when traditional methods struggle to do so.


title: Black phosphorus van der Waals heterostructures light emitting diodes for mid-infrared silicon photonics

abstract: Light-emitting diodes (LEDs) based on III-V/II-VI materials have delivered a compelling performance in the mid-infrared (mid-IR) region, which enabled wide-ranging applications, including environmental monitoring, defense and medical diagnostics. Continued efforts are underway to realize on-chip sensors via heterogeneous integration of mid-IR emitters on a silicon photonic chip. But the uptake of such approach is limited by the high costs and interfacial strains, associated with the process of heterogeneous integrations. Here, the black phosphorus (BP)-based van der Waals (vdW) heterostructures are exploited as room temperature LEDs. The demonstrated devices can emit linearly polarized light, and their spectra cover the technologically important mid-IR atmospheric window (3-4 um). Additionally, the BP LEDs exhibit fast modulation speed as well as exceptional stability, and its peak extrinsic quantum efficiency (QE~0.9%) is comparable to the III-V/II-VI mid-IR LEDs. By leveraging the integrability of vdW heterostructures, we further demonstrate a silicon photonic waveguide-integrated BP LED. The reported hybrid platform holds great promise for mid-IR silicon photonics.


title: Multi-modal Dense Video Captioning

abstract: Dense video captioning is a task of localizing interesting events from an untrimmed video and producing textual description (captions) for each localized event. Most of the previous works in dense video captioning are solely based on visual information and completely ignore the audio track. However, audio, and speech, in particular, are vital cues for a human observer in understanding an environment. In this paper, we present a new dense video captioning approach that is able to utilize any number of modalities for event description. Specifically, we show how audio and speech modalities may improve a dense video captioning model. We apply automatic speech recognition (ASR) system to obtain a temporally aligned textual description of the speech (similar to subtitles) and treat it as a separate input alongside video frames and the corresponding audio track. We formulate the captioning task as a machine translation problem and utilize recently proposed Transformer architecture to convert multi-modal input data into textual descriptions. We demonstrate the performance of our model on ActivityNet Captions dataset. The ablation studies indicate a considerable contribution from audio and speech components suggesting that these modalities contain substantial complementary information to video frames. Furthermore, we provide an in-depth analysis of the ActivityNet Caption results by leveraging the category tags obtained from original YouTube videos. Code is publicly available: github.com/v-iashin/MDVC


title: Dalitz analysis of $D^{0}\to K^{-}\pi^{+}\eta$ decays at Belle

abstract: We present the results of the first Dalitz plot analysis of the decay $D^{0}\to K^{-}\pi^{+}\eta$. The analysis is performed on a data set corresponding to an integrated luminosity of 953 $\rm{fb}^{-1}$ collected by the Belle detector at the asymmetric-energy $e^{+}e^{-}$ KEKB collider. The Dalitz plot is well described by a combination of the six resonant decay channels $\bar{K}^{*}(892)^0\eta$, $K^{-}a_0(980)^+$, $K^{-}a_2(1320)^+$, $\bar{K}^{*}(1410)^0\eta$, $K^{*}(1680)^-\pi^{+}$ and $K_2^{*}(1980)^-\pi^{+}$, together with $K\pi$ and $K\eta$ S-wave components. The decays $K^{*}(1680)^{-}\to K^{-}\eta$ and $K_{2}^{*}(1980)^{-}\to K^{-}\eta$ are observed for the first time. We measure ratio of the branching fractions, $\frac{\mathcal{B}(D^{0}\to K^{-}\pi^{+}\eta)}{\mathcal{B}(D^{0}\to K^{-}\pi^{+})}=0.500\pm0.002{\rm(stat)}\pm0.020{\rm(syst)}\pm0.003{\rm (\mathcal{B}_{PDG})}$. Using the Dalitz fit result, the ratio $\frac{\mathcal{B}(K^{*}(1680)\to K\eta)}{\mathcal{B}(K^{*}(1680)\to K\pi)}$ is measured to be $0.11\pm0.02{\rm(stat)}^{+0.06}_{-0.04}{\rm(syst)}\pm0.04{\rm(\mathcal{B}_{\text{PDG}})}$; this is much lower than the theoretical expectations ($\approx1$) made under the assumption that $K^{*}(1680)$ is a pure $1^{3}D_1$ state. The product branching fraction $\mathcal{B}(D^0\to [K_2^{*}(1980)^-\to K^{-}\eta]\pi^{+})=(2.2^{+1.7}_{-1.9})\times10^{-4}$ is determined. In addition, the $\pi\eta^{\prime}$ contribution to the $a_0(980)^{\pm}$ resonance shape is confirmed with 10.1$\sigma$ statistical significance using the three-channel Flatt\'{e} model. We also measure $\mathcal{B}(D^0\to\bar{K}^{*}(892)^0\eta)=(1.41^{+0.13}_{-0.12})\%$. This is consistent with, and more precise than, the current world average $(1.02\pm0.30)\%$, deviates with a significance of more than $3\sigma$ from the theoretical predictions of (0.51-0.92)%.


title: Scaling Strongly Consistent Replication

abstract: Strong consistency replication helps keep application logic simple and provides significant benefits for correctness and manageability. Unfortunately, the adoption of strongly-consistent replication protocols has been curbed due to their limited scalability and performance. To alleviate the leader bottleneck in strongly-consistent replication protocols, we introduce Pig, an in-protocol communication aggregation and piggybacking technique. Pig employs randomly selected nodes from follower subgroups to relay the leader's message to the rest of the followers in the subgroup, and to perform in-network aggregation of acknowledgments back from these followers. By randomly alternating the relay nodes across replication operations, Pig shields the relay nodes as well as the leader from becoming hotspots and improves throughput scalability. We showcase Pig in the context of classical Paxos protocols employed for strongly consistent replication by many cloud computing services and databases. We implement and evaluate PigPaxos, in comparison to Paxos and EPaxos protocols under various workloads over clusters of size 5 to 25 nodes. We show that the aggregation at the relay has little latency overhead, and PigPaxos can provide more than 3 folds improved throughput over Paxos and EPaxos with little latency deterioration. We support our experimental observations with the analytical modeling of the bottlenecks and show that the rotating of the relay nodes provides the most benefit for reducing the bottlenecks and that the throughput is maximized when employing only 1 randomly rotating relay node.


title: CycleISP: Real Image Restoration via Improved Data Synthesis

abstract: The availability of large-scale datasets has helped unleash the true potential of deep convolutional neural networks (CNNs). However, for the single-image denoising problem, capturing a real dataset is an unacceptably expensive and cumbersome procedure. Consequently, image denoising algorithms are mostly developed and evaluated on synthetic data that is usually generated with a widespread assumption of additive white Gaussian noise (AWGN). While the CNNs achieve impressive results on these synthetic datasets, they do not perform well when applied on real camera images, as reported in recent benchmark datasets. This is mainly because the AWGN is not adequate for modeling the real camera noise which is signal-dependent and heavily transformed by the camera imaging pipeline. In this paper, we present a framework that models camera imaging pipeline in forward and reverse directions. It allows us to produce any number of realistic image pairs for denoising both in RAW and sRGB spaces. By training a new image denoising network on realistic synthetic data, we achieve the state-of-the-art performance on real camera benchmark datasets. The parameters in our model are ~5 times lesser than the previous best method for RAW denoising. Furthermore, we demonstrate that the proposed framework generalizes beyond image denoising problem e.g., for color matching in stereoscopic cinema. The source code and pre-trained models are available at https://github.com/swz30/CycleISP.


title: The Jang equation and the positive mass theorem in the asymptotically hyperbolic setting

abstract: We solve the Jang equation with respect to asymptotically hyperbolic "hyperboloidal" initial data. The results are applied to give a non-spinor proof of the positive mass theorem in the asymptotically hyperbolic setting. This work focuses on the case when the spatial dimension is equal to three.


title: Synchronous whirling of spinning homogeneous elastic cylinders: linear and weakly non-linear analyses

abstract: Stationary whirling of slender and homogeneous (continuous) elastic shafts rotating around their axis, with pin-pin boundary condition at the ends, is revisited by considering the complete deformations in the cross section of the shaft. The stability against a synchronous sinusoidal disturbance of any wave length is investigated and the analytic expression of the buckling amplitude is derived in the weakly non-linear regime by considering both geometric and material (hyper-elastic) non-linearities. The bifurcation is super-critical in the long wave length domain for any elastic constitutive law, and sub-critical in the short wave length limit for a limited range of non-linear material parameters.


title: Polynomial Superlevel Set Representation of the Multistationarity Region of Chemical Reaction Networks

abstract: In this paper we introduce a new representation for the multistationarity region of a reaction network, using polynomial superlevel sets. The advantages of using this polynomial superlevel set representation over the already existing representations (cylindrical algebraic decompositions, numeric sampling, rectangular divisions) is discussed, and algorithms to compute this new representation are provided. The results are given for the general mathematical formalism of a parametric system of equations and so may be applied to other application domains.


title: Kink Instabilities In Relativistic Jets Can Drive Quasi-Periodic Radiation Signatures

abstract: Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate nonthermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have studied the conditions required for the growth of kink instabilities in relativistic jets, the radiation signatures of these instabilities have not been investigated in detail. In this paper, we aim to self-consistently study radiation and polarization signatures from kink instabilities in relativistic jets. We combine large-scale relativistic magnetohydrodynamic (RMHD) simulations with polarized radiation transfer of a magnetized jet, which emerges from the central engine and propagates through the surrounding medium. We observe that a localized region at the central spine of the jet exhibits the strongest kink instabilities, which we identify as the jet emission region. Very interestingly, we find quasi-periodic oscillation (QPO) signatures in the light curve from the emission region. Additionally, the polarization degree appears to be anti-correlated to flares in the light curves. Our analyses show that these QPO signatures are intrinsically driven by kink instabilities, where the period of the QPOs is associated to the kink growth time scale. The latter corresponds to weeks to months QPOs in blazars. The polarization signatures offer unique diagnostics for QPOs driven by kink instabilities.


title: Flag-dipole spinors: On the dual structure derivation and $\mathcal{C}$, $\mathcal{P}$ and $\mathcal{T}$ symmetries

abstract: In this manuscript we report the flag-dipole spinors dual structure direct definition and analyze the properties behind the corresponding operator which generates such structure. This particular construction may be interesting for cosmological, phenomenological and mathematical physics applications. In addition, we analyse the behaviour of the flag-dipole spinors under action of discrete symmetries, facing an \emph{unconventional} property encoded on $(\mathcal{CPT})^2$.


title: Effects of deterministic disorder at deeply subwavelength scales in multilayered dielectric metamaterials

abstract: It is common understanding that multilayered dielectric metamaterials, in the regime of deeply subwavelength layers, are accurately described by simple effective-medium models based on mixing formulas that do not depend on the spatial arrangement. In the wake of recent studies that have shown counterintuitive examples of periodic and aperiodic (orderly or random) scenarios in which this premise breaks down, we study here the effects of deterministic disorder. With specific reference to a model based on Golay-Rudin-Shapiro sequences, we illustrate certain peculiar boundary effects that can occur in finite-size dielectric multilayers, leading to anomalous light-transport properties that are in stark contrast with the predictions from conventional effective-medium theory. Via parametric and comparative studies, we elucidate the underlying physical mechanisms, also highlighting similarities and differences with respect to previously studied geometries. Our outcomes may inspire potential applications to optical sensing, switching and lasing.


title: Predictions for charmed nuclei based on $Y_c N$ forces inferred from lattice QCD simulations

abstract: Charmed nuclei are investigated utilizing $\Lambda_c N$ and $\Sigma_c N$ interactions that have been extrapolated from lattice QCD simulations at unphysical masses of $m_\pi = 410$--$570$ MeV to the physical point using chiral effective field theory as guideline. Calculations of the energies of $\Lambda_c$ single-particle bound states for various charmed nuclei from $^{\ 5}_{\Lambda_c}$Li to $^{209}_{\Lambda_c}$Bi are performed using a perturbative many-body approach. This approach allows one to determine the finite nuclei $\Lambda_c$ self-energy from which the energies of the different bound states can be obtained. Though the $\Lambda_c N$ interaction inferred from the lattice results is only moderately attractive, it supports the existence of charmed nuclei. Already the lightest nucleus considered is found to be bound. The spin-orbit splitting of the p- and d-wave states turns out to be small, as in the case of single $\Lambda$ hypernuclei. Additional calculations based on the Faddeev-Yakubovsky equations suggest that also $A=4$ systems involving a $\Lambda_c$ baryon are likely to be bound, but exclude a bound $^{\, 3}_{\Lambda_c}$He state.


title: Testing the nature of Gauss-Bonnet gravity by four-dimensional rotating black hole shadow

abstract: The recent discovery of the novel four-dimensional static and spherically symmetric Gauss-Bonnet black hole provides a promising bed to test Gauss-Bonnet gravity by using astronomical observations [Phys. Rev. Lett. 124, 081301 (2020)]. In this paper, we first obtain the rotating Gauss-Bonnet black hole solution by using the Newman-Janis algorithm, and then study the shadow cast by the nonrotating and rotating candidate Gauss-Bonnet black holes. The result indicates that positive metric parameter $\alpha$ shrinks the shadow, while negative one enlarges it. Meanwhile, both the distortion and ratio of two diameters of the shadow are found to increase with the metric parameter for certain spin. Comparing with the Kerr black hole, the shadow gets more distorted for $\alpha$, and less distorted for negative one. Furthermore, we calculate the angular diameter of the shadow by making use of the observation of M87*. The result indicates that negative metric parameter $\alpha$ in (-4.5, 0) is more favored. Since the negative energy appears for negative $\alpha$, our results extends the study of Gauss-Bonnet gravity. We believe further study on the four-dimensional rotating black hole may shed new light on Gauss-Bonnet gravity.


title: Hierarchical Models: Intrinsic Separability in High Dimensions

abstract: It has long been noticed that high dimension data exhibits strange patterns. This has been variously interpreted as either a "blessing" or a "curse", causing uncomfortable inconsistencies in the literature. We propose that these patterns arise from an intrinsically hierarchical generative process. Modeling the process creates a web of constraints that reconcile many different theories and results. The model also implies high dimensional data posses an innate separability that can be exploited for machine learning. We demonstrate how this permits the open-set learning problem to be defined mathematically, leading to qualitative and quantitative improvements in performance.


title: From flat to tilted: gradual interfaces in organic thin film growth

abstract: We investigate domain formation and local morphology of thin films of $\alpha$-sexithiophene ($\alpha$-6T) on Au(100) beyond monolayer coverage by combining high resolution scanning tunneling microscopy (STM) experiments with electronic structure theory calculations and computational structure search. We report a layerwise growth of highly-ordered enantiopure domains. For the second and third layer, we show that the molecular orbitals of individual $\alpha$-6T molecules can be well resolved by STM, providing access to detailed information on the molecular orientation. We find that already in the second layer the molecules abandon the flat adsorption structure of the monolayer and adopt a tilted conformation. Although the observed tilted arrangement resembles the orientation of $\alpha$-6T in the bulk, the observed morphology does not yet correspond to a well-defined surface of the $\alpha$-6T bulk structure. A similar behavior is found for the third layer indicating a growth mechanism where the bulk structure is gradually adopted over several layers.


title: On maps which preserve semipositivity and quantifier elimination theory for real numbers

abstract: Assume that $\Phi:\mathbb{M}_{n}(\mathbb{C})\rightarrow\mathbb{M}_{n}(\mathbb{C})$ is a superoperator which preserves hermiticity. We give an algorithm determining whether $\Phi$ preserves semipositivity (we call $\Phi$ positive in this case). Our approach to the problem has a model-theoretic nature, namely, we apply techniques of quantifier elimination theory for real numbers. An approach based on these techniques seems to be the only one that allows to decide whether an arbitrary hermiticity-preserving $\Phi$ is positive. Before we go to detailed analysis of the problem, we argue that quantifier elimination for real numbers (and also for complex numbers) can play a significant role in quantum information theory and other areas as well.


title: NDE 4.0 From Design Thinking to Strategy

abstract: Cyber technologies are offering new horizons for quality control in manufacturing and safety assurance in-service of physical assets. The line between non-destructive evaluation (NDE) and Industry 4.0 is getting blurred since both are sensory data-driven domains. This multidisciplinary approach has led to the emergence of a new capability: NDE 4.0. The NDT community is coming together once again to define the purpose, chart the process, and address the adoption of emerging technologies. In this paper, the authors have taken a design thinking approach to spotlight proper objectives for research on this subject. It begins with qualitative research on twenty different perceptions of stakeholders and misconceptions around the current state of NDE. The interpretation is used to define ten value propositions or use cases under "NDE for Industry 4.0" and "Industry 4.0 for NDE" leading up to the clarity of purpose for NDE 4.0, enhanced safety and economic value for stakeholders. To pursue this worthy cause, the paper delves into some of the top adoption challenges, and proposes a journey of managed innovation, conscious skills development, and a new form of leadership required to succeed in the cyber-physical world.


title: Decomposing Jacobians via Galois covers

abstract: Let $\phi:\,X\rightarrow Y$ be a (possibly ramified) cover between two algebraic curves of positive genus. We develop tools that may identify the Prym variety of $\phi$, up to isogeny, as the Jacobian of a quotient curve $C$ in the Galois closure of the composition of $\phi$ with a well-chosen map $Y\rightarrow \mathbb{P}^1$. This method allows us to recover all previously obtained descriptions of a Prym variety in terms of a Jacobian that are known to us, besides yielding new applications. We also find algebraic equations for some of these new cases, including one where $X$ has genus $3$, $Y$ has genus $1$ and $\phi$ is a degree $3$ map totally ramified over $2$ points.


title: Deep generative models in DataSHIELD

abstract: The best way to calculate statistics from medical data is to use the data of individual patients. In some settings, this data is difficult to obtain due to privacy restrictions. In Germany, for example, it is not possible to pool routine data from different hospitals for research purposes without the consent of the patients. The DataSHIELD software provides an infrastructure and a set of statistical methods for joint analyses of distributed data. The contained algorithms are reformulated to work with aggregated data from the participating sites instead of the individual data. If a desired algorithm is not implemented in DataSHIELD or cannot be reformulated in such a way, using artificial data is an alternative. We present a methodology together with a software implementation that builds on DataSHIELD to create artificial data that preserve complex patterns from distributed individual patient data. Such data sets of artificial patients, which are not linked to real patients, can then be used for joint analyses. We use deep Boltzmann machines (DBMs) as generative models for capturing the distribution of data. For the implementation, we employ the package "BoltzmannMachines" from the Julia programming language and wrap it for use with DataSHIELD, which is based on R. As an exemplary application, we conduct a distributed analysis with DBMs on a synthetic data set, which simulates genetic variant data. Patterns from the original data can be recovered in the artificial data using hierarchical clustering of the virtual patients, demonstrating the feasibility of the approach. Our implementation adds to DataSHIELD the ability to generate artificial data that can be used for various analyses, e. g. for pattern recognition with deep learning. This also demonstrates more generally how DataSHIELD can be flexibly extended with advanced algorithms from languages other than R.


title: Precise deviations for disk counting statistics of invariant determinantal processes

abstract: We consider two-dimensional determinantal processes which are rotation-invariant and study the fluctuations of the number of points in disks. Based on the theory of mod-phi convergence, we obtain Berry-Esseen as well as precise moderate to large deviation estimates for these statistics. These results are consistent with the Coulomb gas heuristic from the physics literature. We also obtain functional limit theorems for the stochastic process $(\# D_r)_{r>0}$ when the radius $r$ of the disk $D_r$ is growing in different regimes. We present several applications to invariant determinantal process, including the polyanalytic Ginibre ensembles, Gaussian analytic function and other hyperbolic models. As a corollary, we compute the precise asymptotics for the entanglement entropy of (integer) Laughlin states for all Landau levels.


title: On a population model in discrete periodic habitat. I. spreading speed and optimal dispersal strategy

abstract: Starting from an age-structured diffusive population growth law for single species in a discrete and periodic habitat, we formulate a stage structured population model with spatially periodic dispersal, mortality and recruitment. With a KPP type setting, after establishing the fundamental solution of a discretized heat equation with spatially periodic dispersal, we apply some recently developed dynamical system theories to obtain the existence of the spreading speed and its coincidence with the minimal speed of pulsating waves, as well as the variational characterization of the speed, by which we further analyze how the habitat periodicity influences the speed. In particular, there is a unique optimal dispersal strategy to maximize the speed.


title: Day Level Forecasting for Coronavirus Disease (COVID-19) Spread: Analysis, Modeling and Recommendations

abstract: In mid of March 2020, Coronaviruses such as COVID-19 is declared as an international epidemic. More than 125000 confirmed cases and 4,607 death cases have been recorded around more than 118 countries. Unfortunately, a coronavirus vaccine is expected to take at least 18 months if it works at all. Moreover, COVID -19 epidemics can mutate into a more aggressive form. Day level information about the COVID -19 spread is crucial to measure the behavior of this new virus globally. Therefore, this study presents a comparison of day level forecasting models on COVID-19 affected cases using time series models and mathematical formulation. The forecasting models and data strongly suggest that the number of coronavirus cases grows exponentially in countries that do not mandate quarantines, restrictions on travel and public gatherings, and closing of schools, universities, and workplaces (Social Distancing).


title: Unified Multi-Domain Learning and Data Imputation using Adversarial Autoencoder

abstract: We present a novel framework that can combine multi-domain learning (MDL), data imputation (DI) and multi-task learning (MTL) to improve performance for classification and regression tasks in different domains. The core of our method is an adversarial autoencoder that can: (1) learn to produce domain-invariant embeddings to reduce the difference between domains; (2) learn the data distribution for each domain and correctly perform data imputation on missing data. For MDL, we use the Maximum Mean Discrepancy (MMD) measure to align the domain distributions. For DI, we use an adversarial approach where a generator fill in information for missing data and a discriminator tries to distinguish between real and imputed values. Finally, using the universal feature representation in the embeddings, we train a classifier using MTL that given input from any domain, can predict labels for all domains. We demonstrate the superior performance of our approach compared to other state-of-art methods in three distinct settings, DG-DI in image recognition with unstructured data, MTL-DI in grade estimation with structured data and MDMTL-DI in a selection process using mixed data.


title: TraLFM: Latent Factor Modeling of Traffic Trajectory Data

abstract: The widespread use of positioning devices (e.g., GPS) has given rise to a vast body of human movement data, often in the form of trajectories. Understanding human mobility patterns could benefit many location-based applications. In this paper, we propose a novel generative model called TraLFM via latent factor modeling to mine human mobility patterns underlying traffic trajectories. TraLFM is based on three key observations: (1) human mobility patterns are reflected by the sequences of locations in the trajectories; (2) human mobility patterns vary with people; and (3) human mobility patterns tend to be cyclical and change over time. Thus, TraLFM models the joint action of sequential, personal and temporal factors in a unified way, and brings a new perspective to many applications such as latent factor analysis and next location prediction. We perform thorough empirical studies on two real datasets, and the experimental results confirm that TraLFM outperforms the state-of-the-art methods significantly in these applications.


title: TTDM: A Travel Time Difference Model for Next Location Prediction

abstract: Next location prediction is of great importance for many location-based applications and provides essential intelligence to business and governments. In existing studies, a common approach to next location prediction is to learn the sequential transitions with massive historical trajectories based on conditional probability. Unfortunately, due to the time and space complexity, these methods (e.g., Markov models) only use the just passed locations to predict next locations, without considering all the passed locations in the trajectory. In this paper, we seek to enhance the prediction performance by considering the travel time from all the passed locations in the query trajectory to a candidate next location. In particular, we propose a novel method, called Travel Time Difference Model (TTDM), which exploits the difference between the shortest travel time and the actual travel time to predict next locations. Further, we integrate the TTDM with a Markov model via a linear interpolation to yield a joint model, which computes the probability of reaching each possible next location and returns the top-rankings as results. We have conducted extensive experiments on two real datasets: the vehicle passage record (VPR) data and the taxi trajectory data. The experimental results demonstrate significant improvements in prediction accuracy over existing solutions. For example, compared with the Markov model, the top-1 accuracy improves by 40% on the VPR data and by 15.6% on the Taxi data.


title: MPE: A Mobility Pattern Embedding Model for Predicting Next Locations

abstract: The wide spread use of positioning and photographing devices gives rise to a deluge of traffic trajectory data (e.g., vehicle passage records and taxi trajectory data), with each record having at least three attributes: object ID, location ID, and time-stamp. In this paper, we propose a novel mobility pattern embedding model called MPE to shed the light on people's mobility patterns in traffic trajectory data from multiple aspects, including sequential, personal, and temporal factors. MPE has two salient features: (1) it is capable of casting various types of information (object, location and time) to an integrated low-dimensional latent space; (2) it considers the effect of ``phantom transitions'' arising from road networks in traffic trajectory data. This embedding model opens the door to a wide range of applications such as next location prediction and visualization. Experimental results on two real-world datasets show that MPE is effective and outperforms the state-of-the-art methods significantly in a variety of tasks.


title: The Shape of the Electric Dipole Function Determines the Sub-Picosecond Dynamics of Anharmonic Vibrational Polaritons

abstract: We develop a fully quantum mechanical methodology to describe the static properties and the dynamics of a single anharmonic vibrational mode interacting with a quantized infrared cavity field in the strong and ultrastrong coupling regimes. By comparing multiconfiguration time-dependent Hartree (MCTDH) simulations for a Morse oscillator in a cavity, with an equivalent formulation of the problem in Hilbert space, we describe for the first time the essential role of permanent dipole moments in the femtosecond dynamics of vibrational polariton wavepackets. We show that depending on the shape of the electric dipole function $d_e(q)$ along the vibrational mode coordinate $q$, molecules can be classified into three general families. For molecules that are polar and have a positive slope of the dipole function at equilibrium, we show that an initial diabatic light-matter product state without vibrational or cavity excitations can evolve into a polariton wavepacket with a large number of intracavity photons, for interaction strengths at the onset of ultrastrong coupling. This build up of intracavity photon amplitude is accompanied by an effective $lengthening$ of the vibrational mode of nearly $10\%$, comparable with a laser-induced vibrational excitation in free space. In contrast, molecules that are also polar at equilibrium but have a negative slope of the dipole function, experience an effective mode $shortening$ under equivalent coupling conditions. Our model predictions are numerically validated using realistic $ab$-$initio$ potentials and dipole functions for HF and CO$_2$ molecules in their ground electronic states. We finally propose a non-adiabatic state preparation scheme to generate vibrational polaritons using nanoscale infrared antennas and UV-vis photochemistry or electron tunneling, to enable the far-field detection of spontaneously generated infrared quantum light.


title: A novel Deep Structure U-Net for Sea-Land Segmentation in Remote Sensing Images

abstract: Sea-land segmentation is an important process for many key applications in remote sensing. Proper operative sea-land segmentation for remote sensing images remains a challenging issue due to complex and diverse transition between sea and lands. Although several Convolutional Neural Networks (CNNs) have been developed for sea-land segmentation, the performance of these CNNs is far from the expected target. This paper presents a novel deep neural network structure for pixel-wise sea-land segmentation, a Residual Dense U-Net (RDU-Net), in complex and high-density remote sensing images. RDU-Net is a combination of both down-sampling and up-sampling paths to achieve satisfactory results. In each down- and up-sampling path, in addition to the convolution layers, several densely connected residual network blocks are proposed to systematically aggregate multi-scale contextual information. Each dense network block contains multilevel convolution layers, short-range connections and an identity mapping connection which facilitates features re-use in the network and makes full use of the hierarchical features from the original images. These proposed blocks have a certain number of connections that are designed with shorter distance backpropagation between the layers and can significantly improve segmentation results whilst minimizing computational costs. We have performed extensive experiments on two real datasets Google Earth and ISPRS and compare the proposed RDUNet against several variations of Dense Networks. The experimental results show that RDUNet outperforms the other state-of-the-art approaches on the sea-land segmentation tasks.


title: Short-term oscillation and falling dynamics for a water drop dripping in quiescent air

abstract: The short-term transient falling dynamics of a dripping water drop in quiescent air has been investigated through both simulation and experiment. The focus is on the short term behavior and the time range considered covers about eight dominant second-mode oscillations of the drop after it is formed. Due to the small fluid inertia the growth of the drop is quasi-static and is well captured by the static pendant drop theory. Nevertheless, the pinching dynamics and the resulting post-formation state of the drop trigger a nonlinear oscillation when the drop falls. The initial shape of the drop when it is just formed is decomposed into spherical harmonic modes. The pinching dynamics such as interface overturning introduces small-scale variation on the drop contour, which in turn contributes to the finite amplitudes of the higher-order modes. Furthermore, the initial kinetic energy when the droplet is just formed is as important as the initial surface energy contained in the drop shape, and is found to amplify the initial oscillation amplitude and to induce a phase shift in the oscillation of all the modes. By incorporating both the initial surface and kinetic energy, the linear model for a free drop oscillation yields very good predictions for the second and third modes. The mode amplitude spectra show both the primary frequencies that are consistent with the Lamb's theory and the secondary frequencies arising from different modes due to nonlinear inter-mode coupling. The complex transient flow inside and outside the drop is induced by the interaction between the falling motion and the nonlinear oscillation. The streamlines indicate that the internal flow is substantially different from the Hill vortex for a falling drop without oscillation. The temporal evolutions of both the internal flow and the wake morphology follow the dominant second oscillation mode.


title: Quantal diffusion approach for multinucleon transfer processes in the $^{58,64}$Ni+$^{208}$Pb reactions: Toward the production of unknown neutron-rich nuclei

abstract: Background: In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This paper aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of $^{64}$Ni+$^{208}$Pb at $E_{c.m.}$=268 MeV and $^{58}$Ni+$^{208}$Pb at $E_{c.m.}$=270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from the time evolution of occupied single-particle orbitals in TDHF. Secondary deexcitation processes, both particle evaporation and fission, are simulated by a statistical model, GEMINI++. Results: We find that the quantal diffusion approach reproduces the measured cross sections for abundant reaction products, including both projectile-like and target-like fragments, as well as fragments of transfer-induced fission. The results underline the importance of beyond mean-field effects in describing multinucleon transfer processes. We emphasize that the quantal diffusion description does not involve any adjustable parameters, once an energy density functional is given. Possible future directions are discussed. (Shortened due to the word limit)


title: Store-Collect in the Presence of Continuous Churn with Application to Snapshots and Lattice Agreement

abstract: We present an algorithm for implementing a store-collect object in an asynchronous crash-prone message-passing dynamic system, where nodes continually enter and leave. The algorithm is very simple and efficient, requiring just one round trip for a store operation and two for a collect. We then show the versatility of the store-collect object for implementing churn-tolerant versions of useful data structures, while shielding the user from the complications of the underlying churn. In particular, we present elegant and efficient implementations of atomic snapshot and generalized lattice agreement objects that use store-collect.


title: Black holes in the four-dimensional Einstein-Lovelock gravity

abstract: A $(3+1)$-dimensional Einstein-Gauss-Bonnet theory of gravity has been recently formulated in [D. Glavan and C. Lin, Phys. Rev. Lett. {\bf 124}, 081301 (2020)] which is different from the pure Einstein theory, i.e., bypasses the Lovelock's theorem and avoids Ostrogradsky instability. The theory was formulated in $D > 4$ dimensions and its action consists of the Einstein-Hilbert term with a cosmological constant, while the Gauss-Bonnet term multiplied by a factor $1/(D-4)$. Then, the four-dimensional theory is defined as the limit $D \to 4$. Here we generalize this approach to the four-dimensional Einstein-Lovelock theory and formulate the most general static $4D$ black-hole solution allowing for a $\Lambda$-term (either positive or negative) and the electric charge $Q$. As metric functions cannot be found in a closed form in the general case, we develop and share publicly the code which constructs the metric functions for every given set of parameters.


title: Thermosolutal Marangoni instability in a viscoelastic liquid film: Effect of heating from the free surface

abstract: We study the Marangoni instability in a thin polymeric liquid film heated from the free surface. Polymeric solutions are usually the binary mixture of a Newtonian solvent with polymeric solute and exhibit viscoelastic behaviour. In the presence of a temperature gradient, stratification of these polymeric solutes can take place via the Soret effect, which may give rise to the solutocapillary effect at the free surface. Considering the Soret effect and incorporating the effects of gravity, here we analyse the stability characteristics of this polymeric thin film bounded between its deformable free surface and a poorly conductive rigid substrate from below. A linear stability analysis around the quiescent base state reveals that under the combined influences of thermo-solutocapillarity, apart from the monotonic disturbance, two different oscillatory instabilities can emerge in the system. The characteristics of each instability mode are discussed, and a complete stability picture is perceived in terms of the phase diagrams, identifying the model parameter space wherein a particular instability mode can get dominant.


title: DistNet: Deep Tracking by displacement regression: application to bacteria growing in the Mother Machine

abstract: The mother machine is a popular microfluidic device that allows long-term time-lapse imaging of thousands of cells in parallel by microscopy. It has become a valuable tool for single-cell level quantitative analysis and characterization of many cellular processes such as gene expression and regulation, mutagenesis or response to antibiotics. The automated and quantitative analysis of the massive amount of data generated by such experiments is now the limiting step. In particular the segmentation and tracking of bacteria cells imaged in phase-contrast microscopy---with error rates compatible with high-throughput data---is a challenging problem. In this work, we describe a novel formulation of the multi-object tracking problem, in which tracking is performed by a regression of the bacteria's displacement, allowing simultaneous tracking of multiple bacteria, despite their growth and division over time. Our method performs jointly segmentation and tracking, leveraging sequential information to increase segmentation accuracy. We introduce a Deep Neural Network architecture taking advantage of a self-attention mechanism which yields extremely low tracking error rate and segmentation error rate. We demonstrate superior performance and speed compared to state-of-the-art methods. Our method is named DiSTNet which stands for DISTance+DISplacement Segmentation and Tracking Network. While this method is particularly well suited for mother machine microscopy data, its general joint tracking and segmentation formulation could be applied to many other problems with different geometries.


title: Twisted conjugacy in fundamental groups of geometric $3$-manifolds

abstract: A group $G$ has property $R_\infty$ if for every $\phi\in Aut(G)$, there are an infinite number of $\phi$-twisted conjugacy classes of elements in $G$. In this note, we determine the $R_\infty$-property for $G=\pi_1(M)$ for all geometric $3$-manifolds $M$.


title: Fast Certification of Collision Probability Bounds with Uncertain Convex Obstacles

abstract: To operate reactively in uncertain environments, robots need to be able to quickly estimate the risk that they will collide with their environment. This ability is important for both planning (to ensure that plans maintain acceptable levels of safety) and execution (to provide real-time warnings when risk exceeds some threshold). Existing methods for estimating this risk are often limited to models with simplified geometry (e.g. point robots); others handle complex geometry but are too slow for many applications. In this paper, we present two algorithms for quickly computing upper bounds on the risk of collision between a robot and uncertain obstacles by searching for certificate regions that capture collision probability mass while avoiding the robot. These algorithms come with strong theoretical guarantees that the true risk does not exceed the estimated value, support arbitrary geometry via convex decomposition, and provide fast query times ($<200\mu$s) in representative scenarios. We characterize the performance of these algorithms in environments of varying complexity, demonstrating at least an order of magnitude speedup over existing techniques.


title: The Parameterized Complexity of Guarding Almost Convex Polygons

abstract: Art Gallery is a fundamental visibility problem in Computational Geometry. The input consists of a simple polygon P, (possibly infinite) sets G and C of points within P, and an integer k; the task is to decide if at most k guards can be placed on points in G so that every point in C is visible to at least one guard. In the classic formulation of Art Gallery, G and C consist of all the points within P. Other well-known variants restrict G and C to consist either of all the points on the boundary of P or of all the vertices of P. Recently, three new important discoveries were made: the above mentioned variants of Art Gallery are all W[1]-hard with respect to k [Bonnet and Miltzow, ESA'16], the classic variant has an O(log k)-approximation algorithm [Bonnet and Miltzow, SoCG'17], and it may require irrational guards [Abrahamsen et al., SoCG'17]. Building upon the third result, the classic variant and the case where G consists only of all the points on the boundary of P were both shown to be \exists R-complete~[Abrahamsen et al., STOC'18]. Even when both G and C consist only of all the points on the boundary of P, the problem is not known to be in NP. Given the first discovery, the following question was posed by Giannopoulos [Lorentz Center Workshop, 2016]: Is Art Gallery FPT with respect to r, the number of reflex vertices? We focus on the variant where G and C are all the vertices of P, called Vertex-Vertex Art Gallery. We show that Vertex-Vertex Art Gallery is solvable in time r^{O(r^2)}n^{O(1)}. Our approach also extends to assert that Vertex-Boundary Art Gallery and Boundary-Vertex Art Gallery are both FPT. We utilize structural properties of "almost convex polygons" to present a two-stage reduction from Vertex-Vertex Art Gallery to a new constraint satisfaction problem (whose solution is also provided in this paper) where constraints have arity 2 and involve monotone functions.


title: Research on a New Convolutional Neural Network Model Combined with Random Edges Adding

abstract: It is always a hot and difficult point to improve the accuracy of convolutional neural network model and speed up its convergence. Based on the idea of small world network, a random edge adding algorithm is proposed to improve the performance of convolutional neural network model. This algorithm takes the convolutional neural network model as a benchmark, and randomizes backwards and cross-layer connections with probability p to form a new convolutional neural network model. The proposed idea can optimize the cross layer connectivity by changing the topological structure of convolutional neural network, and provide a new idea for the improvement of the model. The simulation results based on Fashion-MINST and cifar10 data set show that the model recognition accuracy and training convergence speed are greatly improved by random edge adding reconstructed models with aprobability p = 0.1.


title: Time-dependent entropy of a cooling Bose gas

abstract: Exact analytic solutions of a nonlinear boson diffusion equation with suitable initial conditions that account for evaporative cooling of ultracold atoms, plus boundary conditions at the singularity $\epsilon=\mu<0$ are presented, and used to calculate the time-dependent entropy of a cold quantum gas.


title: Exceptional points in parity-time-symmetric subwavelength metamaterials

abstract: When sources of energy gain and loss are introduced to a wave-scattering system, the underlying mathematical formulation will be non-Hermitian. This paves the way for the existence of exceptional points, where eigenmodes are linearly dependent. The primary goal of this work is to study the existence of exceptional points in high-contrast subwavelength metamaterials. We begin by studying a parity-time-symmetric pair of subwavelength resonators and prove that this system supports asymptotic exceptional points. These are points at which the subwavelength eigenvalues and eigenvectors coincide at leading order in the asymptotic parameters. We then investigate further properties of parity-time-symmetric subwavelength metamaterials. First, we study the exotic scattering behaviour of a metascreen composed of repeating parity-time-symmetric pairs of subwavelength resonators. We prove that the non-Hermitian nature of this structure means that it exhibits asymptotic unidirectional reflectionless transmission at certain frequencies and demonstrate extraordinary transmission close to these frequencies. Thereafter, we consider cavities containing many small resonators and use homogenization theory to show that non-Hermitian behaviour can be replicated at the macroscale.


title: Hyperplane Arrangements of Trained ConvNets Are Biased

abstract: We investigate the geometric properties of the functions learned by trained ConvNets in the preactivation space of their convolutional layers, by performing an empirical study of hyperplane arrangements induced by a convolutional layer. We introduce statistics over the weights of a trained network to study local arrangements and relate them to the training dynamics. We observe that trained ConvNets show a significant statistical bias towards regular hyperplane configurations. Furthermore, we find that layers showing biased configurations are critical to validation performance for the architectures considered, trained on CIFAR10, CIFAR100 and ImageNet.


title: Pressio: Enabling projection-based model reduction for large-scale nonlinear dynamical systems

abstract: This work introduces Pressio, an open-source project aimed at enabling leading-edge projection-based reduced order models (ROMs) for large-scale nonlinear dynamical systems in science and engineering. Pressio provides model-reduction methods that can reduce both the number of spatial and temporal degrees of freedom for any dynamical system expressible as a system of parameterized ordinary differential equations (ODEs). We leverage this simple, expressive mathematical framework as a pivotal design choice to enable a minimal application programming interface (API) that is natural to dynamical systems. The core component of Pressio is a C++11 header-only library that leverages generic programming to support applications with arbitrary data types and arbitrarily complex programming models. This is complemented with Python bindings to expose these C++ functionalities to Python users with negligible overhead and no user-required binding code. We discuss the distinguishing characteristics of Pressio relative to existing model-reduction libraries, outline its key design features, describe how the user interacts with it, and present two test cases -- including one with over 20 million degrees of freedom -- that highlight the performance results of Pressio and illustrate the breath of problems that can be addressed with it.

categories: cs.MS cs.CE physics.comp-ph physics.flu-dyn


title: Marangoni instability in a viscoelastic binary film with cross-diffusive effect

abstract: The viscoelastic fluids are usually the blends of a polymeric solute and a Newtonian solvent. In the presence of a temperature gradient, stratification of these solutes can take place via the Soret effect. Here, we investigate the classical Marangoni instability problem for a thin viscoelastic film considering this binary aspect of the fluid. The film, bounded above by a deformable free surface, is subjected to heating from below by a solid substrate. Linear stability analysis performed numerically for perturbations of finite wavelength (short-wave perturbations) reveals that both monotonic and oscillatory instabilities can emerge in this system. The interaction between the thermocapillary and solutocapillary forces in the presence of Soret diffusion is found to give rise to two different oscillatory instabilities, of which one mode was overlooked previously, even for the Newtonian binary mixtures. As a principal result of this work, we provide a complete picture of the susceptibility to different instability modes based on the model parameter values. Finally, an approximate model is developed under the framework of long-wave analysis, which can qualitatively depict the stability behaviour of the system without numerically solving the problem.


title: When is Ontology-Mediated Querying Efficient?

abstract: In ontology-mediated querying, description logic (DL) ontologies are used to enrich incomplete data with domain knowledge which results in more complete answers to queries. However, the evaluation of ontology-mediated queries (OMQs) over relational databases is computationally hard. This raises the question when OMQ evaluation is efficient, in the sense of being tractable in combined complexity or fixed-parameter tractable. We study this question for a range of ontology-mediated query languages based on several important and widely-used DLs, using unions of conjunctive queries as the actual queries. For the DL ELHI extended with the bottom concept, we provide a characterization of the classes of OMQs that are fixed-parameter tractable. For its fragment EL extended with domain and range restrictions and the bottom concept (which restricts the use of inverse roles), we provide a characterization of the classes of OMQs that are tractable in combined complexity. Both results are in terms of equivalence to OMQs of bounded tree width and rest on a reasonable assumption from parameterized complexity theory. They are similar in spirit to Grohe's seminal characterization of the tractable classes of conjunctive queries over relational databases. We further study the complexity of the meta problem of deciding whether a given OMQ is equivalent to an OMQ of bounded tree width, providing several completeness results that range from NP to 2ExpTime, depending on the DL used. We also consider the DL-Lite family of DLs, including members that admit functional roles.


title: Virtual staining for mitosis detection in Breast Histopathology

abstract: We propose a virtual staining methodology based on Generative Adversarial Networks to map histopathology images of breast cancer tissue from H&E stain to PHH3 and vice versa. We use the resulting synthetic images to build Convolutional Neural Networks (CNN) for automatic detection of mitotic figures, a strong prognostic biomarker used in routine breast cancer diagnosis and grading. We propose several scenarios, in which CNN trained with synthetically generated histopathology images perform on par with or even better than the same baseline model trained with real images. We discuss the potential of this application to scale the number of training samples without the need for manual annotations.


title: The Implicit Regularization of Stochastic Gradient Flow for Least Squares

abstract: We study the implicit regularization of mini-batch stochastic gradient descent, when applied to the fundamental problem of least squares regression. We leverage a continuous-time stochastic differential equation having the same moments as stochastic gradient descent, which we call stochastic gradient flow. We give a bound on the excess risk of stochastic gradient flow at time $t$, over ridge regression with tuning parameter $\lambda = 1/t$. The bound may be computed from explicit constants (e.g., the mini-batch size, step size, number of iterations), revealing precisely how these quantities drive the excess risk. Numerical examples show the bound can be small, indicating a tight relationship between the two estimators. We give a similar result relating the coefficients of stochastic gradient flow and ridge. These results hold under no conditions on the data matrix $X$, and across the entire optimization path (not just at convergence).


title: SAR Tomography at the Limit: Building Height Reconstruction Using Only 3-5 TanDEM-X Bistatic Interferograms

abstract: Multi-baseline interferometric synthetic aperture radar (InSAR) techniques are effective approaches for retrieving the 3-D information of urban areas. In order to obtain a plausible reconstruction, it is necessary to use more than twenty interferograms. Hence, these methods are commonly not appropriate for large-scale 3-D urban mapping using TanDEM-X data where only a few acquisitions are available in average for each city. This work proposes a new SAR tomographic processing framework to work with those extremely small stacks, which integrates the non-local filtering into SAR tomography inversion. The applicability of the algorithm is demonstrated using a TanDEM-X multi-baseline stack with 5 bistatic interferograms over the whole city of Munich, Germany. Systematic comparison of our result with TanDEM-X raw digital elevation models (DEM) and airborne LiDAR data shows that the relative height accuracy of two third buildings is within two meters, which outperforms the TanDEM-X raw DEM. The promising performance of the proposed algorithm paved the first step towards high quality large-scale 3-D urban mapping.


title: Atomic Scale Insights Into The Mechanical Characteristics of Monolayer 1T-Titanium Disulphide: A Molecular Dynamics Study

abstract: In this work, we report on the mechanical responses and fracture behavior of pristine and defected monolayer 1T-Titanium Disulfide using classical molecular dynamics simulation. We investigated the effect of temperature, strain rate and defect ratio on the uniaxial tensile properties in both armchair and zigzag direction. We found that monolayer TiS2 shows isotropic uniaxial tensile properties except for failure strain which is greater in zigzag direction than armchair direction. We also observed a negative correlation of ultimate tensile strength, failure strain and young's modulus with temperature and defect ratio. Results depicts that strain rate has no effect on the young's modulus of monolayer TiS2 but higher strain rate results in higher ultimate tensile strength and failure strain.


title: A Flat-Panel Brightness Model for the Starlink Satellites and Measurement of their Absolute Visual Magnitude

abstract: The Starlink satellites are shaped like flat panels. The flat sides face zenith and nadir during normal operations. Their brightness is determined by the product of the solar illumination on the downward facing side of the panel multiplied by the area of that side projected toward the observer on Earth. This geometry leads to a unique brightness function that is not shared by other satellites. For example, the observed brightness is very sensitive to the solar elevation angle. There are circumstances where sunlight only illuminates the upward facing side of the satellite rendering it invisible to Earth-based observers. A brightness model depending on the solar aspect and the observer aspect of the flat panel, in addition to the satellite distance, is described. Absolute brightness is the only free parameter of the model, and it is taken to be that at a distance of 1,000 km when the solar and observer factors are unity. This model has been successfully fitted to a set of observed magnitudes. The absolute visual magnitude of a Starlink satellite as determined from this fitting is 4.1 +/- 0.1. The model could be used to determine the absolute magnitude of the Starlink satellite known as Dark Sat which has a special low-albedo coating.


title: Semi-abelian spectral data for singular fibres of the $\mathsf{SL}(2,\mathbb{C})$-Hitchin system

abstract: We describe spectral data for singular fibres of the $\mathsf{SL}(2,\mathbb{C})$-Hitchin fibration with irreducible and reduced spectral curve. Using Hecke transformations we give a stratification of these singular spaces by fibre bundles over Prym varieties. By analysing the parameter spaces of Hecke transformations this describes the singular Hitchin fibres as compactifications of abelian group bundles over abelian torsors. We prove that a large class of singular fibres are themselves fibre bundles over Prym varieties. As applications we study irreducible components of singular Hitchin fibres and give a description of $\mathsf{SL}(2,\mathbb{R})$-Higgs bundles in terms of these semi-abelian spectral data.


title: Conserved quantities and regularity in fluid dynamics

abstract: Conserved or dissipated quantities, like energy or entropy, are at the heart of the study of many classes of time-dependent PDEs in connection with fluid mechanics. This is the case, for instance, for the Euler and Navier-Stokes equations, for systems of conservation laws, and for transport equations. In all these cases, a formally conserved quantity may no longer be constant in time for a weak solution at low regularity. The delicate interplay between regularity and conservation of the respective quantity relates to renormalisation in the DiPerna-Lions theory of transport and continuity equations, and to Onsager's conjecture in the realm of ideal incompressible fluids. We will review the classical commutator methods of DiPerna-Lions and Constantin-E-Titi, and then proceed to more recent results.


title: Partial D-operators for the generalized IBP reduction

abstract: Empirical evidence reveals existence of partial D-operators for the generalized IBP (BT) reduction algorithms that are, counterintuitively, much simpler and much easier to find than the complete D-operators from the foundational Bernstein theorem, allowing one to construct first true two-loop examples of generalized IBP identities.


title: Some new rational Gushel fourfolds

abstract: We provide explicit equations and parameterizations of some new rational Gushel-Mukai fourfolds of special type.


title: A Spectral Approach to Network Design

abstract: We present a spectral approach to design approximation algorithms for network design problems. We observe that the underlying mathematical questions are the spectral rounding problems, which were studied in spectral sparsification and in discrepancy theory. We extend these results to incorporate additional non-negative linear constraints, and show that they can be used to significantly extend the scope of network design problems that can be solved. Our algorithm for spectral rounding is an iterative randomized rounding algorithm based on the regret minimization framework. In some settings, this provides an alternative spectral algorithm to achieve constant factor approximation for the classical survivable network design problem, and partially answers a question of Bansal about survivable network design with concentration property. We also show many other applications of the spectral rounding results, including weighted experimental design and additive spectral sparsification.


title: Provably Safe Trajectory Optimization in the Presence of Uncertain Convex Obstacles

abstract: Real-world environments are inherently uncertain, and to operate safely in these environments robots must be able to plan around this uncertainty. In the context of motion planning, we desire systems that can maintain an acceptable level of safety as the robot moves, even when the exact locations of nearby obstacles are not known. In this paper, we solve this chance-constrained motion planning problem using a sequential convex optimization framework. To constrain the risk of collision incurred by planned movements, we employ geometric objects called $\epsilon$-shadows to compute upper bounds on the risk of collision between the robot and uncertain obstacles. We use these $\epsilon$-shadow-based estimates as constraints in a nonlinear trajectory optimization problem, which we then solve by iteratively linearizing the non-convex risk constraints. This sequential optimization approach quickly finds trajectories that accomplish the desired motion while maintaining a user-specified limit on collision risk. Our method can be applied to robots and environments with arbitrary convex geometry; even in complex environments, it runs in less than a second and provides provable guarantees on the safety of planned trajectories, enabling fast, reactive, and safe robot motion in realistic environments.


title: Phase diagram of the $\mathbb{Z}_3$-Fock parafermion chain with pair hopping

abstract: We study a tight binding model of $\mathbb{Z}_3$-Fock parafermions with single-particle and pair-hopping terms. The phase diagram has four different phases: a gapped phase, a gapless phase with central charge $c=2$, and two gapless phases with central charge $c=1$. We characterise each phase by analysing the energy gap, entanglement entropy and different correlation functions. The numerical simulations are complemented by analytical arguments.


title: Enhancing the Monte Carlo Tree Search Algorithm for Video Game Testing

abstract: In this paper, we study the effects of several Monte Carlo Tree Search (MCTS) modifications for video game testing. Although MCTS modifications are highly studied in game playing, their impacts on finding bugs are blank. We focused on bug finding in our previous study where we introduced synthetic and human-like test goals and we used these test goals in Sarsa and MCTS agents to find bugs. In this study, we extend the MCTS agent with several modifications for game testing purposes. Furthermore, we present a novel tree reuse strategy. We experiment with these modifications by testing them on three testbed games, four levels each, that contain 45 bugs in total. We use the General Video Game Artificial Intelligence (GVG-AI) framework to create the testbed games and collect 427 human tester trajectories using the GVG-AI framework. We analyze the proposed modifications in three parts: we evaluate their effects on bug finding performances of agents, we measure their success under two different computational budgets, and we assess their effects on human-likeness of the human-like agent. Our results show that MCTS modifications improve the bug finding performance of the agents.


title: Weight $q$-multiplicities for representations of the exceptional Lie algebra $\mathfrak{g}_2$

abstract: Given a simple Lie algebra $\mathfrak{g}$, Kostant's weight $q$-multiplicity formula is an alternating sum over the Weyl group whose terms involve the $q$-analog of Kostant's partition function. For $\xi$ (a weight of $\mathfrak{g}$), the $q$-analog of Kostant's partition function is a polynomial-valued function defined by $\wp_q(\xi)=\sum c_i q^i$ where $c_i$ is the number of ways $\xi$ can be written as a sum of $i$ positive roots of $\mathfrak{g}$. In this way, the evaluation of Kostant's weight $q$-multiplicity formula at $q = 1$ recovers the multiplicity of a weight in a highest weight representation of $\mathfrak{g}$. In this paper, we give closed formulas for computing weight $q$-multiplicities in a highest weight representation of the exceptional Lie algebra $\mathfrak{g}_2$.