title: Remote Sensing Image Scene Classification with Self-Supervised Paradigm under Limited Labeled Samples

abstract: With the development of deep learning, supervised learning methods perform well in remote sensing images (RSIs) scene classification. However, supervised learning requires a huge number of annotated data for training. When labeled samples are not sufficient, the most common solution is to fine-tune the pre-training models using a large natural image dataset (e.g. ImageNet). However, this learning paradigm is not a panacea, especially when the target remote sensing images (e.g. multispectral and hyperspectral data) have different imaging mechanisms from RGB natural images. To solve this problem, we introduce new self-supervised learning (SSL) mechanism to obtain the high-performance pre-training model for RSIs scene classification from large unlabeled data. Experiments on three commonly used RSIs scene classification datasets demonstrated that this new learning paradigm outperforms the traditional dominant ImageNet pre-trained model. Moreover, we analyze the impacts of several factors in SSL on RSIs scene classification tasks, including the choice of self-supervised signals, the domain difference between the source and target dataset, and the amount of pre-training data. The insights distilled from our studies can help to foster the development of SSL in the remote sensing community. Since SSL could learn from unlabeled massive RSIs which are extremely easy to obtain, it will be a potentially promising way to alleviate dependence on labeled samples and thus efficiently solve many problems, such as global mapping.


title: DIETERpy: a Python framework for The Dispatch and Investment Evaluation Tool with Endogenous Renewables

abstract: DIETER is an open-source power sector model designed to analyze future settings with very high shares of variable renewable energy sources. It minimizes overall system costs, including fixed and variable costs of various generation, flexibility and sector coupling options. Here we introduce DIETERpy that builds on the existing model version, written in the General Algebraic Modeling System (GAMS), and enhances it with a Python framework. This combines the flexibility of Python regarding pre- and post-processing of data with a straightforward algebraic formulation in GAMS and the use of efficient solvers. DIETERpy also offers a browser-based graphical user interface. The new framework is designed to be easily accessible as it enables users to run the model, alter its configuration, and define numerous scenarios without a deeper knowledge of GAMS. Code, data, and manuals are available in public repositories under permissive licenses for transparency and reproducibility.


title: Impact of Short Blocklength Coding on Stability of an AGV Control System in Industry 4.0

abstract: With the advent of 5G and beyond, using wireless communication for closed-loop control and automation processes is one of the main aspects of the envisioned Industry 4.0. In this regard, a major challenge is to ensure a robust and stable control system over an unreliable wireless channel. One of the main use-cases in this context is Automated Guided Vehicle (AGV) control in a future factory. Specifically, we consider a system where an AGV controller is placed in an edge cloud in the factory network infrastructure and the control commands are sent over a time-correlated Rayleigh fading channel. In an industrial control, short packets are exchanged between the controller and the actuator. Therefore, in this case, Shannon's assumption for an infinite block length is not applicable. The objective is to analyse the stability performance of an AGV control system in a Finite Block-Length (FBL) regime. We evaluate the coding rate required to maintain a stable edge cloud based AGV system. The results illustrate that adapting the control parameters can lower the stringent requirements on the coding rate. It reveals that a constant stability performance can be achieved even at higher coding rate by increasing the AGV's velocity. Moreover, this paper also determines the maximum number of AGVs that can be served seamlessly over the available communication resources while maintaining a stable control system.


title: Optimization Landscapes of Wide Deep Neural Networks Are Benign

abstract: We analyze the optimization landscapes of deep learning with wide networks. We highlight the importance of constraints for such networks and show that constraint -- as well as unconstraint -- empirical-risk minimization over such networks has no confined points, that is, suboptimal parameters that are difficult to escape from. Hence, our theories substantiate the common belief that wide neural networks are not only highly expressive but also comparably easy to optimize.


title: Exponential Independence in Subcubic Graphs

abstract: A set $S$ of vertices of a graph $G$ is exponentially independent if, for every vertex $u$ in $S$, $$\sum\limits_{v\in S\setminus \{ u\}}\left(\frac{1}{2}\right)^{{\rm dist}_{(G,S)}(u,v)-1}<1,$$ where ${\rm dist}_{(G,S)}(u,v)$ is the distance between $u$ and $v$ in the graph $G-(S\setminus \{ u,v\})$. The exponential independence number $\alpha_e(G)$ of $G$ is the maximum order of an exponentially independent set in $G$. In the present paper we present several bounds on this parameter and highlight some of the many related open problems. In particular, we prove that subcubic graphs of order $n$ have exponentially independent sets of order $\Omega(n/\log^2(n))$, that the infinite cubic tree has no exponentially independent set of positive density, and that subcubic trees of order $n$ have exponentially independent sets of order $(n+3)/4$.


title: Investigating the magnetospheric accretion process in the young pre-transitional disk system DoAr 44 (V2062~Oph). A multiwavelength interferometric, spectropolarimetric, and photometric observing campaign

abstract: Young stars interact with their accretion disk through their strong magnetosphere. We investigate the magnetospheric accretion process in the young stellar system DoAr 44. We monitored the system over several rotational cycles, combining high-resolution optical and near-IR spectropolarimetry with long-baseline near-IR interferometry and multicolor photometry. DoAr 44 is a young 1.2 solar mass star, moderately accreting from its disk, and seen at a low inclination. We derive a rotational period of 2.96 d from the system's light curve. Several optical and near-IR line profiles probing the accretion funnel flows and the accretion shock are modulated at the stellar rotation period. The most variable line profile, HeI 1083 nm, exhibits modulated redshifted wings a signature of accretion funnel flows, as well as deep blueshifted absorptions indicative of transient outflows. The Zeeman-Doppler analysis suggests the star hosts a mainly dipolar magnetic field, inclined by about 20 deg. onto the spin axis, with an intensity reaching about 800 G at the photosphere, and up to 2 +/- 0.8 kG close to the accretion shock. The magnetic field appears strong enough to disrupt the inner disk close to the corotation radius, at a distance of about 4.6 stellar radii (0.043 au). This supports the upper limit of 5 stellar radii (0.047 au) we derived for the size of the magnetosphere from long baseline interferometry. DoAr 44 is a pre-transitional disk system, exhibiting a 25-30 au gap in its circumstellar disk, with the inner and outer disks being misaligned. On a scale of 0.1 au or less, our results indicate that the system steadily accretes from its inner disk through its tilted dipolar magnetosphere. We conclude that in spite of a highly structured outer disk, perhaps the signature of ongoing planetary formation, the magnetospheric accretion process proceeds unimpeded at the star-disk interaction level.


title: Dynamical generalization of Yetter's model based on a crossed module of discrete groups

abstract: We construct a dynamical lattice model based on a crossed module of possibly non-abelian finite groups. Its degrees of freedom are defined on links and plaquettes, while gauge transformations are based on vertices and links of the underlying lattice. We specify the Hilbert space, define basic observables (including the Hamiltonian) and initiate a~discussion on the model's phase diagram. The constructed model generalizes, and in appropriate limits reduces to, topological theories with symmetries described by groups and crossed modules, lattice Yang-Mills theory and $2$-form electrodynamics. We conclude by reviewing classifying spaces of crossed modules, with an emphasis on the direct relation between their geometry and properties of gauge theories under consideration.

categories: math-ph cond-mat.str-el hep-lat hep-th math.MP


title: Time Matters: Time-Aware LSTMs for Predictive Business Process Monitoring

abstract: Predictive business process monitoring (PBPM) aims to predict future process behavior during ongoing process executions based on event log data. Especially, techniques for the next activity and timestamp prediction can help to improve the performance of operational business processes. Recently, many PBPM solutions based on deep learning were proposed by researchers. Due to the sequential nature of event log data, a common choice is to apply recurrent neural networks with long short-term memory (LSTM) cells. We argue, that the elapsed time between events is informative. However, current PBPM techniques mainly use 'vanilla' LSTM cells and hand-crafted time-related control flow features. To better model the time dependencies between events, we propose a new PBPM technique based on time-aware LSTM (T-LSTM) cells. T-LSTM cells incorporate the elapsed time between consecutive events inherently to adjust the cell memory. Furthermore, we introduce cost-sensitive learning to account for the common class imbalance in event logs. Our experiments on publicly available benchmark event logs indicate the effectiveness of the introduced techniques.


title: OpenTraj: Assessing Prediction Complexity in Human Trajectories Datasets

abstract: Human Trajectory Prediction (HTP) has gained much momentum in the last years and many solutions have been proposed to solve it. Proper benchmarking being a key issue for comparing methods, this paper addresses the question of evaluating how complex is a given dataset with respect to the prediction problem. For assessing a dataset complexity, we define a series of indicators around three concepts: Trajectory predictability; Trajectory regularity; Context complexity. We compare the most common datasets used in HTP in the light of these indicators and discuss what this may imply on benchmarking of HTP algorithms. Our source code is released on Github.


title: Alignment and scattering of colliding active droplets

abstract: Active droplets emit a chemical solute at their surface that modifies their local interfacial tension. They exploit the nonlinear coupling of the convective transport of solute to the resulting Marangoni flows to self-propel. Such swimming droplets are by nature anti-chemotactic and are repelled by their own chemical wake or their neighbours. The rebound dynamics resulting from pairwise droplet interactions was recently analysed in detail for purely head-on collisions using a specific bispherical approach. Here, we extend this analysis and propose a reduced model of a generic collision to characterise the alignment and scattering properties of oblique droplet collisions and their potential impact on collective droplet dynamics. A systematic alignment of the droplets' trajectories is observed for symmetric collisions, when the droplets interact directly, and arises from the finite-time rearrangement of the droplets' chemical wake during the collision. For more generic collisions, complex and diverse dynamical regimes are observed, whether the droplets interact directly or through their chemical wake, resulting in a significant scattering.


title: Variance-Reduced Methods for Machine Learning

abstract: Stochastic optimization lies at the heart of machine learning, and its cornerstone is stochastic gradient descent (SGD), a method introduced over 60 years ago. The last 8 years have seen an exciting new development: variance reduction (VR) for stochastic optimization methods. These VR methods excel in settings where more than one pass through the training data is allowed, achieving a faster convergence than SGD in theory as well as practice. These speedups underline the surge of interest in VR methods and the fast-growing body of work on this topic. This review covers the key principles and main developments behind VR methods for optimization with finite data sets and is aimed at non-expert readers. We focus mainly on the convex setting, and leave pointers to readers interested in extensions for minimizing non-convex functions.


title: Weight Encode Reconstruction Network for Computed Tomography in a Semi-Case-Wise and Learning-Based Way

abstract: Classic algebraic reconstruction technology (ART) for computed tomography requires pre-determined weights of the voxels for projecting pixel values. However, such weight cannot be accurately obtained due to the limitation of the physical understanding and computation resources. In this study, we propose a semi-case-wise learning-based method named Weight Encode Reconstruction Network (WERNet) to tackle the issues mentioned above. The model is trained in a self-supervised manner without the label of a voxel set. It contains two branches, including the voxel weight encoder and the voxel attention part. Using gradient normalization, we are able to co-train the encoder and voxel set numerically stably. With WERNet, the reconstructed result was obtained with a cosine similarity greater than 0.999 with the ground truth. Moreover, the model shows the extraordinary capability of denoising comparing to the classic ART method. In the generalization test of the model, the encoder is transferable from a voxel set with complex structure to the unseen cases without the deduction of the accuracy.


title: Theta characteristics on tropical curves

abstract: We give an explicit description of theta-characteristics on tropical curves and characterize the effective ones. We construct the moduli space for tropical theta-characteristics of genus g as a generalized cone complex. We describe the fibers of the specialization map from the classical moduli space of theta-characteristics on algebraic curves to our moduli space of tropical theta-characteristics.


title: Discontinuous Ground States for NLSE on $\mathbb{R}$ with a F\"{u}l\"{o}p-Tsutsui $\delta$ interaction

abstract: We analyse the existence and the stability of the ground states of the one-dimensional nonlinear Schr\"{o}dinger equation with a focusing power nonlinearity and a defect located at the origin. In this paper a ground state is defined as a global minimizer of the action functional on the Nehari manifold and the defect considered is a F\"{u}l\"{o}p-Tsutsui $\delta$ type, namely a $\delta$ condition that allows discontinuities. The existence of ground states is proved by variational techniques, while the stability results from the Grillakis-Shatah-Strauss theory.


title: Improving performances of MCMC for Nearest Neighbor Gaussian Process models with full data augmentation

abstract: Even though Nearest Neighbor Gaussian Processes (NNGP) alleviate considerably MCMC implementation of Bayesian space-time models, they do not solve the convergence problems caused by high model dimension. Frugal alternatives such as response or collapsed algorithms are an answer.gree Our approach is to keep full data augmentation but to try and make it more efficient. We present two strategies to do so. The first scheme is to pay a particular attention to the seemingly trivial fixed effects of the model. We show empirically that re-centering the latent field on the intercept critically improves chain behavior. We extend this approach to other fixed effects that may interfere with a coherent spatial field. We propose a simple method that requires no tuning while remaining affordable thanks to NNGP's sparsity. The second scheme accelerates the sampling of the random field using Chromatic samplers. This method makes long sequential simulation boil down to group-parallelized or group-vectorized sampling. The attractive possibility to parallelize NNGP likelihood can therefore be carried over to field sampling. We present a R implementation of our methods for Gaussian fields in the public repository https://github.com/SebastienCoube/Improving_NNGP_full_augmentation . An extensive vignette is provided. We run our implementation on two synthetic toy examples along with the state of the art package spNNGP. Finally, we apply our method on a real data set of lead contamination in the United States of America mainland.


title: Shape and dynamics of nonrelativistic vortex strings in parity-breaking media

abstract: The shape and dynamics of the nonrelativistic gauge vortex string in the parity-broken media is considered, upon reducing the problem to finding the extremum of the Abelian Higgs model effective action with the fixed B-type helicity of the gauge field. It is shown that in contrast with the case of the fixed A-type helicity, the static solution of the Ginzburg-Landau energy functional in the London limit is the helix with the specific relation between the curvature and torsion of the vortex line depending on the strength of the space parity violating contribution of the Lifshitz invariant. A nonlinear dynamical equation is linearized in case of small oscillations around the helical contour, and the polarization and dispersion law of the propagated waves are obtained.


title: Edge mode engineering for optimal ultracoherent silicon nitride membranes

abstract: Due to their high force sensitivity, mechanical resonators combining low mechanical dissipation with a small motional mass are highly demanded in fields as diverse as resonant force imaging, mass sensing, or cavity optomechanics. "Soft-clamping" is a phononic engineering technique by which mechanical modes of highly-stressed membranes or strings are localized away from lossy regions, thereby enabling ultrahigh-Q for ng-scale devices. Here, we report on parasitic modes arising from the finite size of the structure which can significantly degrade the performance of these membranes. Through interferometric measurements and finite-element simulations, we show that these parasitic modes can hybridize with the localized modes of our structures, reducing the quality factors by up to one order of magnitude. To circumvent this problem, we engineer the spectral profile of the parasitic modes in order to avoid their overlap with the high-Q defect mode. We verify via a statistical analysis that this modal engineering reproducibly yields higher quality factors in fabricated devices, consistent with theoretically predicted values.


title: On the lower bounds for real double Hurwitz numbers

abstract: As the real counterpart of double Hurwitz number, the real double Hurwitz number depends on the distribution of real branch points. We consider the problem of asymptotic growth of real and complex double Hurwitz numbers. We provide a lower bound for real double Hurwitz numbers based on the tropical computation of real double Hurwitz numbers. By using this lower bound and J. Rau's result ( Math. Ann. 375(1-2): 895-915, 2019), we prove the logarithmic equivalence of real and complex Hurwitz numbers.


title: The Origin of Hadron Masses

abstract: The color confinement can be decently explained by assuming the global $SU(3)$ color symmetry. A hadron is viewed as a bag of a finite size, whose energy is contributed by the color fields within the hadron. In the large momentum frame, the mass of a hadron can be quantized and can be expressed as the sum of the current masses of quarks involved plus a flavor independent term which depends on a universal parameter and can be identified to be the contribution from the QCD trace anomaly. This term depicts the Regge's trajectories of light hadron masses, from which the parameter can be determined and gives the minimal vector meson mass and baryon mass to be 761 MeV and 935 MeV, respectively. The mass formula can explain many features of the hadron spectroscopy.


title: Two-variable logic has weak, but not strong, Beth definability

abstract: We prove that the two-variable fragment of first-order logic has the weak Beth definability property. This makes the two-variable fragment a natural logic separating the weak and the strong Beth properties since it does not have the strong Beth definability property.


title: Simflowny 3: An upgraded platform for scientific modelling and simulation

abstract: Simflowny is an open platform which automatically generates efficient parallel code of scientific dynamical models for different simulation frameworks. Here we present major upgrades on this software to support simultaneously a quite generic family of partial differential equations. These equations can be discretized using: (i) standard finite-difference for systems with derivatives up to any order, (ii) High-Resolution-Shock-Capturing methods to deal with shocks and discontinuities of balance law equations, and (iii) particle-based methods. We have improved the adaptive-mesh-refinement algorithms to preserve the convergence order of the numerical methods, which is a requirement for improving scalability. Finally, we have also extended our graphical user interface (GUI) to accommodate these and future families of equations. This paper summarizes the formal representation and implementation of these new families, providing several validation results.


title: An Evaluation of Classification Methods for 3D Printing Time-Series Data

abstract: Additive Manufacturing presents a great application area for Machine Learning because of the vast volume of data generated and the potential to mine this data to control outcomes. In this paper we present preliminary work on classifying infrared time-series data representing melt-pool temperature in a metal 3D printing process. Our ultimate objective is to use this data to predict process outcomes (e.g. hardness, porosity, surface roughness). In the work presented here we simply show that there is a signal in this data that can be used for the classification of different components and stages of the AM process. In line with other Machine Learning research on time-series classification we use k-Nearest Neighbour classifiers. The results we present suggests that Dynamic Time Warping is an effective distance measure compared with alternatives for 3D printing data of this type.


title: Autoregressive Entity Retrieval

abstract: Entities are at the center of how we represent and aggregate knowledge. For instance, Encyclopedias such as Wikipedia are structured by entities (e.g., one per Wikipedia article). The ability to retrieve such entities given a query is fundamental for knowledge-intensive tasks such as entity linking and open-domain question answering. Current approaches can be understood as classifiers among atomic labels, one for each entity. Their weight vectors are dense entity representations produced by encoding entity meta information such as their descriptions. This approach has several shortcomings: (i) context and entity affinity is mainly captured through a vector dot product, potentially missing fine-grained interactions; (ii) a large memory footprint is needed to store dense representations when considering large entity sets; (iii) an appropriately hard set of negative data has to be subsampled at training time. In this work, we propose GENRE, the first system that retrieves entities by generating their unique names, left to right, token-by-token in an autoregressive fashion. This mitigates the aforementioned technical issues since: (i) the autoregressive formulation directly captures relations between context and entity name, effectively cross encoding both; (ii) the memory footprint is greatly reduced because the parameters of our encoder-decoder architecture scale with vocabulary size, not entity count; (iii) the softmax loss is computed without subsampling negative data. We experiment with more than 20 datasets on entity disambiguation, end-to-end entity linking and document retrieval tasks, achieving new state-of-the-art or very competitive results while using a tiny fraction of the memory footprint of competing systems. Finally, we demonstrate that new entities can be added by simply specifying their names. Code and pre-trained models at https://github.com/facebookresearch/GENRE.


title: Estimating Atmospheric Parameters of DA White Dwarf Stars with Deep Learning

abstract: The determination of atmospheric parameters of white dwarf stars (WDs) is crucial for researches on them. Traditional methodology is to fit the model spectra to observed absorption lines and report the parameters with the lowest $\chi ^2$ error, which strongly relies on theoretical models that are not always publicly accessible. In this work, we construct a deep learning network to model-independently estimate Teff and log g of DA stars (DAs), corresponding to WDs with hydrogen dominated atmospheres. The network is directly trained and tested on the normalized flux pixels of full optical wavelength range of DAs spectroscopically confirmed in the Sloan Digital Sky Survey (SDSS). Experiments in test part yield that the root mean square error (RMSE) for Teff and log g approaches to 900 K and 0.1 dex, respectively. This technique is applicable for those DAs with Teff from 5000 K to 40000 K and log g from 7.0 dex to 9.0 dex. Furthermore, the applicability of this method is verified for the spectra with degraded resolution $\sim 200$. So it is also practical for the analysis of DAs that will be detected by the Chinese Space Station Telescope (CSST).


title: Quantifying Privacy Leakage in Graph Embedding

abstract: Graph embeddings have been proposed to map graph data to low dimensional space for downstream processing (e.g., node classification or link prediction). With the increasing collection of personal data, graph embeddings can be trained on private and sensitive data. For the first time, we quantify the privacy leakage in graph embeddings through three inference attacks targeting Graph Neural Networks. We propose a membership inference attack to infer whether a graph node corresponding to individual user's data was member of the model's training or not. We consider a blackbox setting where the adversary exploits the output prediction scores, and a whitebox setting where the adversary has also access to the released node embeddings. This attack provides an accuracy up to 28% (blackbox) 36% (whitebox) beyond random guess by exploiting the distinguishable footprint between train and test data records left by the graph embedding. We propose a Graph Reconstruction attack where the adversary aims to reconstruct the target graph given the corresponding graph embeddings. Here, the adversary can reconstruct the graph with more than 80% of accuracy and link inference between two nodes around 30% more confidence than a random guess. We then propose an attribute inference attack where the adversary aims to infer a sensitive attribute. We show that graph embeddings are strongly correlated to node attributes letting the adversary inferring sensitive information (e.g., gender or location).


title: Tubular Shape Aware Data Generation for Semantic Segmentation in Medical Imaging

abstract: Chest X-ray is one of the most widespread examinations of the human body. In interventional radiology, its use is frequently associated with the need to visualize various tube-like objects, such as puncture needles, guiding sheaths, wires, and catheters. Detection and precise localization of these tube-like objects in the X-ray images is, therefore, of utmost value, catalyzing the development of accurate target-specific segmentation algorithms. Similar to the other medical imaging tasks, the manual pixel-wise annotation of the tubes is a resource-consuming process. In this work, we aim to alleviate the lack of the annotated images by using artificial data. Specifically, we present an approach for synthetic data generation of the tube-shaped objects, with a generative adversarial network being regularized with a prior-shape constraint. Our method eliminates the need for paired image--mask data and requires only a weakly-labeled dataset (10--20 images) to reach the accuracy of the fully-supervised models. We report the applicability of the approach for the task of segmenting tubes and catheters in the X-ray images, whereas the results should also hold for the other imaging modalities.


title: Exploring the Mass Segregation Effect of X-ray Sources in Globular Clusters. IV. Evidence of Black Hole Burning in $\omega$ Centauri

abstract: Using X-ray sources as sensitive probes of stellar dynamical interactions in globular clusters (GCs), we study the mass segregation and binary burning processes in $\omega$ Cen. We show that the mass segregation of X-ray sources is quenched in $\omega$ Cen, while the X-ray source abundance of $\omega$ Cen is much smaller than other GCs, and the binary hardness ratio (defined as $L_{\rm X}/(L_{\rm K}f_{b})$, with $f_{b}$ the binary fraction, $L_{\rm X}$ and $L_{\rm K}$ the cumulative X-ray and K band luminosity of GCs, respectively) of $\omega$ Cen is located far below the $L_{\rm X}/(L_{\rm K}f_{b})-\sigma_{c}$ correlation line of the dynamically normal GCs. These evidences suggest that the binary burning processes are highly suppressed in $\omega$ Cen, and other "heating mechanisms", very likely a black hole subsystem (BHS), are essential in the dynamical evolution of $\omega$ Cen. Through the "black hole burning" processes (i.e., dynamical hardening of the BH binaries), the BHS is dominating the energy production of $\omega$ Cen, which also makes $\omega$ Cen a promising factory of gravitational-wave sources in the Galaxy.


title: Self-Play Reinforcement Learning for Fast Image Retargeting

abstract: In this study, we address image retargeting, which is a task that adjusts input images to arbitrary sizes. In one of the best-performing methods called MULTIOP, multiple retargeting operators were combined and retargeted images at each stage were generated to find the optimal sequence of operators that minimized the distance between original and retargeted images. The limitation of this method is in its tremendous processing time, which severely prohibits its practical use. Therefore, the purpose of this study is to find the optimal combination of operators within a reasonable processing time; we propose a method of predicting the optimal operator for each step using a reinforcement learning agent. The technical contributions of this study are as follows. Firstly, we propose a reward based on self-play, which will be insensitive to the large variance in the content-dependent distance measured in MULTIOP. Secondly, we propose to dynamically change the loss weight for each action to prevent the algorithm from falling into a local optimum and from choosing only the most frequently used operator in its training. Our experiments showed that we achieved multi-operator image retargeting with less processing time by three orders of magnitude and the same quality as the original multi-operator-based method, which was the best-performing algorithm in retargeting tasks.


title: Continual Learning for Natural Language Generation in Task-oriented Dialog Systems

abstract: Natural language generation (NLG) is an essential component of task-oriented dialog systems. Despite the recent success of neural approaches for NLG, they are typically developed in an offline manner for particular domains. To better fit real-life applications where new data come in a stream, we study NLG in a "continual learning" setting to expand its knowledge to new domains or functionalities incrementally. The major challenge towards this goal is catastrophic forgetting, meaning that a continually trained model tends to forget the knowledge it has learned before. To this end, we propose a method called ARPER (Adaptively Regularized Prioritized Exemplar Replay) by replaying prioritized historical exemplars, together with an adaptive regularization technique based on ElasticWeight Consolidation. Extensive experiments to continually learn new domains and intents are conducted on MultiWoZ-2.0 to benchmark ARPER with a wide range of techniques. Empirical results demonstrate that ARPER significantly outperforms other methods by effectively mitigating the detrimental catastrophic forgetting issue.


title: Proving Highly-Concurrent Traversals Correct

abstract: Modern highly-concurrent search data structures, such as search trees, obtain multi-core scalability and performance by having operations traverse the data structure without any synchronization. As a result, however, these algorithms are notoriously difficult to prove linearizable, which requires identifying a point in time in which the traversal's result is correct. The problem is that traversing the data structure as it undergoes modifications leads to complex behaviors, necessitating intricate reasoning about all interleavings of reads by traversals and writes mutating the data structure. In this paper, we present a general proof technique for proving unsynchronized traversals correct in a significantly simpler manner, compared to typical concurrent reasoning and prior proof techniques. Our framework relies only on sequential properties} of traversals and on a conceptually simple and widely-applicable condition about the ways an algorithm's writes mutate the data structure. Establishing that a target data structure satisfies our condition requires only simple concurrent reasoning, without considering interactions of writes and reads. This reasoning can be further simplified by using our framework. To demonstrate our technique, we apply it to prove several interesting and challenging concurrent binary search trees: the logical-ordering AVL tree, the Citrus tree, and the full contention-friendly tree. Both the logical-ordering tree and the full contention-friendly tree are beyond the reach of previous approaches targeted at simplifying linearizability proofs.


title: GECKO: Reconciling Privacy, Accuracy and Efficiency in Embedded Deep Learning

abstract: Embedded systems demand on-device processing of data using Neural Networks (NNs) while conforming to the memory, power and computation constraints, leading to an efficiency and accuracy tradeoff. To bring NNs to edge devices, several optimizations such as model compression through pruning, quantization, and off-the-shelf architectures with efficient design have been extensively adopted. These algorithms when deployed to real world sensitive applications, requires to resist inference attacks to protect privacy of users training data. However, resistance against inference attacks is not accounted for designing NN models for IoT. In this work, we analyse the three-dimensional privacy-accuracy-efficiency tradeoff in NNs for IoT devices and propose Gecko training methodology where we explicitly add resistance to private inferences as a design objective. We optimize the inference-time memory, computation, and power constraints of embedded devices as a criterion for designing NN architecture while also preserving privacy. We choose quantization as design choice for highly efficient and private models. This choice is driven by the observation that compressed models leak more information compared to baseline models while off-the-shelf efficient architectures indicate poor efficiency and privacy tradeoff. We show that models trained using Gecko methodology are comparable to prior defences against black-box membership attacks in terms of accuracy and privacy while providing efficiency.


title: Static Output Feedback in an Anisotropic Norm Setup Revisited

abstract: The synthesis problem of static output feedback controllers within the anistropic-norm setup is revisited. A tractable synthesis approach involving iterations over a convex optimization problem is suggested, similarly to existing results for the Hinfinity-norm minimization case. The results are formulated by a couple of Linear Matrix Inequalities coupled via a bilinear equality, revealing, as in the Hinfinity case, the duality between the control-type and filtering type LMIs and allowing a tractable iterative method to cope with practical static output feedback synthesis problems. The resulting optimization scheme is then applied to a flight control problem, where the merit of the anisotropic norm setup is shown to provide a useful trade-off between closed loop response and feedback gains.


title: Coded Stochastic ADMM for Decentralized Consensus Optimization with Edge Computing

abstract: Big data, including applications with high security requirements, are often collected and stored on multiple heterogeneous devices, such as mobile devices, drones and vehicles. Due to the limitations of communication costs and security requirements, it is of paramount importance to extract information in a decentralized manner instead of aggregating data to a fusion center. To train large-scale machine learning models, edge/fog computing is often leveraged as an alternative to centralized learning. We consider the problem of learning model parameters in a multi-agent system with data locally processed via distributed edge nodes. A class of mini-batch stochastic alternating direction method of multipliers (ADMM) algorithms is explored to develop the distributed learning model. To address two main critical challenges in distributed networks, i.e., communication bottleneck and straggler nodes (nodes with slow responses), error-control-coding based stochastic incremental ADMM is investigated. Given an appropriate mini-batch size, we show that the mini-batch stochastic ADMM based method converges in a rate of $O(\frac{1}{\sqrt{k}})$, where $k$ denotes the number of iterations. Through numerical experiments, it is revealed that the proposed algorithm is communication-efficient, rapidly responding and robust in the presence of straggler nodes compared with state of the art algorithms.


title: Sharp lower error bounds for strong approximation of SDEs with discontinuous drift coefficient by coupling of noise

abstract: In the past decade, an intensive study of strong approximation of stochastic differential equations (SDEs) with a drift coefficient that has discontinuities in space has begun. In the majority of these results it is assumed that the drift coefficient satisfies piecewise regularity conditions and that the diffusion coefficient is globally Lipschitz continuous and non-degenerate at the discontinuities of the drift coefficient. Under this type of assumptions the best $L_p$-error rate obtained so far for approximation of scalar SDEs at the final time is $3/4$ in terms of the number of evaluations of the driving Brownian motion. In the present article we prove for the first time in the literature sharp lower error bounds for such SDEs. We show that for a huge class of additive noise driven SDEs of this type the $L_p$-error rate $3/4$ can not be improved. For the proof of this result we employ a novel technique by studying equations with coupled noise: we reduce the analysis of the $L_p$-error of an arbitrary approximation based on evaluation of the driving Brownian motion at finitely many times to the analysis of the $L_p$-distance of two solutions of the same equation that are driven by Brownian motions that are coupled at the given time-points and independent, conditioned on their values at these points. To obtain lower bounds for the latter quantity, we prove a new quantitative version of positive association for bivariate normal random variables $(Y,Z)$ by providing explict lower bounds for the covariance $\text{Cov}(f(Y),g(Z))$ in case of piecewise Lipschitz continuous functions $f$ and $g$. In addition it turns out that our proof technique also leads to lower error bounds for estimating occupation time functionals $\int_0^1 f(W_t)\, dt$ of a Brownian motion $W$, which substantially extends known results for the case of $f$ being an indicator function.


title: Gravitational Wave observatories may be able to detect hyperbolic encounters of Black Holes

abstract: Gravitational Wave (GW) astronomy promises to observe different kinds of astrophysical sources. Here we explore the possibility of detection of GWs from hyperbolic interactions of compact stars with ground-based interferometric detectors. It is believed that a bound compact cluster, such as a globular cluster, can be a primary environment for these interactions. We estimate the detection rates for such events by considering local geometry within the cluster, accounting for scattering probability of compact stars at finite distances, and assuming realistic cluster properties guided by available numerical models, their formation times, and evolution of stars inside them. We find that, even in the conservative limit, it may be possible to detect such black hole encounters in the next few years by the present network of observatories with the ongoing sensitivity upgrades and one to few events per year with the next generation observatories. In practice, actual detection rates can significantly surpass the estimated average rates, since the chances of finding outliers in a very large population can be high. Such observations (or, no observation) may provide crucial constraints to estimate the number of isolated compact stars in the universe. These detections will be exciting discoveries on their own and will be complimentary to observations of binary mergers bringing us one step closer to address a fundamental question, how many black holes are there in the observable universe.


title: The Sparse Vector Technique, Revisited

abstract: We revisit one of the most basic and widely applicable techniques in the literature of differential privacy - the sparse vector technique [Dwork et al., STOC 2009]. This simple algorithm privately tests whether the value of a given query on a database is close to what we expect it to be. It allows to ask an unbounded number of queries as long as the answer is close to what we expect, and halts following the first query for which this is not the case. We suggest an alternative, equally simple, algorithm that can continue testing queries as long as any single individual does not contribute to the answer of too many queries whose answer deviates substantially form what we expect. Our analysis is subtle and some of its ingredients may be more widely applicable. In some cases our new algorithm allows to privately extract much more information from the database than the original. We demonstrate this by applying our algorithm to the shifting heavy-hitters problem: On every time step, each of $n$ users gets a new input, and the task is to privately identify all the current heavy-hitters. That is, on time step $i$, the goal is to identify all data elements $x$ such that many of the users have $x$ as their current input. We present an algorithm for this problem with improved error guarantees over what can be obtained using existing techniques. Specifically, the error of our algorithm depends on the maximal number of times that a single user holds a heavy-hitter as input, rather than the total number of times in which a heavy-hitter exists.


title: Are Artificial Dendrites useful in NeuroEvolution?

abstract: The significant role of dendritic processing within neuronal networks has become increasingly clear. This letter explores the effects of including a simple dendrite-inspired mechanism into neuroevolution. The phenomenon of separate dendrite activation thresholds on connections is allowed to emerge under an evolutionary process. It is shown how such processing can be positively selected for, particularly for connections between the hidden and output layer, and increases performance.


title: A General Equivalent Circuit for Lossy Non-Symmetric Reciprocal Two-Ports

abstract: A general equivalent circuit for two-port lossy non-symmetric reciprocal networks is proposed. The equivalent circuit is based on the eigenstate decomposition. It has three complex parameters that are obtained from the eigenvalues and eigenvectors of the admittance or impedance matrix of the network in a straightforward way. The interconnection of the circuit elements is simple and compact. The applicability of the equivalent circuit is general, and the real parts of its immittances are always positive. To verify its behavior, the equivalent circuit of three different structures is extracted: an artificial composite right/left-handed transmission line unit cell, a series-fed coupled patch radiating element, and the complementary strip-slot. These cases offer a wide array of results, obtained from theoretical, simulated, and measured data. Nevertheless, the parameters of the equivalent circuit behave as expected in every case study in a very broad bandwidth. The equivalent circuit can be used advantageously in the analysis and design of non-symmetric two-ports, like unit cells of leaky-wave antennas.


title: Hidden automatic sequences

abstract: An automatic sequence is a letter-to-letter coding of a fixed point of a uniform morphism. More generally, we have morphic sequences, which are letter-to-letter codings of fixed points of arbitrary morphisms. There are many examples where an, a priori, morphic sequence with a \emph{non-uniform} morphism happens to be an automatic sequence. An example is the Lys\"enok morphism $a \to aca$, $b \to d$, $c \to b$, $d \to c$, the fixed point of which is also a 2-automatic sequence. Such an identification is useful for the description of the dynamical systems generated by the fixed point. We give several ways to uncover such hidden automatic sequences, and present many examples. We focus in particular on morphisms associated with Grigorchuk(-like) groups.


title: A straightforward line search approach on the expected empirical loss for stochastic deep learning problems

abstract: A fundamental challenge in deep learning is that the optimal step sizes for update steps of stochastic gradient descent are unknown. In traditional optimization, line searches are used to determine good step sizes, however, in deep learning, it is too costly to search for good step sizes on the expected empirical loss due to noisy losses. This empirical work shows that it is possible to approximate the expected empirical loss on vertical cross sections for common deep learning tasks considerably cheaply. This is achieved by applying traditional one-dimensional function fitting to measured noisy losses of such cross sections. The step to a minimum of the resulting approximation is then used as step size for the optimization. This approach leads to a robust and straightforward optimization method which performs well across datasets and architectures without the need of hyperparameter tuning.


title: Spectral energy distributions of dust and PAHs based on the evolution of grain size distribution in galaxies

abstract: Based on a one-zone evolution model of grain size distribution in a galaxy, we calculate the evolution of infrared spectral energy distribution (SED), considering silicate, carbonaceous dust, and polycyclic aromatic hydrocarbons (PAHs). The dense gas fraction ($\eta_\mathrm{dense}$) of the interstellar medium (ISM), the star formation time-scale ($\tau_\mathrm{SF}$), and the interstellar radiation field intensity normalized to the Milky Way value ($U$) are the main parameters. We find that the SED shape generally has weak mid-infrared (MIR) emission in the early phase of galaxy evolution because the dust abundance is dominated by large grains. At an intermediate stage ($t\sim 1$ Gyr for $\tau_\mathrm{SF}=5$ Gyr), the MIR emission grows rapidly because the abundance of small grains increases drastically by the accretion of gas-phase metals. We also compare our results with observational data of nearby and high-redshift ($z\sim 2$) galaxies taken by \textit{Spitzer}. We broadly reproduce the flux ratios in various bands as a function of metallicity. We find that small $\eta_\mathrm{dense}$ (i.e.\ the ISM dominated by the diffuse phase) is favoured to reproduce the 8 $\mu$m intensity dominated by PAHs both for the nearby and the $z\sim 2$ samples. A long $\tau_\mathrm{SF}$ raises the 8 $\mu$m emission to a level consistent with the nearby low-metallicity galaxies. The broad match between the theoretical calculations and the observations supports our understanding of the grain size distribution, but the importance of the diffuse ISM for the PAH emission implies the necessity of spatially resolved treatment for the ISM.


title: Multiple Infrared Small Targets Detection based on Hierarchical Maximal Entropy Random Walk

abstract: The technique of detecting multiple dim and small targets with low signal-to-clutter ratios (SCR) is very important for infrared search and tracking systems. In this paper, we establish a detection method derived from maximal entropy random walk (MERW) to robustly detect multiple small targets. Initially, we introduce the primal MERW and analyze the feasibility of applying it to small target detection. However, the original weight matrix of the MERW is sensitive to interferences. Therefore, a specific weight matrix is designed for the MERW in principle of enhancing characteristics of small targets and suppressing strong clutters. Moreover, the primal MERW has a critical limitation of strong bias to the most salient small target. To achieve multiple small targets detection, we develop a hierarchical version of the MERW method. Based on the hierarchical MERW (HMERW), we propose a small target detection method as follows. First, filtering technique is used to smooth the infrared image. Second, an output map is obtained by importing the filtered image into the HMERW. Then, a coefficient map is constructed to fuse the stationary dirtribution map of the HMERW. Finally, an adaptive threshold is used to segment multiple small targets from the fusion map. Extensive experiments on practical data sets demonstrate that the proposed method is superior to the state-of-the-art methods in terms of target enhancement, background suppression and multiple small targets detection.


title: Identities among higher genus modular graph tensors

abstract: Higher genus modular graph tensors map Feynman graphs to functions on the Torelli space of genus-$h$ compact Riemann surfaces which transform as tensors under the modular group $Sp(2h , \mathbb Z)$, thereby generalizing a construction of Kawazumi. An infinite family of algebraic identities between one-loop and tree-level modular graph tensors are proven for arbitrary genus and arbitrary tensorial rank. We also derive a family of identities that apply to modular graph tensors of higher loop order.


title: Multi-Resolution 3D Convolutional Neural Networks for Automatic Coronary Centerline Extraction in Cardiac CT Angiography Scans

abstract: We propose a deep learning-based automatic coronary artery tree centerline tracker (AuCoTrack) extending the vessel tracker by Wolterink (arXiv:1810.03143). A dual pathway Convolutional Neural Network (CNN) operating on multi-scale 3D inputs predicts the direction of the coronary arteries as well as the presence of a bifurcation. A similar multi-scale dual pathway 3D CNN is trained to identify coronary artery endpoints for terminating the tracking process. Two or more continuation directions are derived based on the bifurcation detection. The iterative tracker detects the entire left and right coronary artery trees based on only two ostium landmarks derived from a model-based segmentation of the heart. The 3D CNNs were trained on a proprietary dataset consisting of 43 CCTA scans. An average sensitivity of 87.1% and clinically relevant overlap of 89.1% was obtained relative to a refined manual segmentation. In addition, the MICCAI 2008 Coronary Artery Tracking Challenge (CAT08) training and test datasets were used to benchmark the algorithm and to assess its generalization. An average overlap of 93.6% and a clinically relevant overlap of 96.4% were obtained. The proposed method achieved better overlap scores than the current state-of-the-art automatic centerline extraction techniques on the CAT08 dataset with a vessel detection rate of 95%.


title: Simulation of positive streamers in CO$_2$ and in air: the role of photoionization or other electron sources

abstract: Positive streamer discharges have been studied and modelled extensively in air. Here we study positive streamers in CO$_2$ with and without oxygen admixtures; they are relevant for current high voltage technology as well as for discharges in the atmosphere of Venus. We discuss that no efficient photoionization mechanism is known for gases with a large CO$_2$ fraction, as photons in the relevant energy range are rapidly absorbed. Hence positive streamers can propagate only due to some other source of free electrons ahead of the ionization front. Therefore we study positive streamer propagation in CO$_2$ with different levels of background ionization to provide these free electrons. The effect of replacing photoionization by background ionization is studied with simulations in air. Simulating streamers in background fields of 16 to 20~kV/cm at standard temperature and pressure within a gap of 6.4~cm, we find that streamer propagation is rather insensitive to the level of photoionization or background ionization. We also discuss that the results depend not only on the value of breakdown field and applied electric field, and on preionization or photoionization, but also on the electron mobility $\mu(E)$ and the effective ionization coefficient $\alpha_{\rm eff}(E)$, that are gas-dependent functions of the electron energy or the electric field.


title: Lie algebroids. Lie-Rinehart algebras

abstract: After recalling the notion of Lie algebroid, we construct these structures associated with contact forms or systems. We are then interested in particular classes of Lie Rinehart algebras.


title: Image-based underwater 3D reconstruction for Cultural Heritage: from image collection to 3D. Critical steps and considerations

abstract: Underwater Cultural Heritage (CH) sites are widely spread; from ruins in coastlines up to shipwrecks in deep. The documentation and preservation of this heritage is an obligation of the mankind, dictated also by the international treaties like the Convention on the Protection of the Underwater Cultural Her-itage which fosters the use of "non-destructive techniques and survey meth-ods in preference over the recovery of objects". However, submerged CH lacks in protection and monitoring in regards to the land CH and nowadays recording and documenting, for digital preservation as well as dissemination through VR to wide public, is of most importance. At the same time, it is most difficult to document it, due to inherent restrictions posed by the environ-ment. In order to create high detailed textured 3D models, optical sensors and photogrammetric techniques seems to be the best solution. This chapter dis-cusses critical aspects of all phases of image based underwater 3D reconstruc-tion process, from data acquisition and data preparation using colour restora-tion and colour enhancement algorithms to Structure from Motion (SfM) and Multi-View Stereo (MVS) techniques to produce an accurate, precise and complete 3D model for a number of applications.


title: A Deep-Unfolded Reference-Based RPCA Network For Video Foreground-Background Separation

abstract: Deep unfolded neural networks are designed by unrolling the iterations of optimization algorithms. They can be shown to achieve faster convergence and higher accuracy than their optimization counterparts. This paper proposes a new deep-unfolding-based network design for the problem of Robust Principal Component Analysis (RPCA) with application to video foreground-background separation. Unlike existing designs, our approach focuses on modeling the temporal correlation between the sparse representations of consecutive video frames. To this end, we perform the unfolding of an iterative algorithm for solving reweighted $\ell_1$-$\ell_1$ minimization; this unfolding leads to a different proximal operator (a.k.a. different activation function) adaptively learned per neuron. Experimentation using the moving MNIST dataset shows that the proposed network outperforms a recently proposed state-of-the-art RPCA network in the task of video foreground-background separation.


title: The Bernardi formula for non-transitive deformations of the braid arrangement

abstract: Bernardi has given a general formula for the number of regions of a deformation of the braid arrangement as a signed sum over boxed trees. We prove that each set of boxed trees which share an underlying (rooted labeled plane) tree contributes 0, +1, or -1 to this sum, and we give an algorithm for computing this value. For Ish-type arrangements, we further construct a sign-reversing involution which reduces Bernardi's signed sum to the enumeration of a set of (rooted labeled plane) trees. We conclude by explicitly enumerating the trees corresponding to the regions of Ish-type arrangements which are nested, recovering their known counting formula.


title: Splay trees on trees

abstract: Search trees on trees (STTs) are a far-reaching generalization of binary search trees (BSTs), allowing the efficient exploration of tree-structured domains. (BSTs are the special case in which the underlying domain is a path.) Trees on trees have been extensively studied under various guises in computer science and discrete mathematics. Recently Bose, Cardinal, Iacono, Koumoutsos, and Langerman (SODA 2020) considered adaptive STTs and observed that, apart from notable exceptions, the machinery developed for BSTs in the past decades does not readily transfer to STTs. In particular, they asked whether the optimal STT can be efficiently computed or approximated (by analogy to Knuth's algorithm for optimal BSTs), and whether natural self-adjusting BSTs such as Splay trees (Sleator, Tarjan, 1983) can be extended to this more general setting. We answer both questions affirmatively. First, we show that a $(1 + \frac{1}{t})$-approximation of an optimal size-$n$ STT for a given search distribution can be computed in time $O(n^{2t + 1})$ for all integers $t \geq 1$. Second, we identify a broad family of STTs with linear rotation-distance, allowing the generalization of Splay trees to the STT setting. We show that our generalized Splay satisfies a static optimality theorem, asymptotically matching the cost of the optimal STT in an online fashion, i.e. without knowledge of the search distribution. Our results suggest an extension of the dynamic optimality conjecture for Splay trees to the broader setting of trees on trees.


title: Fibonacci-type orbifold data in Ising modular categories

abstract: An orbifold datum is a collection $\mathbb{A}$ of algebraic data in a modular fusion category $\mathcal{C}$. It allows one to define a new modular fusion category $\mathcal{C}_{\mathbb{A}}$ in a construction that is a generalisation of taking the Drinfeld centre of a fusion category. Under certain simplifying assumptions we characterise orbifold data $\mathbb{A}$ in terms of scalars satisfying polynomial equations and give an explicit expression which computes the number of isomorphism classes of simple objects in $\mathcal{C}_{\mathbb{A}}$. In Ising-type modular categories we find new examples of orbifold data which - in an appropriate sense - exhibit Fibonacci fusion rules. The corresponding orbifold modular categories have 11 simple objects, and for a certain choice of parameters one obtains the modular category for $sl(2)$ at level 10. This construction inverts the extension of the latter category by the $E_6$ commutative algebra.


title: "Cellular Network Densification Increases Radio-Frequency Pollution": True or False?

abstract: A very popular theory circulating among non-scientific communities claims that the massive deployment of Base Stations (BSs) over the territory, a.k.a. cellular network densification, always triggers an uncontrolled and exponential increase of human exposure to Radio Frequency "Pollution" (RFP). To face such concern in a way that can be understood by the layman, in this work we develop a very simple model to compute the RFP, based on a set of worst-case and conservative assumptions. We then provide closed-form expressions to evaluate the RFP variation in a pair of candidate 5G deployments, subject to different densification levels. Results, obtained over a wide set of representative 5G scenarios, dispel the myth: cellular network densification triggers an RFP decrease (up to three orders of magnitude) when the radiated power from the BS is adjusted to ensure a minimum sensitivity at the cell edge. Eventually, we analyze the conditions under which the RFP may increase when the network is densified (e.g., when the radiated power does not scale with the cell size), proving that the amount of RFP is always controlled. Finally, the results obtained by simulation confirm the outcomes of the RFP model.


title: Effective Waves for Random Three-dimensional Particulate Materials

abstract: How do you take a reliable measurement of a material whose microstructure is random? When using wave scattering, the answer is often to take an ensemble average (average over time or space). By ensemble averaging we can calculate the average scattered wave and the effective wavenumber. To date, the literature has focused on calculating the effective wavenumber for a plate filled with particles. One clear unanswered question was how to extend this approach to a material of any geometry and for any source. For example, does the effective wavenumber depend on only the microstructure, or also on the material geometry? In this work, we demonstrate that the effective wavenumbers depend on only microstructure and not the geometry, though beyond the long wavelength limit there are multiple effective wavenumbers. We show how to calculate the average wave scattered from a random particulate material of any shape, and for broad frequency ranges. As an example, we show how to calculate the average wave scattered from a sphere filled with particles.


title: Gravitational waves from binary black hole mergers surrounded by scalar field clouds: Numerical simulations and observational implications

abstract: We show how gravitational-wave observations of binary black hole (BBH) mergers can constrain the physical characteristics of a scalar field cloud parameterized by mass $\tilde{\mu}$ and strength $\phi_0$ that may surround them. We numerically study the inspiraling equal-mass, non-spinning BBH systems dressed in such clouds, focusing especially on the gravitational-wave signals emitted by their merger-ringdown phase. These waveforms clearly reveal that larger values of $\tilde{\mu}$ or $\phi_0$ cause bigger changes in the amplitude and frequency of the scalar-field-BBH ringdown signals. We show that the numerical waveforms of scalar-field-BBHs can be modelled as chirping sine-Gaussians, with matches in excess of 95%. This observation enables one to employ computationally expensive Bayesian studies for estimating the parameters of such binaries. Using our chirping sine-Gaussian signal model we establish that observations of BBH mergers at a distance of 450 Mpc will allow to distinguish BBHs without any scalar field from those with a field strength $\phi_0 \gtrsim 5.5\times 10^{-3}$, at any fixed value of $\tilde \mu \in [0.3,0.8]$, with 90% confidence or better, in single detectors with Advanced LIGO/Virgo type sensitivities. This provides hope for the possibility of determining or constraining the mass of ultra-light bosons with gravitational-wave observations of BBH mergers.


title: Turbulence characteristics of switchbacks and non-switchbacks intervals observed by \emph{Parker Solar Probe}

abstract: We use \emph{Parker Solar Probe} (\emph{PSP}) in-situ measurements to analyze the characteristics of solar wind turbulence during the first solar encounter covering radial distances between $35.7R_\odot$ and $41.7R_\odot$. In our analysis we isolate so-called switchback (SB) intervals (folded magnetic field lines) from non-switchback (NSB) intervals, which mainly follow the Parker spiral field. Using a technique based on conditioned correlation functions, we estimate the power spectra of Elsasser, magnetic and bulk velocity fields separately in the SB and NSB intervals. In comparing the turbulent energy spectra of the two types of intervals, we find the following characteristics: 1) The decorrelation length of the backward-propagating Elsasser field $z^-$ is larger in the NSB intervals than the one in the SB intervals; 2) the magnetic power spectrum in SB intervals is steeper, with spectral index close to -5/3, than in NSB intervals, which have a spectral index close to -3/2; 3) both SB and NSB turbulence are imbalanced with NSB having the largest cross-helicity, 4) the residual energy is larger in the SB intervals than in NSB, and 5) the analyzed fluctuations are dominated by Alfv\'enic fluctuations that are propagating in the \emph{sunward} (\emph{anti-sunward}) direction for the SB (NSB) turbulence. These observed features provide further evidence that the switchbacks observed by \emph{PSP} are associated with folded magnetic field lines giving insight into their turbulence nature.


title: Decremental APSP in Directed Graphs Versus an Adaptive Adversary

abstract: Given a directed graph $G = (V,E)$, undergoing an online sequence of edge deletions with $m$ edges in the initial version of $G$ and $n = |V|$, we consider the problem of maintaining all-pairs shortest paths (APSP) in $G$. Whilst this problem has been studied in a long line of research [ACM'81, FOCS'99, FOCS'01, STOC'02, STOC'03, SWAT'04, STOC'13] and the problem of $(1+\epsilon)$-approximate, weighted APSP was solved to near-optimal update time $\tilde{O}(mn)$ by Bernstein [STOC'13], the problem has mainly been studied in the context of oblivious adversaries, which assumes that the adversary fixes the update sequence before the algorithm is started. In this paper, we make significant progress on the problem in the setting where the adversary is adaptive, i.e. can base the update sequence on the output of the data structure queries. We present three new data structures that fit different settings: We first present a deterministic data structure that maintains exact distances with total update time $\tilde{O}(n^3)$. We also present a deterministic data structure that maintains $(1+\epsilon)$-approximate distance estimates with total update time $\tilde O(\sqrt{m} n^2/\epsilon)$ which for sparse graphs is $\tilde O(n^{2+1/2}/\epsilon)$. Finally, we present a randomized $(1+\epsilon)$-approximate data structure which works against an adaptive adversary; its total update time is $\tilde O(m^{2/3}n^{5/3} + n^{8/3}/(m^{1/3}\epsilon^2))$ which for sparse graphs is $\tilde O(n^{2+1/3})$. Our exact data structure matches the total update time of the best randomized data structure by Baswana et al. [STOC'02] and maintains the distance matrix in near-optimal time. Our approximate data structures improve upon the best data structures against an adaptive adversary which have $\tilde{O}(mn^2)$ total update time [JACM'81, STOC'03].


title: Structural parts as quadrics -- elasticity ellipses revisited

abstract: Elasticity ellipses or central ellipses have been long used in graphic statics to capture the elastic behaviour of structural elements. The paper gives a generalisation the concept both in dimensions and in the possibility of degenerate conics / quadrics. The effect of projective transformations of these quadrics is also given, such that the entire mechanical system can be transformed preserving equilibrium and compatibility between its elements.


title: Families of curves orthogonal to the lines $y=mx-2m-m^3$

abstract: The family of lines $y=mx-2m-m^3$, are well known to be normal to the parabola $y^2=4x$. However, this family of lines is normal to a family of curves of which this parabola is just one member. Here, by solving an interesting first order and third degree ODE, we bring out these curves. The resulting one set of curves are "parabola-like" but non-standard ones and the other family is not even "parabola like".


title: Universal evolutionary model for periodical species

abstract: Real-world examples of periods of periodical organisms range from cicadas whose life-cycles are larger prime numbers, like 13 or 17, to bamboos whose periods are large multiples of small primes, like 40 or even 120. The periodicity is caused by interaction of species, be it a predator-prey relationship, symbiosis, commensialism, or competition exclusion principle. We propose a simple mathematical model which explains and models all those principles, including listed extremel cases. This, rather universal, qualitative model is based on the concept of a local fitness function, where a randomly chosen new period is selected if the value of the global fitness function of the species increases. Arithmetically speaking, the different observed interactions are related to only four principles: given a couple of integer periods either (1) their greatest common divisor is one, (2) one of the periods is prime, (3) both periods are equal, or (4) one period is an integer multiple of the other.


title: Airy gas model: From three to reduced dimensions

abstract: By using the propagator of linear potential as a main tool, we extend the Airy gas model, originally developed for the three-dimensional ($d=3$) edge electron gas, to systems in reduced dimensions ($d=2,1$). First, we derive explicit expressions for the edge particle density and the corresponding kinetic energy density (KED) of the Airy gas model in all dimensions. The densities are shown to obey the local virial theorem. We obtain a functional relationship between the positive KED and the particle density and its gradients and analyze the results inside the bulk as a limit of the local-density approximation. We show that in this limit the KED functional reduces to that of the Thomas-Fermi model in $d$ dimensions.


title: Virtual Smartphone: High Fidelity Interaction with Proxy Objects in Virtual Reality

abstract: This workshop paper presents two proxy objects for high fidelity interaction in virtual reality (VR): a paper map and a smartphone. We showcase how our virtual paper map can increase interactivity and orientation, while our virtual smartphone extends the use of a proxy object, as it allows for actual touch input on a real phone leading to an almost infinite set of possible (inter-)actions (e.g. snapping pictures in the virtual world). Observations showed that participants were very precise in holding and interacting with both the paper map and the smartphone even though they did not see their hands in VR. The interaction in general was very intuitive which was mostly attributed to the realistic size of the virtual objects. Using our findings we discuss the trade off between adaptivity and high fidelity of proxy objects in VR.


title: PMSD: Data-Driven Simulation Using System Dynamics and Process Mining

abstract: Process mining extends far beyond process discovery and conformance checking, and also provides techniques for bottleneck analysis and organizational mining. However, these techniques are mostly backward-looking. PMSD is a web application tool that supports forward-looking simulation techniques. It transforms the event data and process mining results into a simulation model which can be executed and validated. PMSD includes log transformation, time window selection, relation detection, interactive model generation, simulating and validating the models in the form of system dynamics, i.e., a technique for an aggregated simulation. The results of the modules are visualized in the tool for a better interpretation


title: Study of autonomous conservative oscillator using an improved perturbation method

abstract: In a recent article \cite{manimegalai2019}, Aboodh transform based homotopy perturbation method ($AT$) has been found to produce approximate analytical solutions in a simple way but with better accuracy in comparison to those obtained from some of the established approximation methods \cite{mehdipour2010application,nofal2013analytical} for some physically relevant anharmonic oscillators such as autonomous conservative oscillator (ACO). In the present article, expansion of frequency ($\omega$) and an auxiliary parameter ($h$) are incorporated in the framework of the homotopy perturbation method (HPM) to improve the accuracy by retaining its simplicity. Laplace transform is used to make the calculation simpler. This improved HPM ($LH$) is simple but provides highly accurate results for ACO in comparison to those obtained from $AT$. The error in the values of frequency and displacement calculated using the $LH$ is found to be one or two order of magnitude less than those obtained from $AT$ for the considered parameter sets.


title: Elementary characters on semigroups: the rational case

abstract: Since polynomials form a subsemigroup of the semigroup of rational functions, every character on rational functions is a character on polynomials. On the other direction, not every character on polynomials is the restriction of a character on rational functions. What are the characters on polynomials that can be extended to rational functions? In this work, we conjecture that the only characters that can be extended are those that depends on the degree, often called elementary. Also, we construct two example of character on polynomials, not elementaries, that cannot be extended to rational functions.


title: $\mathcal{N}=1$ Super-Yang-Mills theory on the lattice with twisted mass fermions

abstract: Super-Yang-Mills theory (SYM) is a central building block for supersymmetric extensions of the Standard Model of particle physics. Whereas the weakly coupled subsector of the latter can be treated within a perturbative setting, the strongly coupled subsector must be dealt with a non-perturbative approach. Such an approach is provided by the lattice formulation. Unfortunately a lattice regularization breaks supersymmetry and consequently the mass degeneracy within a supermultiplet. In this article we investigate the properties of $\mathcal{N}=1$ supersymmetric SU(3) Yang-Mills theory with a lattice Wilson Dirac operator with an additional parity mass, similar as in twisted mass lattice QCD. We show that a special $45^\circ$ twist effectively moves the mass splitting of the chiral partners. Thus, at finite lattice spacing both chiral and supersymmetry are enhanced resulting in an improved continuum extrapolation. Furthermore, we show that for the non-interacting theory at $45^\circ$ twist discretization errors of order $\mathcal{O}(a)$ are suppressed, suggesting that the same happens for the interacting theory as well. As an aside, we demonstrate that the DD$\alpha$AMG multigrid algorithm accelerates the inversion of the Wilson Dirac operator considerably. On a $16^3\times 32$ lattice, speed-up factors of up to 20 are reached if commonly used algorithms are replaced by the DD$\alpha$AMG.


title: MGD-GAN: Text-to-Pedestrian generation through Multi-Grained Discrimination

abstract: In this paper, we investigate the problem of text-to-pedestrian synthesis, which has many potential applications in art, design, and video surveillance. Existing methods for text-to-bird/flower synthesis are still far from solving this fine-grained image generation problem, due to the complex structure and heterogeneous appearance that the pedestrians naturally take on. To this end, we propose the Multi-Grained Discrimination enhanced Generative Adversarial Network, that capitalizes a human-part-based Discriminator (HPD) and a self-cross-attended (SCA) global Discriminator in order to capture the coherence of the complex body structure. A fined-grained word-level attention mechanism is employed in the HPD module to enforce diversified appearance and vivid details. In addition, two pedestrian generation metrics, named Pose Score and Pose Variance, are devised to evaluate the generation quality and diversity, respectively. We conduct extensive experiments and ablation studies on the caption-annotated pedestrian dataset, CUHK Person Description Dataset. The substantial improvement over the various metrics demonstrates the efficacy of MGD-GAN on the text-to-pedestrian synthesis scenario.


title: Causal coupling inference from multivariate time series based on ordinal partition transition networks

abstract: Identifying causal relationships is a challenging yet crucial problem in many fields of science like epidemiology, climatology, ecology, genomics, economics and neuroscience, to mention only a few. Recent studies have demonstrated that ordinal partition transition networks (OPTNs) allow inferring the coupling direction between two dynamical systems. In this work, we generalize this concept to the study of the interactions among multiple dynamical systems and we propose a new method to detect causality in multivariate observational data. By applying this method to numerical simulations of coupled linear stochastic processes as well as two examples of interacting nonlinear dynamical systems (coupled Lorenz systems and a network of neural mass models), we demonstrate that our approach can reliably identify the direction of interactions and the associated coupling delays. Finally, we study real-world observational microelectrode array electrophysiology data from rodent brain slices to identify the causal coupling structures underlying epileptiform activity. Our results, both from simulations and real-world data, suggest that OPTNs can provide a complementary and robust approach to infer causal effect networks from multivariate observational data.


title: Quantum stochastic series expansion methods

abstract: A quantum implementation of the Stochastic Series Expansion (SSE) Monte Carlo method is proposed, and it is shown that quantum SSE offers significant advantages over classical implementations of SSE. In particular, for problems where classical SSE encounters the sign problem, the cost of implementing a Monte Carlo iteration scales only linearly with system size in quantum SSE, while it may scale exponentially with system size in classical SSE. In cases where classical SSE can be efficiently implemented, quantum SSE still offers an advantage by allowing for more general observables to be measured.


title: Regularized K-means through hard-thresholding

abstract: We study a framework of regularized $K$-means methods based on direct penalization of the size of the cluster centers. Different penalization strategies are considered and compared through simulation and theoretical analysis. Based on the results, we propose HT $K$-means, which uses an $\ell_0$ penalty to induce sparsity in the variables. Different techniques for selecting the tuning parameter are discussed and compared. The proposed method stacks up favorably with the most popular regularized $K$-means methods in an extensive simulation study. Finally, HT $K$-means is applied to several real data examples. Graphical displays are presented and used in these examples to gain more insight into the datasets.


title: Coupled Oscillatory Recurrent Neural Network (coRNN): An accurate and (gradient) stable architecture for learning long time dependencies

abstract: Circuits of biological neurons, such as in the functional parts of the brain can be modeled as networks of coupled oscillators. Inspired by the ability of these systems to express a rich set of outputs while keeping (gradients of) state variables bounded, we propose a novel architecture for recurrent neural networks. Our proposed RNN is based on a time-discretization of a system of second-order ordinary differential equations, modeling networks of controlled nonlinear oscillators. We prove precise bounds on the gradients of the hidden states, leading to the mitigation of the exploding and vanishing gradient problem for this RNN. Experiments show that the proposed RNN is comparable in performance to the state of the art on a variety of benchmarks, demonstrating the potential of this architecture to provide stable and accurate RNNs for processing complex sequential data.


title: Distributed Proximal Splitting Algorithms with Rates and Acceleration

abstract: We analyze several generic proximal splitting algorithms well suited for large-scale convex nonsmooth optimization. We derive sublinear and linear convergence results with new rates on the function value suboptimality or distance to the solution, as well as new accelerated versions, using varying stepsizes. In addition, we propose distributed variants of these algorithms, which can be accelerated as well. While most existing results are ergodic, our nonergodic results significantly broaden our understanding of primal-dual optimization algorithms.


title: Absolute X-ray energy measurement using a high-accuracy angle encoder

abstract: This paper presents an absolute X-ray photon energy measurement method that uses a Bond diffractometer. The proposed system enables the prompt and rapid in-situ measurement of photon energies in a wide energy range. The diffractometer uses a reference silicon single crystal plate and a highly accurate angle encoder called SelfA. We evaluate the performance of the system by repeatedly measuring the energy of the first excited state of the potassium-40 nuclide. The excitation energy is determined as 29829.39(6) eV. It is one order of magnitude more precise than the previous measurement. The estimated uncertainty of the photon energy measurement was 0.7 ppm as a standard deviation and the maximum observed deviation was 2 ppm.


title: Optimal optical polarization of nitrogen-vacancy center with arbitrary waveform pulse

abstract: The current work proposes a method for pulsed-light polarization of nitrogen-vacancy (NV) center electron spin. To evaluate the influence of pulsed spin polarization, we establish a polarization evaluation index based on polarizability and polarization time. Master equation model are utilized to theoretically calculate the spin polarization dynamics under light excitation and the optimal polarization conditions for the conventional methods are obtained. A novel pulsed-light polarization method is proposed by changing the optical pumping rate in the master equation from a fixed value to a time variable and an optimal waveform for proposed method is demonstrated through the variational method, which can simultaneously achieve high polarizability and requires a short polarization time. Hence, the polarization evaluation index is improved by ~10%. Moreover, the proposed method is verified by a pulsed-laser experimental system based on an arbitrary waveform generator. The current report shall expand the application horizon of NV center based quantum sensing.


title: Optimization of DMD-based independent amplitude and phase modulation: a spatial resolution and quantization

abstract: This paper presents the results of a comprehensive study on the optimization of the optical system for independent amplitude and phase wavefront manipulation, which is based on a digital micromirror device implementation. The parameters of generated binary fringe patterns to reach the optimal quality of target complex-valued wavefront have been found for several examples, and a general algorithm for DMD pattern optimization is described. It was shown that a trade-off between spatial resolution and quantization of the target amplitude and phase distribution should be achieved. The increase of carrier frequency results in higher spatial resolution of the modulated complex wave but decreases its quantization. The dependence of the best DMD-pattern generation parameters on the type of target complex wavefront is discussed in terms of spatial resolution and quantization degree.


title: On the Disformal Transformation of the Einstein-Hilbert Action

abstract: Disformal transformation is a generalisation of the well-known conformal transformation commonly elaborated in mainstream graduate texts in gravity (relativity) and modern cosmology. This transformation is one of the most important mathematical operations in scalar tensor theories attempting to address pressing problems involving dark energy and dark matter. With this topic yet to penetrate these texts, we present a pedagogically oriented derivation of the disformal transformation of the Einstein-Hilbert action. Along the way of calculation, we encounter apparently problematic terms that could be construed as leading to equations of motion that go beyond second order in derivatives, signalling instability. We demonstrate that these terms can be eliminated and absorbed through the definition of the Riemann curvature tensor. The transformed Einstein-Hilbert action turns out to be a special case of the Horndeski action and the equations of motion for the scalar field that it describes are all up to second order only in derivatives, implying stability.


title: Estimands in Hematologic Oncology Trials

abstract: The estimand framework included in the addendum to the ICH E9 guideline facilitates discussions to ensure alignment between the key question of interest, the analysis, and interpretation. Therapeutic knowledge and drug mechanism play a crucial role in determining the strategy and defining the estimand for clinical trial designs. Clinical trials in patients with hematological malignancies often present unique challenges for trial design due to complexity of treatment options and existence of potential curative but highly risky procedures, e.g. stem cell transplant or treatment sequence across different phases (induction, consolidation, maintenance). Here, we illustrate how to apply the estimand framework in hematological clinical trials and how the estimand framework can address potential difficulties in trial result interpretation. This paper is a result of a cross-industry collaboration to connect the International Conference on Harmonisation (ICH) E9 addendum concepts to applications. Three randomized phase 3 trials will be used to consider common challenges including intercurrent events in hematologic oncology trials to illustrate different scientific questions and the consequences of the estimand choice for trial design, data collection, analysis, and interpretation. Template language for describing estimand in both study protocols and statistical analysis plans is suggested for statisticians' reference.


title: Identification of images of COVID-19 from Chest Computed Tomography (CT) images using Deep learning: Comparing COGNEX VisionPro Deep Learning 1.0 Software with Open Source Convolutional Neural Networks

abstract: For testing patients infected with COVID-19, along with RT-PCR testing, chest radiology images are being used. For the detection of COVID-19 from radiology images, many organizations are proposing the use of Deep Learning. University of Waterloo and DarwinAI, have designed their own Deep Learning model COVIDNet-CT to detect COVID-19 from infected chest CT images. Additionally, they have introduced a CT image dataset COVIDx-CT, from CT images collected by the China National Center for Bioinformation. COVIDx-CT contains 104,009 CT image slices across 1,489 patient cases. After obtaining remarkable results on the identification of COVID-19 from chest X-ray images by using the COGNEX VisionPro Deep Learning Software 1.0 this time we test the performance of the software on the identification of COVID-19 from CT images. COGNEX Deep Learning Software: VisionPro Deep Learning, is a Deep Learning software that is used across various domains ranging from factory automation to life sciences. In this study, we train the classification model on 82,818 chest CT training and validation images from the COVIDx-CT dataset in 3 classes - normal, pneumonia, and COVID-19 and then test the results of the classification on the 21,191 test images are compared with the results of COVIDNet-CT and various other state of the art Deep Learning models from the open-source community. Also, we test how reducing the number of images in the training set effects the results of the software. Overall, VisionPro Deep Learning gives the best results with F-scores over 99%, even as the number of images in the training set is reduced significantly. This software is by no means a stand-alone solution in the detection of COVID-19 but can aid radiologists and clinicians in achieving faster and understandable diagnosis using the full potential of Deep Learning, without the prerequisite of having to code in any programming language.


title: Shielding the vulnerable in an epidemic: a numerical approach

abstract: The death toll for Covid-19 may be reduced by dividing the population into two classes, the vulnerable and the fit, with different lockdown regimes. Instead of one reproduction number there now are four parameters. These make it possible to quantify the effect of the social distancing measures. There is a simple stochastic model for epidemics in a two type population. Apart from the size of the population of the vulnerable and the fit, and the initial number of infected in the two classes, only the four reproduction parameters are needed to run the two type Reed-Frost model. The program is simple and fast. On a pc it takes less than five minutes to do a hundred thousand simulations of the epidemic for a population of the size of the US. Epidemics are non-linear processes. Results may be counterintuitive. The average number of vulnerable persons infected by an infectious fit person is a crucial parameter of the epidemic in the two type population. Intuitively this parameter should be small. However simulations show that even if this parameter is small the death toll may be higher than without shielding. Under certain conditions increasing the value of the parameter may reduce the death toll. The article addresses these blind spots in our intuition.


title: Robust Output Regulation of the Linearized Boussinesq Equations with Boundary Control and Observation

abstract: We study temperature and velocity output tracking problem for a two-dimensional room model with the fluid dynamics governed by the linearized translated Boussinesq equations. Additionally, the room model includes finite-dimensional models for actuation and sensing dynamics, thus the complete model dynamics are governed by an ODE-PDE-ODE system. As the main result, we design a low-dimensional internal model based controller for robust output racking of the room model. Efficiency of the controller is demonstrated through a numerical example of velocity and temperature tracking.


title: Deep Learning Based Computer-Aided Systems for Breast Cancer Imaging : A Critical Review

abstract: This paper provides a critical review of the literature on deep learning applications in breast tumor diagnosis using ultrasound and mammography images. It also summarizes recent advances in computer-aided diagnosis (CAD) systems, which make use of new deep learning methods to automatically recognize images and improve the accuracy of diagnosis made by radiologists. This review is based upon published literature in the past decade (January 2010 January 2020). The main findings in the classification process reveal that new DL-CAD methods are useful and effective screening tools for breast cancer, thus reducing the need for manual feature extraction. The breast tumor research community can utilize this survey as a basis for their current and future studies.


title: Containing a pandemic: Nonpharmaceutical interventions and the "second wave"

abstract: In response to the worldwide outbreak of the coronavirus disease COVID-19, a variety of nonpharmaceutical interventions such as face masks and social distancing have been implemented. A careful assessment of the effects of such containment strategies is required to avoid exceeding social and economical costs as well as a dangerous "second wave" of the pandemic. In this work, we combine a recently developed dynamical density functional theory model and an extended SIRD model with hysteresis to study effects of various measures and strategies using realistic parameters. Depending on intervention thresholds, a variety of phases with different numbers of shutdowns and deaths are found. Spatiotemporal simulations provide further insights into the dynamics of a second wave. Our results are of crucial importance for public health policy.


title: On Differentially Algebraic Generating Series for Walks in the Quarter Plane

abstract: We refine necessary and sufficient conditions for the generating series of a weighted model of a quarter plane walk to be differentially algebraic. In addition, we give algorithms based on the theory of Mordell-Weil lattices, that, for each weighted model, yield polynomial conditions on the weights determining this property of the associated generating series.


title: Augmenting Machine Learning with Information Retrieval to Recommend Real Cloned Code Methods for Code Completion

abstract: Software developers frequently reuse source code from repositories as it saves development time and effort. Code clones accumulated in these repositories hence represent often repeated functionalities and are candidates for reuse in an exploratory or rapid development. In previous work, we introduced DeepClone, a deep neural network model trained by fine tuning GPT-2 model over the BigCloneBench dataset to predict code clone methods. The probabilistic nature of DeepClone output generation can lead to syntax and logic errors that requires manual editing of the output for final reuse. In this paper, we propose a novel approach of applying an information retrieval (IR) technique on top of DeepClone output to recommend real clone methods closely matching the predicted output. We have quantitatively evaluated our strategy, showing that the proposed approach significantly improves the quality of recommendation.


title: Spreading height and critical conditions for the collapse of turbulent fountains in stratified media

abstract: Axisymmetric fountains in stratified environments rise until reaching a maximum height, where the vertical momentum vanishes, and then falls and spread radially as an annular plume following a well-known top-hat profile. Here, firstly, we generalize the model of Morton et al. (Proc. R. Soc. Lond. A \textbf{234}, 1, 1956), in order to correctly determine the dependence of the maximum height and the spreading height with the parameters involved. We obtain the critical conditions for the collapse of the fountain, \textit i.e. when the jet falls up to the source level, and show that the spreading height must be expressed as a function of at least two parameters. To improve the quantitative agreement with the experiments we modify the criterion to take the mixing process in the down flow into account. Numerical simulations were implemented to estimate the parameter values that characterizes this merging. We show that our generalized model agrees very well with the experimental measurements.


title: Surface Treatment of Cu:NiOx Hole-Transporting Layer Using \b{eta}-Alanine for Hysteresis-Free and Thermally Stable Inverted Perovskite Solar Cells

abstract: Inverted perovskite solar cells (PSCs) using a Cu:NiOx hole transporting layer (HTL) often exhibit stability issues and in some cases J/V hysteresis. In this work, we developed a \b{eta}-alanine surface treatment process on Cu:NiOx HTL that provides J/V hysteresis-free, highly efficient, and thermally stable inverted PSCs. The improved device performance due to \b{eta}-alanine-treated Cu:NiOx HTL is attributed to the formation of an intimate Cu:NiOx/perovskite interface and reduced charge trap density in the bulk perovskite active layer. The \b{eta}-alanine surface treatment process on Cu:NiOx HTL eliminates major thermal degradation mechanisms, providing 40 times increased lifetime performance under accelerated heat lifetime conditions. By using the proposed surface treatment, we report optimized devices with high power conversion efficiency (PCE) (up to 15.51%) and up to 1000 h lifetime under accelerated heat lifetime conditions (60 C, N2).


title: Efficient Estimation of Graph Trussness

abstract: A $k$-truss is an edge-induced subgraph $H$ such that each of its edges belongs to at least $k-2$ triangles of $H$. This notion has been introduced around ten years ago in social network analysis and security, as a form of cohesive subgraph that is rich of triangles and less stringent than the clique. The \emph{trussness} of a graph is the maximum $k$ such that a $k$-truss exists. The problem of computing $k$-trusses has been largely investigated from the practical and engineering point of view. On the other hand, the theoretical side of the problem has received much less attention, despite presenting interesting challenges. The existing methods share a common design, based on iteratively removing the edge with smallest support, where the support of an edge is the number of triangles containing it. The aim of this paper is studying algorithmic aspects of graph trussness. While it is possible to show that the time complexity of computing exactly the graph trussness and that of counting/listing all triangles is inherently the same, we provide efficient algorithms for estimating its value, under suitable conditions, with significantly lower complexity than the exact approach. In particular, we provide a $(1 \pm \epsilon)$-approximation algorithm that is asymptotically faster than the exact approach, on graphs which contain $\omega(m \, polylog(n))$ triangles, and has the same running time on graphs that do not. For the latter case, we also show that it is impossible to obtain an approximation algorithm with faster running time than the one of the exact approach when the number of triangles is $O(m)$, unless well known conjectures on triangle-freeness and Boolean matrix multiplication are false.