title: Elucidating Inferential Models with the Cauchy Distribution

abstract: Statistical inference as a formal scientific method to covert experience to knowledge has proven to be elusively difficult. While frequentist and Bayesian methodologies have been accepted in the contemporary era as two dominant schools of thought, it has been a good part of the last hundred years to see growing interests in development of more sound methods, both philosophically, in terms of scientific meaning of inference, and mathematically, in terms of exactness and efficiency. These include Fisher's fiducial argument, the Dempster-Shafe theory of belief functions, generalized fiducial, Confidence Distributions, and the most recently proposed inferential framework, called Inferential Models. Since it is notoriously challenging to make exact and efficient inference about the Cauchy distribution, this article takes it as an example to elucidate different schools of thought on statistical inference. It is shown that the standard approach of Inferential Models produces exact and efficient prior-free probabilistic inference on the location and scale parameters of the Cauchy distribution, whereas all other existing methods suffer from various difficulties.


title: Reaction-diffusive dynamics of number-conserving dissipative quantum state preparation

abstract: The use of dissipation for the controlled creation of nontrivial quantum many-body correlated states is of much fundamental and practical interest. What is the result of imposing number conservation, which, in closed system, gives rise to diffusive spreading? We investigate this question for a paradigmatic model of a two-band system, with dissipative dynamics aiming to empty one band and to populate the other, which had been introduced before for the dissipative stabilization of topological states. Going beyond the mean-field treatment of the dissipative dynamics, we demonstrate the emergence of a diffusive regime for the particle and hole density modes at intermediate length- and time-scales, which, interestingly, can only be excited in nonlinear response to external fields. We also identify processes that limit the diffusive behavior of this mode at the longest length- and time-scales. Strikingly, we find that these processes lead to a reaction-diffusion dynamics governed by the Fisher-Kolmogorov-Petrovsky-Piskunov equation, making the designed dark state unstable towards a state with a finite particle and hole density.

categories: cond-mat.stat-mech cond-mat.mes-hall cond-mat.quant-gas quant-ph


title: Detection of HCN and diverse redox chemistry in the plume of Enceladus

abstract: The Cassini spacecraft observed that Saturn's moon Enceladus possesses a series of jets erupting from its South Polar Terrain. Previous studies of in situ data collected by Cassini's Ion and Neutral Mass Spectrometer (INMS) have identified H$_2$O, CO$_2$, CH$_4$, NH$_3$, and H$_2$ within the plume of ejected material. Identification of minor species in the plume remains an ongoing challenge, owing to the large number of possible combinations that can be used to fit the INMS data. Here, we present the detection of several new compounds of strong importance to the habitability of Enceladus, including HCN, C$_2$H$_2$, C$_3$H$_6$, and C$_2$H$_6$. Our analyses of the low velocity INMS data, coupled with our detailed statistical framework, enable discrimination between previously ambiguous species in the plume by alleviating the effects of high dimensional model fitting. Together with plausible mineralogical catalysts and redox gradients derived from surface radiolysis, these compounds could potentially support extant microbial communities or drive complex organic synthesis leading to the origin of life.


title: On the switching mechanism and optimisation of ion irradiation enabled 2D $MoS_2$ memristors

abstract: Memristors are prominent passive circuit elements with promising futures for energy-efficient in-memory processing and revolutionary neuromorphic computation. State-of-the-art memristors based on two-dimensional (2D) materials exhibit enhanced tunability, scalability and electrical reliability. However, the fundamental of the switching is yet to be clarified before they can meet industrial standards in terms of endurance, variability, resistance ratio, and scalability. This new physical simulator based on the kinetic Monte Carlo (kMC) algorithm reproduces the defect migration process in 2D materials and sheds light on the operation of 2D memristors. The present work employs the simulator to study a two-dimensional $2H-MoS_2$ planar resistive switching (RS) device with an asymmetric defect concentration introduced by ion irradiation. The simulations unveil the non-filamentary RS process and propose practical routes to optimize the device's performance. For instance, the resistance ratio can be increased by 53% by controlling the concentration and distribution of defects, while the variability can be reduced by 55% by increasing 5-fold the device size from 10 to 50 nm. Our simulator also explains the trade-offs between the resistance ratio and variability, resistance ratio and scalability, and variability and scalability. Overall, the simulator may enable an understanding and optimization of devices to expedite cutting-edge applications.


title: Exact emergent higher-form symmetries in bosonic lattice models

abstract: Although condensed matter systems usually do not have higher-form symmetries, we show that, unlike 0-form symmetry, higher-form symmetries can emerge as exact symmetries at low energies and long distances. In particular, emergent higher-form symmetries at zero temperature are robust to arbitrary local UV perturbations in the thermodynamic limit. This result is true for both invertible and non-invertible higher-form symmetries. Therefore, emergent higher-form symmetries are $\textit{exact emergent symmetries}$: they are not UV symmetries but constrain low-energy dynamics as if they were. Since phases of matter are defined in the thermodynamic limit, this implies that a UV theory without higher-form symmetries can have phases characterized by exact emergent higher-form symmetries. We demonstrate this in three lattice models, the quantum clock model and emergent ${\mathbb{Z}_N}$ and ${U(1)}$ ${p}$-gauge theory, finding regions of parameter space with exact emergent (anomalous) higher-form symmetries. Furthermore, we perform a generalized Landau analysis of a 2+1D lattice model that gives rise to $\mathbb{Z}_2$ gauge theory. Using exact emergent 1-form symmetries accompanied by their own energy/length scales, we show that the transition between the deconfined and Higgs/confined phases is continuous and equivalent to the spontaneous symmetry-breaking transition of a $\mathbb{Z}_2$ symmetry, even though the lattice model has no symmetry. Also, we show that this transition line must $\textit{always}$ contain two parts separated by multi-critical points or other phase transitions. We discuss the physical consequences of exact emergent higher-form symmetries and contrast them to emergent ${0}$-form symmetries. Lastly, we show that emergent 1-form symmetries are no longer exact at finite temperatures, but emergent $p$-form symmetries with ${p\geq 2}$ are.


title: Neural Shadow Mapping

abstract: We present a neural extension of basic shadow mapping for fast, high quality hard and soft shadows. We compare favorably to fast pre-filtering shadow mapping, all while producing visual results on par with ray traced hard and soft shadows. We show that combining memory bandwidth-aware architecture specialization and careful temporal-window training leads to a fast, compact and easy-to-train neural shadowing method. Our technique is memory bandwidth conscious, eliminates the need for post-process temporal anti-aliasing or denoising, and supports scenes with dynamic view, emitters and geometry while remaining robust to unseen objects.


title: Channel Estimation for Reconfigurable Intelligent Surface-Assisted Full-Duplex MIMO with Hardware Impairments

abstract: We consider the problem of channel estimation in a multiple-input-multiple-output (MIMO) full-duplex (FD) wireless communication system assisted by a reconfigurable intelligent surface (RIS) with hardware impairments (HI) occurring at the transceivers and RIS elements. We propose an unbiased channel estimator that requires knowledge of only the first and second order statistics of the HI, for which we derive closed form expressions. The proposed estimator reduces to the maximum likelihood estimator (MLE) in the case of ideal hardware. We also describe FD and HD orthogonal pilot schemes that minimize the mean square error of the MLE in the case of ideal hardware. We verify the performance of the estimator under varying conditions of transceiver and RIS HI via numerical simulations.


title: Security-Aware Approximate Spiking Neural Networks

abstract: Deep Neural Networks (DNNs) and Spiking Neural Networks (SNNs) are both known for their susceptibility to adversarial attacks. Therefore, researchers in the recent past have extensively studied the robustness and defense of DNNs and SNNs under adversarial attacks. Compared to accurate SNNs (AccSNN), approximate SNNs (AxSNNs) are known to be up to 4X more energy-efficient for ultra-low power applications. Unfortunately, the robustness of AxSNNs under adversarial attacks is yet unexplored. In this paper, we first extensively analyze the robustness of AxSNNs with different structural parameters and approximation levels under two gradient-based and two neuromorphic attacks. Then, we propose two novel defense methods, i.e., precision scaling and approximate quantization-aware filtering (AQF), for securing AxSNNs. We evaluated the effectiveness of these two defense methods using both static and neuromorphic datasets. Our results demonstrate that AxSNNs are more prone to adversarial attacks than AccSNNs, but precision scaling and AQF significantly improve the robustness of AxSNNs. For instance, a PGD attack on AxSNN results in a 72\% accuracy loss compared to AccSNN without any attack, whereas the same attack on the precision-scaled AxSNN leads to only a 17\% accuracy loss in the static MNIST dataset (4X robustness improvement). Similarly, a Sparse Attack on AxSNN leads to a 77\% accuracy loss when compared to AccSNN without any attack, whereas the same attack on an AxSNN with AQF leads to only a 2\% accuracy loss in the neuromorphic DVS128 Gesture dataset (38X robustness improvement).


title: A Mixed Stirring Mechanism for Debris Discs with Giant and Dwarf Planetary Perturbations

abstract: Debris discs consist of belts of bodies ranging in size from dust grains to planetesimals; these belts are visible markers of planetary systems around other stars that can reveal the influence of extrasolar planets through their shape and structure. Two key stirring mechanisms -- self-stirring by planetesimals and secular perturbation by an external giant planet -- have been identified to explain the dynamics of planetesimal belts; their relative importance has been studied independently, but are yet to be considered in combination. In this work we perform a suite of 286 N-body simulations exploring the evolution of debris discs over 1~Gyr, combining the gravitational perturbations of both dwarf planets embedded in the discs, and an interior giant planet. Our systems were somewhat modeled after the architecture of the outer Solar system: a Solar mass star, a single massive giant planet at 30~au ($M_{\rm GP} =$ 10 to 316~$\mathrm{M}_{\oplus}$), and a debris disc formed by 100 massive dwarf planets and 1000 massless particles ($M_{\rm DD} =$ 3.16 to 31.6~$\mathrm{M}_{\oplus}$). We present the evolution of both the disc and the giant planet after 1~Gyr. The time evolution of the average eccentricity and inclination of the disc is strongly dependent on the giant planet mass as well as on the remaining disc mass. We also found that efficient stirring is achieved even with small disc masses. In general, we find that a mixed mechanism is more efficient in the stirring of cold debris discs than either mechanism acting in isolation.


title: Improving Reliability of Spiking Neural Networks through Fault Aware Threshold Voltage Optimization

abstract: Spiking neural networks have made breakthroughs in computer vision by lending themselves to neuromorphic hardware. However, the neuromorphic hardware lacks parallelism and hence, limits the throughput and hardware acceleration of SNNs on edge devices. To address this problem, many systolic-array SNN accelerators (systolicSNNs) have been proposed recently, but their reliability is still a major concern. In this paper, we first extensively analyze the impact of permanent faults on the SystolicSNNs. Then, we present a novel fault mitigation method, i.e., fault-aware threshold voltage optimization in retraining (FalVolt). FalVolt optimizes the threshold voltage for each layer in retraining to achieve the classification accuracy close to the baseline in the presence of faults. To demonstrate the effectiveness of our proposed mitigation, we classify both static (i.e., MNIST) and neuromorphic datasets (i.e., N-MNIST and DVS Gesture) on a 256x256 systolicSNN with stuck-at faults. We empirically show that the classification accuracy of a systolicSNN drops significantly even at extremely low fault rates (as low as 0.012\%). Our proposed FalVolt mitigation method improves the performance of systolicSNNs by enabling them to operate at fault rates of up to 60\%, with a negligible drop in classification accuracy (as low as 0.1\%). Our results show that FalVolt is 2x faster compared to other state-of-the-art techniques common in artificial neural networks (ANNs), such as fault-aware pruning and retraining without threshold voltage optimization.


title: Polarimetric Reverberation Mapping in Medium-Band Filters

abstract: Earlier, we suggested the "reload" concept of the polarimetric reverberation mapping of active galactic nuclei (AGN), proposed for the first time more than 10 years ago. We have successfully tested this approach of reverberation mapping of the broad emission line on the galaxy Mrk 6. It was shown that such an idea allows one to look at the AGN central parsec structure literally in a new light. However, the method originally assumed the use of spectropolarimetric observations, expensive in terms of telescope time, and implemented on rare large telescopes. Currently, we propose an adaptation of the polarimetric reverberation mapping of broad lines in medium-band filters following the idea of the photometric reverberation mapping, when filters are selected so that their bandwidth is oriented to the broad line and the surrounding continuum near. In this paper, we present the progress status of such monitoring conducted jointly at the Special astrophysical observatory and Asiago Cima Ekar observatory (OAPd/INAF) with support from Rozhen National Astronomical Observatory (NAO), some first results for the most frequently observed AGNs Mrk 335, Mrk 509, and Mrk 817, and the discussion of the future perspectives of the campaign.


title: The non-linear perturbation of a black hole by gravitational waves. III. Newman-Penrose constants

abstract: In this paper we continue our study of the non-linear response of a Schwarzschild black hole to an ingoing gravitational wave by computing the Newman-Penrose (NP) constants. The NP constants are five complex, supertranslation-invariant quantities defined on null infinity $\mathscr{I}^+$ and although put forward in the 60's, they have never been computed in a non-stationary setting. We accomplish this through a numerical implementation of Friedrich's generalized conformal field equations whose semi-global evolution yields direct access to $\mathscr{I}+$. Generalizations of the NP constants' integral expressions are made to allow their computation in a more general gauge that better suits the output of a numerical evolution. Canonical methods of fixing inherent degrees of freedom in their definitions are discussed. The NP constants are then computed for a variety of different ingoing wave profiles in axisymmetry, and then with no symmetry assumptions in 3+1 for which all five are non-zero.


title: Globular Clusters in NGC 4839 Falling into Coma: Evidence for the Second Infall?

abstract: NGC 4839 is the brightest galaxy (cD) of the NGC 4839 group at $R\approx 1$ Mpc in the south-west of the Coma cluster, which is known to be falling into Coma. However, it has been controversial whether it is in the first phase of infall or in the second phase of infall after passing the Coma center. We present a wide field study of globular clusters (GCs) in NGC 4839 and its environment based on Hyper Suprime-Cam $gr$ images in the Subaru archive. We compare the GC system of NGC 4839 with that of NGC 4816, which is the brightest member (S0) of the nearby group and lies at a similar distance in the west from the Coma center. Interestingly the spatial distribution of the GCs in NGC 4839 is significantly more compact than that of the GCs in NGC 4816. In addition, the radial number density profile of the GCs in NGC 4839 shows an abrupt drop at $R_{N4839}\approx 80$ kpc, while that of the GCs in NGC 4816 shows a continuous slow decline even in the outer region at $80<R_{N4816}<500$ kpc. The effective radius of the NGC 4839 GC system is about three times smaller than that of the NGC 4816 GC system. This striking difference can be explained if NGC 4839 lost a significant fraction of the GCs in its outskirt when it passed through Coma. This supports strongly the second infall scenario where the NGC 4839 passed the Coma center about 1.6 Gyr ago, and began the second infall after reaching the apocenter in the south-west recently.


title: On a stratification of positive scalar curvature compact manifolds

abstract: For a compact PSC Riemannian $n$-manifold $(M,g)$, the metric constant $\mathrm {Riem}(g)\in (0, \binom{n}{2}]$ is defined to be the infinimum over $M$ of the spectral scalar curvature $\frac{\sum_{i=1}^N\lambda_i}{\lambda_{\rm max}}$ of $g$, where $\lambda_1, ...,\lambda_N$ are the eigenvalues of the curvature operator of $g$ and $\lambda_{\rm max}$ is the maximal eigenvalue. The functional $g\to \mathrm {Riem}(g)$ is continuous, re-scale invariant and defines a stratification of the space of PSC metrics over $M$. We introduce as well the smooth constant $\mathbf {Riem}(M)\in (0, \binom{n}{2}]$, which is the supremum of $\mathrm {Riem}(g)$ over the set of all psc Riemannian metrics $g$ on $M$. \\ In this paper, we show that in the top layer, compact manifolds with $\mathbf{Riem}=\binom{n}{2}$ are positive space forms. No manifolds have their $\mathbf{Riem}$ in the interval $(\binom{n}{2}-2, \binom{n}{2})$. The manifold $S^{n-1}\times S^1$ and arbitrary connected sums of copies of it with connected sums of positive space forms all have $\mathbf{Riem}=\binom{n-1}{2}$. For $1\leq p\leq n-2\leq 5$, we prove that the manifolds $S^{n-p}\times T^p$ take the intermediate values $\mathbf {Riem}=\binom{n-p}{2}$. From the bottom, we prove that simply connected (resp. $2$-connected, $3$-connected and non-string) compact manifolds of dimension $\geq 5$ (resp. $\geq 7$, $\geq 9$) have $\mathbf{Riem}\geq 1$ (resp. $\geq 3$, $\geq 6$). The proof of these last three results is based on surgery. In fact, we prove that the smooth $\mathbf{Riem}$ constant doesn't decrease after a surgery on the manifold with adequate codimension.


title: Observation of terahertz second harmonic generation from Dirac surface states in the topological insulator Bi$_2$Se$_3$

abstract: We report the observation of second harmonic generation with high conversion efficiency $\sim 0.005\%$ in the terahertz regime from thin films of the topological insulator Bi$_2$Se$_3$ that exhibit the linear photogalvanic effect, measured via time-domain terahertz spectroscopy and terahertz emission, respectively. As neither phenomena is observable from topologically trivial In-doped Bi$_2$Se$_3$, and since no enhancement is observed when subject to band bending, the efficient thickness-independent nonliear responses are attributable to the Dirac fermions of topological surface states of Bi$_2$Se$_3$. This observation of intrinsic terahertz second harmonic generation in an equilibrium system unlocks the full suite of both even and odd harmonic orders in the terahertz regime and opens new pathways to probing quantum geometry via intraband nonlinear processes. We hope our work will motivate the theoretical development of a full treatment of second harmonic generation for probing the quantum geometry in various inversion-breaking topological and twisted materials.


title: Inaccessible Neural Language Models Could Reinvigorate Linguistic Nativism

abstract: Large Language Models (LLMs) have been making big waves in the machine learning community within the past few years. The impressive scalability of LLMs due to the advent of deep learning can be seen as a continuation of empiricist lingusitic methods, as opposed to rule-based linguistic methods that are grounded in a nativist perspective. Current LLMs are generally inaccessible to resource-constrained researchers, due to a variety of factors including closed source code. This work argues that this lack of accessibility could instill a nativist bias in researchers new to computational linguistics, given that new researchers may only have rule-based, nativist approaches to study to produce new work. Also, given that there are numerous critics of deep learning claiming that LLMs and related methods may soon lose their relevancy, we speculate that such an event could trigger a new wave of nativism in the language processing community. To prevent such a dramatic shift and placing favor in hybrid methods of rules and deep learning, we call upon researchers to open source their LLM code wherever possible to allow both empircist and hybrid approaches to remain accessible.


title: Improving the speed of variational quantum algorithms for quantum error correction

abstract: We consider the problem of devising a suitable Quantum Error Correction (QEC) procedures for a generic quantum noise acting on a quantum circuit. In general, there is no analytic universal procedure to obtain the encoding and correction unitary gates, and the problem is even harder if the noise is unknown and has to be reconstructed. The existing procedures rely on Variational Quantum Algorithms (VQAs) and are very difficult to train since the size of the gradient of the cost function decays exponentially with the number of qubits. We address this problem using a cost function based on the Quantum Wasserstein distance of order 1 ($QW_1$). At variance with other quantum distances typically adopted in quantum information processing, $QW_1$ lacks the unitary invariance property which makes it a suitable tool to avoid to get trapped in local minima. Focusing on a simple noise model for which an exact QEC solution is known and can be used as a theoretical benchmark, we run a series of numerical tests that show how, guiding the VQA search through the $QW_1$, can indeed significantly increase both the probability of a successful training and the fidelity of the recovered state, with respect to the results one obtains when using conventional approaches.


title: Convergence for Complex Gaussian Multiplicative Chaos on phase boundaries

abstract: The complex Gaussian Multiplicative Chaos (or complex GMC) is informally defined as a random measure $e^{\gamma X} \mathrm{d} x$ where $X$ is a log correlated Gaussian field on $\mathbb R^d$ and $\gamma=\alpha+i\beta$ is a complex parameter. The correlation function of $X$ is of the form $$ K(x,y)= \log \frac{1}{|x-y|}+ L(x,y),$$ where $L$ is a continuous function. In the present paper, we consider the cases $\gamma\in \mathcal P_{\mathrm{I/II}}$ and $\gamma\in \mathcal{P}'_{\mathrm{II/III}}$ where $$ \mathcal P_{\mathrm{I/II}}:= \{ \alpha+i \beta \ : \alpha,\beta \in \mathbb R \ ; |\alpha|>|\beta| \ ; \ |\alpha|+|\beta|=\sqrt{2d} \}, $$ and $$ \mathcal{P}'_{\mathrm{II/III}}:= \{ \alpha+i \beta \ : \alpha,\beta \in \mathbb R \ ; \ |\alpha|= \sqrt{d/2} \ ; \ |\beta|>\sqrt{2d} \},$$ We prove that if $X$ is replaced by an approximation $X_\epsilon$ obtained via mollification, then $e^{\gamma X_\epsilon} \mathrm{d} x$, when properly rescaled, converges when $\epsilon\to 0$. The limit does not depend on the mollification kernel. When $\gamma\in \mathcal P_{\mathrm{I/II}}$, the convergence holds in probability and in $L^p$ for some value of $p\in [1,\sqrt{2d}/\alpha)$. When $\gamma\in \mathcal{P}'_{\mathrm{II/III}}$ the convergence holds only in law. In this latter case, the limit can be described a complex Gaussian white noise with a random intensity given by a critical real GMC. The regions $\mathcal P_{\mathrm{I/II}}$ and $ \mathcal{P}'_{\mathrm{II/III}}$ correspond to phase boundary between the three different regions of the complex GMC phase diagram. These results complete previous results obtained for the GMC in phase I and III and only leave as an open problem the question of convergence in phase II.


title: Approximate Balancing Weights for Clustered Observational Study Designs

abstract: In a clustered observational study, a treatment is assigned to groups and all units within the group are exposed to the treatment. We develop a new method for statistical adjustment in clustered observational studies using approximate balancing weights, a generalization of inverse propensity score weights that solve a convex optimization problem to find a set of weights that directly minimize a measure of covariate imbalance, subject to an additional penalty on the variance of the weights. We tailor the approximate balancing weights optimization problem to both adjustment sets by deriving an upper bound on the mean square error for each case and finding weights that minimize this upper bound, linking the level of covariate balance to a bound on the bias. We implement the procedure by specializing the bound to a random cluster-level effects model, leading to a variance penalty that incorporates the signal signal-to-noise ratio and penalizes the weight on individuals and the total weight on groups differently according to the the intra-class correlation.


title: Waring identifiability for powers of forms via degenerations

abstract: We discuss an approach to the secant non-defectivity of the varieties parametrizing $k$-th powers of forms of degree $d$. It employs a Terracini type argument along with certain degeneration arguments, some of which are based on toric geometry. This implies a result on the identifiability of the Waring decompositions of general forms of degree kd as a sum of $k$-th powers of degree $d$ forms, for which an upper bound on the Waring rank was proposed by Fr\"oberg, Ottaviani and Shapiro.


title: DriCon: On-device Just-in-Time Context Characterization for Unexpected Driving Events

abstract: Driving is a complex task carried out under the influence of diverse spatial objects and their temporal interactions. Therefore, a sudden fluctuation in driving behavior can be due to either a lack of driving skill or the effect of various on-road spatial factors such as pedestrian movements, peer vehicles' actions, etc. Therefore, understanding the context behind a degraded driving behavior just-in-time is necessary to ensure on-road safety. In this paper, we develop a system called \ourmethod{} that exploits the information acquired from a dashboard-mounted edge-device to understand the context in terms of micro-events from a diverse set of on-road spatial factors and in-vehicle driving maneuvers taken. \ourmethod{} uses the live in-house testbed and the largest publicly available driving dataset to generate human interpretable explanations against the unexpected driving events. Also, it provides a better insight with an improved similarity of $80$\% over $50$ hours of driving data than the existing driving behavior characterization techniques.


title: Mixed volumes of normal complexes

abstract: Normal complexes are orthogonal truncations of polyhedral fans. In this paper, we develop the study of mixed volumes for normal complexes. Our main result is a sufficiency condition that ensures when the mixed volumes of normal complexes associated to a given fan satisfy the Alexandrov-Fenchel inequalities. By specializing to Bergman fans of matroids, we give a new proof of the Heron-Rota-Welsh Conjecture as a consequence of the Alexandrov-Fenchel inequalities for normal complexes.


title: Characterization of Fast Ion Transport via Position-Dependent Optical Deshelving

abstract: Ion transport is an essential operation in some models of quantum information processing, where fast ion shuttling with minimal motional excitation is necessary for efficient, high-fidelity quantum logic. While fast and cold ion shuttling has been demonstrated, the dynamics and specific trajectory of an ion during diabatic transport have not been studied in detail. Here we describe a position-dependent optical deshelving technique useful for sampling an ion's position throughout its trajectory, and we demonstrate the technique on fast linear transport of a $^{40}\text{Ca}^+$ ion in a surface-electrode ion trap. At high speed, the trap's electrode filters strongly distort the transport potential waveform. With this technique, we observe deviations from the intended constant-velocity (100 m/s) transport: we measure an average speed of 83(2) m/s and a peak speed of 251(6) m/s over a distance of 120 $\mu$m


title: Pointwise Bi-Slant Submanifolds in Locally Conformal K\"ahler Manifolds Immersed as Warped Products

abstract: We study immersions of pointwise bi-slant submanifolds of locally conformal K\"ahler manifolds as warped products. In particular, we establish characterisation theorem for a pointwise bi-slant submanifold of a locally conformal K\"ahler manifold to be immersed as a warped product and show that a necessary condition is that the Lee vector field $B$ is orthogonal to the second factor and the warping function $\lambda$ satisfies $\text{grad}(\ln\lambda)=\frac{1}{2}B^T$, where $B^T$ denotes the tangential part of the Lee vector field. We also extend Chen's inequality for the squared length of the second fundamental form to our case and study the corresponding equality case.


title: Time Domain Verification of Differential Transmission Line Modeling Methods

abstract: The advantages and limitations of time-domain pseudo-random binary sequence (PRBS) excitation methods for system identification of individual modes within a multi-conductor transmission system are discussed. We develop the modifications necessary to standard frequency-domain transmission-line models to match time-domain experimental data from several types of transmission systems. We show a variety of experimental results showing very good to excellent agreement with our model's predictions, up to approximately 10 GHz.


title: Non-centrosymmetric Sr$_{2}$IrO$_{4}$ obtained under High Pressure

abstract: Sr$_{2}$IrO$_{4}$ with strong spin-orbit coupling (SOC) and Hubbard repulsion (U) hosts Mott insulating states. The similar crystal structure, magnetic and electronic properties, particularly the $d$-wave gap observed in Sr$_{2}$IrO$_{4}$ enhanced the analogies to cuprate high-$T_{c}$ superconductor, La$_{2}$CuO$_{4}$. The incomplete analogy was due to the lack of broken inversion symmetry phases observed in Sr$_{2}$IrO$_{4}$. Here, under high pressure and high temperature conditions, we report a non-centrosymmetric Sr$_{2}$IrO$_{4}$. The crystal structure and its noncentrosymmetric character were determined by single crystal X-ray diffraction and high-resolution scanning transmission electron microscopy (HR-STEM). The magnetic characterization confirms the Ir$^{4+}$ with $S$ = 1/2 at low temperature in Sr$_{2}$IrO$_{4}$ with magnetic ordering occurred at around 86 K, where a larger moment is observed than the ambient pressure Sr$_{2}$IrO$_{4}$. Moreover, the resistivity measurement shows three-dimensional Mott variable-range hopping existed in the system. This non-centrosymmetric Sr$_{2}$IrO$_{4}$ phase appears to be a unique material to offer further understanding of high-$T_{c}$ superconductivity.


title: Integrable systems with linear periodic integral for the Lie algebra $\mathrm{e}(3)$

abstract: Integrable systems with a linear periodic integral for the Lie algebra $\mathrm{e}(3)$ are considered. One investigates singulariries of the Liouville foliation, bifurcation diagram of the momentum mapping, transformations of Liouville tori, topology of isoenergy surfaces and other topological properties of such systems.


title: Numerical study of the rate of convergence of Chernoff approximations to solutions of the heat equation with full list of illustrations and Python source code

abstract: Chernoff approximations are a flexible and powerful tool of functional analysis, which can be used, in particular, to find numerically approximate solutions of some differential equations with variable coefficients. For many classes of equations such approximations have already been constructed since pioneering papers of Prof. O.G.Somlyanov in 2000, however, the speed of their convergence to the exact solution has not been properly studied. We select the heat equation (because its exact solutions are already known) as a simple yet informative model example for the study of the rate of convergence of Chernoff approximations. Examples illustrating the rate of convergence of Chernoff approximations to the solution of the Cauchy problem for the heat equation are constructed in the paper. Numerically we show that for initial conditions that are smooth enough the order of approximation is equal to the order of Chernoff tangency of the Chernoff function used. We also consider not smooth enough initial conditions and show how H\"older class of initial condition is related to the rate of convergence. This method of study in the future can be applied to general second order parabolic equation with variable coefficients by a slight modification of our Python 3 code. This arXiv version of the text is a supplementary material for our journal article. Here we include all the written text from the article and additionally all illustrations (Appendix A) and full text of the Python 3 code (Appendix B).


title: Laser Inter-Satellite Link Setup Delay: Quantification, Impact, and Tolerable Value

abstract: Dynamic laser inter-satellite links (LISLs) provide the flexibility of connecting a pair of satellites as required (dynamically) while static LISLs need to be active continuously between the energy-constrained satellites. However, due to the LISL establishment time (termed herein as LISL setup delay) being in the order of seconds, realizing dynamic LISLs is currently unfeasible. Towards the realization of dynamic LISLs, we first study the quantification of LISL setup delay; then we calculate the end-to-end latency of a free-space optical satellite network (FSOSN) with the LISL setup delay; subsequently, we analyze the impact of LISL setup delay on the end-to-end latency of the FSOSN. We also provide design guidelines for the laser communication terminal manufacturers in the form of maximum tolerable value of LISL setup delay for which the FSOSN based on Starlink's Phase I satellite constellation will be meaningful to use for low-latency long-distance inter-continental data communications.


title: Demonstration of an AI-driven workflow for autonomous high-resolution scanning microscopy

abstract: With the continuing advances in scientific instrumentation, scanning microscopes are now able to image physical systems with up to sub-atomic-level spatial resolutions and sub-picosecond time resolutions. Commensurately, they are generating ever-increasing volumes of data, storing and analysis of which is becoming an increasingly difficult prospect. One approach to address this challenge is through self-driving experimentation techniques that can actively analyze the data being collected and use this information to make on-the-fly measurement choices, such that the data collected is sparse but representative of the sample and sufficiently informative. Here, we report the Fast Autonomous Scanning Toolkit (FAST) that combines a trained neural network, a route optimization technique, and efficient hardware control methods to enable a self-driving scanning microscopy experiment. The key features of our method are that: it does not require any prior information about the sample, it has a very low computational cost, and that it uses generic hardware controls with minimal experiment-specific wrapping. We test this toolkit in numerical experiments and a scanning dark-field x-ray microscopy experiment of a $WSe_2$ thin film, where our experiments show that a FAST scan of <25% of the sample is sufficient to produce both a high-fidelity image and a quantitative analysis of the surface distortions in the sample. We show that FAST can autonomously identify all features of interest in the sample while significantly reducing the scan time, the volume of data acquired, and dose on the sample. The FAST toolkit is easy to apply for any scanning microscopy modalities and we anticipate adoption of this technique will empower broader multi-level studies of the evolution of physical phenomena with respect to time, temperature, or other experimental parameters.


title: Model selection in atomistic simulation

abstract: There are many atomistic simulation methods with very different costs, accuracies, transferabilities, and numbers of empirical parameters. I show how statistical model selection can compare these methods fairly, even when they are very different. These comparisons are also useful for developing new methods that balance cost and accuracy. As an example, I build a semiempirical model for hydrogen clusters.


title: An Approach to Stochastic Dynamic Games with Asymmetric Information and Hidden Actions

abstract: We consider in discrete time, a general class of sequential stochastic dynamic games with asymmetric information with the following features. The underlying system has Markovian dynamics controlled by the agents' joint actions. Each agent's instantaneous utility depends on the current system state and the agents' joint actions. At each time instant each agent makes a private noisy observation of the current system state and the agents' actions in the previous time instant. In addition, at each time instant all agents have a common noisy observation of the current system state and their actions in the previous time instant. Each agent's actions are part of his private information. The objective is to determine Bayesian Nash Equilibrium (BNE) strategy profiles that are based on a compressed version of the agents' information and can be sequentially computed; such BNE strategy profiles may not always exist. We present an approach/methodology that achieves the above-stated objective, along with an instance of a game where BNE strategy profiles with the above-mentioned characteristics exist. We show that the methodology also works for the case where the agents have no common observations.


title: The covariance metric in the Blaschke locus

abstract: We prove that the Blaschke locus has the structure of a finite dimensional smooth manifold away from the Teichm{\"u}ller space and study its Riemannian manifold structure with respect to the covariance metric introduced by Guillarmou, Knieper and Lefeuvre in \cite{GeodesicStretch}. We also identify some families of geodesics in the Blaschke locus arising from Hitchin representations for orbifolds and show that they have infinite length with respect to the covariance metric.


title: K$^0_\mathrm{S}$ and $\Lambda(\overline\Lambda)$ two-particle femtoscopic correlations in PbPb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

abstract: Two-particle correlations are presented for K$^0_\mathrm{S}$, $\Lambda$, and $\overline\Lambda$ strange hadrons as a function of relative momentum in lead-lead collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The dataset corresponds to an integrated luminosity of 0.607 nb$^{-1}$ and was collected using the CMS detector at the CERN LHC. These correlations are sensitive to quantum statistics and to final-state interactions between the particles. The source size extracted from the K$^0_\mathrm{S}$K$^0_\mathrm{S}$ correlations is found to decrease from 4 to 1 fm in going from central to peripheral collisions. Strong interaction scattering parameters (i.e., scattering length and effective range) are determined from the $\Lambda$K$^0_\mathrm{S}$ and $\Lambda\Lambda$ (including their charge conjugates) correlations using the Lednick$\'y$-Lyuboshitz model and are compared to theoretical and other experimental results.


title: Scalable synthesis and characterization of multilayer $\gamma$-graphyne, new carbon crystals with a small direct bandgap

abstract: $\gamma$-Graphyne is the most symmetric sp2/sp1 allotrope of carbon, which can be viewed as graphene uniformly expanded through insertion of two-carbon acetylenic units between all the aromatic rings. To date, synthesis of bulk $\gamma$-graphyne has remained a challenge. We here report the synthesis of multilayer $\gamma$-graphyne through crystallization-assisted irreversible cross-coupling polymerization. Comprehensive characterization of this new carbon phase is described, including synchrotron X-ray diffraction, electron diffraction, lateral force microscopy, Raman and infrared spectroscopy, and cyclic voltammetry. Experiments indicate that $\gamma$-graphyne is a 0.48 eV bandgap semiconductor, with a hexagonal a-axis spacing of 6.88 {\AA} and an interlayer spacing of 3.48 {\AA}, which is consistent with theoretical predictions. The observed crystal structure has an aperiodic sheet stacking. The material is thermally stable up to 240 $^\circ$C but undergoes a transformation at higher temperatures. While conventional 2D polymerizations and reticular chemistry rely on error correction through reversibility, we demonstrate that a periodic covalent lattice can be synthesized under purely kinetic control. The reported methodology is scalable and inspires extension to other allotropes of the graphyne family.


title: A Novel Framework for Handling Sparse Data in Traffic Forecast

abstract: The ever increasing amount of GPS-equipped vehicles provides in real-time valuable traffic information for the roads traversed by the moving vehicles. In this way, a set of sparse and time evolving traffic reports is generated for each road. These time series are a valuable asset in order to forecast the future traffic condition. In this paper we present a deep learning framework that encodes the sparse recent traffic information and forecasts the future traffic condition. Our framework consists of a recurrent part and a decoder. The recurrent part employs an attention mechanism that encodes the traffic reports that are available at a particular time window. The decoder is responsible to forecast the future traffic condition.


title: HTTE: A Hybrid Technique For Travel Time Estimation In Sparse Data Environments

abstract: Travel time estimation is a critical task, useful to many urban applications at the individual citizen and the stakeholder level. This paper presents a novel hybrid algorithm for travel time estimation that leverages historical and sparse real-time trajectory data. Given a path and a departure time we estimate the travel time taking into account the historical information, the real-time trajectory data and the correlations among different road segments. We detect similar road segments using historical trajectories, and use a latent representation to model the similarities. Our experimental evaluation demonstrates the effectiveness of our approach.


title: Learning to Control and Coordinate Mixed Traffic Through Robot Vehicles at Complex and Unsignalized Intersections

abstract: Intersections are essential road infrastructures for traffic in modern metropolises. However, they can also be the bottleneck of traffic flows as a result of traffic incidents or the absence of traffic coordination mechanisms such as traffic lights. Recently, various control and coordination mechanisms that are beyond traditional control methods have been proposed to improve the efficiency of intersection traffic. Amongst these methods, the control of foreseeable mixed traffic that consists of human-driven vehicles (HVs) and robot vehicles (RVs) has emerged. In this project, we propose a decentralized multi-agent reinforcement learning approach for the control and coordination of mixed traffic at real-world, complex intersections--a topic that has not been previously explored. Comprehensive experiments are conducted to show the effectiveness of our approach. In particular, we show that using 5% RVs, we can prevent congestion formation inside a complex intersection under the actual traffic demand of 700 vehicles per hour. In contrast, without RVs, congestion starts to develop when the traffic demand reaches as low as 200 vehicles per hour. When there exist more than 60% RVs in traffic, our method starts to achieve comparable or even better performance to traffic signals on the average waiting time of all vehicles at the intersection. Our method is also robust against both blackout events and sudden RV percentage drops, and enjoys excellent generalizablility, which is illustrated by its successful deployment in two unseen intersections.


title: Rock Guitar Tablature Generation via Natural Language Processing

abstract: Deep learning has recently empowered and democratized generative modeling of images and text, with additional concurrent works exploring the possibility of generating more complex forms of data, such as audio. However, the high dimensionality, long-range dependencies, and lack of standardized datasets currently makes generative modeling of audio and music very challenging. We propose to model music as a series of discrete notes upon which we can use autoregressive natural language processing techniques for successful generative modeling. While previous works used similar pipelines on data such as sheet music and MIDI, we aim to extend such approaches to the under-studied medium of guitar tablature. Specifically, we develop the first work to our knowledge that models one specific genre as guitar tablature: heavy rock. Unlike other works in guitar tablature generation, we have a freely available public demo at https://huggingface.co/spaces/josuelmet/Metal_Music_Interpolator


title: An efficient hybrid classification approach for COVID-19 based on Harris Hawks Optimization and Salp Swarm Optimization

abstract: Feature selection can be defined as one of the pre-processing steps that decrease the dimensionality of a dataset by identifying the most significant attributes while also boosting the accuracy of classification. For solving feature selection problems, this study presents a hybrid binary version of the Harris Hawks Optimization algorithm (HHO) and Salp Swarm Optimization (SSA) (HHOSSA) for Covid-19 classification. The proposed (HHOSSA) presents a strategy for improving the basic HHO's performance using the Salp algorithm's power to select the best fitness values. The HHOSSA was tested against two well-known optimization algorithms, the Whale Optimization Algorithm (WOA) and the Grey wolf optimizer (GWO), utilizing a total of 800 chest X-ray images. A total of four performance metrics (Accuracy, Recall, Precision, F1) were employed in the studies using three classifiers (Support vector machines (SVMs), k-Nearest Neighbor (KNN), and Extreme Gradient Boosting (XGBoost)). The proposed algorithm (HHOSSA) achieved 96% accuracy with the SVM classifier, and 98% accuracy with two classifiers, XGboost and KNN.


title: Towards Dependable Autonomous Systems Based on Bayesian Deep Learning Components

abstract: As autonomous systems increasingly rely on Deep Neural Networks (DNN) to implement the navigation pipeline functions, uncertainty estimation methods have become paramount for estimating confidence in DNN predictions. Bayesian Deep Learning (BDL) offers a principled approach to model uncertainties in DNNs. However, in DNN-based systems, not all the components use uncertainty estimation methods and typically ignore the uncertainty propagation between them. This paper provides a method that considers the uncertainty and the interaction between BDL components to capture the overall system uncertainty. We study the effect of uncertainty propagation in a BDL-based system for autonomous aerial navigation. Experiments show that our approach allows us to capture useful uncertainty estimates while slightly improving the system's performance in its final task. In addition, we discuss the benefits, challenges, and implications of adopting BDL to build dependable autonomous systems.


title: Open Case Studies: Statistics and Data Science Education through Real-World Applications

abstract: With unprecedented and growing interest in data science education, there are limited educator materials that provide meaningful opportunities for learners to practice statistical thinking, as defined by Wild and Pfannkuch (1999), with messy data addressing real-world challenges. As a solution, Nolan and Speed (1999) advocated for bringing applications to the forefront in undergraduate statistics curriculum with the use of in-depth case studies to encourage and develop statistical thinking in the classroom. Limitations to this approach include the significant time investment required to develop a case study -- namely, to select a motivating question and to create an illustrative data analysis -- and the domain expertise needed. As a result, case studies based on realistic challenges, not toy examples, are scarce. To address this, we developed the Open Case Studies (https://www.opencasestudies.org) project, which offers a new statistical and data science education case study model. This educational resource provides self-contained, multimodal, peer-reviewed, and open-source guides (or case studies) from real-world examples for active experiences of complete data analyses. We developed an educator's guide describing how to most effectively use the case studies, how to modify and adapt components of the case studies in the classroom, and how to contribute new case studies. (https://www.opencasestudies.org/OCS_Guide).


title: Grey area in Embedded WMLES on a transonic nacelle-aircraft configuration

abstract: A scale resolving hybrid RANS-LES technique is applied to an aircraft-nacelle configuration under transonic flow conditions using the unstructured, compressible TAU solver. Therefore, a wall modelled LES methodology is locally applied to the nacelle lower surface in order to examine shock induced separation. In this context a synthetic turbulence generator (STG) is used to shorten the adaption region at the RANS-LES interface. Prior to the actual examinations, fundamental features of the simulation technique are validated by simulations of decaying isotropic turbulence as well as a flat plate flow. For the aircraft-nacelle configuration at a Reynolds number of 3.3 million a sophisticated mesh with 420 million points was designed which refines 32 % of the outer casing surface of the nacelle. The results show a development of a well resolved turbulent boundary layer with a broad spectrum of turbulent scales which demonstrates the applicability of the mesh and method for aircraft configurations. Furthermore, the necessity of a low dissipation low dispersion scheme is demonstrated. However, the distinct adaption region downstream of the STG limits the employment of the method in case of shock buffet for the given flow conditions.


title: Deep Reinforcement Learning for Asset Allocation: Reward Clipping

abstract: Recently, there are many trials to apply reinforcement learning in asset allocation for earning more stable profits. In this paper, we compare performance between several reinforcement learning algorithms - actor-only, actor-critic and PPO models. Furthermore, we analyze each models' character and then introduce the advanced algorithm, so called Reward clipping model. It seems that the Reward Clipping model is better than other existing models in finance domain, especially portfolio optimization - it has strength both in bull and bear markets. Finally, we compare the performance for these models with traditional investment strategies during decreasing and increasing markets.


title: A Minimal Formulation of Session Types

abstract: Session types are a type-based approach to the verification of message-passing programs. They specify communication structures essential for program correctness; a session type says what and when should be exchanged through a channel. Central to session-typed languages are sequencing constructs in types and processes that explicitly specify the order of actions in a protocol. In this paper we study session types without sequencing. The resulting framework of minimal session types is arguably the simplest form of session types one could conceive. In the context of a core process calculus with sessions and higher-order concurrency (abstraction-passing), we establish two main technical results. First, we prove that every process $P$ typable with standard session types can be compiled down into a process $\mathcal{D}(P)$ typable with minimal session types. Second, we prove that $P$ and $\mathcal{D}(P)$ are behaviorally equivalent. These results indicate that having sequencing constructs in processes and session types is convenient but redundant: only sequentiality in processes is truly indispensable, as it can correctly codify sequentiality in types. Our developments draw inspiration from work by Parrow on behavior-preserving decompositions of untyped processes. By casting Parrow's results in the realm of typed processes, our developments reveal a conceptually simple formulation of session types and a principled avenue to the integration of session types into programming languages without sequencing in types.


title: Wavenumber Scattering and Inter-band Targeted Energy Transfer in Phononic Lattices with Local Vibro-Impact Nonlinearities

abstract: We propose a method for manipulating wave propagation in phononic lattices by employing local vibro-impact (VI) nonlinearities to \textit{scatter} energy across the underling linear band structure of the lattice, and \textit{transfer} energy from lower to higher optical bands.Inspired by recent developments in the field of nonlinear targeted energy transfer (TET) using \textit{non-resonant} energy exchanges, we achieve this using spatially localized VI forces that redistribute energy across the linear spectrum of the lattice in a non-resonant fashion.First, a 1-dimensional (1D), 2-band phononic lattice with embedded VI unit cells is computationally studied to demonstrate that energy is scattered in the wavenumber domain, and this nonlinear scattering mechanism depends on the energy of the propagating wave. Next, a 4-band lattice is studied with a similar technique to demonstrate the concept of inter-band targeted energy transfer (IBTET) and to establish analogous scaling relations with respect to energy. To interpret the results of IBTET, we study the nonlinear normal modes (NNMs) of a reduced order model (ROM) of the VI unit cell in the 4-band lattice, using the method of numerical continuation. Interestingly, the slope of the frequency-energy branches of the ROM corresponding to the 1:1 resonance NNM matches remarkably well with the dependence of IBTET to input energy in the 4-band lattice. Moreover, relations between the dynamics of the VI lattice and the NNMs of the underlying Hamiltonian system provide physical interpretations for the relative energy transfers.


title: Probabilistic Constraint Construction for Network-safe Load Coordination

abstract: Distributed Energy Resources (DERs) can provide balancing services to the grid, but their power variations might cause voltage and current constraint violations in the distribution network, compromising network safety. This could be avoided by including network constraints within DER control formulations, but the entities coordinating DERs (e.g., aggregators) may not have access to network information, which typically is known only to the utility. Therefore, it is challenging to develop network-safe DER control algorithms when the aggregator is not the utility; it requires these entities to coordinate with each other. In this paper, we develop an aggregator-utility coordination framework that enables network-safe control of thermostatically-controlled loads to provide frequency regulation. In our framework, the utility sends a network-safe constraint set on the aggregator's command without directly sharing any network information. We propose a constraint set construction algorithm that guarantees satisfaction of a chance constraint on network safety. Assuming monotonicity of the probability of network safety with respect to the aggregator's command, we leverage the bisection method to find the largest possible constraint set, providing maximum flexibility to the aggregator. Simulations show that, compared to two benchmark algorithms, the proposed approach provides a good balance between service quality and network safety.


title: A characterization of the $L^2$-range of the Poisson transforms on a class of vector bundles over the quaternionic hyperbolic spaces

abstract: We study the $L^2$-boundedness of the Poisson transforms associated to the homogeneous vector bundles $ Sp(n,1)\times_{Sp(n)\times Sp(1)} V_\tau$ over the quaternionic hyperbolic spaces $ Sp(n,1)/Sp(n)\times Sp(1)$ associated with irreducible representations $\tau$ of $ Sp(n)\times Sp(1)$ which are trivial on $ Sp(n)$. As a consequence, we describe the image of the section space $L^2(Sp(n,1)\times_{Sp(n)\times Sp(1)} V_\tau)$ under the generalized spectral projections associated to a family of eigensections of the Casimir operator.


title: Reinforcement Learning-based Joint Handover and Beam Tracking in Millimeter-wave Networks

abstract: In this paper, we develop an algorithm for joint handover and beam tracking in millimeter-wave (mmWave) networks. The aim is to provide a reliable connection in terms of the achieved throughput along the trajectory of the mobile user while preventing frequent handovers. We model the association problem as an optimization problem and propose a reinforcement learning-based solution. Our approach learns whether and when beam tracking and handover should be performed and chooses the target base stations. In the case of beam tracking, we propose a tracking algorithm based on measuring a small spatial neighbourhood of the optimal beams in the previous time slot. Simulation results in an outdoor environment show the superior performance of our proposed solution in achievable throughput and the number of handovers needed in comparison to a multi-connectivity baseline and a learning-based handover baseline.


title: A Simultaneous Dual-Frequency Scintillation Arc Survey of Six Bright Canonical Pulsars Using the Upgraded Giant Metrewave Radio Telescope

abstract: We use the upgraded Giant Metrewave Radio Telescope to measure scintillation arc properties in six bright canonical pulsars with simultaneous dual frequency coverage. These observations at frequencies from 300 to 750 MHz allowed for detailed analysis of arc evolution across frequency and epoch. We perform more robust determinations of frequency dependence for arc curvature, scintillation bandwidth, and scintillation timescale, and comparison between arc curvature and pseudo-curvature than allowed by single-frequency-band-per-epoch measurements, which we find to agree with theory and previous literature. We find a strong correlation between arc asymmetry and arc curvature, which we have replicated using simulations, and attribute to a bias in the Hough transform approach to scintillation arc analysis. Possible evidence for an approximately week long timescale over which a given scattering screen dominates signal propagation was found by tracking visible scintillation arcs in each epoch in PSR J1136+1551. The inclusion of a 155 minute observation allowed us to resolve the scale of scintillation variations on short timescales, which we find to be directly tied to the amount of ISM sampled over the observation. Some of our pulsars showed either consistent or emerging asymmetries in arc curvature, indicating instances of refraction across their lines of sight. Significant features in various pulsars, such as multiple scintillation arcs in PSR J1136+1551 and flat arclets in PSR J1509+5531, that have been found in previous works, were also detected. The simultaneous multiple band observing capability of the upgraded GMRT shows excellent promise for future pulsar scintillation work.


title: Partial entropy decomposition reveals higher-order structures in human brain activity

abstract: The standard approach to modeling the human brain as a complex system is with a network, where the basic unit of interaction is a pairwise link between two brain regions. While powerful, this approach is limited by the inability to assess higher-order interactions involving three or more elements directly. In this work, we present a method for capturing higher-order dependencies in discrete data based on partial entropy decomposition (PED). Our approach decomposes the joint entropy of the whole system into a set of strictly non-negative partial entropy atoms that describe the redundant, unique, and synergistic interactions that compose the system's structure. We begin by showing how the PED can provide insights into the mathematical structure of both the FC network itself, as well as established measures of higher-order dependency such as the O-information. When applied to resting state fMRI data, we find robust evidence of higher-order synergies that are largely invisible to standard functional connectivity analyses. This synergistic structure distinct from structural features based on redundancy that have previously dominated FC analyses. Our approach can also be localized in time, allowing a frame-by-frame analysis of how the distributions of redundancies and synergies change over the course of a recording. We find that different ensembles of regions can transiently change from being redundancy-dominated to synergy-dominated, and that the temporal pattern is structured in time. These results provide strong evidence that there exists a large space of unexplored structures in human brain data that have been largely missed by a focus on bivariate network connectivity models. This synergistic "shadow structures" is dynamic in time and, likely will illuminate new and interesting links between brain and behavior.


title: Incremental Dead State Detection in Logarithmic Time

abstract: Identifying live and dead states in an abstract transition system is a recurring problem in formal verification; for example, it arises in our recent work on efficiently deciding regex constraints in SMT. However, state-of-the-art graph algorithms for maintaining reachability information incrementally (that is, as states are visited and before the entire state space is explored) assume that new edges can be added from any state at any time, whereas in many applications, outgoing edges are added from each state as it is explored. To formalize the latter situation, we propose guided incremental digraphs (GIDs), incremental graphs which support labeling closed states (states which will not receive further outgoing edges). Our main result is that dead state detection in GIDs is solvable in $O(\log m)$ amortized time per edge for $m$ edges, improving upon $O(\sqrt{m})$ per edge due to Bender, Fineman, Gilbert, and Tarjan (BFGT) for general incremental directed graphs. We introduce two algorithms for GIDs: one establishing the logarithmic time bound, and a second algorithm to explore a lazy heuristics-based approach. To enable an apples-to-apples experimental comparison, we implemented both algorithms, two simpler baselines, and the state-of-the-art BFGT baseline using a common directed graph interface in Rust. Our evaluation shows $110$-$530$x speedups over BFGT for the largest input graphs over a range of graph classes, random graphs, and graphs arising from regex benchmarks.


title: Planning Visual Inspection Tours for a 3D Dubins Airplane Model in an Urban Environment

abstract: This paper investigates the problem of planning a minimum-length tour for a three-dimensional Dubins airplane model to visually inspect a series of targets located on the ground or exterior surface of objects in an urban environment. Objects are 2.5D extruded polygons representing buildings or other structures. A visibility volume defines the set of admissible (occlusion-free) viewing locations for each target that satisfy feasible airspace and imaging constraints. The Dubins traveling salesperson problem with neighborhoods (DTSPN) is extended to three dimensions with visibility volumes that are approximated by triangular meshes. Four sampling algorithms are proposed for sampling vehicle configurations within each visibility volume to define vertices of the underlying DTSPN. Additionally, a heuristic approach is proposed to improve computation time by approximating edge costs of the 3D Dubins airplane with a lower bound that is used to solve for a sequence of viewing locations. The viewing locations are then assigned pitch and heading angles based on their relative geometry. The proposed sampling methods and heuristics are compared through a Monte-Carlo experiment that simulates view planning tours over a realistic urban environment.


title: Optimization of Hybrid Power Plants: When Is a Detailed Electrolyzer Model Necessary?

abstract: Hybrid power plants comprising renewable power sources and electrolyzers are envisioned to play a key role in accelerating the transition towards decarbonization. It is common in the current literature to use simplified operational models for electrolyzers. It is still an open question whether this is a good practice, and if not, when a more detailed operational model is necessary. This paper answers it by assessing the impact of adding different levels of electrolyzer details, i.e., physics and operational constraints, to the optimal dispatch problem of a hybrid power plant in the day-ahead time stage. Our focus lies on the number of operating states (on, off, standby) as well as the number of linearization segments used for approximating the non-linear hydrogen production curve. For that, we develop several mixed-integer linear models, each representing a different level of operational details. We conduct a thorough comparative ex-post performance analysis under different price conditions, wind farm capacities, and minimum hydrogen demand requirements, and discuss under which operational circumstances a detailed model is necessary. In particular, we provide a case under which a simplified model, compared to a detailed one, results in a decrease in profit of 1.8% and hydrogen production of 13.5% over a year. The key lesson learned is that a detailed model potentially earns a higher profit in circumstances under which the electrolyzer operates with partial loading. This could be the case for a certain range of electricity and hydrogen prices, or limited wind power availability. The detailed model also provides a better estimation of true hydrogen production, facilitating the logistics required.


title: A Framework for Active Haptic Guidance Using Robotic Haptic Proxies

abstract: Haptic feedback is an important component of creating an immersive mixed reality experience. Traditionally, haptic forces are rendered in response to the user's interactions with the virtual environment. In this work, we explore the idea of rendering haptic forces in a proactive manner, with the explicit intention to influence the user's behavior through compelling haptic forces. To this end, we present a framework for active haptic guidance in mixed reality, using one or more robotic haptic proxies to influence user behavior and deliver a safer and more immersive virtual experience. We provide details on common challenges that need to be overcome when implementing active haptic guidance, and discuss example applications that show how active haptic guidance can be used to influence the user's behavior. Finally, we apply active haptic guidance to a virtual reality navigation problem, and conduct a user study that demonstrates how active haptic guidance creates a safer and more immersive experience for users.


title: Coulomb blockade of chiral Majorana and complex fermions far from equilibrium

abstract: We study charge transport in a single-electron transistor implemented as an interferometer such that the Coulomb blockaded middle island contains a circular chiral Majorana or Dirac edge mode. We concentrate on the regime of small conductance and provide an asymptotic solution in the limit of high transport voltage exceeding the charging energy. The solution is achieved using an instanton-like technique. The distinctions between Majorana and Dirac cases appears when the tunnel junctions are unequal. The main difference is in the offset current at high voltages which can be higher up to $50\%$ in Majorana case. It is caused by an additional particle-hole symmetry of the distribution function in the Majorana case. There is also an eye-catching distinction between the oscillations patterns of the current as a function of the gate charge. We conjecture this distinction survives at lower transport voltages as well.

categories: cond-mat.mes-hall cond-mat.str-el cond-mat.supr-con hep-th


title: Algebraic Invariants of Codes on Weighted Projective Planes

abstract: Weighted projective spaces are natural generalizations of projective spaces with a rich structure. Projective Reed-Muller codes are error-correcting codes that played an important role in reliably transmitting information on digital communication channels. In this case study, we explore the power of commutative and homological algebraic techniques to study weighted projective Reed-Muller (WPRM) codes on weighted projective spaces of the form $\mathbb{P}(1,1,a)$. We compute minimal free resolutions and thereby obtain Hilbert series for the vanishing ideal of the $\mathbb{F}_q$-rational points, and compute main parameters for these codes.


title: Light-controlled multi-phase structuring of perovskite crystal enabled by thermoplasmonic metasurface

abstract: Halide perovskites belong to an important family of semiconducting materials with unique electronic properties that enable a myriad of applications, especially in photovoltaics and optoelectronics. Their optical properties, including photoluminescence quantum yield, are affected and notably enhanced at crystal imperfections where the symmetry is broken and the density of states increases. These lattice distortions can be introduced through structural phase transitions, allowing charge gradients to appear near the interfaces between phase structures. In this work, we demonstrate controlled multi-phase structuring in a single perovskite crystal. The concept uses cesium lead bromine (CsPbBr3) placed on a thermoplasmonic TiN/Si metasurface and enables single, double and triple phase structures to form on demand above the room temperature. This approach opens up application horizons of dynamically controlled heterostructures with distinctive electronic and enhanced optical properties.


title: GH-Feat: Learning Versatile Generative Hierarchical Features from GANs

abstract: Recent years witness the tremendous success of generative adversarial networks (GANs) in synthesizing photo-realistic images. GAN generator learns to compose realistic images and reproduce the real data distribution. Through that, a hierarchical visual feature with multi-level semantics spontaneously emerges. In this work we investigate that such a generative feature learned from image synthesis exhibits great potentials in solving a wide range of computer vision tasks, including both generative ones and more importantly discriminative ones. We first train an encoder by considering the pretrained StyleGAN generator as a learned loss function. The visual features produced by our encoder, termed as Generative Hierarchical Features (GH-Feat), highly align with the layer-wise GAN representations, and hence describe the input image adequately from the reconstruction perspective. Extensive experiments support the versatile transferability of GH-Feat across a range of applications, such as image editing, image processing, image harmonization, face verification, landmark detection, layout prediction, image retrieval, etc. We further show that, through a proper spatial expansion, our developed GH-Feat can also facilitate fine-grained semantic segmentation using only a few annotations. Both qualitative and quantitative results demonstrate the appealing performance of GH-Feat.


title: Traffic Steering for 5G Multi-RAT Deployments using Deep Reinforcement Learning

abstract: In 5G non-standalone mode, traffic steering is a critical technique to take full advantage of 5G new radio while optimizing dual connectivity of 5G and LTE networks in multiple radio access technology (RAT). An intelligent traffic steering mechanism can play an important role to maintain seamless user experience by choosing appropriate RAT (5G or LTE) dynamically for a specific user traffic flow with certain QoS requirements. In this paper, we propose a novel traffic steering mechanism based on Deep Q-learning that can automate traffic steering decisions in a dynamic environment having multiple RATs, and maintain diverse QoS requirements for different traffic classes. The proposed method is compared with two baseline algorithms: a heuristic-based algorithm and Q-learningbased traffic steering. Compared to the Q-learning and heuristic baselines, our results show that the proposed algorithm achieves better performance in terms of 6% and 10% higher average system throughput, and 23% and 33% lower network delay, respectively.


title: Next-to-leading-order solution to Kerr-Newman black hole superradiance

abstract: The superradiant instabilities of Kerr-Newman black holes with charged or uncharged massive spin-0 fields are calculated analytically to the next-to-leading order in the limit of $\alpha\sim r_g \mu \ll 1$. A missing factor of $1/2$ in the previous leading-order result is identified. The next-to-leading order result has a compact form and is in good agreement with existing numerical calculations. The percentage error increases with $\alpha$, from a few percent for $\alpha\sim 0.1$ to about $50\%$ for $\alpha\sim 0.4$. Massive neutral scalars too heavy to be produced with Kerr black hole superradiance may exist in the superradiant region of Kerr-Newman black holes.


title: Language-Informed Transfer Learning for Embodied Household Activities

abstract: For service robots to become general-purpose in everyday household environments, they need not only a large library of primitive skills, but also the ability to quickly learn novel tasks specified by users. Fine-tuning neural networks on a variety of downstream tasks has been successful in many vision and language domains, but research is still limited on transfer learning between diverse long-horizon tasks. We propose that, compared to reinforcement learning for a new household activity from scratch, home robots can benefit from transferring the value and policy networks trained for similar tasks. We evaluate this idea in the BEHAVIOR simulation benchmark which includes a large number of household activities and a set of action primitives. For easy mapping between state spaces of different tasks, we provide a text-based representation and leverage language models to produce a common embedding space. The results show that the selection of similar source activities can be informed by the semantic similarity of state and goal descriptions with the target task. We further analyze the results and discuss ways to overcome the problem of catastrophic forgetting.


title: Nuclear level densities and $\gamma$-ray strength functions in $^{120,124}$Sn isotopes: impact of Porter-Thomas fluctuations

abstract: Nuclear level densities (NLDs) and $\gamma$-ray strength functions (GSFs) of $^{120,124}$Sn have been extracted with the Oslo method from proton-$\gamma$ coincidences in the ($p,p^{\prime}\gamma)$ reaction. The functional forms of the GSFs and NLDs have been further constrained with the Shape method by studying primary $\gamma$-transitions to the ground and first excited states.The NLDs demonstrate good agreement with the NLDs of $^{116,118,122}$Sn isotopes measured previously. Moreover, the extracted partial NLD of 1$^{-}$ levels in $^{124}$Sn is shown to be in fair agreement with those deduced from spectra of relativistic Coulomb excitation in forward-angle inelastic proton scattering. The experimental NLDs have been applied to estimate the magnitude of the Porter-Thomas (PT) fluctuations. Within the PT fluctuations, we conclude that the GSFs for both isotopes can be considered to be independent of initial and final excitation energies, in accordance with the generalized Brink-Axel hypothesis. Particularly large fluctuations observed in the Shape-method GSFs present a considerable contribution to the uncertainty of the method, and may be one of the reasons for deviations from the Oslo-method strength at low $\gamma$-ray energies and low values of the NLD (below $\approx1\cdot10^{3}-2\cdot10^{3}$ MeV$^{-1}$).


title: Modeling Strong Lenses from Wide-Field Ground-Based Observations in KiDS and GAMA

abstract: Despite the success of galaxy-scale strong gravitational lens studies with Hubble-quality imaging, the number of well-studied strong lenses remains small. As a result, robust comparisons of the lens models to theoretical predictions are difficult. This motivates our application of automated Bayesian lens modeling methods to observations from public data releases of overlapping large ground-based imaging and spectroscopic surveys: Kilo-Degree Survey (KiDS) and Galaxy and Mass Assembly (GAMA), respectively. We use the open-source lens modeling software PyAutoLens to perform our analysis. We demonstrate the feasibility of strong lens modeling with large-survey data at lower resolution as a complementary avenue to studies that utilize more time-consuming and expensive observations of individual lenses at higher resolution. We discuss advantages and challenges, with special consideration given to determining background source redshifts from single-aperture spectra and to disentangling foreground lens and background source light. High uncertainties in the best-fit parameters for the models due to the limits of optical resolution in ground-based observatories and the small sample size can be improved with future study. We give broadly applicable recommendations for future efforts, and with proper application this approach could yield measurements in the quantities needed for robust statistical inference.


title: Homeostatic regulation of renewing tissue cell populations via crowding control

abstract: To maintain renewing epithelial tissues in a healthy, homeostatic state, (stem) cell divisions and differentiation need to be tightly regulated. Mechanisms of homeostatic control often rely on crowding control: cells are able to sense the cell density in their environment (via various molecular and mechanosensing pathways) and respond by adjusting division, differentiation, and cell state transitions appropriately. Here we determine, via a mathematically rigorous framework, which general conditions for the crowding feedback regulation (i) must be minimally met, and (ii) are sufficient, to allow the maintenance of homeostasis in renewing tissues. We show that those conditions naturally allow for a degree of robustness toward disruption of regulation. Furthermore, intrinsic to this feedback regulation is that stem cell identity is established collectively by the cell population, not by individual cells, which implies the possibility of `quasi-dedifferentiation', in which cells committed to differentiation may reacquire stem cell properties upon depletion of the stem cell pool. These findings can guide future experimental campaigns to identify specific crowding feedback mechanisms.


title: The collision frequency in two unconventional superconductors

abstract: The collision frequency (also known as the inverse scattering lifetime) can be self-consistently calculated from the imaginary part of the zero-temperature elastic scattering cross-section in unconventional superconductors. We find these types of studies helpful to describe a hidden self-consistent damping due to incoherent fermions in two physical spaces: The Phase Space of the Nonequilibrium Statistical Mechanics, and the Configuration Space of Nonrelativistic Quantum Mechanics. The direct relation of the collision frequency with those well-known Physical Spaces is addressed in a singular way this time. Since the use of collisions for different elastic scattering regimes, is a well-developed formalism using retarded, and advanced Green functions in Many Body Physics; in order to describe our findings, we define and characterize a Reduced Phase Space for collision frequencies in the triplet strontium ruthenate compound and the singlet doped with strontium lanthanum cuprate ceramic. Both compounds display different nodal behavior of the superconducting order parameter. In this work, their zero gap behavior is numerically scanned and used to give some illustrative examples. Finally, it intuitively explores the geometrical nonlocality of the collision frequency of this type of hidden self-consistency in the Boltzmann equation, when the zero superconducting gap value drives the physics below the transition temperature, and incoherent fermions quasiparticles govern several nonequilibrium phenomena since the macroscopic behavior remarkably changes with the strontium atomic potential strength, and the concentration inherent to both compounds.


title: Salient Object Detection for Images Taken by People With Vision Impairments

abstract: Salient object detection is the task of producing a binary mask for an image that deciphers which pixels belong to the foreground object versus background. We introduce a new salient object detection dataset using images taken by people who are visually impaired who were seeking to better understand their surroundings, which we call VizWiz-SalientObject. Compared to seven existing datasets, VizWiz-SalientObject is the largest (i.e., 32,000 human-annotated images) and contains unique characteristics including a higher prevalence of text in the salient objects (i.e., in 68\% of images) and salient objects that occupy a larger ratio of the images (i.e., on average, $\sim$50\% coverage). We benchmarked seven modern salient object detection methods on our dataset and found they struggle most with images featuring salient objects that are large, have less complex boundaries, and lack text as well as for lower quality images. We invite the broader community to work on our new dataset challenge by publicly sharing the dataset at https://vizwiz.org/tasks-and-datasets/salient-object .


title: Surface acoustic wave generation and detection in quantum paraelectric regime of SrTiO$_3$-based heterostructure

abstract: Strontium titanate (STO), apart from being a ubiquitous substrate for complex-oxide heterostructures, possesses a multitude of strongly-coupled electronic and mechanical properties. Surface acoustic wave (SAW) generation and detection offers insight into electromechanical couplings that are sensitive to quantum paraelectricity and other structural phase transitions. Propagating SAWs can interact with STO-based electronic nanostructures, in particular LaAlO$_3$/SrTiO$_3$ (LAO/STO). Here we report generation and detection of SAW within LAO/STO heterointerfaces at cryogenic temperatures ($T\ge$~2 K) using superconducting interdigitated transducers (IDTs). The temperature dependence shows an increase in the SAWs quality factor that saturates at $T\approx 8 $ K. The effect of backgate tuning on the SAW resonance frequency shows the possible acoustic coupling with the ferroelastic domain wall evolution. This method of generating SAWs provides a pathway towards dynamic tuning of ferroelastic domain structures, which are expected to influence electronic properties of complex-oxide nanostructures. Devices which incorporate SAWs may in turn help to elucidate the role of ferroelastic domain structures in mediating electronic behavior.


title: A note on properly discontinuous actions

abstract: We compare various notions of proper discontinuity for group actions. We also discuss fundamental domains and criteria for cocompactness.


title: Kinematics and Origin of Gas in the Disk Galaxy NGC 2655

abstract: The new observational data concerning distribution, excitation, and kinematics of the ionized gas in the giant early-type disk galaxy NGC 2655 obtained at the 6m telescope of the Special Astrophysical Observatory (SAO RAS) and at the 2.5m telescope of the Caucasian Mountain Observatory of the Sternberg Astronomical Institute (CMO SAI MSU) are presented in this work. The joint analysis of these and earlier spectral observations has allowed us to make a conclusion about multiple nature of the gas in NGC 2655. Together with a proper large gaseous disk experiencing regular circular rotation in the equatorial plane of the stellar potential of the galaxy for billions years, we observe also remnants of a merged small satellite having striked the central part of NGC 2655 almost vertically for some 10 million years ago.


title: Blind Judgement: Agent-Based Supreme Court Modelling With GPT

abstract: We present a novel Transformer-based multi-agent system for simulating the judicial rulings of the 2010-2016 Supreme Court of the United States. We train nine separate models with the respective authored opinions of each supreme justice active ca. 2015 and test the resulting system on 96 real-world cases. We find our system predicts the decisions of the real-world Supreme Court with better-than-random accuracy. We further find a correlation between model accuracy with respect to individual justices and their alignment between legal conservatism & liberalism. Our methods and results hold significance for researchers interested in using language models to simulate politically-charged discourse between multiple agents.


title: Tetra-penta-deca-hexagonal-graphene (TPDH-graphene) hydrogenation patterns: dynamics and electronic structure

abstract: The advent of graphene has renewed the interest in other 2D carbon-based materials. Bhattacharya and Jana have proposed a new carbon allotrope, composed of different polygonal carbon rings containing 4, 5, 6, and 10 atoms, named Tetra-Penta-Deca-Hexagonal-graphene (TPDH-graphene). This unusual topology created material with interesting mechanical, electronic, and optical properties and several potential applications, including UV protection. Like other 2D carbon structures, chemical functionalizations can be used to tune their TPDH-graphene properties. In this work, we investigated the hydrogenation dynamics of TPDH-graphene and its effects on its electronic structure, combining DFT and fully atomistic reactive molecular dynamics simulations. Our results show that H atoms are mainly incorporated on tetragonal ring sites (up to 80% at 300 K), leading to the appearance of well-delimited pentagonal carbon stripes. The electronic structure of the hydrogenated structures shows the formation of narrow bandgaps with the presence of Dirac cone-like structures, indicative of anisotropic transport properties.


title: Vanishing of Quartic and Sextic Twists of $L$-functions

abstract: Let $E$ be an elliptic curve over $\mathbf{Q}$. We conjecture asymptotic estimates for the number of vanishings of $L(E,1,\chi)$ as $\chi$ varies over all primitive Dirichlet characters of orders 4 and 6, subject to a mild hypothesis on $E$. Our conjectures about these families come from conjectures about random unitary matrices as predicted by the philosophy of Katz-Sarnak. We support our conjectures with numerical evidence. Earlier work by David, Fearnley and Kisilevsky formulates analogous conjectures for characters of any odd prime order. In the composite order case, however, we need to justify our use of random matrix theory heuristics by analyzing the equidistribution of the squares of normalized Gauss sums. Along the way we introduce the notion of totally order $\ell$ characters to quantify how quickly quartic and sextic Gauss sums become equidistributed. Surprisingly, the rate of equidistribution in the full family of quartic (sextic, resp.) characters is much slower than in the sub-family of totally quartic (sextic, resp.) characters. A conceptual explanation for this phenomenon is that the full family of order $\ell$ twisted elliptic curve $L$-functions, with $\ell$ even and composite, is a mixed family with both unitary and orthogonal aspects.


title: Strain-induced Landau levels of Majorana fermions in an anisotropically interacting Kitaev model on a honeycomb lattice

abstract: The energy structure of an anisotropically interacting Kitaev model on a honeycomb lattice under triaxial strain is investigated. A numerical calculation shows that quantized states appear in the low-energy region, even when the anisotropy of the interaction is rather strong. Their energies are proportional to the square root of the quantum number and the quantized state at zero energy appears only on one sublattice. These findings indicate the emergence of the strain-induced Landau levels of Majorana fermions, which is also confirmed by an analytical calculation. These Landau levels are stable, when the direction of triaxial strain is slightly changed from the bond direction.


title: Detection problems in the spiked matrix models

abstract: We study the statistical decision process of detecting the low-rank signal from various signal-plus-noise type data matrices, known as the spiked random matrix models. We first show that the principal component analysis can be improved by entrywise pre-transforming the data matrix if the noise is non-Gaussian, generalizing the known results for the spiked random matrix models with rank-1 signals. As an intermediate step, we find out sharp phase transition thresholds for the extreme eigenvalues of spiked random matrices, which generalize the Baik-Ben Arous-P\'{e}ch\'{e} (BBP) transition. We also prove the central limit theorem for the linear spectral statistics for the spiked random matrices and propose a hypothesis test based on it, which does not depend on the distribution of the signal or the noise. When the noise is non-Gaussian noise, the test can be improved with an entrywise transformation to the data matrix with additive noise. We also introduce an algorithm that estimates the rank of the signal when it is not known a priori.


title: A Levy-driven Ornstein-Uhlenbeck process for the valuation of credit index swaptions

abstract: A Levy-driven Ornstein-Uhlenbeck process is proposed to model the evolution of the risk-free rate and default intensities for the purpose of evaluating option contracts on a credit index. Time evolution in credit markets is assumed to follow a gamma process in order to reflect the different pace at which credit products are exchanged with respect to that of risk-free debt. Formulas for the characteristic function, zero coupon bonds, moments of the process and its stationary distribution are derived. Numerical experiments showing convergence of standard numerical methods for the valuation PIDE to analytical and Montecarlo solutions are shown. Calibration to market prices of options on a credit index is performed, and model and market implied summary statistics for the underlying credit spreads are estimated and compared.


title: Acceptable Bilateral Gamma Parameters

abstract: The purpose of this paper is to utilize statistical methodologies to infer from market prices of assets and their derivatives the magnitude of the set of a measure M that defines acceptance sets of risky future cash flows. Specifically, we estimate upper and lower boundaries of the compensation needed for a given bilateral gamma distributed future cash flow to be acceptable. We show that prospects theory provides a natural interpretation of the behaviors implied by such boundaries, which are not compatible with expected utility theory over terminal wealth. Boundaries for bilateral gamma risk neutral scale parameters for given speed parameters are also estimated and tested against market data and, in particular, comparisons are made with known empirical facts about the magnitude of the acceptance set of a common class of risk measures.


title: TransfQMix: Transformers for Leveraging the Graph Structure of Multi-Agent Reinforcement Learning Problems

abstract: Coordination is one of the most difficult aspects of multi-agent reinforcement learning (MARL). One reason is that agents normally choose their actions independently of one another. In order to see coordination strategies emerging from the combination of independent policies, the recent research has focused on the use of a centralized function (CF) that learns each agent's contribution to the team reward. However, the structure in which the environment is presented to the agents and to the CF is typically overlooked. We have observed that the features used to describe the coordination problem can be represented as vertex features of a latent graph structure. Here, we present TransfQMix, a new approach that uses transformers to leverage this latent structure and learn better coordination policies. Our transformer agents perform a graph reasoning over the state of the observable entities. Our transformer Q-mixer learns a monotonic mixing-function from a larger graph that includes the internal and external states of the agents. TransfQMix is designed to be entirely transferable, meaning that same parameters can be used to control and train larger or smaller teams of agents. This enables to deploy promising approaches to save training time and derive general policies in MARL, such as transfer learning, zero-shot transfer, and curriculum learning. We report TransfQMix's performances in the Spread and StarCraft II environments. In both settings, it outperforms state-of-the-art Q-Learning models, and it demonstrates effectiveness in solving problems that other methods can not solve.


title: HST Low Resolution Stellar Library

abstract: Hubble Space Telescope's (HST) Space Telescope Imaging Spectrograph (STIS) targeted 556 stars in a long-running program called Next Generation Spectral Library (NGSL) via proposals GO9088, GO9786, GO10222, and GO13776. Exposures through three low resolution gratings provide wavelength coverage from 0.2 $< \lambda <$ 1 $\mu$m at $\lambda/\Delta\lambda\sim$ 1000, providing unique coverage in the ultraviolet (UV). The UV grating (G230LB) scatters red light and this results in unwanted flux that becomes especially troubling for cool stars. We applied scattered light corrections based on \cite{2022stis.rept....5W} and flux corrections arising from pointing errors relative to the center of the 0\farcs2 slit. We present 514 fully reduced spectra, fluxed, dereddened, and cross-correlated to zero velocity. Because of the broad spectral range, we can simultaneously study H$\alpha$ and Mg II $\lambda$2800, indicators of chromospheric activity. Their behaviours are decoupled. Besides three cool dwarfs and one giant with mild flares in H$\alpha$, only Be stars show strong H$\alpha$ emission. Mg2800 emission, however, strongly anti-correlates with temperature such that warm stars show absorption and stars cooler than $5000 \: \! \rm{K}$ universally show chromospheric emission regardless of dwarf/giant status or metallicity. Transformed to Mg2800 flux emerging from the stellar surface, we find a correlation with temperature with approximately symmetric astrophysical scatter, in contrast to other workers who find a basal level with asymmetric scatter to strong values. Unsurprisingly, we confirm that Mg2800 activity is variable.


title: Multitask Weakly Supervised Learning for Origin Destination Travel Time Estimation

abstract: Travel time estimation from GPS trips is of great importance to order duration, ridesharing, taxi dispatching, etc. However, the dense trajectory is not always available due to the limitation of data privacy and acquisition, while the origin destination (OD) type of data, such as NYC taxi data, NYC bike data, and Capital Bikeshare data, is more accessible. To address this issue, this paper starts to estimate the OD trips travel time combined with the road network. Subsequently, a Multitask Weakly Supervised Learning Framework for Travel Time Estimation (MWSL TTE) has been proposed to infer transition probability between roads segments, and the travel time on road segments and intersection simultaneously. Technically, given an OD pair, the transition probability intends to recover the most possible route. And then, the output of travel time is equal to the summation of all segments' and intersections' travel time in this route. A novel route recovery function has been proposed to iteratively maximize the current route's co occurrence probability, and minimize the discrepancy between routes' probability distribution and the inverse distribution of routes' estimation loss. Moreover, the expected log likelihood function based on a weakly supervised framework has been deployed in optimizing the travel time from road segments and intersections concurrently. We conduct experiments on a wide range of real world taxi datasets in Xi'an and Chengdu and demonstrate our method's effectiveness on route recovery and travel time estimation.