title: Neuronal Network Inference and Membrane Potential Model using Multivariate Hawkes Processes

abstract: In this work, we propose to catch the complexity of the membrane potential's dynamic of a motoneuron between its spikes, taking into account the spikes from other neurons around. Our approach relies on two types of data: extracellular recordings of multiple spikes trains and intracellular recordings of the membrane potential of a central neuron. Our main contribution is to provide a unified framework and a complete pipeline to analyze neuronal activity from data extraction to statistical inference. The first step of the procedure is to select a subnetwork of neurons impacting the central neuron: we use a multivariate Hawkes process to model the spike trains of all neurons and compare two sparse inference procedures to identify the connectivity graph. Then we infer a jump-diffusion dynamic in which jumps are driven from a Hawkes process, the occurrences of which correspond to the spike trains of the aforementioned subset of neurons that interact with the central neuron. We validate the Hawkes model with a goodness-of-fit test and we show that taking into account the information from the connectivity graph improves the inference of the jump-diffusion process. The entire code has been developed and is freely available on GitHub.


title: Image-based scene recognition for robot navigation considering traversable plants and its manual annotation-free training

abstract: This paper describes a method of estimating the traversability of plant parts covering a path and navigating through them for mobile robots operating in plant-rich environments. Conventional mobile robots rely on scene recognition methods that consider only the geometric information of the environment. Those methods, therefore, cannot recognize paths as traversable when they are covered by flexible plants. In this paper, we present a novel framework of image-based scene recognition to realize navigation in such plant-rich environments. Our recognition model exploits a semantic segmentation branch for general object classification and a traversability estimation branch for estimating pixel-wise traversability. The semantic segmentation branch is trained using an unsupervised domain adaptation method and the traversability estimation branch is trained with label images generated from the robot's traversal experience during the data acquisition phase, coined traversability masks. The training procedure of the entire model is, therefore, free from manual annotation. In our experiment, we show that the proposed recognition framework is capable of distinguishing traversable plants more accurately than a conventional semantic segmentation with traversable plant and non-traversable plant classes, and an existing image-based traversability estimation method. We also conducted a real-world experiment and confirmed that the robot with the proposed recognition method successfully navigated in plant-rich environments.


title: BezierSeg: Parametric Shape Representation for Fast Object Segmentation in Medical Images

abstract: Delineating the lesion area is an important task in image-based diagnosis. Pixel-wise classification is a popular approach to segmenting the region of interest. However, at fuzzy boundaries such methods usually result in glitches, discontinuity, or disconnection, inconsistent with the fact that lesions are solid and smooth. To overcome these undesirable artifacts, we propose the BezierSeg model which outputs bezier curves encompassing the region of interest. Directly modelling the contour with analytic equations ensures that the segmentation is connected, continuous, and the boundary is smooth. In addition, it offers sub-pixel accuracy. Without loss of accuracy, the bezier contour can be resampled and overlaid with images of any resolution. Moreover, a doctor can conveniently adjust the curve's control points to refine the result. Our experiments show that the proposed method runs in real time and achieves accuracy competitive with pixel-wise segmentation models.


title: Monitoring Gamma-Ray Burst VHE emission with the Southern Wide-field-of-view Gamma-ray Observatory

abstract: It has been established that Gamma-Ray Bursts (GRB) can produce Very High Energy radiation (E > 100 GeV), opening a new window on the investigation of particle acceleration and radiation properties in the most energetic domain. We expect that next-generation instruments, such as the Cherenkov Telescope Array (CTA), will mark a huge improvement in their observation. However, constraints on the target visibility and the limited duty cycle of Imaging Atmospheric Cherenkov Telescopes (IACT) reduce their ability to react promptly to transient events and to characterise their general properties. Here we show that an instrument based on the Extensive Air Shower (EAS) array concept, proposed by the Southern Wide Field-of-view Gamma-ray Observatory (SWGO) Collaboration, has promising possibilities to detect and track VHE emission from GRBs. Observations made by the Fermi Large Area Telescope (Fermi-LAT) identified some events with a distinct spectral component, extending above $1\,$GeV or even $10\,$GeV, which can represent a substantial fraction of the emitted energy and also arise in early stages of the process. Using models based on these properties, we estimate the possibilities that a wide field of view and large effective area ground-based monitoring facility has to probe VHE emission from GRBs. We show that the ability to monitor VHE transients with a nearly continuous scanning of the sky grants an opportunity to access simultaneous electromagnetic counterparts to Multi-Messenger triggers up to cosmological scales, in a way that is not available to IACTs.


title: Redundancy Resolution in Kinematic Control of Serial Manipulators in Multi-Obstacle Environment

abstract: The paper focuses on the redundancy resolution in kinematic control of a new type of serial manipulator composed of multiple tensegrity segments, which are moving in a multi-obstacle environment. The general problem is decomposed into two sub-problems, which deal with collision-free path planning for the robot end-effector and collision-free motion planning for the robot body. The first of them is solved via discrete dynamic programming, the second one is worked out using quadratic programming with mixed linear equality/nonequality constraints. Efficiency of the proposed technique is confirmed by simulation.


title: Zigzag nanoribbon of gated bilayer hexagonal crystals with spontaneous edge magnetism

abstract: Zigzag nanoribbons of monolayer graphene-like two-dimensional materials host spontaneous edge magnetism at the zigzag terminations, whose configuration controls the band gap. In this article, the edge magnetism of zigzag nanoribbons of bilayer hexagonal crystals are studied. The specific models of bilayer graphene and bilayer silicene are both studied. As the gated voltage increases, the total energy level, magnetic structure and band structures of the ground state and the first quasi-stable excited state are tuned. For certain region of the gated voltage, the band gaps of spin up and down are opened and closed, respectively, so that the systems are in the spin-polarized metallic phase. The spin-polarized conducting edge band of the bilayer nanoribbons could be applied for gate tunable spintronic devices.


title: More ATP Does Not Equal More Contractility: Power And Remodelling In Reconstituted Actomyosin

abstract: The cytoskeletal component actomyosin is a canonical example of active matter since the powerstroke cycle locally converts chemical energy in the form of adenoside triphosphate (ATP) into mechanical work for remodelling. Observing myosin II minifilaments as they remodel actin {\it in vitro}, we now report that, at high concentrations of ATP, myosin minifilaments form meta-stable swirling patterns that are characterised by recurrent vortex- and spiral-like motifs, whereas at low concentrations of ATP, such structures give way to aster-like patterns. This is shown to be consistent with the (quasi-)steady states of a polar active hydrodynamic theory of actomyosin whose ATP-scaling is obtained from a microscopic description for the ATP-dependent binding statistics of the headgroups of single myosin II minifilaments. The latter codifies the heuristic that, since the powerstroke cycle involves the unbinding of myosin II headgroups from actin, increases in the concentration of ATP reduce the likelihood that a given myosin II minifilament has more than one headgroup bound simultaneously, reducing its ability to generate contractile forces and increasing the relative likelihood of processive motion. As a result, our work demonstrates how ATP not only controls the rate at which work is done -- \textit{i.e.,} the power -- but also the mode by which this occurs.


title: PID-like active control strategy for electroacoustic resonators to design tunable single-degree-of-freedom sound absorbers

abstract: Sound absorption at low frequencies still remains a challenge in both scientific research and engineering practice. Natural porous materials are ineffective in this frequency range, as well as acoustic resonators which present too narrow bandwidth of absorption, thus requiring alternative solutions based on active absorption techniques. In the present work, we propose an active control framework applied on a closed-box loudspeaker to enable the adjustment of the acoustic impedance at the loudspeaker diaphragm. More specifically, based on the proportionality between the pressure inside the enclosure and the axial displacement of the loudspeaker diaphragm at low frequencies, we demonstrate both analytically and experimentally that a PID-like feedback control approach allows tuning independently the compliance, the resistance and the moving mass of the closed-box loudspeaker to implement a prescribed impedance of a single-degree-of-freedom resonator. By considering different control combinations to tailor the resonator characteristics, a perfect absorption (with absorption coefficient equal to 1) is achievable at the target resonance frequency, while enlarging the effective absorption bandwidth. Moreover, the proposed feedback control strategy shows an excellent control accuracy, especially compared to the feedforward-based control formerly reported in the literature. The mismatches between the performance of experimental prototype and the model, likely to result from the control time delay and the inaccuracy in estimating the loudspeaker parameters, can be compensated directly by tuning the control parameters in the control platform. The active resonators implemented through the reported control scheme can be used to build more complex acoustic devices/structures to enable high-efficiency broadband sound absorption or other types of acoustic phenomena such as wavefront shaping.


title: First Law of Black Hole in the Gravitational Electromagnetic System

abstract: After considering the quantum corrections of Einstein-Maxwell theory, the effective theory will contain some higher-curvature terms and nonminimally coupled electromagnetic fields. In this paper, we study the first law of black holes in the gravitational electromagnetic system with the Lagrangian $\math{L}(g_{ab}, R_{abcd}, F_{ab})$. Firstly, we calculate the Noether charge and the variational identity in this theory, and then generically derive the first law of thermodynamics for an asymptotically flat stationary axisymmetrical symmetric black hole without the requirement that the electromagnetic field is smooth on the bifurcation surface. Our results indicate that the first law of black hole thermodynamics might be valid for the Einstein-Maxwell theory with some quantum corrections in the effective region.


title: Projective properties of Divergence-free symmetric tensors, and new dispersive estimates in gas dynamics

abstract: The class of Divergence-free symmetric tensors is ubiquitous in Continuum Mechanics. We show its invariance under projective transformations of the independent variables. This action, which preserves the positiveness, extends Sophus Lie's group analysis of Newtonian dynamics.When applied to models of gas dynamics --~such as Euler system or Boltzmann equation,~-- in combination with Compensated Integrability, this yields new dispersive estimates. The most accurate one is obtained for mono-atomic gases. Then the space-time integral of $t\rho^\frac1d p$ is bounded in terms of the total mass and moment of inertia alone.


title: Cross-cultural Mood Perception in Pop Songs and its Alignment with Mood Detection Algorithms

abstract: Do people from different cultural backgrounds perceive the mood in music the same way? How closely do human ratings across different cultures approximate automatic mood detection algorithms that are often trained on corpora of predominantly Western popular music? Analyzing 166 participants responses from Brazil, South Korea, and the US, we examined the similarity between the ratings of nine categories of perceived moods in music and estimated their alignment with four popular mood detection algorithms. We created a dataset of 360 recent pop songs drawn from major music charts of the countries and constructed semantically identical mood descriptors across English, Korean, and Portuguese languages. Multiple participants from the three countries rated their familiarity, preference, and perceived moods for a given song. Ratings were highly similar within and across cultures for basic mood attributes such as sad, cheerful, and energetic. However, we found significant cross-cultural differences for more complex characteristics such as dreamy and love. To our surprise, the results of mood detection algorithms were uniformly correlated across human ratings from all three countries and did not show a detectable bias towards any particular culture. Our study thus suggests that the mood detection algorithms can be considered as an objective measure at least within the popular music context.


title: Self-Supervised Feature Learning of 1D Convolutional Neural Networks with Contrastive Loss for Eating Detection Using an In-Ear Microphone

abstract: The importance of automated and objective monitoring of dietary behavior is becoming increasingly accepted. The advancements in sensor technology along with recent achievements in machine-learning--based signal-processing algorithms have enabled the development of dietary monitoring solutions that yield highly accurate results. A common bottleneck for developing and training machine learning algorithms is obtaining labeled data for training supervised algorithms, and in particular ground truth annotations. Manual ground truth annotation is laborious, cumbersome, can sometimes introduce errors, and is sometimes impossible in free-living data collection. As a result, there is a need to decrease the labeled data required for training. Additionally, unlabeled data, gathered in-the-wild from existing wearables (such as Bluetooth earbuds) can be used to train and fine-tune eating-detection models. In this work, we focus on training a feature extractor for audio signals captured by an in-ear microphone for the task of eating detection in a self-supervised way. We base our approach on the SimCLR method for image classification, proposed by Chen et al. from the domain of computer vision. Results are promising as our self-supervised method achieves similar results to supervised training alternatives, and its overall effectiveness is comparable to current state-of-the-art methods. Code is available at https://github.com/mug-auth/ssl-chewing .


title: Defect reconstruction in a 2D semi-analytical waveguide model via derivative-based optimization

abstract: This paper considers the reconstruction of a defect in a two-dimensional waveguide during non-destructive ultrasonic inspection using a derivative-based optimization approach. The propagation of the mechanical waves is simulated by the Scaled Boundary Finite Element Method (SBFEM) that builds on a semi-analytical approach. The simulated data is then fitted to a given set of data describing the reflection of a defect to be reconstructed. For this purpose, we apply an iteratively regularized Gauss-Newton method in combination with algorithmic differentiation to provide the required derivative information accurately and efficiently. We present numerical results for three different kinds of defects, namely a crack, a delamination, and a corrosion. These examples show that the parameterization of the defect can be reconstructed efficiently and robustly in the presence of noise.


title: Bite-Weight Estimation Using Commercial Ear Buds

abstract: While automatic tracking and measuring of our physical activity is a well established domain, not only in research but also in commercial products and every-day life-style, automatic measurement of eating behavior is significantly more limited. Despite the abundance of methods and algorithms that are available in bibliography, commercial solutions are mostly limited to digital logging applications for smart-phones. One factor that limits the adoption of such solutions is that they usually require specialized hardware or sensors. Based on this, we evaluate the potential for estimating the weight of consumed food (per bite) based only on the audio signal that is captured by commercial ear buds (Samsung Galaxy Buds). Specifically, we examine a combination of features (both audio and non-audio features) and trainable estimators (linear regression, support vector regression, and neural-network based estimators) and evaluate on an in-house dataset of 8 participants and 4 food types. Results indicate good potential for this approach: our best results yield mean absolute error of less than 1 g for 3 out of 4 food types when training food-specific models, and 2.1 g when training on all food types together, both of which improve over an existing literature approach.


title: One-loop Feynman Integral Reduction by Differential Operators

abstract: For loop integrals, the standard method is reduction. A well-known reduction method for one-loop integrals is the Passarino-Veltman reduction. Inspired by the recent paper [1] where the tadpole reduction coefficients have been solved, in this paper we show the same technique can be used to give a complete integral reduction for any one-loop integrals. The differential operator method is an improved version of the PV-reduction method. Using this method, analytic expressions of all reduction coefficients of the master integrals can be given by algebraic recurrence relation easily. We demonstrate our method explicitly with several examples.


title: Spherically symmetric model atmospheres using approximate lambda operators V. Static inhomogeneous atmospheres of hot dwarf stars

abstract: Context. Clumping is a common property of stellar winds and is being incorporated to a solution of the radiative transfer equation coupled with kinetic equilibrium equations. However, in static hot model atmospheres, clumping and its influence on the temperature and density structures have not been considered and analysed at all to date. This is in spite of the fact that clumping can influence the interpretation of resulting spectra, as many inhomogeneities can appear there; for example, as a result of turbulent motions. Aims. We aim to investigate the effects of clumping on atmospheric structure for the special case of a static, spherically symmetric atmosphere assuming microclumping and a 1-D geometry. Methods. Static, spherically symmetric, non-LTE (local thermodynamic equilibrium) model atmospheres were calculated using the recent version of our code, which includes optically thin clumping. The matter is assumed to consist of dense clumps and a void interclump medium. Clumping is considered by means of clumping and volume filling factors, assuming all clumps are optically thin. Enhanced opacity and emissivity in clumps is multiplied by a volume filling factor to obtain their mean values. These mean values are used in the radiative transfer equation. Equations of kinetic equilibrium and the thermal balance equation use clump values of densities. Equations of hydrostatic and radiative equilibrium use mean values of densities. Results. The atmospheric structure was calculated for selected stellar parameters. Moderate differences were found in temperature structure. However, clumping causes enhanced continuum radiation for the Lyman-line spectral region, while radiation in other parts of the spectrum is lower, depending on the adopted model. The atomic level departure coefficients are influenced by clumping as well.


title: A Random Matrix Perspective on Random Tensors

abstract: Tensor models play an increasingly prominent role in many fields, notably in machine learning. In several applications, such as community detection, topic modeling and Gaussian mixture learning, one must estimate a low-rank signal from a noisy tensor. Hence, understanding the fundamental limits of estimators of that signal inevitably calls for the study of random tensors. Substantial progress has been recently achieved on this subject in the large-dimensional limit. Yet, some of the most significant among these results--in particular, a precise characterization of the abrupt phase transition (with respect to signal-to-noise ratio) that governs the performance of the maximum likelihood (ML) estimator of a symmetric rank-one model with Gaussian noise--were derived based of mean-field spin glass theory, which is not easily accessible to non-experts. In this work, we develop a sharply distinct and more elementary approach, relying on standard but powerful tools brought by years of advances in random matrix theory. The key idea is to study the spectra of random matrices arising from contractions of a given random tensor. We show how this gives access to spectral properties of the random tensor itself. For the aforementioned rank-one model, our technique yields a hitherto unknown fixed-point equation whose solution precisely matches the asymptotic performance of the ML estimator above the phase transition threshold in the third-order case. A numerical verification provides evidence that the same holds for orders 4 and 5, leading us to conjecture that, for any order, our fixed-point equation is equivalent to the known characterization of the ML estimation performance that had been obtained by relying on spin glasses. Moreover, our approach sheds light on certain properties of the ML problem landscape in large dimensions and can be extended to other models, such as asymmetric and non-Gaussian.


title: Self-supervised Answer Retrieval on Clinical Notes

abstract: Retrieving answer passages from long documents is a complex task requiring semantic understanding of both discourse and document context. We approach this challenge specifically in a clinical scenario, where doctors retrieve cohorts of patients based on diagnoses and other latent medical aspects. We introduce CAPR, a rule-based self-supervision objective for training Transformer language models for domain-specific passage matching. In addition, we contribute a novel retrieval dataset based on clinical notes to simulate this scenario on a large corpus of clinical notes. We apply our objective in four Transformer-based architectures: Contextual Document Vectors, Bi-, Poly- and Cross-encoders. From our extensive evaluation on MIMIC-III and three other healthcare datasets, we report that CAPR outperforms strong baselines in the retrieval of domain-specific passages and effectively generalizes across rule-based and human-labeled passages. This makes the model powerful especially in zero-shot scenarios where only limited training data is available.


title: The SMART protocol -- Pulse engineering of a global field for robust and universal quantum computation

abstract: Global control strategies for arrays of qubits are a promising pathway to scalable quantum computing. A continuous-wave global field provides decoupling of the qubits from background noise. However, this approach is limited by variability in the parameters of individual qubits in the array. Here we show that by modulating a global field simultaneously applied to the entire array, we are able to encode qubits that are less sensitive to the statistical scatter in qubit resonance frequency and microwave amplitude fluctuations, which are problems expected in a large scale system. We name this approach the SMART (Sinusoidally Modulated, Always Rotating and Tailored) qubit protocol. We show that there exist optimal modulation conditions for qubits in a global field that robustly provide improved coherence times. We discuss in further detail the example of spins in silicon quantum dots, in which universal one- and two-qubit control is achieved electrically by controlling the spin-orbit coupling of individual qubits and the exchange coupling between spins in neighbouring dots. This work provides a high-fidelity qubit operation scheme in a global field, significantly improving the prospects for scalability of spin-based quantum computer architectures.


title: CPSC: Conformal prediction with shrunken centroids for efficient prediction reliability quantification and data augmentation, a case in alternative herbal medicine classification with electronic nose

abstract: In machine learning applications, the reliability of predictions is significant for assisted decision and risk control. As an effective framework to quantify the prediction reliability, conformal prediction (CP) was developed with the CPKNN (CP with kNN). However, the conventional CPKNN suffers from high variance and bias and long computational time as the feature dimensionality increases. To address these limitations, a new CP framework-conformal prediction with shrunken centroids (CPSC) is proposed. It regularizes the class centroids to attenuate the irrelevant features and shrink the sample space for predictions and reliability quantification. To compare CPKNN and CPSC, we employed them in the classification of 12 categories of alternative herbal medicine with electronic nose as a case and assessed them in two tasks: 1) offline prediction: the training set was fixed and the accuracy on the testing set was evaluated; 2) online prediction with data augmentation: they filtered unlabeled data to augment the training data based on the prediction reliability and the final accuracy of testing set was compared. The result shows that CPSC significantly outperformed CPKNN in both two tasks: 1) CPSC reached a significantly higher accuracy with lower computation cost, and with the same credibility output, CPSC generally achieves a higher accuracy; 2) the data augmentation process with CPSC robustly manifested a statistically significant improvement in prediction accuracy with different reliability thresholds, and the augmented data were more balanced in classes. This novel CPSC provides higher prediction accuracy and better reliability quantification, which can be a reliable assistance in decision support.


title: Analysing digital in-memory computing for advanced finFET node

abstract: Digital In-memory computing improves energy efficiency and throughput of a data-intensive process, which incur memory thrashing and, resulting multiple same memory accesses in a von Neumann architecture. Digital in-memory computing involves accessing multiple SRAM cells simultaneously, which may result in a bit flip when not timed critically. Therefore we discuss the transient voltage characteristics of the bitlines during an SRAM compute. To improve the packaging density and also avoid MOSFET down-scaling issues, we use a 7-nm predictive PDK which uses a finFET node. The finFET process has discrete fins and a lower Voltage supply, which makes the design of in-memory compute SRAM difficult. In this paper, we design a 6T SRAM cell in 7-nm finFET node and compare its SNMs with a UMC 28nm node implementation. Further, we design and simulate the rest of the SRAM peripherals, and in-memory computation for an advanced finFET node.


title: The Homeomorphism Problem For Hyperbolic Manifolds I

abstract: We give a bounded runtime solution to the homeomorphism problem for closed hyperbolic 3-manifolds. This is an algorithm which, given two triangulations of hyperbolic 3-manifolds by at most $t$ tetrahedra, decides if they represent the same hyperbolic 3-manifold with runtime bounded by \[ 2^{2^{t^{O(t)}}}.\] We do this by first finding a hyperbolic structure on each manifold given as a geometric triangulation and then comparing the two as geometric manifolds.


title: Angle Based Feature Learning in GNN for 3D Object Detection using Point Cloud

abstract: In this paper, we present new feature encoding methods for Detection of 3D objects in point clouds. We used a graph neural network (GNN) for Detection of 3D objects namely cars, pedestrians, and cyclists. Feature encoding is one of the important steps in Detection of 3D objects. The dataset used is point cloud data which is irregular and unstructured and it needs to be encoded in such a way that ensures better feature encapsulation. Earlier works have used relative distance as one of the methods to encode the features. These methods are not resistant to rotation variance problems in Graph Neural Networks. We have included angular-based measures while performing feature encoding in graph neural networks. Along with that, we have performed a comparison between other methods like Absolute, Relative, Euclidean distances, and a combination of the Angle and Relative methods. The model is trained and evaluated on the subset of the KITTI object detection benchmark dataset under resource constraints. Our results demonstrate that a combination of angle measures and relative distance has performed better than other methods. In comparison to the baseline method(relative), it achieved better performance. We also performed time analysis of various feature encoding methods.


title: Generalization Bounds using Lower Tail Exponents in Stochastic Optimizers

abstract: Despite the ubiquitous use of stochastic optimization algorithms in machine learning, the precise impact of these algorithms and their dynamics on generalization performance in realistic non-convex settings is still poorly understood. While recent work has revealed connections between generalization and heavy-tailed behavior in stochastic optimization, this work mainly relied on continuous-time approximations; and a rigorous treatment for the original discrete-time iterations is yet to be performed. To bridge this gap, we present novel bounds linking generalization to the lower tail exponent of the transition kernel associated with the optimizer around a local minimum, in both discrete- and continuous-time settings. To achieve this, we first prove a data- and algorithm-dependent generalization bound in terms of the celebrated Fernique-Talagrand functional applied to the trajectory of the optimizer. Then, we specialize this result by exploiting the Markovian structure of stochastic optimizers, and derive bounds in terms of their (data-dependent) transition kernels. We support our theory with empirical results from a variety of neural networks, showing correlations between generalization error and lower tail exponents.


title: Multi-color stimulated Raman scattering with a frame-to-frame wavelength-tunable fiber-based light source

abstract: We present multi-color imaging by stimulated Raman scattering (SRS) enabled by an ultrafast fiber-based light source with integrated amplitude modulation and frame-to-frame wavelength tuning. With a relative intensity noise level of -153.7 dBc/Hz at 20.25 MHz the light source is well suited for SRS imaging and outperforms other fiber-based light source concepts for SRS imaging. The light source is tunable in under 5 ms per arbitrary wavelength step between 700 cm$^{-1}$ and 3200 cm$^{-1}$, which allows for addressing Raman resonances from the fingerprint to the CH-stretch region. Moreover, the compact and environmentally stable system is predestined for fast multi-color assessments of medical or rapidly evolving samples with high chemical specificity, paving the way for diagnostics and sensing outside of specialized laser laboratories.


title: CARLA: A Python Library to Benchmark Algorithmic Recourse and Counterfactual Explanation Algorithms

abstract: Counterfactual explanations provide means for prescriptive model explanations by suggesting actionable feature changes (e.g., increase income) that allow individuals to achieve favorable outcomes in the future (e.g., insurance approval). Choosing an appropriate method is a crucial aspect for meaningful counterfactual explanations. As documented in recent reviews, there exists a quickly growing literature with available methods. Yet, in the absence of widely available opensource implementations, the decision in favor of certain models is primarily based on what is readily available. Going forward - to guarantee meaningful comparisons across explanation methods - we present CARLA (Counterfactual And Recourse LibrAry), a python library for benchmarking counterfactual explanation methods across both different data sets and different machine learning models. In summary, our work provides the following contributions: (i) an extensive benchmark of 11 popular counterfactual explanation methods, (ii) a benchmarking framework for research on future counterfactual explanation methods, and (iii) a standardized set of integrated evaluation measures and data sets for transparent and extensive comparisons of these methods. We have open-sourced CARLA and our experimental results on Github, making them available as competitive baselines. We welcome contributions from other research groups and practitioners.


title: Towards Robust Object Detection: Bayesian RetinaNet for Homoscedastic Aleatoric Uncertainty Modeling

abstract: According to recent studies, commonly used computer vision datasets contain about 4% of label errors. For example, the COCO dataset is known for its high level of noise in data labels, which limits its use for training robust neural deep architectures in a real-world scenario. To model such a noise, in this paper we have proposed the homoscedastic aleatoric uncertainty estimation, and present a series of novel loss functions to address the problem of image object detection at scale. Specifically, the proposed functions are based on Bayesian inference and we have incorporated them into the common community-adopted object detection deep learning architecture RetinaNet. We have also shown that modeling of homoscedastic aleatoric uncertainty using our novel functions allows to increase the model interpretability and to improve the object detection performance being evaluated on the COCO dataset.


title: Bayesian Active Meta-Learning for Few Pilot Demodulation and Equalization

abstract: Two of the main principles underlying the life cycle of an artificial intelligence (AI) module in communication networks are adaptation and monitoring. Adaptation refers to the need to adjust the operation of an AI module depending on the current conditions; while monitoring requires measures of the reliability of an AI module's decisions. Classical frequentist learning methods for the design of AI modules fall short on both counts of adaptation and monitoring, catering to one-off training and providing overconfident decisions. This paper proposes a solution to address both challenges by integrating meta-learning with Bayesian learning. As a specific use case, the problems of demodulation and equalization over a fading channel based on the availability of few pilots are studied. Meta-learning processes pilot information from multiple frames in order to extract useful shared properties of effective demodulators across frames. The resulting trained demodulators are demonstrated, via experiments, to offer better calibrated soft decisions, at the computational cost of running an ensemble of networks at run time. The capacity to quantify uncertainty in the model parameter space is further leveraged by extending Bayesian meta-learning to an active setting. In it, the designer can select in a sequential fashion channel conditions under which to generate data for meta-learning from a channel simulator. Bayesian active meta-learning is seen in experiments to significantly reduce the number of frames required to obtain efficient adaptation procedure for new frames.


title: Aliasing Instabilities in the Numerical Evolution of the Einstein Field Equations

abstract: The Einstein field equations of gravitation are characterized by cross-scale, high-order nonlinear terms, representing a challenge for numerical modeling. In an exact spectral decomposition, high-order nonlinearities correspond to a convolution that numerically might lead to aliasing instabilities. We present a study of this problem, in vacuum conditions, based on the $3+1$ Baumgarte-Shibata-Shapiro-Nakamura (BSSN) formalism. We inspect the emergence of numerical artifacts, in a variety of conditions, by using the Spectral-FIltered Numerical Gravity codE (\texttt{SFINGE}) - a pseudo-spectral algorithm, based on a classical (Cartesian) Fourier decomposition. By monitoring the highest $k-$modes of the dynamical fields, we identify the culprits of the aliasing and propose procedures that cure such instabilities, based on the suppression of the aliased wavelengths. This simple algorithm, together with appropriate treatment of the boundary conditions, can be applied to a variety of gravitational problems, including those related to massive objects dynamics.


title: Convergence rate for the incompressible limit of nonlinear diffusion-advection equations

abstract: The incompressible limit of nonlinear diffusion equations of porous medium type has attracted a lot of attention in recent years, due to its ability to link the weak formulation of cell-population models to free boundary problems of Hele-Shaw type. Although vast literature is available on this singular limit, little is known on the convergence rate of the solutions. In this work, we compute the convergence rate in a negative Sobolev norm and, upon interpolating with BV-uniform bounds, we deduce a convergence rate in appropriate Lebesgue spaces.


title: Dynamic stop pooling for flexible and sustainable ride sharing

abstract: Ride sharing - the bundling of simultaneous trips of several people in one vehicle - may help to reduce the carbon footprint of human mobility. However, standard door-to-door ride sharing services trade reduced route length for increased user travel times and come with the burden of many stops and detours to pick up individual users. Requiring some users to walk to nearby shared stops reduces detours, but could become inefficient if spatio-temporal demand patterns do not well fit the stop locations. Here, we present a simple computational model of dynamic stop pooling with flexible stop positions, and analyze its influence on the performance of ride sharing services. We find that dynamic stop pooling does a-priori not save route length, but occupancy. Intriguingly, it also reduces the travel time although users walk parts of their trip. Together, dynamic stop pooling may break the trade-off between route lengths and travel time in door-to-door ride sharing, thus enabling higher sustainability and service quality.


title: Identifying historical roots in paediatric echocardiography using RPYS

abstract: Echocardiography is a non-invasive diagnostic tool which can be performed on children of all ages including foetuses and new-borns, due to the fact that no radiation is involved. No historical bibliometric study or influential papers in paediatric echocardiography has been performed yet, so our study was aimed to close this gap. References Publication Years Spectrography (RPYS), more precisely CitedReferenceExplorer (CRE) software tool was employed to achieve this aim. We identified 35 influential papers and validation showed that both, the RPYS method and CRE tool performed satisfactorily.


title: Classification of real trivectors in dimension nine

abstract: We classify real trivectors in dimension 9. The corresponding classification over the field C of complex numbers was obtained by Vinberg and Elashvili in 1978. One of the main tools used for their classification was the construction of the representation of SL(9,C) on the space of complex trivectors of C^9 as a theta-representation corresponding to a Z/3Z-grading of the simple complex Lie algebra of type E_8. This divides the trivectors into three groups: nilpotent, semisimple, and mixed trivectors. Our classification follows the same pattern. We use Galois cohomology, first and second, to obtain the classification over R.


title: Dissipative localized states and breathers in phase mismatched singly resonant optical parametric oscillators: Bifurcation structure and stability

abstract: We study the emergence of dissipative localized states in phase mismatched singly resonant optical parametric oscillators. These states arise in two different bistable configurations due to the locking of fronts waves connecting the two coexisting states. In one of these configurations the bistabiity is mediated by the coexistence of two uniform states. Here the localized states are organized in a collapsed snaking bifurcation structure. Moreover, these states undergo oscillatory instabilities which lead to a breathing behavior. When the the bistability is related with the coexistence of an uniform state and a spatially periodic pattern, localized states are organized in a bifurcation structure similar to the standard homoclinic snaking. Performing an exhaustive bifurcation analysis, we characterize in detail the previous structures, their linear stability and the modification of their dynamics as a function of the control parameters of the system.


title: On the importance of wave planet interactions for the migration of two super-Earths embedded in a protoplanetary disk

abstract: We investigate a repulsion mechanism between two low-mass planets migrating in a protoplanetary disk, for which the relative migration switches from convergent to divergent. This mechanism invokes density waves emitted by one planet transferring angular momentum to the coorbital region of the other and then directly to it through the horseshoe drag. We formulate simple analytical estimates, which indicate when the repulsion mechanism is effective. One condition for a planet to be repelled is that it forms a partial gap in the disk and another is that this should contain enough material to support angular momentum exchange with it. Using two-dimensional hydrodynamical simulations we obtain divergent migration of two super-Earths embedded in a protoplanetary disk because of repulsion between them and verify these conditions. To investigate the importance of resonant interaction we study the migration of planet pairs near first-order commensurabilities. It appears that proximity to resonance is significant but not essential. In this context we find repulsion still occurs when the gravitational interaction between the planets is removed sugesting the importance of angular momentum transfer through waves excited by another planet. This may occur through the scattering of coorbital material (the horseshoe drag), or material orbiting close by. Our results indicate that if conditions favor the repulsion between two planets described above, we expect to observe planet pairs with their period ratios greater, often only slightly greater, than resonant values or possibly rarity of commensurability.


title: Ground state energy of the polarized diluted gas of interacting spin $1/2$ fermions

abstract: The effective field theory approach simplifies the perturbative computation of the ground state energy of the diluted gas of fermions allowing in the case of the unpolarized system to easily re-derive the classic results up to the $(k_{\rm F}a_0)^2$ order (where $k_{\rm F}$ is the system's Fermi momentum and $a_0$ the $s$-wave scattering length) and (with more labour) to extend it up to the order $(k_{\rm F}a_0)^4$. The corresponding expansion of the ground state energy of the polarized gas of spin $1/2$ fermions is known analytically (to our best knowledge) only up to the $k_{\rm F}a_0$ (where $k_{\rm F}$ stands for $k_{{\rm F}\uparrow}$ or $k_{{\rm F}\downarrow}$) order. Here we show that the same effective field theory method allows to easily compute also the order $(k_{\rm F}a_0)^2$ correction to this result.


title: Improved Efficiency of Multilevel Monte Carlo for Stochastic PDE through Strong Pairwise Coupling

abstract: Multilevel Monte Carlo (MLMC) has become an important methodology in applied mathematics for reducing the computational cost of weak approximations. For many problems, it is well-known that strong pairwise coupling of numerical solutions in the multilevel hierarchy is needed to obtain efficiency gains. In this work, we show that strong pairwise coupling indeed is also important when (MLMC) is applied to stochastic partial differential equations (SPDE) of reaction-diffusion type, as it can improve the rate of convergence and thus improve tractability. For the (MLMC) method with strong pairwise coupling that was developed and studied numerically on filtering problems in [{\it Chernov et al., Numer. Math., 147 (2021), 71-125}], we prove that the rate of computational efficiency is higher than for existing methods. We also provide numerical comparisons with alternative coupling ideas on linear and nonlinear SPDE to illustrate the importance of this feature.


title: On deformed preprojective algebras

abstract: Deformed preprojective algebras are generalizations of the usual preprojective algebras introduced by Crawley-Boevey and Holland, which have applications to Kleinian singularities, the Deligne-Simpson problem, integrable systems and noncommutative geometry. In this paper we offer three contributions to the study of such algebras: (1) the 2-Calabi-Yau property; (2) the unification of the reflection functors of Crawley-Boevey and Holland with reflection functors for the usual preprojective algebras; and (3) the classification of tilting ideals in 2-Calabi-Yau algebras, and especially in deformed preprojective algebras for extended Dynkin quivers.


title: ricu: R's Interface to Intensive Care Data

abstract: Providing computational infrastructure for handling diverse intensive care unit (ICU) datasets, the R package 'ricu' enables writing dataset-agnostic analysis code, thereby facilitating multi-center training and validation of machine learning models. The package is designed with an emphasis on extensibility both to new datasets as well as clinical data concepts, and currently supports the loading of around 100 patient variables corresponding to a total of 319,402 ICU admissions from 4 data sources collected in Europe and the United States. By allowing for the addition of user-specified medical concepts and data sources the aim of 'ricu' is to foster robust, data-based intensive care research, allowing the user to externally validate their method or conclusion with relative ease, and in turn facilitating reproducible and therefore transparent work in this field.


title: Missing Rung Problem in Vibrational Ladder Climbing

abstract: We observed vanishing of the transition dipole moment, interrupting vibrational ladder climbing (VLC) in molecular systems. We clarified the mechanism of this phenomenon and present a method to use an additional chirped pulse to preserve the VLC. To show the effectiveness of our method, we conducted wavepacket dynamics simulations for LiH dissociations with chirped pulses. The results indicate that the efficiency of LiH dissociation is significantly improved by our method compared to conventional methods. We also revealed the quantum interference effect behind the excitation process of VLC.


title: Quantum Computation Protocol for Dressed Spins in a Global Field

abstract: Spin qubits are contenders for scalable quantum computation because of their long coherence times demonstrated in a variety of materials, but individual control by frequency-selective addressing using pulsed spin resonance creates severe technical challenges for scaling up to many qubits. This individual resonance control strategy requires each spin to have a distinguishable frequency, imposing a maximum number of spins that can be individually driven before qubit crosstalk becomes unavoidable. Here we describe a complete strategy for controlling a large array of spins in quantum dots dressed by an on-resonance global field, namely a field that is constantly driving the spin qubits, to dynamically decouple from the effects of background magnetic field fluctuations. This approach -- previously implemented for the control of single electron spins bound to electrons in impurities -- is here harmonized with all other operations necessary for universal quantum computing with spins in quantum dots. We define the logical states as the dressed qubit states and discuss initialization and readout utilizing Pauli spin blockade, as well as single- and two-qubit control in the new basis. Finally, we critically analyze the limitations imposed by qubit variability and potential strategies to improve performance.


title: Mean Field Game of Optimal Relative Investment with Jump Risk

abstract: This paper studies the n-player game and the mean field game under the CRRA relative performance on terminal wealth, in which the interaction occurs by peer competition. In the model with n agents, the price dynamics of underlying risky assets depend on a common noise and contagious jump risk modelled by a multi-dimensional nonlinear Hawkes process. With a continuum of agents, we formulate the MFG problem and characterize a deterministic mean field equilibrium in an analytical form under some conditions, allowing us to investigate some impacts of model parameters in the limiting model and discuss some financial implications. Moreover, based on the mean field equilibrium, we construct an approximate Nash equilibrium for the n-player game when n is sufficiently large. The explicit order of the approximation error is also derived.


title: Training face verification models from generated face identity data

abstract: Machine learning tools are becoming increasingly powerful and widely used. Unfortunately membership attacks, which seek to uncover information from data sets used in machine learning, have the potential to limit data sharing. In this paper we consider an approach to increase the privacy protection of data sets, as applied to face recognition. Using an auxiliary face recognition model, we build on the StyleGAN generative adversarial network and feed it with latent codes combining two distinct sub-codes, one encoding visual identity factors, and, the other, non-identity factors. By independently varying these vectors during image generation, we create a synthetic data set of fictitious face identities. We use this data set to train a face recognition model. The model performance degrades in comparison to the state-of-the-art of face verification. When tested with a simple membership attack our model provides good privacy protection, however the model performance degrades in comparison to the state-of-the-art of face verification. We find that the addition of a small amount of private data greatly improves the performance of our model, which highlights the limitations of using synthetic data to train machine learning models.


title: LICHEE: Improving Language Model Pre-training with Multi-grained Tokenization

abstract: Language model pre-training based on large corpora has achieved tremendous success in terms of constructing enriched contextual representations and has led to significant performance gains on a diverse range of Natural Language Understanding (NLU) tasks. Despite the success, most current pre-trained language models, such as BERT, are trained based on single-grained tokenization, usually with fine-grained characters or sub-words, making it hard for them to learn the precise meaning of coarse-grained words and phrases. In this paper, we propose a simple yet effective pre-training method named LICHEE to efficiently incorporate multi-grained information of input text. Our method can be applied to various pre-trained language models and improve their representation capability. Extensive experiments conducted on CLUE and SuperGLUE demonstrate that our method achieves comprehensive improvements on a wide variety of NLU tasks in both Chinese and English with little extra inference cost incurred, and that our best ensemble model achieves the state-of-the-art performance on CLUE benchmark competition.


title: Coalitional Control for Self-Organizing Agents

abstract: Coalitional control is concerned with the management of multi-agent systems where cooperation cannot be taken for granted (due to, e.g., market competition, logistics). This paper proposes a model predictive control (MPC) framework aimed at large-scale dynamically-coupled systems whose individual components, possessing a limited model of the system, are controlled independently, pursuing possibly competing objectives. The emergence of cooperating clusters of controllers is contemplated through an autonomous negotiation protocol, based on the characterization as a coalitional game of the benefit derived by a broader feedback and the alignment of the individual objectives. Specific mechanisms for the cooperative benefit redistribution that relax the cognitive requirements of the game are employed to compensate for possible local cost increases due to cooperation. As a result, the structure of the overall MPC feedback can be adapted online to the degree of interaction between different parts of the system, while satisfying the individual interests of the agents. A wide-area control application for the power grid with the objective of minimizing frequency deviations and undesired inter-area power transfers is used as study case.


title: Learn to Match: Automatic Matching Network Design for Visual Tracking

abstract: Siamese tracking has achieved groundbreaking performance in recent years, where the essence is the efficient matching operator cross-correlation and its variants. Besides the remarkable success, it is important to note that the heuristic matching network design relies heavily on expert experience. Moreover, we experimentally find that one sole matching operator is difficult to guarantee stable tracking in all challenging environments. Thus, in this work, we introduce six novel matching operators from the perspective of feature fusion instead of explicit similarity learning, namely Concatenation, Pointwise-Addition, Pairwise-Relation, FiLM, Simple-Transformer and Transductive-Guidance, to explore more feasibility on matching operator selection. The analyses reveal these operators' selective adaptability on different environment degradation types, which inspires us to combine them to explore complementary features. To this end, we propose binary channel manipulation (BCM) to search for the optimal combination of these operators. BCM determines to retrain or discard one operator by learning its contribution to other tracking steps. By inserting the learned matching networks to a strong baseline tracker Ocean, our model achieves favorable gains by $67.2 \rightarrow 71.4$, $52.6 \rightarrow 58.3$, $70.3 \rightarrow 76.0$ success on OTB100, LaSOT, and TrackingNet, respectively. Notably, Our tracker, dubbed AutoMatch, uses less than half of training data/time than the baseline tracker, and runs at 50 FPS using PyTorch. Code and model will be released at https://github.com/JudasDie/SOTS.


title: Relation Aware Semi-autoregressive Semantic Parsing for NL2SQL

abstract: Natural language to SQL (NL2SQL) aims to parse a natural language with a given database into a SQL query, which widely appears in practical Internet applications. Jointly encode database schema and question utterance is a difficult but important task in NL2SQL. One solution is to treat the input as a heterogeneous graph. However, it failed to learn good word representation in question utterance. Learning better word representation is important for constructing a well-designed NL2SQL system. To solve the challenging task, we present a Relation aware Semi-autogressive Semantic Parsing (\MODN) ~framework, which is more adaptable for NL2SQL. It first learns relation embedding over the schema entities and question words with predefined schema relations with ELECTRA and relation aware transformer layer as backbone. Then we decode the query SQL with a semi-autoregressive parser and predefined SQL syntax. From empirical results and case study, our model shows its effectiveness in learning better word representation in NL2SQL.


title: Conservation with moving meshes over orography

abstract: Adaptive meshes have the potential to improve the accuracy and efficiency of atmospheric modelling by increasing resolution where it is most needed. Mesh re-distribution, or r-adaptivity, adapts by moving the mesh without changing the connectivity. This avoids some of the challenges with h-adaptivity (adding and removing points): the solution does not need to be mapped between meshes, which can be expensive and introduces errors, and there are no load balancing problems on parallel computers. A long standing problem with both forms of adaptivity has been changes in volume of the domain as resolution changes at an uneven boundary. We propose a solution to exact local conservation and maintenance of uniform fields while the mesh changes volume as it moves over orography. This is solved by introducing a volume adjustment parameter which tracks the true cell volumes without using expensive conservative mapping. A finite volume solution of the advection equation over orography on moving meshes is described and results are presented demonstrating improved accuracy for cost using moving meshes. Exact local conservation and maintenance of uniform fields is demonstrated and the corrected mesh volume is preserved. We use optimal transport to generate meshes which are guaranteed not to tangle and are equidistributed with respect to a monitor function. This leads to a Monge-Amp\`{e}re equation which is solved with a Newton solver. The superiority of the Newton solver over other techniques is demonstrated in the appendix. However the Newton solver is only efficient if it is applied to the left hand side of the Monge-Amp\`{e}re equation with fixed point iterations for the right hand side.


title: Super-Resolution of Near-Surface Temperature Utilizing Physical Quantities for Real-Time Prediction of Urban Micrometeorology

abstract: The present paper proposes a super-resolution (SR) model based on a convolutional neural network and applies it to the near-surface temperature in urban areas. The SR model incorporates a skip connection, a channel attention mechanism, and separated feature extractors for the inputs of temperature, building height, downward shortwave radiation, and horizontal velocity. We train the SR model with sets of low-resolution (LR) and high-resolution (HR) images from building-resolving large-eddy simulations (LESs) in a city, where the horizontal resolutions of LR and HR are 20 and 5 m, respectively. The generalization capability of the SR model is confirmed with LESs in another city. The estimated HR temperature fields are more accurate than those of the bicubic interpolation and image SR model that takes only the temperature as its input. Except for the temperature input, the building height is the most important to reconstruct the HR temperature and enables the SR model to reduce errors in temperature near building boundaries. The SR model considers the appropriate boundary for each building from its height information. The analysis of attention weights indicates that the importance of the building height increases as the downward shortwave radiation becomes larger. The contrast between sun and shade is strengthened with the increase in solar radiation, which may affect the temperature distribution. The short inference time suggests the potential of the proposed SR model to facilitate a real-time HR prediction in metropolitan areas by combining it with an LR building-resolving LES model.


title: Application of Blockchain in Healthcare and Health Insurance Sector

abstract: Technology has evolved over the years, making our lives easier. It has impacted the healthcare sector, increasing the average life expectancy of human beings. Still, there are gaps that remain unaddressed. There is a lack of transparency in the healthcare system, which results in inherent trust problems between patients and hospitals. In the present day, a patient does not know whether he or she will get the proper treatment from the hospital for the fee charged. A patient can claim reimbursement of the medical bill from any insurance company. However, today there is minimal scope for the Insurance Company to verify the validity of such bills or medical records. A patient can provide fake details to get financial benefits from the insurance company. Again, there are trust issues between the patient (i.e., the insurance claimer) and the insurance company. Blockchain integrated with the smart contract is a well-known disruptive technology that builds trust by providing transparency to the system. In this paper, we propose a blockchain-enabled Secure and Smart HealthCare System. Fairness of all the entities: patient, hospital, or insurance company involved in the system is guaranteed with no one trusting each other. Privacy and security of patients' medical data are ensured as well. We also propose a method for privacy-preserving sharing of aggregated data with the research community for their own purpose. Shared data must not be personally identifiable, i.e, no one can link the acquired data to the identity of any patient or their medical history. We have implemented the prototype in the Ethereum platform and Ropsten test network, and have included the analysis as well.


title: Energy Efficiency Aspects of the AMD Zen 2 Architecture

abstract: In High Performance Computing, systems are evaluated based on their computational throughput. However, performance in contemporary server processors is primarily limited by power and thermal constraints. Ensuring operation within a given power envelope requires a wide range of sophisticated control mechanisms. While some of these are handled transparently by hardware control loops, others are controlled by the operating system. A lack of publicly disclosed implementation details further complicates this topic. However, understanding these mechanisms is a prerequisite for any effort to exploit the full computing capability and to minimize the energy consumption of today's server systems. This paper highlights the various energy efficiency aspects of the AMD Zen 2 microarchitecture to facilitate system understanding and optimization. Key findings include qualitative and quantitative descriptions regarding core frequency transition delays, workload-based frequency limitations, effects of I/O die P-states on memory performance as well as discussion on the built-in power monitoring capabilities and its limitations. Moreover, we present specifics and caveats of idle states, wakeup times as well as the impact of idling and inactive hardware threads and cores on the performance of active resources such as other cores.


title: Cross-Modal Knowledge Transfer via Inter-Modal Translation and Alignment for Affect Recognition

abstract: Multi-modal affect recognition models leverage complementary information in different modalities to outperform their uni-modal counterparts. However, due to the unavailability of modality-specific sensors or data, multi-modal models may not be always employable. For this reason, we aim to improve the performance of uni-modal affect recognition models by transferring knowledge from a better-performing (or stronger) modality to a weaker modality during training. Our proposed multi-modal training framework for cross-modal knowledge transfer relies on two main steps. First, an encoder-classifier model creates task-specific representations for the stronger modality. Then, cross-modal translation generates multi-modal intermediate representations, which are also aligned in the latent space with the stronger modality representations. To exploit the contextual information in temporal sequential affect data, we use Bi-GRU and transformer encoder. We validate our approach on two multi-modal affect datasets, namely CMU-MOSI for binary sentiment classification and RECOLA for dimensional emotion regression. The results show that the proposed approach consistently improves the uni-modal test-time performance of the weaker modalities.


title: Koshliakov zeta functions I: Modular Relations

abstract: We examine an unstudied manuscript of N.~S.~Koshliakov over $150$ pages long and containing the theory of two interesting generalizations $\zeta_p(s)$ and $\eta_p(s)$ of the Riemann zeta function $\zeta(s)$, which we call \emph{Koshliakov zeta functions}. His theory has its genesis in a problem in the analytical theory of heat distribution which was analyzed by him. In this paper, we further build upon his theory and obtain two new modular relations in the setting of Koshliakov zeta functions, each of which gives an infinite family of identities, one for each $p\in\mathbb{R^{+}}$. The first one is a generalization of Ramanujan's famous formula for $\zeta(2m+1)$ and the second is an elegant extension of a modular relation on page $220$ of Ramanujan's Lost Notebook. Several interesting corollaries and applications of these modular relations are obtained including a new representation for $\zeta(4m+3)$.


title: Special values of $L$-functions on regular arithmetic schemes of dimension $1$

abstract: We construct a well-behaved Weil-\'etale complex for a large class of $\mathbb{Z}$-constructible sheaves on a regular irreducible scheme $U$ of finite type over $\mathbb{Z}$ and of dimension $1$. We then give a formula for the special value at $s=0$ of the $L$-function associated to any $\mathbb{Z}$-constructible sheaf on $U$ in terms of Euler characteristics of Weil-\'etale cohomology; for smooth proper curves, we obtain the formula of arXiv:2009.14504. We deduce a special value formula for Artin $L$-functions twisted by a singular irreducible scheme $X$ of finite type over $\mathbb{Z}$ and of dimension $1$. This generalizes and improves all results in arXiv:1611.01720; as a special case, we obtain a special value formula for the arithmetic zeta function of $X$.


title: Isometries of ultrametric normed spaces

abstract: We show that the group of isometries of an ultrametric normed space can be seen as a kind of a fractal. Then, we apply this description to study ultrametric counterparts of some classical problems in Archimedean analysis, such as the so called Probl\`eme des rotations de Mazur or Tingley's problem. In particular, it turns out that, in contrast with the case of real normed spaces, isometries between ultrametric normed spaces can be very far from being linear.


title: The CirCor DigiScope Dataset: From Murmur Detection to Murmur Classification

abstract: Cardiac auscultation is one of the most cost-effective techniques used to detect and identify many heart conditions. Computer-assisted decision systems based on auscultation can support physicians in their decisions. Unfortunately, the application of such systems in clinical trials is still minimal since most of them only aim to detect the presence of extra or abnormal waves in the phonocardiogram signal, i.e., only a binary ground truth variable (normal vs abnormal) is provided. This is mainly due to the lack of large publicly available datasets, where a more detailed description of such abnormal waves (e.g., cardiac murmurs) exists. To pave the way to more effective research on healthcare recommendation systems based on auscultation, our team has prepared the currently largest pediatric heart sound dataset. A total of 5282 recordings have been collected from the four main auscultation locations of 1568 patients, in the process, 215780 heart sounds have been manually annotated. Furthermore, and for the first time, each cardiac murmur has been manually annotated by an expert annotator according to its timing, shape, pitch, grading, and quality. In addition, the auscultation locations where the murmur is present were identified as well as the auscultation location where the murmur is detected more intensively. Such detailed description for a relatively large number of heart sounds may pave the way for new machine learning algorithms with a real-world application for the detection and analysis of murmur waves for diagnostic purposes.


title: What do Firms Gain from Patenting? The Case of the Global ICT Industry

abstract: This study investigates the causal relationship between patent grants and firms' dynamics in the Information and Communication Technology (ICT) industry, as the latter is a peculiar sector of modern economies, often under the lens of antitrust authorities. For our purpose, we exploit matched information about financial accounts and patenting activity in 2009-2017 by 179,660 companies operating in 39 countries. Preliminarily, we show how bigger companies are less than 2% of the sample, although they concentrate about 89% of the grants obtained in the period of analyses. Thus, we test that patent grants in the ICT industry have a significant and large impact on market shares and firm size of smaller companies (31.5% and 30.7%, respectively) in the first year after the grants, while we have no evidence of an impact for bigger companies. After a novel instrumental variable strategy that exploits information at the level of patent offices, we confirm that most of the effects on smaller companies are due to the protection of property rights and not to the innovative content of inventions. Finally, we never observe a significant impact on either profitability or productivity for any firm size category. Eventually, we discuss how our findings support the idea that the ICT industry is a case of endogenous R&D sunk costs, which prevent profit margins from rising in the presence of a relatively high market concentration.


title: Estimating the Peer Degree of Reachable Peers in the Bitcoin P2P Network

abstract: A recent spam wave of IP addresses in the Bitcoin P2P network allowed us to estimate the degree distribution of reachable peers in the network. The resulting distribution shows that about every second reachable peer runs with Bitcoin Core's default setting of a maximum of 125 concurrent connections and nearly all connection slots are taken. We validate this result and, in addition, use our observations of the spam wave to group addresses that belong to the same peer. By doing this grouping, we improve on previous measurements and show that simply counting addresses overestimates the number of reachable peers by 13 %.


title: Progress and opportunities in modelling environmentally assisted cracking

abstract: Environmentally assisted cracking phenomena are widespread across the transport, defence, energy and construction sectors. However, predicting environmentally assisted fractures is a highly cross-disciplinary endeavour that requires resolving the multiple material-environment interactions taking place. In this manuscript, an overview is given of recent breakthroughs in the modelling of environmentally assisted cracking. The focus is on the opportunities created by two recent developments: phase field and multi-physics modelling. The possibilities enabled by the confluence of phase field methods and electro-chemo-mechanics modelling are discussed in the context of three environmental assisted cracking phenomena of particular engineering interest: hydrogen embrittlement, localised corrosion and corrosion fatigue. Mechanical processes such as deformation and fracture can be coupled with chemical phenomena like local reactions, ionic transport and hydrogen uptake and diffusion. Moreover, these can be combined with the prediction of an evolving interface, such as a growing pit or a crack, as dictated by a phase field variable that evolves based on thermodynamics and local kinetics. Suitable for both microstructural and continuum length scales, this new generation of simulation-based, multi-physics phase field models can open new modelling horizons and enable Virtual Testing in harmful environments.


title: Small-Scale air turbulence structure, microphysical time scales and local supersaturation balance at a warm Cloud Top Boundary

abstract: Recent results have shown that there is an acceleration in the spread of the size distribution of droplet populations in the region bordering the cloud and undersaturated ambient. We have analyzed the supersaturation balance in this region, which is typically a highly intermittent shearless turbulent mixing layer, under a condition where there is no mean updraft. We have investigated the evolution of the cloud-clear air interface and of the droplets therein via direct numerical simulations. We have compared horizontal averages of the phase relaxation, evaporation, reaction and condensation times within the cloud-clear air interface for the size distributions of the initial monodisperse and polydisperse droplets. For the monodisperse population, a clustering of the values of the reaction, phase and evaporation times, that is around 20-30 seconds, is observed in the central area of the mixing layer, just before the location where the maximum value of the supersaturation turbulent flux occurs. This clustering of values is similar for the polydisperse population but also includes the condensation time. The mismatch between the time derivative of the supersaturation and the condensation term in the interfacial mixing layer is correlated with the planar covariance of the horizontal longitudinal velocity derivatives of the carrier air flow and the supersaturation field, thus suggesting that a quasi-linear relationship may exist between these quantities.


title: A comparison of data-driven approaches to build low-dimensional ocean models

abstract: We present a comprehensive inter-comparison of linear regression (LR), stochastic, and deep-learning approaches for reduced-order statistical emulation of ocean circulation. The reference dataset is provided by an idealized, eddy-resolving, double-gyre ocean circulation model. Our goal is to conduct a systematic and comprehensive assessment and comparison of skill, cost, and complexity of statistical models from the three methodological classes. The model based on LR is considered as a baseline. Additionally, we investigate its additive white noise augmentation and a multi-level stochastic approach, deep-learning methods, hybrid frameworks (LR plus deep-learning), and simple stochastic extensions of deep-learning and hybrid methods. The assessment metrics considered are: root mean squared error, anomaly cross-correlation, climatology, variance, frequency map, forecast horizon, and computational cost. We found that the multi-level linear stochastic approach performs the best for both short- and long-timescale forecasts. The deep-learning hybrid models augmented by additive state-dependent white noise came second, while their deterministic counterparts failed to reproduce the characteristic frequencies in climate-range forecasts. Pure deep learning implementations performed worse than LR and its noise augmentations. Skills of LR and its white noise extension were similar on short timescales, but the latter performed better on long timescales, while LR-only outputs decay to zero for long simulations. Overall, our analysis promotes multi-level LR stochastic models with memory effects, and hybrid models with linear dynamical core augmented by additive stochastic terms learned via deep learning, as a more practical, accurate, and cost-effective option for ocean emulation than pure deep-learning solutions.


title: Active Learning in Gaussian Process State Space Model

abstract: We investigate active learning in Gaussian Process state-space models (GPSSM). Our problem is to actively steer the system through latent states by determining its inputs such that the underlying dynamics can be optimally learned by a GPSSM. In order that the most informative inputs are selected, we employ mutual information as our active learning criterion. In particular, we present two approaches for the approximation of mutual information for the GPSSM given latent states. The proposed approaches are evaluated in several physical systems where we actively learn the underlying non-linear dynamics represented by the state-space model.


title: The Role of Nanotechnology to Combat Major Recent Worldwide Challenges

abstract: There are two major challenging problems facing humanity today which are the ongoing Coronavirus disease (COVID-19) pandemic started in late 2019 and the environmental crisis due to fossil fuels combustion causing over increase of greenhouse gases in the atmosphere. While the first problem seems to be a temporary remedy for the second one due to less travels by individuals, actual solution of both problems may require a comprehensive scientific, technological, and socioeconomic decisions. In this report we look at the role of nanotechnology and nano-scale materials towards the solution of these two major challenges.


title: The reaction to news in live betting

abstract: Sports betting markets have grown very rapidly recently, with the total European gambling market worth 98.6 billion euro in 2019. Considering a high-resolution (1 Hz) data set provided by a large European bookmaker, we investigate the effect of news on the dynamics of live betting. In particular, we consider stakes placed in a live betting market during football matches. Accounting for the general market activity level within a state-space modelling framework, we focus on the market's response to events such as goals (i.e. major news), but also to the general situation within a match such as the uncertainty about the outcome. Our results indicate that markets might overreact to recent news, confirming cognitive biases known from psychology and behavioural economics.


title: Extremal product-one free sequences over $C_n \rtimes_s C_2$

abstract: Let $G$ be a finite group multiplicatively written. The small Davenport constant of $G$ is the maximum positive integer ${\sf d}(G)$ such that there exists a sequence $S$ of length ${\sf d}(G)$ for which every subsequence of $S$ is product-one free. Let $s^2 \equiv 1 \pmod n$, where $s \not\equiv \pm1 \pmod n$. It has been proven that ${\sf d}(C_n \rtimes_s C_2) = n$ (see Lemma 6 of [Zhuang, Gao; Europ. J. Combin. 26 (2005), 1053-1059]). In this paper, we determine all sequences over $C_n \rtimes_s C_2$ of length $n$ which are product-one free. It completes the classification of all product-one free sequences over every group of the form $C_n \rtimes_s C_2$, including the quasidihedral groups and the modular maximal-cyclic groups.


title: The main zero-sum constants over $D_{2n} \times C_2$

abstract: Let $C_2$ be the cyclic group of order $2$ and $D_{2n}$ be the dihedral group of order $2n$, where $n$ is even. In this paper, we provide the exact values of some zero-sum constants over $D_{2n} \times C_2$, namely small Davenport constant, Gao constant, $\eta$-constant and Erd\H os-Ginzburg-Ziv constant. As a consequence, we prove the Gao's and Zhuang-Gao's Conjectures for this group. These are the first concrete results on zero-sum problems for a family of non-abelian groups of rank greater than $2$.


title: Ising Game on Graphs

abstract: Static and dynamic equilibria in noisy binary choice (Ising) games on complete and random graphs in the annealed approximation are analysed. Two versions, an Ising game with interaction term defined in accordance with the Ising model in statistical physics and a reduced Ising game with a customary definition of interaction term in game theory on graphs, are considered. A detailed analysis of hysteresis phenomenon shaping the pattern of static equilibria based on consideration of elasticity with respect to external influence is conducted. Fokker-Planck equations describing dynamic versions of the games under consideration are written and their asymptotic stationary solutions derived. It is shown that domains of parameters corresponding to the maxima of these probability distributions are identical with the corresponding hysteresis ranges for static equilibria. Same result holds for domains defining local stability of solutions of the evolution equations for the moments. In all the cases considered it is shown that the results for the reduced Ising game coincide with those obtained for the Ising game on complete graphs. It is shown that for s special case of logistic noise the results obtained for static equilibria for the Ising game reproduce those in the Ising model on graphs in statistical physics.


title: A physiology-inspired framework for holistic city simulations

abstract: Life, services and activities within cities have commonly been studied by separate disciplines, each one independent from the others. One such approach is the computer simulation, which enables in-depth modelling and cost-effective evaluation of city phenomena. However, the adoption of integrated city simulations faces several barriers, such as managerial, social, and technical, despite its potential to support city planning and policymaking. This paper introduces the City Physiology: a new conceptual framework to facilitate the integration of city layers when designing holistic simulators. The physiology is introduced and applied through a process of three steps. Firstly, a literature review is offered in order to study the terminology and the progress already made towards integrated modelling of different urban systems. Secondly, interactions between urban systems are extracted from the approaches studied before. Finally, the pipeline to carry out the integration strategy is described. In addition to providing a conceptual tool for holistic simulations, the framework enables the discovery of new research lines generated by previously unseen connections between city layers. Being an open framework, available to all researchers to use and broaden, the authors of this paper envisage that it will be a valuable resource in establishing an exact science of cities.


title: Euler sums of generalized alternating hyperharmonic numbers II

abstract: In this paper, we introduce a new type of generalized alternating hyperharmonic numbers $H_n^{(p,r,s_{1},s_{2})}$, and show that Euler sums of the generalized alternating hyperharmonic numbers $H_n^{(p,r,s_{1},s_{2})}$ can be expressed in terms of linear combinations of classical (alternating) Euler sums.


title: Parameter uncertainty quantification in an idealized GCM with a seasonal cycle

abstract: Climate models are generally calibrated manually by comparing selected climate statistics, such as the global top-of-atmosphere energy balance, to observations. The manual tuning only targets a limited subset of observational data and parameters. Bayesian calibration can estimate climate model parameters and their uncertainty using a larger fraction of the available data and automatically exploring the parameter space more broadly. In Bayesian learning, it is natural to exploit the seasonal cycle, which has large amplitude, compared with anthropogenic climate change, in many climate statistics. In this study, we develop methods for the calibration and uncertainty quantification (UQ) of model parameters exploiting the seasonal cycle, and we demonstrate a proof-of-concept with an idealized general circulation model (GCM). Uncertainty quantification is performed using the calibrate-emulate-sample approach, which combines stochastic optimization and machine learning emulation to speed up Bayesian learning. The methods are demonstrated in a perfect-model setting through the calibration and UQ of a convective parameterization in an idealized GCM with a seasonal cycle. Calibration and UQ based on seasonally averaged climate statistics, compared to annually averaged, reduces the calibration error by up to an order of magnitude and narrows the spread of posterior distributions by factors between two and five, depending on the variables used for UQ. The reduction in the size of the parameter posterior distributions leads to a reduction in the uncertainty of climate model predictions.


title: Closed-form existence conditions for band-gap resonances in a finite periodic chain under general boundary conditions

abstract: Bragg scattering in periodic media generates band gaps, frequency bands where waves attenuate rather than propagate. Yet, a finite periodic structure may exhibit resonance frequencies within these band gaps. This is caused by boundary effects introduced by the truncation of the nominal infinite medium. Previous studies of discrete systems determined existence conditions for band-gap resonances, although the focus has been limited to periodic chains with free-free boundaries. In this paper, we present closed-form existence conditions for band-gap resonances in discrete diatomic chains with general boundary conditions (free-free, free-fixed, fixed-free, or fixed-fixed), odd or even chain parity (contrasting or identical masses at the ends), and the possibility of attaching a unique component (mass and/or spring) at one or both ends. The derived conditions are consistent with those theoretically presented or experimentally observed in prior studies of structures that can be modeled as linear discrete diatomic chains with free-free boundary conditions. An intriguing case is a free-free chain with even parity and an arbitrary additional mass at one end of the chain. Introducing such an arbitrary mass underscores a transition among a set of distinct existence conditions depending on the type of the chain boundaries and parity. The proposed analysis is applicable to linear periodic chains in the form of lumped-parameter models as well as continuous granular media models examined in the low-frequency regime.


title: Information-theoretic symmetry classifications of crystal patterns in the presence of noise and strong Fedorov type pseudosymmetries for an optimal subsequent crystallographic processing of these patterns

abstract: Statistically sound crystallographic symmetry classifications are obtained with information theory based methods in the presence of approximately Gaussian distributed noise. A set of three synthetic images with very strong Fedorov type pseudosymmetries and varying amounts of noise serve as examples. The correct distinctions between genuine symmetries and their Fedorov type pseudosymmetry counterparts failed only for the noisiest image of the series where an inconsistent combination of plane symmetry group and projected Laue class was obtained. Contrary to traditional crystallographic symmetry classifications with an image processing program such as CRISP, the classification process does not need to be supervised by a human being. This enables crystallographic symmetry classification of digital images that are more or less periodic in two dimensions (2D) as recorded with sufficient spatial resolution from a wide range of samples with different types of scanning probe microscopes. Alternatives to the employed objective classification methods as proposed by members of the computational symmetry community and machine learning proponents are briefly discussed in an appendix and are found to be wanting because they ignore Fedorov type pseudosymmetries completely. The information theory based methods are more accurate than visual classifications at first sight by most human experts.


title: First Light And Reionisation Epoch Simulations (FLARES) III: The properties of massive dusty galaxies at cosmic dawn

abstract: Using the First Light And Reionisation Epoch Simulations (\textsc{Flares}) we explore the dust driven properties of massive high-redshift galaxies at $z\in[5,10]$. By post-processing the galaxy sample using the radiative transfer code \textsc{skirt} we obtain the full spectral energy distribution. We explore the resultant luminosity functions, IRX-$\beta$ relations as well as the luminosity-weighted dust temperatures in the Epoch of Reionisation (EoR). We find that most of our results are in agreement with the current set of observations, but under-predict the number densities of bright IR galaxies, which are extremely biased towards the most overdense regions. We see that the \textsc{Flares} IRX-$\beta$ relation (for $5\le z\le8$) predominantly follows the local starburst relation. The IRX shows an increase with stellar mass, plateauing at the high-mass end ($\sim10^{10}$M$_{\odot}$) and shows no evolution in the median normalisation with redshift. We also look at the dependence of the peak dust temperature ($T_{\mathrm{peak}}$) on various galaxy properties including the stellar mass, IR luminosity and sSFR, finding the correlation to be strongest with sSFR. The luminosity-weighted dust temperatures increase towards higher redshifts, with the slope of the $T_{\mathrm{peak}}$ - redshift relation showing a higher slope than the lower redshift relations obtained from previous observational and theoretical works. The results from \textsc{Flares}, which is able to provide a better statistical sample of high-redshift galaxies compared to other simulations, provides a distinct vantage point for the high-redshift Universe.


title: Projective Skip-Connections for Segmentation Along a Subset of Dimensions in Retinal OCT

abstract: In medical imaging, there are clinically relevant segmentation tasks where the output mask is a projection to a subset of input image dimensions. In this work, we propose a novel convolutional neural network architecture that can effectively learn to produce a lower-dimensional segmentation mask than the input image. The network restores encoded representation only in a subset of input spatial dimensions and keeps the representation unchanged in the others. The newly proposed projective skip-connections allow linking the encoder and decoder in a UNet-like structure. We evaluated the proposed method on two clinically relevant tasks in retinal Optical Coherence Tomography (OCT): geographic atrophy and retinal blood vessel segmentation. The proposed method outperformed the current state-of-the-art approaches on all the OCT datasets used, consisting of 3D volumes and corresponding 2D en-face masks. The proposed architecture fills the methodological gap between image classification and ND image segmentation.


title: AI Techniques for Software Requirements Prioritization

abstract: Aspects such as limited resources, frequently changing market demands, and different technical restrictions regarding the implementation of software requirements (features) often demand for the prioritization of requirements. The task of prioritization is the ranking and selection of requirements that should be included in future software releases. In this context, an intelligent prioritization decision support is extremely important. The prioritization approaches discussed in this paper are based on different Artificial Intelligence (AI) techniques that can help to improve the overall quality of requirements prioritization processes


title: Adversarial Attacks Against Deep Reinforcement Learning Framework in Internet of Vehicles

abstract: Machine learning (ML) has made incredible impacts and transformations in a wide range of vehicular applications. As the use of ML in Internet of Vehicles (IoV) continues to advance, adversarial threats and their impact have become an important subject of research worth exploring. In this paper, we focus on Sybil-based adversarial threats against a deep reinforcement learning (DRL)-assisted IoV framework and more specifically, DRL-based dynamic service placement in IoV. We carry out an experimental study with real vehicle trajectories to analyze the impact on service delay and resource congestion under different attack scenarios for the DRL-based dynamic service placement application. We further investigate the impact of the proportion of Sybil-attacked vehicles in the network. The results demonstrate that the performance is significantly affected by Sybil-based data poisoning attacks when compared to adversary-free healthy network scenario.


title: Application of pattern spectra and convolutional neural networks to the analysis of simulated Cherenkov Telescope Array data

abstract: The Cherenkov Telescope Array (CTA) will be the next generation gamma-ray observatory and will be the major global instrument for very-high-energy astronomy over the next decade, offering 5 - 10 x better flux sensitivity than current generation gamma-ray telescopes. Each telescope will provide a snapshot of gamma-ray induced particle showers by capturing the induced Cherenkov emission at ground level. The simulation of such events provides images that can be used as training data for convolutional neural networks (CNNs) to determine the energy of the initial gamma rays. Compared to other state-of-the-art algorithms, analyses based on CNNs promise to further enhance the performance to be achieved by CTA. Pattern spectra are commonly used tools for image classification and provide the distributions of the shapes and sizes of various objects comprising an image. The use of relatively shallow CNNs on pattern spectra would automatically select relevant combinations of features within an image, taking advantage of the 2D nature of pattern spectra. In this work, we generate pattern spectra from simulated gamma-ray events instead of using the raw images themselves in order to train our CNN for energy reconstruction. This is different from other relevant learning and feature selection methods that have been tried in the past. Thereby, we aim to obtain a significantly faster and less computationally intensive algorithm, with minimal loss of performance.


title: Effect of small cosmological constant on electromagnetic memory effect

abstract: We consider a generic scattering process that takes place in a region of size R inside the static patch of the de Sitter spacetime such that R is smaller than the curvature length scale of the background. The effect of curvature can thus be studied perturbatively. We obtain the asymptotic electromagnetic field generated by the scattering process including the leading order correction due to the presence of de Sitter background and discuss its universal aspects. We finally caculate the resultant first order corrections to the flat spacetime velocity memory effect.


title: Implementation of the SMART protocol for global qubit control in silicon

abstract: Quantum computing based on spins in the solid state allows for densely-packed arrays of quantum bits. While high-fidelity operation of single qubits has been demonstrated with individual control pulses, the operation of large-scale quantum processors requires a shift in paradigm towards global control solutions. Here we report the experimental implementation of a new type of qubit protocol - the SMART (Sinusoidally Modulated, Always Rotating and Tailored) protocol. As with a dressed qubit, we resonantly drive a two-level system with a continuous microwave field, but here we add a tailored modulation to the dressing field to achieve increased robustness to detuning noise and microwave amplitude fluctuations. We implement this new protocol to control a single spin confined in a silicon quantum dot and confirm the optimal modulation conditions predicted from theory. Universal control of a single qubit is demonstrated using modulated Stark shift control via the local gate electrodes. We measure an extended coherence time of $2$ ms and an average Clifford gate fidelity $>99$ $\%$ despite the relatively long qubit gate times ($>15$ $\unicode[serif]{x03BC}$s, $20$ times longer than a conventional square pulse gate), constituting a significant improvement over a conventional spin qubit and a dressed qubit. This work shows that future scalable spin qubit arrays could be operated using global microwave control and local gate addressability, while maintaining robustness to relevant experimental inhomogeneities.


title: Discovering sparse control strategies in C. elegans

abstract: Biological circuits such as neural or gene regulation networks use internal states to map sensory input to an adaptive repertoire of behavior. Characterizing this mapping is a major challenge for systems biology, and though experiments that probe internal states are developing rapidly, organismal complexity presents a fundamental obstacle given the many possible ways internal states could map to behavior. Using C. elegans as an example, we propose a protocol for systematic perturbation of neural states that limits experimental complexity but still characterizes collective aspects of the neural-behavioral map. We consider experimentally motivated small perturbations -- ones that are most likely to preserve natural dynamics and are closer to internal control mechanisms -- to neural states and their impact on collective neural behavior. Then, we connect such perturbations to the local information geometry of collective statistics, which can be fully characterized using pairwise perturbations. Applying the protocol to a minimal model of C. elegans neural activity, we find that collective neural statistics are most sensitive to a few principal perturbative modes. Dominant eigenvalues decay initially as a power law, unveiling a hierarchy that arises from variation in individual neural activity and pairwise interactions. Highest-ranking modes tend to be dominated by a few, "pivotal" neurons that account for most of the system's sensitivity, suggesting a sparse mechanism for control of collective behavior.


title: Particle motion in circularly polarized vacuum pp waves

abstract: Bialynicki-Birula and Charzynski argued that a gravitational wave emitted during the merger of a black hole binary may be approximated by a circularly polarized wave which may in turn trap particles [1]. In this paper we consider particle motion in a class of gravitational waves which includes, besides circularly polarized periodic waves (CPP) [2], also the one proposed by Lukash [3] to study anisotropic cosmological models. Both waves have a 7-parameter conformal symmetry which contains, in addition to the generic 5-parameter (broken) Carroll group, also a 6th isometry. The Lukash spacetime can be transformed by a conformal rescaling of time to a perturbed CPP problem. Bounded geodesics, found both analytically and numerically, arise when the Lukash wave is of Bianchi type VI. Their symmetries can also be derived from the Lukash-CPP relation. Particle trapping is discussed.