title: A MUSE Spectro-imaging Study of the Th 28 Jet: Precession in the Inner Jet

abstract: Context: Th 28 is a Classical T Tauri star in the Lupus 3 cloud which drives an extended bipolar jet. Previous studies of the inner jet identified signatures of rotation around the outflow axis, a key result for theories of jet launching. Thus this is an important source in which to investigate the poorly understood jet launching mechanism. We investigate the morphology and kinematics of the Th 28 micro-jets with the aim of characterizing their structure and outflow activity, using optical integral-field spectroscopy observations obtained with VLT/MUSE. We use spectro-imaging and position-velocity maps to investigate the kinematic and morphological features of the jet, and obtain a catalogue of emission lines in which the jet is visible. A Lucy-Richardson deconvolution procedure is used to differentiate the structure of the inner micro-jet region. Spatial profiles extracted perpendicular to the jet axis are fitted to investigate the jet width, opening angle and the evolution of the jet axis. We confirm the previously identified knot HHW$_{2}$ within the red-shifted jet and identify three additional knots in each lobe for the first time. We also find [O III]$\lambda$5007 emission from the blue-shifted micro-jet including the knot closest to the star. Proper motions for the innermost knots on each side are estimated and we show that new knots are ejected on an approximate timescale of 10-15 years. The jet axis centroids show a point-symmetric wiggle within the inner portion of both micro-jets indicating precession. We use the jet shape to measure a precession period of 8 years, with a half-opening angle < 0.6$^{\circ}$. This may provide an alternative explanation for the rotation signatures previously reported. We find the jet shape to be compatible with precession due to a brown dwarf companion orbiting at a separation $\leq$ 0.3 au.


title: LeanML: A Design Pattern To Slash Avoidable Wastes in Machine Learning Projects

abstract: We introduce the first application of the lean methodology to machine learning projects. Similar to lean startups and lean manufacturing, we argue that lean machine learning (LeanML) can drastically slash avoidable wastes in commercial machine learning projects, reduce the business risk in investing in machine learning capabilities and, in so doing, further democratize access to machine learning. The lean design pattern we propose in this paper is based on two realizations. First, it is possible to estimate the best performance one may achieve when predicting an outcome $y \in \mathcal{Y}$ using a given set of explanatory variables $x \in \mathcal{X}$, for a wide range of performance metrics, and without training any predictive model. Second, doing so is considerably easier, faster, and cheaper than learning the best predictive model. We derive formulae expressing the best $R^2$, MSE, classification accuracy, and log-likelihood per observation achievable when using $x$ to predict $y$ as a function of the mutual information $I\left(y; x\right)$, and possibly a measure of the variability of $y$ (e.g. its Shannon entropy in the case of classification accuracy, and its variance in the case regression MSE). We illustrate the efficacy of the LeanML design pattern on a wide range of regression and classification problems, synthetic and real-life.


title: DeformerNet: A Deep Learning Approach to 3D Deformable Object Manipulation

abstract: In this paper, we propose a novel approach to 3D deformable object manipulation leveraging a deep neural network called DeformerNet. Controlling the shape of a 3D object requires an effective state representation that can capture the full 3D geometry of the object. Current methods work around this problem by defining a set of feature points on the object or only deforming the object in 2D image space, which does not truly address the 3D shape control problem. Instead, we explicitly use 3D point clouds as the state representation and apply Convolutional Neural Network on point clouds to learn the 3D features. These features are then mapped to the robot end-effector's position using a fully-connected neural network. Once trained in an end-to-end fashion, DeformerNet directly maps the current point cloud of a deformable object, as well as a target point cloud shape, to the desired displacement in robot gripper position. In addition, we investigate the problem of predicting the manipulation point location given the initial and goal shape of the object.


title: Refined Policy Improvement Bounds for MDPs

abstract: The policy improvement bound on the difference of the discounted returns plays a crucial role in the theoretical justification of the trust-region policy optimization (TRPO) algorithm. The existing bound leads to a degenerate bound when the discount factor approaches one, making the applicability of TRPO and related algorithms questionable when the discount factor is close to one. We refine the results in \cite{Schulman2015, Achiam2017} and propose a novel bound that is "continuous" in the discount factor. In particular, our bound is applicable for MDPs with the long-run average rewards as well.


title: Introducing Introspective Transmission for Reflection Characterization in High Frame-Rate Ultrasound Imaging

abstract: In ultrasound imaging, most of the transmit and receive beamforming schemes assume a homogenous diffuse medium and are evaluated based on contrast, temporal and spatial resolutions. However, most medium are constituted by both diffuse and specular regions and the assumption of a homogeneous medium does not hold good in all cases. Eventhough, there are some adaptive beamforming approaches proposed in literature, they are mostly for receive beamforming. This study is aimed at investigating the relevance of transmission schemes in characterizing the diffuse and specular reflections, particularly at high frame rates. The transmit wavefront interaction behavior on the tissue interfaces for two high frame-rate transmission modalities, i.e. conventional synthetic transmit aperture imaging and multi-angle plane-wave imaging, are analyzed for multiple in-vitro and in-vivo radio-frequency datasets. Two novel visualization perspectives are proposed called contour isolines and directivity variance to understand the wave interaction with different tissue interfaces by considering the scalar and vector aspects of the reflected intensities respectively. We also rationalize the relevance of choosing the appropriate receive beamforming scheme according to the transmission modality through a comparison of delay and sum, filtered delay multiply and sum, minimum variance distortionless response, and specular beamforming algorithms. It is found that a synergistic blend of transmit and receive beamforming schemes adaptive to the tissue is inevitable to avoid any misdiagnosis.


title: An Optimal Polarization-Entanglement Source: Frequency-Converted SPDC with Degeneracy, Indistinguishability, and Ultra-High Purity Configurable Over a Large Spectral Range

abstract: Modeling and simulations of entanglement-based quantum-networking protocols commonly assume perfect entangled states. Some investigations have been performed which show how imperfections cause the efficiency of the protocols to rapidly deteriorate. For polarization-entangled states created by SPDC, the fundamentals of phase matching lead to a trade-off problem for the optimal properties of the state. We present a new SPDC method which circumvents the trade-off problem and allows one to obtain degeneracy, indistinguishability, and heralded-single-photon spectral purity greater than 99% for any target SPDC wavelength in the visible and NIR spectrum. Therefore, our method can in principle generate optimal polarization-entangled states for polarization-entanglement-based quantum-networking protocols.


title: A counterexample to a Gr\"obner approach for noetherianity of the twisted commutative algebra ${\rm Sym}({\rm Sym}^2(\mathbf{C}^\infty))$

abstract: We resolve an open question posed by the authors of arXiv:1501.06925v2 in 2015 concerning a Gr\"obner theoretic approach for the noetherianity of the twisted commutative algebra ${\rm Sym}({\rm Sym}^2(\mathbf{C}^\infty))$. We provide a negative answer to their question by producing an explicit antichain. In doing so, we establish a connection to well studied posets of graphs under the subgraph and induced subgraph relation. We then analyze this connection to suggest future paths of investigation, for example a deeper study of initial ideals in twisted commutative algebras.


title: Accounting for spatial confounding in epidemiological studies with individual-level exposures: An exposure-penalized spline approach

abstract: In the presence of unmeasured spatial confounding, spatial models may actually increase (rather than decrease) bias, leading to uncertainty as to how they should be applied in practice. We evaluated spatial modeling approaches through simulation and application to a big data electronic health record study. Whereas the risk of bias was high for purely spatial exposures (e.g., built environment), we found very limited potential for increased bias for individual-level exposures that cluster spatially (e.g., smoking status). We also proposed a novel exposure-penalized spline approach that selects the degree of spatial smoothing to explain spatial variability in the exposure. This approach appeared promising for efficiently reducing spatial confounding bias.


title: Simulating Quantum Mechanics with a $\theta$-term and an 't Hooft Anomaly on a Synthetic Dimension

abstract: A topological $\theta$-term in gauge theories, including quantum chromodynamics in 3+1 dimensions, gives rise to a sign problem that makes classical Monte Carlo simulations impractical. Quantum simulations are not subject to such sign problems and are a promising approach to studying these theories in the future. In the near term, it is interesting to study simpler models that retain some of the physical phenomena of interest and their implementation on quantum hardware. For example, dimensionally-reducing gauge theories on small spatial tori produces quantum mechanical models which, despite being relatively simple to solve, retain interesting vacuum and symmetry structures from the parent gauge theories. Here we consider quantum mechanical particle-on-a-circle models, related by dimensional reduction to the 1+1d Schwinger model, that possess a $\theta$-term and realize an 't Hooft anomaly or global inconsistency at $\theta = \pi$. These models also exhibit the related phenomena of spontaneous symmetry breaking and instanton-anti-instanton interference in real time. We propose an experimental scheme for the real-time simulation of a particle on a circle with a $\theta$-term and a $\mathbb{Z}_n$ potential using a synthetic dimension encoded in a Rydberg atom. Simulating the Rydberg atom with realistic experimental parameters, we demonstrate that the essential physics can be well-captured by the experiment, with expected behavior in the tunneling rate as a function of $\theta$. Similar phenomena and observables can also arise in more complex quantum mechanical models connected to higher-dimensional nonabelian gauge theories by dimensional reduction.

categories: quant-ph hep-lat hep-ph hep-th physics.atom-ph


title: A comparative study of stochastic and deep generative models for multisite precipitation synthesis

abstract: Future climate change scenarios are usually hypothesized using simulations from weather generators. However, there only a few works comparing and evaluating promising deep learning models for weather generation against classical approaches. This study shows preliminary results making such evaluations for the multisite precipitation synthesis task. We compared two open-source weather generators: IBMWeathergen (an extension of the Weathergen library) and RGeneratePrec, and two deep generative models: GAN and VAE, on a variety of metrics. Our preliminary results can serve as a guide for improving the design of deep learning architectures and algorithms for the multisite precipitation synthesis task.


title: Kpop: A kernel balancing approach for reducing specification assumptions in survey weighting

abstract: With the precipitous decline in response rates, researchers and pollsters have been left with highly non-representative samples, relying on constructed weights to make these samples representative of the desired target population. Though practitioners employ valuable expert knowledge to choose what variables, $X$ must be adjusted for, they rarely defend particular functional forms relating these variables to the response process or the outcome. Unfortunately, commonly-used calibration weights -- which make the weighted mean $X$ in the sample equal that of the population -- only ensure correct adjustment when the portion of the outcome and the response process left unexplained by linear functions of $X$ are independent. To alleviate this functional form dependency, we describe kernel balancing for population weighting (kpop). This approach replaces the design matrix $\mathbf{X}$ with a kernel matrix, $\mathbf{K}$ encoding high-order information about $\mathbf{X}$. Weights are then found to make the weighted average row of $\mathbf{K}$ among sampled units approximately equal that of the target population. This produces good calibration on a wide range of smooth functions of $X$, without relying on the user to decide which $X$ or what functions of them to include. We describe the method and illustrate it by application to polling data from the 2016 U.S. presidential election.


title: Indefinite Causal Orders from Superpositions in Time

abstract: Treating the time of an event as a quantum variable, we derive a scheme in which superpositions in time are used to perform operations in an indefinite causal order. We use some aspects of a recently developed space-time-symmetric formalism of events. We propose a specific implementation of the scheme and recover the Quantum SWITCH, where quantum operations are performed in an order which is entangled with the state of a control qubit. Our scheme does not rely on any exotic quantum gravitational effect, but instead on phenomena which are naturally fuzzy in time, such as the decay of an excited atom.


title: The Convergence Rates of Blockchain Mining Games: A Markovian Approach

abstract: Understanding the strategic behavior of miners in a blockchain is of great importance for its proper operation. A common model for mining games considers an infinite time horizon, with players optimizing asymptotic average objectives. Implicitly, this assumes that the asymptotic behaviors are realized at human-scale times, otherwise invalidating current models. We study the mining game utilizing Markov Decision Processes. Our approach allows us to describe the asymptotic behavior of the game in terms of the stationary distribution of the induced Markov chain. We focus on a model with two players under immediate release, assuming two different objectives: the (asymptotic) average reward per turn and the (asymptotic) percentage of obtained blocks. Using tools from Markov chain analysis, we show the existence of a strategy achieving slow mixing times, exponential in the policy parameters. This result emphasizes the imperative need to understand convergence rates in mining games, validating the standard models. Towards this end, we provide upper bounds for the mixing time of certain meaningful classes of strategies. This result yields criteria for establishing that long-term averaged functions are coherent as payoff functions. Moreover, by studying hitting times, we provide a criterion to validate the common simplification of considering finite states models. For both considered objectives functions, we provide explicit formulae depending on the stationary distribution of the underlying Markov chain. In particular, this shows that both mentioned objectives are not equivalent. Finally, we perform a market share case study in a particular regime of the game. More precisely, we show that an strategic player with a sufficiently large processing power can impose negative revenue on honest players.


title: A Multi-stage Stochastic Programming Model for Adaptive Biomass Processing Operation under Uncertainty

abstract: Variations of physical and chemical characteristics of biomass reduce equipment utilization and increase operational costs of biomass processing. Biomass processing facilities use sensors to monitor the changes in biomass characteristics. Integrating sensory data into the operational decisions in biomass processing will increase its flexibility to the changing biomass conditions. In this paper, we propose a multi-stage stochastic programming model that minimizes the expected operational costs by identifying the initial inventory level and creating an operational decision policy for equipment speed settings. These policies take the sensory information data and the current biomass inventory level as inputs to dynamically adjust inventory levels and equipment settings according to the changes in the biomass' characteristics. We ensure that a prescribed target reactor utilization is consistently achieved by penalizing the violation of the target reactor feeding rate. A case study is developed using real-world data collected at Idaho National Laboratory's biomass processing facility. We show the value of multi-stage stochastic programming from an extensive computational experiment. Our sensitivity analysis indicates that updating the infeed rate of the system, the processing speed of equipment, and bale sequencing based on the moisture level of biomass improves the processing rate of the reactor and reduces operating costs.


title: Entropy exchange and thermal fluctuations in the Jaynes-Cummings model

abstract: The time-dependence of the quantum entropy for a two-level atom interacting with a single-cavity mode is computed using the Jaynes-Cummings model, when the initial state of the radiation field is prepared in a thermal state with temperature fluctuations. In order to describe the out-of-equilibrium situation, the Super Statistics approximation is implemented so that the gamma and the multi-level distribution functions are used to introduce the inverse temperature fluctuations. In the case of the gamma distribution, paralleling the Tsallis non-additive formalism, the entropy for the system is computed with the $q$-logarithm prescription, and the impact of the initial state of the atom is also taken into account. The results show that, in the first distribution, the $q$-parameter (related to the thermal fluctuations) modifies the partial entropies appreciably. In contrast, the way the inverse temperatures are distributed in the second one may lead to changes in the entropy functions.


title: Linear subspaces of minimal codimension in hypersurfaces

abstract: Let $k$ be a perfect field and let $X\subset {\mathbb P}^N$ be a hypersurface of degree $d$ defined over $k$ and containing a linear subspace $L$ defined over an algebraic closure $\overline{k}$ with $\mathrm{codim}_{{\mathbb P}^N}L=r$. We show that $X$ contains a linear subspace $L_0$ defined over $k$ with $\mathrm{codim}_{{\mathbb P}^N}L\le dr$. We conjecture that the intersection of all linear subspaces (over $\overline{k}$) of minimal codimension $r$ contained in $X$, has codimension bounded above only in terms of $r$ and $d$. We prove this when either $d\le 3$ or $r\le 2$.


title: Discriminants et sommes de carr\'es

abstract: In this paper we generalize the notion of discriminant for symmetric matrices and some results about it to the case of symmetric spaces.


title: Isomorphisms and derivations of partial flag incidence algebras

abstract: Let $P$ and $Q$ be finite posets and $R$ a commutative unital ring. In the case where $R$ is indecomposable, we prove that the $R$-linear isomorphisms between partial flag incidence algebras $I^3(P,R)$ and $I^3(Q,R)$ are exactly those induced by poset isomorphisms between $P$ and $Q$. We also show that the $R$-linear derivations of $I^3(P,R)$ are trivial.


title: Robust Risk-Sensitive Reinforcement Learning Agents for Trading Markets

abstract: Trading markets represent a real-world financial application to deploy reinforcement learning agents, however, they carry hard fundamental challenges such as high variance and costly exploration. Moreover, markets are inherently a multiagent domain composed of many actors taking actions and changing the environment. To tackle these type of scenarios agents need to exhibit certain characteristics such as risk-awareness, robustness to perturbations and low learning variance. We take those as building blocks and propose a family of four algorithms. First, we contribute with two algorithms that use risk-averse objective functions and variance reduction techniques. Then, we augment the framework to multi-agent learning and assume an adversary which can take over and perturb the learning process. Our third and fourth algorithms perform well under this setting and balance theoretical guarantees with practical use. Additionally, we consider the multi-agent nature of the environment and our work is the first one extending empirical game theory analysis for multi-agent learning by considering risk-sensitive payoffs.


title: On Some Properties of Irrational Subspaces

abstract: In this paper we discuss some properties of completely irrational subspaces. We prove that there exist completely irrational subspaces that are badly approximable and, moreover, sets of such subspaces are winning in different senses. We get some bounds for Diophantine exponents of vectors that lie in badly approximable subspaces that are completely irrational; in particular, for any vector $\xi$ from two-dimensional badly approximable completely irrational subspace of $\mathbb{R}^d$ one has $\hat{\omega}(\xi) \leq \frac{\sqrt{5} - 1}{2}$. Besides that, some statements about the dimension of subspaces generated by best approximations to completely irrational subspace easily follow from properties that we discuss.


title: Almost invariant subspaces and operators

abstract: We give an elementary proof of an efficient version of the Wagner's theorem on almost invariant subspaces and deduce some consequences in the context of Galois extensions.


title: An Information-state based Approach to the Optimal Output Feedback Control of Nonlinear Systems

abstract: This paper develops a data-based approach to the closed-loop output feedback control of nonlinear dynamical systems with a partial nonlinear observation model. We propose an information state based approach to rigorously transform the partially observed problem into a fully observed problem where the information state consists of the past several observations and control inputs. We further show the equivalence of the transformed and the initial partially observed optimal control problems and provide the conditions to solve for the deterministic optimal solution. We develop a data based generalization of the iterative Linear Quadratic Regulator (iLQR) to partially observed systems using a local linear time varying model of the information state dynamics approximated by an Autoregressive moving average (ARMA) model, that is generated using only the input-output data. This open-loop trajectory optimization solution is then used to design a local feedback control law, and the composite law then provides an optimum solution to the partially observed feedback design problem. The efficacy of the developed method is shown by controlling complex high dimensional nonlinear dynamical systems in the presence of model and sensing uncertainty.


title: Calibration of the liquid argon ionization response to low energy electronic and nuclear recoils with DarkSide-50

abstract: DarkSide-50 has demonstrated the high potential of dual-phase liquid argon time projection chambers in exploring interactions of WIMPs in the GeV/c$^2$ mass range. The technique, based on the detection of the ionization signal amplified via electroluminescence in the gas phase, allows to explore recoil energies down to the sub-keV range. We report here on the DarkSide-50 measurement of the ionization yield of electronic recoils down to $\sim$180~eV$_{er}$, exploiting $^{37}$Ar and $^{39}$Ar decays, and extrapolated to a few ionization electrons with the Thomas-Imel box model. Moreover, we present a model-dependent determination of the ionization response to nuclear recoils down to $\sim$500~eV$_{nr}$, the lowest ever achieved in liquid argon, using \textit{in situ} neutron calibration sources and external datasets from neutron beam experiments.


title: Tracing birth properties of stars with abundance clustering

abstract: To understand the formation and evolution of the Milky Way disk, we must connect its current properties to its past. We explore hydrodynamical cosmological simulations to investigate how the chemical abundances of stars might be linked to their origins. Using hierarchical clustering of abundance measurements in two Milky Way-like simulations with distributed and steady star formation histories, we find that abundance clusters of stars comprise different groups in birth place ($R_\text{birth}$) and time (age). Simulating observational abundance errors (0.05 dex), we find that to trace discrete groups of ($R_\text{birth}$, age) requires a large vector of abundances. Using 15-element abundances (Fe, O, Mg, S, Si, C, P, Mn, Ne, Al, N, V, Ba, Cr, Co), up to $\approx$ 10 clusters can be defined with $\approx$ 25% overlap in ($R_\text{birth}$, age). We build a simple model to show that it is possible to infer a star's age and $R_\text{birth}$ from abundances with precisions of $\pm$0.06 Gyr and $\pm$1.17 kpc respectively. We find that abundance clustering is ineffective for a third simulation, where low-$\alpha$ stars form distributed in the disc and early high-$\alpha$ stars form more rapidly in clumps that sink towards the galactic center as their constituent stars evolve to enrich the interstellar medium. However, this formation path leads to large age-dispersions across the [$\alpha$/Fe]-[Fe/H] plane, which is inconsistent with the Milky Way's observed properties. We conclude that abundance clustering is a promising approach toward charting the history of our Galaxy.


title: Non-Parametric Manifold Learning

abstract: We introduce an estimator for distances in a compact Riemannian manifold based on graph Laplacian estimates of the Laplace-Beltrami operator. We upper bound the error in the estimate of manifold distances, or more precisely an estimate of a spectrally truncated variant of manifold distance of interest in non-commutative geometry (cf. [Connes and Suijelekom, 2020]), in terms of spectral errors in the graph Laplacian estimates and, implicitly, several geometric properties of the manifold. A consequence is a proof of consistency for (untruncated) manifold distances. The estimator resembles, and in fact its convergence properties are derived from, a special case of the Kontorovic dual reformulation of Wasserstein distance known as Connes' Distance Formula.


title: Near-Optimal Algorithms for Linear Algebra in the Current Matrix Multiplication Time

abstract: In the numerical linear algebra community, it was suggested that to obtain nearly optimal bounds for various problems such as rank computation, finding a maximal linearly independent subset of columns (a basis), regression, or low-rank approximation, a natural way would be to resolve the main open question of Nelson and Nguyen (FOCS, 2013). This question is regarding the logarithmic factors in the sketching dimension of existing oblivious subspace embeddings that achieve constant-factor approximation. We show how to bypass this question using a refined sketching technique, and obtain optimal or nearly optimal bounds for these problems. A key technique we use is an explicit mapping of Indyk based on uncertainty principles and extractors, which after first applying known oblivious subspace embeddings, allows us to quickly spread out the mass of the vector so that sampling is now effective. We thereby avoid a logarithmic factor in the sketching dimension that is standard in bounds proven using the matrix Chernoff inequality. For the fundamental problems of rank computation and finding a basis, our algorithms improve Cheung, Kwok, and Lau (JACM, 2013), and are optimal to within a constant factor and a poly(log log(n))-factor, respectively. Further, for constant-factor regression and low-rank approximation we give the first optimal algorithms, for the current matrix multiplication exponent.


title: A Comparison of Methods for OOV-word Recognition on a New Public Dataset

abstract: A common problem for automatic speech recognition systems is how to recognize words that they did not see during training. Currently there is no established method of evaluating different techniques for tackling this problem. We propose using the CommonVoice dataset to create test sets for multiple languages which have a high out-of-vocabulary (OOV) ratio relative to a training set and release a new tool for calculating relevant performance metrics. We then evaluate, within the context of a hybrid ASR system, how much better subword models are at recognizing OOVs, and how much benefit one can get from incorporating OOV-word information into an existing system by modifying WFSTs. Additionally, we propose a new method for modifying a subword-based language model so as to better recognize OOV-words. We showcase very large improvements in OOV-word recognition and make both the data and code available.


title: Wolf 503 b: Characterization of a Sub-Neptune Orbiting a Metal-Poor K Dwarf

abstract: Using radial velocity measurements from four instruments, we report the mass and density of a $2.043\pm0.069 ~\rm{R}_{\oplus}$ sub-Neptune orbiting the quiet K-dwarf Wolf 503 (HIP 67285). In addition, we present improved orbital and transit parameters by analyzing previously unused short-cadence $K2$ campaign 17 photometry and conduct a joint radial velocity-transit fit to constrain the eccentricity at $0.41\pm0.05$. The addition of a transit observation by $Spitzer$ also allows us to refine the orbital ephemeris in anticipation of further follow-up. Our mass determination, $6.26^{+0.69}_{-0.70}~\rm{M}_{\odot}$, in combination with the updated radius measurements, gives Wolf 503 b a bulk density of $\rho = 2.92\pm ^{+0.50}_{-0.44}$ $\rm{g}~\rm{cm}^{-3}$. Using interior composition models, we find this density is consistent with an Earth-like core with either a substantial $\rm{H}_2\rm{O}$ mass fraction ($45^{+19.12}_{-16.15}\%$) or a modest H/He envelope ($0.5\pm0.28\%$). The low H/He mass fraction, along with the old age of Wolf 503 ($11\pm2$ Gyrs), makes this sub-Neptune an opportune subject for testing theories of XUV-driven mass loss while the brightness of its host ($J=8.3$ mag) makes it an attractive target for transmission spectroscopy.


title: Particle Merging-and-Splitting

abstract: Robustly handling collisions between individual particles in a large particle-based simulation has been a challenging problem. We introduce particle merging-and-splitting, a simple scheme for robustly handling collisions between particles that prevents inter-penetrations of separate objects without introducing numerical instabilities. This scheme merges colliding particles at the beginning of the time-step and then splits them at the end of the time-step. Thus, collisions last for the duration of a time-step, allowing neighboring particles of the colliding particles to influence each other. We show that our merging-and-splitting method is effective in robustly handling collisions and avoiding penetrations in particle-based simulations. We also show how our merging-and-splitting approach can be used for coupling different simulation systems using different and otherwise incompatible integrators. We present simulation tests involving complex solid-fluid interactions, including solid fractures generated by fluid interactions.


title: LAORAM: A Look Ahead ORAM Architecture for Training Large Embedding Tables

abstract: Data confidentiality is becoming a significant concern, especially in the cloud computing era. Memory access patterns have been demonstrated to leak critical information such as security keys and a program's spatial and temporal information. This information leak poses an even more significant privacy challenge in machine learning models with embedding tables. Embedding tables are routinely used to learn categorical features from training data. Even knowing the locations of the embedding table entries accessed, not the data within the embedding table, will compromise categorical input data to the model. Embedding entries are privacy-sensitive since they disclose valuable properties about the user. Oblivious RAM (ORAM), and its enhanced variants such as PathORAM have emerged as viable solutions to hide leakage from memory access streams. In this work, we present LAORAM, an ORAM framework explicitly designed to protect user privacy during embedding table training. LAORAM exploits the unique property of training, the training samples used in the future are known beforehand. LAORAM preprocesses the training samples to identify the memory blocks which are accessed together in the near future. The system tries to assign these blocks to as few paths as possible within the PathORAM infrastructure. LAORAM does this operation by combining multiple blocks accessed together as superblocks. To further increase performance, LAORAM uses a fat-tree structure for PathORAM reducing the number of background evictions required, which improves the stash usage. We have evaluated LAORAM using both a recommendation model (DLRM) and a NLP model (XLM-R) embedding table configurations. LAORAM performs 5 times faster than PathORAM on a recommendation dataset (Kaggle) and 5.4x faster on a NLP dataset (XNLI), while guaranteeing the same security guarantees as the original PathORAM.


title: Multiphase CO2 Dispersions in Microfluidics: Formation, Phases, and Mass Transfer

abstract: The dissolution and microfluidic mass transfer of carbon dioxide in water at high-pressure conditions are crucial for a myriad of technological applications, including microreactors, extractions, and carbon capture, utilization, and sequestration (CCUS) processes. In this experimental work, we use a high-pressure microfluidic method to elucidate the mass transfer process of CO2 in water at high pressure. An intriguing multiphase CO2 flow and dispersions are observed when operating at the pressure-temperature ($P$-$T$) condition close to the CO2 gas-liquid phase boundary ($P=6.5$ MPa and $T=23.5\pm 0.5~^{\circ}$C). We propose a series of strategies to unravel this complex multi-phase dynamics by calculating each phase's volume and mass change in a gas-liquid coexistent CO2 dispersion, estimating the possible CO2 concentration change in water, and comparing with the CO2 solubility data. Finally, we quantify the CO2 mass transfer by directly calculating the CO2 dissolution rate in water and estimating the volumetric mass transfer coefficient ($k_{L}a$). The results show that the mass transfer may be influenced by the specific area ($a$), CO2 concentration gradient in the water slug, and the traveling speed of a dispersion.


title: Learning to Limit Data Collection via Scaling Laws: A Computational Interpretation for the Legal Principle of Data Minimization

abstract: Modern machine learning systems are increasingly characterized by extensive personal data collection, despite the diminishing returns and increasing societal costs of such practices. Yet, data minimisation is one of the core data protection principles enshrined in the European Union's General Data Protection Regulation ('GDPR') and requires that only personal data that is adequate, relevant and limited to what is necessary is processed. However, the principle has seen limited adoption due to the lack of technical interpretation. In this work, we build on literature in machine learning and law to propose FIDO, a Framework for Inhibiting Data Overcollection. FIDO learns to limit data collection based on an interpretation of data minimization tied to system performance. Concretely, FIDO provides a data collection stopping criterion by iteratively updating an estimate of the performance curve, or the relationship between dataset size and performance, as data is acquired. FIDO estimates the performance curve via a piecewise power law technique that models distinct phases of an algorithm's performance throughout data collection separately. Empirical experiments show that the framework produces accurate performance curves and data collection stopping criteria across datasets and feature acquisition algorithms. We further demonstrate that many other families of curves systematically overestimate the return on additional data. Results and analysis from our investigation offer deeper insights into the relevant considerations when designing a data minimization framework, including the impacts of active feature acquisition on individual users and the feasability of user-specific data minimization. We conclude with practical recommendations for the implementation of data minimization.


title: Accurate Determination of Hubble Attenuation and Amplification in Expanding and Contracting Cold-Atom Universes

abstract: In the expanding universe, relativistic scalar fields are thought to be attenuated by "Hubble friction", which results from the dilation of the underlying spacetime metric. By contrast, in a contracting universe this pseudo-friction would lead to amplification. Here, we experimentally measure with five-fold better accuracy, both Hubble attenuation and amplification in expanding and contracting toroidally-shaped Bose-Einstein condensates, in which phonons are analogous to cosmological scalar fields. We find that the observed attenuation or amplification depends on the temporal phase of the phonon field, which is only possible for non-adiabatic dynamics. The measured strength of the Hubble friction disagrees with recent theory [J. M. Gomez Llorente and J. Plata, {\it Phys. Rev. A} {\bf 100} 043613 (2019) and S. Eckel and T. Jacobson, {\it SciPost Phys.} {\bf 10} 64 (2021)]; because our experiment probes physics outside the scope of this theory -- with large excitations in rings of intermediate thickness -- this indicates the presence of new physics.


title: End-to-End Tests of the Sensitivity of IceCube to the Neutrino Burst from a Core-Collapse Supernova

abstract: The next galactic supernova presents a once-in-a-lifetime opportunity to obtain detailed information about the explosion of a star and the extreme conditions found within its core. A core-collapse supernova will produce a neutrino burst visible up to half a day before the electromagnetic radiation from the explosion, so the burst will provide an early warning for optical follow-up. Since local supernovae are exceedingly rare, it is critical that neutrino detectors provide prompt alerts after the arrival of a burst. The IceCube Neutrino Observatory operates with $>99\%$ uptime and is sensitive to a variety of supernova models at levels $>10\sigma$ within the Milky Way. IceCube will issue supernova alerts in real time. IceCube's high sensitivity to supernovae, near perfect uptime, and ability to issue prompt alerts makes it a critical component of the worldwide network of detectors known as the SuperNova Early Warning System (SNEWS 2.0). A "Fire Drill" system was designed to inject simulated supernova signals into the IceCube online system. We will discuss IceCube's sensitivity to supernovae near the Milky Way, and describe the data challenges used to ensure the readiness of IceCube and its operators. We will also discuss the coordination of IceCube alerts and data challenges with SNEWS 2.0.


title: Investigating the side-chain structural organization behind the stability of protein folding and binding

abstract: What are the molecular mechanisms that dictate protein-protein binding stability and whether those are related to the ones behind protein fold stability are still largely open questions. Indeed, despite many past efforts, we still lack definitive models to account for experimental quantities like protein melting temperature or complex binding affinity. Here, we investigate and compare chemical and physical features on a dataset of protein with known melting temperature as well as a large dataset of protein-protein complexes with reliable experimental binding affinity. In particular, we probed the aminoacid composition and the organization of the network of intramolecular and intermolecular interaction energies among residues. We found that hydrophobic residues present on the protein surfaces are preferentially located in the binding regions, while charged residues behave oppositely. In addition, the abundance of polar amino acid like Serine and Proline correlates with the binding affinity of the complexes. Analysing the interaction energies we found that distant Coulombic interactions are responsible for thermal stability while the total inter-molecular van der Waals energy correlates with protein-protein binding affinity.


title: Optimality Conditions for Bilevel Imaging Learning Problems with Total Variation Regularization

abstract: We address the problem of optimal scale-dependent parameter learning in total variation image denoising. Such problems are formulated as bilevel optimization instances with total variation denoising problems as lower-level constraints. For the bilevel problem, we are able to derive M-stationarity conditions, after characterizing the corresponding Mordukhovich generalized normal cone and verifying suitable constraint qualification conditions. We also derive B-stationarity conditions, after investigating the Lipschitz continuity and directional differentiability of the lower-level solution operator. A characterization of the Bouligand subdifferential of the solution mapping, by means of a properly defined linear system, is provided as well. Based on this characterization, we propose a two-phase non-smooth trust-region algorithm for the numerical solution of the bilevel problem and test it computationally for two particular experimental settings.


title: Carleson estimates for the Green function on domains with lower dimensional boundaries

abstract: In the present paper, we consider an elliptic divergence form operator in $\mathbb{R}^n\setminus\mathbb{R}^d$ with $d<n-1$ and prove that its Green function is almost affine, in the sense that the normalized difference between the Green function with a sufficiently far away pole and a suitable affine function at every scale satisfies a Carleson measure estimate. The coefficients of the operator can be very oscillatory, and only need to satisfy some condition similar to the traditional quadratic Carleson condition.


title: On the complexity of open shop scheduling with time lags

abstract: The minimization of makespan in open shop with time lags has been shown NP-hard in the strong sense even for the case of two machines and unit-time operations. The minimization of total completion time however has remained open for that case though it has been shown NP-hard in the strong sense for weighted total completion time or for jobs with release dates. This note shows that the minimization of total completion time for two machines and unit-time operations is NP-hard in the strong sense which answers the long standing open question.


title: Signal reflection, quantum non-locality, and delayed choice experiments

abstract: Quantum nonlocality which is conventionally invoked for describing a composite entangled system is shown here to be a possible important characteristic of a single quantum object. To this end, we analyze some interactions of a single photon released from Fabry Perot resonator with environment. The split photon state with oppositely moving parts is shown to obey quantum nonlocality despite the sharp edges truncating each part. Photon post release reflection from a plane mirror is considered. The changing shape of the form during reflection contains moving discontinuities in electric and magnetic components of the pulse. They originate from preexisting edges of the form and move together, first away from and then back to the mirror. At the end of the process, the pulse restores its original shape, with electric component reversed. Altogether, the process demonstrates conservation of moving discontinuities. The considered experimental setup may be used for some new versions of a delayed choice experiment, with various options for insertion of detectors and the respective time delays. In all cases, the delayed insertion does not have any retroactive effect on the process. Key words: Evolution operator, signaling, reflection, discontinuity, quantum nonlocality


title: Physics-constrained intraventricular vector flow mapping by color Doppler

abstract: Color Doppler by transthoracic echocardiography creates 2-D fan-shaped maps of blood velocities in the cardiac cavities. It is a one-component velocimetric technique since it only returns the velocity components parallel to the ultrasound beams. Intraventricular vector flow mapping (iVFM) is a method to recover the blood velocity vectors from the Doppler scalar fields in an echocardiographic three-chamber view. We improved our iVFM numerical scheme by imposing physical constraints. The iVFM consisted in minimizing regularized Doppler residuals subject to the condition that two fluid-dynamics constraints were satisfied, namely planar mass conservation, and free-slip boundary conditions. The optimization problem was solved by using the Lagrange multiplier method. A finite-difference discretization of the optimization problem, written in the polar coordinate system centered on the cardiac ultrasound probe, led to a sparse linear system. The single regularization parameter was determined automatically for non-supervision considerations. The physics-constrained method was validated using realistic intracardiac flow data from a patient-specific CFD model. The numerical evaluations showed that the iVFM-derived velocity vectors were in very good agreement with the CFD-based original velocities, with relative errors ranged between 0.3 and 12%. We calculated two macroscopic measures of flow in the cardiac region of interest, the mean vorticity and mean stream function, and observed an excellent concordance between physics-constrained iVFM and CFD. The capability of physics-constrained iVFM was finally tested with in vivo color Doppler data acquired in patients routinely examined in the echocardiographic laboratory. The vortex that forms during the rapid filling was deciphered. The physics-constrained iVFM algorithm is ready for pilot clinical studies.


title: Irrotational Stokes flow in three system of coordinates

abstract: Irrotational flow is described with the second order elliptic partial differential equation $E^2\psi=0,$ where $\psi$ is the function to be derived and $E^2$ is the Stokes operator. In the present paper we derive the solution of $E^2\psi=0$ in three axisymmetric system of coordinates: the parabolic, the tangent sphere and the cardioid. We prove that Stokes equation separates variables in the parabolic coordinate system, but it R-separates variables in the other two. At the end of this manuscript we summarize all the known results for the 0-eigenspace of $E^2$ in axisymmetric system of coordinates.


title: Local H\"older regularity of minimizers for nonlocal variational problems

abstract: We study the regularity of solutions to a nonlocal variational problem, which is related to the image denoising model, and we show that, in two dimensions, minimizers have the same H\"older regularity as the original image. More precisely, if the datum is (locally) $\beta$-H\"older continuous for some $\beta \in (1-s,\,1]$, where $s \in (0,1)$ is a parameter related to the nonlocal operator, we prove that the solution is also $\beta$-H\"older continuous.


title: Generic projections of the $\mathrm{H}_{4}$ configuration of points

abstract: In the present note we give an example of a finite set of points in $\mathbb{P}^{3}$ which has the so-called geproci property, but it is neither a grid nor a half-grid. This answers a question on the existence of such sets rised by Pokora, Szemberg and Szpond [arXiv:2010.08863].


title: Codimension four regularity of generalized Einstein structures

abstract: We establish codimension 4 regularity of noncollapsed sequences of metrics with bounds on natural generalizations of the Ricci tensor. We obtain a priori L2 curvature estimates on such spaces, with diffeomorphism finiteness results and rigidity theorems as corollaries.


title: Automatic Fatou Property of Law-invariant Risk Measures

abstract: In the paper we investigate automatic Fatou property of law-invariant risk measures on a rearrangement-invariant function space $\mathcal{X}$ other than $L^\infty$. The main result is the following characterization: Every real-valued, law-invariant, coherent risk measure on $\mathcal{X}$ has the Fatou property at every random variable $X\in \mathcal{X}$ whose negative tails have vanishing norm (i.e., $\lim_n\|X\mathbf{1}_{\{X\leq -n\}}\|=0$) if and only if $\mathcal{X}$ satisfies the Almost Order Continuous Equidistributional Average (AOCEA) property, namely, $\mathrm{d}(\mathcal{CL}(X),\mathcal{X}_a) =0$ for any $X\in \mathcal{X}_+$, where $ \mathcal{CL}(X)$ is the convex hull of all random variables having the same distribution as $X$ and $\mathcal{X}_a=\{X\in\mathcal{X}:\lim_n \|X\mathbf{1}_{ \{|X|\geq n\} }\| =0\}$. As a consequence, we show that under the AOCEA property, every real-valued, law-invariant, coherent risk measure on $\mathcal{X}$ admits a tractable dual representation at every $X\in \mathcal{X}$ whose negative tails have vanishing norm. Furthermore, we show that the AOCEA property is satisfied by most classical model spaces, including Orlicz spaces, and therefore the foregoing results have wide applications.


title: Some rigidity results for the Hawking mass and a lower bound for the Bartnik capacity

abstract: We prove rigidity results involving the Hawking mass for surfaces immersed in a $3$-dimensional, complete Riemannian manifold $(M,g)$ with non-negative scalar curvature (resp. with scalar curvature bounded below by $-6$). Roughly, the main result states that if an open subset $\Omega\subset M$ satisfies that every point has a neighbourhood $U\subset \Omega$ such that the supremum of the Hawking mass of surfaces contained in $U$ is non-positive, then $\Omega$ is locally isometric to Euclidean ${\mathbb R}^3$ (resp. locally isometric to the Hyperbolic 3-space ${\mathbb H}^3$). Under mild asymptotic conditions on the manifold $(M,g)$ (which encompass as special cases the standard "asymptotically flat" or, respectively, "asymptotically hyperbolic" assumptions) the previous quasi-local rigidity statement implies a \emph{global rigidity}: if every point in $M$ has a neighbourhood $U$ such that the supremum of the Hawking mass of surfaces contained in $U$ is non-positive, then $(M,g)$ is globally isometric to Euclidean ${\mathbb R}^3$ (resp. globally isometric to the Hyperbolic 3-space ${\mathbb H}^3$). Also, if the space is not flat (resp. not of constant sectional curvature $-1$), the methods give a small yet explicit and strictly positive lower bound on the Hawking mass of suitable spherical surfaces. We infer a small yet explicit and strictly positive lower bound on the Bartnik mass of open sets (non-locally isometric to Euclidean ${\mathbb R}^{3}$) in terms of curvature tensors. Inspired by these results, in the appendix we propose a notion of "sup-Hawking mass" which satisfies some natural properties of a quasi-local mass.


title: Federated Whole Prostate Segmentation in MRI with Personalized Neural Architectures

abstract: Building robust deep learning-based models requires diverse training data, ideally from several sources. However, these datasets cannot be combined easily because of patient privacy concerns or regulatory hurdles, especially if medical data is involved. Federated learning (FL) is a way to train machine learning models without the need for centralized datasets. Each FL client trains on their local data while only sharing model parameters with a global server that aggregates the parameters from all clients. At the same time, each client's data can exhibit differences and inconsistencies due to the local variation in the patient population, imaging equipment, and acquisition protocols. Hence, the federated learned models should be able to adapt to the local particularities of a client's data. In this work, we combine FL with an AutoML technique based on local neural architecture search by training a "supernet". Furthermore, we propose an adaptation scheme to allow for personalized model architectures at each FL client's site. The proposed method is evaluated on four different datasets from 3D prostate MRI and shown to improve the local models' performance after adaptation through selecting an optimal path through the AutoML supernet.


title: Hamiltonian Monte Carlo for Regression with High-Dimensional Categorical Data

abstract: Latent variable models are increasingly used in economics for high-dimensional categorical data like text and surveys. We demonstrate the effectiveness of Hamiltonian Monte Carlo (HMC) with parallelized automatic differentiation for analyzing such data in a computationally efficient and methodologically sound manner. Our new model, Supervised Topic Model with Covariates, shows that carefully modeling this type of data can have significant implications on conclusions compared to a simpler, frequently used, yet methodologically problematic, two-step approach. A simulation study and revisiting Bandiera et al. (2020)'s study of executive time use demonstrate these results. The approach accommodates thousands of parameters and doesn't require custom algorithms specific to each model, making it accessible for applied researchers


title: Single-photon detection with a Josephson junction coupled to a resonator

abstract: We use semiclassical formalism to optimize a microwave single photon detector based on switching events of a current biased Josephson junction coupled to a resonator. In order to detect very rare events, the average time between dark counts $\tau_{\rm dark}$ should be maximized taking into account that the switching time $\tau_{\rm sw}$ should be sufficiently small. We demonstrate that these times can be tuned in the wide range by changing the junction parameters, and the ratios $\tau_{\rm dark}/\tau_{\rm sw} \sim 10^9$ can be achieved. Therefore, a junction-resonator arrangement can be used for detecting extremely low photon fluxes, for instance for searching galactic axions.


title: Decentralized Multi-Agent Reinforcement Learning for Task Offloading Under Uncertainty

abstract: Multi-Agent Reinforcement Learning (MARL) is a challenging subarea of Reinforcement Learning due to the non-stationarity of the environments and the large dimensionality of the combined action space. Deep MARL algorithms have been applied to solve different task offloading problems. However, in real-world applications, information required by the agents (i.e. rewards and states) are subject to noise and alterations. The stability and the robustness of deep MARL to practical challenges is still an open research problem. In this work, we apply state-of-the art MARL algorithms to solve task offloading with reward uncertainty. We show that perturbations in the reward signal can induce decrease in the performance compared to learning with perfect rewards. We expect this paper to stimulate more research in studying and addressing the practical challenges of deploying deep MARL solutions in wireless communications systems.


title: Testing the AGN Radio and Neutrino correlation using the MOJAVE catalog and 10 years of IceCube Data

abstract: On 22 September 2017, IceCube reported a high-energy neutrino event which was found to be coincident with a flaring blazar, TXS 0506+056. This multi-messenger observation hinted at blazars contributing to the observed high-energy astrophysical neutrinos and raised a need for extensive correlation studies. Recent work shows that the internal absorption of gamma rays, and their interactions intrinsic to the source and with the extragalactic background, will cause a lack of energetic gamma-ray and neutrino correlation while hinting towards a correlation between neutrinos and lower photon energy observations in the X-ray and radio bands. Studies based on published IceCube alerts and radio observations report a possible radio-neutrino correlation in both gamma-ray bright and gamma-ray dim active galactic nuclei (AGN). However, they have marginal statistical significance due to limited available data. We present a correlation analysis between 15 GHz radio observations of AGN reported in the MOJAVE XV catalog and 10 years of IceCube detector data and discuss the results derived from a time-averaged stacking analysis.


title: Cell density controls signal propagation waves in a multicellular synthetic gene circuit

abstract: During organismal development, biochemical reaction networks sense and respond to mechanical forces to coordinate embryonic patterning with embryo morphogenesis. Factors such as cortical tension, cell density, and matrix mechanical properties influence differentiation and cell fate decisions by modulating gene regulatory signaling networks. A major goal in synthetic development is to construct gene regulatory circuits that program the patterning and morphogenesis of synthetic multicellular structures. However, in the synthetic context, little is known regarding how the physical properties of the growth environment impact the behavior of synthetic gene circuits. Here, we exploit physical-chemical coupling observed in a synthetic patterning circuit in order to control the size and spatial distribution of patterned synthetic cell sheets. We show that cell density attenuates the propagation of signal between neighboring cells in a multicellular sheet containing a contact-dependent patterning circuit based on the synNotch signaling system. Density-dependent attenuation leads to a signal propagation wave that exhibits distinct qualitative phases of persistent propagation, transient propagation, and no propagation. Through computational modeling, we demonstrate that cell growth parameters determine the phase of propagation observed within a growing cell sheet. Using growth-modulating drugs and spatial density gradients, we control the size of synNotch-activated cell populations and generate tissue-scale activation gradients and kinematic waves. Our study reveals that density-dependent synNotch activity can be exploited to control a synthetic multicellular patterning circuit. More broadly, we show that synthetic gene circuits can be critically impacted by their physical context, providing an alternate means for programming circuit behavior.


title: A skein relation for singular Soergel bimodules

abstract: We study the skein relation that governs the HOMFLYPT invariant of links colored by one-column Young diagrams. Our main result is a categorification of this colored skein relation. This takes the form of a homotopy equivalence between two one-sided twisted complexes constructed from Rickard complexes of singular Soergel bimodules associated to braided webs. Along the way, we prove a conjecture of Beliakova--Habiro relating the colored 2-strand full twist complex with the categorical ribbon element for quantum $\mathfrak{sl}_2$.


title: Inverse transport and diffusion problems in photoacoustic imaging with nonlinear absorption

abstract: Motivated by applications in imaging nonlinear optical absorption by photoacoustic tomography (PAT), we study in this work inverse coefficient problems for a semilinear radiative transport equation and its diffusion approximation with internal data that are functionals of the coefficients and the solutions to the equations. Based on the techniques of first- and second-order linearization, we derive uniqueness and stability results for the inverse problems. For uncertainty quantification purpose, we also establish the stability of the reconstruction of the absorption coefficients with respect to the change in the scattering coefficient.


title: Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport

abstract: We develop a neuroevolution-potential (NEP) framework for generating neural network based machine-learning potentials. They are trained using an evolutionary strategy for performing large-scale molecular dynamics (MD) simulations. A descriptor of the atomic environment is constructed based on Chebyshev and Legendre polynomials. The method is implemented in graphic processing units within the open-source GPUMD package, which can attain a computational speed over $10^7$ atom-step per second using one Nvidia Tesla V100. Furthermore, per-atom heat current is available in NEP, which paves the way for efficient and accurate MD simulations of heat transport in materials with strong phonon anharmonicity or spatial disorder, which usually cannot be accurately treated either with traditional empirical potentials or with perturbative methods.


title: Real-Time Mapping of Tissue Properties for Magnetic Resonance Fingerprinting

abstract: Magnetic resonance Fingerprinting (MRF) is a relatively new multi-parametric quantitative imaging method that involves a two-step process: (i) reconstructing a series of time frames from highly-undersampled non-Cartesian spiral k-space data and (ii) pattern matching using the time frames to infer tissue properties (e.g., T1 and T2 relaxation times). In this paper, we introduce a novel end-to-end deep learning framework to seamlessly map the tissue properties directly from spiral k-space MRF data, thereby avoiding time-consuming processing such as the nonuniform fast Fourier transform (NUFFT) and the dictionary-based Fingerprint matching. Our method directly consumes the non-Cartesian k- space data, performs adaptive density compensation, and predicts multiple tissue property maps in one forward pass. Experiments on both 2D and 3D MRF data demonstrate that quantification accuracy comparable to state-of-the-art methods can be accomplished within 0.5 second, which is 1100 to 7700 times faster than the original MRF framework. The proposed method is thus promising for facilitating the adoption of MRF in clinical settings.


title: Lattice-based simulation of the effects of nutrient concentration and magnetic field exposure on yeast colony growth and morphology

abstract: Yeasts exist in communities that expand over space and time to form complex structures and patterns. We developed a computational lattice-based framework to perform spatial-temporal simulations of budding yeast colonies exposed to different nutrient and magnetic field conditions. The budding patterns of haploid and diploid yeast cells were incorporated into the framework, as well as the filamentous growth that occurs in yeast colonies under nutrient limiting conditions. Simulation of the lattice-based model predicted that magnetic fields decrease colony growth rate, density, and roundness. Magnetic field simulations further predicted that colony elongation and boundary fluctuations increase in a nutrient- and ploidy-dependent manner. These in-silico predictions are an important step towards understanding the effects of the physico-chemical environment on microbial colonies and for informing bioelectromagnetic experiments on yeast colony biofilms and fungal pathogens.


title: Future Intelligent Autonomous Robots, Ethical by Design. Learning from Autonomous Cars Ethics

abstract: Development of the intelligent autonomous robot technology presupposes its anticipated beneficial effect on the individuals and societies. In the case of such disruptive emergent technology, not only questions of how to build, but also why to build and with what consequences are important. The field of ethics of intelligent autonomous robotic cars is a good example of research with actionable practical value, where a variety of stakeholders, including the legal system and other societal and governmental actors, as well as companies and businesses, collaborate bringing about shared view of ethics and societal aspects of technology. It could be used as a starting platform for the approaches to the development of intelligent autonomous robots in general, considering human-machine interfaces in different phases of the life cycle of technology - the development, implementation, testing, use and disposal. Drawing from our work on ethics of autonomous intelligent robocars, and the existing literature on ethics of robotics, our contribution consists of a set of values and ethical principles with identified challenges and proposed approaches for meeting them. This may help stakeholders in the field of intelligent autonomous robotics to connect ethical principles with their applications. Our recommendations of ethical requirements for autonomous cars can be used for other types of intelligent autonomous robots, with the caveat for social robots that require more research regarding interactions with the users. We emphasize that existing ethical frameworks need to be applied in a context-sensitive way, by assessments in interdisciplinary, multi-competent teams through multi-criteria analysis. Furthermore, we argue for the need of a continuous development of ethical principles, guidelines, and regulations, informed by the progress of technologies and involving relevant stakeholders.


title: High performance reduction technique for multiscale finite element modeling (HPR-FE2): towards industrial multiscale FE software

abstract: The authors have shown in previous contributions that reduced order modeling with optimal cubature applied to finite element square (FE2) techniques results in a reliable and affordable multiscale approach, the HPR-FE2 technique. Such technique is assessed here for an industrial case study of a generic 3D reinforced composite whose microstructure is represented by two general microcells accounting for different deformation mechanisms, microstrucural phases and geometry arrangement. Specifically, in this approach the microstrain modes used for building the reduced order model (ROM) are obtained through standard proper orthogonal decomposition (POD) techniques applied over snapshots of a representative sampling strain space. Additionally, a reduced number of integration points is obtained by exactly integrating the main free energy modes resulting from the sampling energy snapshots. The outcome consists of a number of dominant strain modes integrated over a remarkably reduced number of integration points which provide the support to evaluate the constitutive behavior of the microstructural phases. It is emphasized that stresses are computed according to the selected constitutive law at the reduced integration points and, therefore, the strategy inherits advantageous properties such as model completeness and customization of material properties. Overall results are discussed in terms of the consistency of the multiscale analysis, customization of the microscopic material parameters and speedup ratios compared to high-fidelity finite element (HF) simulations.


title: Architectures of Meaning, A Systematic Corpus Analysis of NLP Systems

abstract: This paper proposes a novel statistical corpus analysis framework targeted towards the interpretation of Natural Language Processing (NLP) architectural patterns at scale. The proposed approach combines saturation-based lexicon construction, statistical corpus analysis methods and graph collocations to induce a synthesis representation of NLP architectural patterns from corpora. The framework is validated in the full corpus of Semeval tasks and demonstrated coherent architectural patterns which can be used to answer architectural questions on a data-driven fashion, providing a systematic mechanism to interpret a largely dynamic and exponentially growing field.


title: Precision Electron Measurements in the Solar Wind at 1 au from NASA's Wind Spacecraft

abstract: This work aims to characterize precisely and systematically the non-thermal characteristics of the electron Velocity Distribution Function (eVDF) in the solar wind at 1 au using data from the Wind spacecraft. We present a comprehensive statistical analysis of solar wind electrons at 1 au using the electron analyzers of the 3D-Plasma instrument on board Wind. This work uses a sophisticated algorithm developed to analyze and characterize separately the three populations - core, halo and strahl - of the eVDF up to 2 keV. The eVDF data are calibrated using independent electron parameters obtained from the quasi-thermal noise around the electron plasma frequency measured by the Thermal Noise Receiver. The code determines the respective set of total electron, core, halo and strahl parameters through non-linear least-square fits to the measured eVDF, taking properly into account spacecraft charging and other instrumental effects. We use four years, ~ 280000 independent measurements of core, halo and strahl parameters to investigate the statistical properties of these different populations in the solar wind. We discuss the distributions of their respective densities, drift velocities, temperature, and temperature anisotropies as functions of solar wind speed. We also show distributions with solar wind speed of the total density, temperature, temperature anisotropy and heat flux, as well as those of the proton temperature, proton-to-electron temperature ratio, proton and electron beta. Intercorrelations between some of these parameters are also discussed. The present dataset represents the largest, high-precision, collection of electron measurements in the pristine solar wind at 1~AU. It provides a new wealth of information on electron microphysics. Its large volume will enable future statistical studies of parameter combinations and their dependencies under different plasma conditions.


title: On the Origin of the Pluto System

abstract: The goal of this chapter is to review hypotheses for the origin of the Pluto system in light of observational constraints that have been considerably refined over the 85-year interval between the discovery of Pluto and its exploration by spacecraft. We focus on the giant impact hypothesis currently understood as the likeliest origin for the Pluto-Charon binary, and devote particular attention to new models of planet formation and migration in the outer solar system. We discuss the origins conundrum posed by the system's four small moons. We also elaborate on the implications of these scenarios for the dynamical environment of the early transneptunian disk, the likelihood of finding a Pluto collisional family, and the origin of other binary systems in the Kuiper belt. Finally, we highlight outstanding open issues regarding the origins of the Pluto system and suggest areas of future progress.


title: Stability of superthermal strahl electrons in the solar wind

abstract: We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov-Maxwell linear solver LEOPARD. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfv\'en, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 AU. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modeled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.


title: The Law of Large Documents: Understanding the Structure of Legal Contracts Using Visual Cues

abstract: Large, pre-trained transformer models like BERT have achieved state-of-the-art results on document understanding tasks, but most implementations can only consider 512 tokens at a time. For many real-world applications, documents can be much longer, and the segmentation strategies typically used on longer documents miss out on document structure and contextual information, hurting their results on downstream tasks. In our work on legal agreements, we find that visual cues such as layout, style, and placement of text in a document are strong features that are crucial to achieving an acceptable level of accuracy on long documents. We measure the impact of incorporating such visual cues, obtained via computer vision methods, on the accuracy of document understanding tasks including document segmentation, entity extraction, and attribute classification. Our method of segmenting documents based on structural metadata out-performs existing methods on four long-document understanding tasks as measured on the Contract Understanding Atticus Dataset.


title: Proof of the Jacobi Property for the guiding-center Vlasov-Maxwell Bracket

abstract: The proof of the Jacobi property of the guiding-center Vlasov-Maxwell bracket underlying the Hamiltonian structure of the guiding-center Vlasov-Maxwell equations is presented.


title: A remark on conductor, depth and principal congruence subgroups

abstract: In this paper we study a relation between conductor, depth, and the level of principal congruence subgroups for irreducible admissible representations of ${\rm GL}_n(F)$ for a non-archimedean local field $F$ of characteristic zero. Some global and local applications are also discussed.


title: Analytic gradients in variational quantum algorithms: Algebraic extensions of the parameter-shift rule to general unitary transformations

abstract: Optimization of unitary transformations in Variational Quantum Algorithms benefits highly from efficient evaluation of cost function gradients with respect to amplitudes of unitary generators. We propose several extensions of the parametric-shift-rule to formulating these gradients as linear combinations of expectation values for generators with general eigen-spectrum (i.e. with more than two eigenvalues). Our approaches are exact and do not use any auxiliary qubits, instead they rely on a generator eigen-spectrum analysis. Two main directions in the parametric-shift-rule extensions are 1) polynomial expansion of the exponential unitary operator based on a limited number of different eigenvalues in the generator and 2) decomposition of the generator as a linear combination of low-eigenvalue operators (e.g. operators with only 2 or 3 eigenvalues). These techniques have a range of scalings for the number of needed expectation values with the number of generator eigenvalues from quadratic (for polynomial expansion) to linear and even $\log_2$ (for generator decompositions). This allowed us to propose efficient differentiation schemes superior to previous approaches for commonly used 2-qubit transformations (e.g. match-gates, transmon and fSim gates) and $\hat S^2$-conserving fermionic operators for the variational quantum eigensolver.


title: Loop Transformations using Clang's Abstract Syntax Tree

abstract: OpenMP 5.1 introduced the first loop nest transformation directives unroll and tile, and more are expected to be included in OpenMP 6.0. We discuss the two Abstract Syntax Tree (AST) representations used by Clang's implementation that is currently under development. The first representation is designed for compatibility with the existing implementation and stores the transformed loop nest in a shadow AST next to the syntactical AST. The second representation introduces a new meta AST-node OMPCanonicalLoop that guarantees that the semantic requirements of an OpenMP loop are met, and a CanonicalLoopInfo type that the OpenMPIRBuilder uses to represent literal and transformed loops. This second approach provides a better abstraction of loop semantics, removes the need for shadow AST nodes that are only relevant for code generation, allows sharing the implementation with other front-ends such as flang, but depends on the OpenMPIRBuilder which is currently under development.


title: Measurement of the Bjorken Sum at very low $Q^2$

abstract: We present new data on the Bjorken sum $\overline \Gamma_1^{p-n}(Q^2)$ at 4-momentum transfer $ 0.021 \leq Q^2 \leq 0.496$ GeV$^2$. The data were obtained in two experiments performed at Jefferson Lab: EG4 on polarized protons and deuterons, and E97110 on polarized $^3$He from which neutron data were extracted. The data cover the domain where chiral effective field theory ($\chi$EFT), the leading effective theory of the Strong Force at large distances, is expected to be applicable. We find that our data and the predictions from $\chi$EFT are only in marginal agreement. This is somewhat surprising as the contribution from the $\Delta(1232)$ resonance is suppressed in this observable, which should make it more reliably predicted by $\chi$EFT than quantities in which the $\Delta$ contribution is important. The data are also compared to a number of phenomenological models with various degrees of agreement.


title: A Nobile-like theorem for jet schemes of hypersurfaces

abstract: We prove that, for the jet scheme of a singular hypersurface, the blowup of a certain jet-related module is not an isomorphism. In conjunction with recent developments in the theory of Nash blowups, our result holds over fields of arbitrary characteristic. Our approach is based on explicit presentations given by a higher-order Jacobian matrix combined with a certain jet-related matrix.


title: Mediated Uncoupled Learning: Learning Functions without Direct Input-output Correspondences

abstract: Ordinary supervised learning is useful when we have paired training data of input $X$ and output $Y$. However, such paired data can be difficult to collect in practice. In this paper, we consider the task of predicting $Y$ from $X$ when we have no paired data of them, but we have two separate, independent datasets of $X$ and $Y$ each observed with some mediating variable $U$, that is, we have two datasets $S_X = \{(X_i, U_i)\}$ and $S_Y = \{(U'_j, Y'_j)\}$. A naive approach is to predict $U$ from $X$ using $S_X$ and then $Y$ from $U$ using $S_Y$, but we show that this is not statistically consistent. Moreover, predicting $U$ can be more difficult than predicting $Y$ in practice, e.g., when $U$ has higher dimensionality. To circumvent the difficulty, we propose a new method that avoids predicting $U$ but directly learns $Y = f(X)$ by training $f(X)$ with $S_{X}$ to predict $h(U)$ which is trained with $S_{Y}$ to approximate $Y$. We prove statistical consistency and error bounds of our method and experimentally confirm its practical usefulness.


title: Reflected and Doubly RBSDEs with Irregular Obstacles and a Large Set of Stopping Strategies

abstract: We introduce a new formulation of reflected BSDEs and doubly reflected BSDEs associated with irregular obstacles. In the first part of the paper, we consider an extension of the classical optimal stopping problem over a larger set of stopping systems than the set of stopping times (namely, the set of split stopping times), where the payoff process $\xi$ is irregular and in the case of a general filtration. Split stopping times are a powerful tool for modeling financial contracts and derivatives that depend on multiple conditions or triggers, and for incorporating stochastic processes with jumps and other types of discontinuities. We show that the value family can be aggregated by an optional process $v$, which is characterized as the Snell envelope of the reward process $\xi$ over split stopping times. Using this, we prove the existence and uniqueness of a solution $Y$ to irregular reflected BSDEs. In the second part of the paper, motivated by the classical Dynkin game with completely irregular rewards considered by Grigorova et al. (2018), we generalize the previous equations to the case of two reflecting barrier processes.


title: Delicate Ferromagnetism in MnBi$_6$Te$_{10}$

abstract: Tailoring magnetic orders in topological insulators is critical to the realization of topological quantum phenomena. An outstanding challenge is to find a material where atomic defects lead to tunable magnetic orders while maintaining a nontrivial topology. Here, by combining magnetization measurements, angle-resolved photoemission spectroscopy, and transmission electron microscopy, we reveal disorder-enabled, tunable magnetic ground states in MnBi$_6$Te$_{10}$. In the ferromagnetic phase, an energy gap of 15 meV is resolved at the Dirac point on the MnBi$_2$Te$_4$ termination. In contrast, antiferromagnetic MnBi$_6$Te$_{10}$ exhibits gapless topological surface states on all terminations. Transmission electron microscopy and magnetization measurements reveal substantial Mn vacancies and Mn migration in ferromagnetic MnBi$_6$Te$_{10}$. We provide a conceptual framework where a cooperative interplay of these defects drives a delicate change of overall magnetic ground state energies, and leads to tunable magnetic topological orders. Our work provides a clear pathway for nanoscale defect-engineering towards the realization of topological quantum phases.


title: Direct sampling methods for isotropic and anisotropic scatterers with point source measurements

abstract: In this paper, we consider the inverse scattering problem for recovering either an isotropic or anisotropic scatterer from the measured scattered field initiated by a point source. We propose two new imaging functionals for solving the inverse problem. The first one employs a far-field transform to the data which we then use to derive and provide an explicit decay rate for the imaging functional. In order to analyze the behavior of this imaging functional we use the factorization of the near field operator as well as the Funk-Hecke integral identity. For the second imaging functional the Cauchy data is used to define the functional and its behavior is analyzed using the Green's identities. Numerical experiments are given in two dimensions for both isotropic and anisotropic scatterers.


title: Upper Bounds on Resolvent Degree and Its Growth Rate

abstract: For each $n$, let $\text{RD}(n)$ denote the minimum $d$ for which there exists a formula for the general polynomial of degree $n$ in algebraic functions of at most $d$ variables. In 1945, Segre called for a better understanding of the large $n$ behavior of $\text{RD}(n)$. In this paper, we provide improved thresholds for upper bounds on $\text{RD}(n)$. Our techniques build upon classical algebraic geometry to provide new upper bounds for small $n$ and, in doing so, fix gaps in the proofs of A. Wiman and G.N. Chebotarev in [Wim1927] and [Che1954].


title: Markov Blanket Discovery using Minimum Message Length

abstract: Causal discovery automates the learning of causal Bayesian networks from data and has been of active interest from their beginning. With the sourcing of large data sets off the internet, interest in scaling up to very large data sets has grown. One approach to this is to parallelize search using Markov Blanket (MB) discovery as a first step, followed by a process of combining MBs in a global causal model. We develop and explore three new methods of MB discovery using Minimum Message Length (MML) and compare them empirically to the best existing methods, whether developed specifically as MB discovery or as feature selection. Our best MML method is consistently competitive and has some advantageous features.


title: An Automated Approach to Reasoning About Task-Oriented Insights in Responsive Visualization

abstract: Authors often transform a large screen visualization for smaller displays through rescaling, aggregation and other techniques when creating visualizations for both desktop and mobile devices (i.e., responsive visualization). However, transformations can alter relationships or patterns implied by the large screen view, requiring authors to reason carefully about what information to preserve while adjusting their design for the smaller display. We propose an automated approach to approximating the loss of support for task-oriented visualization insights (identification, comparison, and trend) in responsive transformation of a source visualization. We operationalize identification, comparison, and trend loss as objective functions calculated by comparing properties of the rendered source visualization to each realized target (small screen) visualization. To evaluate the utility of our approach, we train machine learning models on human ranked small screen alternative visualizations across a set of source visualizations. We find that our approach achieves an accuracy of 84% (random forest model) in ranking visualizations. We demonstrate this approach in a prototype responsive visualization recommender that enumerates responsive transformations using Answer Set Programming and evaluates the preservation of task-oriented insights using our loss measures. We discuss implications of our approach for the development of automated and semi-automated responsive visualization recommendation.


title: Autonomy 2.0: Why is self-driving always 5 years away?

abstract: Despite the numerous successes of machine learning over the past decade (image recognition, decision-making, NLP, image synthesis), self-driving technology has not yet followed the same trend. In this paper, we study the history, composition, and development bottlenecks of the modern self-driving stack. We argue that the slow progress is caused by approaches that require too much hand-engineering, an over-reliance on road testing, and high fleet deployment costs. We observe that the classical stack has several bottlenecks that preclude the necessary scale needed to capture the long tail of rare events. To resolve these problems, we outline the principles of Autonomy 2.0, an ML-first approach to self-driving, as a viable alternative to the currently adopted state-of-the-art. This approach is based on (i) a fully differentiable AV stack trainable from human demonstrations, (ii) closed-loop data-driven reactive simulation, and (iii) large-scale, low-cost data collections as critical solutions towards scalability issues. We outline the general architecture, survey promising works in this direction and propose key challenges to be addressed by the community in the future.


title: CoCo: Online Mixed-Integer Control via Supervised Learning

abstract: Many robotics problems, from robot motion planning to object manipulation, can be modeled as mixed-integer convex programs (MICPs). However, state-of-the-art algorithms are still unable to solve MICPs for control problems quickly enough for online use and existing heuristics can typically only find suboptimal solutions that might degrade robot performance. In this work, we turn to data-driven methods and present the Combinatorial Offline, Convex Online (CoCo) algorithm for quickly finding high quality solutions for MICPs. CoCo consists of a two-stage approach. In the offline phase, we train a neural network classifier that maps the problem parameters to a (logical strategy), which we define as the discrete arguments and relaxed big-M constraints associated with the optimal solution for that problem. Online, the classifier is applied to select a candidate logical strategy given new problem parameters; applying this logical strategy allows us to solve the original MICP as a convex optimization problem. We show through numerical experiments how CoCo finds near optimal solutions to MICPs arising in robot planning and control with 1 to 2 orders of magnitude solution speedup compared to other data-driven approaches and solvers.


title: Model for the long and orbital brightness variability of the $\beta$ Lyrae type binary OGLE-BLG-ECL-157529

abstract: Some close binaries of the beta Lyrae type show photometric cycles longer than the orbital one, which are possibly related to changes in their accretion disks. We aim to understand the short- and long-scale changes observed in the light curve of the eclipsing system OGLE-BLG-ECL-157529. In particular, we want to shed light on the contribution of the disk to these changes, especially those related to the long cycle, occurring on timescales of hundreds of days. We studied I-band OGLE photometric times series spanning 18.5 years, constructing disk models by analyzing the orbital light curve at 52 consecutive epochs. An optimized simplex algorithm was used to solve the inverse problem by adjusting the light curve with the best stellar-orbital-disk parameters for the system. We applied principal components analysis to the parameters to evaluate their dependence and variability. We constructed a description of the mass transfer rate in terms of disk parameters. We find that the light variability can be understood in terms of a variable mass transfer rate and variable accretion disk. The system brightness at orbital phase 0.25 follows the long cycle and is correlated with the mass transfer rate and the disk thickness. The long-cycle brightness variations can be understood in terms of differential occultation of the hotter star by a disk of variable thickness. Our model fits the overall light curve during 18.5 years well, including epochs of reversal of main and secondary eclipse depths. The disk radius cyclically change around the tidal radius, decoupled from changes in the mass transfer rate or system brightness, suggesting that viscous delay might explain the non-immediate response. Although the disk is large and fills a large fraction of the hot star Roche lobe, Lindblad resonance are far beyond the disk, excluding viscous dissipation as a major source of photometric variability.


title: Refactoring the MPS/University of Chicago Radiative MHD(MURaM) Model for GPU/CPU Performance Portability Using OpenACC Directives

abstract: The MURaM (Max Planck University of Chicago Radiative MHD) code is a solar atmosphere radiative MHD model that has been broadly applied to solar phenomena ranging from quiet to active sun, including eruptive events such as flares and coronal mass ejections. The treatment of physics is sufficiently realistic to allow for the synthesis of emission from visible light to extreme UV and X-rays, which is critical for a detailed comparison with available and future multi-wavelength observations. This component relies critically on the radiation transport solver (RTS) of MURaM; the most computationally intensive component of the code. The benefits of accelerating RTS are multiple fold: A faster RTS allows for the regular use of the more expensive multi-band radiation transport needed for comparison with observations, and this will pave the way for the acceleration of ongoing improvements in RTS that are critical for simulations of the solar chromosphere. We present challenges and strategies to accelerate a multi-physics, multi-band MURaM using a directive-based programming model, OpenACC in order to maintain a single source code across CPUs and GPUs. Results for a $288^3$ test problem show that MURaM with the optimized RTS routine achieves 1.73x speedup using a single NVIDIA V100 GPU over a fully subscribed 40-core Intel Skylake CPU node and with respect to the number of simulation points (in millions) per second, a single NVIDIA V100 GPU is equivalent to 69 Skylake cores. We also measure parallel performance on up to 96 GPUs and present weak and strong scaling results.


title: Picard understanding Darmok: A Dataset and Model for Metaphor-Rich Translation in a Constructed Language

abstract: Tamarian, a fictional language introduced in the Star Trek episode Darmok, communicates meaning through utterances of metaphorical references, such as "Darmok and Jalad at Tanagra" instead of "We should work together." This work assembles a Tamarian-English dictionary of utterances from the original episode and several follow-on novels, and uses this to construct a parallel corpus of 456 English-Tamarian utterances. A machine translation system based on a large language model (T5) is trained using this parallel corpus, and is shown to produce an accuracy of 76% when translating from English to Tamarian on known utterances.


title: AutoFL: Enabling Heterogeneity-Aware Energy Efficient Federated Learning

abstract: Federated learning enables a cluster of decentralized mobile devices at the edge to collaboratively train a shared machine learning model, while keeping all the raw training samples on device. This decentralized training approach is demonstrated as a practical solution to mitigate the risk of privacy leakage. However, enabling efficient FL deployment at the edge is challenging because of non-IID training data distribution, wide system heterogeneity and stochastic-varying runtime effects in the field. This paper jointly optimizes time-to-convergence and energy efficiency of state-of-the-art FL use cases by taking into account the stochastic nature of edge execution. We propose AutoFL by tailor-designing a reinforcement learning algorithm that learns and determines which K participant devices and per-device execution targets for each FL model aggregation round in the presence of stochastic runtime variance, system and data heterogeneity. By considering the unique characteristics of FL edge deployment judiciously, AutoFL achieves 3.6 times faster model convergence time and 4.7 and 5.2 times higher energy efficiency for local clients and globally over the cluster of K participants, respectively.


title: Declarative Machine Learning Systems

abstract: In the last years machine learning (ML) has moved from a academic endeavor to a pervasive technology adopted in almost every aspect of computing. ML-powered products are now embedded in our digital lives: from recommendations of what to watch, to divining our search intent, to powering virtual assistants in consumer and enterprise settings. Recent successes in applying ML in natural sciences revealed that ML can be used to tackle some of the hardest real-world problems humanity faces today. For these reasons ML has become central in the strategy of tech companies and has gathered even more attention from academia than ever before. Despite these successes, what we have witnessed so far is just the beginning. Right now the people training and using ML models are expert developers working within large organizations, but we believe the next wave of ML systems will allow a larger amount of people, potentially without coding skills, to perform the same tasks. These new ML systems will not require users to fully understand all the details of how models are trained and utilized for obtaining predictions. Declarative interfaces are well suited for this goal, by hiding complexity and favouring separation of interests, and can lead to increased productivity. We worked on such abstract interfaces by developing two declarative ML systems, Overton and Ludwig, that require users to declare only their data schema (names and types of inputs) and tasks rather then writing low level ML code. In this article we will describe how ML systems are currently structured, highlight important factors for their success and adoption, what are the issues current ML systems are facing and how the systems we developed addressed them. Finally we will talk about learnings from the development of ML systems throughout the years and how we believe the next generation of ML systems will look like.


title: Dual Quaternion-Based Visual Servoing for Grasping Moving Objects

abstract: This paper presents a new dual quaternion-based formulation for pose-based visual servoing. Extending our previous work on local contact moment (LoCoMo) based grasp planning, we demonstrate grasping of arbitrarily moving objects in 3D space. Instead of using the conventional axis-angle parameterization, dual quaternions allow designing the visual servoing task in a more compact manner and provide robustness to manipulator singularities. Given an object point cloud, LoCoMo generates a ranked list of grasp and pre-grasp poses, which are used as desired poses for visual servoing. Whenever the object moves (tracked by visual marker tracking), the desired pose updates automatically. For this, capitalising on the dual quaternion spatial distance error, we propose a dynamic grasp re-ranking metric to select the best feasible grasp for the moving object. This allows the robot to readily track and grasp arbitrarily moving objects. In addition, we also explore the robot null-space with our controller to avoid joint limits so as to achieve smooth trajectories while following moving objects. We evaluate the performance of the proposed visual servoing by conducting simulation experiments of grasping various objects using a 7-axis robot fitted with a 2-finger gripper. Obtained results demonstrate the efficiency of our proposed visual servoing.