title: Marginal Effects for Non-Linear Prediction Functions

abstract: Beta coefficients for linear regression models represent the ideal form of an interpretable feature effect. However, for non-linear models and especially generalized linear models, the estimated coefficients cannot be interpreted as a direct feature effect on the predicted outcome. Hence, marginal effects are typically used as approximations for feature effects, either in the shape of derivatives of the prediction function or forward differences in prediction due to a change in a feature value. While marginal effects are commonly used in many scientific fields, they have not yet been adopted as a model-agnostic interpretation method for machine learning models. This may stem from their inflexibility as a univariate feature effect and their inability to deal with the non-linearities found in black box models. We introduce a new class of marginal effects termed forward marginal effects. We argue to abandon derivatives in favor of better-interpretable forward differences. Furthermore, we generalize marginal effects based on forward differences to multivariate changes in feature values. To account for the non-linearity of prediction functions, we introduce a non-linearity measure for marginal effects. We argue against summarizing feature effects of a non-linear prediction function in a single metric such as the average marginal effect. Instead, we propose to partition the feature space to compute conditional average marginal effects on feature subspaces, which serve as conditional feature effect estimates.


title: Electric and Magnetic Moments and Transition Probabilities in $^{208}$Pb $\pm 1$ Nuclei

abstract: We present moments and transition probabilities in the neighboring odd-mass nuclei of $^{208}$Pb calculated fully self-consistently from the s.p. properties of $^{208}$Pb with polarization corrections from its excitations, both given from previous Skyrme-Hartree-Fock and RPA calculations. The electric results agree nicely with the data with two very interesting exceptions. In the magnetic case we obtain similar results. We discuss also polarization contributions to the $l$-forbidden $M1$ transitions, which are, however, much too small compared to the data. With a modified external field operator which accounts effectively for mesonic and many-body effects the description of the data can be substantially improved.


title: Dynamic Infection Spread Model Based Group Testing

abstract: We study a dynamic infection spread model, inspired by the discrete time SIR model, where infections are spread via non-isolated infected individuals. While infection keeps spreading over time, a limited capacity testing is performed at each time instance as well. In contrast to the classical, static, group testing problem, the objective in our setup is not to find the minimum number of required tests to identify the infection status of every individual in the population, but to control the infection spread by detecting and isolating the infections over time by using the given, limited number of tests. In order to analyze the performance of the proposed algorithms, we focus on the mean-sense analysis of the number of individuals that remain non-infected throughout the process of controlling the infection. We propose two dynamic algorithms that both use given limited number of tests to identify and isolate the infections over time, while the infection spreads. While the first algorithm is a dynamic randomized individual testing algorithm, in the second algorithm we employ the group testing approach similar to the original work of Dorfman. By considering weak versions of our algorithms, we obtain lower bounds for the performance of our algorithms. Finally, we implement our algorithms and run simulations to gather numerical results and compare our algorithms and theoretical approximation results under different sets of system parameters.

categories: stat.ME cs.DS cs.IT eess.SP math.IT q-bio.PE


title: Rapid-then-slow migration reproduces mass distribution of TRAPPIST-1 system

abstract: The TRAPPIST-1 system is an iconic planetary system in various aspects (e.g., habitability, resonant relation, and multiplicity) and hence has attracted considerable attention. The mass distribution of the TRAPPIST-1 planets is characterized by two features: the two inner planets are large, and the masses of the four planets in the outer orbit increase with orbital distance. The origin of these features cannot be explained by previous formation models. We investigate whether the mass distribution of the TRAPPIST-1 system can be reproduced by a planet formation model using N-body simulations. We used a gas disk evolution model around a low-mass star constructed by considering disk winds and followed the growth and orbital migration from planetary embryos with the isolation mass, which increases with orbital distance. As a result, we find that from the initial phase, planets in inner orbits undergo rapid orbital migration, and the coalescence growth near the inner disk edge is enhanced. This allows the inner planets to grow larger. Meanwhile, compared with the inner planets, planets in outer orbits migrate more slowly and do not frequently collide with neighboring planets. Therefore, the trend of increasing mass toward the outer orbit, called reversed mass ranking, is maintained. The final mass distribution approximately agrees with the two features of the mass distribution in the TRAPPIST-1 system. We discover that the mass distribution in the TRAPPIST-1 system can be reproduced when embryos experience rapid migration and become trapped near the disk inner edge, and then more massive embryos undergo slower migration. This migration transition can be achieved naturally in a disk evolution model with disk winds.


title: Cosmological dynamics of multifield dark energy

abstract: We numerically and analytically explore the background cosmological dynamics of multifield dark energy with highly nongeodesic or "spinning" field-space trajectories. These extensions of standard single-field quintessence possess appealing theoretical features and observable differences from the cosmological standard model. At the level of the cosmological background, we perform a phase-space analysis and identify approximate attractors with late-time acceleration for a wide range of initial conditions. Focusing on two classes of field-space geometry, we derive bounds on parameter space by demanding viable late-time acceleration and the absence of gradient instabilities, as well as from the de Sitter swampland conjecture.


title: The large-scale structure from non-Gaussian primordial perturbations

abstract: The late-time effect of primordial non-Gaussianity offers a window into the physics of inflation and the very early Universe. In this work we study the consequences of a particular class of primordial non-Gaussianity that is fully characterized by initial density fluctuations drawn from a non-Gaussian probability density function, rather than by construction of a particular form for the primordial bispectrum. We numerically generate multiple realisations of cosmological structure and use the late-time matter power spectrum, bispectrum and trispectrum to determine the effect of these modified initial conditions. We show that the initial non-Gaussianity has only a small imprint on the first three polyspectra, when compared to a standard Gaussian cosmology. Furthermore, some of our models present an interesting scale-dependent deviation from the Gaussian case in the bispectrum and trispectrum, although the signal is at most at the percent level. The majority of our models are consistent with CMB constraints on the trispectrum, while the others are only marginally excluded. Finally, we discuss further possible extensions of our study.


title: The Dark Ages' 21-cm Trispectrum

abstract: We investigate tomography of 21-cm brightness temperature fluctuations during the Dark Ages as a probe for constraining primordial non-Gaussianity. We expand the 21- cm brightness temperature up to cubic order in perturbation theory and improve previous models of the signal by including the effect of the free electron fraction. Using modified standard perturbation theory methods that include baryonic pressure effects we derive an improved secondary bispectrum and for the first time derive the secondary trispectrum of 21-cm brightness temperature fluctuations. We then forecast the amount of information available from the Dark Ages to constrain primordial non-Gaussianity, including the imprints of massive particle exchange during inflation and we determine how much signal is lost due to secondary non-Gaussianity. We find that although secondary non-Gaussianity swamps the primordial signal, primordial non-Gaussianity can still be extracted with signal-to-noise ratios that surpass current and future CMB experiments by several orders of magnitude, depending on the experimental setup. Furthermore, we conclude that for the bi- and trispectra of massive particle exchange marginalizing over other primordial shapes affects signal-to-noise ratios more severely than secondary shapes. Baryonic pressure effects turn out to have a negligible impact on our forecasts, even at scales close to the Jeans scale. The results of this work reinforce the prospects of 21-cm brightness temperature fluctuations from the Dark Ages as the ultimate probe for primordial non-Gaussianity.


title: Bayes' Theorem, Inflation, and the Ekpyrotic Universe

abstract: We calculate the Bayesian evidences for a class of Ekpyrotic universe models, and compare with a model of single field inflation with a Higgs-type potential. Combining parsimony and observational constraints, this gives us a systematic way to evaluate the degree to which Ekpyrotic models are constrained by CMB data from Planck. We integrate the equations of motion numerically to define a likelihood using Planck 2018 data and sample this likelihood to obtain Bayesian evidences. Priors are justified and used to put Ekpyrotic models and inflation on equal footing. We find reasonable preference for one of the considered Ekpyrotic models over the others, but that even this one is disfavored compared with Higgs inflation.


title: Point-NeRF: Point-based Neural Radiance Fields

abstract: Volumetric neural rendering methods like NeRF generate high-quality view synthesis results but are optimized per-scene leading to prohibitive reconstruction time. On the other hand, deep multi-view stereo methods can quickly reconstruct scene geometry via direct network inference. Point-NeRF combines the advantages of these two approaches by using neural 3D point clouds, with associated neural features, to model a radiance field. Point-NeRF can be rendered efficiently by aggregating neural point features near scene surfaces, in a ray marching-based rendering pipeline. Moreover, Point-NeRF can be initialized via direct inference of a pre-trained deep network to produce a neural point cloud; this point cloud can be finetuned to surpass the visual quality of NeRF with 30X faster training time. Point-NeRF can be combined with other 3D reconstruction methods and handles the errors and outliers in such methods via a novel pruning and growing mechanism. The experiments on the DTU, the NeRF Synthetics , the ScanNet and the Tanks and Temples datasets demonstrate Point-NeRF can surpass the existing methods and achieve the state-of-the-art results.


title: (K)not machine learning

abstract: We review recent efforts to machine learn relations between knot invariants. Because these knot invariants have meaning in physics, we explore aspects of Chern-Simons theory and higher dimensional gauge theories. The goal of this work is to translate numerical experiments with Big Data to new analytic results.


title: Parametric Solution to six n degree powers for degree two, three, four, five, six , seven, eight, and nine

abstract: Historically in math literature there are instances where solutions have been arrived at by different authors for equation 1 for six powers on both side of equation, for different degree 2,3,4,5,6,7,8,9. See reference number 1, by A. Bremner & J. Delorme and reference number 10, by Tito Piezas. The difference is that this paper has done systematic analysis of equation 1. While numerical solutions for equation 1, is available on Wolfram math website, search for parametric solutions to equation 1, in various publications for all degree 2,3,4,5,6,7,8,9 did not yield results. The authors of this paper have selected six terms on each side of equation 1, since the difficulty of the problem increases every time a term is deleted on each side of equation 1. The authors have provided parametric solutions for equation 1, for degree 2, 3, 4, 5 & 6 and for degree 7, 8 & 9 solutions using elliptical curve method has been provided. Also we would like to mention that solutions for degree 7, 8 & 9 using elliptic curve method has infinite numerical solutions. Keywords, Pure math, Diophantine equations, Equal sums, parameteric equations.


title: Lensing Machines: Representing Perspective in Latent Variable Models

abstract: Many datasets represent a combination of different ways of looking at the same data that lead to different generalizations. For example, a corpus with examples generated by different people may be mixtures of many perspectives and can be viewed with different perspectives by others. It isnt always possible to represent the viewpoints by a clean separation, in advance, of examples representing each viewpoint and train a separate model for each viewpoint. We introduce lensing, a mixed initiative technique to extract lenses or mappings between machine learned representations and perspectives of human experts, and to generate lensed models that afford multiple perspectives of the same dataset. We apply lensing for two classes of latent variable models: a mixed membership model, a matrix factorization model in the context of two mental health applications, and we capture and imbue the perspectives of clinical psychologists into these models. Our work shows the benefits of the machine learning practitioner formally incorporating the perspective of a knowledgeable domain expert into their models rather than estimating unlensed models themselves in isolation.


title: Multiple shell ejections on a 100 yr timescale from a massive yellow hypergiant

abstract: This contribution focuses on a rare example of the class of post-Red Supergiants, IRAS 17163-3907, the central star of the Fried Egg nebula. In particular, we discuss some of our recently published results in detail. The inner parts of the circumstellar environment of this evolved massive star are probed at milli-arcsec resolution using VLTI's GRAVITY instrument operating in the K-band (2 $\mu$m), while larger, arcsecond, scales are probed by VISIR diffraction limited images around 10 $\mu$m, supplemented by a complete Spectral Energy Distribution. The spectro-interferometric data cover important diagnostic lines (Br$\gamma$, Na I), which we are able to constrain spatially. Both the presence and size of the Na I doublet in emission has been traditionally challenging to explain towards other objects of this class. In this study we show that a two-zone model in Local Thermal Equilibrium can reproduce both the observed sizes and strengths of the emission lines observed in the K-band, without the need of a pseudo-photosphere. In addition, we find evidence for the presence of a third hot inner shell, and demonstrate that the star has undergone at least three mass-loss episodes over roughly the past century. To explain the properties of the observed non-steady mass-loss we explore pulsation-driven and line-driven mass-loss and introduce the bi-stability jump as a possible underlying mechanism to explain mass-loss towards Yellow Hypergiants.


title: Assessment of the effectiveness of Omicron transmission mitigation strategies for European universities using an agent-based network model

abstract: Returning universities to full on-campus operations while the COVID-19 pandemic is ongoing has been a controversial discussion in many countries. The risk of large outbreaks in dense course settings is contrasted by the benefits of in-person teaching. Transmission risk depends on a range of parameters, such as vaccination coverage and efficacy, number of contacts and adoption of non-pharmaceutical intervention measures (NPIs). Due to the generalised academic freedom in Europe, many universities are asked to autonomously decide on and implement intervention measures and regulate on-campus operations. In the context of rapidly changing vaccination coverage and parameters of the virus, universities often lack sufficient scientific insight to base these decisions on. To address this problem, we analyse a calibrated, data-driven agent-based simulation of transmission dynamics of 10755 students and 974 faculty members in a medium-sized European university. We use a co-location network reconstructed from student enrollment data and calibrate transmission risk based on outbreak size distributions in education institutions. We focus on actionable interventions that are part of the already existing decision-making process of universities to provide guidance for concrete policy decisions. Here we show that, with the Omicron variant of the SARS-CoV-2 virus, even a reduction to 25% occupancy and universal mask mandates are not enough to prevent large outbreaks given the vaccination coverage of about 80% recently reported for students in Austria. Our results show that controlling the spread of the virus with available vaccines in combination with NPIs is not feasible in the university setting if presence of students and faculty on campus is required.


title: Emergent strings at infinite distance with broken supersymmetry

abstract: We investigate the infinite-distance properties of families of unstable flux vacua in string theory with broken supersymmetry. To this end, we employ a generalized notion of distance in the moduli space and we build a holographic description for the non-perturbative regime of the tunneling cascade in terms of a renormalization group flow. In one limit we recover an exponentially light tower of Kaluza-Klein states, while in the opposite limit we find a tower of higher-spin excitations of D1-branes, realizing the emergent string proposal. In particular, the holographic description includes a free sector, whose emergent superconformal symmetry resonates with supersymmetric stability, the CFT distance conjecture and S-duality. We compute the anomalous dimensions of scalar vertex operators and single-trace higher-spin currents, finding an exponential suppression with the distance which is not generic from the renormalization group perspective, but appears specific to our settings.


title: Coupled Hydrodynamics in Dipole-Conserving Quantum Systems

abstract: We investigate the coupled dynamics of charge and energy in interacting lattice models with dipole conservation. We formulate a generic hydrodynamic theory for this combination of fractonic constraints and numerically verify its applicability to the late-time dynamics of a specific bosonic quantum system by developing a microscopic non-equilibrium quantum field theory. Employing a self-consistent $1/N$ approximation in the number of field components, we extract all entries of a generalized diffusion matrix and determine their dependence on microscopic model parameters. We discuss the relation of our results to experiments in ultracold atom quantum simulators.

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


title: Co-evolution of the Brightest Cluster Galaxies and their Host Clusters in IllustrisTNG

abstract: We use the IllustrisTNG simulations to explore the dynamic scaling relation between massive clusters and their central galaxies (BCGs). The Illustris TNG300-1 simulation we use includes 280 massive clusters with $M_{200} > 10^{14}$ M$_{\odot}$ enabling a robust statistical analysis. We derive the line-of-sight velocity dispersion of the stellar particles of the BCGs ($\sigma_{*, BCG}$), analogous to the observed BCG stellar velocity dispersion. We also compute the subhalo velocity dispersion to measure the cluster velocity dispersion ($\sigma_{cl}$). Both $\sigma_{*, BCG}$ and $\sigma_{cl}$ are proportional to the cluster halo mass, but the slopes differ slightly. Thus like the observed relation, $\sigma_{*, BCG} / \sigma_{cl}$ declines as a function of $\sigma_{cl}$, but the scatter is large. We explore the redshift evolution of $\sigma_{*, BCG} - \sigma_{cl}$ scaling relation for $z \lesssim 1$ in a way that can be compared directly with observations. The scaling relation has a similar slope at high redshift, but the scatter increases because of the large scatter in $\sigma_{*, BCG}$. The simulations imply that high redshift BCGs are dynamically more complex than their low redshift counterparts.


title: Dynamical friction of black holes in ultralight dark matter

abstract: In this work we derive simple closed-form expressions for the dynamical friction acting on black holes moving through ultralight (scalar field) dark matter, covering both non-relativistic and relativistic black hole speeds. Our derivation is based on long known scattering amplitudes in black hole spacetimes, it includes the effect of black hole spin and can be easily extended to vector and tensor light fields. Our results cover and complement recent numerical and previous non-relativistic treatments of dynamical friction in ultralight dark matter.


title: Extremely massive disc galaxies in the nearby Universe form through gas-rich minor mergers

abstract: In our hierarchical structure-formation paradigm, the observed morphological evolution of massive galaxies -- from rotationally-supported discs to dispersion-dominated spheroids -- is largely explained via galaxy merging. However, since mergers are likely to destroy discs, and the most massive galaxies have the richest merger histories, it is surprising that any discs exist at all at the highest stellar masses. Recent theoretical work by our group has used a cosmological, hydrodynamical simulation to suggest that extremely massive (M* > 10^11.4 MSun) discs form primarily via minor mergers between spheroids and gas-rich satellites, which create new rotational stellar components and leave discs as remnants. Here, we use UV-optical and HI data of massive galaxies, from the SDSS, GALEX, DECaLS and ALFALFA surveys, to test these theoretical predictions. Observed massive discs account for ~13% of massive galaxies, in good agreement with theory (~11%). ~64% of the observed massive discs exhibit tidal features, which are likely to indicate recent minor mergers, in the deep DECaLS images (compared to ~60% in their simulated counterparts). The incidence of these features is at least four times higher than in low-mass discs, suggesting that, as predicted, minor mergers play a significant (and outsized) role in the formation of these systems. The empirical star-formation rates agree well with theoretical predictions and, for a small galaxy sample with HI detections, the HI masses and fractions are consistent with the range predicted by the simulation. The good agreement between theory and observations indicates that extremely massive discs are indeed remnants of recent minor mergers between spheroids and gas-rich satellites.


title: Exact wavefunction dualities and phase diagrams of 3D quantum vertex models

abstract: We construct a broad class of frustration-free quantum vertex models in 3+1D whose ground states are weighted superpositions of classical 3D vertex model configurations. For the simplest members of this class, the corresponding classical vertex models have a $\mathbb{Z}_2$ gauge constraint enriched with a $\mathbb{Z}_2$ global symmetry. We deduce a phase diagram for the quantum vertex models using exact dualities and mappings to classical spin models and height models. We find $\mathbb{Z}_2$ deconfined and $\mathbb{Z}_2$ symmetry-broken phases which are mapped to one another by duality. The transition between these phases is governed by a self-dual gapless point which is related by a different duality to the Rokhsar-Kivelson (RK) point of quantum spin ice. Our results are illustrated for diamond and cubic lattices, but hold for general 3D lattices with even coordination number.


title: Bona-Masso slicing conditions and the lapse close to black-hole punctures

abstract: We consider several families of functions $f(\alpha)$ that appear in the Bona-Masso slicing condition for the lapse function $\alpha$. Focusing on spherically symmetric and time-independent slices we apply these conditions to the Schwarzschild spacetime in order to construct analytical expressions for the lapse $\alpha$ in terms of the areal radius $R$. We then transform to isotropic coordinates and determine the dependence of $\alpha$ on the isotropic radius $r$ in the vicinity of the black-hole puncture. We propose generalizations of previously considered functions $f(\alpha)$ for which, to leading order, the lapse is proportional to $r$ rather than a non-integer power of $r$. We also perform dynamical simulations in spherical symmetry and demonstrate advantages of the above choices in numerical simulations employing spectral methods.


title: The size-luminosity relation of lensed galaxies at $z=6-9$ in the Hubble Frontier Fields

abstract: We measure the size-luminosity relation of photometrically-selected galaxies within the redshift range $z\sim6-9$, using galaxies lensed by six foreground Hubble Frontier Fields (HFF) clusters. The power afforded by strong gravitational lensing allows us to observe fainter and smaller galaxies than in blank fields. We select our sample of galaxies and obtain their properties, e.g., redshift, magnitude, from the photometrically-derived ASTRODEEP catalogues. The intrinsic size is measured with the Lenstruction software, and completeness maps are created as a function of size and luminosity via the GLACiAR2 software. We perform a Bayesian analysis to estimate the intrinsic and incompleteness-corrected size-luminosity distribution, with parameterization $r_e \propto L^\beta$. We find slopes of $\beta\sim0.48\pm0.08$ at $z\sim6-7$ and $\beta\sim0.68\pm0.14$ at $z\sim8.5$, adopting the Bradac lens model. The slope derived by lensed galaxies is steeper than that obtained in blank fields and is consistent with other independent determinations of the size-luminosity relation from the HFF dataset. We also investigate the systematic uncertainties correlated with the choice of lens models, finding that the slopes of size-luminosity relations derived from different models are consistent with each other, i.e. the modeling errors are not a significant source of uncertainty in the size-luminosity relation.


title: A universe field theory for JT gravity

abstract: We present a field theory description for the non-perturbative splitting and joining of baby universes in Euclidean Jackiw-Teitelboim (JT) gravity. We show how the gravitational path integral, defined as a sum over topologies, can be reproduced from the perturbative expansion of a Kodaira-Spencer (KS) field theory for the complex structure deformations of the spectral curve. We use that the Schwinger-Dyson equations for the KS theory can be mapped to the topological recursion relations. We refer to this dual description of JT gravity as a `universe field theory'. By introducing non-compact D-branes in the target space geometry, we can probe non-perturbative aspects of JT gravity. The relevant operators are obtained through a modification of the JT path integral with Neumann boundary conditions. The KS/JT identification suggests that the ensemble average for JT gravity can be understood in terms of a more standard open/closed duality in topological string theory.


title: GreaseLM: Graph REASoning Enhanced Language Models for Question Answering

abstract: Answering complex questions about textual narratives requires reasoning over both stated context and the world knowledge that underlies it. However, pretrained language models (LM), the foundation of most modern QA systems, do not robustly represent latent relationships between concepts, which is necessary for reasoning. While knowledge graphs (KG) are often used to augment LMs with structured representations of world knowledge, it remains an open question how to effectively fuse and reason over the KG representations and the language context, which provides situational constraints and nuances. In this work, we propose GreaseLM, a new model that fuses encoded representations from pretrained LMs and graph neural networks over multiple layers of modality interaction operations. Information from both modalities propagates to the other, allowing language context representations to be grounded by structured world knowledge, and allowing linguistic nuances (e.g., negation, hedging) in the context to inform the graph representations of knowledge. Our results on three benchmarks in the commonsense reasoning (i.e., CommonsenseQA, OpenbookQA) and medical question answering (i.e., MedQA-USMLE) domains demonstrate that GreaseLM can more reliably answer questions that require reasoning over both situational constraints and structured knowledge, even outperforming models 8x larger.


title: Multicore Quantum Computing

abstract: Any architecture for practical quantum computing must be scalable. An attractive approach is to create multiple cores, computing regions of fixed size that are well-spaced but interlinked with communication channels. This exploded architecture can relax the demands associated with a single monolithic device: the complexity of control, cooling and power infrastructure as well as the difficulties of cross-talk suppression and near-perfect component yield. Here we explore interlinked multicore architectures through analytic and numerical modelling. While elements of our analysis are relevant to diverse platforms, our focus is on semiconductor electron spin systems in which numerous cores may exist on a single chip. We model shuttling and microwave-based interlinks and estimate the achievable fidelities, finding values that are encouraging but markedly inferior to intra-core operations. We therefore introduce optimsed entanglement purification to enable high-fidelity communication, finding that $99.5\%$ is a very realistic goal. We then assess the prospects for quantum advantage using such devices in the NISQ-era and beyond: we simulate recently proposed exponentially-powerful error mitigation schemes in the multicore environment and conclude that these techniques impressively suppress imperfections in both the inter- and intra-core operations.


title: Stochastic normalizing flows as non-equilibrium transformations

abstract: Normalizing flows are a class of deep generative models that provide a promising route to sample lattice field theories more efficiently than conventional Monte Carlo simulations. In this work we show that the theoretical framework of stochastic normalizing flows, in which neural-network layers are combined with Monte Carlo updates, is the same that underlies out-of-equilibrium simulations based on Jarzynski's equality, which have been recently deployed to compute free-energy differences in lattice gauge theories. We lay out a strategy to optimize the efficiency of this extended class of generative models and present examples of applications.

categories: hep-lat cond-mat.stat-mech cs.LG stat.ML


title: The central engine of the highest redshift blazar

abstract: We report on a LUCI/Large Binocular Telescope near-infrared (NIR) spectrum of PSO J030947.49+271757.31 (hereafter PSO J0309+27), the highest redshift blazar known to date (z$\sim$6.1). From the C$\rm IV$$\lambda$1549 broad emission line we found that PSO J0309+27 is powered by a 1.45$^{+1.89}_{-0.85}$$\times$10$^9$M$_{\odot}$ supermassive black hole (SMBH) with a bolometric luminosity of $\sim$8$\times$10$^{46}$ erg s$^{-1}$ and an Eddington ratio equal to 0.44$^{+0.78}_{-0.35}$. We also obtained new photometric observations with the Telescopio Nazionale Galileo in J and K bands to better constrain the NIR Spectral Energy Distribution of the source. Thanks to these observations, we were able to model the accretion disk and to derive an independent estimate of the black hole mass of PSO J0309+27, confirming the value inferred from the virial technique. The existence of such a massive SMBH just $\sim$900 million years after the Big Bang challenges models of the earliest SMBH growth, especially if jetted Active Galactic Nuclei are associated to a highly spinning black hole as currently thought. Indeed, in a Eddington-limited accretion scenario and assuming a radiative efficiency of 0.3, typical of a fast rotating SMBH, a seed black hole of more than 10$^6$ M$_{\odot}$ at z = 30 is required to reproduce the mass of PSO J0309+27 at redshift 6. This requirement suggests either earlier periods of rapid black hole growth with super-Eddington accretion and/or that only part of the released gravitational energy goes to heat the accretion disk and feed the black hole.


title: Rapidity dependence of initial state geometry and momentum correlations in p+Pb collisions

abstract: Event geometry and initial state correlations have been invoked as possible explanations of long range azimuthal correlations observed in high multiplicity p+p and p+Pb collisions. We study the rapidity dependence of initial state momentum correlations and event-by-event geometry in $\sqrt{s}=5.02~\rm{TeV}$ p+Pb collisions within the 3+1D IP-Glasma model~\cite{Schenke:2016ksl}, where the longitudinal structure is governed by JIMWLK rapidity evolution of the incoming nuclear gluon distributions. We find that the event geometry is correlated across large rapidity intervals whereas initial state momentum correlations are relatively short range in rapidity. Based on our results, we discuss implications for the relevance of both effects in explaining the origin of collective phenomena in small systems.


title: On the in vivo recognition of kidney stones using machine learning

abstract: Determining the type of kidney stones allows urologists to prescribe a treatment to avoid recurrence of renal lithiasis. An automated in-vivo image-based classification method would be an important step towards an immediate identification of the kidney stone type required as a first phase of the diagnosis. In the literature it was shown on ex-vivo data (i.e., in very controlled scene and image acquisition conditions) that an automated kidney stone classification is indeed feasible. This pilot study compares the kidney stone recognition performances of six shallow machine learning methods and three deep-learning architectures which were tested with in-vivo images of the four most frequent urinary calculi types acquired with an endoscope during standard ureteroscopies. This contribution details the database construction and the design of the tested kidney stones classifiers. Even if the best results were obtained by the Inception v3 architecture (weighted precision, recall and F1-score of 0.97, 0.98 and 0.97, respectively), it is also shown that choosing an appropriate colour space and texture features allows a shallow machine learning method to approach closely the performances of the most promising deep-learning methods (the XGBoost classifier led to weighted precision, recall and F1-score values of 0.96). This paper is the first one that explores the most discriminant features to be extracted from images acquired during ureteroscopies.


title: Vehicle Routing with Stochastic Demands and Partial Reoptimization

abstract: We consider the vehicle routing problem with stochastic demands (VRPSD), a problem in which customer demands are known in distribution at the route planning stage and revealed during route execution upon arrival at each customer. A long-standing open question on the VRPSD concerns the benefits of allowing, during route execution, partial reordering of the planned customer visits. Given the practical importance of this question and the growing interest on the VRPSD under optimal restocking, we study the VRPSD under a recourse policy known as the switch policy. The switch policy is a canonical reoptimization policy that permits only pairs of successive customers to be reordered. We consider this policy jointly with optimal preventive restocking and introduce a branch-cut-and-price algorithm to compute optimal a priori routing plans. This algorithm features pricing routines where value functions represent the expected cost-to-go along planned routes for all possible states and reordering decisions. To ensure pricing tractability, we adopt a strategy that combines elementary pricing with completion bounds of varying complexity, and solve the pricing problem without relying on dominance rules. Our numerical experiments demonstrate the effectiveness of the algorithm for solving instances with up to 50 customers. Notably, they also give us new insights into the value of reoptimization. The switch policy enables significant cost savings over optimal restocking when the planned routes come from an algorithm built on a deterministic approximation of the data, an important scenario given the difficulty of finding optimal VRPSD solutions. The benefits are smaller when comparing optimal a priori VRPSD solutions obtained for both recourse policies. As it appears, further cost savings may require joint reordering and reassignment of customer visits among vehicles when the context permits.


title: LBT search for companions and sub-structures in the (pre)transitional disk of AB Aurigae

abstract: Multi-wavelengths high-resolution imaging of protoplanetary disks has revealed the presence of multiple, varied substructures in their dust and gas components which might be signposts of young, forming planetary systems. AB Aurigae bears an emblematic (pre)transitional disk showing spiral structures observed in the inner cavity of the disk in both the sub-millimeter (ALMA; 1.3mm, $^{12}$CO) and near-infrared (SPHERE; 1.5-2.5$\mu$m) wavelengths which have been claimed to arise from dynamical interactions with a massive companion. In this work, we present new deep $K_s$ (2.16$\mu$m) and L' (3.7$\mu$m) band images of AB Aurigae obtained with LMIRCam on the Large Binocular Telescope, aimed for the detection of both planetary companions and extended disk structures. No point source is recovered, in particular at the outer regions of the disk, where a putative candidate ($\rho = 0.681", PA = 7.6^{\circ}$) had been previously claimed. The nature of a second innermost planet candidate ($\rho = 0.16'', PA = 203.9^{\circ}$) can not be investigated by the new data. We are able to derive 5$\sigma$ detection limits in both magnitude and mass for the system, going from 14 \Mjup at 0.3'' (49 au) down to 3-4 \Mjup at 0.6'' (98 au) and beyond, based on the ATMO 2020 evolutionary models. We detect the inner spiral structures (< 0.5'') resolved in both CO and polarimetric H-band observations. We also recover the ring structure of the system at larger separation (0.5-0.7") showing a clear south-east/north-west asymmetry. This structure, observed for the first time at L'-band, remains interior to the dust cavity seen at ALMA, suggesting an efficient dust trapping mechanism at play in the disk.


title: Detection of Near-Infrared Water Ice at the Surface of the (pre)Transitional Disk of AB Aur: Informing Icy Grain Abundance, Composition, and Size

abstract: We present near-infrared Large Binocular Telescope Interferometer LMIRCam imagery of the disk around the Herbig Ae/Be star AB Aurigae. A comparison of surface brightness at Ks (2.16 ${\mu}$m), H2O narrowband (3.08 ${\mu}$m), and L' (3.7 ${\mu}$m) allows us to probe the presence of icy grains in this (pre)transitional disk environment. By applying Reference Differential Imaging PSF subtraction, we detect the disk at high signal to noise in all three bands. We find strong morphological differences between bands, including asymmetries consistent with observed spiral arms within 100 AU in L'. An apparent deficit of scattered light at 3.08 ${\mu}$m relative to bracketing wavelengths (Ks and L') is evocative of ice absorption at the disk surface layer. However, the ${\Delta}$(Ks-H2O) color is consistent with grains with little to no ice (0-5% by mass). The ${\Delta}$(H2O-L') color, conversely, suggests grains with a much higher ice mass fraction (~0.68), and the two colors cannot be reconciled under a single grain population model. Additionally, we find the extremely red ${\Delta}$(Ks-L') disk color cannot be reproduced under conventional scattered light modeling with any combination of grain parameters or reasonable local extinction values. We hypothesize that the scattering surfaces at the three wavelengths are not co-located, and optical depth effects result in each wavelength probing the grain population at different disk surface depths. The morphological similarity between Ks and H2O suggests their scattering surfaces are near one another, lending credence to the ${\Delta}$(Ks-H2O) disk color constraint of < 5% ice mass fraction for the outermost scattering disk layer.


title: Linear series with $\rho < 0$ via thrifty lego-building

abstract: The moduli space $\mathcal{G}^r_{g,d} \to \mathcal{M}_g$ parameterizing algebraic curves with a linear series of degree $d$ and rank $r$ has expected relative dimension $\rho = g - (r+1)(g-d+r)$. Classical Brill-Noether theory concerns the case $\rho \geq 0$; we consider the non-surjective case $\rho < 0$. We prove the existence of components of this moduli space with the expected relative dimension when $0 > \rho \geq -g+3$, or $0 > \rho \geq -C_r g + \mathcal{O}(g^{5/6})$, where $C_r$ is a constant depending on the rank of the linear series such that $C_r \to 3$ as $r \to \infty$. These results are proved via a two-marked-point generalization suitable for inductive arguments, and the regeneration theorem for limit linear series.


title: Sum Expressions for Kubota-Leopoldt $p$-adic $L$-functions

abstract: When $p$ is an odd prime, Delbourgo observed that any Kubota-Leopoldt $p$-adic $L$-function, when multiplied by an auxiliary Euler factor, can be written as an infinite sum. We shall establish such expressions without restriction on $p$, and without the Euler factor when the character is nontrivial, by computing the periods of appropriate measures. As an application, we will reprove the Ferrero-Greenberg formula for the derivative $L_p'(0,\chi)$. We will also discuss the convergence of sum expressions in terms of elementary $p$-adic analysis, as well as their relation to Stickelberger elements; such discussions in turn give alternative proofs of the validity of sum expressions.


title: The ALICE Experiment Upgrades

abstract: The ALICE experiment profited of the Long Shutdown during 2019-2021 in order to expand its physics capabilities and fully profit from the increased LHC luminosity in Run 3. The Inner Tracking System has been replaced with a new silicon tracker based on MAPS technology, and a new tracking device has been added in front of the Muon Spectrometer to improve its vertexing capabilities. The wire chambers for TPC readout have been replaced with new GEM detectors which will minimize ion backflow and allow for continuous data taking: moreover, a new detector array dedicated to fast triggering has been installed. On the software side, a new first pass reconstruction was added in order to handle and reduce the data flow and storage. These upgrades will be presented together with an outlook of the future ALICE upgrades in view of the LHC Run 4, which will include the replacement of the ITS inner tracking layers with upgraded silicon devices and a high-granularity electromagnetic and hadronic calorimeter in the forward direction (FOCAL)


title: Observation of spin-1 tunneling on a quantum computer

abstract: Spin-1 tunneling and splitting of energy levels as a result of tunneling are observed explicitly on IBM's quantum computer, ibmq-bogota. The spin-1 is realized with two spins-1/2. We detect oscillations of spin-1 between the states $|1\rangle$, $|-1\rangle$ in the result of tunneling on the basis of studies of the time dependence of the mean value of z-component of spin-1 on the quantum device. The energy level splitting is observed quantifying on the IBM's quantum computer the eigenvalues of Hamiltonian which describes the spin tunneling.


title: Computation of weight 1 modular forms with exotic representations

abstract: We present a deterministic algorithm for computing spaces of weight 1 modular forms with exotic representations. This algorithm is an improved version of Schaeffer's Hecke stability method, utilising the author's previous work on the twist-minimal trace formula for weight 2 holomorphic forms, and presenting a method of lifting forms from characteristic p. The algorithm was used to compute all such forms with level at most 10,000. Together with Sutherland's computation of dihedral forms, this allows us to present the dimensions of all weight 1 newform spaces up to level 10,000.


title: Note on $p$-adic Local Functional Equation

abstract: Given primes $\ell\ne p$, we record here a $p$-adic valued Fourier theory on a local field over $\mathbf{Q}_\ell$, which is developed under the perspective of group schemes. As an application, by substituting rigid analysis for complex analysis, it leads naturally to the $p$-adic local functional equation at $\ell$, which strongly resembles the complex one in Tate's thesis.


title: Information-Based Trading

abstract: We consider a pair of traders in a market where the information available to the second trader is a strict subset of the information available to the first trader. The traders make prices based on the information available concerning a security that pays a random cash flow at a fixed time $T$ in the future. Market information is modelled in line with the scheme of Brody, Hughston & Macrina (2007, 2008) and Brody, Davis, Friedman & Hughston (2009). The risk-neutral distribution of the cash flow is known to the traders, who make prices with a fixed multiplicative bid-offer spread and report their prices to a game master who declares that a trade has been made when the bid price of one of the traders crosses the offer price of the other. We prove that the value of the first trader's position is strictly greater than that of the second. The results are analyzed by use of simulation studies and generalized to situations where (a) there is a hierarchy of traders, (b) there are multiple successive trades, and (c) there is inventory aversion.


title: Emittance preservation in advanced accelerators

abstract: Emittance is a beam quality that is vital for many future applications of advanced accelerators, such as compact free-electron lasers and linear colliders. In this paper, we review the challenges of preserving the transverse emittance during acceleration, both inside and outside accelerator stages. Sources of emittance growth range from space charge and instabilities caused by transverse wakefields, which can occur in any advanced accelerator scheme regardless of medium or driver type, to sources more specific to plasma accelerators, such as mismatching, misalignment, ion motion, Coulomb scattering, chromaticity between stages, and more.


title: Variational Autoencoder based Metamodeling for Multi-Objective Topology Optimization of Electrical Machines

abstract: Conventional magneto-static finite element analysis of electrical machine design is time-consuming and computationally expensive. Since each machine topology has a distinct set of parameters, design optimization is commonly performed independently. This paper presents a novel method for predicting Key Performance Indicators (KPIs) of differently parameterized electrical machine topologies at the same time by mapping a high dimensional integrated design parameters in a lower dimensional latent space using a variational autoencoder. After training, via a latent space, the decoder and multi-layer neural network will function as meta-models for sampling new designs and predicting associated KPIs, respectively. This enables parameter-based concurrent multi-topology optimization.


title: Classical Simulation of Variational Quantum Classifiers using Tensor Rings

abstract: In recent times, Variational Quantum Circuits (VQC) have been widely adopted to different tasks in machine learning such as Combinatorial Optimization and Supervised Learning. With the growing interest, it is pertinent to study the boundaries of the classical simulation of VQCs to effectively benchmark the algorithms. Classically simulating VQCs can also provide the quantum algorithms with a better initialization reducing the amount of quantum resources needed to train the algorithm. This manuscript proposes an algorithm that compresses the quantum state within a circuit using a tensor ring representation which allows for the implementation of VQC based algorithms on a classical simulator at a fraction of the usual storage and computational complexity. Using the tensor ring approximation of the input quantum state, we propose a method that applies the parametrized unitary operations while retaining the low-rank structure of the tensor ring corresponding to the transformed quantum state, providing an exponential improvement of storage and computational time in the number of qubits and layers. This approximation is used to implement the tensor ring VQC for the task of supervised learning on Iris and MNIST datasets to demonstrate the comparable performance as that of the implementations from classical simulator using Matrix Product States.


title: A Grid Fault Tolerant Doubly Fed Induction Generator Wind Turbine via Series Connected Grid Side Converter

abstract: With steadily increasing wind turbine penetration, regulatory standards for grid interconnection are evolving to require that wind generation systems ride-through disturbances such as faults and support the grid during such events. Conventional modifications to the doubly fed induction generation (DFIG) architecture for providing ride-through result in limited control of the turbine shaft and grid current during fault events. A DFIG architecture in which the grid side converter is connected in series as opposed to parallel with the grid connection has shown improved low voltage ride through but poor power processing capabilities. In this paper a unified DFIG wind turbine architecture which employs a parallel grid side rectifier and series grid side converter is presented. The combination of these two converters enables unencumbered power processing and improved voltage disturbance ride through. A dynamic model and control structure for this unified architecture is developed. The operation of the system is illustrated using computer simulations.


title: Development and characterization of a fast and low noise readout for the next generation X-ray CCDs

abstract: The broad energy response, low electronic read noise, and good energy resolution have made X-ray Charge-Coupled Devices (CCDs) an obvious choice for developing soft X-ray astronomical instruments over the last half century. They also come in large array formats with small pixel sizes which make them a potential candidate for the next generation astronomical X-ray missions. However, the next generation X-ray telescopic experiments propose for significantly larger collecting area compared to the existing observatories in order to explore the low luminosity and high redshift X-ray universe which requires these detectors to have an order of magnitude faster readout. In this context, the Stanford University (SU) in collaboration with the Massachusetts Institute of Technology (MIT) has initiated the development of fast readout electronics for X-ray CCDs. At SU, we have designed and developed a fast and low noise readout module with the goal of achieving a readout speed of 5 Mpixel/s. We successfully ran a prototype CCD matrix of 512 $\times$ 512 pixels at 4 Mpixels/s. In this paper, we describe the details of the readout electronics and report the performance of the detectors at these readout speeds in terms of read noise and energy resolution. In the future, we plan to continue to improve performance of the readout module and eventually converge to a dedicated ASIC based readout system to enable parallel read out of large array multi-node CCD devices.


title: Probing the Optical Dynamics of Quantum Emitters in Hexagonal Boron Nitride

abstract: Hexagonal boron nitride is a van der Waals material that hosts visible-wavelength quantum emitters at room temperature. However, experimental identification of the quantum emitters' electronic structure is lacking, and key details of their charge and spin properties remain unknown. Here, we probe the optical dynamics of quantum emitters in hexagonal boron nitride using photon emission correlation spectroscopy. Several quantum emitters exhibit ideal single-photon emission with noise-limited photon antibunching, $g^{(2)}(0)=0$. The photoluminescence emission lineshapes are consistent with individual vibronic transitions. However, polarization-resolved excitation and emission suggests the role of multiple optical transitions, and photon emission correlation spectroscopy reveals complicated optical dynamics associated with excitation and relaxation through multiple electronic excited states. We compare the experimental results to quantitative optical dynamics simulations, develop electronic structure models that are consistent with the observations, and discuss the results in the context of ab initio theoretical calculations.


title: Further Evidence for the Minifilament-Eruption Scenario for Solar Polar Coronal Jets

abstract: We examine a sampling of 23 polar-coronal-hole jets. We first identified the jets in soft X-ray (SXR) images from the X-ray telescope (XRT) on the Hinode spacecraft, over 2014-2016. During this period, frequently the polar holes were small or largely obscured by foreground coronal haze, often making jets difficult to see. We selected 23 jets among those adequately visible during this period, and examined them further using Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA) 171, 193, 211, and 304 Ang images. In SXRs we track the lateral drift of the jet spire relative to the jet base's jet bright point (JBP). In 22 of 23 jets the spire either moves away from (18 cases) or is stationary relative to (4 cases) the JBP. The one exception where the spire moved toward the JBP may be a consequence of line-of-sight projection effects at the limb. From the AIA images, we clearly identify an erupting minifilament in 20 of the 23 jets, while the remainder are consistent with such an eruption having taken place. We also confirm that some jets can trigger onset of nearby "sympathetic" jets, likely because eruption of the minifilament field of the first jet removes magnetic constraints on the base-field region of the second jet. The propensity for spire drift away from the JBP, the identification of the erupting minifilament in the majority of jets, and the magnetic-field topological changes that lead to sympathetic jets, all support or are consistent with the minifilament-eruption model for jets.


title: The Rational Selection of Goal Operations and the Integration ofSearch Strategies with Goal-Driven Autonomy

abstract: Intelligent physical systems as embodied cognitive systems must perform high-level reasoning while concurrently managing an underlying control architecture. The link between cognition and control must manage the problem of converting continuous values from the real world to symbolic representations (and back). To generate effective behaviors, reasoning must include a capacity to replan, acquire and update new information, detect and respond to anomalies, and perform various operations on system goals. But, these processes are not independent and need further exploration. This paper examines an agent's choices when multiple goal operations co-occur and interact, and it establishes a method of choosing between them. We demonstrate the benefits and discuss the trade offs involved with this and show positive results in a dynamic marine search task.


title: Triple lines on a cubic threefold

abstract: The present paper deals with lines contained in a smooth complex cubic threefold. It is well-known that the set of lines of the second type on a cubic threefold is a curve on its Fano surface. Here we give a description of the singularities of this curve.


title: Galois scaffolds for cyclic $p^n$-extensions in characteristic $p$

abstract: Let $K$ be a local field of characteristic $p$ and let $L/K$ be a totally ramified Galois extension such that Gal$(L/K)\cong C_{p^n}$. In this paper we find sufficient conditions for $L/K$ to admit a Galois scaffold. This leads to sufficient conditions for the ring of integers $O_L$ to be free of rank 1 over its associated order $A_0$, and to stricter conditions which imply that $A_0$ is a Hopf order in the group ring $K[C_{p^n}]$.


title: Improved unitary uncertainty relations

abstract: We derive strong variance-based uncertainty relations for arbitrary two and more unitary operators by re-examining the mathematical foundation of the uncertainty relation. This is achieved by strengthening the celebrated Cauchy-Schwarz inequality using a method of brackets and convex functions. The unitary uncertainty relations outperform several strong unitary uncertainty relations, notably better than some recent best lower bounds such as [Phys. Rev. Lett. 120, 230402 (2018)] and [Phys. Rev. A. 100, 022116 (2019)].


title: Image-to-Video Re-Identification via Mutual Discriminative Knowledge Transfer

abstract: The gap in representations between image and video makes Image-to-Video Re-identification (I2V Re-ID) challenging, and recent works formulate this problem as a knowledge distillation (KD) process. In this paper, we propose a mutual discriminative knowledge distillation framework to transfer a video-based richer representation to an image based representation more effectively. Specifically, we propose the triplet contrast loss (TCL), a novel loss designed for KD. During the KD process, the TCL loss transfers the local structure, exploits the higher order information, and mitigates the misalignment of the heterogeneous output of teacher and student networks. Compared with other losses for KD, the proposed TCL loss selectively transfers the local discriminative features from teacher to student, making it effective in the ReID. Besides the TCL loss, we adopt mutual learning to regularize both the teacher and student networks training. Extensive experiments demonstrate the effectiveness of our method on the MARS, DukeMTMC-VideoReID and VeRi-776 benchmarks.


title: The Edge-on Galaxies in the Pan-STARRS survey (EGIPS)

abstract: We present a catalogue of 16551 edge-on galaxies created using the public DR2 data of the Pan-STARRS survey. The catalogue covers the three quarters of the sky above Dec.=-30 degrees. The galaxies were selected using a convolutional neural network, trained on a sample of edge-on galaxies identified earlier in the SDSS survey. This approach allows us to dramatically improve the quality of the candidate selection and perform a thorough visual inspection in a reasonable amount of time. The catalogue provides homogeneous information on astrometry, SExtractor photometry, and non-parametric morphological statistics of the galaxies. The photometry is reliably for objects in the 13.8-17.4 r-band magnitude range. According to the HyperLeda database, redshifts are known for about 63 percent of the galaxies in the catalogue. Our sample is well separated into the red sequence and blue cloud galaxy populations. The edge-on galaxies of the red sequence are systematically Delta(g-i)~0.1 mag redder than galaxies oriented at an arbitrary angle to the observer. We found a variation of the galaxy thickness with the galaxy colour. The red sequence galaxies are thicker than the galaxies of the blue cloud. In the blue cloud, on average, thinner galaxies turn out to be bluer. In the future, based on this catalogue it is intended to explore the three-dimensional structure of galaxies of different morphologies, as well as to study the scaling relations for discs and bulges.


title: Automated Catheter Tip Repositioning for Intra-cardiac Echocardiography

abstract: Purpose: Intra-Cardiac Echocardiography (ICE) is a powerful imaging modality for guiding cardiac electrophysiology and structural heart interventions. ICE provides real-time observation of anatomy and devices, while enabling direct monitoring of potential complications. In single operator settings, the physician needs to switch back-and-forth between the ICE catheter and therapy device, making continuous ICE support impossible. Two operators setup are therefore sometimes implemented, with the challenge of increase room occupation and cost. Two operator setups are sometimes implemented, but increase procedural costs and room occupation. Methods: ICE catheter robotic control system is developed with automated catheter tip repositioning (i.e. view recovery) method, which can reproduce important views previously navigated to and saved by the user. The performance of the proposed method is demonstrated and evaluated in a combination of heart phantom and animal experiments. Results: Automated ICE view recovery achieved catheter tip position accuracy of 2.09 +/-0.90 mm and catheter image orientation accuracy of 3.93 +/- 2.07 degree in animal studies, and 0.67 +/- 0.79 mm and 0.37 +/- 0.19 degree in heart phantom studies, respectively. Our proposed method is also successfully used during transeptal puncture in animals without complications, showing the possibility for fluoro-less transeptal puncture with ICE catheter robot. Conclusion: Robotic ICE imaging has the potential to provide precise and reproducible anatomical views, which can reduce overall execution time, labor burden of procedures, and x-ray usage for a range of cardiac procedures. Keywords: Automated View Recovery, Path Planning, Intra-cardiac echocardiography (ICE), Catheter, Tendon-driven manipulator, Cardiac Imaging


title: Neutron transfer reactions on the ground state and isomeric state of a 130Sn beam

abstract: The structure of nuclei around the neutron-rich nucleus 132Sn is of particular interest due to the vicinity of the Z = 50 and N = 82 shell closures and the r-process nucleosynthetic path. Four states in 131Sn with a strong single-particle-like component have previously been studied via the (d,p) reaction, with limited excitation energy resolution. The 130Sn(9Be,8Be)131Sn and 130Sn(13C,12C)131Sn single-neutron transfer reactions were performed in inverse kinematics at the Holifield Radioactive Ion Beam Facility using particle-gamma coincidence spectroscopy. The uncertainties in the energies of the single-particle-like states have been reduced by more than an order of magnitude using the energies of gamma rays. The previous tentative Jpi values have been confirmed. Decays from high-spin states in 131Sn have been observed following transfer on the isomeric component of the 130Sn beam. The improved energies and confirmed spin-parities of the p-wave states important to the r-process lead to direct-semidirect cross-sections for neutron capture on the ground state of 130Sn at 30 keV that are in agreement with previous analyses. A similar assessment of the impact of neutron-transfer on the isomer would require significant nuclear structure and reaction theory input. There are few measurements of transfer reaction on isomers, and this is the first on an isomer in the 132Sn region.


title: Solution of an Acoustic Transmission Inverse Problem by Extended Inversion

abstract: Study of a simple single-trace transmission example shows how an extended source formulation of full-waveform inversion can produce an optimization problem without spurious local minima ("cycle skipping"), hence efficiently solvable via Newton-like local optimization methods. The data consist of a single trace extracted from a causal pressure field, propagating in a homogeneous fluid according to linear acoustics, and recorded at a given distance from a transient point energy source. The source intensity ("wavelet") is presumed quasi-impulsive, with zero energy for time lags greater than a specified maximum lag. The inverse problem is: from the recorded trace, recover both the sound velocity or slowness and source wavelet with specified support, so that the data is fit with prescribed RMS relative error. The least-squares objective function has multiple large residual minimizers. The extended inverse problem permits source energy to spread in time, and replaces the maximum lag constraint by a weighted quadratic penalty. A companion paper shows that for proper choice of weight operator, any stationary point of the extended objective produces a good approximation of the global minimizer of the least squares objective, with slowness error bounded by a multiple of the maximum lag and the assumed noise level. This paper summarizes the theory developed in the companion paper and presents numerical experiments demonstrating the accuracy of the predictions in concrete instances. We also show how to dynamically adjust the penalty scale during iterative optimization to improve the accuracy of the slowness estimate.


title: Computational Modelling of Russia's First 2G-HTS Triaxial Cable

abstract: A better understanding of the interaction between three phases is required when developing superconducting cables for high voltage AC systems. With a particular focus on the energy losses of real power transmission cables, in this paper we utilize the so-called H-formulation of Maxwell equations to devise a 2D model for superconducting triaxial cables. The major aim of this model is to comprehend and reproduce the experimental observations reported on the first triaxial prototype cable developed by SuperOx and VNIIKP. The computationally modelled and prototyped cable is made of up to 87 tapes of 4 mm width SuperOx tape arranged across the three phases. Our computational results are compared to the experimental measurements performed by VNIIKP with the electrical contact method, showing a high degree of accuracy over the outer phase of the cable, whilst revealing technical issues with the experimental measurements at the inner phases. Thus, in consultation with VNIIKP it has been concluded that for the actual experimental measurement of the AC losses at the inner phases, and consequently of the overall cable, a sophisticated calorimetric setup must be built. Still our model is capable to provide an independent assessment of the VNIIKP-SuperOx cable design, by investigating the magnetic profiles per phase in the time domain. In this sense, we confirm that the unbalanced arrange of currents and distancing between the phases affirmatively lead to no magnetic leakages, and therefore to an adequate balance of the cabling inductance.


title: Signal Strength and Noise Drive Feature Preference in CNN Image Classifiers

abstract: Feature preference in Convolutional Neural Network (CNN) image classifiers is integral to their decision making process, and while the topic has been well studied, it is still not understood at a fundamental level. We test a range of task relevant feature attributes (including shape, texture, and color) with varying degrees of signal and noise in highly controlled CNN image classification experiments using synthetic datasets to determine feature preferences. We find that CNNs will prefer features with stronger signal strength and lower noise irrespective of whether the feature is texture, shape, or color. This provides guidance for a predictive model for task relevant feature preferences, demonstrates pathways for bias in machine models that can be avoided with careful controls on experimental setup, and suggests that comparisons between how humans and machines prefer task relevant features in vision classification tasks should be revisited. Code to reproduce experiments in this paper can be found at \url{https://github.com/mwolff31/signal_preference}.


title: Reinforcement Learning for Personalized Drug Discovery and Design for Complex Diseases: A Systems Pharmacology Perspective

abstract: Many multi-genic systemic diseases such as neurological disorders, inflammatory diseases, and the majority of cancers do not have effective treatments yet. Reinforcement learning powered systems pharmacology is a potentially effective approach to design personalized therapies for untreatable complex diseases. In this survey, state-of-the-art reinforcement learning methods and their latest applications to drug design are reviewed. The challenges on harnessing reinforcement learning for systems pharmacology and personalized medicine are discussed. Potential solutions to overcome the challenges are proposed. In spite of successful application of advanced reinforcement learning techniques to target-based drug discovery, new reinforcement learning strategies are needed to address systems pharmacology-oriented personalized de novo drug design.


title: On the Satisfaction Probability of $k$-CNF Formulas

abstract: The satisfaction probability Pr[$\phi$] := Pr$_{\beta:vars(\phi) \to \{0,1\}}[\beta\models \phi]$ of a propositional formula $\phi$ is the likelihood that a random assignment $\beta$ makes the formula true. We study the complexity of the problem $k$SAT-PROB$_{>p}$ = {$\phi$ is a $k$CNF formula | Pr[$\phi$] > p} for fixed $k$ and $p$. While 3SAT-PROB$_{>0}$ = 3SAT is NP-complete and SAT-PROB$_{>1/2}$ is PP-complete, Akmal and Williams recently showed that 3SAT-PROB$_{>1/2}$ lies in P and 4SAT-PROB$_{>1/2}$ is NP-complete; but the methods used to prove these striking results stay silent about, say, 4SAT-PROB$_{>3/4}$, leaving the computational complexity of $k$SAT-PROB$_{>p}$ open for most $k$ and $p$. In the present paper we give a complete characterization in the form of a trichotomy: $k$SAT-PROB$_{>p}$ lies in AC$^0$, is NL-complete, or is NP-complete; and given $k$ and $p$ we can decide which of the three applies. The proof of the trichotomy hinges on a new order-theoretic insight: Every set of $k$CNF formulas contains a formula of maximal satisfaction probability. This deceptively simple statement allows us to (1) kernelize $k$SAT-PROB$_{\ge p}$, (2) show that the variables of the kernel form a backdoor set when the trichotomy states membership in AC$^0$ or NL, and (3) prove a locality property by which Pr[$\phi$] < $p$ implies, for $k$CNF formulas $\phi$, that Pr[$\psi$] < $p$ holds already for some fixed-size subset $\psi$ of $\phi$'s clauses.


title: Environment Generation for Zero-Shot Compositional Reinforcement Learning

abstract: Many real-world problems are compositional - solving them requires completing interdependent sub-tasks, either in series or in parallel, that can be represented as a dependency graph. Deep reinforcement learning (RL) agents often struggle to learn such complex tasks due to the long time horizons and sparse rewards. To address this problem, we present Compositional Design of Environments (CoDE), which trains a Generator agent to automatically build a series of compositional tasks tailored to the RL agent's current skill level. This automatic curriculum not only enables the agent to learn more complex tasks than it could have otherwise, but also selects tasks where the agent's performance is weak, enhancing its robustness and ability to generalize zero-shot to unseen tasks at test-time. We analyze why current environment generation techniques are insufficient for the problem of generating compositional tasks, and propose a new algorithm that addresses these issues. Our results assess learning and generalization across multiple compositional tasks, including the real-world problem of learning to navigate and interact with web pages. We learn to generate environments composed of multiple pages or rooms, and train RL agents capable of completing wide-range of complex tasks in those environments. We contribute two new benchmark frameworks for generating compositional tasks, compositional MiniGrid and gMiniWoB for web navigation.CoDE yields 4x higher success rate than the strongest baseline, and demonstrates strong performance of real websites learned on 3500 primitive tasks.


title: Congruence frames of frames and $\kappa$-frames

abstract: We describe the congruence lattices of frames and $\kappa$-frames. We look at the role that congruence biframes play in the category of strictly zero-dimensional biframes and discuss some reflections and coreflections of congruence frames.


title: The Weak Gram Law for Hecke $L$-functions

abstract: We generalize a theorem by Titchmarsh about the mean value of Hardy's $Z$-function at the Gram points to the Hecke $L$-functions, which in turn implies the weak Gram law for them. Instead of proceeding analogously to Titchmarsh with an approximate functional equation we employ a different method using contour integration.


title: Scaling limits of loop-erased Markov chains on resistance spaces via a partial loop-erasing procedure

abstract: We introduce partial loop-erasing operators. We show that by applying a refinement sequence of partial loop-erasing operators to a finite Markov chain, we get a process equivalent to the chronological loop-erased Markov chain. As an application, we construct loop-erased random paths on bounded domains of resistance spaces as the weak limit of the loop erasure of the Markov chains on a sequence of finite sets approximating the space, and the limit is independent of the approximating sequences. The random paths we constructed are simple paths almost surely, and they can be viewed as the loop-erasure of the paths of the diffusion process. Finally, we show that the scaling limit of the loop-erased random walks on the Sierpi\'nski carpet graphs exists, and is equivalent to the loop-erased random paths on the Sierpi\'nksi carpet.


title: Approximating the discrete and continuous median line segments in $d$ dimensions

abstract: Consider a set $P$ of $n$ points in $\mathbb{R}^d$. In the discrete median line segment problem, the objective is to find a line segment bounded by a pair of points in $P$ such that the sum of the Euclidean distances from $P$ to the line segment is minimized. In the continuous median line segment problem, a real number $\ell>0$ is given, and the goal is to locate a line segment of length $\ell$ in $\mathbb{R}^d$ such that the sum of the Euclidean distances between $P$ and the line segment is minimized. We show how to compute $(1+\epsilon\Delta)$- and $(1+\epsilon)$-approximations to a discrete median line segment in time $O(n\epsilon^{-2d}\log n)$ and $O(n^2\epsilon^{-d})$, respectively, where $\Delta$ is the spread of line segments spanned by pairs of points. While developing our algorithms, by using the principle of pair decomposition, we derive new data structures that allow us to quickly approximate the sum of the distances from a set of points to a given line segment or point. To our knowledge, our utilization of pair decompositions for solving minsum facility location problems is the first of its kind; it is versatile and easily implementable. We prove that it is impossible to construct a continuous median line segment for $n\geq3$ non-collinear points in the plane by using only ruler and compass. In view of this, we present an $O(n^d\epsilon^{-d})$-time algorithm for approximating a continuous median line segment in $\mathbb{R}^d$ within a factor of $1+\epsilon$. The algorithm is based upon generalizing the point-segment pair decomposition from the discrete to the continuous domain. Last but not least, we give an $(1+\epsilon)$-approximation algorithm, whose time complexity is sub-quadratic in $n$, for solving the constrained median line segment problem in $\mathbb{R}^2$ where an endpoint or the slope of the median line segment is given at input.


title: An Ensemble Model for Face Liveness Detection

abstract: In this paper, we present a passive method to detect face presentation attack a.k.a face liveness detection using an ensemble deep learning technique. Face liveness detection is one of the key steps involved in user identity verification of customers during the online onboarding/transaction processes. During identity verification, an unauthenticated user tries to bypass the verification system by several means, for example, they can capture a user photo from social media and do an imposter attack using printouts of users faces or using a digital photo from a mobile device and even create a more sophisticated attack like video replay attack. We have tried to understand the different methods of attack and created an in-house large-scale dataset covering all the kinds of attacks to train a robust deep learning model. We propose an ensemble method where multiple features of the face and background regions are learned to predict whether the user is a bonafide or an attacker.


title: Transmission Estimation at the Fundamental Quantum Cram\'er-Rao Bound with Macroscopic Quantum Light

abstract: The field of quantum metrology seeks to apply quantum techniques and/or resources to classical sensing approaches with the goal of enhancing the precision in the estimation of a parameter beyond what can be achieved with classical resources. Theoretically, the fundamental minimum uncertainty in the estimation of a parameter for a given probing state is bounded by the quantum Cram\'er-Rao bound. From a practical perspective, it is necessary to find physical measurements that can saturate this fundamental limit and to show experimentally that it is possible to perform measurements with the required precision to do so. Here we perform experiments that saturate the quantum Cram\'er-Rao bound for transmission estimation over a wide range of transmissions when probing the system under study with a bright two-mode squeezed state. To properly take into account the imperfections in the generation of the quantum state, we extend our previous theoretical results to incorporate the measured properties of the generated quantum state. For our largest transmission level of 84%, we show a 62% reduction over the optimal classical protocol in the variance in transmission estimation when probing with a bright two-mode squeezed state with 8 dB of intensity-difference squeezing. Given that transmission estimation is an integral part of many sensing protocols, such as plasmonic sensing, spectroscopy, calibration of the quantum efficiency of detectors, etc., the results presented promise to have a significant impact on a number of applications in various fields of research.


title: Universal Online Learning with Unbounded Losses: Memory Is All You Need

abstract: We resolve an open problem of Hanneke on the subject of universally consistent online learning with non-i.i.d. processes and unbounded losses. The notion of an optimistically universal learning rule was defined by Hanneke in an effort to study learning theory under minimal assumptions. A given learning rule is said to be optimistically universal if it achieves a low long-run average loss whenever the data generating process makes this goal achievable by some learning rule. Hanneke posed as an open problem whether, for every unbounded loss, the family of processes admitting universal learning are precisely those having a finite number of distinct values almost surely. In this paper, we completely resolve this problem, showing that this is indeed the case. As a consequence, this also offers a dramatically simpler formulation of an optimistically universal learning rule for any unbounded loss: namely, the simple memorization rule already suffices. Our proof relies on constructing random measurable partitions of the instance space and could be of independent interest for solving other open questions. We extend the results to the non-realizable setting thereby providing an optimistically universal Bayes consistent learning rule.


title: Description-Driven Task-Oriented Dialog Modeling

abstract: Task-oriented dialogue (TOD) systems are required to identify key information from conversations for the completion of given tasks. Such information is conventionally specified in terms of intents and slots contained in task-specific ontology or schemata. Since these schemata are designed by system developers, the naming convention for slots and intents is not uniform across tasks, and may not convey their semantics effectively. This can lead to models memorizing arbitrary patterns in data, resulting in suboptimal performance and generalization. In this paper, we propose that schemata should be modified by replacing names or notations entirely with natural language descriptions. We show that a language description-driven system exhibits better understanding of task specifications, higher performance on state tracking, improved data efficiency, and effective zero-shot transfer to unseen tasks. Following this paradigm, we present a simple yet effective Description-Driven Dialog State Tracking (D3ST) model, which relies purely on schema descriptions and an "index-picking" mechanism. We demonstrate the superiority in quality, data efficiency and robustness of our approach as measured on the MultiWOZ (Budzianowski et al.,2018), SGD (Rastogi et al., 2020), and the recent SGD-X (Lee et al., 2021) benchmarks.


title: Optimal Dynamic Regret in Proper Online Learning with Strongly Convex Losses and Beyond

abstract: We study the framework of universal dynamic regret minimization with strongly convex losses. We answer an open problem in Baby and Wang 2021 by showing that in a proper learning setup, Strongly Adaptive algorithms can achieve the near optimal dynamic regret of $\tilde O(d^{1/3} n^{1/3}\text{TV}[u_{1:n}]^{2/3} \vee d)$ against any comparator sequence $u_1,\ldots,u_n$ simultaneously, where $n$ is the time horizon and $\text{TV}[u_{1:n}]$ is the Total Variation of comparator. These results are facilitated by exploiting a number of new structures imposed by the KKT conditions that were not considered in Baby and Wang 2021 which also lead to other improvements over their results such as: (a) handling non-smooth losses and (b) improving the dimension dependence on regret. Further, we also derive near optimal dynamic regret rates for the special case of proper online learning with exp-concave losses and an $L_\infty$ constrained decision set.


title: Machine learning interatomic potential for high throughput screening and optimization of high-entropy alloys

abstract: We have developed a machine learning-based interatomic potential (MLIP) for the quaternary MoNbTaW (R4) and quinary MoNbTaTiW (R5) high entropy alloys (HEAs). MLIPs enabled accurate high throughput calculations of elastic and mechanical properties of various non-equimolar R4 and R5 alloys, which are otherwise very time consuming calculations when performed using density functional theory (DFT). We demonstrate that the MLIP predicted properties compare well with the DFT results on various test cases and are consistent with the available experimental data. The MLIPs are also utilized for high throughput optimization of non equimolar R4 candidates by guided iterative tuning of R4 compositions to discover candidate materials with promising hardness-ductility combinations. We also used this approach to study the effect of Ti concentration on the elastic and mechanical properties of R4, by statistically averaging the properties of over 100 random structures. MLIP predicted hardness and bulk modulus of equimolar R4 and R5 HEAs are validated using experimentally measured Vickers hardness and modulus. This approach opens a new avenue for employing MLIPs for HEA candidate optimization.


title: Linear lexicographic optimization and preferential bidding system

abstract: Some airlines use the preferential bidding system to construct the schedules of their pilots. In this system, the pilots bid on the different activities and the schedules that lexicographically maximize the scores of the pilots according to their seniority are selected. A sequential approach to solve this maximization problem is natural: the problem is first solved with the bids of the most senior pilot; then it is solved with those of the second most senior without decreasing the score of the most senior, and so on. The literature admits that the structure of the problem somehow imposes such an approach. The problem can be modeled as an integer linear lexicographic program. We propose a new exact method, which relies on column generation for solving its continuous relaxation. To design this column generation, we prove that bounded linear lexicographic programs admit "primal-dual" feasible bases and we show how to compute such bases efficiently. Another contribution on which our exact method relies consists in the extension of standard tools for resource-constrained longest path problems to their lexicographic versions. This is useful in our context since the generation of new columns is modeled as a lexicographic resource-constrained longest path problem. Numerical experiments show that this new method is already able to solve industrial instances provided by Air France, with up to 150 pilots. By adding a last ingredient in the resolution of the longest path problems, which exploits the specificity of the preferential bidding system, the method achieves for these instances computational times that are compatible with operational constraints.


title: The 334-Triangle Graph of $SL_3({\mathbb Z})$

abstract: Long, Reid, and Thistlewaite have shown that some groups generated by representations of the $\Delta 334$ triangle group in $SL_3({\mathbb Z})$ are thin, while the status of others is unknown. In this paper we take a new approach: for each group we introduce a new graph that captures information about representations of $\Delta 334$ in the group. We provide examples of our graph for a variety of groups, and we use information about the graph for $SL_3({\mathbb Z}/2{\mathbb Z})$ to show that the chromatic number of the graph for $SL_3({\mathbb Z})$ is at most eight. By generating a portion of the graph for $SL_3({\mathbb Z})$ we show its chromatic number is at least four; we conjecture it is equal to four.


title: Uncertainty-Cognizant Model Predictive Control for Energy Management of Residential Buildings with PVT and Thermal Energy Storage

abstract: The building sector accounts for almost 40 percent of the global energy consumption. This reveals a great opportunity to exploit renewable energy resources in buildings to achieve the climate target. In this context, this paper offers a building energy system embracing a heat pump, a thermal energy storage system along with grid-connected photovoltaic thermal (PVT) collectors to supply both electric and thermal energy demands of the building with minimum operating cost. To this end, the paper develops a stochastic model predictive control (MPC) strategy to optimally determine the set-point of the whole building energy system while accounting for the uncertainties associated with the PVT energy generation. This system enables the building to 1-shift its electric demand from high-peak to off-peak hours and 2- sell electricity to the grid to make energy arbitrage.


title: A Systematic Exploration of Reservoir Computing for Forecasting Complex Spatiotemporal Dynamics

abstract: A reservoir computer (RC) is a type of simplified recurrent neural network architecture that has demonstrated success in the prediction of spatiotemporally chaotic dynamical systems. A further advantage of RC is that it reproduces intrinsic dynamical quantities essential for its incorporation into numerical forecasting routines such as the ensemble Kalman filter -- used in numerical weather prediction to compensate for sparse and noisy data. We explore here the architecture and design choices for a "best in class" RC for a number of characteristic dynamical systems, and then show the application of these choices in scaling up to larger models using localization. Our analysis points to the importance of large scale parameter optimization. We also note in particular the importance of including input bias in the RC design, which has a significant impact on the forecast skill of the trained RC model. In our tests, the the use of a nonlinear readout operator does not affect the forecast time or the stability of the forecast. The effects of the reservoir dimension, spinup time, amount of training data, normalization, noise, and the RC time step are also investigated. While we are not aware of a generally accepted best reported mean forecast time for different models in the literature, we report over a factor of 2 increase in the mean forecast time compared to the best performing RC model of Vlachas et.al (2020) for the 40 dimensional spatiotemporally chaotic Lorenz 1996 dynamics, and we are able to accomplish this using a smaller reservoir size.


title: The forward and reverse shock dynamics of Cassiopeia A

abstract: We report on proper motion measurements of the forward- and reverse-shock regions of the supernova remnant Cassiopeia A (Cas A), including deceleration/acceleration measurements of the forward shock. The measurements combine 19 years of observations with the Chandra X-ray Observatory, using the 4.2-6 keV continuum band, preferentially targeting X-ray synchrotron radiation. The average expansion rate is $0.218 \pm 0.029$%yr$^{-1}$ for the forward shock, corresponding to a velocity of $\approx 5800$ km/s. The time derivative of the proper motions indicates deceleration in the east, and an acceleration up to $1.1\times 10^{-4}$ yr$^{-2}$ in the western part. The reverse shock moves outward in the East, but in the West it moves toward the center with an expansion rate of $-0.0225\pm 0.0007$ %yr$^{-1}$, corresponding to $-1884\pm 17$ km/s. In the West the reverse shock velocity in the ejecta frame is $\gtrsim 3000$ km/s, peaking at $\sim 8000$ km/s, explaining the presence of X-ray synchrotron emitting filaments there. The backward motion of the reverse shock can be explained by either a scenario in which the forward shock encountered a partial, dense, wind shell, or one in which the shock transgressed initially through a lopsided cavity, created during a brief Wolf-Rayet star phase. Both scenarios are consistent with the local acceleration of the forward shock. Finally we report on the proper motion of the northeastern jet, using both the X-ray continuum band, and the Si XIII K-line emission band. We find expansion rates of respectively 0.21%yr$^{-1}$ and 0.24%yr$^{-1}$, corresponding to velocities at the tip of the X-ray jet of 7830--9200 km/s.


title: Sparse grid implementation of a fixed-point fast sweeping WENO scheme for Eikonal equations

abstract: Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady state solutions of hyperbolic partial differential equations (PDEs). As other types of fast sweeping schemes, fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order. The resulting iterative schemes have fast convergence rate to steady state solutions. Moreover, an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve inverse operation of any nonlinear local system. Hence they are robust and flexible, and have been combined with high order accurate weighted essentially non-oscillatory (WENO) schemes to solve various hyperbolic PDEs in the literature. For multidimensional nonlinear problems, high order fixed-point fast sweeping WENO methods still require quite large amount of computational costs. In this technical note, we apply sparse-grid techniques, an effective approximation tool for multidimensional problems, to fixed-point fast sweeping WENO method for reducing its computational costs. Here we focus on a robust Runge-Kutta (RK) type fixed-point fast sweeping WENO scheme with third order accuracy (Zhang et al. 2006 [33]), for solving Eikonal equations, an important class of static Hamilton-Jacobi (H-J) equations. Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse grid computations of the fixed-point fast sweeping WENO scheme achieve large savings of CPU times on refined meshes, and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.


title: A Computation of the Action of the Morava Stabilizer Group on the Lubin-Tate Deformation Ring

abstract: We compute recursive approximations of the action of the height $h \geq 2$ Morava stabilizer group on the associated Lubin-Tate deformation ring. We then specialize to the case $h=3$ and $p>2$ to calculate the action explicitly. These results are new for $h>2$ and agree with computations by Lader at height $h=2$.


title: Numerical Analysis of a Corrected Smagorinsky Model

abstract: The classical Smagorinsky model's solution is an approximation to a (resolved) mean velocity. Since it is an eddy viscosity model, it cannot represent a flow of energy from unresolved fluctuations to the (resolved) mean velocity. This model has recently been modified to incorporate this flow and still be well-posed. Herein we first develop some basic properties of the modified model. Next, we perform a complete numerical analysis of two algorithms for its approximation. They are tested and proven to be effective.


title: New elements in the center of free alternative algebra

abstract: A new series of central elements is found in the free alternative algebra. More exactly, let $Alt[X]$ and $SMalc[X]\subset Alt[X]$ be the free alternative algebra and the free special Malcev algebra over a field of characteristic 0 on a set of free generators $X$, and let $f(x,y,x_1,\ldots,x_n)\in SMalc[X]$ be a multilinear element which is trivial in the free associative algebra. Then the element $u_n=u_n(x,x_1,\ldots,x_n)=f(x^2,x,x_1,\ldots,x_n)-f(x,x^2,x_1,\ldots,x_n)$ lies in the center of the algebra $Alt[X]$. The elements $u_n(x,x_1,\ldots,x_n)$ are uniquely defined up to a scalar for a given $n$, and they are skew-symmetric on the variables $x_1,\ldots,x_n$. Moreover, $u_n=0$ for $n=4m+2,\,4m+3$. and $u_n\neq 0$ for $n=4m,4m+1$. The ideals generated by the elements $u_{4m},\,u_{4m+1}$ lie in the associative center of the algebra $Alt[X]$ and have trivial multiplication.


title: Enabling Flexibility for Sparse Tensor Acceleration via Heterogeneity

abstract: Recently, numerous sparse hardware accelerators for Deep Neural Networks (DNNs), Graph Neural Networks (GNNs), and scientific computing applications have been proposed. A common characteristic among all of these accelerators is that they target tensor algebra (typically matrix multiplications); yet dozens of new accelerators are proposed for every new application. The motivation is that the size and sparsity of the workloads heavily influence which architecture is best for memory and computation efficiency. To satisfy the growing demand of efficient computations across a spectrum of workloads on large data centers, we propose deploying a flexible 'heterogeneous' accelerator, which contains many 'sub-accelerators' (smaller specialized accelerators) working together. To this end, we propose: (1) HARD TACO, a quick and productive C++ to RTL design flow to generate many types of sub-accelerators for sparse and dense computations for fair design-space exploration, (2) AESPA, a heterogeneous sparse accelerator design template constructed with the sub-accelerators generated from HARD TACO, and (3) a suite of scheduling strategies to map tensor kernels onto heterogeneous sparse accelerators with high efficiency and utilization. AESPA with optimized scheduling achieves 1.96X higher performance, and 7.9X better energy-delay product (EDP) than a Homogeneous EIE-like accelerator with our diverse workload suite.


title: Frequency-dependent Inter-pseudospin Solutions to Superconducting Strontium Ruthenate

abstract: The lasting puzzle of the superconducting order parameter of Sr$_2$RuO$_4$ calls for theoretical studies that include seldom-considered effects. Here we include spin-orbit coupling effects on the electronic structure and then solve the linearized Eliashberg equation in a pseudospin basis, including the possibility that spin and charge fluctuations induce frequency-dependent superconducting order parameters. We find that spin-orbit coupling mixes even and odd contributions in orbital, spin and frequency spaces and that leading inter-pseudospin symmetries, B$_{1g}^+$ and A$_{2g}^-$, have intra-orbital components respectively even and odd in Matsubara frequency. An accidental degeneracy between these could resolve apparent experimental contradictions.

categories: cond-mat.supr-con cond-mat.mtrl-sci cond-mat.str-el quant-ph


title: Superconductivity in Correlated Multi-Orbital Systems with Spin-Orbit Coupling: Coexistence of Even- and Odd-Frequency Pairing and the Case of Strontium Ruthenate

abstract: The superconducting order parameter of strontium ruthenate is the center of a lasting puzzle calling for theoretical studies that include the seldom-considered effects of spin-orbit coupling and the frequency-dependence of the order parameters. Here we generalize the frequency-dependent theory of superconductivity mediated by spin and charge fluctuations to include spin-orbit coupling in multi-orbital systems and we characterize the superconducting states using the spin-parity-orbital-time $SPOT$ quantum numbers, group theory, and phase distributions in the complex plane. We derive a pseudospin formulation that maps the inter-pseudospin sector of the normal state Eliashberg equation to a pseudospin-diagonal one. Possible superconducting order parameters for strontium ruthenate are obtained starting from a realistic density-functional-theory normal state. We find that spin-orbit coupling leads to ubiquitous entanglement of spin and orbital quantum numbers, along with notable mixing between even- and odd-frequency correlations. We propose a phase diagram obtained from the temperature dependence of the leading and subleading symmetries in the pseudospin-orbital basis. An accidental degeneracy between leading inter-pseudospin symmetries in strontium ruthenate, B$_{1g}^+$ and A$_{2g}^-$, could resolve apparent experimental contradictions.