title: Speckle interferometry at SOAR in 2019

abstract: The results of speckle interferometric observations at the 4.1 m Southern Astrophysical Research Telescope (SOAR) in 2019 are given, totaling 2555 measurements of 1972 resolved pairs with separations from 15 mas (median 0.21") and magnitude difference up to 6 mag, and non-resolutions of 684 targets. We resolved for the first time 90 new pairs or subsystems in known binaries. This work continues our long-term speckle program. Its main goal is to monitor orbital motion of close binaries, including members of high-order hierarchies and Hipparcos pairs in the solar neighborhood. We give a list of 127 orbits computed using our latest measurements. Their quality varies from excellent (25 orbits of grades 1 and 2) to provisional (47 orbits of grades 4 and 5).


title: Second largest Eigenpair Statistics for Sparse Graphs

abstract: We develop a formalism to compute the statistics of the second largest eigenpair of weighted sparse graphs with $N\gg 1$ nodes, finite mean connectivity and bounded maximal degree, in cases where the top eigenpair statistics is known. The problem can be cast in terms of optimisation of a quadratic form on the sphere with a fictitious temperature, after a suitable deflation of the original matrix model. We use the cavity and replica methods to find the solution in terms of self-consistent equations for auxiliary probability density functions, which can be solved by an improved population dynamics algorithm enforcing eigenvector orthogonality on-the-fly. The analytical results are in perfect agreement with numerical diagonalisation of large (weighted) adjacency matrices, focussing on the cases of random regular and Erd\H{o}s-R\'enyi graphs. We further analyse the case of sparse Markov transition matrices for unbiased random walks, whose second largest eigenpair describes the non-equilibrium mode with the largest relaxation time. We also show that the population dynamics algorithm with population size $N_P$ does not actually capture the thermodynamic limit $N\to\infty$ as commonly assumed: the accuracy of the population dynamics algorithm has a strongly non-monotonic behaviour as a function of $N_P$, thus implying that an optimal size $N_P^\star=N_P^\star(N)$ must be chosen to best reproduce the results from numerical diagonalisation of graphs of finite size $N$.


title: Molecular Sparse Representation by 3D Ellipsoid Radial Basis Function Neural Networks via $L_1$ Regularization

abstract: In this paper, we have developed an ellipsoid radial basis function neural network (ERBFNN) and algorithm for sparse representing of a molecular shape. To evaluate a sparse representation of the molecular shape model, the Gaussian density map of molecule is approximated by ERBFNN with a relatively small number of neurons. The deep learning models were trained by optimizing a nonlinear loss function with $L_1$ regularization. Experimental results demonstrate that the original molecular shape is able to be represented with good accuracy by much fewer scale of ERBFNN by our algorithm. And our network in principle can be applied to multi-resolution sparse representation of molecular shape and coarse-grained molecular modeling.


title: Lattice-based public key encryption with equality test supporting flexible authorization in standard model

abstract: Public key encryption with equality test (PKEET) supports to check whether two ciphertexts encrypted under different public keys contain the same message or not. PKEET has many interesting applications such as keyword search on encrypted data, encrypted data partitioning for efficient encrypted data management, personal health record systems, spam filtering in encrypted email systems and so on. However, the PKEET scheme lacks an authorization mechanism for a user to control the comparison of its ciphertexts with others. In 2015, Ma et al. introduce the notion of PKEET with flexible authorization (PKEET-FA) which strengthens privacy protection. Since 2015, there are several follow-up works on PKEET-FA. But, all are secure in the random-oracle model. Moreover, all are vulnerable to quantum attacks. In this paper, we provide three constructions of quantum-safe PKEET-FA secure in the standard model. Proposed constructions are secure based on the hardness assumptions of integer lattices and ideal lattices. Finally, we implement the PKEET-FA scheme over ideal lattices.


title: Viscosity Solutions to Second Order Path-Dependent Hamilton-Jacobi-Bellman Equations and Applications

abstract: In this article, a notion of viscosity solutions is introduced for second order path-dependent Hamilton-Jacobi-Bellman (PHJB) equations associated with optimal control problems for path-dependent stochastic differential equations. We identify the value functional of optimal control problems as unique viscosity solution to the associated PHJB equations. We also show that our notion of viscosity solutions is consistent with the corresponding notion of classical solutions, and satisfies a stability property. Applications to backward stochastic Hamilton-Jacobi-Bellman equations are also given.


title: Rational Finance Approach to Behavioral Option Pricing

abstract: When pricing options, there may be different views on the instantaneous mean return of the underlying price process. According to Black (1972), where there exist heterogeneous views on the instantaneous mean return, this will result in arbitrage opportunities. Behavioral finance proponents argue that such heterogenous views are likely to occur and this will not impact option pricing models proposed by rational dynamic asset pricing theory and will not give rise to volatility smiles. To rectify this, a leading advocate of behavioral finance has proposed a behavioral option pricing model. As there may be unexplored links between the behavioral and rational approaches to option pricing, in this paper we revisit Shefrin (2008) option pricing model as an example and suggest one approach to modify this behavioral finance option pricing formula to be consistent with rational dynamic asset pricing theory by introducing arbitrage transaction costs which offset the gains from arbitrage trades.


title: Completeness and injectivity

abstract: We show that for any quantale $\mathcal{Q}$, a $\mathcal{Q}$-category is skeletal and complete if and only if it is injective with respect to fully faithful $\mathcal{Q}$-functors. This is a special case of known theorems due to Hofmann and Stubbe, but we provide a different proof, using the characterisation of the MacNeille completion of a $\mathcal{Q}$-category as its injective envelope. For Lawvere metric spaces, our results yield those of Kemajou, K\"unzi and Otafudu. We point out that their notion of Isbell convexity can be seen as a geometric formulation of categorical completeness for Lawvere metric spaces.


title: Null $\varphi $-Slant Curves in a Main Class of 3-Dimensional Normal Almost Contact B-Metric Manifolds

abstract: We introduce a new type of slant curves in almost contact B-metric manifolds, called $\varphi $-slant curves, by an additional condition which is specific for these manifolds. In this paper we study $\varphi $-slant null curves in a class of 3-dimensional normal almost contact B-metric manifolds and prove that for non-geodesic of them there exists a unique Frenet frame for which the original parameter is distinguished. We investigate some of $\varphi $-slant null curves and with respect to the associated B-metric on the manifold and find relationships between the corresponding Frenet frames and curvatures. We construct the examined curves in a 3-dimensional Lie group and give their matrix representation.


title: Audio and Contact Microphones for Cough Detection

abstract: In the framework of assessing the pathology severity in chronic cough diseases, medical literature underlines the lack of tools for allowing the automatic, objective and reliable detection of cough events. This paper describes a system based on two microphones which we developed for this purpose. The proposed approach relies on a large variety of audio descriptors, an efficient algorithm of feature selection based on their mutual information and the use of artificial neural networks. First, the possible use of a contact microphone (placed on the patient's thorax or trachea) in complement to the audio signal is investigated. This study underlines that this contact microphone suffers from reliability issues, and conveys little new relevant information compared to the audio modality. Secondly, the proposed audio-only approach is compared to a commercially available system using four sensors on a database with different sound categories often misdetected as coughs, and produced in various conditions. With average sensitivity and specificity of 94.7% and 95% respectively, the proposed method achieves better cough detection performance than the commercial system.


title: A class of Integral Operators from Lebesgue spaces into Harmonic Bergman-Besov or Weighted Bloch Spaces

abstract: We consider a class of two-parameter weighted integral operators induced by harmonic Bergman-Besov kernels on the unit ball of $\mathbb{R}^{n}$ and characterize precisely those that are bounded from Lebesgue spaces $L^{p}_{\alpha}$ into Harmonic Bergman-Besov $b^{q}_{\beta}$ or weighted Bloch Spaces $b^{\infty}_{\beta} $, for $1\leq p\leq\infty$, $1\leq q< \infty$ and $\alpha,\beta \in \mathbb{R}$. These operators can be viewed as generalizations of the harmonic Bergman-Besov projections. Also, our results remove the disturbing conditions $\beta>-1$ when $q<\infty$ and $\beta\geq 0$ when $q=\infty$ of Do\u{g}an (A Class of Integral Operators Induced by Harmonic Bergman-Besov kernels on Lebesgue Classes, preprint, 2020) by mapping the operators into these spaces instead of the Lebesgue classes.


title: Polynomial Equations in Subgroups and Applications

abstract: We obtain a new bound for the number of solutions to polynomial equations in cosets of multiplicative subgroups in finite fields, which generalises previous results of P. Corvaja and U. Zannier (2013). We also obtain a conditional improvement of recent results of J. Bourgain, A. Gamburd and P. Sarnak (2016) and S. V. Konyagin, S. V. Makarychev, I. E. Shparlinski and I. V. Vyugin (2019) on the structure of solutions to the reduction of the Markoff equation $x^2 + y^2 + z^2 = 3 x yz$ modulo a prime $p$.


title: Radio frequency emissions from dark-matter-candidate magnetized quark nuggets interacting with matter

abstract: Quark nuggets are theoretical objects composed of approximately equal numbers of up, down, and strange quarks. They are also called strangelets, nuclearites, AQNs, slets, Macros, and MQNs. Quark nuggets are a candidate for dark matter, which has been a mystery for decades despite constituting ${\sim}85\%$ of matter. Most previous models of quark nuggets have assumed no intrinsic magnetic field; however, Tatsumi found that quark nuggets may exist in magnetars as a ferromagnetic liquid with a magnetic field between $10^{ 11}$ T and $10^{ 13}$ T. We apply that result to quark nuggets, a dark-matter candidate consistent with the Standard Model, and report results of analytic calculations and simulations that show they spin up and emit electromagnetic radiation at ${\sim}10^{ 4}$ Hz to ${\sim}10^{ 9}$ Hz after passage through planetary environments. The results depend strongly on the value of $B_{o}$, which is a parameter to guide and interpret observations. A proposed sensor system with three satellites at 51,000 km altitude illustrates the feasibility of using radio-frequency emissions to detect 0.003 to 1600 MQNs, depending on Bo, during a five year mission.


title: Second gradient electrodynamics: Green functions, wave propagation, regularization and self-force

abstract: In this work, the theory of second gradient electrodynamics, which is an important example of generalized electrodynamics, is proposed and investigated. Second gradient electrodynamics is a gradient field theory with up to second-order derivatives of the electromagnetic field strengths in the Lagrangian density. Second gradient electrodynamics possesses a weak nonlocality in space and time. In the framework of second gradient electrodynamics, the retarded Green functions, first-order derivatives of the retarded Green functions, retarded potentials, retarded electromagnetic field strengths, generalized Lienard-Wiechert potentials and the corresponding electromagnetic field strengths are derived for three, two and one spatial dimensions. The behaviour of the electromagnetic fields is investigated on the light cone. In particular, the retarded Green functions and their first-order derivatives show oscillations around the classical solutions inside the forward light cone and it is shown that they are singularity-free and regular on the light cone in three, two and one spatial dimensions. In second gradient electrodynamics, the self-force and the energy release rate are calculated and the equation of motion of a charged point particle, which is an integro-differential equation where the infamous third-order time-derivative of the position does not appear, is determined.


title: Taking advantage of multiplet structure for lineshape analysis in Fourier space

abstract: Lineshape analysis is a recurrent and often computationally intensive task in optics, even more so for multiple peaks in the presence of noise. We demonstrate an algorithm which takes advantage of peak multiplicity (N) to retrieve line shape information. The method is exemplified via analysis of Lorentzian and Gaussian contributions to individual lineshapes for a practical spectroscopic measurement and benefits from a linear increase in sensitivity with the number N. The robustness of the method and its benefits in terms of noise reduction and order of magnitude improvement in run-time performance are discussed.


title: Hardware Implementation of Adaptive Watermarking Based on Local Spatial Disorder Analysis

abstract: With the increasing use of the internet and the ease of exchange of multimedia content, the protection of ownership rights has become a significant concern. Watermarking is an efficient means for this purpose. In many applications, real-time watermarking is required, which demands hardware implementation of low complexity and robust algorithm. In this paper, an adaptive watermarking is presented, which uses embedding in different bit-planes to achieve transparency and robustness. Local disorder of pixels is analyzed to control the strength of the watermark. A new low complexity method for disorder analysis is proposed, and its hardware implantation is presented. An embedding method is proposed, which causes lower degradation in the watermarked image. Also, the performance of proposed watermarking architecture is improved by a pipe-line structure and is tested on an FPGA device. Results show that the algorithm produces transparent and robust watermarked images. The synthesis report from FPGA implementation illustrates a low complexity hardware structure.


title: Simulating the spread of COVID-19 via spatially-resolved susceptible-exposed-infected-recovered-deceased (SEIRD) model with heterogeneous diffusion

abstract: We present an early version of a Susceptible-Exposed-Infected-Recovered-Deceased (SEIRD) mathematical model based on partial differential equations coupled with a heterogeneous diffusion model. The model describes the spatio-temporal spread of the COVID-19 pandemic, and aims to capture dynamics also based on human habits and geographical features. To test the model, we compare the outputs generated by a finite-element solver with measured data over the Italian region of Lombardy, which has been heavily impacted by this crisis between February and April 2020. Our results show a strong qualitative agreement between the simulated forecast of the spatio-temporal COVID-19 spread in Lombardy and epidemiological data collected at the municipality level. Additional simulations exploring alternative scenarios for the relaxation of lockdown restrictions suggest that reopening strategies should account for local population densities and the specific dynamics of the contagion. Thus, we argue that data-driven simulations of our model could ultimately inform health authorities to design effective pandemic-arresting measures and anticipate the geographical allocation of crucial medical resources.


title: Channel-Aware Adversarial Attacks Against Deep Learning-Based Wireless Signal Classifiers

abstract: This paper presents channel-aware adversarial attacks against deep learning-based wireless signal classifiers. There is a transmitter that transmits signals with different modulation types. A deep neural network is used at each receiver to classify its over-the-air received signals to modulation types. In the meantime, an adversary transmits an adversarial perturbation (subject to a power budget) to fool receivers into making errors in classifying signals that are received as superpositions of transmitted signals and adversarial perturbations. First, these evasion attacks are shown to fail when channels are not considered in designing adversarial perturbations. Then, realistic attacks are presented by considering channel effects from the adversary to each receiver. After showing that a channel-aware attack is selective (i.e., it affects only the receiver whose channel is considered in the perturbation design), a broadcast adversarial attack is presented by crafting a common adversarial perturbation to simultaneously fool classifiers at different receivers. The major vulnerability of modulation classifiers to over-the-air adversarial attacks is shown by accounting for different levels of information available about the channel, the transmitter input, and the classifier model. Finally, a certified defense based on randomized smoothing that augments training data with noise is introduced to make the modulation classifier robust to adversarial perturbations.


title: Verlinde gravity effects on the orbits of the planets and the Moon in the Solar System

abstract: In this work, we address the effects of a phenomenon known as Verlinde gravity. Here we show that its effect over the planets and the Moon in our solar system is quite negligible. We find that the Verlinde gravity effects on the orbits of planets are at least 10 times smaller than the precision with which we can determine the Sun's mass, and the one on the orbit of the Moon is about 100 times smaller than the precision with which we can determine the Earth's mass. These results let us infer that statements in the literature that Verlinde gravity is ruled out by the observed motion of planets in our solar system aren't correct.


title: Rate Adaptation in Predictor Antenna Systems

abstract: Predictor antenna (PA) system is referred to as a system with two sets of antennas on the roof of a vehicle, where the PAs positioned in the front of the vehicle are used to predict the channel state observed by the receive antennas (RAs) that are aligned behind the PAs. This letter studies the performance of PA systems in the presence of the mismatching problem, i.e., when the channel observed by the PA is not exactly the same as the channel experienced by the RA. Particularly, we study the effect of spatial mismatching on the accuracy of channel state information estimation and rate adaption. We derive closed-form expressions for instantaneous throughput, outage probability, and the throughput-optimized rate adaptation. Also, we take the temporal evolution of the channel into account and evaluate the system performance in temporally-correlated conditions. The simulation and analytical results show that, while PA-assisted adaptive rate adaptation leads to considerable performance improvement, the throughput and the outage probability are remarkably affected by the spatial mismatch and temporal correlations.


title: Infinite mixtures of multivariate normal-inverse Gaussian distributions for clustering of skewed data

abstract: Mixtures of multivariate normal inverse Gaussian (MNIG) distributions can be used to cluster data that exhibit features such as skewness and heavy tails. However, for cluster analysis, using a traditional finite mixture model framework, either the number of components needs to be known $a$-$priori$ or needs to be estimated $a$-$posteriori$ using some model selection criterion after deriving results for a range of possible number of components. However, different model selection criteria can sometimes result in different number of components yielding uncertainty. Here, an infinite mixture model framework, also known as Dirichlet process mixture model, is proposed for the mixtures of MNIG distributions. This Dirichlet process mixture model approach allows the number of components to grow or decay freely from 1 to $\infty$ (in practice from 1 to $N$) and the number of components is inferred along with the parameter estimates in a Bayesian framework thus alleviating the need for model selection criteria. We provide real data applications with benchmark datasets as well as a small simulation experiment to compare with other existing models. The proposed method provides competitive clustering results to other clustering approaches for both simulation and real data and parameter recovery are illustrated using simulation studies.


title: A Relational Gradient Descent Algorithm For Support Vector Machine Training

abstract: We consider gradient descent like algorithms for Support Vector Machine (SVM) training when the data is in relational form. The gradient of the SVM objective can not be efficiently computed by known techniques as it suffers from the ``subtraction problem''. We first show that the subtraction problem can not be surmounted by showing that computing any constant approximation of the gradient of the SVM objective function is $\#P$-hard, even for acyclic joins. We, however, circumvent the subtraction problem by restricting our attention to stable instances, which intuitively are instances where a nearly optimal solution remains nearly optimal if the points are perturbed slightly. We give an efficient algorithm that computes a ``pseudo-gradient'' that guarantees convergence for stable instances at a rate comparable to that achieved by using the actual gradient. We believe that our results suggest that this sort of stability the analysis would likely yield useful insight in the context of designing algorithms on relational data for other learning problems in which the subtraction problem arises.


title: Orbital precession in the distant solar system; further constraining the Planet Nine hypothesis with numerical simulations

abstract: The longitudes of perihelia and orbital poles of the solar system's dozen or so most remote detected objects are clustered in a manner inconsistent with that of a random sample of uniformly distributed orbits. While small number statistics and observational biases may explain these features, the statistical significance of the clustering has led to the recent development of the "Planet Nine hypothesis." In the proposed scenario, orbits in the distant solar system are shepherded via secular perturbations from an undetected massive planet on an eccentric orbit. However, the precession of perihelia and nodes in the outer Kuiper Belt and inner Oort Cloud are also affected by the the giant planets, passing stars, and the galactic tide. We perform a large suite of numerical simulations designed to study the orbital alignment of Extreme Trans-Neptunian Objects (ETNOs) and Inner Oort Cloud Objects (IOCOs). In our various integrations that include Planet Nine, we consistently find that >60% of ETNOs and IOCOs that are detectable after 4 Gyr are also anti-aligned in perihelia with the distant massive perturber. However, when we randomly select 17 objects from this sample of remaining orbits, there is significant scatter in the degree of longitude of perihelion and orbital pole clustering that might be observed. Furthermore, we argue that, in the absence of Planet Nine, 17 randomly drawn orbits should still exhibit some clustering even if the underlying distribution is uniform. Thus, we find that still more ETNO and IOCO detections are required to confidently infer the presence of Planet Nine.


title: Inevitable imprints of patchy reionization on the cosmic microwave background anisotropy

abstract: Reionization of the cosmic neutral hydrogen by the first stars in the Universe is an inhomogeneous process which produces spatial fluctuations in free electron density. These fluctuations lead to observable signatures in cosmological probes like the cosmic microwave background (CMB). We explore the effect of the electron density fluctuations on CMB using photon-conserving semi-numerical simulations of reionization named \texttt{SCRIPT}. We show that the amplitude of the kinematic Sunyaev-Zeldovich (kSZ) and $B$-mode polarization signal depends on the patchiness in the spatial distribution of electrons along with the dependence on mid-point and extent of the reionization history. Motivated by this finding, we provide new scaling relations for the amplitude of kSZ and $B$-mode polarization signal which can capture the effects arising from the mean optical depth, width of reionization, and spatial fluctuations in the electron density arising from patchy reionization. We show that the amplitude of the kSZ and $B$-mode polarization signal exhibits different dependency on the width of reionization and the patchiness of reionization, and hence a joint study of these CMB probes will be able to break the degeneracy. By combining external datasets from 21~cm measurements, the degeneracy can be further lifted by directly exploring the sizes of the ionized regions.


title: Dynamic Zoom Simulations: a fast, adaptive algorithm for simulating lightcones

abstract: The advent of a new generation of large-scale galaxy surveys is pushing cosmological numerical simulations in an uncharted territory. The simultaneous requirements of high resolution and very large volume pose serious technical challenges, due to their computational and data storage demand. In this paper, we present a novel approach dubbed Dynamic Zoom Simulations -- or DZS -- developed to tackle these issues. Our method is tailored to the production of lightcone outputs from N-body numerical simulations, which allow for a more efficient storage and post-processing compared to standard comoving snapshots, and more directly mimic the format of survey data. In DZS, the resolution of the simulation is dynamically decreased outside the lightcone surface, reducing the computational work load, while simultaneously preserving the accuracy inside the lightcone and the large-scale gravitational field. We show that our approach can achieve virtually identical results to traditional simulations at half of the computational cost for our largest box. We also forecast this speedup to increase up to a factor of 5 for larger and/or higher-resolution simulations. We assess the accuracy of the numerical integration by comparing pairs of identical simulations run with and without DZS. Deviations in the lightcone halo mass function, in the sky-projected lightcone, and in the 3D matter lightcone always remain below 0.1%. In summary, our results indicate that the DZS technique may provide a highly-valuable tool to address the technical challenges that will characterise the next generation of large-scale cosmological simulations.


title: Loop amplitudes monodromy relations and color-kinematics duality

abstract: Color-kinematics duality is a remarkable conjectured property of gauge theory which, together with double copy, is at the heart of a wealth of new developments in scattering amplitudes. So far, its validity has been verified in most cases only empirically, lacking a first-principle understanding beyond tree-level. In this paper we provide initial steps in a first-principle understanding of color-kinematics duality and double-copy at loop level, through a detailed analysis of the field-theory limit of the monodromy relations of string theory at one loop. In this limit, we dissect the type of Feynman graphs generated and the relations they obey. We find that graphs with contact-terms are unavoidable and are generated in the field theory limit of "bulk" contours which do not have a standard physical interpretation in string perturbation theory. We show how they are related to ambiguities in the definition of the loop momentum and that their role is precisely to cancel those ambiguities.


title: Bandlimited Entanglement Harvesting

abstract: There are many reasons to believe that there is a fundamental minimum length scale below which distances cannot be reliably resolved. One method of constructing a quantum field with a finite minimum length scale is to use bandlimited quantum field theory, where the spacetime is mathematically both continuous and discrete. This is a modification to the field, which has been shown to have many consequences at the level of the field. We consider an operational approach and use a pair of particle detectors (two-level qubits) as a local probe of the field, which are coupled to the vacuum of the bandlimited massless scalar field in a time dependent way through a switching function. We show that, mathematically, the bandlimit modifies the spatial profile of the detectors so that they are only quasi-local. We explore two different types of switching functions, Gaussian and Dirac delta. We find that with Gaussian switching, the bandlimit exponentially suppresses the de-excitation of the detectors when the energy gap between the two levels is larger than the bandlimit. If the detectors are prepared in ground state, in certain regions of the parameter space they are able to extract more entanglement from the field than if there was no bandlimit. When the detectors couple with Dirac-delta switching, we show that a particle detector is most sensitive to the bandlimit when it couples to a small but finite region of spacetime. We find that the effects of a bandlimit are detectable using local probes. This work is important because it illustrates the possible observable consequences of a fundamental bandlimit in a quantum field.


title: The $n$-point streaming model: how velocities shape correlation functions in redshift space

abstract: Starting from first principles, we derive the fundamental equations that relate the $n$-point correlation functions in real and redshift space. Our result generalises the so-called `streaming model' to higher-order statistics: the full $n$-point correlation in redshift-space is obtained as an integral of its real-space counterpart times the joint probability density of $n-1$ relative line-of-sight peculiar velocities. Equations for the connected $n$-point correlation functions are obtained by recursively applying the generalised streaming model for decreasing $n$. Our results are exact within the distant-observer approximation and completely independent of the nature of the tracers for which the correlations are evaluated. Focusing on 3-point statistics, we use an $N$-body simulation to study the joint probability density function of the relative line-of-sight velocities of pairs of particles in a triplet. On large scales, we find that this distribution is approximately Gaussian and that its moments can be accurately computed with standard perturbation theory. We use this information to formulate a phenomenological 3-point Gaussian streaming model. A practical implementation is obtained by using perturbation theory at leading order to approximate several statistics in real space. In spite of this simplification, the resulting predictions for the matter 3-point correlation function in redshift space are in rather good agreement with measurements performed in the simulation. We discuss the limitations of the simplified model and suggest a number of possible improvements. Our results find direct applications in the analysis of galaxy clustering but also set the basis for studying 3-point statistics with future peculiar-velocity surveys and experiments based on the kinetic Sunyaev-Zel'dovich effect.


title: Resonant Neutrino Self-Interactions

abstract: If neutrinos have self-interactions, these will induce scatterings between astrophysical and cosmic neutrinos. Prior work proposed to look for possible resulting resonance features in astrophysical neutrino spectra in order to seek a neutrino self-interaction which can be either diagonal in the neutrino flavor space or couple different neutrino flavors. The calculation of the astrophysical spectra involves either a Monte Carlo simulation or a computationally intensive numerical integration of an integro-partial-differential equation. As a result only limited regions of the neutrino self-interaction parameter space have been explored, and only flavor-diagonal self-interactions have been considered. Here, we present a fully analytic form for the astrophysical neutrino spectra for arbitrary neutrino number and arbitrary self-coupling matrix that accurately obtains the resonance features in the observable neutrino spectra. The results can be applied to calculations of the diffuse supernova neutrino background and of the spectrum from high-energy astrophysical neutrino sources. We illustrate with a few examples.


title: On the Paradoxical Impact of Blending by Red Clump Giants

abstract: The impact of blending by RCGs (red clump giants, or relatively metal-rich red horizontal branch stars) is discussed as it relates to RRab and classical Cepheids, and invariably establishing an improved distance scale. An analysis of OGLE Magellanic Cloud variables reaffirms that blending with RCGs may advantageously thrust remote extragalactic stars into the range of detectability. Specifically, simulations of Magellanic Cloud RRab and RCG blends partly reproduce bright non-canonical trends readily observed in amplitude-magnitude space ($I_c$ vs.~$A_{I_c}$). Conversely, the larger magnitude offset between classical Cepheids and RCGs causes the latter's influence to be challenging to address. The relative invariance of a Wesenheit function's slope to metallicity (e.g., $W_{VI_c}$) implies that a deviation from the trend could reveal blending and photometric inaccuracies (e.g., standardization), as blending by RCGs (a proxy of an evolved red stellar demographic) can flatten period-Wesenehit relations owing to the increased impact on less-luminous shorter-period Cepheids. That could partly explain both a shallower inferred Wesenheit function and over-estimated $H_0$ values. A consensus framework to identify and exploit blending is desirable, as presently $H_0$ estimates from diverse teams are unwittingly leveraged without homogenizing the disparate approaches (e.g., no blending correction to a sizable $\simeq 0.^{\rm m}3$).


title: Getting High: High Fidelity Simulation of High Granularity Calorimeters with High Speed

abstract: Accurate simulation of physical processes is crucial for the success of modern particle physics. However, simulating the development and interaction of particle showers with calorimeter detectors is a time consuming process and drives the computing needs of large experiments at the LHC and future colliders. Recently, generative machine learning models based on deep neural networks have shown promise in speeding up this task by several orders of magnitude. We investigate the use of a new architecture -- the Bounded Information Bottleneck Autoencoder -- for modelling electromagnetic showers in the central region of the Silicon-Tungsten calorimeter of the proposed International Large Detector. Combined with a novel second post-processing network, this approach achieves an accurate simulation of differential distributions including for the first time the shape of the minimum-ionizing-particle peak compared to a full GEANT4 simulation for a high-granularity calorimeter with 27k simulated channels. The results are validated by comparing to established architectures. Our results further strengthen the case of using generative networks for fast simulation and demonstrate that physically relevant differential distributions can be described with high accuracy.

categories: physics.ins-det hep-ex hep-ph physics.data-an


title: Long-term evolution of a magnetic massive merger product

abstract: About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $\tau$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars.


title: Finite-temperature spectroscopy of dirty helical Luttinger liquids

abstract: We develop a theory of finite-temperature momentum-resolved tunneling spectroscopy (MRTS) for disordered, interacting two-dimensional topological-insulator edges. The MRTS complements conventional electrical transport measurement in characterizing the properties of the helical Luttinger liquid edges. Using standard bosonization technique, we study low-energy spectral function and the MRTS tunneling current, providing a detailed description controlled by disorder, interaction, and temperature, taking into account Rashba spin orbit coupling, interedge interaction and distinct edge velocities. Our theory provides a systematic description of the spectroscopic signals in the MRTS measurement and we hope will stimulate future experimental studies on the two-dimensional time-reversal invariant topological insulator.


title: Quantum dynamics of the corner of the Bianchi IX model in WKB approximation

abstract: In this paper we analyse the Bianchi IX Universe dynamics within the corner region associated to the potential term which the spatial curvature induces in the Minisuperspace. The study is done in the vacuum and in the presence of a massless scalar field $\phi$ and a cosmological constant term $\Lambda$. We investigate the dynamics in terms of WKB scenario for which the isotropic Misner variables (the volume) and one of the two anisotropic ones (and $\phi$ when present) are treated on a semi-classical level, while the remaining anisotropy degree of freedom, the one trapped in the corner, is described on a pure quantum level. The quantum dynamics always reduces to the one of a time-dependent Schr\"{o}edinger equation for a harmonic potential with a time dependent frequency. The vacuum case is treated in the limits of a collpasing and an expanding Universe, while the dynamics in presence of $\phi$ and $\Lambda$ is studied for $t \rightarrow \infty $. In both analysis the quantum dynamics of the anisotropy variable is associated to a decaying standard deviation of its probability density, corresponding to a suppression of the quantum anisotropy associated. In the vacuum case, the corner configuration becomes an attractor for the dynamics and the evolution resembles that one of a Taub cosmology in the limit of a non-singular initial Universe. This suggests that if the Bianchi dynamics enters enough the potential corner then the initial singularity is removed and a Taub picture emerges. The case when $\phi$ is present well mimics the De-Sitter phase of an inflationary Universe. Here we show that both the classical and quantum anisotropies are exponentially suppressed, so that the resulting dynamics corresponds to an isotropic closed Robertson-Walker geometry.


title: A hybrid post-Newtonian -- effective-one-body scheme for spin-precessing compact-binary waveforms

abstract: We introduce \texttt{TEOBResumSP}: an efficient, accurate hybrid scheme for generating gravitational waveforms from spin-precessing compact binaries. The precessing waveforms are generated via the established technique of Euler rotating the non-precessing \texttt{TEOBResumS} waveforms from a precessing frame to an inertial frame. We obtain the Euler angles by solving the post-Newtonian precession equations expanded to second post-Newtonian order. Current version of \texttt{TEOBResumSP} produces precessing waveforms through the inspiral phase up to the onset of the merger. We compare \texttt{TEOBResumSP} to current state-of-the-art precessing approximants \texttt{NRSur7dq4}, \texttt{SEOBNRv4PHM}, and \texttt{IMRPhenomPv3HM} for 200 cases of precessing compact binary inspirals with orbital inclinations up to 90 degrees, mass ratios up to four, and the effective precession parameter $\chi_p$ up to 0.75. We further provide an extended comparison with \texttt{SEOBNRv4PHM} involving 1030 more inspirals with $\chi_p\le 1$ and mass ratios up to 10. We find that 91\% of the \texttt{TEOBResumSP}-\texttt{NRSur7dq4} matches, 85\% of the \texttt{TEOBResumSP}-\texttt{SEOBNRv4PHM} matches, and 77\% of the \texttt{TEOBResumSP}-\texttt{IMRPhenomPv3HM} matches are greater than $0.965$. Most disagreements occur for large mass ratios and $\chi_p \gtrsim 0.6$. We identify the mismatch of the \emph{non}-precessing $(2,1)$ mode as one of the leading causes of disagreements. We also introduce a new parameter, $\chi_{\perp,\text{max}}$, to measure the strength of precession and hint that the mismatch between the above approximants shows an exponential dependence on $\chi_{\perp,\text{max}}$ though this requires further study. Our results indicate that \texttt{TEOBResumSP} is on its way to becoming a robust precessing approximant to be employed in the parameter estimation of generic-spin compact binaries.


title: Enabling Language Models to Fill in the Blanks

abstract: We present a simple approach for text infilling, the task of predicting missing spans of text at any position in a document. While infilling could enable rich functionality especially for writing assistance tools, more attention has been devoted to language modeling---a special case of infilling where text is predicted at the end of a document. In this paper, we aim to extend the capabilities of language models (LMs) to the more general task of infilling. To this end, we train (or fine-tune) off-the-shelf LMs on sequences containing the concatenation of artificially-masked text and the text which was masked. We show that this approach, which we call infilling by language modeling, can enable LMs to infill entire sentences effectively on three different domains: short stories, scientific abstracts, and lyrics. Furthermore, we show that humans have difficulty identifying sentences infilled by our approach as machine-generated in the domain of short stories.


title: A tidal disruption event coincident with a high-energy neutrino

abstract: Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently reported the likely association of one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the likely association of a radio-emitting tidal disruption event, AT2019dsg, with a second high-energy neutrino. AT2019dsg was identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility. The probability of finding any coincident radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. Assuming that the association is genuine, our observations suggest that tidal disruption events with mildly-relativistic outflows contribute to the cosmic neutrino flux.


title: Phase-Space Spectral Line De-confusion in Intensity Mapping

abstract: Line intensity mapping (LIM) is a promising tool to efficiently probe the three-dimensional large-scale structure by mapping the aggregate emission of a spectral line from all sources that trace the matter density field. Spectral lines from different redshifts can fall in the same observed frequency and be confused, however, which is a major challenge in LIM. In this work, we develop a line de-confusion technique in map space capable of reconstructing the three-dimensional spatial distribution of line-emitting sources. If multiple spectral lines of a source population are observable in multiple frequencies, using the sparse approximation, our technique iteratively extracts sources along a given line of sight by fitting the LIM data to a set of spectral templates. We demonstrate that the technique successfully extracts sources with emission lines present at a few $\sigma$ above the noise level, taking into account uncertainties in the source modeling and presence of continuum foreground contamination and noise fluctuations. As an example, we consider a TIME/CONCERTO-like survey targeting [C II] at the epoch of reionization, and reliably reconstruct the 3D spatial distribution of the CO interlopers at $0.5\lesssim z\lesssim 1.5$. We also demonstrate a successful de-confusion for the SPHEREx mission in the near-infrared wavelengths. Potentially, the reconstructed maps can be further cross-correlated with a (galaxy) tracer population to estimate the total interloper power in the linear clustering regime. This technique is a general framework to extract the phase-space distribution of low-redshift interlopers, without the need of external information, for any line de-confusion problem.


title: Where did the globular clusters of the Milky Way form? Insights from the E-MOSAICS simulations

abstract: Globular clusters (GCs) are typically old, with most having formed at z >~ 2. This makes understanding their birth environments difficult, as they are typically too distant to observe with sufficient angular resolution to resolve GC birth sites. Using 25 cosmological zoom-in simulations of Milky Way-like galaxies from the E-MOSAICS project, with physically-motivated models for star formation, feedback, and the formation, evolution, and disruption of GCs, we identify the birth environments of present-day GCs. We find roughly half of GCs in these galaxies formed in-situ (52.0 +/- 1.0 per cent) between z ~ 2 - 4, in turbulent, high-pressure discs fed by gas that was accreted without ever being strongly heated through a virial shock or feedback. A minority of GCs form during mergers (12.6 +/- 0.6 per cent in major mergers, and 7.2 +/- 0.5 per cent in minor mergers), but we find that mergers are important for preserving the GCs seen today by ejecting them from their natal, high density interstellar medium (ISM), where proto-GCs are rapidly destroyed due to tidal shocks from ISM substructure. This chaotic history of hierarchical galaxy assembly acts to mix the spatial and kinematic distribution of GCs formed through different channels, making it difficult to use observable GC properties to distinguish GCs formed in mergers from ones formed by smooth accretion, and similarly GCs formed in-situ from those formed ex-situ. These results suggest a simple picture of GC formation, in which GCs are a natural outcome of normal star formation in the typical, gas-rich galaxies that are the progenitors of present-day galaxies.


title: Intrinsic sign problems in topological quantum field theories

abstract: The sign problem is a widespread numerical hurdle preventing us from simulating the equilibrium behavior of various problems at the forefront of physics. Focusing on an important sub-class of such problems, bosonic $(2+1)$-dimensional topological quantum field theories, here we provide a simple criterion to diagnose intrinsic sign problems---that is, sign problems that are inherent to that phase of matter and cannot be removed by any local unitary transformation. Explicitly, \textit{if the exchange statistics of the anyonic excitations do not form complete sets of roots of unity, then the model has an intrinsic sign problem}. This establishes a concrete connection between the statistics of anyons, contained in the modular $S$ and $T$ matrices, and the presence of a sign problem in a microscopic Hamiltonian. Furthermore, it places constraints on the phases that can be realised by stoquastic Hamiltonians. We prove this and a more restrictive criterion for the large set of gapped bosonic models described by an abelian topological quantum field theory at low-energy, and offer evidence that it applies more generally with analogous results for non-abelian and chiral theories.

categories: cond-mat.str-el cond-mat.stat-mech hep-lat quant-ph


title: Accretion of Galaxy Groups into Galaxy Clusters

abstract: We study the role of group infall in the assembly and dynamics of galaxy clusters in $\Lambda$CDM. We select $10$ clusters with virial mass $M_{\rm 200} \sim 10^{14} \, M_{\odot}$ from the cosmological hydrodynamical simulation Illustris and follow their galaxies with stellar mass $M_{\star} \geq 1.5 \times 10^8 \, M_{\odot}$. A median of $\sim 38\%$ of surviving galaxies at $z=0$ are accreted as part of groups and did not infall directly from the field, albeit with significant cluster-to-cluster scatter. The evolution of these galaxy associations is quick, with observational signatures of their common origin eroding rapidly in $1$-$3$ Gyr after infall. Substructure plays a dominant role in fostering the conditions for galaxy mergers to happen, even within the cluster environment. Integrated over time, we identify (per cluster) an average of $17 \pm 6$ mergers that occur in infalling galaxy associations, of which $7 \pm 3$ occur well within the virial radius of their cluster hosts. The number of mergers show large dispersion from cluster to cluster, with our most massive system having $42$ mergers above our mass cut-off. These mergers, which are typically gas rich for dwarfs and a combination of gas rich and gas poor for $M_{\star} \sim 10^{11} \, M_{\odot}$, may contribute significantly within $\Lambda$CDM to the formation of specific morphologies, such as lenticulars (S0) and blue compact dwarfs in groups and clusters.


title: Non-ideal magnetohydrodynamics vs turbulence I: Which is the dominant process in protostellar disc formation?

abstract: Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of circumstellar discs that form nearly simultaneously with the formation of the protostar. This is done by modelling the gravitational collapse of a 1~M$_\odot$ gas cloud that is threaded with a magnetic field and imposed with both rotational and turbulent velocities. We investigate magnetic fields that are parallel/anti-parallel and perpendicular to the rotation axis, two rotation rates and four Mach numbers. Disc formation occurs preferentially in the models that include non-ideal MHD where the magnetic field is anti-parallel or perpendicular to the rotation axis. This is independent of the initial rotation rate and level of turbulence, suggesting that subsonic turbulence plays a minimal role in influencing the formation of discs. Aside from first core outflows which are influenced by the initial level of turbulence, non-ideal MHD processes are more important than turbulent processes during the formation of discs around low-mass stars.


title: Three- and two-point spatial correlations of IGM at $z\sim 2$: Cloud based analysis using simulations

abstract: Ly$\alpha$ forest absorption spectra decomposed into multiple Voigt profile components (clouds) allow us to study clustering of intergalactic medium (IGM) as a function of HI column density ($N_{\rm HI}$). Here, we explore the transverse three-point correlation ($\zeta$) of these Ly$\alpha$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We find $\zeta$ to depend strongly on $N_{\rm HI}$ and scale and weakly on angle ($\theta$) of the triplet configuration. We show that the "hierarchical ansatz" is applicable for scales $\ge~ 3h^{-1}$cMpc, and obtain a median reduced three-point correlation (Q) in the range 0.2-0.7. We show, $\zeta$ is influenced strongly by the thermal and ionization state of the gas. As found in the case of galaxies, the influence of physical parameters on Q is weaker compared to that of $\zeta$. We show difference in $\zeta$ and Q between different simulations are minimized if we use appropriate $N_{\rm HI}$ cut-offs corresponding to a given baryon over-density ($\Delta$) using the measured $N_{\rm HI}~vs~\Delta$ relationship obtained from individual simulations. Additionally, we see the effect of pressure broadening on $\zeta$ in a model with artificially boosted heating rates. However, for models with realistic thermal and ionization histories the effect of pressure broadening on $\zeta$ is weak and sub-dominant compared to other local effects. We find strong redshift evolution shown by $\zeta$, mainly originating from the redshift evolution of thermal and ionization state of the IGM. We discuss the observational requirements for the detection of three-point correlation, specifically, in small intervals of configuration parameters and redshift.


title: Rozansky-Witten geometry of Coulomb branches and logarithmic knot invariants

abstract: By studying Rozansky-Witten theory with non-compact target spaces we find new connections with knot invariants whose physical interpretation was not known. This opens up several new avenues, which include a new formulation of $q$-series invariants of 3-manifolds in terms of affine Grassmannians and a generalization of Akutsu-Deguchi-Ohtsuki knot invariants.


title: Non-ideal magnetohydrodynamics vs turbulence II: Which is the dominant process in stellar core formation?

abstract: Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of protostars by following the gravitational collapse of 1~M$_\odot$ gas clouds through the first hydrostatic core to stellar densities. The clouds are imposed with both rotational and turbulent velocities, and are threaded with a magnetic field that is parallel/anti-parallel or perpendicular to the rotation axis; we investigate two rotation rates and four Mach numbers. The initial radius and mass of the stellar core are only weakly dependent on the initial parameters. In the models that include ideal MHD, the magnetic field strength implanted in the protostar at birth is much higher than observed, independent of the initial level of turbulence; only non-ideal MHD can reduce this strength to near or below the observed levels. This suggests that not only is ideal MHD an incomplete picture of star formation, but that the magnetic fields in low mass stars are implanted later in life by a dynamo process. Non-ideal MHD suppresses magnetically launched stellar core outflows, but turbulence permits thermally launched outflows to form a few years after stellar core formation.


title: Origin of cosmological temperature

abstract: A classical solution of the Standard Model + General Relativity is given by an elliptic function whose periodicity in imaginary time is the origin of cosmological temperature. Nothing beyond the Standard Model is assumed. The solution is a $\mathrm{Spin}(4)$-symmetric universe expanding prior to the electroweak transition. A rapidly oscillating $\mathrm{SU}(2)$ gauge field holds the Higgs field to $0$ with strength inversely proportional to the scale factor $a$. When $a$ reaches $a_{\scriptscriptstyle\mathrm{EW}}$ the solution becomes unstable and the electoweak transition begins. $a_{\scriptscriptstyle\mathrm{EW}}$ is the only free parameter in the solution. The temperature at $a_{\scriptscriptstyle\mathrm{EW}}$ is $m_{H}/{(6\pi)^{1/2}} = 28.8\,\text{GeV} = 3.34 \times 10^{14}\, \text{K}$ whatever the value of $a_{\scriptscriptstyle\mathrm{EW}}$.


title: MOTS in Schwarzschild: multiple self-intersections and extreme mass ratio mergers

abstract: We study the open and closed axisymmetric marginally outer trapped surfaces contained in leaves of constant Painlev\'e-Gullstrand time for Schwarzschild spacetimes. We identify a family of closed MOTS in the black hole interior characterized by an arbitrary number of self-intersections. This suggests that the self-intersecting behaviour reported in [Phys.Rev.D 100, 084044 (2019)] may be a far more generic phenomenon than expected. We also consider open surfaces, finding that their behaviour is highly constrained but includes surfaces with multiple self-intersections inside the horizon. We argue that the behaviour of open MOTS identifies and constrains the possible local behaviour of MOTS during extreme mass ratio mergers.


title: Atomic frustrated impurity states in Weyl metals

abstract: We theoretically analyze the effect of the inversion symmetry breaking on the structure of the impurity molecular states in Weyl metals. We show that for the case of a highly noncentrosymmetric Weyl metallic host, the standard picture of the alternating bonding and antibonding orbitals breaks down, and a qualitatively different frustrated atomic state emerges. This is a consequence of the pseudogap closing and related delicate Fano interplay between intra- and inter-impurity scattering channels.


title: Aharonov-Bohm Oscillations in Minimally Twisted Bilayer Graphene

abstract: We investigate transport in the network of valley Hall states that emerges in minimally twisted bilayer graphene under interlayer bias. To this aim, we construct a scattering theory that captures the network physics. In the absence of forward scattering, symmetries constrain the network model to a single parameter that interpolates between one-dimensional chiral zigzag modes and pseudo-Landau levels. Moreover, we show how the coupling of zigzag modes affects magnetotransport. In particular, we find that scattering between parallel zigzag channels gives rise to Aharonov-Bohm oscillations that are robust against temperature, while coupling between zigzag modes propagating in different directions leads to Shubnikov-de Haas oscillations that are smeared out at finite temperature.


title: Dark Matter Subhalos, Strong Lensing and Machine Learning

abstract: We investigate the possibility of applying machine learning techniques to images of strongly lensed galaxies to detect a low mass cut-off in the spectrum of dark matter sub-halos within the lens system. We generate lensed images of systems containing substructure in seven different categories corresponding to lower mass cut-offs ranging from $10^9M_\odot$ down to $10^6M_\odot$. We use convolutional neural networks to perform a multi-classification sorting of these images and see that the algorithm is able to correctly identify the lower mass cut-off within an order of magnitude to better than 93% accuracy.


title: Electromagnetic dual Einstein-Maxwell-scalar models

abstract: Electromagnetic duality is discussed in the context of Einstein-Maxwell-scalar (EMS) models including axionic-type couplings. This family of models introduces two non-minimal coupling functions $f(\phi)$ and $g(\phi)$, depending on a real scalar field $\phi$. Interpreting the scalar field as a medium, one naturally defines constitutive relations as in relativistic non-linear electrodynamics. Requiring these constitutive relations to be invariant under the $SO(2)$ electromagnetic duality rotations of Maxwell's theory, defines 1-parameter, closed $\textit{duality orbits}$ in the space of EMS models, connecting different electromagnetic fields in "dual" models with different coupling functions, but leaving both the scalar field and the spacetime geometry invariant. This mapping works as a solution generating technique, extending any given solution of a specific model to a (different) solution for any of the dual models along the whole duality orbit. We illustrate this technique by considering the duality orbits seeded by specific EMS models wherein solitonic and black hole solutions are known. For dilatonic models, specific rotations are equivalent to $S$-duality.


title: Light bending by a slowly rotating source in quadratic theories of gravity

abstract: In this paper we study the light bending caused by a slowly rotating source in the context of quadratic theories of gravity, in which the Einstein--Hilbert action is extended by additional terms quadratic in the curvature tensors. The deflection angle is computed employing the method based on the Gauss--Bonnet theorem and working in the approximation of a weak lens; also, we assume that the source and observer are at an infinite distance. The formalism presented is very general and applies to any spacetime metric in the limit of weak gravitational field and slow rotation. We find the explicit formula for the deflection angle for several local and nonlocal theories, and also discuss some phenomenological implications.


title: Well-conditioned eigenvalue problems that overflow

abstract: In this note we present a parameterized class of lower triangular matrices. The components of the eigenvectors grow rapidly and will exceed the representational range of any finite number system. The eigenvalues and the eigenvectors are well-conditioned with respect to componentwise relative perturbations of the matrix. This class of matrices is well suited for testing software for computing eigenvectors as these routines must be able to handle overflow successfully.


title: Singularity resolution depends on the clock

abstract: We study the quantum cosmology of a flat Friedmann-Lema\^{i}tre-Robertson-Walker universe filled with a (free) massless scalar field and a perfect fluid that represents radiation or a cosmological constant whose value is not fixed by the action, as in unimodular gravity. We study two versions of the quantum theory: the first is based on a time coordinate conjugate to the radiation/dark energy matter component, i.e., conformal time (for radiation) or unimodular time. As shown by Gryb and Th\'ebault, this quantum theory achieves a type of singularity resolution; we illustrate this and other properties of this theory. The theory is then contrasted with a second type of quantisation in which the logarithm of the scale factor serves as time, which has been studied in the context of the "perfect bounce" for quantum cosmology. Unlike the first quantum theory, the second one contains semiclassical states that follow classical trajectories and evolve into the singularity without obstruction, thus showing no singularity resolution. We discuss how a complex scale factor best describes the semiclassical dynamics. This cosmological model serves as an illustration of the problem of time in quantum cosmology.


title: Tunable coupling of two mechanical resonators by a graphene membrane

abstract: Coupled nanomechanical resonators are interesting for both fundamental studies and practical applications as they offer rich and tunable oscillation dynamics. At present, the mechanical coupling in such systems is often mediated by a fixed geometry, such as a joint clamping point of the resonators or a displacement-dependent force. Here we show a graphene-integrated electromechanical system consisting of two physically separated mechanical resonators -- a comb-drive actuator and a suspended silicon beam -- that are tunably coupled by a graphene membrane. The graphene membrane, moreover, provides a sensitive electrical read-out for the two resonating systems silicon structures showing 16 different modes in the frequency range from 0.4~to 24~MHz. In addition, by pulling on the graphene membrane with an electrostatic potential applied to one of the silicon resonators, we control the mechanical coupling, quantified by the $g$-factor, from 20 kHz to 100 kHz. Our results pave the way for coupled nanoelectromechanical systems requiring controllable mechanically coupled resonators.


title: A Pattern-based Approach to Detect and Improve Non-descriptive Test Names

abstract: Unit tests are an important artifact that supports the software development process in several ways. For example, when a test fails, its name can provide the first step towards understanding the purpose of the test. Unfortunately, unit tests often lack descriptive names. In this paper, we propose a new, pattern-based approach that can help developers improve the quality of test names of JUnit tests by making them more descriptive. It does this by detecting non-descriptive test names and in some cases, providing additional information about how the name can be improved. Our approach was assessed using an empirical evaluation on 34352 JUnit tests. The results of the evaluation show that the approach is feasible, accurate, and useful at discriminating descriptive and non-descriptive names with a 95% true-positive rate.


title: What is interdisciplinarity? Views from students and professors in a non-major intro physics course

abstract: We present an investigation into the interdisciplinary role of physics in a physics-for-non-physicists course at Pomona College. This work is guided by prior research into introductory physics for life-science (IPLS) courses, but attends to significant differences in the scope and context of this course. We interviewed enrolled students, physics professors, and professors from non-physics disciplines to explore the function of this course and the role of physics in the education of non-physics-science students. Interviews were audio recorded and transcribed, then analyzed to identify emergent themes. These themes outline the authentic physics, including content knowledge and other, broader learning objectives, that play an important and distinct role in the science education of enrolled students. Stakeholders generally align in their emphasis of interdisciplinary relevance with some divergence in the specific articulation of that idea. The differences can be understood through the stakeholders' distinct areas of expertise, with non-physics professors expressing value through relevance to their discipline and physics professors focusing on essential aspects of physics.


title: Multiplexed Single Photons from Deterministically Positioned Nanowire Quantum Dots

abstract: Solid-state quantum emitters are excellent sources of on-demand indistinguishable or entangled photons and can host long-lived spin memories, crucial resources for photonic quantum information applications. However, their scalability remains an outstanding challenge. Here we present a scalable technique to multiplex streams of photons from multiple independent quantum dots, on-chip, into a fiber network for use off-chip. Multiplexing is achieved by incorporating a multi-core fiber into a confocal microscope and spatially matching the multiple foci, seven in this case, to quantum dots in an array of deterministically positioned nanowires. First, we report the coherent control of the emission of biexciton-exciton cascade from a single nanowire quantum dot under resonant two-photon excitation. Then, as a proof-of-principle demonstration, we perform parallel spectroscopy on the nanowire array to identify two nearly identical quantum dots at different positions which are subsequently tuned into resonance with an external magnetic field. Multiplexing of background-free single photons from these two quantum dots is then achieved. Our approach, applicable to all types of quantum emitters, can readily be scaled up to multiplex $>100$ quantum light sources, providing a breakthrough in hardware for photonic based quantum technologies. Immediate applications include quantum communication, quantum simulation, and quantum computation.


title: Minimal Model for Fast Scrambling

abstract: We study quantum information scrambling in spin models with both long-range all-to-all and short-range interactions. We argue that a simple global, spatially homogeneous interaction together with local chaotic dynamics is sufficient to give rise to fast scrambling, which describes the spread of quantum information over the entire system in a time that is logarithmic in the system size. This is illustrated in two tractable models: (1) a random circuit with Haar random local unitaries and a global interaction and (2) a classical model of globally coupled non-linear oscillators. We use exact numerics to provide further evidence by studying the time evolution of an out-of-time-order correlator and entanglement entropy in spin chains of intermediate sizes. Our results pave the way towards experimental investigations of fast scrambling and aspects of quantum gravity with quantum simulators.


title: Machine Learning for Physics-Informed Generation of Dispersed Multiphase Flow Using Generative Adversarial Networks

abstract: Fluid flow around a random distribution of stationary spherical particles is a problem of substantial importance in the study of dispersed multiphase flows. In this paper we present a machine learning methodology using Generative Adversarial Network framework and Convolutional Neural Network architecture to recreate particle-resolved fluid flow around a random distribution of monodispersed particles. The model was applied to various Reynolds number and particle volume fraction combinations spanning over a range of [2.69, 172.96] and [0.11, 0.45] respectively. Test performance of the model for the studied cases is very promising.


title: A Repo Model of Fire Sales with VWAP and LOB Pricing Mechanisms

abstract: We consider a network of banks that optimally choose a strategy of asset liquidations and borrowing in order to cover short term obligations. The borrowing is done in the form of collateralized repurchase agreements, the haircut level of which depends on the total liquidations of all the banks. Similarly the fire-sale price of the asset obtained by each of the banks depends on the amount of assets liquidated by the bank itself and by other banks. By nature of this setup, banks' behavior is considered as a Nash equilibrium. This paper provides two forms for market clearing to occur: through a common closing price and through an application of the limit order book. The main results of this work are providing the existence of maximal and minimal clearing solutions (i.e., liquidations, borrowing, fire sale prices, and haircut levels) as well as sufficient conditions for uniqueness of the clearing solutions.


title: Solar Bayesian Analysis Toolkit -- a new Markov chain Monte Carlo IDL code for Bayesian parameter inference

abstract: We present the Solar Bayesian Analysis Toolkit (SoBAT) which is a new easy to use tool for Bayesian analysis of observational data, including parameter inference and model comparison. SoBAT is aimed (but not limited) to be used for the analysis of solar observational data. We describe a new Interactive Data Language (IDL) code designed to facilitate the comparison of user-supplied model with data. Bayesian inference allows prior information to be taken into account. The use of Markov chain Monte Carlo (MCMC) sampling allows efficient exploration of large parameter spaces and provides reliable estimation of model parameters and their uncertainties. The Bayesian evidence for different models can be used for quantitative comparison. The code is tested to demonstrate its ability to accurately recover a variety of parameter probability distributions. Its application to practical problems is demonstrated using studies of the structure and oscillation of coronal loops.


title: Flavour symmetries in the SMEFT

abstract: We analyse how $U(3)^5$ and $U(2)^5$ flavour symmetries act on the Standard Model Effective Field Theory, providing an organising principle to classify the large number of dimension-six operators involving fermion fields. A detailed counting of such operators, at different order in the breaking terms of both these symmetries, is presented. A brief discussion about possible deviations from these two reference cases, and a simple example of the usefulness of this classification scheme for high-$p_T$ analyses at the LHC, are also presented.


title: Stability of clusters in the second-order Kuramoto model on random graphs

abstract: The Kuramoto model of coupled phase oscillators with inertia on Erdos-Renyi graphs is analyzed in this work. For a system with intrinsic frequencies sampled from a bimodal distribution we identify a variety of two cluster patterns and study their stability. To this end, we decompose the description of the cluster dynamics into two systems: one governing the (macro) dynamics of the centers of mass of the two clusters and the second governing the (micro) dynamics of individual oscillators inside each cluster. The former is a low-dimensional ODE whereas the latter is a system of two coupled Vlasov PDEs. Stability of the cluster dynamics depends on the stability of the low-dimensional group motion and on coherence of the oscillators in each group. We show that the loss of coherence in one of the clusters leads to the loss of stability of a two-cluster state and to formation of chimera states. The analysis of this paper can be generalized to cover states with more than two clusters and to coupled systems on W-random graphs. Our results apply to a model of a power grid with fluctuating sources.


title: Testing the violation of the equivalence principle in the electromagnetic sector and its consequences in $f(T)$ gravity

abstract: A violation of the distance-duality relation is directly linked with a temporal variation of the electromagnetic fine-structure constant. We consider a number of well-studied $f(T)$ gravity models and we revise the theoretical prediction of their corresponding induced violation of the distance-duality relationship. We further extract constraints on the involved model parameters through fine-structure constant variation data, alongside with supernovae data, and Hubble parameter measurements. Moreover, we constrain the evolution of the effective $f(T)$ gravitational constant. Finally, we compare with revised constraints on the phenomenological parametrisations of the violation of the equivalence principle in the electromagnetic sector.


title: Micro Cold Traps on the Moon

abstract: Water ice is thought to be trapped in large permanently shadowed regions (PSRs) in the Moon's polar regions, due to their extremely low temperatures. Here, we show that many unmapped cold traps exist on small spatial scales, substantially augmenting the areas where ice may accumulate. Using theoretical models and data from the Lunar Reconnaissance Orbiter, we estimate the contribution of shadows on scales from 1 km down to 1 cm, the smallest distance over which we find cold-trapping to be effective for water ice. Approximately 10-20\% of the permanent cold trap area for water is found to be contained in these "micro cold traps," which are the most numerous cold traps on the Moon. Consideration of all spatial scales therefore substantially increases the number of cold traps over previous estimates, for a total area of ~40,000 km^2. A majority of cold traps for water ice is found at latitudes >80{\deg} because permanent shadows equatorward of 80{\deg} are typically too warm to support ice accumulation. Our results show that water trapped at the lunar poles may be more accessible as a resource for future missions than previously thought.


title: Electro-Quasistatic Animal Body Communication for Chronic Untethered Rodent Biopotential Recording

abstract: Continuous multi-channel monitoring of biopotential signals is vital in understanding the body as a whole, facilitating accurate models and predictions in neural research. The current state of the art in wireless technologies for untethered biopotential recordings rely on radiative electromagnetic (EM) fields. In such transmissions, only a small fraction of this energy is received since the EM fields are widely radiated resulting in lossy inefficient systems. Using the body as a communication medium (similar to a 'wire') allows for the containment of the energy within the body, yielding order(s) of magnitude lower loss than radiative EM communication. In this work, we introduce Animal Body Communication (ABC), which utilizes the concept of using the body as a medium into the domain of chronic animal biopotential recording. This work, for the first time, develops the theory and models for animal body communication circuitry and channel loss. Using this theoretical model, a sub-inch$^3$, custom-designed sensor node is built using off the shelf components which is capable of sensing and transmitting biopotential signals, through the body of the rat at significantly lower powers compared to traditional wireless transmissions. In-vivo experimental analysis proves that ABC successfully transmits acquired electrocardiogram (EKG) signals through the body with correlation accuracy >99% when compared to traditional wireless communication modalities, with a 50x reduction in power consumption.


title: A Parallel Hybrid Technique for Multi-Noise Removal from Grayscale Medical Images

abstract: Medical imaging is the technique used to create images of the human body or parts of it for clinical purposes. Medical images always have large sizes and they are commonly corrupted by single or multiple noise type at the same time, due to various reasons, these two reasons are the triggers for moving toward parallel image processing to find alternatives of image de-noising techniques. This paper presents a parallel hybrid filter implementation for gray scale medical image de-noising. The hybridization is between adaptive median and wiener filters. Parallelization is implemented on the adaptive median filter to overcome the latency of neighborhood operation, parfor implicit parallelism powered by MatLab 2013a is used. The implementation is tested on an image of 2.5 MB size, which is divided into 2, 4 and 8 partitions; a comparison between the proposed implementation and sequential implementation is given, in terms of time. Thus, each case has the best time when assigned to number of threads equal to the number of its partitions. Moreover, Speed up and efficiency are calculated for the algorithm and they show a measured enhancement.


title: The string C-group representations of the Suzuki and Rudvalis sporadic groups

abstract: We present new algorithms to classify all string C-group representations of a given group $G$. We use these algorithms to classify all string C-group representations of the sporadic groups of Suzuki and Rudvalis.


title: Correlated electron-hole State in Twisted Double Bilayer Graphene

abstract: When twisted to angles near 1{\deg}, graphene multilayers provide a new window on electron correlation physics by hosting gate-tuneable strongly-correlated states, including insulators, superconductors, and unusual magnets. Here we report the discovery of a new member of the family, density-wave states, in double bilayer graphene twisted to 2.37{\deg}. At this angle the moir\'e states retain much of their isolated bilayer character, allowing their bilayer projections to be separately controlled by gates. We use this property to generate an energetic overlap between narrow isolated electron and hole bands with good nesting properties. Our measurements reveal the formation of ordered states with reconstructed Fermi surfaces, consistent with density-wave states, for equal electron and hole densities. These states can be tuned without introducing chemical dopants, thus opening the door to a new class of fundamental studies of density-waves and their interplay with superconductivity and other types of order, a central issue in quantum matter physics.


title: Approaching the upper limits of the local density of states via optimized metallic cavities

abstract: By computational optimization of air-void cavities in metallic substrates, we show that the local density of states (LDOS) can reach within a factor of $\approx 10$ of recent theoretical upper limits, and within a factor $\approx 4$ for the single-polarization LDOS, demonstrating that the theoretical limits are nearly attainable. Optimizing the total LDOS results in a spontaneous symmetry breaking where it is preferable to couple to a specific polarization. Moreover, simple shapes such as optimized cylinders attain nearly the performance of complicated many-parameter optima, suggesting that only one or two key parameters matter in order to approach the theoretical LDOS bounds for metallic resonators.


title: Computational Investigation of Copper Phosphides as Conversion Anodes for Lithium-Ion Batteries

abstract: Using first principles structure searching with density-functional theory (DFT) we identify a novel $Fm\bar{3}m$ phase of Cu$_2$P and two low-lying metastable structures, an $I\bar{4}3d$--Cu$_3$P phase, and a $Cm$--Cu$_3$P$_{11}$ phase. The computed pair distribution function of the novel $Cm$--Cu$_3$P$_{11}$ phase shows its structural similarity to the experimentally identified $Cm$--Cu$_2$P$_7$ phase. The relative stability of all Cu--P phases at finite temperatures is determined by calculating the Gibbs free energy using vibrational effects from phonon modes at 0 K. From this, a finite-temperature convex hull is created, on which $Fm\bar{3}m$--Cu$_2$P is dynamically stable and the Cu$_{3-x}$P ($x < 1$) defect phase $Cmc2_1$--Cu$_8$P$_3$ remains metastable (within 20 meV/atom of the convex hull) across a temperature range from 0 K to 600 K. Both CuP$_2$ and Cu$_3$P exhibit theoretical gravimetric capacities higher than contemporary graphite anodes for Li-ion batteries; the predicted Cu$_2$P phase has a theoretical gravimetric capacity of 508 mAh/g as a Li-ion battery electrode, greater than both Cu$_3$P (363 mAh/g) and graphite (372 mAh/g). Cu$_2$P is also predicted to be both non-magnetic and metallic, which should promote efficient electron transfer in the anode. Cu$_2$P's favorable properties as a metallic, high-capacity material suggest its use as a future conversion anode for Li-ion batteries; with a volume expansion of 99% during complete cycling, Cu$_2$P anodes could be more durable than other conversion anodes in the Cu--P system with volume expansions greater than 150%.


title: Provably insecure group authentication: Not all security proofs are what they claim to be

abstract: A paper presented at the ICICS 2019 conference describes what is claimed to be a `provably secure group authentication [protocol] in the asynchronous communication model'. We show here that this is far from being the case, as the protocol is subject to serious attacks. To try to explain this troubling case, an earlier (2013) scheme on which the ICICS 2019 protocol is based was also examined and found to possess even more severe flaws - this latter scheme was previously known to be subject to attack, but not in quite as fundamental a way as is shown here. Examination of the security theorems provided in both the 2013 and 2019 papers reveals that in neither case are they exactly what they seem to be at first sight; the issues raised by this are also briefly discussed.


title: Dimensionless equations in non-relativistic quantum mechanics

abstract: We discuss the numerous advantages of using dimensionless equations in non-relativistic quantum mechanics. Dimensionless equations are considerably simpler and reveal the number of relevant parameters in the models. They are less prone to round-off errors when applying numerical methods because all the quantities are of the other of unity. A dimensionless equation facilitates the application of perturbation theory and provides a glimpse of the sort of solution we are going to obtain beforehand.


title: Classification of finite irreducible conformal modules over $N=2$ Lie conformal superalgebras of Block type

abstract: We introduce the $N=2$ Lie conformal superalgebras ${\frak {K}}(p)$ of Block type, and classify their finite irreducible conformal modules for any nonzero parameter $p$. where $p$ is a nonzero complex number. In particular, we show that such a conformal module admits a nontrivial extension of a finite conformal module $M$ over $K_2$ if $p=-1$ and $M$ has rank $(2+2)$, where $K_2$ is an $N=2$ conformal subalgebra of ${\frak {K}}(p)$. As a byproduct, we obtain the classification of finite irreducible conformal modules over a series of finite Lie conformal superalgebras ${\frak k}(n)$ for $n\ge1$. Composition factors of all the involved reducible conformal modules are also determined.


title: Gap Sets for the Spectra of Cubic Graphs

abstract: We study gaps in the spectra of the adjacency matrices of large finite cubic graphs. It is known that the gap intervals $(2 \sqrt{2},3)$ and $[-3,-2)$ achieved in cubic Ramanujan graphs and line graphs are maximal. We give constraints on spectra in [-3,3] which are maximally gapped and construct examples which achieve these bounds. These graphs yield new instances of maximally gapped intervals. We also show that every point in $[-3,3)$ can be gapped by cubic graphs, even by planar ones. Our results show that the study of spectra of cubic, and even planar cubic, graphs is subtle and very rich.


title: Roles of Dynamic State Estimation in Power System Modeling, Monitoring and Operation

abstract: Power system dynamic state estimation (DSE) remains an active research area. This is driven by the absence of accurate models, the increasing availability of fast-sampled, time-synchronized measurements, and the advances in the capability, scalability, and affordability of computing and communications. This paper discusses the advantages of DSE as compared to static state estimation, and the implementation differences between the two, including the measurement configuration, modeling framework and support software features. The important roles of DSE are discussed from modeling, monitoring and operation aspects for today's synchronous machine dominated systems and the future power electronics-interfaced generation systems. Several examples are presented to demonstrate the benefits of DSE on enhancing the operational robustness and resilience of 21st century power system through time critical applications. Future research directions are identified and discussed, paving the way for developing the next generation of energy management systems.


title: Clasper Concordance, Whitney towers and repeating Milnor invariants

abstract: We show that for each $k\in \mathbb{N}$, a link $L\subset S^3$ bounds a degree $k$ Whitney tower in the 4-ball if and only if it is \emph{$C_k$-concordant} to the unlink. This means that $L$ is obtained from the unlink by a finite sequence of concordances and degree $k$ clasper surgeries. In our construction the trees associated to the Whitney towers coincide with the trees associated to the claspers. As a corollary to our previous obstruction theory for Whitney towers in the 4-ball, it follows that the $C_k$-concordance filtration of links is classified in terms of Milnor invariants, higher-order Sato-Levine and Arf invariants. Using a new notion of $k$-repeating twisted Whitney towers, we also classify a natural generalization of the notion of link homotopy, called twisted \emph{self $C_k$-concordance}, in terms of $k$-repeating Milnor invariants and $k$-repeating Arf invariants.


title: A search for the $Z\gamma$ decay mode of the Higgs boson in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

abstract: A search for the $Z\gamma$ decay of the Higgs boson, with $Z$ boson decays into pairs of electrons or muons is presented. The analysis uses proton$-$proton collision data at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS detector at the Large Hadron Collider. The observed data are consistent with the expected background with a $p$-value of 1.3%. An upper limit at 95% confidence level on the production cross-section times the branching ratio for $pp\to H\to Z\gamma$ is set at 3.6 times the Standard Model prediction while 2.6 times is expected in the presence of the Standard Model Higgs boson. The best-fit value for the signal yield normalised to the Standard Model prediction is $2.0^{+1.0}_{-0.9}$ where the statistical component of the uncertainty is dominant.


title: Breaking the Selection Rules of Spin-Forbidden Molecular Absorption in Plasmonic Nanocavities

abstract: Controlling absorption and emission of organic molecules is crucial for efficient light-emitting diodes, organic solar cells and single-molecule spectroscopy. Here, a new molecular absorption is activated inside a gold plasmonic nanocavity, and found to break selection rules via spin-orbit coupling. Photoluminescence excitation scans reveal absorption from a normally spin-forbidden singlet to triplet state transition, while drastically enhancing the emission rate by several thousand fold. The experimental results are supported by density functional theory, revealing the manipulation of molecular absorption by nearby metallic gold atoms.


title: Algebraically cofibrant and fibrant objects revisited

abstract: We extend all known results about transferred model structures on algebraically cofibrant and fibrant objects by working with weak model categories. We show that for an accessible weak model category there are always Quillen equivalent transferred weak model structures on both the categories of algebraically cofibrant and algebraically fibrant objects. Under additional assumptions, these transferred weak model structures are shown to be left, right or Quillen model structures. By combining both constructions, we show that each combinatorial weak model category is connected, via a chain of Quillen equivalences, to a combinatorial Quillen model category in which all objects are fibrant.


title: Process Knowledge Driven Change Point Detection for Automated Calibration of Discrete Event Simulation Models Using Machine Learning

abstract: Initial development and subsequent calibration of discrete event simulation models for complex systems require accurate identification of dynamically changing process characteristics. Existing data driven change point methods (DD-CPD) assume changes are extraneous to the system, thus cannot utilize available process knowledge. This work proposes a unified framework for process-driven multi-variate change point detection (PD-CPD) by combining change point detection models with machine learning and process-driven simulation modeling. The PD-CPD, after initializing with DD-CPD's change point(s), uses simulation models to generate system level outputs as time-series data streams which are then used to train neural network models to predict system characteristics and change points. The accuracy of the predictive models measures the likelihood that the actual process data conforms to the simulated change points in system characteristics. PD-CPD iteratively optimizes change points by repeating simulation and predictive model building steps until the set of change point(s) with the maximum likelihood is identified. Using an emergency department case study, we show that PD-CPD significantly improves change point detection accuracy over DD-CPD estimates and is able to detect actual change points.


title: Design and visualization of Riemannian metrics

abstract: Local and global illumination were recently defined in Riemannian manifolds to visualize classical Non-Euclidean spaces. This work focuses on Riemannian metric construction in $\mathbb{R}^3$ to explore special effects like warping, mirages, and deformations. We investigate the possibility of using graphs of functions and diffeomorphism to produce such effects. For these, their Riemannian metrics and geodesics derivations are provided, and ways of accumulating such metrics. We visualize, in "real-time", the resulting Riemannian manifolds using a ray tracing implemented on top of Nvidia RTX GPUs.


title: Computationally Inequivalent Summations and Their Parenthetic Forms

abstract: Floating-point addition on a finite-precision machine is not associative, so not all mathematically equivalent summations are computationally equivalent. Making this assumption can lead to numerical error in computations. Proper ordering and parenthesizing is a low-overhead way of mitigating such error in a floating point summation. Ordered and parenthesized summations fall into equivalence classes. We describe these classes, and the parenthetic forms summations in these classes take. We provide summation-related interpretations for sequences known in other contexts, and give new recursive and closed formulas for sequences not previously related to summation. We also introduce a data structure that facilitates understanding of these objects, and use it to consider certain forms of summation used by default in widely used computer languages. Finally, we relate this data structure to other mathematical constructs from the fields of mathematical analysis and algorithmic analysis.


title: Effective theories for a nonrelativistic field in an expanding universe: Induced self-interaction, pressure, sound speed, and viscosity

abstract: A massive, nonrelativistic scalar field in an expanding spacetime is usually approximated by a pressureless perfect fluid, which leads to the standard conclusion that such a field can play the role of cold dark matter. In this paper, we systematically study these approximations, incorporating subleading corrections. We provide two equivalent effective descriptions of the system, each of which offers its own advantages and insights: (i) A nonrelativistic effective field theory (EFT) with which we show that the relativistic corrections induce an effective self-interaction for the nonrelativistic field. As a byproduct, our EFT also allows one to construct the exact solution, including oscillatory behavior, which is often difficult to achieve from the exact equations. (ii) An effective (imperfect) fluid description, with which we demonstrate that, for a perturbed Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) universe: (a) The pressure is small but nonzero (and positive), even for a free theory with no tree-level self-interactions. (b) The sound speed of small fluctuations is also nonzero (and positive), reproducing already known leading-order results, correcting a subdominant term, and identifying a new contribution that had been omitted in previous analyses. (c) The fluctuations experience a negative effective bulk viscosity. The positive sound speed and the negative bulk viscosity act in favor of and against the growth of overdensities, respectively. The net effect may be considered a smoking gun for ultra-light dark matter.


title: Citations versus expert opinions: Citation analysis of Featured Reviews of the American Mathematical Society

abstract: Peer review and citation metrics are two means of gauging the value of scientific research, but the lack of publicly available peer review data makes the comparison of these methods difficult. Mathematics can serve as a useful laboratory for considering these questions because as an exact science, there is a narrow range of reasons for citations. In mathematics, virtually all published articles are post-publication reviewed by mathematicians in Mathematical Reviews (MathSciNet) and so the data set was essentially the Web of Science mathematics publications from 1993 to 2004. For a decade, especially important articles were singled out in Mathematical Reviews for featured reviews. In this study, we analyze the bibliometrics of elite articles selected by peer review and by citation count. We conclude that the two notions of significance described by being a featured review article and being highly cited are distinct. This indicates that peer review and citation counts give largely independent determinations of highly distinguished articles. We also consider whether hiring patterns of subfields and mathematicians' interest in subfields reflect subfields of featured review or highly cited articles. We reexamine data from two earlier studies in light of our methods for implications on the peer review/citation count relationship to a diversity of disciplines.


title: New bounds on the periodicity transient of the powers of a tropical matrix: using cyclicity and factor rank

abstract: Building on the weak CSR approach developed in a previous paper by Merlet, Nowak and Sergeev, we establish new bounds for the periodicity threshold of the powers of a tropical matrix. According to that approach, bounds on the ultimate periodicity threshold take the form of T=max(T_1,T_2), where T_1 is a bound on the time after which the weak CSR expansion starts to hold and T_2 is a bound on the time after which the first CSR term starts to dominate. The new bounds on T_1 and T_2 established in this paper make use of the cyclicity of the associated graph and the (tropical) factor rank of the matrix, which leads to much improved bounds in favorable cases. For T_1, in particular, we obtain new extensions of bounds of Schwarz, Kim and Gregory-Kirkland-Pullman, previously known as bounds on exponents of digraphs. For similar bounds on T_2, we introduce the novel concept of walk reduction threshold and establish bounds on it that use both cyclicity and factor rank.