title: Atomistic mechanisms of binary alloy surface segregation from nanoseconds to seconds using accelerated dynamics

abstract: Although the equilibrium composition of many alloy surfaces is well understood, the rate of transient surface segregation during annealing is not known, despite its crucial effect on alloy corrosion and catalytic reactions occurring on overlapping timescales. In this work, CuNi bimetallic alloys representing (100) surface facets are annealed in vacuum using atomistic simulations to observe the effect of vacancy diffusion on surface separation. We employ multi-timescale methods to sample the early transient, intermediate, and equilibrium states of slab surfaces during the separation process, including standard MD as well as three methods to perform atomistic, long-time dynamics: parallel trajectory splicing (ParSplice), adaptive kinetic Monte Carlo (AKMC), and kinetic Monte Carlo (KMC). From nanosecond (ns) to second timescales, our multiscale computational methodology can observe rare stochastic events not typically seen with standard MD, closing the gap between computational and experimental timescales for surface segregation. Rapid diffusion of a vacancy to the slab is resolved by all four methods in tens of ns. Stochastic re-entry of vacancies into the subsurface, however, is only seen on the microsecond timescale in the two KMC methods. Kinetic vacancy trapping on the surface and its effect on the segregation rate are discussed. The equilibrium composition profile of CuNi after segregation during annealing is estimated to occur on a timescale of seconds as determined by KMC, a result directly comparable to nanoscale experiments.


title: Beyond $L_p$ clipping: Equalization-based Psychoacoustic Attacks against ASRs

abstract: Automatic Speech Recognition (ASR) systems convert speech into text and can be placed into two broad categories: traditional and fully end-to-end. Both types have been shown to be vulnerable to adversarial audio examples that sound benign to the human ear but force the ASR to produce malicious transcriptions. Of these attacks, only the "psychoacoustic" attacks can create examples with relatively imperceptible perturbations, as they leverage the knowledge of the human auditory system. Unfortunately, existing psychoacoustic attacks can only be applied against traditional models, and are obsolete against the newer, fully end-to-end ASRs. In this paper, we propose an equalization-based psychoacoustic attack that can exploit both traditional and fully end-to-end ASRs. We successfully demonstrate our attack against real-world ASRs that include DeepSpeech and Wav2Letter. Moreover, we employ a user study to verify that our method creates low audible distortion. Specifically, 80 of the 100 participants voted in favor of all our attack audio samples as less noisier than the existing state-of-the-art attack. Through this, we demonstrate both types of existing ASR pipelines can be exploited with minimum degradation to attack audio quality.


title: IRS-Aided Radar: Enhanced Target Parameter Estimation via Intelligent Reflecting Surfaces

abstract: The intelligent reflecting surface (IRS) technology has recently attracted a lot of interest in wireless communications research. An IRS consists of passive reflective elements capable of tuning the phase, amplitude, frequency and polarization of the impinging waveforms. Given such desirable properties, the wireless channel characteristics can be controlled and optimized for specific signal design and processing needs -- thus promising significant potential in radar applications. In this paper, we establish the theoretical foundations for introducing IRS into a radar system and study the potential to improve target parameter estimation. More specifically, we will investigate the deployment of IRS in cases where the line-of-sight (LOS) link is weak or blocked by obstructions. We demonstrate that the IRS can provide a virtual or non-line-of-sight (NLOS) link between the radar and target leading to an enhanced radar performance. The effectiveness of such an IRS-provided virtual link in estimating the moving target parameters is illustrated under both optimized and non-optimized IRS scenarios. Numerical simulations indicate that the IRS can enhance the target parameter estimation when the LOS link is weaker than $\sim 10^{-1}$ in relative strength in comparison with the NLOS link.


title: VAC-CNN: A Visual Analytics System for Comparative Studies of Deep Convolutional Neural Networks

abstract: The rapid development of Convolutional Neural Networks (CNNs) in recent years has triggered significant breakthroughs in many machine learning (ML) applications. The ability to understand and compare various CNN models available is thus essential. The conventional approach with visualizing each model's quantitative features, such as classification accuracy and computational complexity, is not sufficient for a deeper understanding and comparison of the behaviors of different models. Moreover, most of the existing tools for assessing CNN behaviors only support comparison between two models and lack the flexibility of customizing the analysis tasks according to user needs. This paper presents a visual analytics system, VAC-CNN (Visual Analytics for Comparing CNNs), that supports the in-depth inspection of a single CNN model as well as comparative studies of two or more models. The ability to compare a larger number of (e.g., tens of) models especially distinguishes our system from previous ones. With a carefully designed model visualization and explaining support, VAC-CNN facilitates a highly interactive workflow that promptly presents both quantitative and qualitative information at each analysis stage. We demonstrate VAC-CNN's effectiveness for assisting novice ML practitioners in evaluating and comparing multiple CNN models through two use cases and one preliminary evaluation study using the image classification tasks on the ImageNet dataset.


title: Second Timescale Photometry of the Very Fast Nova V1674 Her with Palomar Gattini-IR

abstract: We report second-timescale infrared photometry of the nova V1674 Her using Palomar Gattini-IR. These observations constitute the first infrared and highest temporal resolution data (resolution of ~ 0.84 s) of the nova reported to date. PGIR observed in this fast readout mode for more than an hour on three nights between 3 and 6 days after discovery. We searched for periodic variability using a Lomb-Scargle periodogram and did not detect anything down to a three sigma upper limit of 0.093 mag. This suggests that the periodic variability detected in the nova by Patterson et al. 2021 was lower by at least a factor of about 1.65 in the first week of the eruption.


title: Pediatric Otoscopy Video Screening with Shift Contrastive Anomaly Detection

abstract: Ear related concerns and symptoms represents the leading indication for seeking pediatric healthcare attention. Despite the high incidence of such encounters, the diagnostic process of commonly encountered disease of the middle and external presents significant challenge. Much of this challenge stems from the lack of cost effective diagnostic testing, which necessitating the presence or absence of ear pathology to be determined clinically. Research has however demonstrated considerable variation among clinicians in their ability to accurately diagnose and consequently manage ear pathology. With recent advances in computer vision and machine learning, there is an increasing interest in helping clinicians to accurately diagnose middle and external ear pathology with computer-aided systems. It has been shown that AI has the capacity to analyse a single clinical image captured during examination of the ear canal and eardrum from which it can determine the likelihood of a pathognomonic pattern for a specific diagnosis being present. The capture of such an image can however be challenging especially to inexperienced clinicians. To help mitigate this technical challenge we have developed and tested a method using video sequences. We present a two stage method that first, identifies valid frames by detecting and extracting ear drum patches from the video sequence, and second, performs the proposed shift contrastive anomaly detection to flag the otoscopy video sequences as normal or abnormal. Our method achieves an AUROC of 88.0% on the patient-level and also outperforms the average of a group of 25 clinicians in a comparative study, which is the largest of such published to date. We conclude that the presented method achieves a promising first step towards automated analysis of otoscopy video.


title: Lower bounds for the cyclicity of centers of quadratic three-dimensional systems

abstract: We consider quadratic three-dimensional differential systems having a Hopf singular point. We study the cyclicity when the singular point is a center on the center manifold using higher order developments of the Lyapunov constants. As a result, we make a chart of the cyclicity by establishing the lower bounds for several known systems in the literature, among them the Rossler, Lorenz and Moon-Rand systems. Moreover, we obtain an example of a jerk system for which is possible to bifurcate 12 limit-cycles from the center, which is a new lower bound for three-dimensional quadratic systems.


title: Symbolic Substitutions, Bratteli Diagrams and Operator Algebras

abstract: In this paper we look at symbolic substitutions and their relationship to Bratteli diagrams and their associated operator algebras. In particular, we consider the equivalence relation on substitutions induced by telescope equivalence of Bratteli diagrams. Such an equivalence preserves pure aperiodicity and primitivity but fails to preserve rank, order, and number of letters. In a similar manner, we consider the equivalence relation on substitutions induced by telescope equivalence of ordered Bratteli diagrams. This results in a finer equivalence but fails to provide a complete invariant. An application to Fibonacci-like substitutions is developed.


title: Large sieve estimate for multivariate polynomial moduli and applications

abstract: We prove large sieve inequalities with multivariate polynomial moduli and deduce a general Bombieri--Vinogradov type theorem for a class of polynomial moduli having a sufficient number of variables compared to its degree. This sharpens previous results of the first author in two aspects: the range of the moduli as well as the class of polynomials which can be handled. As a consequence, we deduce that there exist infinitely many primes $p$such that $p-1$ has a prime divisor of size $\gg p^{2/5+o(1)}$ that is the value of an incomplete norm form polynomial.


title: Do female physics students benefit from informal physics programs they facilitate?

abstract: Gender bias, reduced sense of belonging, and lower physics self-efficacy are among the challenges faced by female students who choose to study physics. Prior studies focusing on this underrepresented group have examined the experiences and impacts of formal educational settings, leaving the impact of informal physics programs as a relatively overlooked area. Existing research on the impact of informal physics programs indicates that student facilitators, who help run the programs, can experience positive impacts on their learning and sense of community beyond the formal setting of a classroom or laboratory. In this study we took a first step, narrowing our focus to explore the relation between facilitation of informal physics programs and female students' physics identity, persistence, mindset, and worldview. We analyzed survey responses (32) and interviews (11) collected from undergraduate and graduate female students at a large, land-grant university. Our results, based on self-reported data, showed a statistically significant shift in confidence of choice of major after facilitating informal physics programs. Analysis of interviews indicated a positive effect of facilitation of informal programs on female student interest and motivation with regards to the field of physics, improved performance and competence beliefs, and the development of characteristics indicative of a growth mindset. A semantic network analysis showed statistically significant interdependencies between positive persistence and constructs including growth mindset, important early undergraduate experiences, gender stereotype threat, external recognition, and confidence.


title: Active Learning for Deep Visual Tracking

abstract: Convolutional neural networks (CNNs) have been successfully applied to the single target tracking task in recent years. Generally, training a deep CNN model requires numerous labeled training samples, and the number and quality of these samples directly affect the representational capability of the trained model. However, this approach is restrictive in practice, because manually labeling such a large number of training samples is time-consuming and prohibitively expensive. In this paper, we propose an active learning method for deep visual tracking, which selects and annotates the unlabeled samples to train the deep CNNs model. Under the guidance of active learning, the tracker based on the trained deep CNNs model can achieve competitive tracking performance while reducing the labeling cost. More specifically, to ensure the diversity of selected samples, we propose an active learning method based on multi-frame collaboration to select those training samples that should be and need to be annotated. Meanwhile, considering the representativeness of these selected samples, we adopt a nearest neighbor discrimination method based on the average nearest neighbor distance to screen isolated samples and low-quality samples. Therefore, the training samples subset selected based on our method requires only a given budget to maintain the diversity and representativeness of the entire sample set. Furthermore, we adopt a Tversky loss to improve the bounding box estimation of our tracker, which can ensure that the tracker achieves more accurate target states. Extensive experimental results confirm that our active learning-based tracker (ALT) achieves competitive tracking accuracy and speed compared with state-of-the-art trackers on the seven most challenging evaluation benchmarks.


title: Color-Kinematics Relation from the Feynman Diagram Perspective

abstract: Feynman diagrams for gluon tree amplitudes are studied in the Feynman gauge and in any number of spacetime dimensions. The color-kinematics combinations $\Delta=n_s-n_t-n_u$ of numerators are explicitly calculated for $N=4,5,6$ gluons to see whether the color-kinematics relation $\Delta=0$ is satisfied. This is a tedious task because of the presence of four-gluon vertices, and the large number of Feynman diagrams, numerators, and $\Delta$ combinations involved, especially when $N=6$. For on-shell amplitudes, it is found that $\Delta=0$ for $N=4$, but $\Delta\not=0$ for $N=5$ and $N=6$ owing to the presence of the four-gluon vertex. However, a {\it local} generalized gauge transformation can bring about $\Delta=0$ for $N=5$, but not for $N=6$. This raises the question whether gluon amplitudes satisfying the color-kinematics relation contain non-local interactions.


title: SWAP Test for an Arbitrary Number of Quantum States

abstract: We develop a recursive algorithm to generalize the quantum SWAP test for an arbitrary number $m$ of quantum states requiring $O(m)$ controlled-swap (CSWAP) gates and $O(\log m)$ ancillary qubits. We construct a quantum circuit able to simultaneously measure overlaps of $m$ arbitrary pure states. Our construction relies on a pairing unitary that generates a superposition state where every pair of input states is labelled by a basis state formed by the ancillaries.


title: Covariate Balancing Methods for Randomized Controlled Trials Are Not Adversarially Robust

abstract: The first step towards investigating the effectiveness of a treatment via a randomized trial is to split the population into control and treatment groups then compare the average response of the treatment group receiving the treatment to the control group receiving the placebo. In order to ensure that the difference between the two groups is caused only by the treatment, it is crucial that the control and the treatment groups have similar statistics. Indeed, the validity and reliability of a trial are determined by the similarity of two groups' statistics. Covariate balancing methods increase the similarity between the distributions of the two groups' covariates. However, often in practice, there are not enough samples to accurately estimate the groups' covariate distributions. In this paper, we empirically show that covariate balancing with the Standardized Means Difference (SMD) covariate balancing measure, as well as Pocock's sequential treatment assignment method, are susceptible to worst-case treatment assignments. Worst-case treatment assignments are those admitted by the covariate balance measure, but result in highest possible ATE estimation errors. We developed an adversarial attack to find adversarial treatment assignment for any given trial. Then, we provide an index to measure how close the given trial is to the worst-case. To this end, we provide an optimization-based algorithm, namely Adversarial Treatment ASsignment in TREatment Effect Trials (ATASTREET), to find the adversarial treatment assignments.


title: Construction of Fuchsian Schottky group with conformal boundary at infinity

abstract: In this article, we have constructed an interesting type of generalized Schottky group, named as Fuchsian Schottky group of arbitrary finite rank, in the context of the classical Schottky group (i.e., Schottky curves which are Euclidean circles). After that, we initiated the construction of the generalized Fuchsian Schottky group of any finite rank by including orientation-reversing isometries of the hyperbolic plane as side-pairing transformations. Further, we have investigated the hyperbolic ends for any arbitrary finite rank Fuchsian Schottky groups from the point of view of the Euler characteristic in the hyperbolic surface. Finally, we have shown that the compact core of the conformally compact Riemann surface can be decomposed into non-tight pairs of pants by using suitable twist parameters with some fixed Bers' constant. The Fenchel-Nielsen coordinates for Teichm\"uller space corresponding to any finite rank Fuchsian Schottky groups are also obtained.


title: Memory visualization tool for training neural network

abstract: Software developed helps world a better place ranging from system software, open source, application software and so on. Software engineering does have neural network models applied to code suggestion, bug report summarizing and so on to demonstrate their effectiveness at a real SE task. Software and machine learning algorithms combine to make software give better solutions and understanding of environment. In software, there are both generalized applications which helps solve problems for entire world and also some specific applications which helps one particular community. To address the computational challenge in deep learning, many tools exploit hardware features such as multi-core CPUs and many-core GPUs to shorten the training time. Machine learning algorithms have a greater impact in the world but there is a considerable amount of memory utilization during the process. We propose a new tool for analysis of memory utilized for developing and training deep learning models. Our tool results in visual utilization of memory concurrently. Various parameters affecting the memory utilization are analysed while training. This tool helps in knowing better idea of processes or models which consumes more memory.


title: On the Second-order Convergence Properties of Random Search Methods

abstract: We study the theoretical convergence properties of random-search methods when optimizing non-convex objective functions without having access to derivatives. We prove that standard random-search methods that do not rely on second-order information converge to a second-order stationary point. However, they suffer from an exponential complexity in terms of the input dimension of the problem. In order to address this issue, we propose a novel variant of random search that exploits negative curvature by only relying on function evaluations. We prove that this approach converges to a second-order stationary point at a much faster rate than vanilla methods: namely, the complexity in terms of the number of function evaluations is only linear in the problem dimension. We test our algorithm empirically and find good agreements with our theoretical results.


title: Image Quality Assessment using Contrastive Learning

abstract: We consider the problem of obtaining image quality representations in a self-supervised manner. We use prediction of distortion type and degree as an auxiliary task to learn features from an unlabeled image dataset containing a mixture of synthetic and realistic distortions. We then train a deep Convolutional Neural Network (CNN) using a contrastive pairwise objective to solve the auxiliary problem. We refer to the proposed training framework and resulting deep IQA model as the CONTRastive Image QUality Evaluator (CONTRIQUE). During evaluation, the CNN weights are frozen and a linear regressor maps the learned representations to quality scores in a No-Reference (NR) setting. We show through extensive experiments that CONTRIQUE achieves competitive performance when compared to state-of-the-art NR image quality models, even without any additional fine-tuning of the CNN backbone. The learned representations are highly robust and generalize well across images afflicted by either synthetic or authentic distortions. Our results suggest that powerful quality representations with perceptual relevance can be obtained without requiring large labeled subjective image quality datasets. The implementations used in this paper are available at \url{https://github.com/pavancm/CONTRIQUE}.


title: Reconstruction and Stability Analysis of Some Cosmological Bouncing Solutions in $F(\mathcal{R},T)$ Theory

abstract: The present article investigates the possibility of reconstruction of the generic function in $F(\mathcal{R},T)$ gravitational theory by considering some well-known cosmological bouncing models namely exponential evaluation, oscillatory, power law and matter bounce model, where $\mathcal{R}$ and $T$ are Ricci scalar and trace of energy-momentum tensor, respectively. Due to the complexity of dynamical field equations, we propose some ansatz forms of function $F(\mathcal{R},T)$ in perspective models and examine that which type of Lagrangian is capable to reproduce bouncing solution via analytical expression. It is seen that for some cases of exponential, oscillatory and matter bounce models, it is possible to get analytical solution while in other cases, it is not possible to achieve exact solutions so only complementary solutions can be discussed. However, for power law model, all forms of generic function can be reconstructed analytically. Further, we analyze the energy conditions and stability of these reconstructed cosmological bouncing models which have analytical forms. It is found that these models are stable for linear forms of Lagrangian only but the reconstructed solutions for power law are unstable for some non-linear forms of Lagrangian.


title: eXtreme19: when art and science make the front page

abstract: Back in mid 2018, we were organizing eXtreme19, a conference on astro-particle physics held in Padova on the topic of extremely energetic emission from galaxies. For the preparation of the graphical material in support of the conference we seeked for a collaboration with talented art students. To this purpose, we joined the Italian programme PCTO (Percorsi per le Competenze Trasversali e per l Orientamento) of high school student stages in job centers. Emily, Beatrice, and Chiara from the high school Liceo Artistico Statale Amedeo Modigliani in Padova accepted our invitation and started a 6-month stage at the Padova University in close contact with us. The challenge was to interbreed our scientific description of a relativistic jet of a powerful galaxy and their artistic assimilation and subsequent representation of it. During this period, they elaborated outstanding and innovative graphical material used for the webpage as well as the conference poster. The quality of the graphics was indeed excellent: one of their drawings became the cover of the February 2020 issue of the prestigious Nature Astronomy journal.


title: Facial Recognition in Collaborative Learning Videos

abstract: Face recognition in collaborative learning videos presents many challenges. In collaborative learning videos, students sit around a typical table at different positions to the recording camera, come and go, move around, get partially or fully occluded. Furthermore, the videos tend to be very long, requiring the development of fast and accurate methods. We develop a dynamic system of recognizing participants in collaborative learning systems. We address occlusion and recognition failures by using past information about the face detection history. We address the need for detecting faces from different poses and the need for speed by associating each participant with a collection of prototype faces computed through sampling or K-means clustering. Our results show that the proposed system is proven to be very fast and accurate. We also compare our system against a baseline system that uses InsightFace [2] and the original training video segments. We achieved an average accuracy of 86.2% compared to 70.8% for the baseline system. On average, our recognition rate was 28.1 times faster than the baseline system.


title: $B$ Flavour Anomalies: 2021 Theoretical Status Report

abstract: At the present time, there are discrepancies with the predictions of the SM in several observables involving $b \to s \ell^+ \ell^-$ and $b \to c \ell^- {\bar\nu}_\ell$ decays. These are the $B$ flavour anomalies. In this review, we summarize the data as of Moriond 2021 and present theoretical new-physics explanations from both a model-independent effective-field-theory point of view and through the building of explicit models. Throughout, we stress the complementarity of these two approaches. We also discuss combined explanations of both $B$ anomalies, and present models that also explain other problems, such as dark matter, $(g-2)_\mu$, neutrino properties, and hadronic anomalies.


title: On the underlying drag-reduction mechanisms of flow-control strategies in a transitional channel flow: temporal approach

abstract: The underlying mechanisms of three different flow-control strategies on drag reduction in a channel flow are investigated by direct numerical simulations at friction Reynolds numbers ranging from 65 to 85. These strategies include the addition of long-chain polymers, the incorporation of slip surfaces, and the application of an external body force. While it has been believed that such methods lead to a skin-friction reduction by controlling near-wall flow structures, the underlying mechanisms at play are still not as clear. In this study, a temporal analysis is employed to elucidate underlying drag-reduction mechanisms among these methods. The analysis is based on the lifetime of intermittent phases represented by the active and hibernating phases of a minimal turbulent channel flow (Xi & Graham, Phy. Rev. Lett. 2010). At a similar amount of drag reduction, the polymer and slip methods show a similar mechanism, while the body force method is different. The polymers and slip surfaces cause hibernating phases to happen more frequently, while the duration of active phases is decreased. However, the body forces cause hibernating phases to happen less frequently but prolong its duration to achieve a comparable amount of drag reduction. A possible mechanism behind the body force method is associated with its unique roller-like vortical structures formed near the wall. These structures appear to prevent interactions between inner and outer regions by which hibernating phases are prolonged. It should motivate adaptive flow-control strategies to exploit the distinct underlying mechanisms for robust control of turbulent drag at low Reynolds numbers.


title: Learning Neural Transmittance for Efficient Rendering of Reflectance Fields

abstract: Recently neural volumetric representations such as neural reflectance fields have been widely applied to faithfully reproduce the appearance of real-world objects and scenes under novel viewpoints and lighting conditions. However, it remains challenging and time-consuming to render such representations under complex lighting such as environment maps, which requires individual ray marching towards each single light to calculate the transmittance at every sampled point. In this paper, we propose a novel method based on precomputed Neural Transmittance Functions to accelerate the rendering of neural reflectance fields. Our neural transmittance functions enable us to efficiently query the transmittance at an arbitrary point in space along an arbitrary ray without tedious ray marching, which effectively reduces the time-complexity of the rendering. We propose a novel formulation for the neural transmittance function, and train it jointly with the neural reflectance fields on images captured under collocated camera and light, while enforcing monotonicity. Results on real and synthetic scenes demonstrate almost two order of magnitude speedup for renderings under environment maps with minimal accuracy loss.


title: Testing the Universality of Free Fall using correlated $^{39}$K -- $^{87}$Rb interferometers

abstract: We demonstrate how simultaneously-operated $^{39}$K -- $^{87}$Rb interferometers exhibiting a high level of correlation can be used to make competitive tests of the university of free fall. This work provides an overview of our experimental apparatus and data analysis procedure, including a detailed study of systematic effects. With a total interrogation time of $2T = 40$ ms in a compact apparatus, we reach a statistical uncertainty on the measurement of the E\"{o}tv\"{o}s parameter of $7.8 \times 10^{-8}$ after $2.4 \times 10^4$ s of integration. The main limitations to our measurement arise from a combination of wavefront aberrations, the quadratic Zeeman effect in $^{39}$K, parasitic interferometers in $^{87}$Rb, and the velocity sensitivity of our detection system. These systematic errors limit the accuracy of our measurement to $\eta = 0.9(1.6) \times 10^{-6}$. We discuss prospects for improvements using ultracold atoms at extended interrogation times.


title: Fingerprints of hot-phonon physics in time-resolved correlated quantum lattice dynamics

abstract: The time dynamics of the energy flow from electronic to lattice degrees of freedom in pump-probe setups could be strongly affected by the presence of a hot-phonon bottleneck, which can sustain longer coherence of the optically excited electronic states. Recently, hot-phonon physics has been experimentally observed and theoretically described in MgB$_2$, the electron-phonon based superconductor with $T_{\rm c}\approx 39$ K. By employing a combined ab-initio and quantum-field-theory approach and by taking MgB$_2$ as an example, here we propose a novel path for revealing the presence and characterizing the properties of hot phonons through a direct analysis of the information encoded in the lattice inter-atomic correlations. Such method exploits the underlying symmetry of the $E_{2g}$ hot modes characterized by a out-of-phase in-plane motion of the two boron atoms. Since hot phonons occur typically at high-symmetry points of the Brillouin zone, with specific symmetries of the lattice displacements, the present analysis is quite general and it could aid in revealing the hot-phonon physics in other promising materials, such as graphene, boron nitride, or black phosphorus.


title: Who's Afraid of the Supersymmetric Dark? The Standard Model vs Low-Energy Supergravity

abstract: Use of supergravity equations in astronomy and late-universe cosmology is often criticized on three grounds: (i) phenomenological success usually depends on the supergravity form for the scalar potential applying at the relevant energies; (ii) the low-energy scalar potential is extremely sensitive to quantum effects involving very massive particles and so is rarely well-approximated by classical calculations of its form; and (iii) almost all Standard Model particles count as massive for these purposes and none of these are supersymmetric. Why should Standard Model loops preserve the low-energy supergravity form even if supersymmetry is valid at energies well above the electroweak scale? We use recently developed tools for coupling supergravity to non-supersymmetric matter to estimate the loop effects of heavy non-supersymmetric particles on the low-energy effective action, and provide evidence that the supergravity form is stable against integrating out such particles (and so argues against the above objection). This suggests an intrinsically supersymmetric picture of Nature where supersymmetry survives to low energies within the gravity sector but not the visible sector (for which supersymmetry is instead non-linearly realized). We explore the couplings of both sectors in this picture and find that the presence of auxiliary fields in the gravity sector makes the visible sector share many features usually attributed to linearly realized supersymmetry although (unlike for the MSSM) a second Higgs doublet is not required for all Yukawa couplings to be non-vanishing and changes the dimension of the operator generating the Higgs mass. We discuss the naturalness of this picture and some of the implications it might have when searching for dark-sector physics.


title: Meta-analysis of Uniform Scaling Factors for Harmonic Frequency Calculations

abstract: Vibrational frequency calculations performed under the harmonic approximation are widespread across chemistry. However, it is well-known that the calculated harmonic frequencies tend to systematically overestimate experimental fundamental frequencies; a limitation commonly overcome with multiplicative scaling factors. In practice, multiplicative scaling factors are derived for each individual model chemistry choice (i.e., a level of theory and basis set pair), where performance is judged by, for example, the root-mean square error (RMSE) between the predicted scaled and experimental frequencies. However, despite the overwhelming number of scaling factors reported in the literature and model chemistry approximations available, there is little guidance for users on appropriate model chemistry choices for harmonic frequency calculations. Here, we compile and analyse the data for 1495 scaling factors calculated using 141 levels of theory and 109 basis sets. Our meta-analysis of this data shows that scaling factors and RMSE approach convergence with only hybrid functionals and double-zeta basis sets, with anharmonicity error already dominating model chemistry errors. Noting inconsistent data and the lack of independent testing, we can nevertheless conclude that a minimum error of 25cm-1 -- arising from insufficiently accurate treatment of anharmonicity -- is persistent regardless of the model chemistry choice. Based on the data we compiled and cautioning the need for a future systematic benchmarking study, we recommend wB97X-D/def2-TZVP for most applications and B2PLYP/def2-TZVPD for superior intensity predictions. With a smaller benchmark set, direct comparison strongly prefers wB97X-D/6-31G* to B3LYP/6-31G*.


title: Large Phonon Thermal Hall Conductivity in a Simple Antiferromagnetic Insulator

abstract: Phonons are known to generate a thermal Hall effect in certain insulators, including oxides with rare-earth impurities, quantum paraelectrics, multiferroic materials and cuprate Mott insulators. In each case, a special feature of the material is presumed relevant for the underlying mechanism that confers chirality to phonons in a magnetic field. The question is whether a phonon Hall effect is an unusual occurrence - linked to special characteristics such as skew scattering off rare-earth impurities, structural domains, ferroelectricity, ferromagnetism - or a much more common property of insulators than hitherto believed. To help answer this question, we have turned to a simple insulator, with none of the previously encountered special features: the cubic antiferromagnet Cu$_3$TeO$_6$. We find that it has the largest thermal Hall conductivity $\kappa_{\rm{xy}}$ of any insulator so far. We show that this record-high $\kappa_{\rm{xy}}$ signal is due to phonons and it does not require the presence of magnetic order, as it persists above the ordering temperature. We conclude that the phonon Hall effect is likely to be a fairly common property of solids.


title: Optomechanical quantum entanglement mediated by acoustic phonon fields

abstract: We present exact solutions for the quantum time evolution of two spatially separated, local inductor-capacitor (LC) oscillators that are coupled optomechanically to a long elastic strip that functions as a quantum thermal acoustic field environment. We show that the optomechanical coupling to the acoustic environment gives rise to causal entanglement dynamics between the two LC oscillators in the absence of resonant photon exchange between them, and that significant entanglement develops regardless of the environment temperature. Such a process establishes that distributed entanglement may be generated between superconducting qubits via a connected phonon bus bar, without the need for resonant phonon release and capture.


title: Three-Loop Corrections to Lamb Shift in Muonium and Positronium

abstract: We calculate hard spin-independent three-loop radiative corrections to energy levels in muonium and positronium which are due to radiative insertions in two-photon exchange diagrams. These corrections could be relevant for the new generation of precise $1S-2S$ and $2S-2P$ measurements in muonium and positronium.


title: A Variational Graph Autoencoder for Manipulation Action Recognition and Prediction

abstract: Despite decades of research, understanding human manipulation activities is, and has always been, one of the most attractive and challenging research topics in computer vision and robotics. Recognition and prediction of observed human manipulation actions have their roots in the applications related to, for instance, human-robot interaction and robot learning from demonstration. The current research trend heavily relies on advanced convolutional neural networks to process the structured Euclidean data, such as RGB camera images. These networks, however, come with immense computational complexity to be able to process high dimensional raw data. Different from the related works, we here introduce a deep graph autoencoder to jointly learn recognition and prediction of manipulation tasks from symbolic scene graphs, instead of relying on the structured Euclidean data. Our network has a variational autoencoder structure with two branches: one for identifying the input graph type and one for predicting the future graphs. The input of the proposed network is a set of semantic graphs which store the spatial relations between subjects and objects in the scene. The network output is a label set representing the detected and predicted class types. We benchmark our new model against different state-of-the-art methods on two different datasets, MANIAC and MSRC-9, and show that our proposed model can achieve better performance. We also release our source code https://github.com/gamzeakyol/GNet.


title: A resolution criterion based on characteristic timescales for MHD simulations of molecular clouds

abstract: We investigate the effect of numerical magnetic diffusion in magnetohydrodynamic simulations of magnetically supported molecular clouds. To this end, we have performed numerical studies on adaptive mesh isothermal simulations of marginally sub-critical molecular clouds. We find that simulations with low and intermediate resolutions collapse, contrary to what is theoretically expected. However, the simulation with the highest numerical resolution oscillates around an equilibrium state without collapsing. In order to quantify the numerical diffusion of the magnetic field, we ran a second suit of current-sheet simulations in which the numerical magnetic diffusion coefficient can be directly measured, and computed the corresponding diffusion times at various numerical resolutions. On this basis, we propose a criterion for the resolution of magnetic fields in MHD simulations based on requiring that the diffusion time to be larger than the characteristic timescale of the physical process responsible for the dynamic evolution of the structure.


title: The Pareto Frontier of model selection for general Contextual Bandits

abstract: Recent progress in model selection raises the question of the fundamental limits of these techniques. Under specific scrutiny has been model selection for general contextual bandits with nested policy classes, resulting in a COLT2020 open problem. It asks whether it is possible to obtain simultaneously the optimal single algorithm guarantees over all policies in a nested sequence of policy classes, or if otherwise this is possible for a trade-off $\alpha\in[\frac{1}{2},1)$ between complexity term and time: $\ln(|\Pi_m|)^{1-\alpha}T^\alpha$. We give a disappointing answer to this question. Even in the purely stochastic regime, the desired results are unobtainable. We present a Pareto frontier of up to logarithmic factors matching upper and lower bounds, thereby proving that an increase in the complexity term $\ln(|\Pi_m|)$ independent of $T$ is unavoidable for general policy classes. As a side result, we also resolve a COLT2016 open problem concerning second-order bounds in full-information games.


title: Generating GitHub Repository Descriptions: A Comparison of Manual and Automated Approaches

abstract: Given the vast number of repositories hosted on GitHub, project discovery and retrieval have become increasingly important for GitHub users. Repository descriptions serve as one of the first points of contact for users who are accessing a repository. However, repository owners often fail to provide a high-quality description; instead, they use vague terms, the purpose of the repository is poorly explained, or the description is omitted entirely. In this work, we examine the current practice of writing GitHub repository descriptions. Our investigation leads to the proposal of the LSP (Language, Software technology, and Purpose) template to formulate good descriptions for GitHub repositories that are clear, concise, and informative. To understand the extent to which current automated techniques can support generating repository descriptions, we compare the performance of state-of-the-art text summarization methods on this task. Finally, our user study with GitHub users reveals that automated summarization can adequately be used for default description generation for GitHub repositories, while the descriptions which follow the LSP template offer the most effective instrument for communicating with GitHub users.


title: A New Method for Estimating Starspot Lifetimes Based on Autocorrelation Functions

abstract: We present a method that utilizes autocorrelation functions from long-term precision broadband differential light curves to estimate the average lifetimes of starspot groups for two large sample of Kepler stars: stars with and without previously known rotation periods. Our method is calibrated by comparing the strengths of the first few normalized autocorrelation peaks using ensembles of models that have various starspot lifetimes. We find that we must mix models of short and long lifetimes together (in heuristically determined ratios) to align the models with the Kepler data. Our fundamental result is that short starspot group lifetimes (1-4 rotations) are implied when the first normalized peak is weaker than about 0.4, long lifetimes (15 or greater) are implied when it is greater than about 0.7, and in between are the intermediate cases. Rotational lifetimes can be converted to physical lifetimes if the rotation period is known. Stars with shorter rotation periods tend to have longer rotational (but not physical) spot lifetimes, and cooler stars tend to have longer physical spot lifetimes than warmer stars with the same rotation period. The distributions of the physical lifetimes are lognormal for both samples and generally longer in the first sample. The shorter lifetimes in the stars without known periods probably explain why their periods are difficult to measure. Some stars exhibit longer than average physical starspot lifetimes; their percentage drops with increasing temperature from nearly half at 3000K to nearly zero for hotter than 6000K.


title: Generative Flows as a General Purpose Solution for Inverse Problems

abstract: Due to the success of generative flows to model data distributions, they have been explored in inverse problems. Given a pre-trained generative flow, previous work proposed to minimize the 2-norm of the latent variables as a regularization term. The intuition behind it was to ensure high likelihood latent variables that produce the closest restoration. However, high-likelihood latent variables may generate unrealistic samples as we show in our experiments. We therefore propose a solver to directly produce high-likelihood reconstructions. We hypothesize that our approach could make generative flows a general purpose solver for inverse problems. Furthermore, we propose 1 x 1 coupling functions to introduce permutations in a generative flow. It has the advantage that its inverse does not require to be calculated in the generation process. Finally, we evaluate our method for denoising, deblurring, inpainting, and colorization. We observe a compelling improvement of our method over prior works.


title: Dynamic bulk-boundary correspondence for anomalous Floquet topology

abstract: Periodically driven systems with internal and spatial symmetries can exhibit a variety of anomalous boundary behaviors at both the zero and $\pi$ quasienergies despite the trivial bulk Floquet bands. These phenomena are called anomalous Floquet topology (AFT) as they are unconnected from their static counterpart, arising from the winding of the time evolution unitary rather than the bulk Floquet bands at the end of the driving period. In this paper, we systematically derive the first and inversion-symmetric second-order AFT bulk-boundary correspondence for Altland-Zirnbauer (AZ) classes BDI, D, DIII, AII. For each AZ class, we start a dimensional hierarchy with a parent dimension having $\mathbb{Z}$ classification, then use it as an interpolating map to classify the lower-dimensional descendants. From the Atiyah-Hirzebruch spectral sequence (AHSS), we identify the subspace that contains topological information and faithfully derive the AFT bulk-boundary correspondence for both the parent and descendants. Our theory provides analytic tools for out-of-equilibrium topological phenomena.


title: Towards Realistic Market Simulations: a Generative Adversarial Networks Approach

abstract: Simulated environments are increasingly used by trading firms and investment banks to evaluate trading strategies before approaching real markets. Backtesting, a widely used approach, consists of simulating experimental strategies while replaying historical market scenarios. Unfortunately, this approach does not capture the market response to the experimental agents' actions. In contrast, multi-agent simulation presents a natural bottom-up approach to emulating agent interaction in financial markets. It allows to set up pools of traders with diverse strategies to mimic the financial market trader population, and test the performance of new experimental strategies. Since individual agent-level historical data is typically proprietary and not available for public use, it is difficult to calibrate multiple market agents to obtain the realism required for testing trading strategies. To addresses this challenge we propose a synthetic market generator based on Conditional Generative Adversarial Networks (CGANs) trained on real aggregate-level historical data. A CGAN-based "world" agent can generate meaningful orders in response to an experimental agent. We integrate our synthetic market generator into ABIDES, an open source simulator of financial markets. By means of extensive simulations we show that our proposal outperforms previous work in terms of stylized facts reflecting market responsiveness and realism.


title: Controlling Smart Propagation Environments: Long-Term versus Short-Term Phase Shift Optimization

abstract: Reconfigurable intelligent surfaces (RISs) have recently gained significant interest as an emerging technology for future wireless networks. This paper studies an RIS-assisted propagation environment, where a single-antenna source transmits data to a single-antenna destination in the presence of a weak direct link. We analyze and compare RIS designs based on long-term and short-term channel statistics in terms of coverage probability and ergodic rate. For the considered optimization designs, closed-form expressions for the coverage probability and ergodic rate are derived. We use numerical simulations to analyze and compare against analytic results in finite samples. Also, we show that the considered optimal phase shift designs outperform several heuristic benchmarks.


title: Uncertainty quantification in non-rigid image registration via stochastic gradient Markov chain Monte Carlo

abstract: We develop a new Bayesian model for non-rigid registration of three-dimensional medical images, with a focus on uncertainty quantification. Probabilistic registration of large images with calibrated uncertainty estimates is difficult for both computational and modelling reasons. To address the computational issues, we explore connections between the Markov chain Monte Carlo by backpropagation and the variational inference by backpropagation frameworks, in order to efficiently draw samples from the posterior distribution of transformation parameters. To address the modelling issues, we formulate a Bayesian model for image registration that overcomes the existing barriers when using a dense, high-dimensional, and diffeomorphic transformation parametrisation. This results in improved calibration of uncertainty estimates. We compare the model in terms of both image registration accuracy and uncertainty quantification to VoxelMorph, a state-of-the-art image registration model based on deep learning.


title: Exploring System Performance of Continual Learning for Mobile and Embedded Sensing Applications

abstract: Continual learning approaches help deep neural network models adapt and learn incrementally by trying to solve catastrophic forgetting. However, whether these existing approaches, applied traditionally to image-based tasks, work with the same efficacy to the sequential time series data generated by mobile or embedded sensing systems remains an unanswered question. To address this void, we conduct the first comprehensive empirical study that quantifies the performance of three predominant continual learning schemes (i.e., regularization, replay, and replay with examples) on six datasets from three mobile and embedded sensing applications in a range of scenarios having different learning complexities. More specifically, we implement an end-to-end continual learning framework on edge devices. Then we investigate the generalizability, trade-offs between performance, storage, computational costs, and memory footprint of different continual learning methods. Our findings suggest that replay with exemplars-based schemes such as iCaRL has the best performance trade-offs, even in complex scenarios, at the expense of some storage space (few MBs) for training examples (1% to 5%). We also demonstrate for the first time that it is feasible and practical to run continual learning on-device with a limited memory budget. In particular, the latency on two types of mobile and embedded devices suggests that both incremental learning time (few seconds - 4 minutes) and training time (1 - 75 minutes) across datasets are acceptable, as training could happen on the device when the embedded device is charging thereby ensuring complete data privacy. Finally, we present some guidelines for practitioners who want to apply a continual learning paradigm for mobile sensing tasks.


title: Optimal cooling of an internally heated disc

abstract: Motivated by the search for sharp bounds on turbulent heat transfer as well as the design of optimal heat exchangers, we consider incompressible flows that most efficiently cool an internally heated disc. Heat enters via a distributed source, is passively advected and diffused, and exits through the boundary at a fixed temperature. We seek an advecting flow to optimize this exchange. Previous work on energy-constrained cooling with a constant source has conjectured that global optimizers should resemble convection rolls; we prove one-sided bounds on energy-constrained cooling corresponding to, but not resolving, this conjecture. In the case of an enstrophy constraint, our results are more complete: we construct a family of self-similar, tree-like ``branching flows'' whose cooling we prove is within a logarithm of globally optimal. These results hold for general space- and time-dependent source-sink distributions that add more heat than they remove. Our main technical tool is a non-local Dirichlet-like variational principle for bounding solutions of the inhomogeneous advection-diffusion equation with a divergence-free velocity.


title: Self-aware Social Learning over Graphs

abstract: In this paper we study the problem of social learning under multiple true hypotheses and self-interested agents which exchange information over a graph. In this setup, each agent receives data that might be generated from a different hypothesis (or state) than the data other agents receive. In contrast to the related literature in social learning, which focuses on showing that the network achieves consensus, here we study the case where every agent is self-interested and wants to find the hypothesis that generates its own observations. However, agents do not know which ones of their peers wants to find the same state with them and as a result they do not know which agents they should cooperate with. To this end, we propose a scheme with adaptive combination weights and study the consistency of the agents' learning process. The scheme allows each agent to identify and collaborate with neighbors that observe the same hypothesis, while excluding others, thus resulting in improved performance compared to both non-cooperative learning and cooperative social learning solutions. We analyze the asymptotic behavior of agents' beliefs under the proposed social learning algorithm and provide sufficient conditions that enable all agents to correctly identify their true hypotheses. The theoretical analysis is corroborated by numerical simulations.


title: Emulation of physical processes with Emukit

abstract: Decision making in uncertain scenarios is an ubiquitous challenge in real world systems. Tools to deal with this challenge include simulations to gather information and statistical emulation to quantify uncertainty. The machine learning community has developed a number of methods to facilitate decision making, but so far they are scattered in multiple different toolkits, and generally rely on a fixed backend. In this paper, we present Emukit, a highly adaptable Python toolkit for enriching decision making under uncertainty. Emukit allows users to: (i) use state of the art methods including Bayesian optimization, multi-fidelity emulation, experimental design, Bayesian quadrature and sensitivity analysis; (ii) easily prototype new decision making methods for new problems. Emukit is agnostic to the underlying modeling framework and enables users to use their own custom models. We show how Emukit can be used on three exemplary case studies.


title: When cardinals determine the power set: inner models and H\"artig quantifier logic

abstract: We make use of some observations on the core model, for example assuming $V=L [ E ]$, and that there is no inner model with a Woodin cardinal, and $M$ is an inner model with the same cardinals as $V$, then $V=M$. We conclude in this latter situation that "$x=\mathcal{P} ( y )$" is $\Sigma_{1} ( Card )$ where $Card$ is a predicate true of just the infinite cardinals. It is known that this implies the validities of second order logic are reducible to $V_I$ the set of validities of the H\"artig quantifier logic. We draw some further conclusions on the L\"owenheim number, $\ell_{I}$ of the latter logic: that if no $L[E]$ model has a cardinal strong up to an $\aleph$-fixed point, and $\ell_{I}$ is less than the least weakly inaccessible $\delta$, then (i) $\ell_I$ is a limit of measurable cardinals of $K$; (ii) the Weak Covering Lemma holds at $\delta$.


title: Germanium monosulfide as a natural platform for highly anisotropic THz polaritons

abstract: Terahertz (THz) electromagnetic radiation is key to optically access collective excitations such as magnons (spins), plasmons (electrons), or phonons (atomic vibrations), thus bridging between optics and solid-state physics. Confinement of THz light to the nanometer length scale is desirable for local probing of such excitations in low dimensional systems, thereby inherently circumventing the large footprint and low spectral density of far-field THz optics. For that purpose, phonon polaritons (PhPs, i.e., light coupled to lattice vibrations in polar crystals) in anisotropic van der Waals (vdW) materials have recently emerged as a promising platform for THz nanooptics; yet the amount of explored, viable materials is still exiguous. Hence, there is a demand for the exploration of novel materials that feature not only THz PhPs at different spectral regimes, but also exhibit unique anisotropic (directional) electrical, thermoelectric, and vibronic properties. To that end, we introduce here the semiconducting alpha-germanium(II) sulfide (GeS) as an intriguing candidate. By employing THz nano-spectroscopy supported by theoretical analysis, we provide a thorough characterization of the different in-plane hyperbolic and elliptical PhP modes in GeS. We find not only PhPs with long life times ($\tau$ > 2 ps) and excellent THz light confinement ($\lambda_0/\lambda$ > 45), but also an intrinsic, phonon-induced anomalous dispersion as well as signatures of naturally occurring PhP canalization within one single GeS slab.


title: Efficient ISDA Initial Margin Calculations Using Least Squares Monte-Carlo

abstract: Non-cleared bilateral OTC derivatives between two financial firms or systemically important non-financial entities are subject to regulations that require the posting of initial and variation margin. The ISDA standard approach (SIMM) provides a way for computing the initial margin. It involves computing sensitivities of the contracts with respect to several market factors. In this paper, the authors extend the well known LSMC technique to efficiently estimate the sensitivities required in the ISDA SIMM methodology.


title: Fast PDE-constrained optimization via self-supervised operator learning

abstract: Design and optimal control problems are among the fundamental, ubiquitous tasks we face in science and engineering. In both cases, we aim to represent and optimize an unknown (black-box) function that associates a performance/outcome to a set of controllable variables through an experiment. In cases where the experimental dynamics can be described by partial differential equations (PDEs), such problems can be mathematically translated into PDE-constrained optimization tasks, which quickly become intractable as the number of control variables and the cost of experiments increases. In this work we leverage physics-informed deep operator networks (DeepONets) -- a self-supervised framework for learning the solution operator of parametric PDEs -- to build fast and differentiable surrogates for rapidly solving PDE-constrained optimization problems, even in the absence of any paired input-output training data. The effectiveness of the proposed framework will be demonstrated across different applications involving continuous functions as control or design variables, including time-dependent optimal control of heat transfer, and drag minimization of obstacles in Stokes flow. In all cases, we observe that DeepONets can minimize high-dimensional cost functionals in a matter of seconds, yielding a significant speed up compared to traditional adjoint PDE solvers that are typically costly and limited to relatively low-dimensional control/design parametrizations.


title: Scalable Bayesian divergence time estimation with ratio transformations

abstract: Divergence time estimation is crucial to provide temporal signals for dating biologically important events, from species divergence to viral transmissions in space and time. With the advent of high-throughput sequencing, recent Bayesian phylogenetic studies have analyzed hundreds to thousands of sequences. Such large-scale analyses challenge divergence time reconstruction by requiring inference on highly-correlated internal node heights that often become computationally infeasible. To overcome this limitation, we explore a ratio transformation that maps the original N - 1 internal node heights into a space of one height parameter and N - 2 ratio parameters. To make analyses scalable, we develop a collection of linear-time algorithms to compute the gradient and Jacobian-associated terms of the log-likelihood with respect to these ratios. We then apply Hamiltonian Monte Carlo sampling with the ratio transform in a Bayesian framework to learn the divergence times in four pathogenic virus phylogenies: West Nile virus, rabies virus, Lassa virus and Ebola virus. Our method both resolves a mixing issue in the West Nile virus example and improves inference efficiency by at least 5-fold for the Lassa and rabies virus examples. Our method also makes it now computationally feasible to incorporate mixed-effects molecular clock models for the Ebola virus example, confirms the findings from the original study and reveals clearer multimodal distributions of the divergence times of some clades of interest.


title: Squeezing multiple soft particles into a constriction: transition to clogging

abstract: We study numerically how multiple deformable capsules squeeze into a constriction. This situation is largely encountered in microfluidic chips designed to manipulate living cells, which are soft entities. We use fully three-dimensional simulations based on the lattice Boltzmann method to compute the flow of the suspending fluid, and on the immersed boundary method to achieve the two-way fluid-structure interaction. The mechanics of the capsule membrane elasticity is computed with the finite element method. We obtain two main states: continuous passage of the particles, and their blockage that leads to clogging the constriction. The transition from one state to another is dictated by the ratio between the size of the capsules and the constriction width, and by the capsule membrane deformability. This latter is found to enhance particle passage through narrower constrictions, where rigid particles with similar diameter are blocked and lead to clogging.

categories: cond-mat.soft physics.bio-ph physics.comp-ph physics.flu-dyn


title: Adaptive Probabilistic Model for Energy-Efficient Distance-based Clustering in WSNs (Adapt-P): A LEACH-based Analytical Study

abstract: Network lifetime and energy consumption of data transmission have been primary Quality of Service (QoS) obligations in Wireless Sensor Networks (WSNs). The environment of a WSN is often organized into clusters to mitigate the management complexity of such obligations. However, the distance between Sensor Nodes (SNs) and the number of clusters per round are vital factors that affect QoS performance of a WSN. A designer's conundrum resolves around the desire to sustain a balance between the limited residual energy of SNs and the demand for prolonged network lifetime. Any imbalance in controlling such objectives results in either QoS penalties due to draining SN energies, or an over-cost environment that is significantly difficult to distribute and operate. Low-Energy Adaptive Clustering Hierarchy (LEACH) is a distributed algorithm proposed to tackle such difficulties. Proposed LEACH-based algorithms focus on residual energies of SNs to compute a probability function that selects cluster-heads and an optimal energy-efficient path toward a destination SN. Nevertheless, these algorithms do not consider variations in network's state at run-time. Such a state changes in an adaptive manner according to existing network structures and conditions. Thus, cluster-heads per round are not elected adaptively depending on the state and distances between SNs. This paper proposes an energy-efficient adaptive distance-based clustering called Adapt-P, in which an adaptive probability function is developed to formulate clusters. A near-optimal distance between each cluster-head and its cluster-members is formulated so that energy consumption of the network is mitigated and network lifetime is maximized. The cluster-head selection probability is adapted at the end of each round based on the maximum number of cluster-heads permitted per round found a priori and the number of alive SNs in the network.


title: A new Wilson Line-based classical action for gluodynamics

abstract: We develop a new classical action that in addition to $\mathrm{MHV}$ vertices contains also $\mathrm{N^kMHV}$ vertices, where $1\leq k \leq n-4$ with $n$ the number of external legs. The lowest order vertex is the four-point MHV vertex -- there is no three point vertex and thus the amplitude calculation involves fewer vertices than in the CSW method and, obviously, considerably fewer than in the standard Yang-Mills action. The action is obtained by a canonical transformation of the Yang-Mills action in the light-cone gauge, where the field transformations are based on the Wilson line functionals.


title: Oscillating Wandering Domains for p-Adic Transcendental Entire Maps

abstract: We give examples of transcendental entire maps over $\mathbb{C}_p$ having an oscillating wandering Fatou component.


title: Negotiating Networks in Oligopoly Markets for Price-Sensitive Products

abstract: We present a novel framework to learn functions that estimate decisions of sellers and buyers simultaneously in an oligopoly market for a price-sensitive product. In this setting, the aim of the seller network is to come up with a price for a given context such that the expected revenue is maximized by considering the buyer's satisfaction as well. On the other hand, the aim of the buyer network is to assign probability of purchase to the offered price to mimic the real world buyers' responses while also showing price sensitivity through its action. In other words, rejecting the unnecessarily high priced products. Similar to generative adversarial networks, this framework corresponds to a minimax two-player game. In our experiments with simulated and real-world transaction data, we compared our framework with the baseline model and demonstrated its potential through proposed evaluation metrics.


title: Observational constraints on the maximum energies of accelerated particles in supernova remnants: low maximum energies and a large variety

abstract: Supernova remnants (SNRs) are thought to be the most promising sources of Galactic cosmic rays. One of the principal questions is whether they are accelerating particles up to the maximum energy of Galactic cosmic rays ($\sim$ PeV). In this work, a systematic study of gamma-ray emitting SNRs is conducted as an advanced study of Suzuki et al. 2021. Our purpose is to newly measure the evolution of maximum particle energies with increased statistics and better age estimates. We model their gamma-ray spectra to constrain the particle-acceleration parameters. Two candidates of the maximum energy of freshly accelerated particles, the gamma-ray cutoff and break energies, are found to be well below PeV. We also test a spectral model that includes both the freshly accelerated and escaping particles to estimate the maximum energies more reliably, but no tighter constraints are obtained with current statistics. The average time dependences of the cutoff energy ($\propto t^{-0.81 \pm 0.24}$) and break energy ($\propto t^{-0.77 \pm 0.23}$) cannot be explained with the simplest acceleration condition (Bohm limit), and requires shock-ISM (interstellar medium) interaction. The average maximum energy during lifetime is found to be $\lesssim 20$ TeV $(t_{\rm M}/1~{\rm kyr})^{-0.8}$ with $t_{\rm M}$ being the age at the maximum, which reaches PeV if $t_{\rm M} \lesssim 10$ yr. The maximum energies during lifetime are suggested to have a variety of 1.1--1.8 dex from object to object. Although we cannot isolate the cause of this variety, this work provides an important clue to understand the microphysics of particle acceleration in SNRs.


title: Recurrence equations involving different orthogonal polynomial sequences and applications

abstract: Consider $\{p_n\}_{n=0}^{\infty}$, a sequence of polynomials orthogonal with respect to $w(x)>0$ on $(a,b)$, and polynomials $\{g_{n,k}\}_{n=0}^{\infty},k \in \mathbb{N}_0$, orthogonal with respect to $c_k(x)w(x)>0$ on $(a,b)$, where $c_{k}(x)$ is a polynomial of degree $k$ in $x$. We show how Christoffel's formula can be used to obtain mixed three-term recurrence equations involving the polynomials $p_n$, $p_{n-1}$ and $g_{n-m,k},m\in\{2,3,\dots, n-1\}$. In order for the zeros of $p_n$ and $G_{m-1}g_{n-m,k}$ to interlace (assuming $p_n$ and $g_{n-m,k}$ are co-prime), the coefficient of $p_{n-1}$, namely $G_{m-1}$, should be of exact degree $m-1$, in which case restrictions on the parameter $k$ are necessary. The zeros of $G_{m-1}$ can be considered to be inner bounds for the extreme zeros of the (classical or $q$-classical) orthogonal polynomial $p_n$ and we give examples to illustrate the accuracy of these bounds. Because of the complexity the mixed three-term recurrence equations in each case, algorithmic tools, mainly Zeilberger's algorithm and its $q$-analogue, are used to obtain them.


title: Reconciling Risk Allocation and Prevalence Estimation in Public Health Using Batched Bandits

abstract: In many public health settings, there is a perceived tension between allocating resources to known vulnerable areas and learning about the overall prevalence of the problem. Inspired by a door-to-door Covid-19 testing program we helped design, we combine multi-armed bandit strategies and insights from sampling theory to demonstrate how to recover accurate prevalence estimates while continuing to allocate resources to at-risk areas. We use the outbreak of an infectious disease as our running example. The public health setting has several characteristics distinguishing it from typical bandit settings, such as distribution shift (the true disease prevalence is changing with time) and batched sampling (multiple decisions must be made simultaneously). Nevertheless, we demonstrate that several bandit algorithms are capable out-performing greedy resource allocation strategies, which often perform worse than random allocation as they fail to notice outbreaks in new areas.


title: Institutional Incentives for the Evolution of Committed Cooperation: Ensuring Participation is as Important as Enhancing Compliance

abstract: Both conventional wisdom and empirical evidence suggests that arranging a prior commitment or agreement before an interaction enhances the chance of reaching mutual cooperation. Yet it is not clear what mechanisms can promote the participation in and compliance with such a commitment, especially when the former is costly and deviating from the latter is profitable. Prior work either considers regimented commitments where compensation is assumed enforceable from dishonest committers, or assume implicit commitments from every individual (so they are all in and thus being treated as such). Here we develop a theory of participation and compliance with respect to an explicit prior commitment under institutional incentives where individuals, at first, decide whether or not to join a cooperative agreement to play a one-shot social dilemma game. Using a mathematical model, we determine when participating in a costly commitment and complying with it, is an evolutionary stable strategy (ESS) when playing against all other possible strategies, and results in high levels of cooperation in the population. We show that, given a sufficient budget for providing incentives, reward of commitment compliant behaviours better promotes cooperation than punishment of non-compliant ones. Moreover, by sparing part of this budget for rewarding those who are willing to participate in a commitment, the overall frequency of cooperation can be significantly enhanced, for both reward and punishment. Finally, we find that, surprisingly, the presence of errors in a participation decision favours evolutionary stability of commitment compliant strategies and higher levels of cooperation.


title: Calibrating DFT formation enthalpy calculations by multi-fidelity machine learning

abstract: Machine learning materials properties measured by experiments is valuable yet difficult due to the limited amount of experimental data. In this work, we use a multi-fidelity random forest model to learn the experimental formation enthalpy of materials with prediction accuracy higher than the empirically corrected PBE functional (PBEfe) and meta-GGA functional (SCAN), and it outperforms the hotly studied deep neural-network based representation learning and transfer learning. We then use the model to calibrate the DFT formation enthalpy in the Materials Project database, and discover materials with underestimated stability. The multi-fidelity model is also used as a data-mining approach to find how DFT deviates from experiments by the explaining the model output.


title: History Aware Multimodal Transformer for Vision-and-Language Navigation

abstract: Vision-and-language navigation (VLN) aims to build autonomous visual agents that follow instructions and navigate in real scenes. To remember previously visited locations and actions taken, most approaches to VLN implement memory using recurrent states. Instead, we introduce a History Aware Multimodal Transformer (HAMT) to incorporate a long-horizon history into multimodal decision making. HAMT efficiently encodes all the past panoramic observations via a hierarchical vision transformer (ViT), which first encodes individual images with ViT, then models spatial relation between images in a panoramic observation and finally takes into account temporal relation between panoramas in the history. It, then, jointly combines text, history and current observation to predict the next action. We first train HAMT end-to-end using several proxy tasks including single step action prediction and spatial relation prediction, and then use reinforcement learning to further improve the navigation policy. HAMT achieves new state of the art on a broad range of VLN tasks, including VLN with fine-grained instructions (R2R, RxR), high-level instructions (R2R-Last, REVERIE), dialogs (CVDN) as well as long-horizon VLN (R4R, R2R-Back). We demonstrate HAMT to be particularly effective for navigation tasks with longer trajectories.


title: Reconstruction of Mutation Laws in Heterogeneous Tumours with Local and Nonlocal Dynamics

abstract: Cancer cell mutations occur when cells undergo multiple cell divisions, and these mutations can be spontaneous or environmentally-induced. The mechanisms that promote and sustain these mutations are still not fully understood. This study deals with the identification (or reconstruction) of the usually unknown cancer cell mutation law, which lead to the transformation of a primary tumour cell population into a secondary, more aggressive cell population. We focus on local and nonlocal mathematical models for cell dynamics and movement, and identify these mutation laws from macroscopic tumour snapshot data collected at some later stage in the tumour evolution. In a local cancer invasion model, we first reconstruct the mutation law when we assume that the mutations depend only on the surrounding cancer cells (i.e., the ECM plays no role in mutations). Second, we assume that the mutations depend on the ECM only, and we reconstruct the mutation law in this case. Third, we reconstruct the mutation when we assume that there is no prior knowledge about the mutations. Finally, for the nonlocal cancer invasion model, we reconstruct the mutation law that depends on the cancer cells and on the ECM. For these numerical reconstructions, our approximations are based on the finite difference method combined with the finite elements method. As the inverse problem is ill-posed, we use the Tikhonov regularisation technique in order to regularise the solution. Stability of the solution is examined by adding additive noise into the measurements.


title: Physics informed machine learning with Smoothed Particle Hydrodynamics: Hierarchy of reduced Lagrangian models of turbulence

abstract: Building efficient, accurate and generalizable reduced order models of developed turbulence remains a major challenge. This manuscript approaches this problem by developing a hierarchy of parameterized reduced Lagrangian models for turbulent flows, and investigates the effects of enforcing physical structure through Smoothed Particle Hydrodynamics (SPH) versus relying on neural networks (NN)s as universal function approximators. Starting from Neural Network (NN) parameterizations of a Lagrangian acceleration operator, this hierarchy of models gradually incorporates a weakly compressible and parameterized SPH framework, which enforces physical symmetries, such as Galilean, rotational and translational invariances. Within this hierarchy, two new parameterized smoothing kernels are developed in order to increase the flexibility of the learn-able SPH simulators. For each model we experiment with different loss functions which are minimized using gradient based optimization, where efficient computations of gradients are obtained by using Automatic Differentiation (AD) and Sensitivity Analysis (SA). Each model within the hierarchy is trained on two data sets associated with weekly compressible Homogeneous Isotropic Turbulence (HIT): (1) a validation set using weakly compressible SPH; and (2) a high fidelity set from Direct Numerical Simulations (DNS). Numerical evidence shows that encoding more SPH structure improves generalizability to different turbulent Mach numbers and time shifts, and that including the novel parameterized smoothing kernels improves the accuracy of SPH at the resolved scales.


title: Commutators of multilinear strongly singular integrals on non-homogeneous metric measure spaces

abstract: Let $(X,d,\mu)$ denotes non-homogeneous metric measure space satisfying geometrically doubling and the upper doubling measure condition. In this paper, the boundedness in Lebesgue spaces for two kinds of commutators, which are iterated commutators and commutators in summation form, generated by multilinear strongly singular integral operators with RBMO$(\mu)$ function on non-homogeneous metric measure spaces $(X,d,\mu)$ is obtained.


title: A Novel AQC Factoring Algorithm

abstract: Due to recent technological advances, actual quantum devices are being constructed and used to perform computations. As a result, many classical problems are being restated so as to be solved on quantum computers. Some examples include satisfiability problems; clustering and classification; protein folding; and simulating many-body systems. Converting these classical problems to a quantum framework is not always straightforward. As such, instances where researchers explicitly elucidate the conversion process are not only valuable in their own right, but are likely to spawn new ideas and creative ways in regards to problem solving. In this paper, we propose a classical factoring algorithm, which we then convert into a quantum framework. Along the way, we discuss the subtle similarities and differences between the approaches, and provide a general comparison of their performance. It is our desire to not only introduce an interesting approach to factoring, but to hopefully promote more creative ways to solving problems using quantum computers. The key to our algorithm is that we convert the factoring problem to a graph theory problem using elements from group theory. The move to a graph-theoretic approach ultimately eases the transition to a quantum setting.


title: Products of reflections in smooth Bruhat intervals

abstract: A permutation is called smooth if the corresponding Schubert variety is smooth. Gilboa and Lapid prove that in the symmetric group, multiplying the reflections below a smooth element $w$ in Bruhat order in a compatible order yields back the element $w$. We strengthen this result by showing that such a product in fact determines a saturated chain $e \to w$ in Bruhat order, and that this property characterizes smooth elements.


title: EarthGAN: Can we visualize the Earth's mantle convection using a surrogate model?

abstract: Scientific simulations are often used to gain insight into foundational questions. However, many potentially useful simulation results are difficult to visualize without powerful computers. In this research, we seek to build a surrogate model, using a generative adversarial network, to allow for the visualization of the Earth's Mantle Convection data set on readily accessible hardware. We present our preliminary method and results, and all code is made publicly available. The preliminary results show that a surrogate model of the Earth's Mantle Convection data set can generate useful results. A comparison to the "ground-truth" is provided.


title: Cosmography with Standard Sirens from the Einstein Telescope

abstract: We discuss the power of third-generation gravitational wave detectors to constrain cosmographic parameters in the case of electromagnetically bright standard sirens focusing on the specific case of the Einstein Telescope. We analyze the impact that the redshift source distribution, the number of detections and the observational error in the luminosity distance have on the inference of the first cosmographic parameters, and show that with a few hundreds detections the Hubble constant can be recovered at sub-percent level whereas the deceleration parameter at a few percent level, both with negligible bias.


title: Exposure of occupations to technologies of the fourth industrial revolution

abstract: The fourth industrial revolution (4IR) is likely to have a substantial impact on the economy. Companies need to build up capabilities to implement new technologies, and automation may make some occupations obsolete. However, where, when, and how the change will happen remain to be determined. Robust empirical indicators of technological progress linked to occupations can help to illuminate this change. With this aim, we provide such an indicator based on patent data. Using natural language processing, we calculate patent exposure scores for more than 900 occupations, which represent the technological progress related to them. To provide a lens on the impact of the 4IR, we differentiate between traditional and 4IR patent exposure. Our method differs from previous approaches in that it both accounts for the diversity of task-level patent exposures within an occupation and reflects work activities more accurately. We find that exposure to 4IR patents differs from traditional patent exposure. Manual tasks, and accordingly occupations such as construction and production, are exposed mainly to traditional (non-4IR) patents but have low exposure to 4IR patents. The analysis suggests that 4IR technologies may have a negative impact on job growth; this impact appears 10 to 20 years after patent filing. Further, we compared the 4IR exposure to other automation and AI exposure scores. Whereas many measures refer to theoretical automation potential, our patent-based indicator reflects actual technology diffusion. Our work not only allows analyses of the impact of 4IR technologies as a whole, but also provides exposure scores for more than 300 technology fields, such as AI and smart office technologies. Finally, the work provides a general mapping of patents to tasks and occupations, which enables future researchers to construct individual exposure measures.


title: Generalizing the Lehmer's totient problem

abstract: An important unsolved question in number theory is the Lehmer's totient problem that asks whether there exists any composite number $n$ such that $\varphi(n)\mid n-1$, where $\varphi$ is the Euler's totient function. It is known that if any such $n$ exists, it must be odd, square-free, greater that $10^{30}$, and divisible by at least $15$ distinct primes. Such a number must be also a Carmichael number. In this short note, we discuss a group-theoretical analogous problem involving the function that counts the number of automorphisms of a finite group. Another way to generalize the Lehmer's totient problem is also proposed.


title: Did Goryachev et al. detect megahertz gravitational waves?

abstract: Goryachev et al. [1] recently announced the detection of "two strongly significant events" in their Bulk Acoustic Wave High Frequency Gravitational Wave Antenna. They claim many possibilities for the cause of these events, including high-frequency megahertz gravitational waves. We demonstrate these events are not due to gravitational waves for two reasons. 1) The inferred stochastic gravitational-wave background from these events implies the gravitational-wave energy density of the Universe is $\Omega_{\rm gw}\approx 10^8$, approximately $10^8$ times the closure density of the Universe. 2) The low-frequency gravitational-wave memory signal that accompanies any high-frequency gravitational-wave source visible by the current generation of high-frequency detectors would have been visible by LIGO/Virgo as a transient burst with signal-to-noise ratio $\gtrsim10^6$. The non-detection of such loud memory bursts throughout the operation of LIGO/Virgo rules out the gravitational-wave explanation for the high-frequency events detected by Goryachev et al. We discuss broader implications of this work for the ongoing experimental search for ultra high-frequency (MHz-GHz) gravitational waves.


title: On a group-theoretical generalization of the Euler's totient function

abstract: Let $G$ be a finite group and $\varphi(G)=|\{a\in G \mid o(a)=\exp(G)\}|$, where $o(a)$ denotes the order of $a$ in $G$ and $\exp(G)$ denotes the exponent of $G$. Under a natural hypothesis, in this note we determine the groups $G$ such that $\forall\, H,K\leq G$, $H\subseteq K$ implies $\varphi(H)\mid\varphi(K)$. This partially answers Problem 5.4 in \cite{6}.


title: Formation of graphene nanoribbons on the macrofacets of vicinal 6H-SiC(0001) surfaces

abstract: Thermal decomposition of vicinal 6H-SiC(0001) surfaces with off-angles toward the $[1\bar{1}00]$ direction results in the appearance of pairs of (0001) macroterraces and $(1\bar{1}0n)$ macrofacets covered with graphene, as follows. A carpet-like carbon layer grows on the surface, covering both the macroterraces and macrofacets; it forms $(6\sqrt{3} \times 6\sqrt{3})$ buffer layer on the former ones, whereas its partial periodic bonding with the SiC steps on the latter ones generates a pseudo-graphene nanoribbon (pseudo-GNR) array. The nanoribbons have a width of 2 nm and are aligned in the $[1\bar{2}10]$ direction with a spatial periodicity of 3.3 nm. Here, the Raman spectroscopy analysis of the pseudo-GNR array showed the absence of the 2D peak and the polarization dependence of the $G$ and $D$ peaks, which is typical of the armchair edge nanoribbon.


title: Submegahertz spectral width photon pair source based on fused silica microspheres

abstract: High efficiency, sub-MHz bandwidth photon pair generators will enable the field of quantum technology to transition from laboratory demonstrations to transformational applications involving information transfer from photons to atoms. While spontaneous parametric processes are able to achieve high efficiency photon pair generation, the spectral bandwidth tends to be relatively large, as defined by phase-matching constraints. To solve this fundamental limitation, we use an ultra-high quality factor (Q) fused silica microsphere resonant cavity to form a photon pair generator. We present the full theory for the SWM process in these devices, fully taking into account all relevant source characteristics in our experiments. The exceptionally narrow (down to kHz-scale) linewidths of these devices in combination with the device size results in a reduction in the bandwidth of the photon pair generation, allowing sub-MHz spectral bandwidth to be achieved. Specifically, using a pump source centered around 1550nm, photon pairs with the signal and idler modes at wavelengths close to 1540nm and 1560nm, respectively, are demonstrated. We herald a single idler-mode photon by detecting the corresponding signal photon, filtered via transmission through a wavelength division multiplexing channel of choice. We demonstrate the extraction of the spectral profile of a single peak in the single-photon frequency comb from a measurement of the signal-idler time of emission distribution. These improvements in device design and experimental methods enabled the narrowest spectral width (366kHz) to date in a heralded single photon source basedon SFWM.


title: Deep Learning Tools for Audacity: Helping Researchers Expand the Artist's Toolkit

abstract: We present a software framework that integrates neural networks into the popular open-source audio editing software, Audacity, with a minimal amount of developer effort. In this paper, we showcase some example use cases for both end-users and neural network developers. We hope that this work fosters a new level of interactivity between deep learning practitioners and end-users.


title: Sampling Multiple Nodes in Large Networks: Beyond Random Walks

abstract: Sampling random nodes is a fundamental algorithmic primitive in the analysis of massive networks, with many modern graph mining algorithms critically relying on it. We consider the task of generating a large collection of random nodes in the network assuming limited query access (where querying a node reveals its set of neighbors). In current approaches, based on long random walks, the number of queries per sample scales linearly with the mixing time of the network, which can be prohibitive for large real-world networks. We propose a new method for sampling multiple nodes that bypasses the dependence in the mixing time by explicitly searching for less accessible components in the network. We test our approach on a variety of real-world and synthetic networks with up to tens of millions of nodes, demonstrating a query complexity improvement of up to $\times 20$ compared to the state of the art.


title: Construction of a set of circulant matrix submatrices for faster MDS matrix verification

abstract: The present paper focuses on the construction of a set of submatrices of a circulant matrix such that it is a smaller set to verify that the circulant matrix is an MDS (maximum distance separable) one, comparing to the complete set of square submatrices needed in general case. The general MDS verification method requires to test for singular submatrices: if at least one square submatrix is singular the matrix is not MDS. However, the complexity of the general method dramatically increases for matrices of a greater dimension. We develop an algorithm that constructs a smaller subset of submatrices thanks to a simple structure of circulant matrices. The algorithm proposed in the paper reduces the size of the testing set by approximately two matrix orders.


title: Highly nonlinear dipolar exciton-polaritons in bilayer MoS$_2$

abstract: Realizing nonlinear optical response in the low photon density limit in solid-state systems has been a long-standing challenge. Semiconductor microcavities in the strong coupling regime hosting exciton-polaritons have emerged as attractive candidates in this context. However, the weak interaction between these quasiparticles has been a hurdle in this quest. Dipolar excitons provide an attractive strategy to overcome this limitation but are often hindered by their weak oscillator strength. The interlayer dipolar excitons in naturally occurring homobilayer MoS$_2$ alleviates this issue owing to their hybridization of interlayer charge transfer exciton and intralayer B exciton. Here we demonstrate the formation of dipolar exciton polaritons in bilayer MoS$_2$ resulting in unprecedented nonlinear interaction strengths. A ten-fold increase in nonlinearity is observed for the interlayer dipolar excitons compared to the conventional A excitons. These highly nonlinear dipolar polaritons will likely be a frontrunner in the quest for solid-state quantum nonlinear devices.


title: Enhancing Ground-Based Observations of Trans-Neptunian Objects using a Single-Epoch Parallax Measurement from L2

abstract: Space-based observatories at the second Sun-Earth Lagrange Point (L2) offer a unique opportunity to efficiently determine the orbits of distant solar system objects by taking advantage of the parallax effect that arises from nearly simultaneous ground- and space-based observations. Given the typical orbit of an observatory about L2, the observational baseline of $\sim1.5\times10^6$ km between L2 and Earth results in an instantaneous parallax of $\sim$10''-100'', even for the most distant of detectable trans-Neptunian objects (TNOs) in our solar system. Current ground-based strategies for measuring the orbits of TNOs are very expensive and require multiple years of observations. We show that the direct constraint on the distance to a TNO, afforded by near-simultaneous ground- and space-based observations, allows us to confidently determine orbits with as few as three ground-based observations spanning 24 hr combined with a single observational epoch from an observatory orbiting L2.


title: Eigenvalue Bounds for Double Saddle-Point Systems

abstract: We derive bounds on the eigenvalues of a generic form of double saddle-point matrices. The bounds are expressed in terms of extremal eigenvalues and singular values of the associated block matrices. Inertia and algebraic multiplicity of eigenvalues are considered as well. The analysis includes bounds for preconditioned matrices based on block diagonal preconditioners using Schur complements, and it is shown that in this case the eigenvalues are clustered within a few intervals bounded away from zero. Analysis for approximations of Schur complements is included. Some numerical experiments validate our analytical findings.


title: From SKA to SKAO: Early Progress in the SKAO Construction

abstract: The Square Kilometre Array telescopes have recently started their construction phase, after years of pre-construction effort. The new SKA Observatory (SKAO) intergovernmental organisation has been created, and the start of construction ($\mathrm{T_0}$) has already happened. In this talk, we summarise the construction progress of our facility, and the role that agile software development and open-source collaboration, and in particular the development of our TANGO-based control system, is playing.


title: Recipes for when Physics Fails: Recovering Robust Learning of Physics Informed Neural Networks

abstract: Physics-informed Neural Networks (PINNs) have been shown to be effective in solving partial differential equations by capturing the physics induced constraints as a part of the training loss function. This paper shows that a PINN can be sensitive to errors in training data and overfit itself in dynamically propagating these errors over the domain of the solution of the PDE. It also shows how physical regularizations based on continuity criteria and conservation laws fail to address this issue and rather introduce problems of their own causing the deep network to converge to a physics-obeying local minimum instead of the global minimum. We introduce Gaussian Process (GP) based smoothing that recovers the performance of a PINN and promises a robust architecture against noise/errors in measurements. Additionally, we illustrate an inexpensive method of quantifying the evolution of uncertainty based on the variance estimation of GPs on boundary data. Robust PINN performance is also shown to be achievable by choice of sparse sets of inducing points based on sparsely induced GPs. We demonstrate the performance of our proposed methods and compare the results from existing benchmark models in literature for time-dependent Schr\"odinger and Burgers' equations.