title: Entanglement Induced by Noncommutativity: Anisotropic Harmonic Oscillator in Noncommutative space

abstract: Quantum entanglement, induced by spatial noncommutativity, is investigated for an anisotropic harmonic oscillator. Exact solutions for the system are obtained after the model is re-expressed in terms of canonical variables, by performing a particular Bopp's shift to the noncommuting degrees of freedom. Employing Simon's separability criterion, we find that the states of the system are entangled provided a unique function of the (mass and frequency) parameters obeys an inequality. Entanglement of Formation for this system is also computed and its relation to the degree of anisotropy is discussed. It is worth mentioning that, even in a noncommutative space, entanglement is generated only if the harmonic oscillator is anisotropic. Interestingly, the Entanglement of Formation saturates for higher values of the deformation parameter $\theta$, that quantifies spatial noncommutativity.


title: Lachesis: Automatic Partitioning for UDF-Centric Analytics

abstract: Persistent partitioning is effective in avoiding expensive shuffling operations. However it remains a significant challenge to automate this process for Big Data analytics workloads that extensively use user defined functions (UDFs), where sub-computations are hard to be reused for partitionings compared to relational applications. In addition, functional dependency that is widely utilized for partitioning selection is often unavailable in the unstructured data that is ubiquitous in UDF-centric analytics. We propose the Lachesis system, which represents UDF-centric workloads as workflows of analyzable and reusable sub-computations. Lachesis further adopts a deep reinforcement learning model to infer which sub-computations should be used to partition the underlying data. This analysis is then applied to automatically optimize the storage of the data across applications to improve the performance and users' productivity.


title: Moduli spaces of semistable pairs on projective Deligne-Mumford stacks

abstract: We generalize the construction of a moduli space of semistable pairs parametrizing isomorphism classes of morphisms from a fixed coherent sheaf to any sheaf with fixed Hilbert polynomial under a notion of stability to the case of projective Deligne-Mumford stacks. We study the deformation and obstruction theories of stable pairs, and then prove the existence of virtual fundamental classes for some cases of dimension two and three. This leads to a definition of Pandharipande-Thomas invariants on three-dimensional smooth projective Deligne-Mumford stacks.


title: Sliced Kernelized Stein Discrepancy

abstract: Kernelized Stein discrepancy (KSD), though being extensively used in goodness-of-fit tests and model learning, suffers from the curse-of-dimensionality. We address this issue by proposing the sliced Stein discrepancy and its scalable and kernelized variants, which employ kernel-based test functions defined on the optimal one-dimensional projections. When applied to goodness-of-fit tests, extensive experiments show the proposed discrepancy significantly outperforms KSD and various baselines in high dimensions. For model learning, we show its advantages over existing Stein discrepancy baselines by training independent component analysis models with different discrepancies. We further propose a novel particle inference method called sliced Stein variational gradient descent (S-SVGD) which alleviates the mode-collapse issue of SVGD in training variational autoencoders.


title: Integrating Out New Fermions at One Loop

abstract: We present the fermionic universal one--loop effective action obtained by integrating out heavy vector--like fermions at one loop using functional techniques. Even though previous approaches are able to handle integrating out heavy fermions with non--chiral interactions, i.e. vanishing $\gamma^5$ interaction terms, the computations proceed in a tedious manner that obscures a physical interpretation. We show how directly tackling the fermionic functional determinant not only allows for a much simpler and transparent computation, but is also able to account for chiral interaction terms in a simple, algorithmic way. Finally, we apply the obtained results to integrate out at one loop the vector--like fermions appearing in a toy model and in a fermionic model that exhibits strong cosmological phase transitions.


title: Actionable Attribution Maps for Scientific Machine Learning

abstract: The scientific community has been increasingly interested in harnessing the power of deep learning to solve various domain challenges. However, despite the effectiveness in building predictive models, fundamental challenges exist in extracting actionable knowledge from the deep neural network due to their opaque nature. In this work, we propose techniques for exploring the behavior of deep learning models by injecting domain-specific actionable concepts as tunable ``knobs'' in the analysis pipeline. By incorporating the domain knowledge with generative modeling, we are not only able to better understand the behavior of these black-box models, but also provide scientists with actionable insights that can potentially lead to fundamental discoveries.


title: Next Generation Wi-Fi and 5G NR-U in the 6 GHz Bands: Opportunities & Challenges

abstract: The ever-increasing demand for unlicensed spectrum has prompted regulators in the US and Europe to consider opening up the 6 GHz bands for unlicensed access. These bands will open up 1.2 GHz of additional spectrum for unlicensed radio access technologies (RATs), such as Wi-Fi and 5G New Radio Unlicensed (NR-U), in the US and if permitted, 500 MHz of additional spectrum in Europe. The abundance of spectrum in these bands creates new opportunities for the design of mechanisms and features that can support the emerging bandwidth-intensive and latency-sensitive applications. However, coexistence of unlicensed devices both with the bands' incumbent users and across different unlicensed RATs present significant challenges. In this paper, we provide a comprehensive survey of the existing literature on various issues surrounding the operations of unlicensed RATs in the 6 GHz bands. In particular, we discuss how key features in next-generation Wi-Fi are being designed to leverage these additional unlicensed bands. We also shed light on the foreseeable challenges that designers of unlicensed RATs might face in the near future. Our survey encompasses key research papers, contributions submitted to standardization bodies and regulatory agencies, and documents presented at various other venues. Finally, we highlight a few key research problems that are likely to arise due to unlicensed operations in the 6 GHz bands. Tackling these research challenges effectively will be critical in ensuring that the new unlicensed bands are efficiently utilized while guaranteeing the interference-free operation of the bands' incumbent users.


title: RCP: A Low-overhead Reversible Coherence Protocol

abstract: This paper proposes RCP, a new reversible coherence protocol that ensures invisible speculative load execution (ISLE) with low overhead. RCP can be combined with processor mechanisms that eliminate the effects of speculative instructions on other instructions to achieve low overhead invisible speculative execution (ISE). ISE provides protection that is at least as strong as speculative privacy tracking (SPT) and stronger than speculative taint tracking (STT). RCP is designed by systematically extending the existing coherence protocol to incorporate speculative loads and states. The protocol is implemented in Gem5 and verified with Murphi. The results show that RCP based ISE incurs lower overhead than STT/SDO/SPT while providing similar/stronger protection.


title: A counterexample to vanishing conjectures for negative $K$-theory

abstract: In a 2006 article Schlichting conjectured that the negative {\it K--}theory of any abelian category must vanish. This conjecture was generalized in a 2019 article by Antieau, Gepner and Heller, who hypothesized that the negative {\it K--}theory of any category with a bounded {\it t--}structure must vanish. Both conjectures will be shown to be false.


title: Theory-Inspired Path-Regularized Differential Network Architecture Search

abstract: Despite its high search efficiency, differential architecture search (DARTS) often selects network architectures with dominated skip connections which lead to performance degradation. However, theoretical understandings on this issue remain absent, hindering the development of more advanced methods in a principled way. In this work, we solve this problem by theoretically analyzing the effects of various types of operations, e.g. convolution, skip connection and zero operation, to the network optimization. We prove that the architectures with more skip connections can converge faster than the other candidates, and thus are selected by DARTS. This result, for the first time, theoretically and explicitly reveals the impact of skip connections to fast network optimization and its competitive advantage over other types of operations in DARTS. Then we propose a theory-inspired path-regularized DARTS that consists of two key modules: (i) a differential group-structured sparse binary gate introduced for each operation to avoid unfair competition among operations, and (ii) a path-depth-wise regularization used to incite search exploration for deep architectures that often converge slower than shallow ones as shown in our theory and are not well explored during the search. Experimental results on image classification tasks validate its advantages.


title: Near-field cavity optomechanical coupling in a compound semiconductor nanowire

abstract: A III-V compound semiconductor nanowire is an attractive material for a novel hybrid quantum interface that interconnects photons, electrons, and phonons through a wavelength-tunable quantum structure embedded in its free-standing structure. In such a nanomechanical element, however, a challenge is how to detect and manipulate a small number of phonons via its tiny mechanical motion. A solution would be to couple an optical cavity to a nanowire by introducing the ``cavity optomechanics'' framework, but the typical size difference between them becomes a barrier to achieving this. Here, we demonstrate near-field coupling of a silica microsphere cavity and an epitaxially grown InP/InAs free-standing nanowire. The evanescent optomechanical coupling enables not only fine probing of the mechanical motion by balanced homodyne interferometry but also tuning of the resonance frequency, linewidth, Duffing nonlinearity, and vibration axis in the nanowire. Combining this cavity optomechanics with epitaxial nanowire engineering opens the way to novel quantum metrology and information processing.


title: Autoregressive Mixture Models for Serial Correlation Clustering of Time Series Data

abstract: Clustering time series into similar groups can improve models by combining information across like time series. While there is a well developed body of literature for clustering of time series, these approaches tend to generate clusters independently of model training which can lead to poor model fit. We propose a novel distributed approach that simultaneously clusters and fits autoregression models for groups of similar individuals. We apply a Wishart mixture model so as to cluster individuals while modeling the corresponding autocovariance matrices at the same time. The fitted Wishart scale matrices map to cluster-level autoregressive coefficients through the Yule-Walker equations, fitting robust parsimonious autoregressive mixture models. This approach is able to discern differences in underlying autocorrelation variation of time series in settings with large heterogeneous datasets. We prove consistency of our cluster membership estimator, asymptotic distributions of coefficients and compare our approach against competing methods through simulation as well as by fitting a COVID-19 forecast model.


title: Associative Memory in Iterated Overparameterized Sigmoid Autoencoders

abstract: Recent work showed that overparameterized autoencoders can be trained to implement associative memory via iterative maps, when the trained input-output Jacobian of the network has all of its eigenvalue norms strictly below one. Here, we theoretically analyze this phenomenon for sigmoid networks by leveraging recent developments in deep learning theory, especially the correspondence between training neural networks in the infinite-width limit and performing kernel regression with the Neural Tangent Kernel (NTK). We find that overparameterized sigmoid autoencoders can have attractors in the NTK limit for both training with a single example and multiple examples under certain conditions. In particular, for multiple training examples, we find that the norm of the largest Jacobian eigenvalue drops below one with increasing input norm, leading to associative memory.


title: AdaSGD: Bridging the gap between SGD and Adam

abstract: In the context of stochastic gradient descent(SGD) and adaptive moment estimation (Adam),researchers have recently proposed optimization techniques that transition from Adam to SGD with the goal of improving both convergence and generalization performance. However, precisely how each approach trades off early progress and generalization is not well understood; thus, it is unclear when or even if, one should transition from one approach to the other. In this work, by first studying the convex setting, we identify potential contributors to observed differences in performance between SGD and Adam. In particular,we provide theoretical insights for when and why Adam outperforms SGD and vice versa. We ad-dress the performance gap by adapting a single global learning rate for SGD, which we refer to as AdaSGD. We justify this proposed approach with empirical analyses in non-convex settings. On several datasets that span three different domains,we demonstrate how AdaSGD combines the benefits of both SGD and Adam, eliminating the need for approaches that transition from Adam to SGD.


title: Stress-induced ordering of two-dimensional packings of elastic spheres

abstract: Packing of particles in confined environments is a common problem in multiple fields. Here, based on the two-dimensional Hertzian particle model, we study the packing of deformable spherical particles under compression, and reveal the crucial role of stress as an ordering field in regulating particle arrangement. Specifically, under increasing compression, the squeezed particles spontaneously organize into regular packings in the sequence of triangular and square lattices, pentagonal tessellation, and the reentrant triangular lattice. The rich ordered patterns and complex structures revealed in this work suggest a fruitful organizational strategy based on the interplay of external stress and intrinsic elastic instability of particle arrays.


title: Coherent detection schemes for subcarrier wave continuous variable quantum key distribution

abstract: We examine different methods to implement coherent detection in the subcarrier wave quantum key distribution (SCW QKD) systems. For classical wavefields, we present the models describing homodyne-type and heterodyne-type coherent detection schemes needed to extract information from the quadrature phase-coded multimode signals used in SCW QKD. Practical feasibility of the proposed schemes is corroborated by the experiments.


title: Properly colored Hamilton cycles in Dirac-type hypergraphs

abstract: We consider a robust variant of Dirac-type problems in $k$-uniform hypergraphs. For instance, we prove that if $H$ is a $k$-uniform hypergraph with minimum codegree at least $(1/2 + \gamma )n$, $\gamma >0$, and $n$ is sufficiently large, then any edge coloring $\phi$ satisfying appropriate local constraints yields a properly colored tight Hamilton cycle in $H$. Similar results for loose cycles are also shown.


title: Adversarial Deep Ensemble: Evasion Attacks and Defenses for Malware Detection

abstract: Malware remains a big threat to cyber security, calling for machine learning based malware detection. While promising, such detectors are known to be vulnerable to evasion attacks. Ensemble learning typically facilitates countermeasures, while attackers can leverage this technique to improve attack effectiveness as well. This motivates us to investigate which kind of robustness the ensemble defense or effectiveness the ensemble attack can achieve, particularly when they combat with each other. We thus propose a new attack approach, named mixture of attacks, by rendering attackers capable of multiple generative methods and multiple manipulation sets, to perturb a malware example without ruining its malicious functionality. This naturally leads to a new instantiation of adversarial training, which is further geared to enhancing the ensemble of deep neural networks. We evaluate defenses using Android malware detectors against 26 different attacks upon two practical datasets. Experimental results show that the new adversarial training significantly enhances the robustness of deep neural networks against a wide range of attacks, ensemble methods promote the robustness when base classifiers are robust enough, and yet ensemble attacks can evade the enhanced malware detectors effectively, even notably downgrading the VirusTotal service.


title: Hand-drawn Symbol Recognition of Surgical Flowsheet Graphs with Deep Image Segmentation

abstract: Perioperative data are essential to investigating the causes of adverse surgical outcomes. In some low to middle income countries, these data are computationally inaccessible due to a lack of digitization of surgical flowsheets. In this paper, we present a deep image segmentation approach using a U-Net architecture that can detect hand-drawn symbols on a flowsheet graph. The segmentation mask outputs are post-processed with techniques unique to each symbol to convert into numeric values. The U-Net method can detect, at the appropriate time intervals, the symbols for heart rate and blood pressure with over 99 percent accuracy. Over 95 percent of the predictions fall within an absolute error of five when compared to the actual value. The deep learning model outperformed template matching even with a small size of annotated images available for the training set.


title: Escape of a passive particle from activity-induced energy landscape: Emergence of slow and fast effective diffusion

abstract: Spontaneous persistent motions driven by active processes play a central role to maintain the living cells far from equilibrium. In the majority of the research works, the steady state dynamics of an active system has been described in terms of an effective temperature. In the majority of the research works, the steady state dynamics of an active system has been described in terms of an effective temperature. By contrast, we have examined a prototype model for diffusion in an activity-induced rugged energy landscape to describe the slow dynamics of a tagged particle in a dense active environment. The expression for the mean escape time from the active rugged energy landscape holds only in the limit of low activity and the mean escape time from the rugged energy landscape increases with activity. The precise form of the active correlation will determine whether the mean escape time will depend on the persistence time or not. The active rugged energy landscape approach also allows an estimate of non-equilibrium effective diffusivity characterizing the slow diffusive motion of the tagged particle due to activity. On the other hand, in a dilute environment, high activity augments the diffusion of the tagged particle. The enhanced diffusion can be attributed to an effective temperature, higher than the ambient temperature and is used to calculate the Kramers' mean escape time, which decreases with activity. Our results have direct relevance to recent experiments on tagged particle diffusion in condensed phases.


title: Sinkhorn EM: An Expectation-Maximization algorithm based on entropic optimal transport

abstract: We study Sinkhorn EM (sEM), a variant of the expectation maximization (EM) algorithm for mixtures based on entropic optimal transport. sEM differs from the classic EM algorithm in the way responsibilities are computed during the expectation step: rather than assign data points to clusters independently, sEM uses optimal transport to compute responsibilities by incorporating prior information about mixing weights. Like EM, sEM has a natural interpretation as a coordinate ascent procedure, which iteratively constructs and optimizes a lower bound on the log-likelihood. However, we show theoretically and empirically that sEM has better behavior than EM: it possesses better global convergence guarantees and is less prone to getting stuck in bad local optima. We complement these findings with experiments on simulated data as well as in an inference task involving C. elegans neurons and show that sEM learns cell labels significantly better than other approaches.


title: Systems of parameters and the Cohen--Macaulay property

abstract: We recall a numerical criteria for Cohen--Macaulayness related to system of parameters, and introduce monomial ideals of K\"onig type which include the edge ideals of K\"onig graphs. We show that a monomial ideal is of K\"onig type if and only if its corresponding residue class ring admits a system of parameters whose elements are of the form $x_i-x_j$. This provides an algebraic characterization of K\"onig graphs. We use this special parameter systems for the study of the edge ideal of K\"onig graphs and the study of the order complex of a certain family of posets. Finally, for any simplicial complex $\Delta$ we introduce a system of parameters for $K[\Delta]$ with a universal construction principle, independent of the base field and only dependent on the faces of $\Delta$. This system of parameters is an efficient tool to test Cohen--Macaulayness of the Stanley--Reisner ring of a simplicial complex.


title: Parameter Estimation of Fire Propagation Models Using Level Set Methods

abstract: The availability of wildland fire propagation models with parameters estimated in an accurate way starting from measurements of fire fronts is crucial to predict the evolution of fire and allocate resources for firefighting. Thus, we propose an approach to estimate the parameters of a wildland fire propagation model combining an empirical fire spread rate and level set methods to describe the evolution of the fire front over time and space. After validating the model, the estimation of parameters in the spread rate is performed by using fire front shapes measured at different time instants as well as wind velocity and direction, landscape elevation, and vegetation distribution. Parameter estimation is performed by solving an optimization problem, where the objective function to be minimized is the symmetric difference between predicted and measured fronts. Numerical results obtained by the application of the proposed method are reported in two simulated scenarios and in an application case study using real data of the 2002 Troy fire in Southern California, thus showing the effectiveness of the proposed approach.


title: Hierarchical Graph Matching Network for Graph Similarity Computation

abstract: Graph edit distance / similarity is widely used in many tasks, such as graph similarity search, binary function analysis, and graph clustering. However, computing the exact graph edit distance (GED) or maximum common subgraph (MCS) between two graphs is known to be NP-hard. In this paper, we propose the hierarchical graph matching network (HGMN), which learns to compute graph similarity from data. HGMN is motivated by the observation that two similar graphs should also be similar when they are compressed into more compact graphs. HGMN utilizes multiple stages of hierarchical clustering to organize a graph into successively more compact graphs. At each stage, the earth mover distance (EMD) is adopted to obtain a one-to-one mapping between the nodes in two graphs (on which graph similarity is to be computed), and a correlation matrix is also derived from the embeddings of the nodes in the two graphs. The correlation matrices from all stages are used as input for a convolutional neural network (CNN), which is trained to predict graph similarity by minimizing the mean squared error (MSE). Experimental evaluation on 4 datasets in different domains and 4 performance metrics shows that HGMN consistently outperforms existing baselines in the accuracy of graph similarity approximation.


title: Soft X-ray Vortex Beam detected by In-line Holography

abstract: We demonstrate the in-line holography for soft x-ray vortex beam having an orbital angular momentum. A hologram is recorded as an interference between a Bragg diffraction wave from a fork grating and a divergence wave generated by a Fresnel zone plate. The obtained images exhibit fork-shaped interference fringes, which confirms the formation of the vortex beam. By analyzing the interference image, we successfully obtained the spiral phase distribution. The results demonstrate that the in-line holography technique is promising for the characterization of topological magnets, such as magnetic skyrmions.


title: On Ulrich bundles on projective bundles

abstract: In this article, the existence of Ulrich bundles on projective bundles $\mathbb P(E) \to X$ is discussed. In the case, that the base variety $X$ is a curve or surface, a close relationship between Ulrich bundles on $X$ and those on $\mathbb P(E)$ is established for specific polarisations. This yields the existence of Ulrich bundles on a wide range of projective bundles over curves and surfaces.


title: Security Issues of Low Power Wide Area Networks in the Context of LoRa Networks

abstract: Low Power Wide Area Networks (LPWAN) have been used to support low cost and mobile bi-directional communications for the Internet of Things (IoT), smart city and a wide range of industrial applications. A primary security concern of LPWAN technology is the attacks that block legitimate communication between nodes resulting in scenarios like loss of packets, delayed packet arrival, and skewed packet reaching the reporting gateway. LoRa (Long Range) is a promising wireless radio access technology that supports long-range communication at low data rates and low power consumption. LoRa is considered as one of the ideal candidates for building LPWANs. We use LoRa as a reference technology to review the IoT security threats on the air and the applicability of different countermeasures that have been adopted so far. LoRa nodes that are close to the gateway use a small SF than the nodes which are far away. But it also implies long in-the-air transmission time, which makes the transmitted packets vulnerable to different kinds of malicious attacks, especially in the physical and the link layer. Therefore, it is not possible to enforce a fixed set of rules for all LoRa nodes since they have different levels of vulnerabilities. Our survey reveals that there is an urgent need for secure and uninterrupted communication between an end-device and the gateway, especially when the threat models are unknown in advance. We explore the traditional countermeasures and find that most of them are ineffective now, such as frequency hopping and spread spectrum methods. In order to adapt to new threats, the emerging countermeasures using game-theoretic approaches and reinforcement machine learning methods can effectively identify threats and dynamically choose the corresponding actions to resist threats, thereby making secured and reliable communications.


title: Analysis of Close Stellar Encounters with the Solar System

abstract: We consider 36 candidates for close (within 1 pc) encounters with the Solar system. These stars have been selected in accordance with the results of an analysis of their motion obtained by various authors. For most stars from this list the kinematic characteristics have been taken from the Gaia DR2 catalogue. The parameters of the encounters of these stars with the Solar system have been calculated using three methods: (1) the linear one, (2) the epicyclic one, and (3) by integrating the orbits in an axisymmetric potential. We have concluded that the epicyclic method works quite well only on a time interval no longer than $\pm$1 Myr. Based on the third method, in good agreement with the first method, for the first time we have obtained the following estimates of the encounter parameters for the star Gaia DR2 3130033734235815424: $t_{min}=-0.62\pm0.12$ Myr and $d_{min}=0.30\pm0.10$ pc.


title: Graph Neural Networks for Leveraging Industrial Equipment Structure: An application to Remaining Useful Life Estimation

abstract: Automated equipment health monitoring from streaming multisensor time-series data can be used to enable condition-based maintenance, avoid sudden catastrophic failures, and ensure high operational availability. We note that most complex machinery has a well-documented and readily accessible underlying structure capturing the inter-dependencies between sub-systems or modules. Deep learning models such as those based on recurrent neural networks (RNNs) or convolutional neural networks (CNNs) fail to explicitly leverage this potentially rich source of domain-knowledge into the learning procedure. In this work, we propose to capture the structure of a complex equipment in the form of a graph, and use graph neural networks (GNNs) to model multi-sensor time-series data. Using remaining useful life estimation as an application task, we evaluate the advantage of incorporating the graph structure via GNNs on the publicly available turbofan engine benchmark dataset. We observe that the proposed GNN-based RUL estimation model compares favorably to several strong baselines from literature such as those based on RNNs and CNNs. Additionally, we observe that the learned network is able to focus on the module (node) with impending failure through a simple attention mechanism, potentially paving the way for actionable diagnosis.


title: On the Fundamental Limits of Coded Caching Systems with Restricted Demand Types

abstract: Caching is a technique to reduce the communication load in peak hours by prefetching contents during off-peak hours. An information-theoretic framework for coded caching was introduced by Maddah-Ali and Niesen in a recent work, where it was shown that significant improvement can be obtained compared to uncoded caching. Considerable efforts have been devoted to identify the precise information-theoretic fundamental limits of the coded caching systems, however the difficulty of this task has also become clear. One of the reasons for this difficulty is that the original coded caching setting allows all possible multiple demand types during delivery, which in fact introduces tension in the coding strategy. In this paper, we seek to develop a better understanding of the fundamental limits of coded caching by investigating systems with certain demand type restrictions. We first consider the canonical three-user three-file system, and show that, contrary to popular beliefs, the worst demand type is not the one in which all three files are requested. Motivated by these findings, we focus on coded caching systems where every file must be requested by at least one user. A novel coding scheme is proposed, which can provide new operating points that are not covered by any previously known schemes.


title: Enabling Continual Learning with Differentiable Hebbian Plasticity

abstract: Continual learning is the problem of sequentially learning new tasks or knowledge while protecting previously acquired knowledge. However, catastrophic forgetting poses a grand challenge for neural networks performing such learning process. Thus, neural networks that are deployed in the real world often struggle in scenarios where the data distribution is non-stationary (concept drift), imbalanced, or not always fully available, i.e., rare edge cases. We propose a Differentiable Hebbian Consolidation model which is composed of a Differentiable Hebbian Plasticity (DHP) Softmax layer that adds a rapid learning plastic component (compressed episodic memory) to the fixed (slow changing) parameters of the softmax output layer; enabling learned representations to be retained for a longer timescale. We demonstrate the flexibility of our method by integrating well-known task-specific synaptic consolidation methods to penalize changes in the slow weights that are important for each target task. We evaluate our approach on the Permuted MNIST, Split MNIST and Vision Datasets Mixture benchmarks, and introduce an imbalanced variant of Permuted MNIST -- a dataset that combines the challenges of class imbalance and concept drift. Our proposed model requires no additional hyperparameters and outperforms comparable baselines by reducing forgetting.


title: Complex-birefringent dielectric metasurfaces for arbitrary polarization-pair transformations

abstract: Birefringent materials or nanostructures that introduce phase differences between two linear polarizations underpin the operation of wave plates for polarization control of light. Here we develop metasurfaces realizing a distinct class of complex-birefringent wave plates, which combine polarization transformation with a judiciously tailored polarization-dependent phase retardance and amplitude filtering via diffraction. We prove that the presence of loss enables the mapping from any chosen generally non-orthogonal pair of polarizations to any other pair at the output. We establish an optimal theoretical design-framework based on pairwise nanoresonator structures and experimentally demonstrate unique properties of metasurfaces in the amplification of small polarization differences and polarization coupling with unconventional phase control. Furthermore, we reveal that these metasurfaces can perform arbitrary transformations of biphoton polarization-encoded quantum states, including the modification of the degree of entanglement. Thereby, such flat devices can facilitate novel types of multi-functional polarization optics for classical and quantum applications.


title: PeTar: a high-performance N-body code for modeling massive collisional stellar systems

abstract: The numerical simulations of massive collisional stellar systems, such as globular clusters (GCs), are very time-consuming. Until now, only a few realistic million-body simulations of GCs with a small fraction of binaries (5%) have been performed by using the NBODY6++GPU code. Such models took half a year computational time on a GPU based super-computer. In this work, we develop a new N-body code, PeTar, by combining the methods of Barnes-Hut tree, Hermite integrator and slow-down algorithmic regularization (SDAR). The code can accurately handle an arbitrary fraction of multiple systems (e.g. binaries, triples) while keeping a high performance by using the hybrid parallelization methods with MPI, OpenMP, SIMD instructions and GPU. A few benchmarks indicate that PeTar and NBODY6++GPU have a very good agreement on the long-term evolution of the global structure, binary orbits and escapers. On a highly configured GPU desktop computer, the performance of a million-body simulation with all stars in binaries by using PeTar is 11 times faster than that of NBODY6++GPU. Moreover, on the Cray XC50 supercomputer, PeTar well scales when number of cores increase. The ten million-body problem, which covers the region of ultra compact dwarfs and nuclearstar clusters, becomes possible to be solved.


title: From Poincar\'e Inequalities to Nonlinear Matrix Concentration

abstract: This paper deduces exponential matrix concentration from a Poincar\'e inequality via a short, conceptual argument. Among other examples, this theory applies to matrix-valued functions of a uniformly log-concave random vector. The proof relies on the subadditivity of Poincar\'e inequalities and a chain rule inequality for the trace of the matrix Dirichlet form. It also uses a symmetrization technique to avoid difficulties associated with a direct extension of the classic scalar argument.


title: Nonlinear Matrix Concentration via Semigroup Methods

abstract: Matrix concentration inequalities provide information about the probability that a random matrix is close to its expectation with respect to the $l_2$ operator norm. This paper uses semigroup methods to derive sharp nonlinear matrix inequalities. In particular, it is shown that the classic Bakry-\'Emery curvature criterion implies subgaussian concentration for "matrix Lipschitz" functions. This argument circumvents the need to develop a matrix version of the log-Sobolev inequality, a technical obstacle that has blocked previous attempts to derive matrix concentration inequalities in this setting. The approach unifies and extends much of the previous work on matrix concentration. When applied to a product measure, the theory reproduces the matrix Efron-Stein inequalities due to Paulin et al. It also handles matrix-valued functions on a Riemannian manifold with uniformly positive Ricci curvature.


title: Alignment of irregular grains by radiative torques: efficiency study

abstract: We study the efficiency of grain alignment by radiative torques (RATs) for an ensemble of irregular grains. The grains are modeled as ensembles of oblate and prolate spheroids, deformed as Gaussian random ellipsoids, and their scattering interactions are solved using numerically exact methods. We define the fraction of the grains that both rotate fast and demonstrate perfect alignment with grain long axes perpendicular to the magnetic field. We demonstrate that for typical interstellar conditions the degree of alignment arising from the RAT mechanism is significantly larger than that arising from the Davis-Greenstein (DG) mechanism based on paramagnetic relaxation. We quantify a factor related to the efficacy of alignment and show that it is related to a $q_\mathrm{max}$-factor of analytical model (AMO) of the RAT theory. Our results indicate that the RAT alignment can potentially be sufficiently strong and to explain observations even if grains do not have magnetic inclusions.


title: Low--Temperature Plasma Technology: Impact on Indian Rural Life

abstract: The low temperature plasma technology has been emerged as an alternative environmental friendly technology that has a potential to replace the traditional technologies in various sectors such as textile industries, agriculture, food industries and waste water treatment. The future prospective of the low temperature plasma technology on the life of rural Indians has been discussed in this brief report.


title: Erd\H{o}s distinct distances in hyperbolic surfaces

abstract: In this paper, we introduce the notion of "geodesic cover" for Fuchsian groups, which summons copies of fundamental polygons in the hyperbolic plane to cover pairs of representatives realizing distances in the corresponding hyperbolic surface. Then we use estimates of geodesic-covering numbers to study the distinct distances problem in hyperbolic surfaces. Especially, for $Y$ from a large class of hyperbolic surfaces, we establish the nearly optimal bound $\geq c(Y)N/\log N$ for distinct distances determined by any $N$ points in $Y$, where $c(Y)>0$ is some constant depending only on $Y$. In particular, for $Y$ being modular surface or standard regular of genus $g\geq 2$, we evaluate $c(Y)$ explicitly. We also derive new sum-product type estimates.


title: Adaptive and Pressure-Robust Discretization of Incompressible Pressure-Driven Phase-Field Fracture

abstract: In this work, we consider pressurized phase-field fracture problems in nearly and fully incompressible materials. To this end, a mixed form for the solid equations is proposed. To enhance the accuracy of the spatial discretization, a residual-type error estimator is developed. Our algorithmic advancements are substantiated with several numerical tests that are inspired from benchmark configurations. Therein, a primal-based formulation is compared to our newly developed mixed phase-field fracture method for Poisson ratios approaching $\nu \to 0.5$. Finally, for $\nu = 0.5$, we compare the numerical results of the mixed formulation with a pressure robust modification.


title: $\Lambda\Lambda N$-$\Xi NN$ $S$ wave resonance

abstract: We use an existing model of the $\Lambda\Lambda N - \Xi NN$ three-body system based in two-body separable interactions to study the $(I,J^P)=(1/2,1/2^+)$ three-body channel. For the $\Lambda\Lambda$, $\Xi N$, and $\Lambda\Lambda - \Xi N$ amplitudes we have constructed separable potentials based on the most recent results of the HAL QCD Collaboration. They are characterized by the existence of a resonance just below or above the $\Xi N$ threshold in the so-called $H$-dibaryon channel, $(i,j^p)=(0,0^+)$. A three-body resonance appears {2.3} MeV above the $\Xi d$ threshold. We show that if the $\Lambda\Lambda - \Xi N$ $H$-dibaryon channel is not considered, the $\Lambda\Lambda N - \Xi NN$ $S$ wave resonance disappears. Thus, the possible existence of a $\Lambda\Lambda N - \Xi NN$ resonance would be sensitive to the $\Lambda\Lambda - \Xi N$ interaction. The existence or nonexistence of this resonance could be evidenced by measuring, for example, the $\Xi d$ cross section.


title: Ferromagnetic contamination of Ultra-Low-Field-NMR sample containers. Quantification of the problem and possible solutions

abstract: The presence of a weak remanence in Ultra-Low-Field (ULF) NMR sample containers is investigated on the basis of proton precession. The high-sensitivity magnetometer used for the NMR detection, enables simultaneously the measurement of the static field produced in the sample proximity by ferromagnetic contaminants. The presence of the latter is studied by high resolution chemical analyses of the surface, based on X-ray fluorescence spectroscopy and secondary ions mass spectroscopy. Methodologies to reduce the contamination are explored and characterized. This study is of relevance in any ULF-NMR experiment, as in the ULF regime spurious ferromagnetism becomes easily a dominant cause of artefacts.


title: Forced-exploration free Strategies for Unimodal Bandits

abstract: We consider a multi-armed bandit problem specified by a set of Gaussian or Bernoulli distributions endowed with a unimodal structure. Although this problem has been addressed in the literature (Combes and Proutiere, 2014), the state-of-the-art algorithms for such structure make appear a forced-exploration mechanism. We introduce IMED-UB, the first forced-exploration free strategy that exploits the unimodal-structure, by adapting to this setting the Indexed Minimum Empirical Divergence (IMED) strategy introduced by Honda and Takemura (2015). This strategy is proven optimal. We then derive KLUCB-UB, a KLUCB version of IMED-UB, which is also proven optimal. Owing to our proof technique, we are further able to provide a concise finite-time analysis of both strategies in an unified way. Numerical experiments show that both IMED-UB and KLUCB-UB perform similarly in practice and outperform the state-of-the-art algorithms.


title: Visualizing an adjustable WO3/p-GaN heterojunction

abstract: The p-n junctions based on typical semiconductors are the elementary units for the modern electronic devices and chip industry. While the rectification property of those p-n junction is usually fixed once the unit is fabricated. Here, we proposed an adjustable n-WO3/p-GaN heterojunction with controllable electronic properties. For the prepared n-WO3/p-GaN heterojunction, it is almost transparent and shows typical p-n junction rectification. While if gradually doping some hydrogen atoms into WO3 layer by a facile electron-proton synergistic route, the heterojunction can be turned dynamically from the typical p-n junction (n-WO3/p-GaN) to standard Schottky contact (HxWO3/p-GaN) step by step. More importantly, this evolution can be directly visualized by eyesight due to the pronounced electrochromic characteristic of WO3 layer. By connecting two HxWO3/p-GaN heterojunctions, the controllable bi-functional rectification can be achieved. In addition, the HxWO3/p-GaN heterojunction can recovered to the original p-n jucntion just by annealing at ambient, demonstrating the heterojunction is controllable and reusable. The current study will open up tremendous opportunities for dynamic electronic devices in the future.


title: Tackling Occlusion in Siamese Tracking with Structured Dropouts

abstract: Occlusion is one of the most difficult challenges in object tracking to model. This is because unlike other challenges, where data augmentation can be of help, occlusion is hard to simulate as the occluding object can be anything in any shape. In this paper, we propose a simple solution to simulate the effects of occlusion in the latent space. Specifically, we present structured dropout to mimick the change in latent codes under occlusion. We present three forms of dropout (channel dropout, segment dropout and slice dropout) with the various forms of occlusion in mind. To demonstrate its effectiveness, the dropouts are incorporated into two modern Siamese trackers (SiamFC and SiamRPN++). The outputs from multiple dropouts are combined using an encoder network to obtain the final prediction. Experiments on several tracking benchmarks show the benefits of structured dropouts, while due to their simplicity requiring only small changes to the existing tracker models.


title: Mitigating undesirable emergent behavior arising between driver and semi-automated vehicle

abstract: Emergent behavior arising in a joint human-robot system cannot be fully predicted based on an understanding of the individual agents. Typically, robot behavior is governed by algorithms that optimize a reward function that should quantitatively capture the joint system's goal. Although reward functions can be updated to better match human needs, this is no guarantee that no misalignment with the complex and variable human needs will occur. Algorithms may learn undesirable behavior when interacting with the human and the intrinsically unpredictable human-inhabited world, thereby producing further misalignment with human users or bystanders. As a result, humans might behave differently than anticipated, causing robots to learn differently and undesirable behavior to emerge. With this short paper, we state that to design for Human-Robot Interaction that mitigates such undesirable emergent behavior, we need to complement advancements in human-robot interaction algorithms with human factors knowledge and expertise. More specifically, we advocate a three-pronged approach that we illustrate using a particularly challenging example of safety-critical human-robot interaction: a driver interacting with a semi-automated vehicle. Undesirable emergent behavior should be mitigated by a combination of 1) including driver behavioral mechanisms in the vehicle's algorithms and reward functions, 2) model-based approaches that account for interaction-induced driver behavioral adaptations and 3) driver-centered interaction design that promotes driver engagement with the semi-automated vehicle, and the transparent communication of each agent's actions that allows mutual support and adaptation. We provide examples from recent empirical work in our group, in the hope this proves to be fruitful for discussing emergent human-robot interaction.


title: Subspace approximation with outliers

abstract: The subspace approximation problem with outliers, for given $n$ points in $d$ dimensions $x_{1},\ldots, x_{n} \in R^{d}$, an integer $1 \leq k \leq d$, and an outlier parameter $0 \leq \alpha \leq 1$, is to find a $k$-dimensional linear subspace of $R^{d}$ that minimizes the sum of squared distances to its nearest $(1-\alpha)n$ points. More generally, the $\ell_{p}$ subspace approximation problem with outliers minimizes the sum of $p$-th powers of distances instead of the sum of squared distances. Even the case of robust PCA is non-trivial, and previous work requires additional assumptions on the input. Any multiplicative approximation algorithm for the subspace approximation problem with outliers must solve the robust subspace recovery problem, a special case in which the $(1-\alpha)n$ inliers in the optimal solution are promised to lie exactly on a $k$-dimensional linear subspace. However, robust subspace recovery is Small Set Expansion (SSE)-hard. We show how to extend dimension reduction techniques and bi-criteria approximations based on sampling to the problem of subspace approximation with outliers. To get around the SSE-hardness of robust subspace recovery, we assume that the squared distance error of the optimal $k$-dimensional subspace summed over the optimal $(1-\alpha)n$ inliers is at least $\delta$ times its squared-error summed over all $n$ points, for some $0 < \delta \leq 1 - \alpha$. With this assumption, we give an efficient algorithm to find a subset of $poly(k/\epsilon) \log(1/\delta) \log\log(1/\delta)$ points whose span contains a $k$-dimensional subspace that gives a multiplicative $(1+\epsilon)$-approximation to the optimal solution. The running time of our algorithm is linear in $n$ and $d$. Interestingly, our results hold even when the fraction of outliers $\alpha$ is large, as long as the obvious condition $0 < \delta \leq 1 - \alpha$ is satisfied.


title: Ergodic property for Galton-Watson processes in which individuals have variable life times

abstract: This paper is concerned with an extended Galton-Watson process so as to allow individuals to live and reproduce for more than one unit time. We assume that each individual can live $k$ seasons (time-units) with probability $h_k$, and produce $m$ offspring with probability $p_m$ during each season. These can be seen as Galton-Watson processes with countably infinitely many types in which particles of type $i$ may only have offspring of type $i+1$ and type $1$. Let $\textbf{M}$ be its mean progeny matrix and $\gamma$ be the convergence radius of the power series $\sum_{k\geq 0}r^k(\textbf{M}^k)_{ij}$. We first derive formula of calculating $\gamma$ and show that $\gamma$, in supercritical case, is actually the extinction probability of a Galton-Watson process. Next, we give clear criteria for $\textbf{M}$ to be $\gamma$-transient, $\gamma$-positive and $\gamma$-null recurrent from which the ergodic property of the process is discussed. The criteria for $\gamma$ and $\gamma$-recurrence of $\textbf{M}$ rely on the properties of lifetime distribution which are easier to be verified than current results. Finally, we show the asymptotic behavior of the total population size of each type of individuals under certain conditions which illustrates the evolution of Galton-Watson process in which individuals have variable lifetimes.


title: Supersymmetric spin-phonon coupling prevents odd integer spins from quantum tunneling

abstract: Quantum tunneling of the magnetization is a phenomenon that impedes the use of small anisotropic spin systems for storage purposes even at the lowest temperatures. Phonons, usually considered for temperature dependent relaxation of magnetization over the anisotropy barrier, also contribute to magnetization tunneling for integer spin quantum numbers. Here we demonstrate that certain spin-phonon Hamiltonians are unexpectedly robust against the opening of a tunneling gap, even for strong spin-phonon coupling. The key to understanding this phenomenon is provided by an underlying supersymmetry that involves both spin and phonon degrees of freedom.


title: Higher order Levi forms on homogeneous CR manifolds

abstract: We investigate the nondegeneracy of higher order Levi forms on weakly nondegenerate homogeneous $CR$ manifolds. Improving previous results, we prove that general orbits of real forms in complex flag manifolds have order less or equal $3$ and the compact ones less or equal~$2$. Finally we construct by Lee extensions weakly nondegenerate $CR$ vector bundles with arbitrary orders of nondegeneracy.


title: On the Maxey-Riley equation of motion and its extension to high Reynolds numbers

abstract: The inertial response of a particle to turbulent flows is a problem of relevance to a wide range of environmental and engineering problems. The equation most often used to describe the force balance is the Maxey-Riley equation, which includes in addition to buoyancy and steady drag forces, an unsteady Basset drag force related to past particle acceleration. Here we provide a historical review of how the Maxey-Riley equation was developed and how it is only suited for studies where the Reynolds number is less than unity. Revisiting the innovative mathematical methods employed by Basset (1888), we introduce an alternative formulation for the unsteady drag for application to a broader range of particle motions. While the Basset unsteady drag is negligible at higher Reynolds numbers, the revised unsteady drag is not.


title: Accelerating Binarized Neural Networks via Bit-Tensor-Cores in Turing GPUs

abstract: Despite foreseeing tremendous speedups over conventional deep neural networks, the performance advantage of binarized neural networks (BNNs) has merely been showcased on general-purpose processors such as CPUs and GPUs. In fact, due to being unable to leverage bit-level-parallelism with a word-based architecture, GPUs have been criticized for extremely low utilization (1%) when executing BNNs. Consequently, the latest tensorcores in NVIDIA Turing GPUs start to experimentally support bit computation. In this work, we look into this brand new bit computation capability and characterize its unique features. We show that the stride of memory access can significantly affect performance delivery and a data-format co-design is highly desired to support the tensorcores for achieving superior performance than existing software solutions without tensorcores. We realize the tensorcore-accelerated BNN design, particularly the major functions for fully-connect and convolution layers -- bit matrix multiplication and bit convolution. Evaluations on two NVIDIA Turing GPUs show that, with ResNet-18, our BTC-BNN design can process ImageNet at a rate of 5.6K images per second, 77% faster than state-of-the-art. Our BNN approach is released on https://github.com/pnnl/TCBNN.


title: Counting independent sets in Riordan graphs

abstract: The notion of a Riordan graph was introduced recently, and it is a far-reaching generalization of the well-known Pascal graphs and Toeplitz graphs. However, apart from a certain subclass of Toeplitz graphs, nothing was known on independent sets in Riordan graphs. In this paper, we give exact enumeration and lower and upper bounds for the number of independent sets for various classes of Riordan graphs. Remarkably, we offer a variety of methods to solve the problems that range from the structural decomposition theorem to methods in combinatorics on words. Some of our results are valid for any graph.


title: Baskets and fibred links realizing $A_{n}$

abstract: We prove that basket links, whose symmetrized Seifert form is congruent to the Cartan matrix of the simply laced Dynkin diagram $A_{n}$, are isotopic to the torus link $T(2,n+1)$. In addition, we provide examples of links, constructed by plumbing $n$ positive Hopf bands the core curves of which intersect at most once, with symmetrized Seifert form congruent to the Cartan matrix $A_{n}$, that are not isotopic to $T(2,n+1)$.


title: Early Exit or Not: Resource-Efficient Blind Quality Enhancement for Compressed Images

abstract: Lossy image compression is pervasively conducted to save communication bandwidth, resulting in undesirable compression artifacts. Recently, extensive approaches have been proposed to reduce image compression artifacts at the decoder side; however, they require a series of architecture-identical models to process images with different quality, which are inefficient and resource-consuming. Besides, it is common in practice that compressed images are with unknown quality and it is intractable for existing approaches to select a suitable model for blind quality enhancement. In this paper, we propose a resource-efficient blind quality enhancement (RBQE) approach for compressed images. Specifically, our approach blindly and progressively enhances the quality of compressed images through a dynamic deep neural network (DNN), in which an early-exit strategy is embedded. Then, our approach can automatically decide to terminate or continue enhancement according to the assessed quality of enhanced images. Consequently, slight artifacts can be removed in a simpler and faster process, while the severe artifacts can be further removed in a more elaborate process. Extensive experiments demonstrate that our RBQE approach achieves state-of-the-art performance in terms of both blind quality enhancement and resource efficiency. The code is available at https://github.com/RyanXingQL/RBQE.


title: Spin and density self-ordering in dynamic polarization gradients fields

abstract: We study the zero-temperature quantum phase diagram for a two-component Bose-Einstein condensate in an optical cavity. The two atomic spin states are Raman coupled by two transverse orthogonally-polarized, blue detuned plane-wave lasers inducing a repulsive cavity potential. For weak pump the lasers favor a state with homogeneous density and predefined uniform spin direction. When one pump laser is polarized parallel to the cavity mode polarization, the photons coherently scattered into the resonator induce a polarization gradient along the cavity axis, which mediates long-range density-density, spin-density, and spin-spin interactions. We show that the coupled atom-cavity system implements central aspects of the $t$-$J$-$V$-$W$ model with a rich phase diagram. At the mean-field limit we identify at least four qualitatively distinct density- and spin-ordered phases including ferro- and anti-ferromagnetic order along the cavity axis, which can be controlled via the pump strength and detuning. A real time observation of amplitude and phase of the emitted fields bears strong signatures of the realized phase and allows for real-time determination of phase transition lines. Together with measurements of the population imbalance most properties of the phase diagram can be reconstructed.


title: Pan-Sharpening with Color-Aware Perceptual Loss and Guided Re-Colorization

abstract: We present a novel color-aware perceptual (CAP) loss for learning the task of pan-sharpening. Our CAP loss is designed to focus on the deep features of a pre-trained VGG network that are more sensitive to spatial details and ignore color information to allow the network to extract the structural information from the PAN image while keeping the color from the lower resolution MS image. Additionally, we propose "guided re-colorization", which generates a pan-sharpened image with real colors from the MS input by "picking" the closest MS pixel color for each pan-sharpened pixel, as a human operator would do in manual colorization. Such a re-colorized (RC) image is completely aligned with the pan-sharpened (PS) network output and can be used as a self-supervision signal during training, or to enhance the colors in the PS image during test. We present several experiments where our network trained with our CAP loss generates naturally looking pan-sharpened images with fewer artifacts and outperforms the state-of-the-arts on the WorldView3 dataset in terms of ERGAS, SCC, and QNR metrics.


title: Survey of Gravitationally lensed Objects in HSC Imaging (SuGOHI) -- VII. Discovery and Confirmation of Three Strongly Lensed Quasars

abstract: We present spectroscopic confirmation of three new two-image gravitationally lensed quasars, compiled from existing strong lens and X-ray catalogs. Images of HSC J091843.27$-$022007.5 show a red galaxy with two blue point sources at either side, separated by 2.26 arcsec. This system has a source and a lens redshifts $z_s=0.804$ and $z_{\ell}=0.459$, respectively, as obtained by our follow-up spectroscopic data. CXCO J100201.50$+$020330.0 shows two point sources separated by 0.85 arcsec on either side of an early-type galaxy. The follow-up spectroscopic data confirm the fainter quasar has the same redshift with the brighter quasar from the SDSS fiber spectrum at $z_s=2.016$. The deflecting foreground galaxy is a typical early-type galaxy at a redshift of $z_{\ell}=0.439$. SDSS J135944.21$+$012809.8 has two point sources with quasar spectra at the same redshift $z_s=1.096$, separated by 1.05 arcsec, and fits to the HSC images confirm the presence of a galaxy between these. These discoveries demonstrate the power of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)'s deep imaging and wide sky coverage. Combined with existing X-ray source catalogues and follow-up spectroscopy, the HSC-SSP provides us unique opportunities to find multiple-image quasars lensed by a foreground galaxy.


title: Territories of Parrondo's paradox and its relation with entanglement in quantum walks

abstract: Parrondo's paradox is a well-known counterintuitive phenomenon, where the combination of unfavorable situations can establish favorable ones. In this paper, we study one-dimensional discrete-time quantum walks, manipulating two different coins (two-state) operators representing two losing games A and B, respectively, to create the Parrondo effect in the quantum domain. We exhibit that games A and B are losing games when played individually but could produce a winning expectation when played alternatively for a particular sequence of different periods for distinct choices of the relative phase. Furthermore, we investigate the regimes of the relative phase of the initial state of coins where Parrondo games exist. Moreover, we also analyze the relationships between Parrondo's game and quantum entanglement and show regimes where the Parrondo sequence may generate a maximal entangler state in our scheme. Along with the applications of different kinds of quantum walks, our outcomes potentially encourage the development of new quantum algorithms.


title: Engineering light absorption at critical coupling via bound states in the continuum

abstract: Recent progress in nanophotonics is driven by the desire to engineer light-matter interaction in two-dimensional (2D) materials using high-quality resonances in plasmonic and dielectric structures. Here, we demonstrate a link between the radiation control at critical coupling and the metasurface-based bound states in the continuum (BIC) physics, and develop a generalized theory to engineer light absorption of 2D materials in coupling resonance metasurfaces. In a typical example of hybrid graphene-dielectric metasurfaces, we present the manipulation of absorption bandwidth by more than one order of magnitude by simultaneously adjusting the asymmetry parameter of silicon resonators governed by BIC and the graphene surface conductivity while the absorption efficiency maintains maximum. This work reveals the generalized role of BIC in the radiation control at critical coupling and provides promising strategies in engineering light absorption of 2D materials for high-efficiency optoelectronics device applications, e.g., light emission, detection and modulation.


title: Growth mechanism of polymer membranes obtained by H-bonding across immiscible liquid interfaces

abstract: Complexation of polymers at liquid interfaces is an emerging technique to produce all-liquid printable and self-healing devices and membranes. It is crucial to control the assembly process but the mechanisms at play remain unclear. Using two different reflectometric methods, we investigate the spontaneous growth of H-bonded PPO-PMAA membranes at a flat liquid-liquid interface. We find that the membrane thickness h grows with time t as h~t^(1/2), which is reminiscent of a diffusion-limited process. However, counter-intuitively, we observe that this process is faster as the PPO molar mass increases. We are able to rationalize these results with a model which considers the diffusion of the PPO chains within the growing membrane. The architecture of the latter is described as a gel-like porous network, with a pore size much smaller than the radius of the diffusing PPO chains, thus inducing entropic barriers that hinder the diffusion process. From the comparison between the experimental data and the result of the model, we extract some key piece of information about the microscopic structure of the membrane. This study opens the route toward the rational design of self-assembled membranes and capsules with optimal properties.


title: Minimal phase-field crystal modeling of vapor-liquid-solid coexistence and transitions

abstract: A phase-field crystal model based on the density-field approach incorporating high-order interparticle direct correlations is developed to study vapor-liquid-solid coexistence and transitions within a single continuum description. Conditions for the realization of the phase coexistence and transition sequence are systematically analyzed and shown to be satisfied by a broad range of model parameters, demonstrating the high flexibility and applicability of the model. Both temperature-density and temperature-pressure phase diagrams are identified, while structural evolution and coexistence among the three phases are examined through dynamical simulations. The model is also able to produce some temperature and pressure related material properties, including effects of thermal expansion and pressure on equilibrium lattice spacing, and temperature dependence of saturation vapor pressure. This model can be used as an effective approach for investigating a variety of material growth and deposition processes based on vapor-solid, liquid-solid, and vapor-liquid-solid growth.


title: On the Demystification of Knowledge Distillation: A Residual Network Perspective

abstract: Knowledge distillation (KD) is generally considered as a technique for performing model compression and learned-label smoothing. However, in this paper, we study and investigate the KD approach from a new perspective: we study its efficacy in training a deeper network without any residual connections. We find that in most of the cases, non-residual student networks perform equally or better than their residual versions trained on raw data without KD (baseline network). Surprisingly, in some cases, they surpass the accuracy of baseline networks even with the inferior teachers. After a certain depth of non-residual student network, the accuracy drop, coming from the removal of residual connections, is substantial, and training with KD boosts the accuracy of the student up to a great extent; however, it does not fully recover the accuracy drop. Furthermore, we observe that the conventional teacher-student view of KD is incomplete and does not adequately explain our findings. We propose a novel interpretation of KD with the Trainee-Mentor hypothesis, which provides a holistic view of KD. We also present two viewpoints, loss landscape, and feature reuse, to explain the interplay between residual connections and KD. We substantiate our claims through extensive experiments on residual networks.


title: Robust Kernel Density Estimation with Median-of-Means principle

abstract: In this paper, we introduce a robust nonparametric density estimator combining the popular Kernel Density Estimation method and the Median-of-Means principle (MoM-KDE). This estimator is shown to achieve robustness to any kind of anomalous data, even in the case of adversarial contamination. In particular, while previous works only prove consistency results under known contamination model, this work provides finite-sample high-probability error-bounds without a priori knowledge on the outliers. Finally, when compared with other robust kernel estimators, we show that MoM-KDE achieves competitive results while having significant lower computational complexity.


title: A Novel Bistatic Joint Radar-Communication System in Multi-path Environments

abstract: Radar detection and communication can be operated simultaneously in joint radar-communication (JRC) system. In this paper, we propose a bistatic JRC system which is applicable in multi-path environments. Basing on a novel joint waveform, a joint detection process is designed for both target detection and channel estimation. Meanwhile, a low-cost channel equalization method that utilizes the channel state information acquired from the detection process is proposed. The numerical results show that the symbol error rate (SER) of the proposed system is similar to that of the binary frequency shift keying system, and the signal to noise ratio requirement in multi-path environments is less than 2 dB higher compared with that in single-path environment to reach a SER of 10-5. Besides, the knowledge of the embedded information is not required for the joint detection process and the detection performance is robust to unknown information.


title: Comment on "Measurability of nonequilibrium thermodynamics in terms of the Hamiltonian of mean force"

abstract: In a recent paper [P. Strasberg and M. Esposito, Phys. Rev. E {\bf 101}, 050101(R) (2020)] an attempt is presented to formulate the nonequilibrium thermodynamics of an open system in terms of the Hamiltonian of mean force. The purpose of the present comment is to clarify severe restrictions of this approach and also to stress that recently noted ambiguities [P. Talkner and P. H\"anggi, Phys. Rev. E {\bf 94}, 022143 (2016)] of fluctuating thermodynamic potentials cannot be removed in the suggested way.


title: QCD sum rules with spectral densities solved in inverse problems

abstract: We construct QCD sum rules for nonperturbative studies without assuming the quark-hadron duality for the spectral density at low energy on the hadron side. Instead, both resonance and continuum contributions to the spectral density are solved with the operator-product-expansion input on the quark side by treating sum rules as an inverse problem. This new formalism does not involve the continuum threshold, does not require the Borel transformation and stability analysis, and can be extended to extract properties of excited states. Taking the two-current correlator as an example, we demonstrate that the series of $\rho$ resonances can emerge in our formalism, and the decay constants $f_{\rho(770)} (f_{\rho(1450)},f_{\rho(1700)},f_{\rho(1900)})\approx$ 0.22 (0.19, 0.14, 0.14) GeV for the masses $m_{\rho(770)} (m_{\rho(1450)},m_{\rho(1700)},m_{\rho(1900)})\approx$ 0.78 (1.46, 1.70, 1.90) GeV are determined. We also show that the decay width $\Gamma_{\rho(770)}\approx 0.17$ GeV can be obtained by substituting a Breit-Wigner parametrization for the $\rho(770)$ pole on the hadron side. It is observed that quark condensates of dimension-six on the quark side are crucial for establishing those $\rho$ resonances. Handling the conventional sum rules with the duality assumption as an inverse problem, we find that the multiple pole sum rules widely adopted in the literature do not describe the $\rho$ excitations reasonably. The precision of our theoretical outcomes can be improved systematically by including higher-order and higher-power corrections on the quark side. Broad applications of this formalism to abundant low energy QCD observables are expected.


title: Neutrino amplitude decomposition: Toward observing the atmospheric - solar wave interference

abstract: Observation of the interference between the atmospheric and solar oscillation waves with the correct magnitude would provide another manifestation of the three-generation structure of leptons. As a prerequisite for such analyses we develop a method for decomposing the oscillation $S$ matrix into the atmospheric and solar amplitudes. Though the similar method was recently proposed successfully in vacuum, once an extension into the matter environment is attempted, it poses highly nontrivial problems. Even for an infinitesimal matter potential, inherent mixture of the atmospheric and solar oscillation waves occurs, rendering a simple extension of the vacuum definition untenable. We utilize general kinematic structure as well as analyses of the five perturbative frameworks, in which the nature of matter-dressed atmospheric and solar oscillations are known, to understand the origin of the trouble, how to deal with the difficulty, and to grasp the principle of decomposition. Then, we derive the amplitude decomposition formulas in these frameworks, and discuss properties of the decomposed probabilities. We mostly discuss the $\nu_{\mu} \rightarrow \nu_{e}$ channel, but a comparison with the $\nu_{\mu} \rightarrow \nu_{\tau}$ channel reveals an interesting difference.


title: Polynomial stabilization of non-smooth direct/indirect elastic/viscoelastic damping problem involving Bresse system

abstract: We consider an elastic/viscoelastic transmission problem for the Bresse system with fully Dirichlet or Dirichlet-Neumann-Neumann boundary conditions. The physical model consists of three wave equations coupled in certain pattern. The system is damped directly or indirectly by global or local Kelvin-Voigt damping. Actually, the number of the dampings, their nature of distribution (locally or globally) and the smoothness of the damping coefficient at the interface play a crucial role in the type of the stabilization of the corresponding semigroup. Indeed, using frequency domain approach combined with multiplier techniques and the construction of a new multiplier function, we establish different types of energy decay rate (see the table of stability results below). Our results generalize and improve many earlier ones in the literature and in particular some studies done on the Timoshenko system with Kelvin-Voigt damping.


title: On the Potential of Dynamic Substructuring Methods for Model Updating

abstract: While purely data-driven assessment is feasible for the first levels of the Structural Health Monitoring (SHM) process, namely damage detection and arguably damage localization, this does not hold true for more advanced processes. The tasks of damage quantification and eventually residual life prognosis are invariably linked to availability of a representation of the system, which bears physical connotation. In this context, it is often desirable to assimilate data and models, into what is often termed a digital twin of the monitored system. One common take to such an end lies in exploitation of structural mechanics models, relying on use of Finite Element approximations. proper updating of these models, and their incorporation in an inverse problem setting may allow for damage quantification and localization, as well as more advanced tasks, including reliability analysis and fatigue assessment. However, this may only be achieved by means of repetitive analyses of the forward model, which implies considerable computational toll, when the model used is a detailed FE representation. In tackling this issue, reduced order models can be adopted, which retain the parameterisation and link to the parameters regulating the physical properties, albeit greatly reducing the computational burden. In this work a detailed FE model of a wind turbine tower is considered, reduced forms of this model are found using both the Craig Bampton and Dual Craig Bampton methods. These reduced order models are then used and compared in a Transitional Markov Chain Monte Carlo procedure to localise and quantify damage which is introduced to the system.


title: Regression with reject option and application to kNN

abstract: We investigate the problem of regression where one is allowed to abstain from predicting. We refer to this framework as regression with reject option as an extension of classification with reject option. In this context, we focus on the case where the rejection rate is fixed and derive the optimal rule which relies on thresholding the conditional variance function. We provide a semi-supervised estimation procedure of the optimal rule involving two datasets: a first labeled dataset is used to estimate both regression function and conditional variance function while a second unlabeled dataset is exploited to calibrate the desired rejection rate. The resulting predictor with reject option is shown to be almost as good as the optimal predictor with reject option both in terms of risk and rejection rate. We additionally apply our methodology with kNN algorithm and establish rates of convergence for the resulting kNN predictor under mild conditions. Finally, a numerical study is performed to illustrate the benefit of using the proposed procedure.


title: Accretion bursts in low-metallicity protostellar disks

abstract: The early evolution of protostellar disks with metallicities in the $Z=1.0-0.01~Z_\odot$ range was studied with a particular emphasis on the strength of gravitational instability and the nature of protostellar accretion in low-metallicity systems. Numerical hydrodynamics simulations in the thin-disk limit were employed that feature separate gas and dust temperatures, and disk mass-loading from the infalling parental cloud cores. Models with cloud cores of similar initial mass and rotation pattern, but distinct metallicity were considered to distinguish the effect of metallicity from that of initial conditions. The early stages of disk evolution in low-metallicity models are characterized by vigorous gravitational instability and fragmentation. Disk instability is sustained by continual mass-loading from the collapsing core. The time period that is covered by this unstable stage is much shorter in the $Z=0.01~Z_\odot$ models as compared to their higher metallicity counterparts thanks to the higher mass infall rates caused by higher gas temperatures (that decouple from lower dust temperatures) in the inner parts of collapsing cores. Protostellar accretion rates are highly variable in the low-metallicity models reflecting a highly dynamical nature of the corresponding protostellar disks. The low-metallicity systems feature short, but energetic episodes of mass accretion caused by infall of inward-migrating gaseous clumps that form via gravitational fragmentation of protostellar disks. These bursts seem to be more numerous and last longer in the $Z=0.1~Z_\odot$ models in comparison to the $Z=0.01~Z_\odot$ case. Variable protostellar accretion with episodic bursts is not a particular feature of solar metallicity disks. It is also inherent to gravitationally unstable disks with metallicities up to 100 times lower than solar.


title: Terahertz Pulse Shaping Using Diffractive Surfaces

abstract: Recent advances in deep learning have been providing non-intuitive solutions to various inverse problems in optics. At the intersection of machine learning and optics, diffractive networks merge wave-optics with deep learning to design task-specific elements to all-optically perform various tasks such as object classification and machine vision. Here, we present a diffractive network, which is used to shape an arbitrary broadband pulse into a desired optical waveform, forming a compact pulse engineering system. We experimentally demonstrate the synthesis of square pulses with different temporal-widths by manufacturing passive diffractive layers that collectively control both the spectral amplitude and the phase of an input terahertz pulse. Our results constitute the first demonstration of direct pulse shaping in terahertz spectrum, where a complex-valued spectral modulation function directly acts on terahertz frequencies. Furthermore, a Lego-like physical transfer learning approach is presented to illustrate pulse-width tunability by replacing part of an existing network with newly trained diffractive layers, demonstrating its modularity. This learning-based diffractive pulse engineering framework can find broad applications in e.g., communications, ultra-fast imaging and spectroscopy.


title: Exponential inequalities for sampling designs

abstract: In this work we introduce a general approach, based on the mar-tingale representation of a sampling design and Azuma-Hoeffding's inequality , to derive exponential inequalities for the difference between a Horvitz-Thompson estimator and its expectation. Applying this idea, we establish such inequalities for Chao's procedure, Till{\'e}'s elimination procedure, the generalized Midzuno method as well as for Brewer's method. As a by-product, we prove that the first three sampling designs are (conditionally) negatively associated. For such sampling designs, we show that that the inequality we obtain is usually sharper than the one obtained by applying known results for negatively associated random variables.


title: Symbolic Execution and Debugging Synchronization

abstract: In this thesis, we introduce the idea of combining symbolic execution with dynamic analysis for reverse engineering. Differently from DSE, we devise an approach where the reverse engineer can use a debugger to drive and inspect a concrete execution engine of the application code and then, when needed, transfer the execution into a symbolic executor in order to automatically identify the input values required to reach a target point in the code. After that, the user can also transfer back the correct input values found with symbolic execution in order to continue the debugging. The synchronization between a debugger and a symbolic executor can enhance manual dynamic analysis and allow a reverser to easily solve small portions of code without leaving the debugger. We implemented a synchronization mechanism on top of the binary analysis framework angr, allowing for transferring the state of the debugged process to the angr environment and back. The backend library is debugger agnostic and can be extended to work with various frontends. We implemented a frontend for the IDA Pro debugger and one for the GNU Debugger, which are both widely popular among reverse engineers.


title: Denjoy sub-systems and horseshoes

abstract: We introduce a notion of weak Denjoy subsystem (WDS) that generalizes the Aubry-Mather Cantor sets to diffeomorphisms of manifolds. We explain how a rotation number can be associated to such a WDS. Then we build in any horseshoe a continuous one parameter family of such WDS that is indexed by its rotation number. Looking at the inverse problem in the setting of Aubry-Mather theory, we also prove that for a generic conservative twist map of the annulus, the majority of the Aubry-Mather sets are contained in some horseshoe that is associated to a Aubry-Mather set with a rational rotation number.


title: Synchrotron radiation from ultrahigh-intensity laser-plasma interactions and competition with Bremsstrahlung in thin foil targets

abstract: By means of particle-in-cell numerical simulations, we investigate the emission of high-energy photons in laser-plasma interactions under ultrahigh-intensity conditions relevant to multi-petawatt laser systems. We first examine the characteristics of synchrotron radiation from laser-driven plasmas of varying density and size. In particular, we show and explain the dependence of the angular distribution of the radiated photons on the transparency or opacity of the plasma. We then study the competition of the synchrotron and Bremsstrahlung emissions in copper foil targets irradiated by $10^{22}\,\rm W\,cm^{-2}$, $50 \, \rm fs$ laser pulses. Synchrotron emission is observed to be maximized for target thicknesses of a few $10 \, \rm nm$, close to the relativistic transparency threshold, and to be superseded by Bremsstrahlung in targets a few $\mu$m thick. At their best efficiency, both mechanisms are found to radiate about one percent of the laser energy into photons with energies above $10\,\rm keV$. Their energy and angular spectra are thoroughly analyzed in light of the ultrafast target dynamics.


title: Ultrasonic Actuation of a Fine-Needle Improves Biopsy Yield

abstract: Despite the ubiquitous use over the past 150 years, the functions of the current medical needle are facilitated only by mechanical shear and cutting by the needle tip,i.e.the lancet. In this study, we demonstrate how nonlinear ultrasonics (NLU) extends the functionality of the medical needle far beyond its present capability. The NLU actions were found to be localized to the proximity of the needle tip, the SonoLancet, but the effects extend several millimeters from the physical needle boundary. The observed nonlinear phenomena, transient cavitation, fluid streams, translation of micro- and nanoparticles and atomization, were quantitatively characterized. In the fine-needle biopsy application, the SonoLancet contributed to obtaining tissue cores with increase in tissue yield by 3-6x in different tissue types compared to conventional needle biopsy technique using the same 21G needle. In conclusion, the SonoLancet could be of interest to several other medical applications, including drug or gene delivery, cell modulation, and minimally invasive surgical procedures.


title: Conformal mapping of the Borel plane: going beyond perturbative QCD

abstract: The power corrections in the Operator Product Expansion (OPE) of QCD correlators can be viewed mathematically as an illustration of the transseries concept, which allows to recover a function from its asymptotic divergent expansion. Alternatively, starting from the divergent behavior of the perturbative QCD encoded in the singularities in the Borel plane, a modified expansion can be defined by means of the conformal mapping of this plane. A comparison of the two approaches concerning their ability to recover nonperturbative properties of the true correlator was not explored up to now. In the present paper, we make a first attempt to investigate this problem. We use for illustration the Adler function and observables expressed as integrals of this function along contours in the complex energy plane. We show that the expansions based on the conformal mapping of the Borel plane go beyond finite-order perturbation theory, containing an infinite number of terms when reexpanded in powers of the coupling. Moreover, the expansion functions exhibit nonperturbative features of the true function, while the expansions have a tamed behavior at large orders and are expected even to be convergent. Using these properties, we argue that there are no mathematical reasons for supplementing the expansions based on the conformal mapping of the Borel plane by additional arbitrary power corrections. Therefore, we make the conjecture that they provide an alternative to the standard OPE in approximating the QCD correlator. This conjecture allows to slightly improve the accuracy of the strong coupling extracted from the hadronic $\tau$ decay width. Using the optimal expansions based on conformal mapping and the contour-improved prescription of renormalization-group resummation, we obtain $\alpha_s(m_\tau^2)=0.314 \pm 0.006$, which implies $\alpha_s(m_Z^2)=0.1179 \pm 0.0008$.


title: Spatiotemporal Multi-Resolution Approximations for Analyzing Global Environmental Data

abstract: Technological developments and open data policies have made large, global environmental datasets accessible to everyone. For analysing such datasets, including spatiotemporal correlations using traditional models based on Gaussian processes does not scale with data volume and requires strong assumptions about stationarity, separability, and distance measures of covariance functions that are often unrealistic for global data. Only very few modeling approaches suitably model spatiotemporal correlations while addressing both computational scalability as well as flexible covariance models. In this paper, we provide an extension to the multi-resolution approximation (MRA) approach for spatiotemporal modeling of global datasets. MRA has been shown to be computationally scalable in distributed computing environments and allows for integrating arbitrary user-defined covariance functions. Our extension adds a spatiotemporal partitioning, and fitting of complex covariance models including nonstationarity with kernel convolutions and spherical distances. We evaluate the effect of the MRA parameters on estimation and spatiotemporal prediction using simulated data, where computation times reduced around two orders of magnitude with an increase of the root-mean-square prediction error of around five percent. This allows for trading off computation times against prediction errors, and we derive a practical strategy for selecting the MRA parameters. We demonstrate how the approach can be practically used for analyzing daily sea surface temperature and precipitation data on global scale and compare models with different complexities in the covariance function.


title: A Framework for Learning Invariant Physical Relations in Multimodal Sensory Processing

abstract: Perceptual learning enables humans to recognize and represent stimuli invariant to various transformations and build a consistent representation of the self and physical world. Such representations preserve the invariant physical relations among the multiple perceived sensory cues. This work is an attempt to exploit these principles in an engineered system. We design a novel neural network architecture capable of learning, in an unsupervised manner, relations among multiple sensory cues. The system combines computational principles, such as competition, cooperation, and correlation, in a neurally plausible computational substrate. It achieves that through a parallel and distributed processing architecture in which the relations among the multiple sensory quantities are extracted from time-sequenced data. We describe the core system functionality when learning arbitrary non-linear relations in low-dimensional sensory data. Here, an initial benefit rises from the fact that such a network can be engineered in a relatively straightforward way without prior information about the sensors and their interactions. Moreover, alleviating the need for tedious modelling and parametrization, the network converges to a consistent description of any arbitrary high-dimensional multisensory setup. We demonstrate this through a real-world learning problem, where, from standard RGB camera frames, the network learns the relations between physical quantities such as light intensity, spatial gradient, and optical flow, describing a visual scene. Overall, the benefits of such a framework lie in the capability to learn non-linear pairwise relations among sensory streams in an architecture that is stable under noise and missing sensor input.


title: The CARMENES search for exoplanets around M dwarfs: Convective shift and starspot constraints from chromatic radial velocities

abstract: Context. Variability caused by stellar activity represents a challenge to the discovery and characterization of terrestrial exoplanets and complicates the interpretation of atmospheric planetary signals. Aims. We aim to use a detailed modeling tool to reproduce the effect of active regions on radial velocity measurements, which aids the identification of the key parameters that have an impact on the induced variability. Methods. We analyzed the effect of stellar activity on radial velocities as a function of wavelength by simulating the impact of the properties of spots, shifts induced by convective motions, and rotation. We focused our modeling effort on the active star YZ CMi (GJ 285), which was photometrically and spectroscopically monitored with CARMENES and the Telescopi Joan Or\'o. Results. We demonstrate that radial velocity curves at different wavelengths yield determinations of key properties of active regions, including spot filling factor, temperature contrast, and location, thus solving the degeneracy between them. Most notably, our model is also sensitive to convective motions. Results indicate a reduced convective shift for M dwarfs when compared to solar-type stars (in agreement with theoretical extrapolations) and points to a small global convective redshift instead of blueshift. Conclusions. Using a novel approach based on simultaneous chromatic radial velocities and light curves, we can set strong constraints on stellar activity, including an elusive parameter such as the net convective motion effect.


title: Calibrating the surface brightness - color relation for late-type red giants stars in the visible domain using VEGA/CHARA interferometric observations

abstract: The surface brightness - color relationship (SBCR) is a poweful tool for determining the angular diameter of stars from photometry. It was for instance used to derive the distance of eclipsing binaries in the Large Magellanic Cloud (LMC), which led to its distance determination with an accuracy of 1%. We calibrate the SBCR for red giant stars in the 2.1 < V-K < 2.5 color range using homogeneous VEGA/CHARA interferometric data secured in the visible domain, and compare it to the relation based on infrared interferometric observations, which were used to derive the distance to the LMC. Observations of eight G-K giants were obtained with the VEGA/CHARA instrument. The derived limb-darkened angular diameters were combined with a homogeneous set of infrared magnitudes in order to constrain the SBCR. The average precision we obtain on the limb-darkened angular diameters of the eight stars in our sample is 2.4%. For the four stars in common observed by both VEGA/CHARA and PIONIER/VLTI, we find a 1 sigma agreement for the angular diameters. The SBCR we obtain in the visible has a dispersion of 0.04 magnitude and is consistent with the one derived in the infrared (0.018 magnitude). The consistency of the infrared and visible angular diameters and SBCR reinforces the result of 1% precision and accuracy recently achieved on the distance of the LMC using the eclipsing-binary technique. It also indicates that it is possible to combine interferometric observations at different wavelengths when the SBCR is calibrated.


title: Predicted masses of Galactic Cepheids in the Gaia Data Release 2

abstract: On the basis of recently computed nonlinear convective pulsation models of Galactic Cepheids, spanning wide ranges of input stellar parameters, we derive theoretical mass-dependent Period-Wesenheit relations in the Gaia bands, namely $G$, $G_{BP}$ and $G_{RP}$, that are found to be almost independent of the assumed efficiency of super-adiabatic convection. The application to a selected sub-sample of Gaia Data Release 2 Galactic Cepheids database allows us to derive mass-dependent estimates of their individual distances. By imposing their match with the astrometric values inferred from Gaia, we are able to evaluate the individual mass of each pulsator. The inferred mass distribution is peaked around 5.6$M_{\odot}$ and 5.4$M_{\odot}$ for the F and FO pulsators, respectively. If the estimated Gaia parallax offset $<\Delta\varpi>$=0.046 mas is applied to Gaia parallaxes before imposing their coincidence with the theoretical ones, the inferred mass distribution is found to shift towards lower masses, namely $\sim$5.2$M_{\odot}$ and 5.1$M_{\odot}$ for the F and FO pulsators, respectively. The comparison with independent evaluations of the stellar masses, for a subset of binary Cepheids in our sample, seems to support the predictive capability of current theoretical scenario. By forcing the coincidence of our mass determinations with these literature values we derive an independent estimate of the mean offset to be applied to Gaia DR2 parallaxes, $<\Delta\varpi>$=0.053 $\pm$ 0.029 mas, slightly higher but in agreement within the errors with Riess2018 value.


title: Social Distancing 2.0 with Privacy-Preserving Contact Tracing to Avoid a Second Wave of COVID-19

abstract: How to avoid a second wave of COVID-19 after reopening the economy is a pressing question. The extremely high basic reproductive number $R_0$ (5.7 to 6.4, shown in new studies) of SARS-CoV-2 further complicates the challenge. Here we assess effects of Social distancing 2.0, i.e. proximity alert (to maintain inter-personal distance) plus privacy-preserving contact tracing. To solve the dual task, we developed an open source mobile app. The app uses a Bluetooth-based, decentralized contact tracing platform over which the anonymous user ID cannot be linked by the government or a third party. Modelling results show that a 50\% adoption rate of Social distancing 2.0, with privacy-preserving contact tracing, would suffice to decrease the $R_0$ to less than 1 and prevent the resurgence of COVID-19 epidemic.


title: On Dynamic Substructuring of Systems with Localised Nonlinearities

abstract: Dynamic substructuring (DS) methods encompass a range of techniques to decompose large structural systems into multiple coupled subsystems. This decomposition has the principle benefit of reducing computational time for dynamic simulation of the system. In this context, DS methods may form an essential component of hybrid simulation, wherein they can be used to couple physical and numerical substructures at reduced computational cost. Since most engineered systems are inherently nonlinear, particular potential lies in incorporating nonlinear methods in existing substructuring schemes which are largely linear methods. The most widely used and studied DS methods are classical linear techniques such as the Craig-Bampton (CB) method. However, as linear methods they naturally break down in the presence of nonlinearities. Recent advancements in substructuring have involved the development of enrichments to linear methods, which allow for some degree of nonlinearity to be captured. The use of mode shape derivatives has been shown to be able to capture geometrically non-linear effects as an extension to the CBmethod. Other candidates include the method of Finite Element Tearing and Interconnecting. In this work, a virtual hybrid simulation is presented in which a linear elastic vehicle frame supported on four nonlinear spring damper isolators is decomposed into separate domains. One domain consisting of the finite element model of the vehicle frame, which is reduced using the CB method. The second domain consists of the nonlinear isolators whose restoring forces are characterised by nonlinear spring and damper forces. Coupling between the models is carried out using a Lagrange multiplier method and time series simulations of the system are conducted and compared to the full global system with regards to simulation time and accuracy.