title: Black holes seeding cosmological phase transitions

abstract: We consider a generic first-order phase transition at finite temperature and investigate to what extent a population of primordial black holes, of variable masses, can affect the rate of bubble nucleation. Using a thin-wall approximation, we construct the Euclidean configurations that describe transition at finite temperature. After the transition, the remnant black hole mass is dictated dynamically by the equations of motion. The transition exponent is computed and displays an explicit dependence on temperature. We find the configuration with the lowest Euclidean action to be static and $O(3)$ symmetric; therefore, the transition takes place via thermal excitation. The transition exponent exhibits a strong dependence on the seed mass black hole, $M_+$, being almost directly proportional. A new nucleation condition in the presence of black holes is derived and the nucleation temperature is compared to the familiar flat-space result, i.e. $S_3/T$. For an electroweak-like phase transition it is possible to enhance the nucleation rate if $M_+ \lesssim 10^{15} M_{\rm P}$. Finally, we outline the possible transition scenarios and the consequences for the power spectrum of stochastic gravitational waves produced due to the first-order phase transition.


title: Bright high-purity quantum emitters in aluminium nitride integrated photonics

abstract: Solid-state quantum emitters (QEs) are fundamental in photonic-based quantum information processing. There is strong interest to develop high-quality QEs in III-nitride semiconductors because of their sophisticated manufacturing driven by large and growing applications in optoelectronics, high voltage power transistors, and microwave amplifiers. Here, we report the generation and direct integration of QEs in an aluminium nitride-based photonic integrated circuit platform. For individual waveguide-integrated QEs, we measure an off-chip count rate exceeding $6 \times 10^{4}$ counts per second (cps) (saturation rate > $8.6 \times 10^{4}$ cps). In an unpatterned thin-film sample, we measure antibunching with $g^{(2)}(0) \sim 0.05$ and photon count rates exceeding $8 \times 10^{5}$ cps (saturation rate > $1 \times 10^{6}$ cps). Although spin and detailed optical linewidth measurements are left for future work, these results already show the potential for high-quality QEs monolithically integrated in a wide range of III-nitride device technologies that would enable new quantum device opportunities and industrial scalability.


title: A Guaranteed Discovery at Future Muon Colliders

abstract: The longstanding muon g-2 anomaly may indicate the existence of new particles that couple to muons, which could either be light (< GeV) and weakly coupled, or heavy (>> 100 GeV) with large couplings. If light new states are responsible, upcoming intensity frontier experiments will discover further evidence of new physics. However, if heavy particles are responsible, many candidates are beyond the reach of existing colliders. We show that, if the g-2 anomaly is confirmed and no explanation is found at low-energy experiments, a high-energy muon collider program is guaranteed to make fundamental discoveries about our universe. New physics scenarios that account for the anomaly can be classified as either "Singlet" or "Electroweak" (EW) models, involving only EW singlets or new EW-charged states respectively. We argue that a TeV-scale future muon collider will discover all possible singlet model solutions to the anomaly. If this does not yield a discovery, the next step would be a O(10 TeV) muon collider. Such a machine would either discover new particles associated with high-scale EW model solutions to the anomaly, or empirically prove that nature is fine-tuned, both of which would have profound consequences for fundamental physics.


title: A nonlocal isoperimetric problem with density perimeter

abstract: We consider the minimization of an energy functional given by the sum of a density perimeter and a nonlocal interaction of Riesz type with exponent $\alpha$, under volume constraint, where the strength of the nonlocal interaction is controlled by a parameter $\gamma$. We show that for a wide class of density functions the energy admits a minimizer for any value of $\gamma$. Moreover these minimizers are bounded. For monomial densities of the form $|x|^p$ we prove that when $\gamma$ is sufficiently small the unique minimizer is given by the ball of fixed volume. In contrast with the constant density case, here the $\gamma\to 0$ limit corresponds, under a suitable rescaling, to a small mass $m=|\Omega|\to 0$ limit when $p<d-\alpha+1$, but to a large mass $m\to\infty$ for powers $p>d-\alpha+1$.


title: Comparison of improved TMD and CGC frameworks in forward quark dijet production

abstract: For studying small-$x$ gluon saturation in forward dijet production in high-energy dilute-dense collisions, the improved TMD (ITMD) factorization formula was recently proposed. In the Color Glass Condensate (CGC) framework, it represents the leading term of an expansion in inverse powers of the hard scale. It contains the leading-twist TMD factorization formula relevant for small gluon's transverse momentum $k_t$, but also incorporates an all-order resummation of kinematical twists, resulting in a proper matching to high-energy factorization at large $k_t$. In this paper, we evaluate the accuracy of the ITMD formula quantitatively, for the case of quark dijet production in high-energy proton-proton($p+p$) and proton-nucleus ($p+A$) collisions at LHC energies. We do so by comparing the quark-antiquark azimuthal angle $\Delta\phi$ distribution to that obtained with the CGC formula. For a dijet with each quark momentum $p_t$ much larger than the target saturation scale, $Q_s$, the ITMD formula is a good approximation to the CGC formula in a wide range of azimuthal angle. It becomes less accurate as the jet $p_t$'s are lowered, as expected, due to the presence of genuine higher-twists contributions in the CGC framework, which represent multi-body scattering effects absent in the ITMD formula. We find that, as the hard jet momenta are lowered, the accuracy of ITMD start by deteriorating at small angles, in the high-energy-factorization regime, while in the TMD regime near $\Delta\phi=\pi$, very low values of $p_t$ are needed to see differences between the CGC and the ITMD formula. In addition, the genuine twists corrections to ITMD become visible for higher values of $p_t$ in $p+A$ collisions, compared to $p+p$ collisions, signaling that they are enhanced by the target saturation scale.


title: The Assembly History of M87 Through Radial Variations in Chemical Abundances of its Field Star And Globular Cluster Populations

abstract: We present an extensive study of spectroscopically-derived chemical abundances for M87 and its globular cluster (GC) system. Using observations from the Mitchell spectrograph at McDonald, LRIS at Keck, and Hectospec on the MMT, we derive new metallicity gradients from $\sim 2$ to $140$ kpc. We use a novel hierarchical statistical framework to simultaneously separate the GC system into subpopulations while measuring the metallicity gradients of those subpopulations. We create physically-motivated spectral stacks of the GC subpopulations by leveraging the output of this statistical framework to perform the first application of abundance tagging in a massive ETG to better constrain the origins of the GC subpopulations and, thus, the assembly history of M87. We find a metal-poor, $\alpha$-enhanced population of GCs in both in the inner and outer halo unanticipated by current cosmological simulations of galaxy evolution. We use the remarkably flat metallicity gradients we find for both the metal-rich and metal-poor GC subpopulations in the inner halo as tentative evidence that some amount of the metal-poor GCs formed directly in the halo of M87 at high redshift.


title: TinyRadarNN: Combining Spatial and Temporal Convolutional Neural Networks for Embedded Gesture Recognition with Short Range Radars

abstract: This work proposes a low-power high-accuracy embedded hand-gesture recognition algorithm targeting battery-operated wearable devices using low power short-range RADAR sensors. A 2D Convolutional Neural Network (CNN) using range frequency Doppler features is combined with a Temporal Convolutional Neural Network (TCN) for time sequence prediction. The final algorithm has a model size of only 46 thousand parameters, yielding a memory footprint of only 92 KB. Two datasets containing 11 challenging hand gestures performed by 26 different people have been recorded containing a total of 20,210 gesture instances. On the 11 hand gesture dataset, accuracies of 86.6% (26 users) and 92.4% (single user) have been achieved, which are comparable to the state-of-the-art, which achieves 87% (10 users) and 94% (single user), while using a TCN-based network that is 7500x smaller than the state-of-the-art. Furthermore, the gesture recognition classifier has been implemented on a Parallel Ultra-Low Power Processor, demonstrating that real-time prediction is feasible with only 21 mW of power consumption for the full TCN sequence prediction network, while a system-level power consumption of less than 100 mW is achieved. We provide open-source access to all the code and data collected and used in this work on tinyradar.ethz.ch.


title: Irksome: Automating Runge--Kutta time-stepping for finite element methods

abstract: While implicit Runge--Kutta methods possess high order accuracy and important stability properties, implementation difficulties and the high expense of solving the coupled algebraic system at each time step are frequently cited as impediments. We present IIrksome, a high-level library for manipulating UFL (Unified Form Language) expressions of semidiscrete variational forms to obtain UFL expressions for the coupled Runge--Kutta stage equations at each time step. Irksome works with the Firedrake package to enable the efficient solution of the resulting coupled algebraic systems. Numerical examples confirm the efficacy of the software and our solver techniques for various problems.


title: Determination of the escape velocity of the Milky Way using a proper motion selected halo sample

abstract: The {\it Gaia} mission has provided the largest catalogue ever of sources with tangential velocity information. However, using this catalogue for dynamical studies is difficult because most of the stars lack line-of-sight velocity measurements. Recently, we presented a selection of $\sim 10^7$ halo stars with accurate distances that have been selected based on their photometry and proper motions. Using this sample, we model the tail of the velocity distribution with a power-law distribution, a commonly used approach first established by \cite{Leonard1990THESPEED}. For the first time ever we use tangential velocities measured accurately for an unprecedented number of halo stars to estimate the escape velocity. In the solar neighbourhood, we obtain a very precise estimate of the escape velocity which is $497^{+8}_{-8}~{\rm km/s}$. This estimate is most likely biased low, our best guess is by 10\%. As a result, the true escape velocity most likely is closer to $550~{\rm km/s}$. The escape velocity directly constrains the total mass of the Milky Way. To find the best fitting halo mass and concentration parameter we adjusted the dark (spherical NFW) halo of a realistic Milky Way potential while keeping the circular velocity at the solar radius fixed at $v_c(R_\odot) = 232.8~{\rm km/s}$. The resulting halo parameters are $M_{200}^{+10\%} = 1.11^{+0.08}_{-0.07} \cdot10^{12} ~{\rm M}_\odot$ and concentration parameter $c^{+10\%} = 11.8^{+0.3}_{-0.3}$, where we use the explicit notation to indicate that these have been corrected for the 10\% bias. The slope of the escape velocity with galactocentric distance is as expected in the inner Galaxy based on Milky Way models. Curiously, we find a disagreement beyond the solar radius which is likely an effect of a change in the shape of the velocity distribution and could be related to the presence of velocity clumps.


title: Transactions on Red-black and AVL trees in NVRAM

abstract: Byte-addressable non-volatile memory (NVRAM) supports persistent storage with low latency and high bandwidth. Complex data structures in it ought to be updated transactionally, so that they remain recoverable at all times. Traditional database technologies such as keeping a separate log, a journal, or shadow data work on a coarse-grained level, where the whole transaction is made visible using a final atomic update operation. These methods typically need significant additional space overhead and induce non-trivial overhead for log pruning, state maintenance, and resource (de-)allocation. Thus, they are not necessarily the best choice for NVRAM, which supports fine-grained, byte-addressable access. We present a generic transaction mechanism to update dynamic complex data structures `in-place' with a constant memory overhead. It is independent of the size of the data structure. We demonstrate and evaluate our approach on Red-Black Trees and AVL Trees with a redo log of constant size (4 resp. 2 cache lines). The redo log guarantees that each accepted (started) transaction is executed eventually despite arbitrary many system crashes and recoveries in the meantime. We update complex data structures in local and remote NVRAM providing exactly once semantics and durable linearizability for multi-reader single-writer access. To persist data, we use the available processor instructions for NVRAM in the local case and remote direct memory access (RDMA) combined with a software agent in the remote case.


title: An Introduction to Effective Field Theories

abstract: These notes are an overview of effective field theory (EFT) methods. I discuss toy model EFTs, chiral perturbation theory, Fermi liquid theory, and non-relativistic QED, and use these examples to introduce a variety of EFT concepts, including: matching a tree and loop level; summation of large logarithms; naturalness problems; spurion fields; coset construction; Wess-Zumino-Witten terms; chiral and gauge anomalies; field rephasings; and method of regions. These lecture notes were prepared for the 2nd Joint ICTP-Trieste/ICTP-SAIFR school on Particle Physics, Sao Paulo, Brazil, June 22 - July 3, 2020.


title: Diffusion approximation for noise-induced evolution of first integrals in multifrequency systems

abstract: We consider fast oscillating perturbations of dynamical systems in regions where one can introduce action-angle type coordinates. In an appropriate time scale, a diffusion approximation of the first-integrals evolution is described under the assumption that the set of resonance tori is small enough. If the action-angle coordinates can be introduced just piece-wise , the limiting diffusion process should be considered on an open book space. Such a process can be described by differential operators, one on each page, supplemented by some gluing conditions on the binding of the open book.


title: Multifrequency study of HH 137 and HH 138: Discovering new knots and molecular outflows with Gemini and APEX

abstract: We present a multi-wavelength study of two HH objects (137 and 138) that may be associated. We use Gemini H$_2$ (2.12 $\mu$m) and K (2.2 $\mu$m) images, as well as APEX molecular line observations and Spitzer image archives. Several H$_2$ knots, linked to the optical chain of knots of HH 137, are identified in the Gemini and Spitzer 4.5 $\mu$m images. New shock excited regions related to the optical knots delineating HH 138 are also reported. In addition, a bright 4.5 $\mu$m 0.09 pc-long arc-shaped structure, roughly located mid-way between HH 137 and HH 138, is found to be associated with two Spitzer Class I/II objects, which are likely to be the exciting stars. These sources are almost coincident with a high-density molecular clump detected in $^{12}$CO(3-2), $^{13}$CO(3-2), C$^{18}$O(3-2), HCO$^{+}$(3-2) and HCN(3-2) molecular lines with an LTE mass of 36 M$_{\odot}$. The $^{12}$CO(3-2) emission distribution over the observed region reveals molecular material underlying three molecular outflows. Two of them (outflows 1 and 2) are linked to all optical knots of HH 137 and HH 138 and to the H$_2$ and 4.5 $\mu$m shock emission knots. In fact, the outflow 2 shows an elongated $^{12}$CO blue lobe that coincides with all the H$_2$ knots of HH 137 which end at a terminal H$_2$ bow shock. We propose a simple scenario that connects the outflows to the dust clumps detected in the region. A third possible outflow is located to the north-east projected towards a secondary weak and cold dust clump.


title: Affine Deligne-Lusztig varieties and folded galleries governed by chimneys

abstract: We characterize the nonemptiness and dimension problems for an affine Deligne-Lusztig variety $X_x(b)$ in the affine flag variety in terms of galleries that are positively folded with respect to a chimney. If the parabolic subgroup associated to the Newton point of b has rank 1, we then prove nonemptiness for a certain class of Iwahori-Weyl group elements x by explicitly constructing such galleries.


title: Gravity/Ensemble Duality

abstract: For the first time, a gravitational calculation was recently shown to yield the Page curve for the entropy of Hawking radiation, consistent with unitary evolution. However, the calculation takes as essential input Hawking's result that the radiation entropy becomes large at late times. We call this apparent contradiction the state paradox. We exhibit its manifestations in standard and doubly-holographic settings, with and without an external bath. We clarify which version(s) of the Ryu-Takayanagi prescription apply in each setting. We show that the two possible homology rules in the presence of a braneworld generate a bulk dual of the state paradox. The paradox is resolved if the gravitational path integral computes averaged quantities in a suitable ensemble of unitary theories, a possibility supported independently by several recent developments.


title: All states are universal catalysts in quantum thermodynamics

abstract: Quantum catalysis is a fascinating concept which demonstrates that certain transformations can only become possible when given access to a specific resource that has to be returned unaffected. It was first discovered in the context of entanglement theory and since then applied in a number of resource-theoretic frameworks, including quantum thermodynamics. Although in that case the necessary (and sometimes also sufficient) conditions on the existence of a catalyst are known, almost nothing is known about the precise form of the catalyst state required by the transformation. In particular, it is not clear whether it has to have some special properties or be finely tuned to the desired transformation. In this work we describe a surprising property of multi-copy states: we show that in resource theories governed by majorization all resourceful states are catalysts for all allowed transformations. In quantum thermodynamics this means that the so-called "second laws of thermodynamics" do not require a fine-tuned catalyst but rather any state, given sufficiently many copies, can serve as a useful catalyst. These analytic results are accompanied by several numerical investigations that indicate that neither a multi-copy form nor a very large dimension catalyst are required to activate most allowed transformations catalytically.


title: $H_0$ Ex Machina: Vacuum Metamorphosis and Beyond $H_0$

abstract: We do not solve tensions with concordance cosmology; we do obtain $H_0\approx 74\,$km/s/Mpc from CMB+BAO+SN data in our model, but that is not the point. Discrepancies in Hubble constant values obtained by various astrophysical probes should not be viewed in isolation. While one can resolve at least some of the differences through either an early time transition or late time transition in the expansion rate, these introduce other changes. We advocate a holistic approach, using a wide variety of cosmic data, rather than focusing on one number, $H_0$. Vacuum metamorphosis, a late time transition physically motivated by quantum gravitational effects and with the same number of parameters as \lcdm, can successfully give a high $H_0$ value from cosmic microwave background data but fails when combined with multiple distance probes. We also explore the influence of spatial curvature, and of a conjoined analysis of cosmic expansion and growth.


title: A Gap Between the Hypergraph and Stabilizer Entropy Cones

abstract: It was recently found that the stabilizer and hypergraph entropy cones coincide for four parties, leading to a conjecture of their equivalence at higher party numbers. In this note, we show this conjecture to be false by proving new inequalities obeyed by all hypergraph entropy vectors that exclude particular stabilizer states on six qubits. By further leveraging this connection, we improve the characterization of stabilizer entropies and show that all linear rank inequalities at five parties, except for classical monotonicity, form facets of the stabilizer cone. Additionally, by studying minimum cuts on hypergraphs, we prove some structural properties of hypergraph representations of entanglement and generalize the notion of entanglement wedge nesting in holography.


title: Probabilistic definition of the perturbative theoretical uncertainty from missing higher orders

abstract: We consider the problem of quantifying the uncertainty on theoretical predictions based on perturbation theory due to missing higher orders. The most widely used approach, scale variation, is largely arbitrary and it has no probabilistic foundation, making it not suitable for robust data analysis. In 2011, Cacciari and Houdeau proposed a model based on a Bayesian approach to provide a probabilistic definition of the theory uncertainty from missing higher orders. In this work, we propose an improved version of the Cacciari-Houdeau model, that overcomes some limitations. In particular, it performs much better in case of perturbative expansions with large high-order contributions (as it often happens in QCD). In addition, we propose an alternative model based on the same idea of scale variation, which overcomes some of the shortcomings of the canonical approach, on top of providing a probabilistically-sound result. Moreover, we address the problem of the dependence of theoretical predictions on unphysical scales (such as the renormalization scale), and propose a solution to obtain a scale-independent result within the probabilistic framework. We validate these methods on expansions with known sums, and apply them to a number of physical observables in particle physics. We also investigate some variations, improvements and combinations of the models. We believe that these methods provide a powerful tool to reliably estimate theory uncertainty from missing higher orders that can be used in any physics analysis. The results of this work are easily accessible through a public code named THunc.


title: Reductions of $2$-dimensional semi-stable representations with large $\mathcal L$-invariant

abstract: We determine reductions of 2-dimensional, irreducible, semi-stable, and non-crystalline representations of $\mathrm{Gal}(\overline{\mathbb Q}_p/\mathbb Q_p)$ with Hodge--Tate weights $0 < k-1$ and with $\mathcal L$-invariant whose $p$-adic norm is sufficiently large, depending on $k$. Our main result provides the first systematic examples of the reductions for $k \geq p$.


title: Rapid Transitions with Robust Accelerated Delayed Self Reinforcement for Consensus-Based Networks

abstract: Rapid transitions are important for quick response of consensus-based, multi-agent networks to external stimuli. While high-gain can increase response speed, potential instability tends to limit the maximum possible gain, and therefore, limits the maximum convergence rate to consensus during transitions. Since the update law for multi-agent networks with symmetric graphs can be considered as the gradient of its Laplacian-potential function, Nesterov-type accelerated-gradient approaches from optimization theory, can further improve the convergence rate of such networks. An advantage of the accelerated-gradient approach is that it can be implemented using accelerated delayed-self-reinforcement (A-DSR), which does not require new information from the network nor modifications in the network connectivity. However, the accelerated-gradient approach is not directly applicable to general directed graphs since the update law is not the gradient of the Laplacian-potential function. The main contribution of this work is to extend the accelerated-gradient approach to general directed graph networks, without requiring the graph to be strongly connected. Additionally, while both the momentum term and outdated-feedback term in the accelerated-gradient approach are important in general, it is shown that the momentum term alone is sufficient to achieve balanced robustness and rapid transitions without oscillations in the dominant mode, for networks whose graph Laplacians have real spectrum. Simulation results are presented to illustrate the performance improvement with the proposed Robust A-DSR of 40% in structural robustness and 50% in convergence rate to consensus, when compared to the case without the A-DSR. Moreover, experimental results are presented that show a similar 37% faster convergence with the Robust A-DSR when compared to the case without the A-DSR.


title: Neutron Star Equation of State in light of GW190814

abstract: The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. We construct heavy neutron stars by introducing non-trivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.


title: Optimization Landscape of Tucker Decomposition

abstract: Tucker decomposition is a popular technique for many data analysis and machine learning applications. Finding a Tucker decomposition is a nonconvex optimization problem. As the scale of the problems increases, local search algorithms such as stochastic gradient descent have become popular in practice. In this paper, we characterize the optimization landscape of the Tucker decomposition problem. In particular, we show that if the tensor has an exact Tucker decomposition, for a standard nonconvex objective of Tucker decomposition, all local minima are also globally optimal. We also give a local search algorithm that can find an approximate local (and global) optimal solution in polynomial time.


title: Measurements of the Time-Dependent Cosmic-Ray Sun Shadow with Seven Years of IceCube Data -- Comparison with the Solar Cycle and Magnetic Field Models

abstract: Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. In this paper, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the Moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the Moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4sigma level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of $3\sigma$ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-gamma Experiment, a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.


title: Spatial Separation of the Conformers of Methyl Vinyl Ketone

abstract: Methyl vinyl ketone (C$_4$H$_6$O) is a volatile, labile organic compound of importance in atmospheric chemistry. We prepared a molecular beam of methyl vinyl ketone with a rotational temperature of 1.2(2)~K and demonstrated the spatial separation of the \emph{s-cis} and \emph{s-trans} conformers of methyl vinyl ketone using the electrostatic deflector. The resulting sample density was $1.5(2)\times10^{8}~\text{cm}^{-3}$ for the direct beam in the laser ionization region. These conformer-selected methyl vinyl ketone samples are well suited for conformer-specific chemical reactivity studies such as in Diels-Alder cycloaddition reactions.


title: pyCallisto: A Python Library To Process The CALLISTO Spectrometer Data

abstract: CALLISTO is a radio spectrometer designed to monitor the transient radio emissions/bursts originated from the solar corona in the frequency range $45-870$ MHz. At present, there are $\gtrsim 150$ stations (together forms an e-CALLISTO network) around the globe continuously monitoring the Sun 24 hours a day. We have developed a pyCallisto, a python library to process the CALLISTO data observed by all stations of the e-CALLISTO network. In this article, we demonstrate various useful functions that are routinely used to process the CALLISTO data with suitable examples. This library is not only efficient in processing the data but plays a significant role in developing automatic classification algorithms of different types of solar radio bursts.


title: Quantum well states in fractured crystals of the heavy fermion material CeCoIn$_5$

abstract: Quantum well states appear in metallic thin films due to the confinement of the wave function by the film interfaces. Using angle-resolved photoemission spectroscopy, we unexpectedly observe quantum well states in fractured single crystals of CeCoIn$_5$. We confirm that confinement occurs by showing that these states' binding energies are photon-energy independent and are well described with a phase accumulation model, commonly applied to quantum well states in thin films. This indicates that atomically flat thin films can be formed by fracturing hard single crystals. For the two samples studied, our observations are explained by free-standing flakes with thicknesses of 206 and 101 \r{A}. We extend our analysis to extract bulk properties of CeCoIn$_5$. Specifically, we obtain the dispersion of a three-dimensional band near the zone center along in-plane and out-of-plane momenta. We establish part of its Fermi surface, which corresponds to a hole pocket centered at $\Gamma$. We also reveal a change of its dispersion with temperature, a signature that may be caused by the Kondo hybridization.


title: A Compact Gated-Synapse Model for Neuromorphic Circuits

abstract: This work reports a compact behavioral model for gated-synaptic memory. The model is developed in Verilog-A for easy integration into computer-aided design of neuromorphic circuits using emerging memory. The model encompasses various forms of gated synapses within a single framework and is not restricted to only a single type. The behavioral theory of the model is described in detail along with a full list of the default parameter settings. The model includes parameters such as a device's ideal set time, threshold voltage, general evolution of the conductance with respect to time, decay of the device's state, etc. Finally, the model's validity is shown via extensive simulation and fitting to experimentally reported data on published gated-synapses.


title: Correlation functions of composite Ramond fields in deformed D1-D5 orbifold SCFT$_2$

abstract: We study two families of composite twisted Ramond fields (made by products of two operators) in the $\cal {N}=(4,4)$ supersymmetric D1-D5 SCFT$_2$ deformed by a marginal modulus operator away from its $(T^4)^N/ S_N$ free orbifold point. We construct the large-$N$ contributions to the four-point functions with two composite operators and two deformation fields. These functions allow us to derive short-distance OPE limits and to calculate the anomalous dimensions of the composite operators. We demonstrate that one can distinguish two sets of composite Ramond states with twists $m_1$ and $m_2$: protected states, for which $m_1+m_2=N$, and "lifted" states for which $m_1+m_2<N$. The latter require an appropriate renormalisation. We also derive the leading order corrections to their two-point functions, and to their three-point functions with the deformation operator.


title: Realizing the Frenkel-Kontorova model with Rydberg-dressed atoms

abstract: We propose a method to realize the Frenkel-Kontorova model using an array of Rydberg dressed atoms. Our platform can be used to study this model with a range of realistic interatomic potentials. In particular, we concentrate on two types of interaction potentials: a springlike potential and a repulsive long-range potential. We numerically calculate the phase diagram for such systems and characterize the Aubry-like and commensurate-incommensurate phase transitions. Experimental realizations of this system that are possible to achieve using current technology are discussed.


title: Exponential multivalued forbidden configurations

abstract: The forbidden number $\mathrm{forb}(m,F)$, which denotes the maximum number of unique columns in an $m$-rowed $(0,1)$-matrix with no submatrix that is a row and column permutation of $F$, has been widely studied in extremal set theory. Recently, this function was extended to $r$-matrices, whose entries lie in $\{0,1,\dots,r-1\}$. The combinatorics of the generalized forbidden number is less well-studied. In this paper, we provide exact bounds for many $(0,1)$-matrices $F$, including all $2$-rowed matrices when $r > 3$. We also prove a stability result for the $2\times 2$ identity matrix. Along the way, we expose some interesting qualitative differences between the cases $r=2$, $r = 3$, and $r > 3$.


title: Imaging and Spectroscopy of Domains of the Cellular Membrane by Photothermal-Induced Resonance

abstract: We use photothermal induced resonance (PTIR) imaging and spectroscopy, in resonant and non-resonant mode, to study the cytoplasmic membrane and surface of intact cells. Non-resonant PTIR images apparently provide rich details of the cell surface. However we show that non-resonant image contrast does not arise from the infrared absorption of surface molecules and is instead dominated by the mechanics of tip-sample contact. In contrast, spectra and images of the cellular surface can be selectively obtained by tuning the pulsing structure of the laser to restrict thermal wave penetration to the surface layer. Resonant PTIR images reveal surface structures and domains that range in size from about 20 nm to 1 um and are associated to the cytoplasmic membrane and its proximity. Resonant PTIR spectra of the cell surface are comparable to far-field IR spectra and provide the first selective measurement of the IR absorption spectrum of the cellular membrane of an intact cell. In resonant PTIR images, signal intensity, and therefore contrast, can be ascribed to a variety of factors, including mechanical, thermodynamic and spectroscopic properties of the cellular surface. While PTIR images are difficult to interpret in terms of spectroscopic absorption, they are easy to collect and provide unique contrast mechanisms without any exogenous labeling. As such they provide a new paradigm in cellular imaging and membrane biology and can be used to address a range of critical questions, from the nature of membrane lipid domains to viral membrane fusion.


title: The environmental dependence of rapidly-quenching and rejuvenating galaxies

abstract: By combining H$\alpha$ flux measurements from the Sloan Digital Sky Survey (SDSS) with UV flux observations from the Galaxy Evolution Explorer (GALEX), we examine the environmental dependence (through central/satellite distinction) of the rapid quenching and rejuvenation of galaxies. H$\alpha$ emissions trace the most massive stars, thereby indicating star-formation on timescales of $\sim 10$ Myr, while UV emission traces star-formation on timescales of $\sim 100$ Myr. These varying timescales are exploited to probe the most recent star-formation histories of galaxies. In this work, we define a class of transient galaxies which have UV emission typical of star-formation but negligible H$\alpha$ emission. We find that the occurrence of these transients has a strong stellar mass dependence in both the satellite and central population. However, while at stellar masses greater than $\sim 10^{10} M_\odot$ they occur with equal frequency regardless of environmental class, at lower stellar masses they are more common in satellites only, with an excess of about 1 percentage point across all low stellar mass galaxies. These satellite transients also have a strong halo mass and group-centric radial dependence suggesting they are driven by an environmental process. Finally, we select a sample of galaxies with H$\alpha$ emission but not UV emission which could contain short-timescale rejuvenating galaxies. These rejuvenating candidates are few in number and do not have a strong difference in their occurrence rate in centrals or satellites. These unique probes point to an environmental quenching mechanism which occurs on short timescales after the satellite has been in the group environment for a significant time - consistent with 'delayed-then-rapid' quenching.


title: Heterogeneous Optoelectronic Characteristics of Si Micropillar Arrays Fabricated by Metal-Assisted Chemical Etching

abstract: Recent progress achieved in metal-assisted chemical etching (MACE) has enabled the production of high-quality micropillar arrays for various optoelectronic applications. Si micropillars produced by MACE often show a porous Si/SiOx shell on crystalline pillar cores introduced by local electrochemical reactions. In this paper, we report the distinct optoelectronic characteristics of the porous Si/SiOx shell correlated to their chemical compositions. Local photoluminescent (PL) images obtained with an immersion oil objective lens in confocal microscopy show a red emission peak (about 650 nm) along the perimeter of the pillars that is threefold stronger compared to their center. On the basis of our analysis, we find an unexpected PL increase (about 540 nm) at the oil/shell interface. We suggest that both PL enhancements are mainly attributed to the porous structures, a similar behavior observed in previous MACE studies. Surface potential maps simultaneously recorded with topography reveal a significantly high surface potential on the sidewalls of MACE-synthesized pillars (+0.5 V), which is restored to the level of planar Si control (-0.5 V) after removing SiOx in hydrofluoric acid. These distinct optoelectronic characteristics of the Si/SiOx shell can be beneficial for various sensor architectures.


title: Adversarial Learning for Debiasing Knowledge Graph Embeddings

abstract: Knowledge Graphs (KG) are gaining increasing attention in both academia and industry. Despite their diverse benefits, recent research have identified social and cultural biases embedded in the representations learned from KGs. Such biases can have detrimental consequences on different population and minority groups as applications of KG begin to intersect and interact with social spheres. This paper aims at identifying and mitigating such biases in Knowledge Graph (KG) embeddings. As a first step, we explore popularity bias -- the relationship between node popularity and link prediction accuracy. In case of node2vec graph embeddings, we find that prediction accuracy of the embedding is negatively correlated with the degree of the node. However, in case of knowledge-graph embeddings (KGE), we observe an opposite trend. As a second step, we explore gender bias in KGE, and a careful examination of popular KGE algorithms suggest that sensitive attribute like the gender of a person can be predicted from the embedding. This implies that such biases in popular KGs is captured by the structural properties of the embedding. As a preliminary solution to debiasing KGs, we introduce a novel framework to filter out the sensitive attribute information from the KG embeddings, which we call FAN (Filtering Adversarial Network). We also suggest the applicability of FAN for debiasing other network embeddings which could be explored in future work.


title: Bosonic Bott Index and Disorder-Induced Topological Transitions of Magnons

abstract: We investigate the role of disorder on the various topological magnonic phases present in deformed honeycomb ferromagnets. To this end, we introduce a bosonic Bott index to characterize the topology of magnon spectra in finite, disordered systems. The consistency between the Bott index and Chern number is numerically established in the clean limit. We demonstrate that topologically protected magnon edge states are robust to moderate disorder and, as anticipated, localized in the strong regime. We predict a disorder-driven topological phase transition, a magnonic analog of the "topological Anderson insulator" in electronic systems, where the disorder is responsible for the emergence of the nontrivial topology. Combining the results for the Bott index and transport properties, we show that bulk-boundary correspondence holds for disordered topological magnons. Our results open the door for research on topological magnonics as well as other bosonic excitations in finite and disordered systems.


title: Controlled creation of a singular spinor vortex by circumventing the Dirac belt trick

abstract: Persistent topological defects and textures are particularly dramatic consequences of superfluidity. Among the most fascinating examples are the singular vortices arising from the rotational symmetry group SO(3), with surprising topological properties illustrated by Dirac's famous belt trick. Despite considerable interest, controlled preparation and detailed study of vortex lines with complex internal structure in fully three-dimensional spinor systems remains an outstanding experimental challenge. Here, we propose and implement a reproducible and controllable method for creating and detecting a singular SO(3) line vortex from the decay of a non-singular spin texture in a ferromagnetic spin-1 Bose--Einstein condensate. Our experiment explicitly demonstrates the SO(3) character and the unique spinor properties of the defect. Although the vortex is singular, its core fills with atoms in the topologically distinct polar magnetic phase. The resulting stable, coherent topological interface has analogues in systems ranging from condensed matter to cosmology and string theory.


title: Dynamic Knapsack Optimization Towards Efficient Multi-Channel Sequential Advertising

abstract: In E-commerce, advertising is essential for merchants to reach their target users. The typical objective is to maximize the advertiser's cumulative revenue over a period of time under a budget constraint. In real applications, an advertisement (ad) usually needs to be exposed to the same user multiple times until the user finally contributes revenue (e.g., places an order). However, existing advertising systems mainly focus on the immediate revenue with single ad exposures, ignoring the contribution of each exposure to the final conversion, thus usually falls into suboptimal solutions. In this paper, we formulate the sequential advertising strategy optimization as a dynamic knapsack problem. We propose a theoretically guaranteed bilevel optimization framework, which significantly reduces the solution space of the original optimization space while ensuring the solution quality. To improve the exploration efficiency of reinforcement learning, we also devise an effective action space reduction approach. Extensive offline and online experiments show the superior performance of our approaches over state-of-the-art baselines in terms of cumulative revenue.


title: Effective conductivity of the multidimensional chessboard

abstract: An algebraic formula for the effective conductivity of a d-dimensional, two-component chessboard (checkerboard) is proposed. The derivation relies on the self-duality of the square bond lattice, the principle of universality, and the analytical capabilities of the Walker Diffusion Method.


title: First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

abstract: Physical and chemical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.


title: First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations

abstract: Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.


title: First results from SMAUG: Uncovering the Origin of the Multiphase Circumgalactic Medium with a Comparative Analysis of Idealized and Cosmological Simulations

abstract: We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.


title: First results from SMAUG: The need for preventative stellar feedback and improved baryon cycling in semi-analytic models of galaxy formation

abstract: Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.


title: Learning and Planning in Average-Reward Markov Decision Processes

abstract: We introduce learning and planning algorithms for average-reward MDPs, including 1) the first general proven-convergent off-policy model-free control algorithm without reference states, 2) the first proven-convergent off-policy model-free prediction algorithm, and 3) the first off-policy learning algorithm that converges to the actual value function rather than to the value function plus an offset. All of our algorithms are based on using the temporal-difference error rather than the conventional error when updating the estimate of the average reward. Our proof techniques are a slight generalization of those by Abounadi, Bertsekas, and Borkar (2001). In experiments with an Access-Control Queuing Task, we show some of the difficulties that can arise when using methods that rely on reference states and argue that our new algorithms can be significantly easier to use.


title: Estimation and Decomposition of Rack Force for Driving on Uneven Roads

abstract: The force transmitted from the front tires to the steering rack of a vehicle, called the rack force, plays an important role in the function of electric power steering (EPS) systems. Estimates of rack force can be used by EPS to attenuate road feedback and reduce driver effort. Further, estimates of the components of rack force (arising, for example, due to steering angle and road profile) can be used to separately compensate for each component and thereby enhance steering feel. In this paper, we present three vehicle and tire model-based rack force estimators that utilize sensed steering angle and road profile to estimate total rack force and individual components of rack force. We test and compare the real-time performance of the estimators by performing driving experiments with non-aggressive and aggressive steering maneuvers on roads with low and high frequency profile variations. The results indicate that for aggressive maneuvers the estimators using non-linear tire models produce more accurate rack force estimates. Moreover, only the estimator that incorporates a semi-empirical Rigid Ring tire model is able to capture rack force variation for driving on a road with high frequency profile variation. Finally, we present results from a simulation study to validate the component-wise estimates of rack force.


title: Connected sums of sphere products and minimally non-Golod complexes

abstract: We show that if the moment-angle complex $\mathcal{Z}_K$ associated to a simplicial complex $K$ is homotopy equivalent to a connected sum of sphere products with two spheres in each product, then $K$ decomposes as the simplicial join of an $n$-simplex $\Delta^n$ and a minimally non-Golod complex. In particular, we prove that $K$ is minimally non-Golod for every moment-angle complex $\mathcal{Z}_K$ homeomorphic to a connected sum of two-fold products of spheres, answering a question of Grbi\'c, Panov, Theriault and Wu.


title: Mechanisms of dynamic near-wake modulation of a utility-scale wind turbine

abstract: The current study uses large eddy simulations to investigate the transient response of a utility-scale wind turbine wake to dynamic changes in atmospheric and operational conditions, as observed in previous field-scale measurements. Most wind turbine wake investigations assume quasi-steady conditions, but real wind turbines operate in a highly stochastic atmosphere, and their operation (e.g., blade pitch, yaw angle) changes constantly in response. Furthermore, dynamic control strategies have been recently proposed to optimize wind farm power generation and longevity. Therefore, improved understanding of dynamic wake behaviors is essential. First, changes in blade pitch are investigated and the wake expansion response is found to display hysteresis as a result of flow inertia. The timescales of the wake response to different pitch rates are quantified. Next, changes in wind direction with different timescales are explored. Under short timescales, the wake deflection is in the opposite direction of that observed under quasi-steady conditions. Finally, yaw changes are implemented at different rates, and the maximum inverse wake deflection and timescale are quantified, showing a clear dependence on yaw rate. To gain further physical understanding of the mechanism behind the inverse wake deflection, the streamwise vorticity in different parts of the wake is quantified. The results of this study provide guidance for the design of advanced wake flow control algorithms. The lag in wake response observed for both blade pitch and yaw changes shows that proposed dynamic control strategies must implement turbine operational changes with a timescale on the order of the rotor timescale or slower.


title: Scaling Symbolic Methods using Gradients for Neural Model Explanation

abstract: Symbolic techniques based on Satisfiability Modulo Theory (SMT) solvers have been proposed for analyzing and verifying neural network properties, but their usage has been fairly limited owing to their poor scalability with larger networks. In this work, we propose a technique for combining gradient-based methods with symbolic techniques to scale such analyses and demonstrate its application for model explanation. In particular, we apply this technique to identify minimal regions in an input that are most relevant for a neural network's prediction. Our approach uses gradient information (based on Integrated Gradients) to focus on a subset of neurons in the first layer, which allows our technique to scale to large networks. The corresponding SMT constraints encode the minimal input mask discovery problem such that after masking the input, the activations of the selected neurons are still above a threshold. After solving for the minimal masks, our approach scores the mask regions to generate a relative ordering of the features within the mask. This produces a saliency map which explains "where a model is looking" when making a prediction. We evaluate our technique on three datasets - MNIST, ImageNet, and Beer Reviews, and demonstrate both quantitatively and qualitatively that the regions generated by our approach are sparser and achieve higher saliency scores compared to the gradient-based methods alone. Code and examples are at - https://github.com/google-research/google-research/tree/master/smug_saliency


title: Recursion operators and the hierarchies of MKdV equations related to $D_4^{(1)}$, $D_4^{(2)}$ and $D_4^{(3)}$ Kac-Moody algebras

abstract: We constructed the three nonequivalent gradings in the algebra $D_4 \simeq so(8)$. The first one is the standard one obtained with the Coxeter automorphism $C_1=S_{\alpha_2} S_{\alpha_1}S_{\alpha_3}S_{\alpha_4}$ using its dihedral realization. In the second one we use $C_2 = C_1R$ where $R$ is the mirror automorphism. The third one is $C_3 = S_{\alpha_2}S_{\alpha_1}T$ where $T$ is the external automorphism of order 3. For each of these gradings we constructed the basis in the corresponding linear subspaces $\mathfrak{g}^{(k)}$, the orbits of the Coxeter automorphisms and the related Lax pairs generating the corresponding mKdV hierarchies. We found compact expressions for each of the hierarchies in terms of the recursion operators. At the end we wrote explicitly the first nontrivial mKdV equations and their Hamiltonians. For $D_4^{(1)}$ these are in fact two mKdV systems, due to the fact that in this case the exponent $3$ has multiplicity 2. Each of these mKdV systems consist of 4 equations of third order with respect to $\partial_x$. For $D_4^{(2)}$ this is a system of three equations of third order with respect to $\partial_x$. Finally, for $D_4^{(3)}$ this is a system of two equations of fifth order with respect to $\partial_x$.


title: Universal linguistic inductive biases via meta-learning

abstract: How do learners acquire languages from the limited data available to them? This process must involve some inductive biases - factors that affect how a learner generalizes - but it is unclear which inductive biases can explain observed patterns in language acquisition. To facilitate computational modeling aimed at addressing this question, we introduce a framework for giving particular linguistic inductive biases to a neural network model; such a model can then be used to empirically explore the effects of those inductive biases. This framework disentangles universal inductive biases, which are encoded in the initial values of a neural network's parameters, from non-universal factors, which the neural network must learn from data in a given language. The initial state that encodes the inductive biases is found with meta-learning, a technique through which a model discovers how to acquire new languages more easily via exposure to many possible languages. By controlling the properties of the languages that are used during meta-learning, we can control the inductive biases that meta-learning imparts. We demonstrate this framework with a case study based on syllable structure. First, we specify the inductive biases that we intend to give our model, and then we translate those inductive biases into a space of languages from which a model can meta-learn. Finally, using existing analysis techniques, we verify that our approach has imparted the linguistic inductive biases that it was intended to impart.


title: General decay and blow-up of solutions for a nonlinear wave equation with memory and fractional boundary damping terms

abstract: The paper studies the global existence and general decay of solutions using Lyaponov functional for a nonlinear wave equation, taking into account the fractional derivative boundary condition and memory term. In addition, we establish the blow up of solutions with nonpositive initial energy.


title: Impact of contamination factors on the COVID-19 evolution in Senegal

abstract: In this article, we perform an analysis of COVID-19 on one of the South Saharan countries (hot zone), the Senegal (West Africa). Many questions remain unanswered: why the African continent is not very contaminated compared to other continents. Factors of cross immunity, temperature, population density, youth, etc. are taken into account for an analysis of the contamination factors. Numerical simulations are carried out for a prediction over the coming week.


title: Link Prediction Using Supervised Machine Learning based on Aggregated and Topological Features

abstract: Link prediction is an important task in social network analysis. There are different characteristics (features) in a social network that can be used for link prediction. In this paper, we evaluate the effectiveness of aggregated features and topological features in link prediction using supervised learning. The aggregated features, in a social network, are some aggregation functions of the attributes of the nodes. Topological features describe the topology or structure of a social network, and its underlying graph. We evaluated the effectiveness of these features by measuring the performance of different supervised machine learning methods. Specifically, we selected five well-known supervised methods including J48 decision tree, multi-layer perceptron (MLP), support vector machine (SVM), logistic regression and Naive Bayes (NB). We measured the performance of these five methods with different sets of features of the DBLP Dataset. Our results indicate that the combination of aggregated and topological features generates the best performance. For evaluation purposes, we used accuracy, area under the ROC curve (AUC) and F-Measure. Our selected features can be used for the analysis of almost any social network. This is because these features provide the important characteristics of the underlying graph of the social networks. The significance of our work is that the selected features can be very effective in the analysis of big social networks. In such networks we usually deal with big data sets, with millions or billions of instances. Using fewer, but more effective, features can help us for the analysis of big social networks.


title: Veering triangulations and the Thurston norm: homology to isotopy

abstract: We show that a veering triangulation $\tau$ specifies a face $\sigma$ of the Thurston norm ball of a closed three-manifold, and computes the Thurston norm in the cone over $\sigma$. Further, we show that $\tau$ collates exactly the taut surfaces representing classes in the cone over $\sigma$ up to isotopy. The analysis includes nonlayered veering triangulations and nonfibered faces.


title: Simple halo model formalism for the cosmic infrared background and its correlation with the thermal Sunyaev Zel'dovich effect

abstract: Modelling the anisotropies in the CIB on all the scales is a challenging task due to the complex nature of the galaxy evolution and thus often requires too many parameters in order to fit the observational data. In this paper, we present a new halo model for the anisotropies of the CIB using only four parameters. Our model connects the mass accretion on the dark matter halos to the star formation rate. Using this model, we find that the halo mass with the maximum efficiency for converting the accreted baryons into stars is $\log_{10}M_\mathrm{max} = {12.94}^{+0.02}_{-0.02} \: M_\odot$, consistent with other studies. Accounting for the mass loss through stellar evolution, we find, for an intermediate age galaxy, that the star formation efficiency defined as $M_\star(z)/M_b(z)$ is equal to 0.19 and 0.21 at redshift 0.1 and 2 respectively, in good agreement with the values obtained by previous studies. It is the first time that a CIB model is used to fit Planck and Herschel CIB power spectra simultaneously. However we find that the large angular scale Planck and Herschel data are not fully compatible with the small scale Herschel data (for $\ell>3000$). The CIB is expected to be correlated with the tSZ signal of the galaxy clusters. Using this halo model for the CIB, and a halo model for the tSZ with a single parameter, we also provide a consistent framework to calculate this cross-correlation, which requires no extra parameter. The CIB$\times$tSZ correlation is found to be higher when inferred with a combination of two widely spaced frequency channels (e.g. 143x857 GHz). The CIB, tSZ, and CIB$\times$tSZ act as foregrounds while measuring the kSZ power spectrum from the CMB power spectrum and need to be removed. Due to its simplistic nature and low number of parameters, the halo model formalism presented here is quite useful for such an analysis to measure the kSZ power spectrum accurately.


title: Yetter-Drinfel'd algebras and coideals of Weak Hopf $C^*$-Algebras

abstract: We characterize braided commutative Yetter-Drinfeld $C^*$-algebras over weak Hopf $C^*$-algebras in categorical terms. Using this, we then study quotient type coideal subalgebras of a given weak Hopf $C^*$-algebra $\mathcal G$ and coideal subalgebras invariant with respect to the adjoint action of $\mathcal G$. Finally, as an example, we explicitly describe quotient type coideal subalgebras of the weak Hopf $C^*$-algebras associated with Tambara-Yamagami categories.


title: Asymmetric metric learning for knowledge transfer

abstract: Knowledge transfer from large teacher models to smaller student models has recently been studied for metric learning, focusing on fine-grained classification. In this work, focusing on instance-level image retrieval, we study an asymmetric testing task, where the database is represented by the teacher and queries by the student. Inspired by this task, we introduce asymmetric metric learning, a novel paradigm of using asymmetric representations at training. This acts as a simple combination of knowledge transfer with the original metric learning task. We systematically evaluate different teacher and student models, metric learning and knowledge transfer loss functions on the new asymmetric testing as well as the standard symmetric testing task, where database and queries are represented by the same model. We find that plain regression is surprisingly effective compared to more complex knowledge transfer mechanisms, working best in asymmetric testing. Interestingly, our asymmetric metric learning approach works best in symmetric testing, allowing the student to even outperform the teacher.


title: Massive photons and electrically charged black holes

abstract: The characteristic time of disappearance of electric field in massive electrodynamics during the capture of electric charge by black holes is calculated. It is shown that this time does not depend upon the photon mass. The electric field at large distances disappears with the speed of light.


title: Inference in Bayesian Additive Vector Autoregressive Tree Models

abstract: Vector autoregressive (VAR) models assume linearity between the endogenous variables and their lags. This assumption might be overly restrictive and could have a deleterious impact on forecasting accuracy. As a solution, we propose combining VAR with Bayesian additive regression tree (BART) models. The resulting Bayesian additive vector autoregressive tree (BAVART) model is capable of capturing arbitrary non-linear relations between the endogenous variables and the covariates without much input from the researcher. Since controlling for heteroscedasticity is key for producing precise density forecasts, our model allows for stochastic volatility in the errors. We apply our model to two datasets. The first application shows that the BAVART model yields highly competitive forecasts of the US term structure of interest rates. In a second application, we estimate our model using a moderately sized Eurozone dataset to investigate the dynamic effects of uncertainty on the economy.


title: A high precision falling-ball viscometer using a fast camera

abstract: This paper describes a simple and inexpensive method of measuring viscosity of a Newtonian fluid using the ball drop technique and an inexpensive point and shoot ~1000 frame per second camera. We successfully measured the viscosity of glycerol and glycerol-water mixture with high precision. We used three different size copper balls of diameters 0.8 mm, 1.59 mm, and 2.38 mm to check the accuracy of the measured viscosity in different concentrations of glycerol-water mixer solutions ranging from 50% to 100% (pure glycerol). Our measurements are in excellent agreement with the measurements conducted by other standard techniques. The simple and inexpensive techniques and physics we present in this manuscript can be employed to create a simple viscosity measurement setup for learning about complex fluid mechanics even at the undergraduate laboratory and high school teaching laboratory.


title: A Generative Neural Network Framework for Automated Software Testing

abstract: Search Based Software Testing (SBST) is a popular automated testing technique which uses a feedback mechanism to search for faults in software. Despite its popularity, it has fundamental challenges related to the design, construction and interpretation of the feedback. Neural Networks (NN) have been hugely popular in recent years for a wide range of tasks. We believe that they can address many of the issues inherent to common SBST approaches. Unfortunately, NNs require large and representative training datasets. In this work we present an SBST framework based on a deconvolutional generative neural network. Not only does it retain the beneficial qualities that make NNs appropriate for SBST tasks, it also produces its own training data which circumvents the problem of acquiring a training dataset that limits the use of NNs. We demonstrate through a series of experiments that this architecture is possible and practical. It generates diverse, sensible program inputs, while exploring the space of program behaviours. It also creates a meaningful ordering over program behaviours and is able to find crashing executions. This is all done without any prior knowledge of the program. We believe this proof of concept opens new directions for future work at the intersection of SBST and neural networks.


title: Learning Sparse Prototypes for Text Generation

abstract: Prototype-driven text generation uses non-parametric models that first choose from a library of sentence "prototypes" and then modify the prototype to generate the output text. While effective, these methods are inefficient at test time as a result of needing to store and index the entire training corpus. Further, existing methods often require heuristics to identify which prototypes to reference at training time. In this paper, we propose a novel generative model that automatically learns a sparse prototype support set that, nonetheless, achieves strong language modeling performance. This is achieved by (1) imposing a sparsity-inducing prior on the prototype selection distribution, and (2) utilizing amortized variational inference to learn a prototype retrieval function. In experiments, our model outperforms previous prototype-driven language models while achieving up to a 1000x memory reduction, as well as a 1000x speed-up at test time. More interestingly, we show that the learned prototypes are able to capture semantics and syntax at different granularity as we vary the sparsity of prototype selection, and that certain sentence attributes can be controlled by specifying the prototype for generation.


title: Nonlinear coupling in an asymmetric pendulum

abstract: We investigate the nonlinear effect of a pendulum with the upper end fixed to an elastic rod which is only allowed to vibrate horizontally. The pendulum will start rotating and trace a delicate stationary pattern when released without initial angular momentum. We explain it as amplitude modulation due to nonlinear coupling between the two degrees of freedom. Though the phenomenon of conversion between radial and azimuthal oscillations is common for asymmetric pendulums, nonlinear coupling between the two oscillations is usually overlooked. In this paper, we build a theoretical model and obtain the pendulum's equations of motion. The pendulum's motion patterns are solved numerically and analytically using the method of multiple scales. In the analytical solution, the modulation period not only depends on the dynamical parameters, but also on the pendulum's initial releasing positions, which is a typical nonlinear behavior. The analytical approximate solutions are supported by numerical results. This work provides a good demonstration as well as a research project of nonlinear dynamics on different levels from high school to undergraduate students.


title: The atomic structure of ancient grain boundaries

abstract: Democritus and the early atomists held that "the material cause of all things that exist is the coming together of atoms and void. Atoms are eternal and have many different shapes, and they can cluster together to create things that are perceivable. Differences in shape, arrangement, and position of atoms produce different phenomena". Like the atoms of Democritus, the Grim Reaper solution to curve shortening flow is eternal and indivisible -- it does not split off a line, and is itself its only "asymptotic translator". Confirming the heuristic described by Huisken and Sinestrari [J. Differential Geom. 101, 2 (2015), 267-287], we show that it gives rise to a great diversity of convex ancient and translating solutions to mean curvature flow, through the evolution of families of Grim hyperplanes in suitable configurations. We construct, in all dimensions $n\ge 2$, a large family of new examples, including both symmetric and asymmetric examples, as well as many eternal examples that do not evolve by translation. The latter resolve a conjecture of White [J. Amer. Math. Soc. 16, 1 (2003), 123-138]. We also provide a detailed asymptotic analysis of convex ancient solutions in slab regions in general. Roughly speaking, we show that they decompose "backwards in time" into a canonical configuration of Grim hyperplanes which satisfies certain necessary conditions. An analogous decomposition holds "forwards in time" for eternal solutions. One consequence is a new rigidity result for translators. Another is that, in dimension two, solutions are necessarily reflection symmetric across the mid-plane of their slab.


title: Parallel Betweenness Computation in Graph Database for Contingency Selection

abstract: Parallel betweenness computation algorithms are proposed and implemented in a graph database for power system contingency selection. Principles of the graph database and graph computing are investigated for both node and edge betweenness computation. Experiments on the 118-bus system and a real power system show that speed-up can be achieved for both node and edge betweenness computation while the speeding effect on the latter is more remarkable due to the data retrieving advantages of the graph database on the power network data.


title: A Nonmonotone Matrix-Free Algorithm for Nonlinear Equality-Constrained Least-Squares Problems

abstract: Least squares form one of the most prominent classes of optimization problems, with numerous applications in scientific computing and data fitting. When such formulations aim at modeling complex systems, the optimization process must account for nonlinear dynamics by incorporating constraints. In addition, these systems often incorporate a large number of variables, which increases the difficulty of the problem, and motivates the need for efficient algorithms amenable to large-scale implementations. In this paper, we propose and analyze a Levenberg-Marquardt algorithm for nonlinear least squares subject to nonlinear equality constraints. Our algorithm is based on inexact solves of linear least-squares problems, that only require Jacobian-vector products. Global convergence is guaranteed by the combination of a composite step approach and a nonmonotone step acceptance rule. We illustrate the performance of our method on several test cases from data assimilation and inverse problems: our algorithm is able to reach the vicinity of a solution from an arbitrary starting point, and can outperform the most natural alternatives for these classes of problems.


title: Handling Missing Data in Decision Trees: A Probabilistic Approach

abstract: Decision trees are a popular family of models due to their attractive properties such as interpretability and ability to handle heterogeneous data. Concurrently, missing data is a prevalent occurrence that hinders performance of machine learning models. As such, handling missing data in decision trees is a well studied problem. In this paper, we tackle this problem by taking a probabilistic approach. At deployment time, we use tractable density estimators to compute the "expected prediction" of our models. At learning time, we fine-tune parameters of already learned trees by minimizing their "expected prediction loss" w.r.t.\ our density estimators. We provide brief experiments showcasing effectiveness of our methods compared to few baselines.


title: Nuclear X-ray Activity in Low-Surface-Brightness Galaxies: Prospects for Constraining the Local Black Hole Occupation Fraction with a Chandra Successor Mission

abstract: About half of nearby galaxies have a central surface brightness >1 magnitude below that of the sky. The overall properties of these low-surface-brightness galaxies (LSBGs) remain understudied, and in particular we know very little about their massive black hole population. This gap must be closed to determine the frequency of massive black holes at z=0 as well as to understand their role in regulating galaxy evolution. Here we investigate the incidence and intensity of nuclear, accretion-powered X-ray emission in a sample of 32 nearby LSBGs with the Chandra X-ray Observatory. A nuclear X-ray source is detected in 4 galaxies (12.5%). Based on an X-ray binary contamination assessment technique developed for normal galaxies, we conclude that the detected X-ray nuclei indicate low-level accretion from massive black holes. The active fraction is consistent with that expected from the stellar mass distribution of the LSBGs, but not their total baryonic mass, when using a scaling relation from an unbiased X-ray survey of normal galaxies. This suggests that their black holes co-evolved with their stellar population. In addition, the apparent agreement nearly doubles the number of galaxies available within ~100 Mpc for which a measurement of nuclear activity can efficiently constrain the frequency of black holes as a function of stellar mass. We conclude by discussing the feasibility of measuring this occupation fraction to a few percent precision below 1e10 solar masses with high-resolution, wide-field X-ray missions currently under consideration.


title: Fourier DiffuserScope: Single-shot 3D Fourier light field microscopy with a diffuser

abstract: Light field microscopy (LFM) uses a microlens array (MLA) near the sensor plane of a microscope to achieve single-shot 3D imaging of a sample without any moving parts. Unfortunately, the 3D capability of LFM comes with a significant loss of lateral resolution at the focal plane. Placing the MLA near the pupil plane of the microscope, instead of the image plane, can mitigate the artifacts and provide an efficient forward model, at the expense of field-of-view (FOV). Here, we demonstrate improved resolution across a large volume with Fourier DiffuserScope, which uses a diffuser in the pupil plane to encode 3D information, then computationally reconstructs the volume by solving a sparsity-constrained inverse problem. Our diffuser consists of randomly placed microlenses with varying focal lengths; the random positions provide a larger FOV compared to a conventional MLA, and the diverse focal lengths improve the axial depth range. To predict system performance based on diffuser parameters, we for the first time establish a theoretical framework and design guidelines, which are verified by numerical simulations, then build an experimental system that achieves $< 3$ um lateral and $4$ um axial resolution over a $1000 \times 1000 \times 280$ um$^3$ volume. Our diffuser design outperforms the MLA used in LFM, providing more uniform resolution over a larger volume, both laterally and axially.


title: Material Recognition for Automated Progress Monitoring using Deep Learning Methods

abstract: Recent advancements in Artificial intelligence, especially deep learning, has changed many fields irreversibly by introducing state of the art methods for automation. Construction monitoring has not been an exception; as a part of construction monitoring systems, material classification and recognition have drawn the attention of deep learning and machine vision researchers. However, to create production-ready systems, there is still a long path to cover. Real-world problems such as varying illuminations and reaching acceptable accuracies need to be addressed in order to create robust systems. In this paper, we have addressed these issues and reached a state of the art performance, i.e., 97.35% accuracy rate for this task. Also, a new dataset containing 1231 images of 11 classes taken from several construction sites is gathered and publicly published to help other researchers in this field.


title: SeMPE: Secure Multi Path Execution Architecture for Removing Conditional Branch Side Channels

abstract: One of the most prevalent source of side channel vulnerabilities is the secret-dependent behavior of conditional branches (SDBCB). The state-of-the-art solution relies on Constant-Time Expressions, which require high programming effort and incur high performance overheads. In this paper, we propose SeMPE, an approach that relies on architecture support to eliminate SDBCB without requiring much programming effort while incurring low performance overheads. The key idea is that when a secret-dependent branch is encountered, the SeMPE microarchitecture fetches, executes, and commits both paths of the branch, preventing the adversary from inferring secret values from the branching behavior of the program. To enable that, SeMPE relies on an architecture that is capable of safely executing both branch paths sequentially. Through microbenchmarks and an evaluation of a real-world library, we show that SeMPE incurs near ideal execution time overheads, which is the sum of the execution time of all branch paths of secret-dependent branches. SeMPE outperforms code generated by FaCT, a constant-time expression language, by up to a factor of 18x.


title: Data-Driven and Online Estimation of Linear Sensitivity Distribution Factors: A Low-rank Approach

abstract: Estimation of sensitivity matrices in electrical transmission systems allows grid operators to evaluate in real-time how changes in power injections reflect into changes in power flows. In this paper, we propose a robust low-rank minimization approach to estimate sensitivity matrices based on measurements of power injections and power flows. An online proximal-gradient method is proposed to estimate sensitivities on-the-fly from real-time measurements. The proposed method obtains meaningful estimates with fewer measurements when the regression model is underdetermined, in contrast with existing methods based on least-squares approaches. In addition, our method can also identify faulty measurements and handle missing data. In this work, convergence results in terms of dynamic regret are presented. Numerical tests corroborate the effectiveness of the novel approach and the robustness of missing measurements and outliers.


title: Random nearest neighbor graphs: the translation invariant case

abstract: If $(\omega(e))$ is a family of random variables (weights) assigned to the edges of $\mathbb{Z}^d$, the nearest neighbor graph is the directed graph induced by all edges $\langle x,y \rangle$ such that $\omega(\{x,y\})$ is minimal among all neighbors $y$ of $x$. That is, each vertex points to its closest neighbor, if the weights are viewed as edge-lengths. Nanda-Newman introduced nearest neighbor graphs when the weights are i.i.d. and continuously distributed and proved that a.s., all components of the undirected version of the graph are finite. We study the case of translation invariant, distinct weights, and prove that nearest neighbor graphs do not contain doubly-infinite directed paths. In contrast to the i.i.d. case, we show that in this stationary case, the graphs can contain either one or two infinite components (but not more) in dimension two, and $k$ infinite components for any $k \in [1,\infty]$ in dimension $\geq 3$. The latter constructions use a general procedure to exhibit a certain class of directed graphs as nearest neighbor graphs with distinct weights, and thereby characterize all translation invariant nearest neighbor graphs. We also discuss relations to geodesic graphs from first-passage percolation and implications for the coalescing walk model of Chaika-Krishnan.


title: On Cohen-Macaulay Hopf monoids in species

abstract: We study Cohen-Macaulay Hopf monoids in the category of species. The goal is to apply techniques from topological combinatorics to the study of polynomial invariants arising from combinatorial Hopf algebras. Given a polynomial invariant arising from a linearized Hopf monoid, we show that under certain conditions it is the Hilbert polynomial of a relative simplicial complex. If the Hopf monoid is Cohen-Macaulay, we give necessary and sufficient conditions for the corresponding relative simplicial complex to be relatively Cohen-Macaulay, which implies that the polynomial has a nonnegative $h$-vector. We apply our results to the weak and strong chromatic polynomials of acyclic mixed graphs, and the order polynomial of a double poset.


title: New developer metrics: Are comments as crucial as code contributions?

abstract: Open-source code development has become widespread in recent years. As a result, open-source software platforms have also become popular, and millions of developers from diverse locations are able to contribute to the same projects. On these platforms, various knowledge about them is obtained from user activity. This information is used in the form of developer metrics to solve a variety of challenges. In this study, we proposed new developer metrics, including commenting and issue-related activity, that require less information. We concluded that commenting on any feature of a project can be equally as valuable as code contribution. In addition, besides the quantitative ones, metrics based on only the existence of the activity have been shown to offer also considerable results. We saw that issues were crucial in identifying user contributions. Even if a developer makes a contribution to only one issue on a project, the relation between the developer and the project is tight. The hit scores are relatively lower because of the sparsity problem of our dataset; even so, we believe that we have presented improvable and remarkable new developer metrics.


title: Functoriality of Quantum Resource Theory and Variable-Domain Modal Logic

abstract: Quantum resource theory is a cutting-edge tool used to study practical implementations of quantum mechanical principles under realistic operational constraints. It does this by modelling quantum systems as restricted classes of possible or permissible experimental operations. Modal logic provides a formal tool for studying possibility and impossibility is a completely general logical setting. Here, I show that quantum resource theories may be functorially translated into models of variable-domain S4 modal logic in a way that provides a new class of formal techniques for exploring quantum resource-theoretic problems. I then extend this functorial relationship to an injective one by adding structure to these logical models to reflect the convertibility preorder of resources in the underlying resource theory. I conclude by discussing how this viewpoint may be deployed concretely.


title: Almost Existence From the Feral Perspective and Some Questions

abstract: This paper is concerned with the "almost existence" phenomenon for periodic orbits of Hamiltonian dynamical systems. In particular, we recover this result in both some standard and some novel cases via feral curves and an adiabatic degeneration.


title: Cameron-Liebler line classes

abstract: New examples of Cameron-Liebler line classes in $\mathrm{PG}(3,q)$ are given with parameter $\frac{1}{2}(q^2 -1)$. These examples have been constructed for many odd values of $q$ using a computer search, by forming a union of line orbits from a cyclic collineation group acting on the space. While there are many equivalent characterizations of these objects, perhaps the most significant is that a set of lines $\mathcal{L}$ in $\mathrm{PG}(3,q)$ is a Cameron-Liebler line class with parameter $x$ if and only if every spread $\mathcal{S}$ of the space shares precisely $x$ lines with $\mathcal{L}$. These objects are related to generalizations of symmetric tactical decompositions of $\mathrm{PG}(3,q)$, as well as to subgroups of $\mathrm{P\Gamma L}(4,q)$ having equally many orbits on points and lines of $\mathrm{PG}(3,q)$. Furthermore, in some cases the line classes we construct are related to two-intersection sets in $\mathrm{AG}(2,q)$. Since there are very few known examples of these sets for $q$ odd, any new results in this direction are of particular interest.


title: Human Trust-based Feedback Control: Dynamically varying automation transparency to optimize human-machine interactions

abstract: Human trust in automation plays an essential role in interactions between humans and automation. While a lack of trust can lead to a human's disuse of automation, over-trust can result in a human trusting a faulty autonomous system which could have negative consequences for the human. Therefore, human trust should be calibrated to optimize human-machine interactions with respect to context-specific performance objectives. In this article, we present a probabilistic framework to model and calibrate a human's trust and workload dynamics during his/her interaction with an intelligent decision-aid system. This calibration is achieved by varying the automation's transparency---the amount and utility of information provided to the human. The parameterization of the model is conducted using behavioral data collected through human-subject experiments, and three feedback control policies are experimentally validated and compared against a non-adaptive decision-aid system. The results show that human-automation team performance can be optimized when the transparency is dynamically updated based on the proposed control policy. This framework is a first step toward widespread design and implementation of real-time adaptive automation for use in human-machine interactions.


title: RLWE/PLWE equivalence for totally real cyclotomic subextensions via quasi-Vandermonde matrices

abstract: We propose and justify a generalised approach to prove the polynomial reduction of the RLWE to the PLWE problem attached to the ring of integers of a monogenic number field. We prove such equivalence in the case of the maximal totally real subextension of the 4p-th cyclotomic field, with p arbitrary prime.


title: Gaia and Hubble unveil the kinematics of stellar populations in the Type II globular clusters {\omega} Centauri and M 22

abstract: The origin of multiple stellar populations in Globular Clusters (GCs) is one of the greatest mysteries of modern stellar astrophysics. N-body simulations suggest that the present-day dynamics of GC stars can constrain the events that occurred at high redshift and led to the formation of multiple populations. Here, we combine multi-band photometry from the Hubble Space Telescope (HST) and ground-based facilities with HST and Gaia Data Release 2 proper motions to investigate the spatial distributions and the motions in the plane of the sky of multiple populations in the type II GCs NGC 5139 ($\omega\,$Centauri) and NGC 6656 (M 22). We first analyzed stellar populations with different metallicities. Fe-poor and Fe-rich stars in M 22 share similar spatial distributions and rotation patterns and exhibit similar isotropic motions. Similarly, the two main populations with different iron abundance in $\omega\,$Centauri share similar ellipticities and rotation patterns. When analyzing different radial regions, we find that the rotation amplitude decreases from the center towards the external regions. Fe-poor and Fe-rich stars of $\omega\,$Centauri are radially anisotropic in the central region and show similar degrees of anisotropy. We also investigate the stellar populations with different light-element abundances and find that their N-rich stars exhibit higher ellipticity than N-poor stars. In $\omega\,$Centauri Centauri both stellar groups are radially anisotropic. Interestingly, N-rich, Fe-rich stars exhibit different rotation patterns than N-poor stars with similar metallicities. The stellar populations with different nitrogen of M 22 exhibit similar rotation patterns and isotropic motions. We discuss these findings in the context of the formation of multiple populations.


title: High-Fidelity Machine Learning Approximations of Large-Scale Optimal Power Flow

abstract: The AC Optimal Power Flow (AC-OPF) is a key building block in many power system applications. It determines generator setpoints at minimal cost that meet the power demands while satisfying the underlying physical and operational constraints. It is non-convex and NP-hard, and computationally challenging for large-scale power systems. Motivated by the increased stochasticity in generation schedules and increasing penetration of renewable sources, this paper explores a deep learning approach to deliver highly efficient and accurate approximations to the AC-OPF. In particular, the paper proposes an integration of deep neural networks and Lagrangian duality to capture the physical and operational constraints. The resulting model, called OPF-DNN, is evaluated on real case studies from the French transmission system, with up to 3,400 buses and 4,500 lines. Computational results show that OPF-DNN produces highly accurate AC-OPF approximations whose costs are within 0.01% of optimality. OPF-DNN generates, in milliseconds, solutions that capture the problem constraints with high fidelity.


title: Data integration in high dimension with multiple quantiles

abstract: This article deals with the analysis of high dimensional data that come from multiple sources (experiments) and thus have different possibly correlated responses, but share the same set of predictors. The measurements of the predictors may be different across experiments. We introduce a new regression approach with multiple quantiles to select those predictors that affect any of the responses at any quantile level and estimate the nonzero parameters. Our estimator is a minimizer of a penalized objective function, which aggregates the data from the different experiments. We establish model selection consistency and asymptotic normality of the estimator. In addition we present an information criterion, which can also be used for consistent model selection. Simulations and two data applications illustrate the advantages of our method, which takes the group structure induced by the predictors across experiments and quantile levels into account.