title: Gems of the Galaxy Zoos -- a Wide-Ranging Hubble Space Telescope Gap-Filler Program

abstract: We describe the Gems of the Galaxy Zoos (Zoo Gems) project, a gap-filler project using short windows in the Hubble Space Telescope's schedule. As with previous snapshot programs, targets are taken from a pool based on position; we combine objects selected by volunteers in both the Galaxy Zoo and Radio Galaxy Zoo citizen-science projects. Zoo Gems uses exposures with the Advanced Camera for Surveys (ACS) to address a broad range of topics in galaxy morphology, interstellar-medium content, host galaxies of active galactic nuclei, and galaxy evolution. Science cases include studying galaxy interactions, backlit dust in galaxies, post-starburst systems, rings and peculiar spiral patterns, outliers from the usual color-morphology relation, Green Pea compact starburst systems, double radio sources with spiral host galaxies, and extended emission-line regions around active galactic nuclei. For many of these science categories, final selection of targets from a larger list used public input via a voting process. Highlights to date include the prevalence of tightly-wound spiral structure in blue, apparently early-type galaxies, a nearly complete Einstein ring from a group lens, redder components at lower surface brightness surrounding compact Green Pea starbursts, and high-probability examples of spiral galaxies hosting large double radio sources.


title: On Lyapunov Stability of Positive and Conservative Time Integrators and Application to Second Order Modified Patankar--Runge--Kutta Schemes

abstract: Since almost twenty years, modified Patankar--Runge--Kutta (MPRK) methods have proven to be efficient and robust numerical schemes that preserve positivity and conservativity of the production-destruction system irrespectively of the time step size chosen. Due to these advantageous properties they are used for a wide variety of applications. Nevertheless, until now, an analytic investigation of the stability of MPRK schemes is still missing, since the usual approach by means of Dahlquist's equation is not feasible. Therefore, we consider a positive and conservative 2D test problem and provide statements usable for a stability analysis of general positive and conservative time integrator schemes based on the center manifold theory. We use this approach to investigate the Lyapunov stability of the second order MPRK22($\alpha$) and MPRK22ncs($\alpha$) schemes. We prove that MPRK22($\alpha$) schemes are unconditionally stable and derive the stability regions of MPRK22ncs($\alpha$) schemes. Finally, numerical experiments are presented, which confirm the theoretical results.


title: Exact Privacy Analysis of the Gaussian Sparse Histogram Mechanism

abstract: Sparse histogram methods can be useful for returning differentially private counts of items in large or infinite histograms, large group-by queries, and more generally, releasing a set of statistics with sufficient item counts. We consider the Gaussian version of the sparse histogram mechanism and study the exact $\epsilon,\delta$ differential privacy guarantees satisfied by this mechanism. We compare these exact $\epsilon,\delta$ parameters to the simpler overestimates used in prior work to quantify the impact of their looser privacy bounds.


title: Growth optimization and device integration of narrow-bandgap graphene nanoribbons

abstract: The electronic, optical and magnetic properties of graphene nanoribbons (GNRs) can be engineered by controlling their edge structure and width with atomic precision through bottom-up fabrication based on molecular precursors. This approach offers a unique platform for all-carbon electronic devices but requires careful optimization of the growth conditions to match structural requirements for successful device integration, with GNR length being the most critical parameter. In this work, we study the growth, characterization, and device integration of 5-atom wide armchair GNRs (5-AGNRs), which are expected to have an optimal band gap as active material in switching devices. 5-AGNRs are obtained via on-surface synthesis under ultra-high vacuum conditions from Br- and I-substituted precursors. We show that the use of I-substituted precursors and the optimization of the initial precursor coverage quintupled the average 5-AGNR length. This significant length increase allowed us to integrate 5-AGNRs into devices and to realize the first field-effect transistor based on narrow bandgap AGNRs that shows switching behavior at room temperature. Our study highlights that optimized growth protocols can successfully bridge between the sub-nanometer scale, where atomic precision is needed to control the electronic properties, and the scale of tens of nanometers relevant for successful device integration of GNRs.


title: System Identification with Variance Minimization via Input Design

abstract: The subspace method is one of the mainstream system identification method of linear systems, and its basic idea is to estimate the system parameter matrices by projecting them into a subspace related to input and output. However, most of the existing subspace methods cannot have the statistic performance guaranteed since the lack of closed-form expression of the estimation. Meanwhile, traditional subspace methods cannot deal with the uncertainty of the noise, and thus stable identification results cannot be obtained. In this paper, we propose a novel improved subspace method from the perspective of input design, which guarantees the consistent and stable identification results with the minimum variance. Specifically, we first obtain a closed-form estimation of the system matrix, then analyze the statistic performance by deriving the maximum identification deviation. This identification deviation maximization problem is non-convex, and is solved by splitting it into two sub-problems with the optimal solution guaranteed. Next, an input design method is proposed to deal with the uncertainty and obtain stable identification results by minimizing the variance. This problem is formulated as a constrained min-max optimization problem. The optimal solution is obtained from transforming the cost function into a convex function while ensuring the safety constraints through the method of predictive control. We prove the consistency and the convergence of the proposed method. Simulation demonstrates the effectiveness of our method.


title: A New Temporal Interpretation of Cluster Editing

abstract: The NP-complete graph problem Cluster Editing seeks to transform a static graph into a disjoint union of cliques by making the fewest possible edits to the edges. We introduce a natural interpretation of this problem in temporal graphs, whose edge sets change over time. This problem is NP-complete even when restricted to temporal graphs whose underlying graph is a path, but we obtain two polynomial-time algorithms for restricted cases. In the static setting, it is well-known that a graph is a disjoint union of cliques if and only if it contains no induced copy of $P_3$; we demonstrate that no general characterisation involving sets of at most four vertices can exist in the temporal setting, but obtain a complete characterisation involving forbidden configurations on at most five vertices. This characterisation gives rise to an FPT algorithm parameterised simultaneously by the permitted number of modifications and the lifetime of the temporal graph.


title: Optimal fenestration of the Fontan circulation

abstract: In this paper, we develop a pulsatile compartmental model of the Fontan circulation and use it to explore the effects of a fenestration added to this physiology. A fenestration is a shunt between the systemic and pulmonary veins that is added either at the time of Fontan conversion or at a later time for the treatment of complications. This shunt increases cardiac output and decreases systemic venous pressure. However, these hemodynamic benefits are achieved at the expense of a decrease in the arterial oxygen saturation. The model developed this paper incorporates fenestration size as a parameter and describes both blood flow and oxygen transport. It is calibrated to clinical data from Fontan patients, and we use it to study the impact of a fenestration on several hemodynamic variables. In certain scenarios corresponding to high-risk Fontan physiology, we demonstrate the existence of an optimal fenestration size that maximizes oxygen delivery to the systemic tissues.


title: Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives

abstract: This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program.


title: The string little algebra

abstract: We consider the string, like point particles and branes, to be an irreducible representation of the semi-direct product of the Cartan involution invariant subalgebra of E11 and its vector representation. We show that the subalgebra that preserves the string charges, the string little algebra, is essentially the Borel subalgebra of E9. We also show that the known string physical states carry a representation of parts of this algebra.


title: Keyword localisation in untranscribed speech using visually grounded speech models

abstract: Keyword localisation is the task of finding where in a speech utterance a given query keyword occurs. We investigate to what extent keyword localisation is possible using a visually grounded speech (VGS) model. VGS models are trained on unlabelled images paired with spoken captions. These models are therefore self-supervised -- trained without any explicit textual label or location information. To obtain training targets, we first tag training images with soft text labels using a pretrained visual classifier with a fixed vocabulary. This enables a VGS model to predict the presence of a written keyword in an utterance, but not its location. We consider four ways to equip VGS models with localisations capabilities. Two of these -- a saliency approach and input masking -- can be applied to an arbitrary prediction model after training, while the other two -- attention and a score aggregation approach -- are incorporated directly into the structure of the model. Masked-based localisation gives some of the best reported localisation scores from a VGS model, with an accuracy of 57% when the system knows that a keyword occurs in an utterance and need to predict its location. In a setting where localisation is performed after detection, an $F_1$ of 25% is achieved, and in a setting where a keyword spotting ranking pass is first performed, we get a localisation P@10 of 32%. While these scores are modest compared to the idealised setting with unordered bag-of-word-supervision (from transcriptions), these models do not receive any textual or location supervision. Further analyses show that these models are limited by the first detection or ranking pass. Moreover, individual keyword localisation performance is correlated with the tagging performance from the visual classifier. We also show qualitatively how and where semantic mistakes occur, e.g. that the model locates surfer when queried with ocean.


title: Learning to reason about and to act on physical cascading events

abstract: Reasoning and interacting with dynamic environments is a fundamental problem in AI, but it becomes extremely challenging when actions can trigger cascades of cross-dependent events. We introduce a new supervised learning setup called {\em Cascade} where an agent is shown a video of a physically simulated dynamic scene, and is asked to intervene and trigger a cascade of events, such that the system reaches a "counterfactual" goal. For instance, the agent may be asked to "Make the blue ball hit the red one, by pushing the green ball". The agent intervention is drawn from a continuous space, and cascades of events makes the dynamics highly non-linear. We combine semantic tree search with an event-driven forward model and devise an algorithm that learns to search in semantic trees in continuous spaces. We demonstrate that our approach learns to effectively follow instructions to intervene in previously unseen complex scenes. It can also reason about alternative outcomes, when provided an observed cascade of events.


title: Experimental analysis of energy transfers between a quantum emitter and light fields

abstract: Energy transfer between quantum systems can either be achieved through an effective unitary interaction or through the generation of entanglement. This observation defines two types of energy exchange: unitary and correlation energy. Here we propose and implement experimental protocols to access these energy transfers in interactions between a quantum emitter and light fields. Upon spontaneous emission, we measure the unitary energy transfer from the emitter to the optical field and show that it never exceeds half of the total energy and is reduced when introducing decoherence. We then study the interference of the emitted field and a laser field at a beam splitter and show that the energy transfers quantitatively depend on the quantum purity of the emitted field.


title: A Survey on Retrieval-Augmented Text Generation

abstract: Recently, retrieval-augmented text generation attracted increasing attention of the computational linguistics community. Compared with conventional generation models, retrieval-augmented text generation has remarkable advantages and particularly has achieved state-of-the-art performance in many NLP tasks. This paper aims to conduct a survey about retrieval-augmented text generation. It firstly highlights the generic paradigm of retrieval-augmented generation, and then it reviews notable approaches according to different tasks including dialogue response generation, machine translation, and other generation tasks. Finally, it points out some important directions on top of recent methods to facilitate future research.


title: Hypergraph Fuss-Catalan Numbers

abstract: The Catalan numbers $C_n$ are an extremely well-studied sequence of numbers that appear as the answer to many combinatorial problems. Two generalizations of these numbers that have been studied are the Fuss-Catalan numbers and the Hypergraph Catalan numbers. In this paper, we study the combination of these, the Hypergraph Fuss-Catalan numbers. We provide some combinatorial interpretations of these numbers, as well as describe their generating function.


title: On Joint Communication and Channel Discrimination

abstract: We consider a basic communication and sensing setup comprising a transmitter, a receiver and a sensor. The transmitter sends an encoded sequence to the receiver through a discrete memoryless channel, and the receiver is interested in decoding the sequence. On the other hand, the sensor picks up a noisy version of the transmitted sequence through one of two possible discrete memoryless channels. The sensor knows the transmitted sequence and wishes to discriminate between the two possible channels, i.e. to identify the channel that has generated the output given the input. We study the trade-off between communication and sensing in the asymptotic regime, captured in terms of the coding rate to the receiver against the discrimination error exponent at the sensor. We characterize the optimal rate-exponent trade-off for general discrete memoryless channels with an input cost constraint.


title: Tailoring Gradient Methods for Differentially-Private Distributed Optimization

abstract: Decentralized optimization is gaining increased traction due to its widespread applications in large-scale machine learning and multi-agent systems. The same mechanism that enables its success, i.e., information sharing among participating agents, however, also leads to the disclosure of individual agents' private information, which is unacceptable when sensitive data are involved. As differential privacy is becoming a de facto standard for privacy preservation, recently results have emerged integrating differential privacy with distributed optimization. However, directly incorporating differential privacy design in existing distributed optimization approaches significantly compromises optimization accuracy. In this paper, we propose to redesign and tailor gradient methods for differentially-private distributed optimization, and propose two differential-privacy oriented gradient methods that can ensure both rigorous epsilon-differential privacy and optimality. The first algorithm is based on static-consensus based gradient methods, and the second algorithm is based on dynamic-consensus (gradient-tracking) based distributed optimization methods and, hence, is applicable to general directed interaction graph topologies. Both algorithms can simultaneously ensure almost sure convergence to an optimal solution and a finite privacy budget, even when the number of iterations goes to infinity. To our knowledge, this is the first time that both goals are achieved simultaneously. Numerical simulations using a distributed estimation problem and experimental results on a benchmark dataset confirm the effectiveness of the proposed approaches.


title: Sumsets and Veronese varieties

abstract: In this paper, to any subset $\mathcal{A} \subset \mathbb{Z}^{n}$ we explicitly associate a unique monomial projection $Y_{n,d_{\mathcal{A}}}$ of a Veronese variety, whose Hilbert function coincides with the cardinality of the $t$-fold sumsets $t\mathcal{A}$. This link allows us to tackle the classical problem of determining the polynomial $p_{\mathcal{A}} \in \mathbb{Q}[t]$ such that $|t\mathcal{A}| = p_{\mathcal{A}}(t)$ for all $t \geq t_0$ and the minimum integer $n_0(\mathcal{A}) \leq t_0$ for which this condition is satisfied, i.e. the so-called {\em phase transition} of $|t\mathcal{A}|$. We use the Castelnuovo--Mumford regularity and the geometry of $Y_{n,d_{\mathcal{A}}}$ to describe the polynomial $p_{\mathcal{A}}(t)$ and to derive new bounds for $n_0(\mathcal{A})$ under some technical assumptions on the convex hull of $\mathcal{A}$; and vice versa we apply the theory of sumsets to obtain geometric information of the varieties $Y_{n,d_{\mathcal{A}}}$.


title: NeuRegenerate: A Framework for Visualizing Neurodegeneration

abstract: Recent advances in high-resolution microscopy have allowed scientists to better understand the underlying brain connectivity. However, due to the limitation that biological specimens can only be imaged at a single timepoint, studying changes to neural projections is limited to general observations using population analysis. In this paper, we introduce NeuRegenerate, a novel end-to-end framework for the prediction and visualization of changes in neural fiber morphology within a subject, for specified age-timepoints.To predict projections, we present neuReGANerator, a deep-learning network based on cycle-consistent generative adversarial network (cycleGAN) that translates features of neuronal structures in a region, across age-timepoints, for large brain microscopy volumes. We improve the reconstruction quality of neuronal structures by implementing a density multiplier and a new loss function, called the hallucination loss.Moreover, to alleviate artifacts that occur due to tiling of large input volumes, we introduce a spatial-consistency module in the training pipeline of neuReGANerator. We show that neuReGANerator has a reconstruction accuracy of 94% in predicting neuronal structures. Finally, to visualize the predicted change in projections, NeuRegenerate offers two modes: (1) neuroCompare to simultaneously visualize the difference in the structures of the neuronal projections, across the age timepoints, and (2) neuroMorph, a vesselness-based morphing technique to interactively visualize the transformation of the structures from one age-timepoint to the other. Our framework is designed specifically for volumes acquired using wide-field microscopy. We demonstrate our framework by visualizing the structural changes in neuronal fibers within the cholinergic system of the mouse brain between a young and old specimen.


title: An Optimal Transport Perspective on Unpaired Image Super-Resolution

abstract: Real-world image super-resolution (SR) tasks often do not have paired datasets, which limits the application of supervised techniques. As a result, the tasks are usually approached by unpaired techniques based on Generative Adversarial Networks (GANs), which yield complex training losses with several regularization terms, e.g., content or identity losses. We theoretically investigate optimization problems which arise in such models and find two surprizing observations. First, the learned SR map is always an optimal transport (OT) map. Second, we theoretically prove and empirically show that the learned map is biased, i.e., it does not actually transform the distribution of low-resolution images to high-resolution ones. Inspired by these findings, we propose an algorithm for unpaired SR which learns an unbiased OT map for the perceptual transport cost. Unlike the existing GAN-based alternatives, our algorithm has a simple optimization objective reducing the need for complex hyperparameter selection and an application of additional regularizations. At the same time, it provides a nearly state-of-the-art performance on the large-scale unpaired AIM19 dataset.


title: An Eye for an Eye: Defending against Gradient-based Attacks with Gradients

abstract: Deep learning models have been shown to be vulnerable to adversarial attacks. In particular, gradient-based attacks have demonstrated high success rates recently. The gradient measures how each image pixel affects the model output, which contains critical information for generating malicious perturbations. In this paper, we show that the gradients can also be exploited as a powerful weapon to defend against adversarial attacks. By using both gradient maps and adversarial images as inputs, we propose a Two-stream Restoration Network (TRN) to restore the adversarial images. To optimally restore the perturbed images with two streams of inputs, a Gradient Map Estimation Mechanism is proposed to estimate the gradients of adversarial images, and a Fusion Block is designed in TRN to explore and fuse the information in two streams. Once trained, our TRN can defend against a wide range of attack methods without significantly degrading the performance of benign inputs. Also, our method is generalizable, scalable, and hard to bypass. Experimental results on CIFAR10, SVHN, and Fashion MNIST demonstrate that our method outperforms state-of-the-art defense methods.


title: On Linear Separability under Linear Compression with Applications to Hard Support Vector Machine

abstract: This paper investigates the theoretical problem of maintaining linear separability of the data-generating distribution under linear compression. While it has been long known that linear separability may be maintained by linear transformations that approximately preserve the inner products between the domain points, the limit to which the inner products are preserved in order to maintain linear separability was unknown. In this paper, we show that linear separability is maintained as long as the distortion of the inner products is smaller than the squared margin of the original data-generating distribution. The proof is mainly based on the geometry of hard support vector machines (SVM) extended from the finite set of training examples to the (possibly) infinite domain of the data-generating distribution. As applications, we derive bounds on the (i) compression length of random sub-Gaussian matrices; and (ii) generalization error for compressive learning with hard-SVM.


title: Diffusion-mediated surface reactions and stochastic resetting

abstract: In this paper, we investigate the effects of stochastic resetting on diffusion in $\R^d\backslash \calU$, where $\calU$ is a bounded obstacle with a partially absorbing surface $\partial \calU$. We begin by considering a Robin boundary condition with a constant reactivity $\kappa_0$, and show how previous results are recovered in the limits $\kappa_0\rightarrow 0,\infty$. We then generalize the Robin boundary condition to a more general probabilistic model of diffusion-mediated surface reactions using an encounter-based approach. The latter considers the joint probability density or propagator $P(\x,\ell,t|\x_0)$ for the pair $(\X_t,\ell_t)$ in the case of a perfectly reflecting surface, where $\X_t$ and $\ell_t$ denote the particle position and local time, respectively. The local time determines the amount of time that a Brownian particle spends in a neighborhood of the boundary. The effects of surface reactions are then incorporated via an appropriate stopping condition for the boundary local time. We construct the boundary value problem (BVP) satisfied by the propagator in the presence of resetting, and use this to derive implicit equations for the marginal density of particle position and the survival probability. We highlight the fact that these equations are difficult to solve in the case of non-constant reactivities, since resetting is not governed by a renewal process. We then consider a simpler problem in which both the position and local time are reset. In this case, the survival probability with resetting can be expressed in terms of the survival probability without resetting, which allows us to explore the dependence of the MFPT on the resetting rate $r$ and the type of surface reactions. The theory is illustrated using the example of a spherically symmetric surface.


title: Brain Source Localization by Alternating Projection

abstract: We present a novel solution to the problem of localizing magnetoencephalography (MEG) and electroencephalography (EEG) brain signals. The solution is sequential and iterative, and is based on minimizing the least-squares criterion by the Alternating Projection algorithm. Results from simulated and experimental MEG data from a human subject demonstrated robust performance, with consistently superior localization accuracy than scanning methods belonging to the beamformer and multiple-signal classification (MUSIC) families. Importantly, the proposed solution is more robust to forward model errors resulting from head rotations and translations, with a significant advantage in highly correlated sources.


title: Distributed Ride-Matching for Shared Ridehailing Service with Intelligent City Infrastructure

abstract: High computational time is one of the most important operational issues in centralized dynamic shared ridehailing services. To resolve this issue, we propose a distributed ride-matching system that is based on vehicle to infrastructure (V2I) and infrastructure to infrastructure (I2I) communication. The application on downtown Toronto road network demonstrated that the distributed system resulted in a speed-up of 125 times in terms of computational time and showed high scalability. Moreover, the service rate in the proposed system improved by 7% compared to the centralized. However, the centralized system showed 29% and 17% improvement in wait time and detour time, respectively.


title: The ocean fine spray

abstract: A major fraction of the atmospheric aerosols come from the ocean spray originated by the bursting of bubbles from breaking waves. A theoretical framework that incorporates the latest knowledge on film and jet droplets from bubble bursting is proposed. Assuming that their relics constitute the ultimate origin of primary and secondary sea aerosols through a diversity of physicochemical routes, the model can be reduced to a single controlling parameter to predict the global probability density distribution (pdf) of the ocean spray. The bursting and collapse of small bubbles on the sea surface from about 10 to 100 microns produces an extreme energy focusing and the ejection of a rapid liquid spout whose size reaches the free molecular regime of the gaseous environment. In these rarefied conditions, simulations show that this spout yields a jet of sub-micrometer and nanometric scale droplets whose number and speed can be far beyond any previous estimation, overcoming by orders of magnitude alternative mechanisms recently proposed. The one-parameter model fits remarkably well published experimental measurements along five orders of magnitude of spray size, from about 5 nm to about 0.5 mm. According to this proposal, the majority of aerosols determining the life on our planet would have their extremely elusive birth in the uterus-like nano-shape of small bursting bubbles on the ocean surface at the very latest instants of collapse.


title: An ASP approach for reasoning on neural networks under a finitely many-valued semantics for weighted conditional knowledge bases

abstract: Weighted knowledge bases for description logics with typicality have been recently considered under a "concept-wise" multipreference semantics (in both the two-valued and fuzzy case), as the basis of a logical semantics of MultiLayer Perceptrons (MLPs). In this paper we consider weighted conditional ALC knowledge bases with typicality in the finitely many-valued case, through three different semantic constructions. For the boolean fragment LC of ALC we exploit ASP and "asprin" for reasoning with the concept-wise multipreference entailment under a phi-coherent semantics, suitable to characterize the stationary states of MLPs. As a proof of concept, we experiment the proposed approach for checking properties of trained MLPs. The paper is under consideration for acceptance in TPLP.


title: Time delay interferometry without clock synchronisation

abstract: Time-delay interferometry (TDI) is a data processing technique for LISA designed to suppress the otherwise overwhelming laser noise by several orders of magnitude. It is widely believed that TDI can only be applied once all phase or frequency measurements from each spacecraft have been synchronized to a common time frame. We demonstrate analytically, using as an example the commonly-used Michelson combination X, that TDI can be computed using the raw, unsynchronized data, thereby avoiding the need for an initial synchronization processing step and significantly simplifying the initial noise reduction pipeline. Furthermore, the raw data is free of any potential artifacts introduced by clock synchronization and reference frame transformation algorithms, which allows to operate directly on the MHz beatnotes. As a consequence, in-band clock noise is directly suppressed as part of TDI, in contrast to the approach previously proposed in the literature (in which large trends in the beatnotes are removed before the main laser-noise reduction step, and clock noise is suppressed in an extra processing step). We validate our algorithm with full-scale numerical simulations that use LISA Instrument and PyTDI and show that we reach the same performance levels as the previously proposed methods, ultimately limited by the clock sideband stability.


title: GLISp-r: A preference-based optimization algorithm with convergence guarantees

abstract: Preference-based optimization algorithms are iterative procedures that seek the optimal calibration of a decision vector based only on comparisons between couples of different tunings. At each iteration, a human decision-maker expresses a preference between two calibrations (samples), highlighting which one, if any, is better than the other. The optimization procedure must use the observed preferences to find the tuning of the decision vector that is most preferred by the decision-maker, while also minimizing the number of comparisons. In this work, we formulate the preference-based optimization problem from a utility theory perspective. Then, we propose GLISp-r, an extension of a recent preference-based optimization procedure called GLISp. The latter uses a Radial Basis Function surrogate to describe the tastes of the decision-maker. Iteratively, GLISp proposes new samples to compare with the best calibration available by trading off exploitation of the surrogate model and exploration of the decision space. In GLISp-r, we propose a different criterion to use when looking for new candidate samples that is inspired by MSRS, a popular procedure in the black-box optimization framework. Compared to GLISp, GLISp-r is less likely to get stuck on local optima of the preference-based optimization problem. We motivate this claim theoretically, with a proof of global convergence, and empirically, by comparing the performances of GLISp and GLISp-r on several benchmark optimization problems.


title: Neutrino point source searches for dark matter spikes

abstract: Any dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation, leading to a potentially detectable neutrino signal. In this paper we examine $10-10^5 M_\odot$ black holes associated with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today, in light of neutrino data from the ANTARES and IceCube detectors. In various regions of the sky, we determine the minimum distance away from the solar system that a dark matter spike must be in order to have not been detected as a neutrino point source for a variety of representative dark matter annihilation channels. Given these constraints on the distribution of dark matter spikes in the Galaxy, we place significant limits on the formation of the first generation of stars in early minihalos -- stronger than previous limits from gamma-ray searches in Fermi Gamma-Ray Space Telescope data. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted; thus neutrino observations may be used to constrain the properties of Dark Stars. The limits are particularly strong for heavier WIMPs. For WIMP masses $\sim 5 \,$TeV, we show that $\lesssim 10 \%$ of minihalos can host first stars that collapse into BHs larger than $10^3 M_\odot$.


title: A first order description of a nonlinear SPDE in the spirit of rough paths

abstract: We consider a nonlinear stochastic partial differential equation (SPDE) in divergence form where the forcing term is a Gaussian noise, that is white in time and colored in space such that the gradient of the solution is H\"older-continuous, but not differentiable. Then, we prove a generalized Taylor expansion of the difference between the solution to the SPDE and the solution to its linearization around a fixed basepoint. The result is reminiscent of the theory of (controlled) rough paths and agrees with the general observation, that, in settings with a rough driver, subtracting the solution to the linearized equation yields a more regular object.


title: Language Models Explain Word Reading Times Better Than Empirical Predictability

abstract: Though there is a strong consensus that word length and frequency are the most important single-word features determining visual-orthographic access to the mental lexicon, there is less agreement as how to best capture syntactic and semantic factors. The traditional approach in cognitive reading research assumes that word predictability from sentence context is best captured by cloze completion probability (CCP) derived from human performance data. We review recent research suggesting that probabilistic language models provide deeper explanations for syntactic and semantic effects than CCP. Then we compare CCP with (1) Symbolic n-gram models consolidate syntactic and semantic short-range relations by computing the probability of a word to occur, given two preceding words. (2) Topic models rely on subsymbolic representations to capture long-range semantic similarity by word co-occurrence counts in documents. (3) In recurrent neural networks (RNNs), the subsymbolic units are trained to predict the next word, given all preceding words in the sentences. To examine lexical retrieval, these models were used to predict single fixation durations and gaze durations to capture rapidly successful and standard lexical access, and total viewing time to capture late semantic integration. The linear item-level analyses showed greater correlations of all language models with all eye-movement measures than CCP. Then we examined non-linear relations between the different types of predictability and the reading times using generalized additive models. N-gram and RNN probabilities of the present word more consistently predicted reading performance compared with topic models or CCP.


title: Structure-preserving GANs

abstract: Generative adversarial networks (GANs), a class of distribution-learning methods based on a two-player game between a generator and a discriminator, can generally be formulated as a minmax problem based on the variational representation of a divergence between the unknown and the generated distributions. We introduce structure-preserving GANs as a data-efficient framework for learning distributions with additional structure such as group symmetry, by developing new variational representations for divergences. Our theory shows that we can reduce the discriminator space to its projection on the invariant discriminator space, using the conditional expectation with respect to the sigma-algebra associated to the underlying structure. In addition, we prove that the discriminator space reduction must be accompanied by a careful design of structured generators, as flawed designs may easily lead to a catastrophic "mode collapse" of the learned distribution. We contextualize our framework by building symmetry-preserving GANs for distributions with intrinsic group symmetry, and demonstrate that both players, namely the equivariant generator and invariant discriminator, play important but distinct roles in the learning process. Empirical experiments and ablation studies across a broad range of data sets, including real-world medical imaging, validate our theory, and show our proposed methods achieve significantly improved sample fidelity and diversity -- almost an order of magnitude measured in Fr\'echet Inception Distance -- especially in the small data regime.


title: Well-posedness and dynamics of solutions to the generalized KdV with low power nonlinearity

abstract: We consider two types of the generalized Korteweg - de Vries equation, where the nonlinearity is given with or without absolute values, and, in particular, including the low powers of nonlinearity, an example of which is the Schamel equation. We first prove the local well-posedness of both equations in a weighted subspace of $H^1$ that includes functions with polynomial decay, extending the result of Linares et al [39] to fractional weights. We then investigate solutions numerically, confirming the well-posedness and extending it to a wider class of functions that includes exponential decay. We include a comparison of solutions to both types of equations, in particular, we investigate soliton resolution for the positive and negative data with different decay rates. Finally, we study the interaction of various solitary waves in both models, showing the formation of solitons, dispersive radiation and even breathers, all of which are easier to track in nonlinearities with lower power.


title: NANOGrav Signal from the End of Inflation and the LIGO Mass and Heavier Primordial Black Holes

abstract: Releasing the 12.5-year pulsar timing array data, the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has recently reported the evidence for a stochastic common-spectrum which would herald the detection of a stochastic gravitational wave background (SGWB) for the first time. We investigate if the signal could be generated from the end of a $\sim 10$ MeV but still phenomenologically viable double-field inflation when the field configuration settles to its true vacuum. During the double-field inflation at such scales, bubbles of true vacuum that can collapse to LIGO mass and heavier primordial black holes form. We show that only when this process happens with a first-order phase transition, the produced gravitational wave spectrum can match with the NANOGrav acclaimed SGWB signal. We show that the produced gravitational wave spectrum matches the NANOGrav SGWB signal only when this process happens through a first-order phase transition. Using LATTICEEASY, we also examine the previous observation in the literature that by lowering the scale of preheating, despite the shift of the peak frequency of the gravitational wave profile to smaller values, the amplitude of the SGWB could be kept almost constant. We notice that this observation breaks down at the preheating scale, $M\lesssim 10^{-14}~m_{{}_{\rm Pl}}$.


title: Exact general solutions for cosmological scalar field evolution in a background-dominated expansion

abstract: We derive exact general solutions (as opposed to attractor particular solutions) and corresponding first integrals for the evolution of a scalar field $\phi$ in a universe dominated by a background fluid with equation of state parameter $w_B$. In addition to the previously-examined linear [$V(\phi) = V_0 \phi$] and quadratic [$V(\phi) = V_0 \phi^2$] potentials, we show that exact solutions exist for the power law potential $V(\phi) = V_0 \phi^n$ with $n = 4(1+w_B)/(1-w_B) + 2$ and $n = 2(1+w_B)/(1-w_B)$. These correspond to the potentials $V(\phi) = V_0 \phi^6$ and $V(\phi) = V_0 \phi^2$ for matter domination and $V(\phi) = V_0 \phi^{10}$ and $V(\phi) = V_0 \phi^4$ for radiation domination. The $\phi^6$ and $\phi^{10}$ potentials can yield either oscillatory or non-oscillatory evolution, and we use the first integrals to determine how the initial conditions map onto each form of evolution. The exponential potential yields an exact solution for a stiff/kination ($w_B = 1$) background. We use this exact solution to derive an analytic expression for the evolution of the equation of state parameter, $w_\phi$, for this case.

categories: gr-qc astro-ph.CO hep-th math-ph math.MP


title: A combined inelastic neutron scattering and \textit{ab initio} lattice dynamics study of FeSi

abstract: The phonon renormalization across the semiconductor-to-metal crossover in FeSi is investigated by inelastic neutron scattering combined with \textit{ab-initio} lattice dynamical calculations. A significant part of reciprocal space with a particular focus on the 110$-$001 scattering plane is mapped by the time-of-flight inelastic neutron scattering data taken below and above the crossover. Individual momentum values are investigated in more detail as a function of temperature. The data reveal that the anomalous phonon softening upon metallization is not exclusive to the high symmetry $R$ and $\Gamma$ points. Several other phonon modes around the $R$-point as well as the phonon modes at the $M$ and $X$ points of the Brillouin zone exhibit anomalous phonon softening with magnitudes comparable to that observed at the $R$-point. The momentum dependence of the phonon softening is reproduced by the lattice dynamical calculation based on the density functional perturbation theory. We discuss our findings with respect to the nature of the semiconductor-to-metal crossover in FeSi, for which different microscopic origins have been proposed, i.e., lattice thermal disorder and electronic correlation effects.


title: Ultra-Wideband Teach and Repeat

abstract: Autonomously retracing a manually-taught path is desirable for many applications, and Teach and Repeat (T&R) algorithms present an approach that is suitable for long-range autonomy. In this paper, ultra-wideband (UWB) ranging-based T&R is proposed for vehicles with limited resources. By fixing single UWB transceivers at far-apart unknown locations in an indoor environment, a robot with 3 UWB transceivers builds a locally consistent map during the teach pass by fusing the range measurements under a custom ranging protocol with an on-board IMU and height measurements. The robot then uses information from the teach pass to retrace the same trajectory autonomously. The proposed ranging protocol and T&R algorithm are validated in simulation, where it is shown that the robot can successfully retrace the taught trajectory with sub-metre tracking error.


title: THz radiation induced circular Hall effect in graphene

abstract: We report on the observation of the circular transversal terahertz photoconductivity in monolayer graphene supplied by a back gate. The photoconductivity response is caused by the free carrier absorption and reverses its sign upon switching the radiation helicity. The observed dc Hall effect manifests the time inversion symmetry breaking induced by circularly polarized terahertz radiation in the absence of a magnetic field. For low gate voltages, the photosignal is found to be proportional to the radiation intensity and can be ascribed to the alignment of electron momenta by the combined action of THz and static electric fields as well as by the dynamic heating and cooling of the electron gas. Strikingly, at high gate voltages, we observe that the linear-in-intensity Hall photoconductivity vanishes; the photoresponse at low intensities becomes superlinear and varies with the square of the radiation intensity. We attribute this behavior to the interplay of the second- and fourth-order effects in the radiation electric field which has not been addressed theoretically so far and requires additional studies.


title: Probabilistically Robust Learning: Balancing Average- and Worst-case Performance

abstract: Many of the successes of machine learning are based on minimizing an averaged loss function. However, it is well-known that this paradigm suffers from robustness issues that hinder its applicability in safety-critical domains. These issues are often addressed by training against worst-case perturbations of data, a technique known as adversarial training. Although empirically effective, adversarial training can be overly conservative, leading to unfavorable trade-offs between nominal performance and robustness. To this end, in this paper we propose a framework called probabilistic robustness that bridges the gap between the accurate, yet brittle average case and the robust, yet conservative worst case by enforcing robustness to most rather than to all perturbations. From a theoretical point of view, this framework overcomes the trade-offs between the performance and the sample-complexity of worst-case and average-case learning. From a practical point of view, we propose a novel algorithm based on risk-aware optimization that effectively balances average- and worst-case performance at a considerably lower computational cost relative to adversarial training. Our results on MNIST, CIFAR-10, and SVHN illustrate the advantages of this framework on the spectrum from average- to worst-case robustness.


title: Filter Bubble effect in the multistate voter model

abstract: Social media influence online activity by recommending to users content strongly correlated with what they have preferred in the past. In this way they constrain users within filter bubbles that strongly limit their exposure to new or alternative content. We investigate this type of dynamics by considering a multistate voter model where, with a given probability $\lambda$, a user interacts with a "personalized information" suggesting the opinion most frequently held in the past. By means of theoretical arguments and numerical simulations, we show the existence of a nontrivial transition between a region (for small $\lambda$) where consensus is reached and a region (above a threshold $\lambda_c$) where the system gets polarized and clusters of users with different opinions persist indefinitely. The threshold always vanishes for large system size $N$, showing that consensus becomes impossible for a large number of users. This finding opens new questions about the side effects of the widespread use of personalized recommendation algorithms.


title: On Igusa local zeta functions of Hauser hybrid polynomials

abstract: Let $K$ be a local field and $f(x)\in K[x]$ be a non-constant polynomial. When ${\rm char}K=0$, Igusa showed the local zeta function is a rational function. However, when ${\rm char}K>0$, the rationality of the local zeta function is unknown in general. In this paper, we study the local zeta functions for the so-called hybrid polynomials in three variables with coefficients in a non-archimedean local field of positive characteristic. These hybrid polynomials were first introduced by Hauser in 2003 to study the resolution of singularities in positive characteristic. We establish the rationality theorem for these local zeta functions and list explicitly all the candidate poles. Our result generalizes the work of Le$\acute{o}$n-Cardenal, Ibadula and Segers and that of Yin and Hong.


title: Surrogate Modeling for Physical Systems with Preserved Properties and Adjustable Tradeoffs

abstract: Determining the proper level of details to develop and solve physical models is usually difficult when one encounters new engineering problems. Such difficulty comes from how to balance the time (simulation cost) and accuracy for the physical model simulation afterwards. We propose a framework for automatic development of a family of surrogate models of physical systems that provide flexible cost-accuracy tradeoffs to assist making such determinations. We present both a model-based and a data-driven strategy to generate surrogate models. The former starts from a high-fidelity model generated from first principles and applies a bottom-up model order reduction (MOR) that preserves stability and convergence while providing a priori error bounds, although the resulting reduced-order model may lose its interpretability. The latter generates interpretable surrogate models by fitting artificial constitutive relations to a presupposed topological structure using experimental or simulation data. For the latter, we use Tonti diagrams to systematically produce differential equations from the assumed topological structure using algebraic topological semantics that are common to various lumped-parameter models (LPM). The parameter for the constitutive relations are estimated using standard system identification algorithms. Our framework is compatible with various spatial discretization schemes for distributed parameter models (DPM), and can supports solving engineering problems in different domains of physics.


title: Low-Rank and Row-Sparse Decomposition for Joint DOA Estimation and Distorted Sensor Detection

abstract: Distorted sensors could occur randomly and may lead to the breakdown of a sensor array system. We consider an array model within which a small number of sensors are distorted by unknown sensor gain and phase errors. With such an array model, the problem of joint direction-of-arrival (DOA) estimation and distorted sensor detection is investigated and the problem is formulated under the framework of low-rank and row-sparse decomposition. We derive an iteratively reweighted least squares (IRLS) algorithm to solve the resulting problem in both noiseless and noisy cases. The convergence property of the IRLS algorithm is analyzed by means of the monotonicity and boundedness of the objective function. Extensive simulations are conducted regarding parameter selection, convergence speed, computational complexity, and performances of DOA estimation as well as distorted sensor detection. Even though the IRLS algorithm is slightly worse than the alternating direction method of multipliers in detecting the distorted sensors, the results show that our approach outperforms several state-of-the-art techniques in terms of convergence speed, computational cost, and DOA estimation performance.


title: Federated Reinforcement Learning for Collective Navigation of Robotic Swarms

abstract: The recent advancement of Deep Reinforcement Learning (DRL) contributed to robotics by allowing automatic controller design. The automatic controller design is a crucial approach for designing swarm robotic systems, which require more complex controllers than a single robot system to lead a desired collective behaviour. Although the DRL-based controller design method showed its effectiveness, the reliance on the central training server is a critical problem in real-world environments where robot-server communication is unstable or limited. We propose a novel Federated Learning (FL) based DRL training strategy (FLDDPG) for use in swarm robotic applications. Through the comparison with baseline strategies under a limited communication bandwidth scenario, it is shown that the FLDDPG method resulted in higher robustness and generalisation ability into a different environment and real robots, while the baseline strategies suffer from the limitation of communication bandwidth. This result suggests that the proposed method can benefit swarm robotic systems operating in environments with limited communication bandwidth, e.g., in high-radiation, underwater, or subterranean environments.


title: Pop Quiz! Can a Large Language Model Help With Reverse Engineering?

abstract: Large language models (such as OpenAI's Codex) have demonstrated impressive zero-shot multi-task capabilities in the software domain, including code explanation. In this work, we examine if this ability can be used to help with reverse engineering. Specifically, we investigate prompting Codex to identify the purpose, capabilities, and important variable names or values from code, even when the code is produced through decompilation. Alongside an examination of the model's responses in answering open-ended questions, we devise a true/false quiz framework to characterize the performance of the language model. We present an extensive quantitative analysis of the measured performance of the language model on a set of program purpose identification and information extraction tasks: of the 136,260 questions we posed, it answered 72,754 correctly. A key takeaway is that while promising, LLMs are not yet ready for zero-shot reverse engineering.


title: Improved Integrality Gap in Max-Min Allocation: or Topology at the North Pole

abstract: In the max-min allocation problem a set $P$ of players are to be allocated disjoint subsets of a set $R$ of indivisible resources, such that the minimum utility among all players is maximized. We study the restricted variant, also known as the Santa Claus problem, where each resource has an intrinsic positive value, and each player covets a subset of the resources. Bez\'akov\'a and Dani showed that this problem is NP-hard to approximate within a factor less than $2$, consequently a great deal of work has focused on approximate solutions. The principal approach for obtaining approximation algorithms has been via the Configuration LP (CLP) of Bansal and Sviridenko. Accordingly, there has been much interest in bounding the integrality gap of this CLP. The existing algorithms and integrality gap estimations are all based one way or another on the combinatorial augmenting tree argument of Haxell for finding perfect matchings in certain hypergraphs. Our main innovation in this paper is to introduce the use of topological methods for the restricted max-min allocation problem, to replace the combinatorial argument. This approach yields substantial improvements in the integrality gap of the CLP. In particular we improve the previously best known bound of $3.808$ to $3.534$. We also study the $(1,\varepsilon)$-restricted version, in which resources can take only two values, and improve the integrality gap in most cases.


title: An Adiabatic Capacitive Artificial Neuron with RRAM-based Threshold Detection for Energy-Efficient Neuromorphic Computing

abstract: In the quest for low power, bio-inspired computation both memristive and memcapacitive-based Artificial Neural Networks (ANN) have been the subjects of increasing focus for hardware implementation of neuromorphic computing. One step further, regenerative capacitive neural networks, which call for the use of adiabatic computing, offer a tantalising route towards even lower energy consumption, especially when combined with `memimpedace' elements. Here, we present an artificial neuron featuring adiabatic synapse capacitors to produce membrane potentials for the somas of neurons; the latter implemented via dynamic latched comparators augmented with Resistive Random-Access Memory (RRAM) devices. Our initial 4-bit adiabatic capacitive neuron proof-of-concept example shows 90% synaptic energy saving. At 4 synapses/soma we already witness an overall 35% energy reduction. Furthermore, the impact of process and temperature on the 4-bit adiabatic synapse shows a maximum energy variation of 30% at 100 degree Celsius across the corners without any functionality loss. Finally, the efficacy of our adiabatic approach to ANN is tested for 512 & 1024 synapse/neuron for worst and best case synapse loading conditions and variable equalising capacitance's quantifying the expected trade-off between equalisation capacitance and range of optimal power-clock frequencies vs. loading (i.e. the percentage of active synapses).


title: Relative Position Prediction as Pre-training for Text Encoders

abstract: Meaning is defined by the company it keeps. However, company is two-fold: It's based on the identity of tokens and also on their position (topology). We argue that a position-centric perspective is more general and useful. The classic MLM and CLM objectives in NLP are easily phrased as position predictions over the whole vocabulary. Adapting the relative position encoding paradigm in NLP to create relative labels for self-supervised learning, we seek to show superior pre-training judged by performance on downstream tasks.


title: Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

abstract: The high computational scaling with the number of correlated electrons and the size of the basis set is a bottleneck which limits applications of coupled cluster (CC) algorithms. This is particularly so for calculations based on 2- or 4-component relativistic Hamiltonians, which often employ uncontracted basis sets and lead to large virtual molecular orbital (VMO) spaces. This problem may be alleviated by employing a more compact set of virtual spinors than those provided by the canonical Hartree-Fock (HF) set, such as those based on natural orbitals (NOs). In this paper we describe the implementation of a module for generating NOs for correlated wavefunctions, and in particular MP2 frozen natural orbitals (MP2FNOs), as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [J. Pototschnig et. al., J. Chem. Theory Comput. 17, 5509, 2021]. Our implementation is capable of manipulating both complex and quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, so that the code also makes it possible to generate CCSD natural orbitals in AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties at CCSD level such as the electric field gradients at the nuclei (EFGs), electric dipole and quadrupole moments for hydrogen halides HX (X=F-Ts), and parity-violating energy differences (PV) for H$_2$Y$_2$ (Y=O-Se), we find that MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the periodic table and that, with VMO spaces truncated to around half the size of the full spaces ones, it is possible to obtain reliable estimates for both energies and all molecular properties considered.


title: Improving Screening Processes via Calibrated Subset Selection

abstract: Many selection processes such as finding patients qualifying for a medical trial or retrieval pipelines in search engines consist of multiple stages, where an initial screening stage focuses the resources on shortlisting the most promising candidates. In this paper, we investigate what guarantees a screening classifier can provide, independently of whether it is constructed manually or trained. We find that current solutions do not enjoy distribution-free theoretical guarantees -- we show that, in general, even for a perfectly calibrated classifier, there always exist specific pools of candidates for which its shortlist is suboptimal. Then, we develop a distribution-free screening algorithm -- called Calibrated Subset Selection (CSS) -- that, given any classifier and some amount of calibration data, finds near-optimal shortlists of candidates that contain a desired number of qualified candidates in expectation. Moreover, we show that a variant of CSS that calibrates a given classifier multiple times across specific groups can create shortlists with provable diversity guarantees. Experiments on US Census survey data validate our theoretical results and show that the shortlists provided by our algorithm are superior to those provided by several competitive baselines.


title: Canceling and inverting normal and anomalous group-velocity dispersion using space-time wave packets

abstract: Angular dispersion can counterbalance normal group-velocity dispersion (GVD) that increases the wave-vector length in a dispersive medium. By tilting the wave vector, angular dispersion reduces the axial wave number in this case to match the pre-GVD value. By the same token, however, angular dispersion fails to counterbalance anomalous GVD, which in contrast reduces the wave-vector length. Consequently, GVD-cancellation via angular dispersion has not been demonstrated to date in the anomalous dispersion regime. We synthesize here structured femtosecond pulsed beams known as `space-time' wave packets designed to realize dispersion-cancellation symmetrically in either the normal- or anomalous-GVD regimes by virtue of non-differentiable angular dispersion inculcated into the pulsed field. Furthermore, we also verify GVD-inversion: reversing the GVD sign experienced by the field with respect to that dictated by the chromatic dispersion of the medium itself.


title: Phonon scattering and vibrational localization in 2D embedded nanoparticle composites

abstract: In this work, a Landauer approach enabled by the Frequency Domain Perfectly Matched Layer Method (FDPML) is used to study phonon transport in a series of large 2D domains with randomly embedded nanoparticles over a wide range of nanoparticle loadings and wavelengths. The effect of nanoparticle packing density on the mean free path and localization length is characterized. We observe that in the Mie scattering regime, the independent scattering approximation is valid up to volume fractions exceeding 10% and often higher depending on scattering parameter, indicating the mean free path can usually be calculated much less expensively using the number density and the scattering cross-section of a single scatterer. We also study localization lengths and their dependence on particle loading. In the case of heavy particles in a lighter matrix, we have been able to observe localization only at volume fractions >30% using the Landauer approach and only for high frequency modes, exceeding the vibration frequencies of the embedded nanoparticles. Using modal analysis we show that localization in nanoparticle laden materials is primarily due to energetic confinement rather than Anderson localization. Subsequently, we show that by using light particles in a heavy matrix the fraction of confined modes can be substantially increased.


title: Fusions of the generalized Hamming scheme on a strongly-regular graph

abstract: In this paper we show that for any fusion $\mathcal{B}$ of an association scheme $\mathcal{A}$, the generalized Hamming scheme $H(n,\mathcal{B})$ is a nontrivial fusion of $H(n,\mathcal{A})$. We analyze the case where $\mathcal{A}$ is the association scheme on a strongly-regular graph, and determine the parameters of all strongly-regular graphs for which the generalized Hamming scheme, $H(2,\mathcal{A})$, has extra fusions, in addition to the one arising from the trivial fusion of $\mathcal{A}$.


title: Investigating the detectability of hydrocarbons in exoplanet atmospheres with JWST

abstract: We investigate at what abundances various hydrocarbon molecules (e.g. acetylene (C$_2$H$_2$), ethylene (C$_2$H$_4$), and methane (CH$_4$)) become detectable when observing the atmospheres of various planets using JWST. We focused on atmospheric models based on the parameters of a small sample of planets: HD 189733b, HD 209458b, HD 97658b, and Kepler-30c. We computed model transmission spectra using the Bayesian retrieval package ARCiS. We simulated observed spectra using the PandExo package. We subsequently ran retrievals on these spectra to determine whether the molecular abundances can be accurately retrieved from these simulated spectra. We find that generally we can detect and retrieve abundances of the hydrocarbon species as long as they have a VMR above approximately 1x10$^{-7}$-1x10$^{-6}$, at least for the brighter targets. There are variations based on planet type and instrument(s) used, and these limits will likely change depending on the abundance of other strong absorbers. We also find scenarios where the presence of one hydrocarbon is confused with another; this is often improved when two instruments are combined. C$_2$H$_2$, CH$_4$, and C$_2$H$_4$ will all be detectable with JWST, provided they are present in high enough abundances, and that the optimal instruments are chosen. A combination of two instruments, either NIRSpec G395M and MIRI LRS, or NIRCam F322W2 and MIRI LRS, is best for observing these species in bright exoplanet systems, with NIRSpec G395M and MIRI LRS the best option for the HD 189733b-like atmosphere with clouds included. The use of NIRSpec Prism is tentatively found to be best for fainter targets, potentially in combination with the MIRI LRS slit mode, although the target we test is too faint to draw any strong conclusions. Instrument sensitivity, noise, and wavelength range are all thought to play a role in being able to distinguish spectral features.


title: Defect Controlled Ferromagnetic Ordering in Au Implanted TiSe$_2$ Nanocrystals

abstract: Layered transition metal dichalcogenides (TMDs) are attracting increasing attention because they exhibit unconventional magnetic properties due to crystal imperfections in their usually non-magnetic 2D structure. This work aims to investigate the magnetic response of self-engineered Se deficient TiSe$_2$ thin films, synthesized using chemical vapour deposition. We demonstrate tunability of the ferromagnetic order with the introduction of Au atoms using low energy Au ion implantation, which works as a controlling knob to vary the stoichiometry of Se in TiSe$_{2-x}$. The corresponding isothermal field-magnetization curves fit well with a modified Brillouin J function with J value of 1.5 for Ti$^{3+}$, and 4 for Au$^{3+}$, accounting for the diamagnetism that arises from Au implantation. We propose a qualitative model for the experimentally observed magnetization as a function of ion fluence, corroborated with high-resolution transmission electron microscopy. Depending on the Au nanoparticle size in the implanted samples, magnetization saturates faster at a much lower applied magnetic field than the pristine sample. Our findings hold potential to expand the range of 2D ferromagnetic materials for spintronics and magnetic sensing applications.


title: Analogies and Feature Attributions for Model Agnostic Explanation of Similarity Learners

abstract: Post-hoc explanations for black box models have been studied extensively in classification and regression settings. However, explanations for models that output similarity between two inputs have received comparatively lesser attention. In this paper, we provide model agnostic local explanations for similarity learners applicable to tabular and text data. We first propose a method that provides feature attributions to explain the similarity between a pair of inputs as determined by a black box similarity learner. We then propose analogies as a new form of explanation in machine learning. Here the goal is to identify diverse analogous pairs of examples that share the same level of similarity as the input pair and provide insight into (latent) factors underlying the model's prediction. The selection of analogies can optionally leverage feature attributions, thus connecting the two forms of explanation while still maintaining complementarity. We prove that our analogy objective function is submodular, making the search for good-quality analogies efficient. We apply the proposed approaches to explain similarities between sentences as predicted by a state-of-the-art sentence encoder, and between patients in a healthcare utilization application. Efficacy is measured through quantitative evaluations, a careful user study, and examples of explanations.


title: Softly higher-derivative massive gravity

abstract: In this work, we study the higher-derivative massive gravity in $D$-dimensional space-time with the mass term arisen due to a 2-brane embedded in a $4D$ Minkowski spacetime. We consider the effect of a resonance mass term from the DGP braneworld model for the higher-derivative massive gravity. Moreover, we add the gravitational Chern-Simons and Ricci-Cotton terms to this model and evaluate the graviton propagator using a basis of Barnes-Rivers spin projectors. Using the obtained propagator, we proceed with discussing the consistency of the model, writing the dispersion relations, and analyzing causality and unitarity.


title: The slurk Interaction Server Framework: Better Data for Better Dialog Models

abstract: This paper presents the slurk software, a lightweight interaction server for setting up dialog data collections and running experiments. Slurk enables a multitude of settings including text-based, speech and video interaction between two or more humans or humans and bots, and a multimodal display area for presenting shared or private interactive context. The software is implemented in Python with an HTML and JS frontend that can easily be adapted to individual needs. It also provides a setup for pairing participants on common crowdworking platforms such as Amazon Mechanical Turk and some example bot scripts for common interaction scenarios.


title: Antiferromagnetic parametric resonance driven by voltage-controlled magnetic anisotropy

abstract: Voltage controlled magnetic anisotropy (VCMA) is a low-energy alternative to manipulate the ferromagnetic state, which has been recently considered also in antiferromagnets (AFMs). Here, we theoretically demonstrate that VCMA can be used to excite linear and parametric resonant modes in easy-axis AFMs with perpendicular anisotropy, thus opening the way for an efficient electrical control of the Neel vector, and for detection of high-frequency dynamics. Our work leads to two key results: (i) VCMA parametric pumping experiences the so-called exchange enhancement of the coupling efficiency and, thus, is 1-2 orders of magnitude more efficient than microwave magnetic fields or spin-orbit-torques, and (ii) it also allows for zero-field parametric resonance, which cannot be achieved by other parametric pumping mechanisms in AFMs with out-of-plane easy axis. Therefore, we demonstrate that the VCMA parametric pumping is the most promising method for coherent excitation and manipulation of AFM order in perpendicular easy-axis AFMs.


title: Error Correction in ASR using Sequence-to-Sequence Models

abstract: Post-editing in Automatic Speech Recognition (ASR) entails automatically correcting common and systematic errors produced by the ASR system. The outputs of an ASR system are largely prone to phonetic and spelling errors. In this paper, we propose to use a powerful pre-trained sequence-to-sequence model, BART, further adaptively trained to serve as a denoising model, to correct errors of such types. The adaptive training is performed on an augmented dataset obtained by synthetically inducing errors as well as by incorporating actual errors from an existing ASR system. We also propose a simple approach to rescore the outputs using word level alignments. Experimental results on accented speech data demonstrate that our strategy effectively rectifies a significant number of ASR errors and produces improved WER results when compared against a competitive baseline. We also highlight a negative result obtained on the related grammatical error correction task in Hindi language showing the limitation in capturing wider context by our proposed model.


title: GADGET: Online Resource Optimization for Scheduling Ring-All-Reduce Learning Jobs

abstract: Fueled by advances in distributed deep learning (DDL), recent years have witnessed a rapidly growing demand for resource-intensive distributed/parallel computing to process DDL computing jobs. To resolve network communication bottleneck and load balancing issues in distributed computing, the so-called ``ring-all-reduce'' decentralized architecture has been increasingly adopted to remove the need for dedicated parameter servers. To date, however, there remains a lack of theoretical understanding on how to design resource optimization algorithms for efficiently scheduling ring-all-reduce DDL jobs in computing clusters. This motivates us to fill this gap by proposing a series of new resource scheduling designs for ring-all-reduce DDL jobs. Our contributions in this paper are three-fold: i) We propose a new resource scheduling analytical model for ring-all-reduce deep learning, which covers a wide range of objectives in DDL performance optimization (e.g., excessive training avoidance, energy efficiency, fairness); ii) Based on the proposed performance analytical model, we develop an efficient resource scheduling algorithm called GADGET (greedy ring-all-reduce distributed graph embedding technique), which enjoys a provable strong performance guarantee; iii) We conduct extensive trace-driven experiments to demonstrate the effectiveness of the GADGET approach and its superiority over the state of the art.


title: L3Cube-MahaCorpus and MahaBERT: Marathi Monolingual Corpus, Marathi BERT Language Models, and Resources

abstract: We present L3Cube-MahaCorpus a Marathi monolingual data set scraped from different internet sources. We expand the existing Marathi monolingual corpus with 24.8M sentences and 289M tokens. We further present, MahaBERT, MahaAlBERT, and MahaRoBerta all BERT-based masked language models, and MahaFT, the fast text word embeddings both trained on full Marathi corpus with 752M tokens. We show the effectiveness of these resources on downstream Marathi sentiment analysis, text classification, and named entity recognition (NER) tasks. We also release MahaGPT, a generative Marathi GPT model trained on Marathi corpus. Marathi is a popular language in India but still lacks these resources. This work is a step forward in building open resources for the Marathi language. The data and models are available at https://github.com/l3cube-pune/MarathiNLP .


title: Saving Brian's Privacy: the Perils of Privacy Exposure through Reverse DNS

abstract: Given the importance of privacy, many Internet protocols are nowadays designed with privacy in mind (e.g., using TLS for confidentiality). Foreseeing all privacy issues at the time of protocol design is, however, challenging and may become near impossible when interaction out of protocol bounds occurs. One demonstrably not well understood interaction occurs when DHCP exchanges are accompanied by automated changes to the global DNS (e.g., to dynamically add hostnames for allocated IP addresses). As we will substantiate, this is a privacy risk: one may be able to infer device presence and network dynamics from virtually anywhere on the Internet -- and even identify and track individuals -- even if other mechanisms to limit tracking by outsiders (e.g., blocking pings) are in place. We present a first of its kind study into this risk. We identify networks that expose client identifiers in reverse DNS records and study the relation between the presence of clients and said records. Our results show a strong link: in 9 out of 10 cases, records linger for at most an hour, for a selection of academic, enterprise and ISP networks alike. We also demonstrate how client patterns and network dynamics can be learned, by tracking devices owned by persons named Brian over time, revealing shifts in work patterns caused by COVID-19 related work-from-home measures, and by determining a good time to stage a heist.


title: Strategies for variational quantum compiling of a zero-phase beamsplitter on the Xanadu X8 processor

abstract: In the context of variational compiling of a continuous-variable (CV) unitary operation, the architecture and parameter space of the Xanadu X8 processor constrain both the set of feasible compiling problems and the allowed cost functions. In this paper, we motivate a faithful cost function for variational compiling of a two-mode, continuous-variable beamsplitter gate with real matrix elements (i.e., a "zero-phase" beamsplitter) that complies with the constraints of the X8 processor. This cost function is then computed on the X8. Despite the noise in the processor, we find that the cost function exhibits optimum parameter resilience and, therefore, that this variational compiling problem is feasible on the X8. The intent of the paper is partly to report a proof-of-principle cost function calculation on near-term CV hardware, and partly to present methods that may be relevant for CV variational compiling problems in more complicated, large scale settings.


title: The Role of Second Law of Thermodynamics in Continuum Physics: A Muschik and Ehrentraut Theorem Revisited

abstract: Second law of thermodynamics imposes that in any thermodynamic process the entropy production must be nonnegative. In continuum physics such a requirement is fulfilled by postulating the constitutive equations which represent the material properties of the bodies in such a way that second law of thermodynamics is satisfied in arbitrary processes. Such an approach, first assumed in some pioneering papers by Coleman and Noll \cite{ColNol} and Coleman and Mizel \cite{ColMiz}, in practice regards second law of thermodynamics as a restriction on the constitutive equations, which must guarantee that any solution of the balance laws satisfies also the entropy inequality. As observed by Muschik and Ehrentraut \cite{MusEhr}, this is a useful operative assumption, but not a consequence of general physical laws. Indeed, a different point of view, which regards second law of thermodynamics as a restriction on the thermodynamic processes, i.e., on the solutions of the system of balance laws, is possible. This is tantamount to assume that there are solutions of the balance laws which satisfy the entropy inequality, and solutions which do not satisfy it. In order to decide what is the correct approach, Muschik and Erhentraut postulated an amendment to the second law, which makes explicit the evident but rather hidden assumption that in any point of the body the entropy production is zero if, and only if, this point is thermodynamic equilibrium. Then they proved that, given the amendment, second law of thermodynamics is necessarily a restriction on the constitutive equations and not on the thermodynamic processes. In the present paper we revisit their proof, lighting up some geometric aspects which were hidden in Ref. \cite{MusEhr}. Moreover, we propose an alternative formulation of second law of thermodynamics which incorporates the amendment.


title: A Recommender System Based on a Double Feature Allocation Model

abstract: A collaborative filtering recommender system predicts user preferences by discovering common features among users and items. We implement such inference using a Bayesian double feature allocation model, that is, a model for random pairs of subsets. We use an Indian buffet process (IBP) to link users and items to features. Here a feature is a subset of users and a matching subset of items. By training feature-specific rating effects, we predict ratings. We use MovieLens Data to demonstrate posterior inference in the model and prediction of user preferences for unseen items compared to items they have previously rated. Part of the implementation is a novel semi-consensus Monte Carlo method to accomodate large numbers of users and items, as is typical for related applications. The proposed approach implements parallel posterior sampling in multiple shards of users while sharing item-related global parameters across shards.


title: Unprecedented Daylight Display of Kreutz Sungrazers in AD 363?

abstract: In the context of the recently proposed contact-binary model (Sekanina 2021), I investigate the circumstances of the first perihelion passage of the Kreutz sungrazers in orbits with barycentric periods near 735 yr, following the initial near-aphelion splitting of the presumed progenitor, Aristotle's comet of 372 BC. Given favorable conditions at this breakup and at episodes of secondary fragmentation in its aftermath, the fragments should have arrived at their first perihelion nearly simultaneously, reminiscent of the anticipated outcome for the two-superfragment model's perihelion return of AD 356 (Sekanina & Chodas 2004). The relevant case of a swarm of Kreutz sungrazers is examined to appraise possible scientific ramifications of the brief remark by Ammianus Marcellinus, a Roman historian, that "in broad daylight comets were seen" in late AD 363, only seven years later. The tested scenario, which does not contradict Ammianus' narrative and is consistent with the contact-binary model, involves a set of ten sungrazers visible in the daytime, all reaching perihelion over a period of 4.6 days. As part of this work, I comment on the role of the rapidly developing, brilliant post-perihelion tail; revise the apparent magnitude typical for the first and last naked-eye sightings; compare the visibility conditions in full daylight, in twilight, and at night; and, for the first time, present circumstantial evidence that favors comet X/1106 C1 as the parent to C/1843 D1 rather than to C/1882 R1 and C/1965 S1.


title: Axisymmetric gravity currents in anisotropic porous media

abstract: We explore the motion of an axisymmetric gravity current in an anisotropic porous medium in which the horizontal permeability is larger than the vertical permeability. It is well known that the classical axisymmetric gravity current supplied by a constant point source of fluid has an unphysical singularity near the origin. We address this by considering a pressure-dominated region near the origin which allows for vertical flow from the source, such that the current remains of finite depth, whilst beyond this region the flow is gravity-dominated. At early times the inner pressure-driven region controls the spreading of the current, but at late times the inner region occupies a progressively smaller fraction of the current such that the radius increases as $\sim t^{3/7}$, while the depth near the origin increases approximately as $\sim t^{1/7}$. The presence of anisotropy highlights this phenomenon, since the vertical permeability maintains an effect on the flow at late times through the pressure-driven flow near the origin. Using these results we provide some quantitative insights into the dominant dynamics which control CO$_2$ migration through permeable aquifers, as occurs in the context of carbon capture and storage.


title: General linear stability properties of monoclinal shallow waves

abstract: We analyze the linear stability of monoclinal traveling waves on a constant incline, which connect uniform flowing regions of differing depths. The classical shallow-water equations are employed, subject to a general resistive drag term. This approach incorporates many flow rheologies into a single setting and enables us to investigate the features that set different systems apart. We derive simple formulae for the onset of linear instability, the corresponding linear growth rates and related properties including the existence of monoclinal waves, development of shocks and whether instability is initially triggered up- or downstream of the wave front. Also included within our framework is the presence of shear in the flow velocity profile, which is often neglected in depth-averaged studies. We find that it can significantly modify the threshold for instability. Constant corrections to the governing equations to account for sheared profiles via a 'momentum shape factor' act to stabilize traveling waves. More general correction terms are found to have a nontrivial and potentially important quantitative effect on the properties explored. Finally, we have investigated the spatial properties of the dominant (fastest growing) linear modes. We derive equations for their amplitude and frequency and find that both features can become severely amplified near the front of the traveling wave. For flood waves that propagate into a dry downstream region, this amplification is unbounded in the limit of high disturbance frequency. We show that the rate of divergence is a function of the spatial dependence of the wave depth profile at the front, which may be determined straightforwardly from the drag law.


title: Exclusive nonleptonic Bc-meson decays to S-wave charmonium states

abstract: We study the exclusive two-body nonleptonic $B_c\to X_{c\bar{c}}M$ decays, where $X_{c\bar{c}}$ is either a ground(1S) or a radially excited (2S or 3S) charmonium and $M$ is a pseudoscalar $(P)$ or a vector $(V)$ meson. We consider here three categories of decays: $B_c \to PP, PV, VP$ decays within the framework of relativistic independent quark(RIQ) model based on a flavor-independent interaction potential in scalar-vector harmonic form. Using the factorization approximation, we calculate the weak form factors from the overlapping integrals of meson wave functions obtained in the RIQ model and predict the branching fractions for a set of exclusive nonleptonic $B_c$-decays in reasonable agreement with other model predictions. Some of the decays of interest are found to have branching fractions $\sim (10^{-3}-10^{-4})$ within the detection ability of the current experiments and can be precisely measured at LHCb in near future. In the wake of the recent measurement of $B_c\to J/\psi\pi(K)$, $B_c\to J/\psi\pi(D_s)$, $B_c\to \pi(J/\psi,\psi(2S))$ and $B_c\to J/\psi(\pi,\mu \nu)$ reported by the LHCb Collaborations, we predict the ratios: ${\cal R}_{K/\pi}$, ${\cal R}_{D_s/\pi}$ and ${\cal R}_{\psi(2S)/{J/\psi}}$ in broad agreement with the LHCb data though our predicted ratio ${\cal R}_{\pi /{\mu\nu}}$ is found to be underestimated. The results indicate that the present approach works well in the description of exclusive nonleptonic $B_c$-decays within the framework of the RIQ model.


title: Winding number and homotopy for quaternionic curves

abstract: Following a recent approach to quaternionic curves, we defined the quaternionic polar angle that enabled us to define global properties of quaternionic curves, namely the winding number and the homotopy concept. The results admit various applications, including further analogies to plane curves, and physical applications.


title: Unified Scaling Laws for Routed Language Models

abstract: The performance of a language model has been shown to be effectively modeled as a power-law in its parameter count. Here we study the scaling behaviors of Routing Networks: architectures that conditionally use only a subset of their parameters while processing an input. For these models, parameter count and computational requirement form two independent axes along which an increase leads to better performance. In this work we derive and justify scaling laws defined on these two variables which generalize those known for standard language models and describe the performance of a wide range of routing architectures trained via three different techniques. Afterwards we provide two applications of these laws: first deriving an Effective Parameter Count along which all models scale at the same rate, and then using the scaling coefficients to give a quantitative comparison of the three routing techniques considered. Our analysis derives from an extensive evaluation of Routing Networks across five orders of magnitude of size, including models with hundreds of experts and hundreds of billions of parameters.


title: Kernel-Function Based Quantum Algorithms for Finite Temperature Quantum Simulation

abstract: Computing finite temperature properties of a quantum many-body system is key to describing a broad range of correlated quantum many-body physics from quantum chemistry and condensed matter to thermal quantum field theories. Quantum computing with rapid developments in recent years has a huge potential to impact the computation of quantum thermodynamics. To fulfill the potential impacts, it is crucial to design quantum algorithms that utilize the computation power of the quantum computing devices. Here we present a quantum kernel function expansion (QKFE) algorithm for solving thermodynamic properties of quantum many-body systems. In this quantum algorithm, the many-body density of states is approximated by a kernel-Fourier expansion, whose expansion moments are obtained by random state sampling and quantum interferometric measurements. As compared to its classical counterpart, namely the kernel polynomial method (KPM), QKFE has an exponential advantage in the cost of both time and memory. In computing low temperature properties, QKFE becomes inefficient, as similar to classical KPM. To resolve this difficulty, we further construct a thermal ensemble and approaches the low temperature regime step-by-step. For quantum Hamiltonians, whose ground states are preparable with polynomial quantum circuits, THEI has an overall polynomial complexity. We demonstrate its efficiency with applications to one and two-dimensional quantum spin models, and a fermionic lattice. With our analysis on the realization with digital and analogue quantum devices, we expect the quantum algorithm is accessible to current quantum technology.


title: Low regularity integrators via decorated trees

abstract: We introduce a general framework of low regularity integrators which allows us to approximate the time dynamics of a large class of equations, including parabolic and hyperbolic problems, as well as dispersive equations, up to arbitrary high order on general domains. The structure of the local error of the new schemes is driven by nested commutators which in general require (much) lower regularity assumptions than classical methods do. Our main idea lies in embedding the central oscillations of the nonlinear PDE into the numerical discretisation. The latter is achieved by a novel decorated tree formalism inspired by singular SPDEs with Regularity Structures and allows us to control the nonlinear interactions in the system up to arbitrary high order on the infinite dimensional (continuous) as well as finite dimensional (discrete) level.


title: Carrollian superconformal theories and super BMS

abstract: We investigate the supersymmetric versions of Bondi-Metzner-Sachs or, equivalently, conformal Carroll symmetry in boundary dimensions $d>3$, with applications of flat space holography in mind. We identify the contraction of the relativistic symmetry relevant for our purposes and construct a finite-dimensional Carrollian superconformal algebra (CSA) before proposing an infinite-dimensional lift. We provide the superspace formulation for $\mathcal{N}=1$ CSA and work towards an understanding of the representation theory of the algebra. We conclude with some aspects of $\mathcal{N}>1$ CSA.


title: Deviation from maximal entanglement for mid-spectrum eigenstates of local Hamiltonians

abstract: In a spin chain governed by a local Hamiltonian, we consider a microcanonical ensemble in the middle of the energy spectrum and a contiguous subsystem whose length is a constant fraction of the system size. We prove that if the bandwidth of the ensemble is greater than a certain constant, then the average entanglement entropy (between the subsystem and the rest of the system) of eigenstates in the ensemble deviates from the maximum entropy by at least a positive constant. This result highlights the difference between the entanglement entropy of mid-spectrum eigenstates of (chaotic) local Hamiltonians and that of random states. We also prove that the former deviates from the thermodynamic entropy at the same energy by at least a positive constant.


title: Evitable iterates of the consistency operator

abstract: Let's fix a reasonable subsystem $T$ of arithmetic; why are natural extensions of $T$ pre-well-ordered by consistency strength? In previous work, an approach to this question was proposed. The goal of this work was to classify the recursive functions that are monotone with respect to the Lindenabum algebra of $T$. According to an optimistic conjecture, roughly, every such function must be equivalent to an iterate $\mathsf{Con}_T^\alpha$ of the consistency operator in the limit. In previous work the author established the first case of this optimistic conjecture; roughly, every recursive monotone function is either as weak as the identity operator in the limit or as strong as $\mathsf{Con}_T$ in the limit. Yet in this note we prove that this optimistic conjecture fails already at the next step; there are recursive monotone functions that are neither as weak as $\mathsf{Con}_T$ in the limit nor as strong as $\mathsf{Con}_T^2$ in the limit. In fact, for every $\alpha$, we produce a function that is cofinally as strong as $\mathsf{Con}^\alpha_T$ yet cofinally as weak as $\mathsf{Con}_T$.


title: Surface acoustic waves inside polystyrene microparticles through photoacoustic microscopy

abstract: We demonstrate surface acoustic waves inside polystyrene microspheres of different sizes experimentally through photoacoustic microscopy and validate the experimental result with simulation. A novel method for sample preparation of a lifted sample is also presented where the microparticles are suspended in agarose above the surface of the petridish. Another objective of this study was to investigate the results with different laser focus hitting on the microparticles and their impact on photoacoustic images and photoacoustic signals. An absorbing microsphere is excited with a pulsed laser of wavelength 532 nm and the photoacoustic signal is detected using a 40 MHz transducer. On analyzing the photoacoustic signals from microspheres, we find the signature of surface acoustic waves.