title: Timing the earliest quenching events with a robust sample of massive quiescent galaxies at 2 < z < 5

abstract: We present a sample of 151 massive ($M_* > 10^{10}\mathrm{M_\odot}$) quiescent galaxies at $2 < z < 5$, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at $z>3$, of which 2 are at $z>4$. We report formation redshifts, demonstrating that the oldest objects formed at $z > 6$, however individual ages from our photometric data have significant uncertainties, typically $\sim0.5$ Gyr. We demonstrate that the UVJ colours of the quiescent population evolve with redshift at $z>3$, becoming bluer and more similar to post-starburst galaxies at lower redshift. Based upon this we construct a model for the time-evolution of quiescent galaxy UVJ colours, concluding that the oldest objects are consistent with forming the bulk of their stellar mass at $z\sim6-7$ and quenching at $z\sim5$. We report spectroscopic redshifts for two of our objects at $z=3.440$ and $3.396$, which exhibit extremely weak Ly$\alpha$ emission in ultra-deep VANDELS spectra. We calculate star-formation rates based on these line fluxes, finding that these galaxies are consistent with our quiescent selection criteria, provided their Ly$\alpha$ escape fractions are $>3$ and $>10$ per cent respectively. We finally report that our highest-redshift robust object exhibits a continuum break at $\lambda\sim7000$A in a spectrum from VUDS, consistent with our photometric redshift of $z_\mathrm{phot}=4.72^{+0.06}_{-0.04}$. If confirmed as quiescent this object would be the highest-redshift known quiescent galaxy. To obtain stronger constraints on the times of the earliest quenching events, high-SNR spectroscopy must be extended to $z\gtrsim3$ quiescent objects.


title: Continuous Emotion Recognition via Deep Convolutional Autoencoder and Support Vector Regressor

abstract: Automatic facial expression recognition is an important research area in the emotion recognition and computer vision. Applications can be found in several domains such as medical treatment, driver fatigue surveillance, sociable robotics, and several other human-computer interaction systems. Therefore, it is crucial that the machine should be able to recognize the emotional state of the user with high accuracy. In recent years, deep neural networks have been used with great success in recognizing emotions. In this paper, we present a new model for continuous emotion recognition based on facial expression recognition by using an unsupervised learning approach based on transfer learning and autoencoders. The proposed approach also includes preprocessing and post-processing techniques which contribute favorably to improving the performance of predicting the concordance correlation coefficient for arousal and valence dimensions. Experimental results for predicting spontaneous and natural emotions on the RECOLA 2016 dataset have shown that the proposed approach based on visual information can achieve CCCs of 0.516 and 0.264 for valence and arousal, respectively.


title: Macroscopic loops in the loop O(n) model via the XOR trick

abstract: The loop $O(n)$ model is a family of probability measures on collections of non-intersecting loops on the hexagonal lattice, parameterized by a loop-weight $n$ and an edge-weight $x$. Nienhuis predicts that, for $0 \leq n \leq 2$, the model exhibits two regimes separated by $x_c(n) = 1/\sqrt{2 + \sqrt{2-n}}$: when $x < x_c(n)$, the loop lengths have exponential tails, while, when $x \geq x_c(n)$, the loops are macroscopic. In this paper, we prove three results regarding the existence of long loops in the loop $O(n)$ model: - In the regime $(n,x) \in [1,1+\delta) \times (1- \delta, 1]$ with $\delta >0$ small, a configuration sampled from a translation-invariant Gibbs measure will either contain an infinite path or have infinitely many loops surrounding every face. In the subregime $n \in [1,1+\delta)$ and $x \in (1-\delta,1/\sqrt{n}]$ our results further imply Russo--Seymour--Welsh theory. This is the first proof of the existence of macroscopic loops in a positive area subset of the phase diagram. - Existence of loops whose diameter is comparable to that of a finite domain whenever $n=1, x \in (1,\sqrt{3}]$; this regime is equivalent to part of the antiferromagnetic regime of the Ising model on the triangular lattice. - Existence of non-contractible loops on a torus when $n \in [1,2], x=1$. The main ingredients of the proof are: (i) the `XOR trick': if $\omega$ is a collection of short loops and $\Gamma$ is a long loop, then the symmetric difference of $\omega$ and $\Gamma$ necessarily includes a long loop as well; (ii) a reduction of the problem of finding long loops to proving that a percolation process on an auxiliary planar graph, built using the Chayes--Machta and Edwards--Sokal geometric expansions, has no infinite connected components; and (iii) a recent result on the percolation threshold of Benjamini--Schramm limits of planar graphs.


title: An augmented filled function for global nonlinear integer optimization

abstract: The problem of finding global minima of nonlinear discrete functions arises in many fields of practical matters. In recent years, methods based on discrete filled functions become popular as ways of solving these sort of problems. However, they rely on the steepest descent method for local searches. Here we present an approach that does not depend on a particular local optimization method, and a new discrete filled function with the useful property that a good continuous global optimization algorithm applied to it leads to an approximation of the solution of the nonlinear discrete problem. Numerical results are given showing the efficiency of the new approach.


title: Complete Leibniz Algebras

abstract: Leibniz algebras are certain generalization of Lie algebras. It is natural to generalize concepts in Lie algebras to Leibniz algebras and investigate whether the corresponding results still hold. In this paper we introduce the notion of complete Leibniz algebras as generalization of complete Lie algebras. Then we study properties of complete Leibniz algebras and their holomorphs.


title: A cold atom quantum simulator to explore pairing, condensation, and pseudogaps in extended Hubbard--Holstein models

abstract: We describe a quantum simulator for the Hubbard--Holstein model (HHM), comprising two dressed Rydberg atom species held in a monolayer by independent painted potentials, predicting that boson-mediated preformed pairing, and Berezinskii--Kosterlitz--Thouless (BKT) transition temperatures are experimentally accessible. The HHM is important for modeling the essential physics of unconventional superconductors. Experimentally realizable quantum simulators for HHMs are needed: (1) since HHMs are difficult to solve numerically and analytically, (2) to explore how competition between electron-phonon interactions and strong repulsion affects pairing in unconventional superconductors, (3) to understand the role of boson-mediated local pairing in pseudogaps and fermion condensates. We propose and study a quantum simulator for the HHM, using optical lattices, painted using zeros in the AC stark shift, to control two Rydberg atom species independently within a monolayer. We predict that interactions are sufficiently tunable to probe: (1) both HHMs and highly unconventional phonon-mediated repulsions, (2) the competition between intermediate-strength phonon and Coulomb mediated interactions, (3) BKT transitions, and preformed pairing that could be used to examine key hypotheses related to the pseudogap. We discuss how the quantum simulator can be used to investigate boson-mediated pairing and condensation of fermions in unconventional superconductors.


title: Lifted Reed-Solomon Codes with Application to Batch Codes

abstract: Guo, Kopparty and Sudan have initiated the study of error-correcting codes derived by lifting of affine-invariant codes. Lifted Reed-Solomon (RS) codes are defined as the evaluation of polynomials in a vector space over a field by requiring their restriction to every line in the space to be a codeword of the RS code. In this paper, we investigate lifted RS codes and discuss their application to batch codes, a notion introduced in the context of private information retrieval and load-balancing in distributed storage systems. First, we improve the estimate of the code rate of lifted RS codes for lifting parameter $m\ge 3$ and large field size. Second, a new explicit construction of batch codes utilizing lifted RS codes is proposed. For some parameter regimes, our codes have a better trade-off between parameters than previously known batch codes.


title: Electric field control of interfacial Dzyaloshinskii-Moriya interaction in Pt/Co/AlO$_x$ thin films

abstract: We studied electric field modification of magnetic properties in a Pt/Co/AlO$_x$ trilayer via magneto-optical Kerr microscopy. We observed the spontaneous formation of labyrinthine magnetic domain structure due to thermally activated domain nucleation and propagation under zero applied magnetic field. A variation of the period of the labyrinthine structure under electric field is observed as well as saturation magnetization and magnetic anisotropy variations. Using an analytical formula of the stripe equilibrium width we estimate the variation of the interfacial Dzyaloshinskii-Moriya interaction under electric field as function of the exchange stiffness constant.


title: Cooperative Energy Scheduling for Microgrids under Peak Demand Energy Plans

abstract: A cooperative energy scheduling method is proposed that allows joint energy optimization for a group of microgrids to achieve cost savings that the microgrids could not achieve individually. The discussed microgrids may be commercial entities in a distribution network under utility electricity rate plans comprising both Time of Use (ToU) and peak demand charge. Defining a stable operation as a situation where all microgrids would be willing to participate, it is shown that under such rate plans and in particular due to the peak demand charge, a cost distribution that is seemingly fair does not necessarily result in a stable cooperation. These results are derived in this paper using concepts from cooperative games. It is therefore sought to devise a stable cost distribution algorithm that, while maximizing some measure of fairness among the participating microgrids, ensures they all benefit from their participation. A simple case study is presented that demonstrates fairness and stability aspects of the cooperation.


title: Contractibility of space of stability conditions on the projective plane via global dimension function

abstract: We compute the global dimension function $\mathrm{gldim}$ on the principal component $\mathrm{Stab}^{\dag}(\mathbb{P}^2)$ of the space of Bridgeland stability conditions on $\mathbb{P}^2$. It admits $2$ as the minimum value and the preimage $\mathrm{gldim}^{-1}(2)$ is contained in the closure $\bar{\mathrm{Stab}^{\mathrm{Geo}}(\mathbb{P}^2)}$ of the subspace consisting of geometric stability conditions. We show that $\mathrm{gldim}^{-1}[2,x)$ contracts to $\mathrm{gldim}^{-1}(2)$ for any real number $x\geq 2$ and that $\mathrm{gldim}^{-1}(2)$ is contractible.


title: Pretrained Transformers for Simple Question Answering over Knowledge Graphs

abstract: Answering simple questions over knowledge graphs is a well-studied problem in question answering. Previous approaches for this task built on recurrent and convolutional neural network based architectures that use pretrained word embeddings. It was recently shown that finetuning pretrained transformer networks (e.g. BERT) can outperform previous approaches on various natural language processing tasks. In this work, we investigate how well BERT performs on SimpleQuestions and provide an evaluation of both BERT and BiLSTM-based models in datasparse scenarios.


title: Capacity-achieving Polar-based LDGM Codes with Crowdsourcing Applications

abstract: In this paper we study codes with sparse generator matrices. More specifically, codes with a certain constraint on the weight of all the columns in the generator matrix are considered. The end result is the following. For any binary-input memoryless symmetric (BMS) channel and any epsilon > 2 epsilon*, where epsilon^* = \frac{1}{6}-\frac{5}{3}\log{\frac{4}{3}} \approx 0.085, we show an explicit sequence of capacity-achieving codes with all the column wights of the generator matrix upper bounded by (\log N)^{1+epsilon}, where N is the code block length. The constructions are based on polar codes. Applications to crowdsourcing are also shown.


title: How $Not$ to Compute a Fourier Transform

abstract: We revisit the Fourier transform of a Hankel function, of considerable importance in the theory of knife edge diffraction. Our approach is based directly upon the underlying Bessel equation, which admits manipulation into an alternate second order differential equation, one of whose solutions is precisely the desired transform, apart from an $\textit{a priori}$ unknown constant, and a second, undesired solution of logarithmic type. A modest amount of analysis is then required to exhibit that constant as having its proper value, and to purge the logarithmic accompaniment. The intervention of this analysis, which relies upon an interplay of asymptotic and close-in functional behaviors, prompts our somewhat ironic, mildly puckish caveat, our negation $\textit{not}$ in the title. In a concluding section we show that this same transform is still more readily exhibited as an easy byproduct of the inhomogeneous wave equation in two dimensions satisfied by the Green's function $\textit{G}$, itself proportional to a Hankel function. This latter discussion lapses of course into the argot of physicists and, in its r$\hat{{\rm{o}}}$le of a mere afterthought, makes no claim whatsoever to any originality.


title: Consensus-Based Optimization on the Sphere: Convergence to Global Minimizers and Machine Learning

abstract: We investigate the implementation of a new stochastic Kuramoto-Vicsek-type model for global optimization of nonconvex functions on the sphere. This model belongs to the class of Consensus-Based Optimization. In fact, particles move on the sphere driven by a drift towards an instantaneous consensus point, which is computed as a convex combination of particle locations, weighted by the cost function according to Laplace's principle, and it represents an approximation to a global minimizer. The dynamics is further perturbed by a random vector field to favor exploration, whose variance is a function of the distance of the particles to the consensus point. In particular, as soon as the consensus is reached the stochastic component vanishes. The main results of this paper are about the proof of convergence of the numerical scheme to global minimizers provided conditions of well-preparation of the initial datum. The proof combines previous results of mean-field limit with a novel asymptotic analysis, and classical convergence results of numerical methods for SDE. We present several numerical experiments, which show that the algorithm proposed in the present paper scales well with the dimension and is extremely versatile. To quantify the performances of the new approach, we show that the algorithm is able to perform essentially as good as ad hoc state of the art methods in challenging problems in signal processing and machine learning, namely the phase retrieval problem and the robust subspace detection.

categories: cs.LG cs.NA math.AP math.NA math.OC stat.ML


title: Hierarchy and co-evolution processes in urban systems

abstract: The concept of hierarchy in complex systems is tightly linked to co-evolutionary processes. We propose here to explore it in the case of the co-evolution between transportation networks and territories. More precisely, we extend a co-evolution model for systems of cities and infrastructure networks, and systematically study its behavior following specific hierarchy indicators we introduce. We show that population hierarchy and network hierarchy are tightly linked, but that a broad range of regimes can exist. Model exploration furthermore yields non-trivial stylized facts which can be taken into account for territorial planning on such long time scales with co-evolutionary processes.


title: Deontological Ethics By Monotonicity Shape Constraints

abstract: We demonstrate how easy it is for modern machine-learned systems to violate common deontological ethical principles and social norms such as "favor the less fortunate," and "do not penalize good attributes." We propose that in some cases such ethical principles can be incorporated into a machine-learned model by adding shape constraints that constrain the model to respond only positively to relevant inputs. We analyze the relationship between these deontological constraints that act on individuals and the consequentialist group-based fairness goals of one-sided statistical parity and equal opportunity. This strategy works with sensitive attributes that are Boolean or real-valued such as income and age, and can help produce more responsible and trustworthy AI.


title: Chiral Abelian gauge theories with few fermions

abstract: We construct chiral theories with the smallest number $n_\chi$ of Weyl fermions that form an anomaly-free set under various Abelian gauge groups. For the $U(1)$ group, where $n_\chi = 5$, we show that the general solution to the anomaly equations is a set of charges given by cubic polynomials in three integer parameters. For the $U(1) \times U(1)$ gauge group we find $n_\chi = 6$, and derive the general solution to the anomaly equations, in terms of 6 parameters. For $U(1) \times U(1) \times U(1)$ we show that $n_\chi = 8$, and present some families of solutions. These chiral gauge theories have potential applications to dark matter models, right-handed neutrino interactions, and other extensions of the Standard Model. As an example, we present a simple dark sector with a natural mass hierarchy between three dark matter components.


title: On the Holton-Lindzen-Plumb model for mean flow reversals in stratified fluids

abstract: The Holton-Lindzen-Plumb model describes the spontaneous emergence of mean flow reversals in stratified fluids. It has played a central role in understanding the quasi-biennial oscillation of equatorial winds in Earth's stratosphere and has arguably become a linchpin of wave-mean flow interaction theory in geophysical and astrophysical fluid dynamics. The derivation of the model's equation from primitive equations follows from several assumptions, including quasi-linear approximations, WKB expansion of the wavefield, simplifications of boundary layer terms, among others. Starting from the two-dimensional, non-rotating, Boussinesq equations, we present in this paper a self-consistent derivation of the Holton-Lindzen-Plumb model and show the existence of a distinguished limit for which all approximations remains valid. We furthermore discuss the important role of boundary conditions, and the relevance of this model to describe secondary bifurcations associated with a quasi-periodic route to chaos.


title: Iterative Motion Compensation reconstruction ultra-short TE(iMoCo UTE) for high resolution free breathing pulmonary MRI

abstract: Purpose, To develop a high scanning efficiency, motion corrected imaging strategy for free-breathing pulmonary MRI by combining a motion compensation reconstruction with a UTE acquisition, called iMoCo UTE. Methods, An optimized golden angle ordering radial UTE sequence was used to continuously acquire data for 5 minutes. All readouts were grouped to different respiratory motion states based on self-navigator signals, then motion resolved data was reconstructed by XD Golden angle RAdial Sparse Parallel reconstruction (XD GRASP). One state from the motion resolved images was selected as a reference, and motion fields from the other states to the reference were derived via non-rigid registration. Finally, all motion resolved data and motion fields were reconstructed by using an iterative motion compensation reconstruction with a total generalized variation sparse constraint. Results, The iMoCo UTE strategy was evaluated in volunteers and non-sedated pediatric patient(4-6 y/o) studies. Images reconstructed with iMoCo UTE provided sharper anatomical lung structures, and higher apparent SNR and CNR, compared to using other motion correction strategies, such as soft-gating, motion resolved reconstruction, and non-rigid motion compensation(MoCo). iMoCo UTE also showed promising results in an infant study. Conclusions, The proposed iMoCo UTE combines self-navigation, motion modeling, and a compressed sensing reconstruction to increase scan efficiency, SNR, and reduce respiratory motion in lung MRI. This proposed strategy shows improvements in free breathing lung MRI scans, especially in very challenging application situations, such as pediatric MRI studies.


title: Consensus-Based Optimization on Hypersurfaces: Well-Posedness and Mean-Field Limit

abstract: We introduce a new stochastic differential model for global optimization of nonconvex functions on compact hypersurfaces. The model is inspired by the stochastic Kuramoto-Vicsek system and belongs to the class of Consensus-Based Optimization methods. In fact, particles move on the hypersurface driven by a drift towards an instantaneous consensus point, computed as a convex combination of the particle locations weighted by the cost function according to Laplace's principle. The consensus point represents an approximation to a global minimizer. The dynamics is further perturbed by a random vector field to favor exploration, whose variance is a function of the distance of the particles to the consensus point. In particular, as soon as the consensus is reached, then the stochastic component vanishes. In this paper, we study the well-posedness of the model and we derive rigorously its mean-field approximation for large particle limit.


title: Persistence and periodic solutions in systems of delay differential equations

abstract: We study semi-dynamical systems associated to delay differential equations. We give a simple criteria to obtain weak and strong persistence and provide sufficient conditions to guarantee uniform persistence. Moreover, we show the existence of non-trivial $T$-periodic solutions via topological degree techniques. Finally, we prove that, in some sense, the conditions are also necessary.


title: PT-Symmetric Non-Hermitian Quantum Field Theories with Supersymmetry

abstract: We formulate supersymmetric non-Hermitian quantum field theories with PT symmetry, starting with free chiral boson/fermion models and then including trilinear superpotential interactions. We consider models with both Dirac and Majorana fermions, analyzing them in terms of superfields and at the component level. We also discuss the relation between the equations of motion, the (non-)invariance of the Lagrangian and the (non-)conservation of the supercurrents in the two models. We exhibit a similarity transformation that maps the free-field supersymmetric PT-symmetric Dirac model to a supersymmetric Hermitian theory, but there is, in general, no corresponding similarity transformation for the Majorana model. In this model, we find generically a mass splitting between bosons and fermions, even though its construction is explicitly supersymmetric, offering a novel non-Hermitian mechanism for soft supersymmetry breaking.


title: Ultimate greedy approximation of independent sets in subcubic graphs

abstract: We study the approximability of the maximum size independent set (MIS) problem in bounded degree graphs. This is one of the most classic and widely studied NP-hard optimization problems. We focus on the well known minimum degree greedy algorithm for this problem. This algorithm iteratively chooses a minimum degree vertex in the graph, adds it to the solution and removes its neighbors, until the remaining graph is empty. The approximation ratios of this algorithm have been very widely studied, where it is augmented with an advice that tells the greedy which minimum degree vertex to choose if it is not unique. Our main contribution is a new mathematical theory for the design of such greedy algorithms with efficiently computable advice and for the analysis of their approximation ratios. With this new theory we obtain the ultimate approximation ratio of 5/4 for greedy on graphs with maximum degree 3, which completely solves the open problem from the paper by Halldorsson and Yoshihara (1995). Our algorithm is the fastest currently known algorithm with this approximation ratio on such graphs. We apply our new algorithm to the minimum vertex cover problem on graphs with maximum degree 3 to obtain a substantially faster 6/5-approximation algorithm than the one currently known. We complement our positive, upper bound results with negative, lower bound results which prove that the problem of designing good advice for greedy is computationally hard and even hard to approximate on various classes of graphs. These results significantly improve on such previously known hardness results. Moreover, these results suggest that obtaining the upper bound results on the design and analysis of greedy advice is non-trivial.


title: Quasicrystalline and crystalline rare-gas clusters produced in supersonic jets: Impact of the jet clustering level on cathodoluminescence spectra

abstract: The paper proposes a new approach to studying cathodoluminescence spectra of substrate-free rare-gas clusters produced in supersonic jets exhausting into a vacuum. The approach, which takes into account the fraction of the clustered substance in the jet, is applied to quantitatively analyze integrated intensities of the luminescence bands of the neutral and charged excimer complexes (Rg2)* and (Rg4+)* measured for nanoclusters of three rare gases (Rg = Ar, Kr, and Xe) with average sizes ranging from 100 to 18000 atoms per clusters (diameters varying from 2 to 13 nm). The amount of the clustered substance, which affects the absolute values of integrated intensity of the bands, is shown to be proportional to the logarithm of the average size of clusters in the jet. Analysis of normalized intensities allowed us to spectroscopically find two ranges of average sizes of Ar, Kr, and Xe nanoclusters which, in accordance with the electron diffraction studies, can be assigned to quasicrystalline icosahedral and crystalline fcc structures in clusters, as well as to find the cluster size range in which both structures coexist. We show that in fcc clusters the luminescence of the neutral molecules (Rg2)* comes from within the volume of the cluster, while the charged excimer complexes (Rg4+)* emit mostly from subsurface layers.


title: Combining realization space models of polytopes

abstract: In this paper we examine four different models for the realization space of a polytope: the classical model, the Grassmannian model, the Gale transform model, and the slack variety. Respectively, they identify realizations of the polytopes with the matrix whose columns are the coordinates of their vertices, the column space of said matrix, their Gale transforms, and their slack matrices. Each model has been used to study realizations of polytopes. In this paper we establish very explicitly the maps that allow us to move between models, study their precise relationships, and combine the strengths of different viewpoints. As an illustration, we combine the compact nature of the Grassmannian model with the slack variety to obtain a reduced slack model that allows us to perform slack ideal calculations that were previously out of computational reach. These calculations allow us to answer the question of [Criado, Santos. Topological prismatoids and small simplicial spheres of large diameter. Experimental Mathematics, 1-13, 2019], about the realizability of a family of prismatoids, in general in the negative by proving the non-realizability of one of them.


title: Modeling synchronization in forced turbulent oscillator flows

abstract: Periodically forced, oscillatory fluid flows have been the focus of intense research for decades due to their richness as a nonlinear dynamical system and their relevance to applications in transportation, aeronautics, and energy conversion. Recently, it has been observed that turbulent bluff-body wakes exhibit a subharmonic resonant response when excited with specific spatial symmetries at twice the natural vortex shedding frequency, which is hypothesized to be caused by triadic interactions. The focus of this paper is to provide new physical insight into the dynamics of turbulent oscillator flows, based on improved mechanistic models informed by a comprehensive experimental study of the turbulent wake behind a D-shaped body under periodic forcing. We confirm for the first time the role of resonant triadic interactions in the forced flow by studying the dominant components in the power spectra across multiple excitation frequencies and amplitudes. We then develop an extended Stuart-Landau model for the forced global wake mode, incorporating parametric and non-harmonic forcing. This model captures the system dynamics and reveals the boundaries of multiple synchronization regions. Further, it is possible to identify model coefficients from sparse measurement data, making it applicable to a wide range of turbulent oscillator flows. We believe these generalized synchronization models will be valuable for prediction, control, and understanding of the underlying physics in this ubiquitous class of flows.


title: Compositions with 3 Pairwise Coprime Parts

abstract: How many ways can we write $n$ as a sum of $3$ positive integers, no pair of which share a common factor? We express this quantity in terms of the number of solutions to a certain class of linear Diophantine equations. This allows us to show that there are $$ \prod_{p \mid n} \left( 1- \frac{1}{p^2} \right) \prod_{q \nmid n} \left( 1- \frac{3}{q^2} \right) \frac{n^2}{2} + O(n^{3/2+o(1)}) $$ such compositions, where the products are over primes that respectively do and don't divide $n$. This strengthens the previous result of Bubbolini, Luca, and Spiga (arXiv:1202.1670)


title: Efficient Langevin dynamics for "noisy" forces

abstract: Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Stochastic approaches provide a gentler system-size dependency at the cost of introducing "noisy" forces, which serve to limit the efficiency of the sampling. In the first-order Langevin dynamics (FOLD), efficient sampling is achievable by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator (Mazzola et al., Phys. Rev. Lett., 118, 015703 (2017)). However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure which severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.


title: Photon-Photon Scattering in Very Special Relativity

abstract: Starting from the Mandelstam-Leibbrandt prescription, we introduce a general rule for the null vector $\bar{n}$ to compute any $SIM(2)$ integral and diagram with an arbitrary number of external legs. Using the new prescription, the computation of the low energy limit of photon-photon scattering under the Very Special Relativity (VSR) framework is presented. The prescription preserves automatically the Ward identities corresponding to the gauge symmetry. Within the low momentum approximation we get the standard unpolarized differential cross section for photon-photon scattering. The result suggests that loops with any external photon legs on-shell will be zero in VSR.


title: Neural MMO v1.3: A Massively Multiagent Game Environment for Training and Evaluating Neural Networks

abstract: Progress in multiagent intelligence research is fundamentally limited by the number and quality of environments available for study. In recent years, simulated games have become a dominant research platform within reinforcement learning, in part due to their accessibility and interpretability. Previous works have targeted and demonstrated success on arcade, first person shooter (FPS), real-time strategy (RTS), and massive online battle arena (MOBA) games. Our work considers massively multiplayer online role-playing games (MMORPGs or MMOs), which capture several complexities of real-world learning that are not well modeled by any other game genre. We present Neural MMO, a massively multiagent game environment inspired by MMOs and discuss our progress on two more general challenges in multiagent systems engineering for AI research: distributed infrastructure and game IO. We further demonstrate that standard policy gradient methods and simple baseline models can learn interesting emergent exploration and specialization behaviors in this setting.


title: Reflection principles, GCH and the uniformization properties

abstract: Reflection principles (or dually speaking, compactness principles) often give rise to combinatorial guessing principles. Uniformization properties, on the other hand, are examples of anti-guessing principles. We discuss the tension and the compatibility between reflection principles and uniformization properties at the level of the second uncountable cardinal.


title: Theory inspired deep network for instantaneous-frequency extraction and signal components recovery from discrete blind-source data

abstract: This paper is concerned with the inverse problem of recovering the unknown signal components, along with extraction of their instantaneous frequencies (IFs), governed by the adaptive harmonic model (AHM), from discrete (and possibly non-uniform) samples of the blind-source composite signal. None of the existing decomposition methods and algorithms, including the most popular empirical mode decomposition (EMD) computational scheme and its current modifications, is capable of solving this inverse problem. In order to meet the AHM formulation and to extract the IFs of the decomposed components, called intrinsic mode functions (IMFs), each IMF of EMD is extended to an analytic function in the upper half of the complex plane via the Hilbert transform, followed by taking the real part of the polar form of the analytic extension. Unfortunately, this approach most often fails to resolve the inverse problem satisfactorily. More recently, to resolve the inverse problem, the notion of synchrosqueezed wavelet transform (SST) was proposed by Daubechies and Maes, and further developed in many other papers, while a more direct method, called signal separation operation (SSO), was proposed and developed in our previous work published in the journal, Applied and Computational Harmonic Analysis, vol. 30(2):243-261, 2016. In the present paper, we propose a synthesis of SSO using a deep neural network, based directly on a discrete sample set, that may be non-uniformly sampled, of the blind-source signal. Our method is localized, as illustrated by a number of numerical examples, including components with different signal arrival and departure times. It also yields short-term prediction of the signal components, along with their IFs. Our neural networks are inspired by theory, designed so that they do not require any training in the traditional sense.


title: Tracing the anemic stellar halo of M101

abstract: Models of galaxy formation in a cosmological context predict that massive disk galaxies should have structured extended stellar halos. Recent studies in integrated light, however, report a few galaxies, including the nearby disk galaxy M101, that have no measurable stellar halos to the detection limit. We aim to quantify the stellar content and structure of M101's outskirts by resolving its stars. We present the photometry of its stars based on deep F606W and F814W images taken with Hubble Space Telescope as part of the GHOSTS survey. The constructed CMDs of stars reach down to two magnitudes below the tip of the red giant branch. We derived radial number density profiles of the bright red giant branch (RGB) stars. The mean color of the RGB stars at $R \sim$ 40 -- 60 kpc is similar to those of metal-poor globular clusters in the Milky Way. We also derived radial surface brightness profiles using the public image data provided by the Dragonfly team. Both the radial number density and surface brightness profiles were converted to radial mass density profiles and combined. We find that the mass density profiles show a weak upturn at the very outer region, where surface brightness is as faint as $\mu_g\approx 34$ mag arcsec$^{-1}$. An exponential disk + power-law halo model on the mass density profiles finds the total stellar halo mass of $M_{halo}=8.2_{-2.2}^{+3.5}\times 10^7M_\odot$. The total stellar halo mass does not exceed $M_{halo} = 3.2 \times 10^8$ $M_{\odot}$ when strongly truncated disk models are considered. Combining the halo mass with the total stellar mass of M101, we obtain the stellar halo mass fraction of $M_{halo}/M_{gal} = 0.20_{-0.08}^{+0.10}\%$ with an upper limit of 0.78\%. We compare the halo properties of M101 with those of six GHOSTS survey galaxies as well as the Milky Way and M31 and find that M101 has an anemic stellar halo.


title: Conformal Skorokhod embeddings and related extremal problems

abstract: The conformal Skorokhod embedding problem (CSEP) is a planar variant of the classical problem where the solution is now a simply connected domain $D\subset\mathbb{C}$ whose exit time embeds a given probability distribution $\mu$ by projecting the stopped Brownian motion onto the real axis. In this paper we explore two new research directions for the CSEP by proving general bounds on the principal Dirichlet eigenvalue of a solution domain in terms of the corresponding $\mu$ and by proposing related extremal problems. Moreover, we give a new and nontrivial example of an extremal domain $\mathbb{U}$ that attains the lowest possible principal Dirichlet eigenvalue over all domains solving the CSEP for the uniform distribution on $[-1,1]$. Remarkably, the boundary of $\mathbb{U}$ is related to the Grim Reaper translating solution to the curve shortening flow in the plane. The novel tool used in the proof of the sharp lower bound is a precise relationship between the widths of the orthogonal projections of a simply connected planar domain and the support of its harmonic measure that we develop in the paper. The upper bound relies on spectral bounds for the torsion function which have recently appeared in the literature.


title: Distance mapping applied to four well-known planetary nebulae and a nova shell

abstract: Accurate distance estimates of astrophysical objects such as planetary nebulae (PNe), and nova and supernova remnants, among others, allow us to constrain their physical characteristics, such as size, mass, luminosity, and age. An innovative technique based on the expansion parallax method, the so-called distance mapping technique (DMT), provides distance maps of expanding nebulae as well as an estimation of their distances. The DMT combines the tangential velocity vectors obtained from 3D morpho-kinematic models and the observed proper motion vectors to estimate the distance. We applied the DMT to four PNe (NGC 6702, NGC 6543, NGC 6302, and BD+30 3639) and one nova remnant (GK Persei) and derived new distances in good agreement with previous studies. New simple morpho-kinematic shape models were generated for NGC 6543, NGC 6302, and NGC 6702, whereas for BD+30 3639 and GK Persei published models were used. We demonstrate that the DMT is a useful tool to obtain distance values of PNe, in addition to revealing kinematically peculiar regions within the nebulae. Distances are also derived from the trigonometric Gaia parallaxes. The effect of the non-negligible parallax offset in the second Gaia data release is also discussed.


title: Learning Deep Analysis Dictionaries for Image Super-Resolution

abstract: Inspired by the recent success of deep neural networks and the recent efforts to develop multi-layer dictionary models, we propose a Deep Analysis dictionary Model (DeepAM) which is optimized to address a specific regression task known as single image super-resolution. Contrary to other multi-layer dictionary models, our architecture contains L layers of analysis dictionary and soft-thresholding operators to gradually extract high-level features and a layer of synthesis dictionary which is designed to optimize the regression task at hand. In our approach, each analysis dictionary is partitioned into two sub-dictionaries: an Information Preserving Analysis Dictionary (IPAD) and a Clustering Analysis Dictionary (CAD). The IPAD together with the corresponding soft-thresholds is designed to pass the key information from the previous layer to the next layer, while the CAD together with the corresponding soft-thresholding operator is designed to produce a sparse feature representation of its input data that facilitates discrimination of key features. DeepAM uses both supervised and unsupervised setup. Simulation results show that the proposed deep analysis dictionary model achieves better performance compared to a deep neural network that has the same structure and is optimized using back-propagation when training datasets are small.


title: The Ballet of Triangle Centers on the Elliptic Billiard

abstract: The dynamic geometry of the family of 3-periodics in the Elliptic Billiard is mystifying. Besides conserving perimeter and the ratio of inradius-to-circumradius, it has a stationary point. Furthermore, its triangle centers sweep out mesmerizing loci including ellipses, quartics, circles, and a slew of other more complex curves. Here we explore a bevy of new phenomena relating to (i) the shape of 3-periodics and (ii) the kinematics of certain Triangle Centers constrained to the Billiard boundary, specifically the non-monotonic motion some can display with respect to 3-periodics. Hypnotizing is the joint motion of two such non-monotonic Centers, whose many stops-and-gos are akin to a Ballet.


title: Algorithmic Aspects of Some Variants of Domination in Graphs

abstract: A set $S \subseteq V$ is a dominating set in G if for every u \in V \ S, there exists $v \in S$ such that $(u,v) \in E$, i.e., $N[S] = V$. A dominating set $S$ is an Isolate Dominating Set} (IDS) if the induced subgraph $G[S]$ has at least one isolated vertex. It is known that Isolate Domination Decision problem (IDOM) is NP-complete for bipartite graphs. In this paper, we extend this by showing that the IDOM is NP-complete for split graphs and perfect elimination bipartite graphs, a subclass of bipartite graphs. A set $S \subseteq V$ is an independent set if G[S] has no edge. A set S \subseteq V is a secure dominating set of $G$ if, for each vertex $u \in V \setminus S$, there exists a vertex $v \in S$ such that $ (u,v) \in E $ and $(S \ \{v\}) \cup \{u\}$ is a dominating set of $G$. In addition, we initiate the study of a new domination parameter called, independent secure domination. A set $S\subseteq V$ is an Independent Secure Dominating Set (InSDS) if $S$ is an independent set and a secure dominating set of $G$. The minimum size of an InSDS in $G$ is called the independent secure domination number of $G$ and is denoted by $\gamma_{is}(G)$. Given a graph $ G$ and a positive integer $ k,$ the InSDM problem is to check whether $ G $ has an independent secure dominating set of size at most $ k.$ We prove that InSDM is NP-complete for bipartite graphs and linear time solvable for bounded tree-width graphs and threshold graphs, a subclass of split graphs. The MInSDS problem is to find an independent secure dominating set of minimum size, in the input graph. Finally, we prove that the MInSDS problem is APX-hard for graphs with maximum degree $5.$


title: Temperature dependence of spin pinning and spin-wave dispersion in nanoscopic ferromagnetic waveguides

abstract: The field of magnonics attracts significant attention due to the possibility of utilizing information coded into the spin-wave phase or amplitude to perform computation operations on the nanoscale. Recently, spin waves were investigated in Yttrium Iron Garnet (YIG) waveguides with widths ranging down to 50 nm and aspect ratios thickness over width approaching unity. A critical width was found, below which the exchange interaction suppresses the dipolar pinning phenomenon and the system becomes unpinned. Here we continue these investigations and analyse the pinning phenomenon and spin-wave dispersions as a function of temperature, thickness and material of choice. Higher order modes, the influence of a finite wavevector along the waveguide and the impact of the pinning phenomenon on the spin-wave lifetime are discussed as well as the influence of a trapezoidal cross section and edge roughness of the waveguides. The presented results are of particular interest for potential applications in magnonic devices and the incipient field of quantum magnonics at cryogenic temperatures.


title: Optimal upper bound of entropic uncertainty relation for mutually unbiased bases

abstract: We have obtained the optimal upper bound of entropic uncertainty relation for $N$ Mutually Unbiased Bases (MUBs). We have used the methods of variational calculus for the states that can be written in terms of $N$ MUBs. Our result is valid for any state when $N$ is $d+1$, where $d$ is the dimension of the related system. We provide a quantitative criterion for the extendibilty of MUBs. In addition, we have applied our result to the mutual information of $d+1$ observables conditioned with a classical memory.


title: Investigating elastic constants across diverse strain-matrix sets

abstract: Elastic constants and mechanical properties play a pivotal role across multiple disciplines and engineering applications. We introduced the optimized high-efficient strain-matrix set (OHESS) that determines the second-order elastic constants of materials using the stress-strain method. Herein, we systematically investigate the computational efficiency of OHESS across a broad range of crystal systems and compare it with other notable stress-strain approaches, such as the single-element strain-matrix sets and the universal linear-independent coupling strains. Notably, our data affirm the superior efficacy of OHESS among the strain-matrix sets under consideration. We believe OHESS will markedly improve computational efficiency in determining the elastic constants and mechanical properties, becoming an indispensable tool for material research, design, and high-throughput screening.


title: Nonuniform sampling and approximation in Sobolev space from the perturbation of framelet system

abstract: The Sobolev space $H^{\varsigma}(\mathbb{R}^{d})$, where $\varsigma > d/2$, is an important function space that has many applications in various areas of research. Attributed to the inertia of a measurement instrument, it is desirable in sampling theory to recover a function by its nonuniform sampling. In the present paper, based on dual framelet systems for the Sobolev space pair $(H^{s}(\mathbb{R}^{d}), H^{-s}(\mathbb{R}^{d}))$, where $d/2<s<\varsigma$, we investigate the problem of constructing the approximations to all the functions in $H^{\varsigma}(\mathbb{R}^{d})$ by nonuniform sampling. We first establish the convergence rate of the framelet series in $(H^{s}(\mathbb{R}^{d}), H^{-s}(\mathbb{R}^{d}))$, and then construct the framelet approximation operator that acts on the entire space $H^{\varsigma}(\mathbb{R}^{d})$. We examine the stability property for the framelet approximation operator with respect to the perturbations of shift parameters, and obtain an estimate bound for the perturbation error. Our result shows that under the condition $d/2<s<\varsigma$, the approximation operator is robust to shift perturbations. Motivated by some recent work on nonuniform sampling and approximation in Sobolev space (e.g., [20]), we don't require the perturbation sequence to be in $\ell^{\alpha}(\mathbb{Z}^{d})$. Our results allow us to establish the approximation for every function in $H^{\varsigma}(\mathbb{R}^{d})$ by nonuniform sampling. In particular, the approximation error is robust to the jittering of the samples.


title: Ultra-Discretization of $D_6^{(1)}$- Geometric Crystal at the spin node

abstract: Let $\mathfrak g$ be an affine Lie algebra with index set $I = \{0, 1, 2, \cdots , n\}$. It is conjectured in \cite{KNO} that for each Dynkin node $k \in I \setminus \{0\}$ the affine Lie algebra $\mathfrak g$ has a positive geometric crystal whose ultra-discretization is isomorphic to the limit of a coherent family of perfect crystals for the Langland dual ${\mathfrak g} ^L$. In this paper we show that at the spin node $k=6$, the family of perfect crystals given in \cite{KMN2} form a coherent family and show that its limit $B^{6,\infty}$ is isomorphic to the ultra-discretization of the positive geometric crystal we constructed in \cite{MP} for the affine Lie algebra $D_6^{(1)}$ which proves the conjecture in this case.


title: On the Information Bottleneck Problems: Models, Connections, Applications and Information Theoretic Views

abstract: This tutorial paper focuses on the variants of the bottleneck problem taking an information theoretic perspective and discusses practical methods to solve it, as well as its connection to coding and learning aspects. The intimate connections of this setting to remote source-coding under logarithmic loss distortion measure, information combining, common reconstruction, the Wyner-Ahlswede-Korner problem, the efficiency of investment information, as well as, generalization, variational inference, representation learning, autoencoders, and others are highlighted. We discuss its extension to the distributed information bottleneck problem with emphasis on the Gaussian model and highlight the basic connections to the uplink Cloud Radio Access Networks (CRAN) with oblivious processing. For this model, the optimal trade-offs between relevance (i.e., information) and complexity (i.e., rates) in the discrete and vector Gaussian frameworks is determined. In the concluding outlook, some interesting problems are mentioned such as the characterization of the optimal inputs ("features") distributions under power limitations maximizing the "relevance" for the Gaussian information bottleneck, under "complexity" constraints.


title: Probabilistic Complexity Classes through Semantics

abstract: In a recent paper, the author has shown how Interaction Graphs models for linear logic can be used to obtain implicit characterisations of non-deterministic complexity classes. In this paper, we show how this semantic approach to Implicit Complexity Theory (ICC) can be used to characterise deterministic and probabilistic models of computation. In doing so, we obtain correspondences between group actions and both deterministic and probabilistic hierarchies of complexity classes. As a particular case, we provide the first implicit characterisations of the classes PLogspace (un-bounded error probabilistic logarithmic space) and PPtime (unbounded error probabilistic polynomial time)


title: Classification of Computer Models with Labelled Outputs

abstract: Classification is a vital tool that is important for modelling many complex numerical models. A model or system may be such that, for certain areas of input space, the output either does not exist, or is not in a quantifiable form. Here, we present a new method for classification where the model outputs are given distinct classifying labels, which we model using a latent Gaussian process (GP). The latent variable is estimated using MCMC sampling, a unique likelihood and distinct prior specifications. Our classifier is then verified by calculating a misclassification rate across the input space. Comparisons are made with other existing classification methods including logistic regression, which models the probability of being classified into one of two regions. To make classification predictions we draw from an independent Bernoulli distribution, meaning that distance correlation is lost from the independent draws and so can result in many misclassifications. By modelling the labels using a latent GP, this problem does not occur in our method. We apply our novel method to a range of examples including a motivating example which models the hormones associated with the reproductive system in mammals, where the two labelled outputs are high and low rates of reproduction.


title: Age-Conditioned Synthesis of Pediatric Computed Tomography with Auxiliary Classifier Generative Adversarial Networks

abstract: Deep learning is a popular and powerful tool in computed tomography (CT) image processing such as organ segmentation, but its requirement of large training datasets remains a challenge. Even though there is a large anatomical variability for children during their growth, the training datasets for pediatric CT scans are especially hard to obtain due to risks of radiation to children. In this paper, we propose a method to conditionally synthesize realistic pediatric CT images using a new auxiliary classifier generative adversarial network (ACGAN) architecture by taking age information into account. The proposed network generated age-conditioned high-resolution CT images to enrich pediatric training datasets.


title: The Birth of a Massive First-Star Binary

abstract: We study the formation of massive Population III binary stars using a newly developed radiation hydrodynamics code with the adaptive mesh refinement and adaptive ray-tracing methods. We follow the evolution of a typical primordial star-forming cloud obtained from a cosmological hydrodynamics simulation. Several protostars form as a result of disk fragmentation and grow in mass by the gas accretion, which is finally quenched by the radiation feedback from the protostars. Our code enables us, for the first time, to consider the feedback by both the ionizing and dissociating radiation from the multiple protostars, which is essential for self-consistently determining their final masses. At the final step of the simulation, we observe a very wide ($\gtrsim 10^4\,\mathrm{au}$) binary stellar system consisting of $60$ and $70\,M_\odot$ stars. One of the member stars also has two smaller mass ($10\,M_\odot$) companion stars orbiting at $200$ and $800\,\mathrm{au}$, making up a mini-triplet system. Our results suggest that massive binary or multiple systems are common among Population III stars.


title: Anisotropy-mediated reentrant localization

abstract: We consider a 2d dipolar system, $d=2$, with the generalized dipole-dipole interaction $\sim r^{-a}$, and the power $a$ controlled experimentally in trapped-ion or Rydberg-atom systems via their interaction with cavity modes. We focus on the dilute dipolar excitation case when the problem can be effectively considered as single-particle with the interaction providing long-range dipolar-like hopping. We show that the spatially homogeneous tilt $\beta$ of the dipoles giving rise to the anisotropic dipole exchange leads to the non-trivial reentrant localization beyond the locator expansion, $a<d$, unlike the models with random dipole orientation. The Anderson transitions are found to occur at the finite values of the tilt parameter $\beta = a$, $0<a<d$, and $\beta = a/(a-d/2)$, $d/2<a<d$, showing the robustness of the localization at small and large anisotropy values. Both extensive numerical calculations and analytical methods show power-law localized eigenstates in the bulk of the spectrum, obeying recently discovered duality $a\leftrightarrow 2d-a$ of their spatial decay rate, on the localized side of the transition, $a>a_{AT}$. This localization emerges due to the presence of the ergodic extended states at either spectral edge, which constitute a zero fraction of states in the thermodynamic limit, decaying though extremely slowly with the system size.


title: Close Binary Companions to APOGEE DR16 Stars: 20,000 Binary-star Systems Across the Color-Magnitude Diagram

abstract: Many problems in contemporary astrophysics---from understanding the formation of black holes to untangling the chemical evolution of galaxies---rely on knowledge about binary stars. This, in turn, depends on discovery and characterization of binary companions for large numbers of different kinds of stars in different chemical and dynamical environments. Current stellar spectroscopic surveys observe hundreds of thousands to millions of stars with (typically) few observational epochs, which allows binary discovery but makes orbital characterization challenging. We use a custom Monte Carlo sampler (The Joker) to perform discovery and characterization of binary systems through radial-velocities, in the regime of sparse, noisy, and poorly sampled multi-epoch data. We use it to generate posterior samplings in Keplerian parameters for 232,531 sources released in APOGEE Data Release 16. Our final catalog contains 19,635 high-confidence close-binary (P < few years, a < few AU) systems that show interesting relationships between binary occurrence rate and location in the color-magnitude diagram. We find notable faint companions at high masses (black-hole candidates), at low masses (substellar candidates), and at very close separations (mass-transfer candidates). We also use the posterior samplings in a (toy) hierarchical inference to measure the long-period binary-star eccentricity distribution. We release the full set of posterior samplings for the entire parent sample of 232,531 stars. This set of samplings involves no heuristic "discovery" threshold and therefore can be used for myriad statistical purposes, including hierarchical inferences about binary-star populations and sub-threshold searches.


title: More accurate simulations with separate initial conditions for baryons and dark matter

abstract: We revisit techniques for performing cosmological simulations with both baryons and cold dark matter when each fluid has different initial conditions, as is the case at the end of the radiation era. Most simulations do not reproduce the linear prediction for the difference between the cold dark matter and baryon perturbations. We show that this is due to the common use of offset regular grids when setting up the particle initial conditions. The correct behaviour can be obtained without any loss of simulation resolution by using a Lagrangian glass for the baryon particles. We further show that the difference between cold dark matter and baryons may affect predictions for the Lyman-alpha forest flux power spectrum at the 5% level, potentially impacting current cosmological constraints.


title: A proposal for reconciling diverse experiments on the superconducting state in Sr2RuO4

abstract: A variety of precise experiments have been carried out to establish the character of the superconducting state in Sr2RuO4. Many of these appear to imply contradictory conclusions concerning the symmetries of this state. Here, we propose that these results can be reconciled if we assume that there is a near-degeneracy between a d_{x^2-y^2} (B_{1g} in group theory nomenclature) and a g_{xy(x^2-y^2)} (A_{2g}) superconducting state. From a weak-coupling perspective, such an accidental degeneracy can occur at a point at which a balance between the on-site and nearest-neighbor repulsions triggers a d-wave to g-wave transition.


title: Functoriality of HKR isomorphisms

abstract: For a closed embedding of smooth schemes $X\hookrightarrow S$ with a fixed first order splitting, one can construct HKR isomorphisms between the derived scheme $X\times^R_S X$ and the total space of the shifted normal bundle $\mathbb{N}_{X/S}[-1]$, due to Arinkin-C\u{a}ld\u{a}raru, Arinkin-C\u{a}ld\u{a}raru-Hablicsek, and Grivaux. In this paper, we study functoriality property of the HKR isomorphisms for a sequence of closed embeddings $X\hookrightarrow Y\hookrightarrow S$. The HKR isomorphism is functorial when a certain cohomology class, which we call the Bass-Quillen class, vanishes. We obtain Lie theoretic interpretations for the HKR isomorphisms and for the Bass-Quillen class as well.


title: Modular Hamiltonians for Euclidean Path Integral States

abstract: We study half-space/Rindler modular Hamiltonians for excited states created by turning on sources for local operators in the Euclidean path integral in relativistic quantum field theories. We derive a simple, manifestly Lorentzian formula for the modular Hamiltonian to all orders in perturbation theory in the sources. We apply this formula to the case of shape-deformed half spaces in the vacuum state, and obtain the corresponding modular Hamiltonian to all orders in the shape deformation in terms of products of half-sided null energy operators, i.e., stress tensor components integrated along the future and past Rindler horizons. In the special case where the shape deformation is purely null, our perturbation series can be resummed, and agrees precisely with the known formula for vacuum modular Hamiltonians for null cuts of the Rindler horizon. Finally, we study some universal properties of modular flow (corresponding to Euclidean path integral states) of local operators inside correlation functions in conformal field theories. In particular, we show how the flow becomes the local boost in the limit where the operator being flowed approaches the entanglement cut.


title: A large deviation principle linking lineage statistics to fitness in microbial populations

abstract: In exponentially proliferating populations of microbes, the population typically doubles at a rate less than the average doubling time of a single-cell due to variability at the single-cell level. It is known that the distribution of generation times obtained from a single lineage is, in general, insufficient to determine a population's growth rate. Is there an explicit relationship between observables obtained from a single lineage and the population growth rate? We show that a population's growth rate can be represented in terms of averages over isolated lineages. This lineage representation is related to a large deviation principle that is a generic feature of exponentially proliferating populations. Due to the large deviation structure of growing populations, the number of lineages needed to obtain an accurate estimate of the growth rate depends exponentially on the duration of the lineages, leading to a non-monotonic convergence of the estimate, which we verify in both synthetic and experimental data sets.


title: Das ist der HAMMER: Consistent new physics interpretations of semileptonic decays

abstract: Precise measurements of $b\to c\tau\bar\nu$ decays require large resource-intensive Monte Carlo (MC) samples, which incorporate detailed simulations of detector responses and physics backgrounds. Extracted parameters may be highly sensitive to the underlying theoretical models used in the MC generation. Because new physics (NP) can alter decay distributions and acceptances, the standard practice of fitting NP Wilson coefficients to SM-based measurements of the $R(D^{(*)})$ ratios can be biased. The newly developed HAMMER software tool enables efficient reweighting of MC samples to arbitrary NP scenarios or to any hadronic matrix elements. We demonstrate how HAMMER allows avoidance of biases through self-consistent fits directly to the NP Wilson coefficients. We also present example analyses that demonstrate the sizeable biases that can otherwise occur from naive NP interpretations of SM-based measurements. The HAMMER library is presently interfaced with several existing experimental analysis frameworks and we provide an overview of its structure.


title: Embedding Hard Physical Constraints in Neural Network Coarse-Graining of 3D Turbulence

abstract: In the recent years, deep learning approaches have shown much promise in modeling complex systems in the physical sciences. A major challenge in deep learning of PDEs is enforcing physical constraints and boundary conditions. In this work, we propose a general framework to directly embed the notion of an incompressible fluid into Convolutional Neural Networks, and apply this to coarse-graining of turbulent flow. These physics-embedded neural networks leverage interpretable strategies from numerical methods and computational fluid dynamics to enforce physical laws and boundary conditions by taking advantage the mathematical properties of the underlying equations. We demonstrate results on three-dimensional fully-developed turbulence, showing that this technique drastically improves local conservation of mass, without sacrificing performance according to several other metrics characterizing the fluid flow.


title: Learning Deep Analysis Dictionaries -- Part II: Convolutional Dictionaries

abstract: In this paper, we introduce a Deep Convolutional Analysis Dictionary Model (DeepCAM) by learning convolutional dictionaries instead of unstructured dictionaries as in the case of deep analysis dictionary model introduced in the companion paper. Convolutional dictionaries are more suitable for processing high-dimensional signals like for example images and have only a small number of free parameters. By exploiting the properties of a convolutional dictionary, we present an efficient convolutional analysis dictionary learning approach. A L-layer DeepCAM consists of L layers of convolutional analysis dictionary and element-wise soft-thresholding pairs and a single layer of convolutional synthesis dictionary. Similar to DeepAM, each convolutional analysis dictionary is composed of a convolutional Information Preserving Analysis Dictionary (IPAD) and a convolutional Clustering Analysis Dictionary (CAD). The IPAD and the CAD are learned using variations of the proposed learning algorithm. We demonstrate that DeepCAM is an effective multilayer convolutional model and, on single image super-resolution, achieves performance comparable with other methods while also showing good generalization capabilities.


title: Consistency checks for two-body finite-volume matrix elements: II. Perturbative systems

abstract: Using the general formalism presented in Refs. [1,2], we study the finite-volume effects for the $\mathbf{2}+\mathcal{J}\to\mathbf{2}$ matrix element of an external current coupled to a two-particle state of identical scalars with perturbative interactions. Working in a finite cubic volume with periodicity $L$, we derive a $1/L$ expansion of the matrix element through $\mathcal O(1/L^5)$ and find that it is governed by two universal current-dependent parameters, the scalar charge and the threshold two-particle form factor. We confirm the result through a numerical study of the general formalism and additionally through an independent perturbative calculation. We further demonstrate a consistency with the Feynman-Hellmann theorem, which can be used to relate the $1/L$ expansions of the ground-state energy and matrix element. The latter gives a simple insight into why the leading volume corrections to the matrix element have the same scaling as those in the energy, $1/L^3$, in contradiction to earlier work, which found a $1/L^2$ contribution to the matrix element. We show here that such a term arises at intermediate stages in the perturbative calculation, but cancels in the final result.


title: Limit theorems of stochastic differential equations with jumps

abstract: In this paper, we study the convergence for solutions to a sequence of (possibly degenerate) stochastic differential equations with jumps, when the coefficients converge in some appropriate sense. Our main tools are the superposition principles. And then we analyze some special cases and give some concrete and verifiable conditions.


title: Gating creates slow modes and controls phase-space complexity in GRUs and LSTMs

abstract: Recurrent neural networks (RNNs) are powerful dynamical models for data with complex temporal structure. However, training RNNs has traditionally proved challenging due to exploding or vanishing of gradients. RNN models such as LSTMs and GRUs (and their variants) significantly mitigate these issues associated with training by introducing various types of gating units into the architecture. While these gates empirically improve performance, how the addition of gates influences the dynamics and trainability of GRUs and LSTMs is not well understood. Here, we take the perspective of studying randomly initialized LSTMs and GRUs as dynamical systems, and ask how the salient dynamical properties are shaped by the gates. We leverage tools from random matrix theory and mean-field theory to study the state-to-state Jacobians of GRUs and LSTMs. We show that the update gate in the GRU and the forget gate in the LSTM can lead to an accumulation of slow modes in the dynamics. Moreover, the GRU update gate can poise the system at a marginally stable point. The reset gate in the GRU and the output and input gates in the LSTM control the spectral radius of the Jacobian, and the GRU reset gate also modulates the complexity of the landscape of fixed-points. Furthermore, for the GRU we obtain a phase diagram describing the statistical properties of fixed-points. We also provide a preliminary comparison of training performance to the various dynamical regimes realized by varying hyperparameters. Looking to the future, we have introduced a powerful set of techniques which can be adapted to a broad class of RNNs, to study the influence of various architectural choices on dynamics, and potentially motivate the principled discovery of novel architectures.

categories: cs.LG cond-mat.dis-nn cond-mat.stat-mech stat.ML


title: Magnetic properties of rare earth and transition metal based perovskite type high entropy oxides

abstract: High entropy oxides (HEO) are a recently introduced class of oxide materials, which are characterized by a large number of elements (i.e. five or more) sharing one lattice site which crystallize in a single phase structure. One complex example of the rather young HEO family are the rare-earth transition metal perovskite high entropy oxides. In this comprehensive study, we provide an overview over the magnetic properties of three perovskite type high entropy oxides. The compounds have a rare-earth site which is occupied by five different rare-earth elements, while the transition metal site is occupied by a single transition metal. In this way a comparison to the parent binary oxides, namely the orthocobaltites, -chromites and -ferrites is possible. X-ray absorption near edge spectroscopy (XANES), magnetometry and M\"ossbauer spectroscopy are employed to characterize these complex materials. In general, we find surprising similarities to the magnetic properties of the binary oxides, despite the chemical disorder on the rare-earth site. However distinct differences and interesting magnetic properties are also observed such as noncollinearity, spin reorientation transitions as well as large coercive fields of up to 2\,T at ambient temperature. Both the chemical disorder on the RE A-site, and the nature of the TM on the B-site play an important role in the physical properties of these high entropy oxides.


title: Hypercomplex-Valued Recurrent Correlation Neural Networks

abstract: Recurrent correlation neural networks (RCNNs), introduced by Chiueh and Goodman as an improved version of the bipolar correlation-based Hopfield neural network, can be used to implement high-capacity associative memories. In this paper, we extend the bipolar RCNNs for processing hypercomplex-valued data. Precisely, we present the mathematical background for a broad class of hypercomplex-valued RCNNs. Then, we provide the necessary conditions which ensure that a hypercomplex-valued RCNN always settles at an equilibrium using either synchronous or asynchronous update modes. Examples with bipolar, complex, hyperbolic, quaternion, and octonion-valued RCNNs are given to illustrate the theoretical results. Finally, computational experiments confirm the potential application of hypercomplex-valued RCNNs as associative memories designed for the storage and recall of gray-scale images.


title: Valence instability across magnetostructural transition in USb$_2$

abstract: We have performed pressure dependent X-ray diffraction and resonant X-ray emission spectroscopy experiments on USb$_2$ to further characterize the AFM-FM transition occurring near 8 GPa. We have found the magnetic transition coincides with a tetragonal to orthorhombic transition resulting in a 17% volume collapse as well as a transient $\textit{f}$-occupation enhancement. Compared to UAs$_2$ and UAsS, USb$_2$ shows a reduced bulk modulus and transition pressure and an increased volume collapse at the structural transition. Except for an enhancement across the transition region, the $\textit{f}$-occupancy decreases steadily from 1.96 to 1.75.


title: Multiple string LED driver with flexible and high performance PWM dimming control

abstract: The main objectives in driving multiple LED strings include achieving uniform current control and high performance PWM dimming for all strings. This work proposes a new multiple string LED driver to achieve not only current balance, but also flexible and wide range PWM dimming ratio for each string. A compact single-inductor multiple-output topology is adopted in the driver, accompanied by synchronous integrators and variable dimming frequency, to achieve both high efficiency and high performance dimming. By using the proposed variable dimming frequency scheme, high dimming frequency is applied to a string with high dimming ratio, which helps to maintain the deviation of LED string current in an acceptable range, while low dimming frequency is applied to a string with low dimming ratio, which helps to achieve rectangular LED current waveform. Meanwhile, the new time multiplexing control scheme automatically optimizes the LED string's bus voltages, thus minimizes each string's power loss. A three string LED driver prototype is constructed to validate the effectiveness of the peoposed control scheme, where the three strings can have different dimming ratios between 4% and 100%.


title: Non-linear extension of interval arithmetic and exact resolution of interval equations over square regions

abstract: The interval numbers is the set of compact intervals of $\mathbb{R}$ with addition and multiplication operation, which are very useful for solving calculations where there are intervals of error or uncertainty, however, it lacks an algebraic structure with an inverse element, both additive and multiplicative This fundamental disadvantage results in overestimation of solutions in an interval equation or also overestimation of the image of a function over square regions. In this article we will present a solution to this problem, through a morphism that preserves both the addiction and the multiplication between the space of the interval numbers to the space of square diagonal matrices.


title: Optimal Transport Based Seismic Inversion: Beyond Cycle Skipping

abstract: Full-waveform inversion (FWI) is today a standard process for the inverse problem of seismic imaging. PDE-constrained optimization is used to determine unknown parameters in a wave equation that represent geophysical properties. The objective function measures the misfit between the observed data and the calculated synthetic data, and it has traditionally been the least-squares norm. In a sequence of papers, we introduced the Wasserstein metric from optimal transport as an alternative misfit function for mitigating the so-called cycle skipping, which is the trapping of the optimization process in local minima. In this paper, we first give a sharper theorem regarding the convexity of the Wasserstein metric as the objective function. We then focus on two new issues. One is the necessary normalization of turning seismic signals into probability measures such that the theory of optimal transport applies. The other, which is beyond cycle skipping, is the inversion for parameters below reflecting interfaces. For the first, we propose a class of normalizations and prove several favorable properties for this class. For the latter, we demonstrate that FWI using optimal transport can recover geophysical properties from domains where no seismic waves travel through. We finally illustrate these properties by the realistic application of imaging salt inclusions, which has been a significant challenge in exploration geophysics.


title: Local nucleon-nucleon and three-nucleon interactions within chiral effective field theory

abstract: To obtain an understanding of the structure and reactions of nuclear systems from first principles has been a long-standing goal of nuclear physics. In this respect, few- and many-body systems provide a unique laboratory for studying nuclear interactions. During the past decades, the development of accurate representations of the nuclear force has undergone substantial progress. Particular emphasis has been devoted to chiral effective field theory (EFT), a low-energy effective representation of quantum chromodynamics (QCD). Within chiral EFT, many studies have been carried out dealing with the construction of both the nucleon-nucleon ($N\!N$) and three-nucleon ($3N$) interactions. The aim of the present article is to give a detailed overview of the chiral interaction models that are local in configuration space, and show recent results for nuclear systems obtained by employing these local chiral forces.


title: Semi-Exact Control Functionals From Sard's Method

abstract: The numerical approximation of posterior expected quantities of interest is considered. A novel control variate technique is proposed for post-processing of Markov chain Monte Carlo output, based both on Stein's method and an approach to numerical integration due to Sard. The resulting estimators are proven to be polynomially exact in the Gaussian context, while empirical results suggest the estimators approximate a Gaussian cubature method near the Bernstein-von-Mises limit. The main theoretical result establishes a bias-correction property in settings where the Markov chain does not leave the posterior invariant. Empirical results are presented across a selection of Bayesian inference tasks. All methods used in this paper are available in the R package ZVCV.


title: Exotic Phases of a Higgs-Yukawa Model with Reduced Staggered Fermions

abstract: We investigate the phase structure of a four dimensional SO(4) invariant lattice Higgs-Yukawa model comprising four reduced staggered fermions interacting with a real scalar field. The fermions belong to the fundamental representation of the symmetry group while the three scalar field components transform in the self-dual representation of SO(4). We explore the phase diagram and find evidence of a continuous transition between a phase where the fermions are massless to one where the fermions acquire mass. This transition is not associated with symmetry breaking and there is no obvious local order parameter.


title: Simulating quantum chemistry in the seniority-zero space on qubit-based quantum computers

abstract: Accurate quantum chemistry simulations remain challenging on classical computers for problems of industrially relevant sizes and there is reason for hope that quantum computing may help push the boundaries of what is technically feasible. While variational quantum eigensolver (VQE) algorithms may already turn noisy intermediate scale quantum (NISQ) devices into useful machines, one has to make all efforts to use the scarce quantum resources as efficiently as possible. We combine the so-called seniority-zero, or paired-electron, approximation of computational quantum chemistry with techniques for simulating molecular chemistry on gate-based quantum computers and obtain a very resource efficient quantum simulation algorithm. While some accuracy is lost through the paired-electron approximation, we show that using the freed-up quantum resources for increasing the basis set can lead to more accurate results and reductions in the necessary number of quantum computing runs by several orders of magnitude, already for a simple system like lithium hydride. We also discuss an error mitigation scheme based on post-selection which shows an attractive scaling when the given Hamiltonian format is considered, increasing the viability of its NISQ implementation.


title: Degenerate Sub-keV Fermion Dark Matter from a Solution to the Hubble Tension

abstract: We present a dark sector model addressing both the Hubble tension and the core-cusp problem. The model is based on a hidden Abelian gauge symmetry group with some chiral fermions required by the anomaly cancellation conditions, producing a candidate for the decaying fermion dark matter as a solution to the Hubble tension. Moreover, the sub-keV mass regime and the thermal history of the dark sector help the dark matter candidate resolve the core-cusp problem occurring in the standard $\Lambda$CDM cosmology.


title: Unsupervised Bilingual Lexicon Induction Across Writing Systems

abstract: Recent embedding-based methods in unsupervised bilingual lexicon induction have shown good results, but generally have not leveraged orthographic (spelling) information, which can be helpful for pairs of related languages. This work augments a state-of-the-art method with orthographic features, and extends prior work in this space by proposing methods that can learn and utilize orthographic correspondences even between languages with different scripts. We demonstrate this by experimenting on three language pairs with different scripts and varying degrees of lexical similarity.


title: Fibration theorems for TQ-completion of structured ring spectra

abstract: The aim of this short paper is to establish a spectral algebra analog of the Bousfield-Kan "fibration lemma" under appropriate conditions. We work in the context of algebraic structures that can be described as algebras over an operad $\mathcal{O}$ in symmetric spectra. Our main result is that completion with respect to topological Quillen homology (or TQ-completion, for short) preserves homotopy fibration sequences provided that the base and total $\mathcal{O}$-algebras are connected. Our argument essentially boils down to proving that the natural map from the homotopy fiber to its TQ-completion tower is a pro-$\pi_*$ isomorphism. More generally, we also show that similar results remain true if we replace "homotopy fibration sequence" with "homotopy pullback square."


title: Model for the electrical conductivity in dense plasma mixtures

abstract: A new model for the electrical conductivity of dense plasmas with a mixture of ion species, containing no adjustable parameters, is presented. The model takes the temperature, mass density and relative abundances of the species as input. It takes into account partial ionization, ionic structure, and core-valence orthogonality, and uses quantum mechanical calculations of cross sections. Comparison to an existing high fidelity but computationally expensive method reveals good agreement. The new model is computationally efficient and can reach high temperatures. A new mixing rule is also presented that gives reasonably accurate conductivities for high temperature plasma mixtures.


title: Resilient Finite Time Consensus: A Discontinuous Systems Perspective

abstract: Many algorithms have been proposed in prior literature to guarantee resilient multi-agent consensus in the presence of adversarial attacks or faults. The majority of prior work present excellent results that focus on discrete-time or discretized continuous-time systems. Fewer authors have explored applying similar resilient techniques to continuous-time systems without discretization. These prior works typically consider asymptotic convergence and make assumptions such as continuity of adversarial signals, the existence of a dwell time between switching instances for the system dynamics, or the existence of trusted agents that do not misbehave. In this paper, we expand the study of resilient continuous-time systems by removing many of these assumptions and using discontinuous systems theory to provide conditions for normally-behaving agents with nonlinear dynamics to achieve consensus in finite time despite the presence of adversarial agents.


title: On Implicit Regularization in $\beta$-VAEs

abstract: While the impact of variational inference (VI) on posterior inference in a fixed generative model is well-characterized, its role in regularizing a learned generative model when used in variational autoencoders (VAEs) is poorly understood. We study the regularizing effects of variational distributions on learning in generative models from two perspectives. First, we analyze the role that the choice of variational family plays in imparting uniqueness to the learned model by restricting the set of optimal generative models. Second, we study the regularization effect of the variational family on the local geometry of the decoding model. This analysis uncovers the regularizer implicit in the $\beta$-VAE objective, and leads to an approximation consisting of a deterministic autoencoding objective plus analytic regularizers that depend on the Hessian or Jacobian of the decoding model, unifying VAEs with recent heuristics proposed for training regularized autoencoders. We empirically verify these findings, observing that the proposed deterministic objective exhibits similar behavior to the $\beta$-VAE in terms of objective value and sample quality.


title: Gravitational Lensing in Rotating and Twisting Universes

abstract: Gravitational lensing caused by the gravitational field of massive objects has been studied and acknowledged for a long period of time. In this paper, however, we propose a different mechanism where the bending of light stems from the non-linear interaction of gravitational, electromagnetic and axion waves that creates the high curvature zone in the space-time fabric. The striking distinction in the present study is that in contrast to the convex lensing in the gravitational field of a massive object, hyperbolic nature of the high curvature zone of the background space-time may give rise to concave lensing. Expectedly detection of this kind of lensing becomes possible through satellite detectors.


title: Polariton fluids for analogue gravity physics

abstract: Analogue gravity enables the study of fields on curved spacetimes in the laboratory. There are numerous experimental platforms in which amplification at the event horizon or the ergoregion has been observed. Here, we demonstrate how optically generating a defect in a polariton microcavity enables the creation of one- and two-dimensional, transsonic fluid flows. We show that this highly tuneable method permits the creation of sonic horizons. Furthermore, we present a rotating geometry akin to the water-wave bathtub vortex. These experiments usher-in the possibility of observing stimulated as well as spontaneous amplification by the Hawking, Penrose and Zeld'ovich effects in fluids of light.


title: Design Principles and Clinician Preferences for Pharmacogenomic Clinical Decision Support Alerts

abstract: OBJECTIVE: To better understand clinician needs and preferences for the display of pharmacogenomic (PGx) information in clinical decision support (CDS) tools. MATERIALS AND METHODS: We developed a semi-structured interview to collect feedback and preferences in six key areas of PGx CDS design, from clinicians who had prior experience with live PGx CDS tools. Eight clinicians from Northwestern Medicine's (NM) General Internal Medicine clinic participated in the study. RESULTS: Clinicians expressed preference for interruptive pop-up alerts during order entry, brief descriptions of relevant drug-gene interactions, and a clear and specific recommended alternative course of action when a medication is contraindicated. They did not wish to see detailed genetic data, preferring phenotypic information predicted from the genotype. Nor did they wish to be interrupted when genetic test results do not indicate a change in treatment plan. Clinicians reported little familiarity with Clinical Pharmacogenetic Implementation Consortium prescribing recommendations but reported trusting recommendations of their professional societies and resources like UpToDate. Analysis of unstructured comments concurred with structured results, indicating a general uncertainty among participants around how to interpret and apply PGx information in practice. DISCUSSION: Results point to several underlying principles that can inform future PGx CDS alert designs: Be Specific and Actionable; Be Brief; Display Phenotypes not Genotypes; Rely on Sources Clinicians Already Trust; and, Be Adaptable to Learning Effects. CONCLUSION: This study is part of a broader socio-technical design approach to PGx CDS design underway at NM and provides a baseline for future PGx CDS development. Designs based on these results have the potential to improve clinician education and adherence levels, and to improve patient outcomes.


title: Open IT tools for inquiry-based physics education

abstract: In case of inquiry-based education, the priority is the pupil's activity, developing his/her practical and research skills. We can use the knowledge and skills obtained by the pupil from other subjects. Applied informatics into physics seems to be a suitable application of digital literacy. Open IT tools are the ones of the available IT tools that bring several benefits to promote freedom of inquiry. An overview of open IT tools for physics education includes microprocessor kits with a variety of analog and digital electronic components, a large number of application software, simulation software, programming languages, libraries and development environments, web tools, virtual laboratories, data catalogs, hypertext documents and other auxiliary resources. Such open IT tools together with suitable open educational resources will be applicable in physics, mathematics, informatics and other subjects.


title: Black hole, neutron star and white dwarf merger rates in AGN disks

abstract: Advanced LIGO \& Advanced Virgo are detecting a large number of binary stellar origin black hole (BH) mergers. A promising channel for accelerated BH merger lies in active galactic nucleus (AGN) disks of gas around super-masssive black holes. Here we investigate the relative number of compact object mergers in AGN disk models, including BH, neutron stars (NS) and white dwarfs, via Monte Carlo simulations. We find the number of all merger types in the bulk disk grows $\propto t^{1/3}$ which is driven by the Hill sphere of the more massive merger component. Median mass ratios of NS-BH mergers in AGN disks are $\tilde{q}=0.07\pm 0.06(0.14\pm 0.07)$ for mass functions (MF) $M^{-1(-2)}$. If a fraction $f_{AGN}$ of the observed rate of BH-BH mergers (${\cal}{R}_{BH-BH}$) come from AGN, we expect a rate of NS-BH (NS-NS) mergers from the AGN channel of ${\cal}{R}_{BH-NS} \sim f_{AGN}[10,300]\rm{Gpc}^{-3} \rm{yr}^{-1}$,(${\cal}{R}_{NS-NS} \leq f_{AGN}400\rm{Gpc}^{-3} \rm{yr}^{-1}$). Allowing for the ratio of NS-NS/BH-BH LIGO search volumes, from preliminary O3 results we can rule out the AGN channel as a dominant contribution to observed NS-NS mergers. The number of lower mass gap events in this channel is a strong function of the nuclear MF and mass segregation efficiency. Compact object merger ratios derived from LIGO can restrict models of MF, mass segregation and populations embedded in AGN disks. The expected number of EM counterparts to NS-BH mergers in AGN disks at $z<1$ is $\sim [30,900]{\rm{yr}}^{-1}(f_{AGN}/0.1)$. EM searches for flaring due to embedded events in large AGN surveys will complement LIGO constraints on AGN models and the embedded populations that must live in them.


title: Design Principles Developed through User-Centered and Socio-Technical Methods Improve Clinician Satisfaction, Speed, and Confidence in Pharmacogenomic Clinical Decision Support

abstract: OBJECTIVE: To design and evaluate new pharmacogenomic (PGx) clinical decision support (CDS) alerts, built to adhere to PGx CDS design principles developed through socio-technical approaches. MATERIALS AND METHODS: Based on previously identified design principles, we created 11 new PGx CDS alert designs and developed an interactive web application containing realistic clinical scenarios and user workflows that mimicked a real-world EHR system. We recruited General Internal Medicine and Cardiology clinicians from Northwestern Medicine and recorded their interactions with the original and new designs. We measured clinician response, satisfaction, speed, and confidence through questionnaires and analysis of the recordings. RESULTS: The study included 12 clinicians. Participants were significantly more satisfied (p=0.0000001), faster (p=0.009), and more confident (p<.05) with the new designs than the original ones. The study lacked statistical power to determine whether prescribing accuracy was improved, but participants were no less accurate, and clinical actions were more concordant with alert interactions (p=0.004) with the new designs. We found a significant learning curve associated with the original designs, which was eliminated with the new designs. DISCUSSION: This study successfully demonstrates that socio-technical and user-centered design techniques can improve PGx CDS alert designs. Best practices for PGx CDS design are limited in the literature, with few effectiveness studies available. These results can help guide future PGx CDS implementations to be more clinician friendly and less time-consuming. CONCLUSION: The results of this study support the PGx CDS design principles we proposed in previous work. As a next step, the new designs should be implemented in a live setting for further validation.


title: Observation of Collisions between Two Ultracold Ground-State CaF Molecules

abstract: We measure inelastic collisions between ultracold CaF molecules by combining two optical tweezers, each containing a single molecule. We observe collisions between $^2\Sigma$ CaF molecules in the absolute ground state $|X,v=0, N=0,F=0\rangle$, and in excited hyperfine and rotational states. In the absolute ground state, we find a two-body loss rate of $7(4) \times 10^{-11} \text{cm}^{3}/\text{s}$, which is below, but close to the predicted universal loss rate.


title: Studying Lagrangian theories with machine learning: a toy model

abstract: The existence or not of pathologies in the context of Lagrangian theory is studied with the aid of Machine Learning algorithms. Using an example in the framework of classical mechanics, we make a proof of concept, that the construction of new physical theories using machine learning is possible. Specifically, we utilize a fully-connected, feed-forward neural network architecture, aiming to discriminate between ``healthy'' and ``non-healthy'' Lagrangians, without explicitly extracting the relevant equations of motion. The network, after training, is used as a fitness function in the concept of a genetic algorithm and new healthy Lagrangians are constructed. These new Lagrangians are different from the Lagrangians contained in the initial data set. Hence, searching for Lagrangians possessing a number of pre-defined properties is significantly simplified within our approach. The framework employed in this work can be used to explore more complex physical theories, such as generalizations of General Relativity in gravitational physics, or constructions in solid state physics, in which the standard procedure can be laborious.


title: A direct proof of APN-ness of the Kasami functions

abstract: Using recent results on solving the equation $X^{2^k+1}+X+a=0$ over a finite field $\mathbb{F}_{2^n}$, we address an open question raised by the first author in WAIFI 2014 concerning the APN-ness of the Kasami functions $x\mapsto x^{2^{2k}-2^k+1}$ with $gcd(k,n)=1$, $x\in\mathbb{F}_{2^n}$.


title: Exact solution for the Anisotropic Ornstein-Uhlenbeck Process

abstract: Active Matter models commonly consider particles with overdamped dynamics subject to a force (speed) with constant modulus and random direction. Some models include also random noise in particle displacement (Wiener process) resulting in a diffusive motion at short time scales. On the other hand, Ornstein-Uhlenbeck processes consider Langevin dynamics for the particle velocity and predict a motion that is not diffusive at short time scales. However, experiments show that migrating cells may present a varying speed as well as a short-time diffusive behavior. While Ornstein-Uhlenbeck processes can describe the varying speed, Active Mater models can explain the short-time diffusive behavior. Isotropic models cannot explain both: short-time diffusion renders instantaneous velocity ill-defined, hence impeding dynamical equations that consider velocity time-derivatives. On the other hand, both models apply for migrating biological cells and must, in some limit, yield the same observable predictions. Here we propose and analytically solve an Anisotropic Ornstein-Uhlenbeck process that considers polarized particles, with a Langevin dynamics for the particle movement in the polarization direction while following a Wiener process for displacement in the orthogonal direction. Our characterization provides a theoretically robust way to compare movement in dimensionless simulations to movement in dimensionful experiments, besides proposing a procedure to deal with inevitable finite precision effects in experiments or simulations.


title: Symplectic Manifolds and Isomonodromic Deformations

abstract: We study moduli spaces of meromorphic connections (with arbitrary order poles) over Riemann surfaces together with the corresponding spaces of monodromy data (involving Stokes matrices). Natural symplectic structures are found and described both explicitly and from an infinite dimensional viewpoint (generalising the Atiyah-Bott approach). This enables us to give an intrinsic symplectic description of the isomonodromic deformation equations of Jimbo, Miwa and Ueno, thereby putting the existing results for the six Painleve equations and Schlesinger's equations into a uniform framework.