title: On the Ternary Ohta-Kawasaki Free Energy and Its One-dimensional Global Minimizers

abstract: We study the ternary Ohta-Kawasaki free energy that has been used to model triblock copolymer systems. Its one-dimensional global minimizers are conjectured to have cyclic patterns. However, some physical experiments and computer simulations found triblock copolymers forming noncyclic lamellar patterns. In this work, by comparing the free energies of the cyclic pattern and some noncyclic candidates, we show that the conjecture does not hold for some choices of parameters. Our results suggest that even in one dimension, the global minimizers may take on very different patterns in different parameter regimes. To unify the existing choices of the long range coefficient matrix, we present a reformulation of the long range term using a generalized charge interpretation, and thereby propose conditions on the matrix in order for the global minimizers to reproduce physically relevant nanostructures of block copolymers.


title: Nonlocal conductivity, continued fractions and current vortices in electron fluids

abstract: Vortices in electron fluids are a key indicator of electron hydrodynamics. However, a comprehensive framework linking macroscopic vorticity measurements with microscopic interactions and scattering mechanisms has been lacking. We employ wavenumber-dependent conductivity $\sigma(k)$ incorporating realistic microscopic scattering processes, aiming to clarify the relationship between nonlocal response and vortices across ballistic and hydrodynamic phases. Vorticity is found to take similar values in both phases but feature very different sensitivity to momentum-relaxing scattering, with ballistic vortical flows being orders-of-magnitude more resilient than the hydrodynamic ones. This behavior can serve as a simple diagnostic of the microscopic origin of vorticity in electron fluids.


title: On the size of subsets of $\mathbb{F}_q^n$ avoiding solutions to linear systems with repeated columns

abstract: Consider a system of $m$ balanced linear equations in $k$ variables with coefficients in $\mathbb{F}_q$. If $k \geq 2m + 1$, then a routine application of the slice rank method shows that there are constants $\beta,\gamma \geq 1$ with $\gamma < q$ such that, for every subset $S \subseteq \mathbb{F}_q^n$ of size at least $\beta \cdot \gamma^n$, the system has a solution $(x_1,\ldots,x_k) \in S^k$ with $x_1,\ldots,x_k$ not all equal. Building on a series of papers by Mimura and Tokushige and on a paper by Sauermann, this paper investigates the problem of finding a solution of higher non-degeneracy; that is, a solution where $x_1,\ldots,x_k$ are pairwise distinct, or even a solution where $x_1,\ldots,x_k$ do not satisfy any balanced linear equation that is not a linear combination of the equations in the system. In this paper, we focus on linear systems with repeated columns. For a large class of systems of this type, we prove that there are constants $\beta,\gamma \geq 1$ with $\gamma < q$ such that every subset $S \subseteq \mathbb{F}_q^n$ of size at least $\beta \cdot \gamma^n$ contains a solution that is non-degenerate (in one of the two senses described above). This class is disjoint from the class covered by Sauermann's result, and captures the systems studied by Mimura and Tokushige into a single proof. Moreover, a special case of our results shows that, if $S \subseteq \mathbb{F}_p^n$ is a subset such that $S - S$ does not contain a non-trivial $k$-term arithmetic progression (with $p$ prime and $3 \leq k \leq p$), then $S$ must have exponentially small density.


title: Optimal control of PDEs using physics-informed neural networks

abstract: Physics-informed neural networks (PINNs) have recently become a popular method for solving forward and inverse problems governed by partial differential equations (PDEs). By incorporating the residual of the PDE into the loss function of a neural network-based surrogate model for the unknown state, PINNs can seamlessly blend measurement data with physical constraints. Here, we extend this framework to PDE-constrained optimal control problems, for which the governing PDE is fully known and the goal is to find a control variable that minimizes a desired cost objective. We provide a set of guidelines for obtaining a good optimal control solution; first by selecting an appropriate PINN architecture and training parameters based on a forward problem, second by choosing the best value for a critical scalar weight in the loss function using a simple but effective two-step line search strategy. We then validate the performance of the PINN framework by comparing it to adjoint-based nonlinear optimal control, which performs gradient descent on the discretized control variable while satisfying the discretized PDE. This comparison is carried out on several distributed control examples based on the Laplace, Burgers, Kuramoto-Sivashinsky, and Navier-Stokes equations. Finally, we discuss the advantages and caveats of using the PINN and adjoint-based approaches for solving optimal control problems constrained by nonlinear PDEs.


title: Restormer: Efficient Transformer for High-Resolution Image Restoration

abstract: Since convolutional neural networks (CNNs) perform well at learning generalizable image priors from large-scale data, these models have been extensively applied to image restoration and related tasks. Recently, another class of neural architectures, Transformers, have shown significant performance gains on natural language and high-level vision tasks. While the Transformer model mitigates the shortcomings of CNNs (i.e., limited receptive field and inadaptability to input content), its computational complexity grows quadratically with the spatial resolution, therefore making it infeasible to apply to most image restoration tasks involving high-resolution images. In this work, we propose an efficient Transformer model by making several key designs in the building blocks (multi-head attention and feed-forward network) such that it can capture long-range pixel interactions, while still remaining applicable to large images. Our model, named Restoration Transformer (Restormer), achieves state-of-the-art results on several image restoration tasks, including image deraining, single-image motion deblurring, defocus deblurring (single-image and dual-pixel data), and image denoising (Gaussian grayscale/color denoising, and real image denoising). The source code and pre-trained models are available at https://github.com/swz30/Restormer.


title: Long-range, Non-local Switching of Spin Textures in a Frustrated Antiferromagnet

abstract: Antiferromagnetic spintronics is an emerging area of quantum technologies that leverage the coupling between spin and orbital degrees of freedom in exotic materials. Spin-orbit interactions allow spin or angular momentum to be injected via electrical stimuli to manipulate the spin texture of a material, enabling the storage of information and energy. In general, the physical process is intrinsically local: spin is carried by an electrical current, imparted into the magnetic system, and the spin texture then rotates. The collective excitations of complex spin textures have rarely been utilized in this context, even though they can in principle transport spin over much longer distances, using much lower power. In this study, we show that spin information can be transported and stored non-locally in the material Fe$_x$NbS$_2$. We propose that collective modes leverage the strong magnetoelastic coupling in the system to achieve this, revealing a novel way to store spin information in complex magnetic systems


title: Swin Transformer V2: Scaling Up Capacity and Resolution

abstract: Large-scale NLP models have been shown to significantly improve the performance on language tasks with no signs of saturation. They also demonstrate amazing few-shot capabilities like that of human beings. This paper aims to explore large-scale models in computer vision. We tackle three major issues in training and application of large vision models, including training instability, resolution gaps between pre-training and fine-tuning, and hunger on labelled data. Three main techniques are proposed: 1) a residual-post-norm method combined with cosine attention to improve training stability; 2) A log-spaced continuous position bias method to effectively transfer models pre-trained using low-resolution images to downstream tasks with high-resolution inputs; 3) A self-supervised pre-training method, SimMIM, to reduce the needs of vast labeled images. Through these techniques, this paper successfully trained a 3 billion-parameter Swin Transformer V2 model, which is the largest dense vision model to date, and makes it capable of training with images of up to 1,536$\times$1,536 resolution. It set new performance records on 4 representative vision tasks, including ImageNet-V2 image classification, COCO object detection, ADE20K semantic segmentation, and Kinetics-400 video action classification. Also note our training is much more efficient than that in Google's billion-level visual models, which consumes 40 times less labelled data and 40 times less training time. Code is available at \url{https://github.com/microsoft/Swin-Transformer}.


title: Assisted Robust Reward Design

abstract: Real-world robotic tasks require complex reward functions. When we define the problem the robot needs to solve, we pretend that a designer specifies this complex reward exactly, and it is set in stone from then on. In practice, however, reward design is an iterative process: the designer chooses a reward, eventually encounters an "edge-case" environment where the reward incentivizes the wrong behavior, revises the reward, and repeats. What would it mean to rethink robotics problems to formally account for this iterative nature of reward design? We propose that the robot not take the specified reward for granted, but rather have uncertainty about it, and account for the future design iterations as future evidence. We contribute an Assisted Reward Design method that speeds up the design process by anticipating and influencing this future evidence: rather than letting the designer eventually encounter failure cases and revise the reward then, the method actively exposes the designer to such environments during the development phase. We test this method in a simplified autonomous driving task and find that it more quickly improves the car's behavior in held-out environments by proposing environments that are "edge cases" for the current reward.


title: Rate-optimal Bayesian Simple Regret in Best Arm Identification

abstract: We consider best arm identification in the multi-armed bandit problem. Assuming certain continuity conditions of the prior, we characterize the rate of the Bayesian simple regret. Differing from Bayesian regret minimization (Lai, 1987), the leading term in the Bayesian simple regret derives from the region where the gap between optimal and suboptimal arms is smaller than $\sqrt{\frac{\log T}{T}}$. We propose a simple and easy-to-compute algorithm with its leading term matching with the lower bound up to a constant factor; simulation results support our theoretical findings.


title: SimMIM: A Simple Framework for Masked Image Modeling

abstract: This paper presents SimMIM, a simple framework for masked image modeling. We simplify recently proposed related approaches without special designs such as block-wise masking and tokenization via discrete VAE or clustering. To study what let the masked image modeling task learn good representations, we systematically study the major components in our framework, and find that simple designs of each component have revealed very strong representation learning performance: 1) random masking of the input image with a moderately large masked patch size (e.g., 32) makes a strong pre-text task; 2) predicting raw pixels of RGB values by direct regression performs no worse than the patch classification approaches with complex designs; 3) the prediction head can be as light as a linear layer, with no worse performance than heavier ones. Using ViT-B, our approach achieves 83.8% top-1 fine-tuning accuracy on ImageNet-1K by pre-training also on this dataset, surpassing previous best approach by +0.6%. When applied on a larger model of about 650 million parameters, SwinV2-H, it achieves 87.1% top-1 accuracy on ImageNet-1K using only ImageNet-1K data. We also leverage this approach to facilitate the training of a 3B model (SwinV2-G), that by $40\times$ less data than that in previous practice, we achieve the state-of-the-art on four representative vision benchmarks. The code and models will be publicly available at https://github.com/microsoft/SimMIM.


title: PyTorchVideo: A Deep Learning Library for Video Understanding

abstract: We introduce PyTorchVideo, an open-source deep-learning library that provides a rich set of modular, efficient, and reproducible components for a variety of video understanding tasks, including classification, detection, self-supervised learning, and low-level processing. The library covers a full stack of video understanding tools including multimodal data loading, transformations, and models that reproduce state-of-the-art performance. PyTorchVideo further supports hardware acceleration that enables real-time inference on mobile devices. The library is based on PyTorch and can be used by any training framework; for example, PyTorchLightning, PySlowFast, or Classy Vision. PyTorchVideo is available at https://pytorchvideo.org/


title: Simple but Effective: CLIP Embeddings for Embodied AI

abstract: Contrastive language image pretraining (CLIP) encoders have been shown to be beneficial for a range of visual tasks from classification and detection to captioning and image manipulation. We investigate the effectiveness of CLIP visual backbones for Embodied AI tasks. We build incredibly simple baselines, named EmbCLIP, with no task specific architectures, inductive biases (such as the use of semantic maps), auxiliary tasks during training, or depth maps -- yet we find that our improved baselines perform very well across a range of tasks and simulators. EmbCLIP tops the RoboTHOR ObjectNav leaderboard by a huge margin of 20 pts (Success Rate). It tops the iTHOR 1-Phase Rearrangement leaderboard, beating the next best submission, which employs Active Neural Mapping, and more than doubling the % Fixed Strict metric (0.08 to 0.17). It also beats the winners of the 2021 Habitat ObjectNav Challenge, which employ auxiliary tasks, depth maps, and human demonstrations, and those of the 2019 Habitat PointNav Challenge. We evaluate the ability of CLIP's visual representations at capturing semantic information about input observations -- primitives that are useful for navigation-heavy embodied tasks -- and find that CLIP's representations encode these primitives more effectively than ImageNet-pretrained backbones. Finally, we extend one of our baselines, producing an agent capable of zero-shot object navigation that can navigate to objects that were not used as targets during training. Our code and models are available at https://github.com/allenai/embodied-clip


title: Multiwavelength View of the M87 Black Hole Captured by the Event Horizon Telescope

abstract: In 2017, the first image of the center of the M87 galaxy was captured by the Event Horizon Telescope (EHT). It revealed a ring morphology and a size consistent with theoretical expectations for the light pattern around a weakly accreting supermassive black hole of ~6.5 billion solar masses. In parallel to the EHT measurements, an extensive multi-wavelength (MWL) campaign with ground- and space-based facilities from radio all the way up to the TeV energy range was organized. In this note we will give an overview of the results from this campaign. M87 was found to be in a historically low state. At X-ray energies the emission from the core dominates over HST-1. We present the most complete simultaneous, MWL spectrum of the active nucleus to date.


title: Near-Chandrasekhar-Mass Type Ia Supernovae from the Double-Degenerate Channel

abstract: Recent observational evidence has demonstrated that white dwarf (WD) mergers are a highly efficient mechanism for mass accretion onto WDs in the galaxy. In this paper, we show that WD mergers naturally produce highly-magnetized, uniformly-rotating WDs, including a substantial population within a narrow mass range close to the Chandrasekhar mass ($M_{\rm Ch}$). These near-$M_{\rm Ch}$ WD mergers subsequently undergo rapid spin up and compression on a $\sim 10^2$ yr timescale, either leading to central ignition and a normal SN Ia via the DDT mechanism, or alternatively to a failed detonation and SN Iax through pure deflagration. The resulting SNe Ia and SNe Iax will have spectra, light curves, polarimetry, and nucleosynthetic yields similar to those predicted to arise through the canonical near-$M_{\rm Ch}$ single degenerate (SD) channel, but with a $t^{-1}$ delay time distribution characteristic of the double-degenerate (DD) channel. Furthermore, in contrast to the SD channel, WD merger near-$M_{\rm Ch}$ SNe Ia and SNe Iax will not produce observable companion signatures. We discuss a range of implications of these findings, from SNe Ia explosion mechanisms, to galactic nucleosynthesis of iron peak elements including manganese.


title: Effective dimensions of infinite-dimensional Hilbert spaces: A phase-space approach

abstract: By employing Husimi quasiprobability distributions, we show that a bounded portion of an unbounded phase space induces a finite effective dimension in an infinite dimensional Hilbert space. We compare our general expressions with numerical results for the spin-boson Dicke model in the chaotic energy regime, restricting its unbounded four-dimensional phase space to a classically chaotic energy shell. This effective dimension can be employed to characterize quantum phenomena in infinite dimensional systems, such as localization and scarring.


title: Upper Limit on the QCD Axion Mass from Isolated Neutron Star Cooling

abstract: The quantum chromodynamics (QCD) axion may modify the cooling rates of neutron stars (NSs). The axions are produced within the NS cores from nucleon bremsstrahlung and, when the nucleons are in superfluid states, Cooper pair breaking and formation processes. We show that four of the nearby isolated Magnificent Seven NSs along with PSR J0659 are prime candidates for axion cooling studies because they are coeval, with ages of a few hundred thousand years known from kinematic considerations, and they have well-measured surface luminosities. We compare these data to dedicated NS cooling simulations incorporating axions, profiling over uncertainties related to the equation of state, NS masses, surface compositions, and superfluidity. Our calculations of the axion and neutrino emissivities include high-density suppression factors that also affect SN 1987A and previous NS cooling limits on axions. We find no evidence for axions in the isolated NS data, and within the context of the KSVZ QCD axion model we constrain $m_a \lesssim 16$ meV at 95% confidence. An improved understanding of NS cooling and nucleon superfluidity could further improve these limits or lead to the discovery of the axion at weaker couplings.


title: Nonrelativistic Conformal Invariance in Mesoscopic Two-Dimensional Fermi Gases

abstract: Two-dimensional Fermi gases with universal short-range interactions are known to exhibit a quantum anomaly, where a classical scale and conformal invariance is broken by quantum effects at strong coupling. We argue that in a quasi two-dimensional geometry, a conformal window remains at weak interactions. Using degenerate perturbation theory, we verify the conformal symmetry by computing the energy spectrum of mesoscopic particle ensembles in a harmonic trap, which separates into conformal towers formed by so-called primary states and their center-of-mass and breathing-mode excitations, the latter having excitation energies at precisely twice the harmonic oscillator energy. In addition, using Metropolis importance sampling, we compute the hyperradial distribution function of the many-body wave functions, which are predicted by the conformal symmetry in closed analytical form. The weakly interacting Fermi gas constitutes a system where the nonrelativistic conformal symmetry can be revealed using elementary methods, and our results are testable in current experiments on mesoscopic Fermi gases.


title: Theory of a continuous bandwidth-tuned Wigner-Mott transition

abstract: We develop a theory for a continuous bandwidth-tuned transition at fixed \textit{fractional} electron filling from a metal with a generic Fermi surface to a `Wigner-Mott' insulator that spontaneously breaks crystalline space-group symmetries. Across the quantum critical point, (i) the entire electronic Fermi surface disappears abruptly upon approaching from the metallic side, and (ii) the insulating charge gap and various order parameters associated with the spontaneously broken space-group symmetries vanish continuously upon approaching from the insulating side. Additionally, the insulating side hosts a Fermi surface of neutral spinons. We present a framework for describing such continuous metal-insulator transitions (MITs) and analyze the example of a bandwidth-tuned transition at a filling, $\nu=1/6$, for spinful electrons on the triangular lattice. By extending the theory to a certain large-$N$ limit, we provide a concrete example of such a continuous MIT and discuss numerous experimental signatures near the critical point. We place our results in the context of recent experiments in moir\'e transition metal dichalcogenide materials.


title: Unitarity Bounds on Effective Field Theories at the LHC

abstract: Effective Field Theory (EFT) extensions of the Standard Model are tools to compute observables $\big(e.g.$ cross sections with partonic center-of-mass energy $\sqrt{\hat{s}}\,\big)$ as a systematically improvable expansion suppressed by a new physics scale $M$. If one is interested in EFT predictions in the parameter space where $M<\sqrt{\hat{s}}$, concerns of self-consistency emerge, which can manifest as a violation of perturbative partial-wave unitarity. However, when we search for the effects of an EFT at a hadron collider with center-of-mass energy $\sqrt{s}$ using an inclusive strategy, we typically do not have access to the event-by-event value of $\sqrt{\hat{s}}$. This motivates the need for a formalism that incorporates parton distribution functions into the perturbative partial-wave unitarity analysis. Developing such a framework and initiating an exploration of its implications is the goal of this work. Our approach opens up a potentially valid region of the EFT parameter space where $M \ll \sqrt{s}$. We provide evidence that there exist valid EFTs in this parameter space. The perturbative unitarity bounds are sensitive to the details of a given search, an effect we investigate by varying kinematic cuts.


title: Stochastic optimization for learning quantum state feedback control

abstract: High fidelity state preparation represents a fundamental challenge in the application of quantum technology. While the majority of optimal control approaches use feedback to improve the controller, the controller itself often does not incorporate explicit state dependence. Here, we present a general framework for training deep feedback networks for open quantum systems with quantum nondemolition measurement that allows a variety of system and control structures that are prohibitive by many other techniques and can in effect react to unmodeled effects through nonlinear filtering. We demonstrate that this method is efficient due to inherent parallelizability, robust to open system interactions, and outperforms landmark state feedback control results in simulation.


title: Chronostar. II. Kinematic age and substructure of the Scorpius-Centaurus OB2 association

abstract: The nearest region of massive star formation - the Scorpius-Centaurus OB2 association (Sco-Cen) - is a local laboratory ideally suited to the study of a wide range of astrophysical phenomena. Precision astrometry from the Gaia mission has expanded the census of this region by an order of magnitude. However, Sco-Cen's vastness and complex substructure make kinematic analysis of its traditional three regions, Upper Scorpius, Upper Centaurus-Lupus and Lower Centaurus-Crux, challenging. Here we use Chronostar, a Bayesian tool for kinematic age determination, to carry out a new kinematic decomposition of Sco-Cen using full 6-dimensional kinematic data. Our model identifies 8 kinematically distinct components consisting of 8,185 stars distributed in dense and diffuse groups, each with an independently-fit kinematic age; we verify that these kinematic estimates are consistent with isochronal ages. Both Upper Centaurus-Lupus and Lower Centaurus-Crux are split into two parts. The kinematic age of the component that includes PDS 70, one of the most well studied systems currently forming planets, is 15$\pm$3 Myr.


title: Cosmological constraints from the tomographic cross-correlation of DESI Luminous Red Galaxies and Planck CMB lensing

abstract: We use luminous red galaxies selected from the imaging surveys that are being used for targeting by the Dark Energy Spectroscopic Instrument (DESI) in combination with CMB lensing maps from the Planck collaboration to probe the amplitude of large-scale structure over $0.4\le z\le 1$. Our galaxy sample, with an angular number density of approximately $500\,\mathrm{deg}^{-2}$ over 18,000 sq.deg., is divided into 4 tomographic bins by photometric redshift and the redshift distributions are calibrated using spectroscopy from DESI. We fit the galaxy autospectra and galaxy-convergence cross-spectra using models based on cosmological perturbation theory, restricting to large scales that are expected to be well described by such models. Within the context of $\Lambda$CDM, combining all 4 samples and using priors on the background cosmology from supernova and baryon acoustic oscillation measurements, we find $S_8=\sigma_8(\Omega_m/0.3)^{0.5}=0.73\pm 0.03$. This result is lower than the prediction of the $\Lambda$CDM model conditioned on the Planck data. Our data prefer a slower growth of structure at low redshift than the model predictions, though at only modest significance.


title: The eMERLIN and EVN view of FR0 radio galaxies

abstract: We present the results from high-resolution observations carried out with the eMERLIN UK-array and the European VLBI network (EVN) for a sample of 15 FR0s, i.e. compact core-dominated radio sources associated with nearby early-type galaxies (ETGs) which represent the bulk of the local radio galaxy population. The 5-GHz eMERLIN observations available for 5 objects exhibit sub-mJy core components, and reveal pc-scale twin jets for 4 out of 5 FR0s once the eMERLIN and JVLA archival visibilities data are combined. The 1.66-GHz EVN observations available for 10 FR0s display one- and two-sided jetted morphologies and compact cores. The pc-scale core emission contributes, on average, to about one tenth of the total extended radio emission, although we note a increasing core contribution for flat/inverted-spectrum sources. We found an unprecedented linear correlation between the pc-scale core luminosity ($\sim$10$^{21.3}$-10$^{23.6}$ W Hz$^{-1}$) and [O III] line luminosity, generally considered as proxy of the accretion power, for a large sample of LINER-type radio-loud low-luminosity active nuclei, all hosted in massive ETGs, which include FR0s and FRIs. This result represents further evidence of a common jet-disc coupling in FR0s and FRIs, despite they differ in kpc-scale radio structure. For our objects and for other FR0 samples reported in the literature, we estimate the jet brightness sidedness ratios, which typically range between 1 and 3. This parameter roughly gauges the jet bulk Lorentz factor $\Gamma$, which turns out to range between 1 and 2.5 for most of the sample. This corroborates the scenario that FR0s are characterized by mildly-relativistic jets, possibly as a result of lower spinning black holes (BHs) than the highly-spinning BHs of relativistic-jetted radio galaxies, FRIs.


title: Fragmenting Active Asteroid 331P/Gibbs

abstract: We describe active asteroid 331P/Gibbs (2012 F5) using archival Hubble Space Telescope data taken between 2015 and 2018. 331P is an outer main-belt active asteroid with a long-lived debris trail that formed in 2011. Embedded in the debris trail we identify 19 fragments with radii between 0.04 and 0.11 km (albedo 0.05 assumed) containing about 1 percent of the mass of the primary nucleus. The largest shows a photometric range (1.5 magnitudes), a V-shaped minimum and a two-peaked lightcurve period near 9 hours, consistent with a symmetric contact binary (Drahus and Waniak 2016). Less convincing explanations are that 331P-A is a monolithic, elongated splinter or that its surface shows hemispheric 4:1 albedo variations. The debris trail is composed of centimeter sized and larger particles ejected with characteristic 10 cm s$^{-1}$ speeds following a size distribution with index q = 3.7+/-0.1 to 4.1+/-0.2. The HST data show that earlier, ground-based measurements of the nucleus were contaminated by near-nucleus debris, which cleared by 2015. We find that the primary nucleus has effective radius 0.8+/-0.1 km and is in rapid rotation (3.26+/-0.01 hour, range 0.25 magnitudes, minimum density 1600 kg/m3 if strengthless. The properties of 331P are consistent with a) formation about 1.5 Myr ago by impact shattering of a precursor body (Novakovic et al. 2014) b) spin-up by radiation torques to critical rotation c) ejection of about 1 percent of the nucleus mass in mid-2011 by rotational instability and d) subsequent evolution of the fragments and dispersal of the debris by radiation pressure.


title: Crosscap States in Integrable Field Theories and Spin Chains

abstract: We study crosscap states in integrable field theories and spin chains in 1+1 dimensions. We derive an exact formula for overlaps between the crosscap state and any excited state in integrable field theories with diagonal scattering. We then compute the crosscap entropy, i.e. the overlap for the ground state, in some examples. In the examples we analyzed, the result turns out to decrease monotonically along the renormalization group flow except in cases where the discrete symmetry is spontaneously broken in the infrared. We next introduce crosscap states in integrable spin chains, and obtain determinant expressions for the overlaps with energy eigenstates. These states are long-range entangled and their entanglement entropy grows linearly until the size of the subregion reaches half the system size. This property is reminiscent of pure-state black holes in holography and makes them interesting for use as initial conditions of quantum quench. As side results, we propose a generalization of Zamolodchikov's staircase model to flows between D-series minimal models, and discuss the relation to fermionic minimal models and the GSO projection.


title: A transformer-based model for default prediction in mid-cap corporate markets

abstract: In this paper, we study mid-cap companies, i.e. publicly traded companies with less than US $10 billion in market capitalisation. Using a large dataset of US mid-cap companies observed over 30 years, we look to predict the default probability term structure over the medium term and understand which data sources (i.e. fundamental, market or pricing data) contribute most to the default risk. Whereas existing methods typically require that data from different time periods are first aggregated and turned into cross-sectional features, we frame the problem as a multi-label time-series classification problem. We adapt transformer models, a state-of-the-art deep learning model emanating from the natural language processing domain, to the credit risk modelling setting. We also interpret the predictions of these models using attention heat maps. To optimise the model further, we present a custom loss function for multi-label classification and a novel multi-channel architecture with differential training that gives the model the ability to use all input data efficiently. Our results show the proposed deep learning architecture's superior performance, resulting in a 13% improvement in AUC (Area Under the receiver operating characteristic Curve) over traditional models. We also demonstrate how to produce an importance ranking for the different data sources and the temporal relationships using a Shapley approach specific to these models.


title: Surface Growth in Deformable Solids using an Eulerian Formulation

abstract: Accretion and ablation, i.e. the addition and removal of mass at the surface, is important in a wide range of physical processes including solidification, growth of biological tissues, environmental processes, and additive manufacturing. The description of accretion requires the addition of new continuum particles to the body, and is therefore challenging for standard continuum formulations for solids that require a reference configuration. Recent work has proposed an Eulerian approach to this problem, enabling the side-stepping of the issue of constructing the reference configuration. However, this raises the complementary challenge of determining the stress response of the solid, which typically requires the deformation gradient that is not immediately available in the Eulerian formulation. To resolve this, the approach introduced the elastic deformation as an additional kinematic descriptor of the added material, and its evolution has been shown to be governed by a transport equation. In this work, the method of characteristics is applied to solve concrete simplified problems motivated by biomechanics and manufacturing. Specifically, (1) for a problem with both ablation and accretion in a fixed domain, and (2) for a problem with a time-varying domain, the closed-form solution is obtained in the Eulerian framework using the method of characteristics without explicit construction of the reference configuration.


title: Operator backflow and the classical simulation of quantum transport

abstract: Tensor product states have proved extremely powerful for simulating the area-law entangled states of many-body systems, such as the ground states of gapped Hamiltonians in one dimension. The applicability of such methods to the \emph{dynamics} of many-body systems is less clear: the memory required grows exponentially in time in most cases, quickly becoming unmanageable. New methods reduce the memory required by selectively discarding/dissipating parts of the many-body wavefunction which are expected to have little effect on the hydrodynamic observables typically of interest: for example, some methods discard fine-grained correlations associated with $n$-point functions, with $n$ exceeding some cutoff $\ell_*$. In this work, we present a theory for the sizes of `backflow corrections', i.e., systematic errors due to discarding this fine-grained information. In particular, we focus on their effect on transport coefficients. Our results suggest that backflow corrections are exponentially suppressed in the size of the cutoff $\ell_*$. Moreover, the backflow errors themselves have a hydrodynamical expansion, which we elucidate. We test our predictions against numerical simulations run on random unitary circuits and ergodic spin-chains. These results lead to the conjecture that transport coefficients in ergodic diffusive systems can be captured to a given precision $\epsilon$ with an amount of memory scaling as $\exp[\mathcal{O}(\log(\epsilon)^2)]$, significantly better than the naive estimate of memory $\exp[\mathcal{O}(\mathrm{poly}(\epsilon^{-1}))]$ required by more brute-force methods.


title: Smart readout of nondestructive image sensors with single-photon sensitivity

abstract: Image sensors with nondestructive charge readout provide single-photon or single-electron sensitivity, but at the cost of long readout times. We present a smart readout technique to allow the use of these sensors in visible-light and other applications that require faster readout times. The method optimizes the readout noise and time by changing the number of times pixels are read out either statically, by defining an arbitrary number of regions of interest (ROI) in the array, or dynamically, depending on the charge or energy of interest (EOI) in the pixel. This technique is tested in a Skipper CCD showing that it is possible to obtain deep sub-electron noise, and therefore, high resolution of quantized charge, while dynamically changing the readout noise of the sensor. These faster, low noise readout techniques show that the skipper CCD is a competitive technology even where other technologies such as Electron Multiplier Charge Coupled Devices (EMCCD), silicon photo multipliers, etc. are currently used. This technique could allow skipper CCDs to benefit new astronomical instruments, quantum imaging, exoplanet search and study, and quantum metrology.


title: Solvable model for a charge-$4e$ superconductor

abstract: A charge-$4e$ superconductor forms due to the condensation of quartets of electrons. While in previous works the mechanism for the formation of charge-$4e$ superconductivity has been analyzed in terms of the binding of Cooper pairs in unconventional superconductors, its properties in the fermionic sector have not been studied systematically due to its inherently interacting nature even at the mean-field level. Here we propose a solvable model for a charge-$4e$ superconductor -- a spinful version of the Sachdev-Ye-Kitaev model with an anomalous quartic term. We show that the ground state is gapless and resembles a heavy Fermi liquid. We analytically solve for the superfluid density and show that it is perturbative in the strength of the charge-$4e$ order parameter, in sharp contrast with a regular (charge-$2e$) superconductor. Upon lowering temperature, we show that the correlation between charge-$4e$ order and regular interaction terms can drive a first-order phase transition to a charge-$2e$ superconducting state.


title: An Abstract Model for Branch and Cut

abstract: Branch and cut is the dominant paradigm for solving a wide range of mathematical programming problems -- linear or nonlinear -- combining efficient search (via branch and bound) and relaxation-tightening procedures (via cutting planes, or cuts). While there is a wealth of computational experience behind existing cutting strategies, there is simultaneously a relative lack of theoretical explanations for these choices, and for the tradeoffs involved therein. Recent papers have explored abstract models for branching and for comparing cuts with branch and bound. However, to model practice, it is crucial to understand the impact of jointly considering branching and cutting decisions. In this paper, we provide a framework for analyzing how cuts affect the size of branch-and-cut trees, as well as their impact on solution time. Our abstract model captures some of the key characteristics of real-world phenomena in branch-and-cut experiments, regarding whether to generate cuts only at the root or throughout the tree, how many rounds of cuts to add before starting to branch, and why cuts seem to exhibit nonmonotonic effects on the solution process.


title: Visual Goal-Directed Meta-Learning with Contextual Planning Networks

abstract: The goal of meta-learning is to generalize to new tasks and goals as quickly as possible. Ideally, we would like approaches that generalize to new goals and tasks on the first attempt. Toward that end, we introduce contextual planning networks (CPN). Tasks are represented as goal images and used to condition the approach. We evaluate CPN along with several other approaches adapted for zero-shot goal-directed meta-learning. We evaluate these approaches across 24 distinct manipulation tasks using Metaworld benchmark tasks. We found that CPN outperformed several approaches and baselines on one task and was competitive with existing approaches on others. We demonstrate the approach on a physical platform on Jenga tasks using a Kinova Jaco robotic arm.


title: Two Schr\"odinger-like Equations for hadrons

abstract: In this talk, based on [1,2], I argue that the holographic Schr\"odinger Equation of $(3+1)$-dim, conformal light-front QCD and the 't Hooft Equation of $(1+1)$-dim, large $N_c$ QCD, can be complementary to each other in providing a first approximation to hadron spectroscopy. Together, the two equations play a role in hadronic physics analogous that of the ordinary Schr\"odinger Equation in atomic physics.


title: Obstacles to Redistribution Through Markets and One Solution

abstract: Dworczak et al. (2021) study when certain market structures are optimal in the presence of heterogeneous preferences. A key assumption is that the social planner knows the joint distribution of the value of the good and marginal value of money. This paper studies whether relevant features of this distribution are identified from choice data. We show that the features of the distribution needed to characterize optimal market structure cannot be identified when demand is known for all prices. While this is a negative result, we show that the distribution of good value and marginal utility of money is fully identified when there is an observed measure of quality that varies. Thus, while Dworczak et al. (2021) abstract from quality, we show how including quality is crucial for potential applications.


title: Birational rigidity and K-stability of Fano hypersurfaces with ordinary double points

abstract: Extending previous results, we prove that for $n \ge 5$ all hypersurfaces of degree $n+1$ in ${\mathbb P}^{n+1}$ with isolated ordinary double points are birational superrigid and K-stable, hence admit a weak K\"ahler--Einstein metric.


title: Efficiently Transforming Tables for Joinability

abstract: Data from different sources rarely conform to a single formatting even if they describe the same set of entities, and this raises concerns when data from multiple sources must be joined or cross-referenced. Such a formatting mismatch is unavoidable when data is gathered from various public and third-party sources. Commercial database systems are not able to perform the join when there exist differences in data representation or formatting, and manual reformatting is both time consuming and error-prone. We study the problem of efficiently joining textual data under the condition that the join columns are not formatted the same and cannot be equi-joined, but they become joinable under some transformations. The problem is challenging simply because the number of possible transformations explodes with both the length of the input and the number of rows, even if each transformation is formed using very few basic units. We show that an efficient algorithm can be developed based on the common characteristics of the joined columns, and develop one such algorithm over a rich set of basic operations that can be composed to form transformations. We compare both the coverage and the running time of our algorithm to a state-of-the-art approach, and show that our algorithm covers every transformation that is covered in the state-of-the-art approach but is a few orders of magnitude faster, as evaluated on various real and synthetic data.


title: Min-max construction of minimal surfaces with a fixed angle at the boundary

abstract: We prove the existence of minimal surfaces in a bounded convex subset of $\mathbb R^3$, $\mathcal M$, intersecting the boundary of $\mathcal M$ with a fixed contact angle. The proof is based on a min-max construction in the spirit of Almgren-Pitts for the capillarity functional.


title: Jet energy spectrum and substructure in $e^+e^-$ collisions at 91.2 GeV with ALEPH Archived Data

abstract: The first measurements of energy spectra and substructure of anti-$k_{T}$ jets in hadronic $Z^0$ decays in $e^+e^-$ collisions are presented. The archived $e^+e^-$ annihilation data at a center-of-mass energy of 91.2 GeV were collected with the ALEPH detector at LEP in 1994. In addition to inclusive jet and leading dijet energy spectra, various jet substructure observables are analyzed as a function of jet energy which includes groomed and ungroomed jet mass to jet energy ratios, groomed momentum sharing, and groomed jet radius. The results are compared with perturbative QCD calculations and predictions from the SHERPA, HERWIG v7.1.5, PYTHIA 6, PYTHIA 8, and PYQUEN event generators. The jet energy spectra agree with perturbative QCD calculations which include the treatment of logarithms of the jet radius and threshold logarithms. None of the event generators give a fully satisfactory description of the data.


title: Quantum-Logic Gate between Two Optical Photons with an Average Efficiency above 40%

abstract: Optical qubits uniquely combine information transfer in optical fibers with a good processing capability and are therefore attractive tools for quantum technologies. A large challenge, however, is to overcome the low efficiency of two-qubit logic gates. The experimentally achieved efficiency in an optical controlled NOT (CNOT) gate reached approximately 11% in 2003 and has seen no increase since. Here we report on a new platform that was designed to surpass this long-standing record. The new scheme avoids inherently probabilistic protocols and, instead, combines aspects of two established quantum nonlinear systems: atom-cavity systems and Rydberg electromagnetically induced transparency. We demonstrate a CNOT gate between two optical photons with an average efficiency of 41.7(5)% at a postselected process fidelity of 81(2)%. Moreover, we extend the scheme to a CNOT gate with multiple target qubits and produce entangled states of presently up to five photons. All these achievements are promising and have the potential to advance optical quantum information processing in which almost all advanced protocols would profit from high-efficiency logic gates.


title: An Atomic Spectral Survey of WASP-76b: Resolving Chemical Gradients and Asymmetries

abstract: Ultra-hot Jupiters are gas giants that orbit so close to their host star that they are tidally locked, causing a permanent hot dayside and a cooler nightside. Signatures of their nonuniform atmospheres can be observed with high-resolution transit transmission spectroscopy by resolving time-dependent velocity shifts as the planet rotates and varying areas of the evening and morning terminator are probed. These asymmetric shifts were seen for the first time in iron absorption in WASP-76b. Here, we search for other atoms/ions in the planet's transmission spectrum and study the asymmetries in their signals. We detect Li I, Na I, Mg I, Ca II, V I, Cr I, Mn I, Fe I, Ni I, and Sr II, and tentatively detect H I, K I, and Co I, of which V, Cr, Ni, Sr II, and Co have not been reported before. We notably do not detect Ti or Al, even though these species should be readily observable, and hypothesize this could be due to condensation or cold trapping. We find that the observed signal asymmetries in the detected species can be explained in different ways. We find a relation between the expected condensation or ionization temperatures and the strength of the observed asymmetry, which could indicate rain-out or recombination on the nightside. However, we also find a dependence on the signal broadening, which could imply a two-zoned atmospheric model, in which the lower atmosphere is dominated by a day-to-night wind, while the upper atmosphere is dominated by a vertical wind or outflow. These observations provide a new level of modeling constraint and will aid our understanding of atmospheric dynamics in highly irradiated planets.


title: Interactive Set Discovery

abstract: We study the problem of set discovery where given a few example tuples of a desired set, we want to find the set in a collection of sets. A challenge is that the example tuples may not uniquely identify a set, and a large number of candidate sets may be returned. Our focus is on interactive exploration to set discovery where additional example tuples from the candidate sets are shown and the user either accepts or rejects them as members of the target set. The goal is to find the target set with the least number of user interactions. The problem can be cast as an optimization problem where we want to find a decision tree that can guide the search to the target set with the least number of questions to be answered by the user. We propose a general algorithm that is capable of reaching an optimal solution and two variations of it that strike a balance between the quality of a solution and the running time. We also propose a novel pruning strategy that safely reduces the search space without introducing false negatives. We evaluate the efficiency and the effectiveness of our algorithms through an extensive experimental study using both real and synthetic datasets and comparing them to previous approaches in the literature. We show that our pruning strategy reduces the running time of the search algorithms by 2-5 orders of magnitude.


title: Exploring leakage in dielectric films via automated experiment in scanning probe microscopy

abstract: Electronic conduction pathways in dielectric thin films are explored using automated experiments in scanning probe microscopy (SPM). Here, we use large field of view scanning to identify the position of localized conductive spots and develop a SPM workflow to probe their dynamic behavior at higher spatial resolution as a function of time, voltage, and scanning process in an automated fashion. Using this approach, we observe the variable behaviors of the conductive spots in a 20 nm thick ferroelectric Hf0.54Zr0.48O2 film, where conductive spots disappear and reappear during continuous scanning. There are also new conductive spots that appear during scanning. The automated workflow is universal and can be integrated into a wide range of microscopy techniques, including SPM, electron microscopy, optical microscopy, and chemical imaging.


title: Microwave studies of the three chiral ensembles in chains of coupled dielectric resonators

abstract: Random matrix theory has proven very successful in the understanding of the spectra of chaotic systems. Depending on symmetry with respect to time reversal and the presence or absence of a spin 1/2 there are three ensembles, the Gaussian orthogonal (GOE), Gaussian unitary (GUE), and Gaussian symplectic (GSE) one. With a further particle-antiparticle symmetry the chiral variants of these ensembles, the chiral orthogonal, unitary, and symplectic ensembles (the BDI, AIII, and CII in Cartan's notation) appear which are the main point of interest in this paper. Following a recently published work on chiral random matrix ensembles and their experimental realizations (Rehemanjiang et al., PRL 124, 116801 (2020)), this is achieved by using dielectric cylinders placed between two parallel aluminium plates. These cylinders act as microwave resonators which are used to create tight-binding chains of finite length up to N=5. The different ensembles are achieved by using different types of couplings: For the orthogonal case spatial proximity is used, for the unitary case microwave circulators are used, and for the symplectic case a combination of circulators and cables is used to create the necessary symmetry. In all cases the predicted repulsion behavior between positive and negative eigenvalues for energies close to zero are verified by a comparison with theory taking the finite size of the systems into account. We will show that the difference to the expected universal behavior is given by logarithmic corrections only. These corrections stem from the Hamiltonians having zero entries in their off-diagonal blocks.


title: BLUE: A 3D Dynamic Bipedal Robot

abstract: The objective of this work is to design a mechatronic bipedal robot with mobility in 3D environments. The designed robot has a total of six actuated degrees of freedom (DoF), each leg has two DoF located at the hip: one for abduction/adduction and another for thigh flexion/extension, and a third DoF at the knee for the shin flexion/extension. This robot is designed with point-feet legs to achieve a dynamic underactuated walking. Each actuator in the robot includes a DC gear motor, an encoder for position measurement, a flexible joint to form a series flexible actuator, and a feedback controller to ensure trajectory tracking. In order to reduce the total mass of the robot, the shin is designed using topology optimization. The resulting design is fabricated using 3D printed parts, which allows to get a robot's prototype to validate the selection of actuators. The preliminary experiments confirm the robot's ability to maintain a stand-up position and let us drawn future works in dynamic control and trajectory generation for periodic stable walking.


title: An outlier-resistant $\kappa$-generalized approach for robust physical parameter estimation

abstract: In this work we propose a robust methodology to mitigate the undesirable effects caused by outliers to generate reliable physical models. In this way, we formulate the inverse problems theory in the context of Kaniadakis statistical mechanics (or $\kappa$-statistics), in which the classical approach is a particular case. In this regard, the errors are assumed to be distributed according to a finite-variance $\kappa$-generalized Gaussian distribution. Based on the probabilistic maximum-likelihood method we derive a $\kappa$-objective function associated with the finite-variance $\kappa$-Gaussian distribution. To demonstrate our proposal's outlier-resistance, we analyze the robustness properties of the $\kappa$-objective function with help of the so-called influence function. In this regard, we discuss the role of the entropic index ($\kappa$) associated with the Kaniadakis $\kappa$-entropy in the effectiveness in inferring physical parameters by using strongly noisy data. In this way, we consider a classical geophysical data-inverse problem in two realistic circumstances, in which the first one refers to study the sensibility of our proposal to uncertainties in the input parameters, and the second is devoted to the inversion of a seismic data set contaminated by outliers. The results reveal an optimum $\kappa$-value at the limit $\kappa \rightarrow 2/3$, which is related to the best results.


title: A cut finite element method for the Darcy problem

abstract: We present and analyze a cut finite element method for the weak imposition of the Neumann boundary conditions of the Darcy problem. The Raviart-Thomas mixed element on both triangular and quadrilateral meshes is considered. Our method is based on the Nitsche formulation studied in [10.1515/jnma-2021-0042] and can be considered as a first attempt at extension in the unfitted case. The key feature is to add two ghost penalty operators to stabilize both the velocity and pressure fields. We rigorously prove our stabilized formulation to be well-posed and derive a priori error estimates for the velocity and pressure fields. We show that an upper bound for the condition number of the stiffness matrix holds as well. Numerical examples corroborating the theory are included.


title: Hybrid bounds for the sup-norm of automorphic forms in higher rank

abstract: Let $A$ be a central division algebra of prime degree $p$ over $\mathbb{Q}$. We obtain subconvex hybrid bounds, uniform in both the eigenvalue and the discriminant, for the sup-norm of Hecke-Maass forms on the compact quotients of $\operatorname{SL}_p(\mathbb{R})/\operatorname{SO}(p)$ by unit groups of orders in $A$. The exponents in the bounds are explicit and polynomial in $p$. We also prove subconvex hybrid bounds in the case of certain Eichler-type orders in division algebras of arbitrary odd degree.


title: Min-max theory for capillary surfaces

abstract: We develop a min-max theory for the construction of capillary surfaces in 3-manifolds with smooth boundary. In particular, for a generic set of ambient metrics, we prove the existence of nontrivial, smooth, almost properly embedded surfaces with any given constant mean curvature $c$, and with smooth boundary contacting at any given constant angle $\theta$. Moreover, if $c$ is nonzero and $\theta$ is not $\frac{\pi}{2}$, then our min-max solution always has multiplicity one. We also establish a stable Bernstein theorem for minimal hypersurfaces with certain contact angles in higher dimensions.


title: Ribbon Homology Cobordisms and Link Floer Homology

abstract: We make use of link Floer homology to study cobordisms between links embedded in 4-dimensional ribbon homology cobordisms. Combining results of Daemi--Lidman--Vela-Vick--Wong and Zemke, we show that ribbon homology concordances induce split injections on $\mathcal{HFL}^-$. We also make use of the torsion submodule of $\text{HFL}^-$ to give restrictions on the number of critical points in ribbon homology concordances.


title: Floquet dynamical quantum phase transitions under synchronized periodic driving

abstract: We study a generic class of fermionic two-band models under synchronized periodic driving, i.e., with the different terms in a Hamiltonian subject to periodic drives with the same frequency and phase. With all modes initially in a maximally mixed state, the synchronized drive is found to produce nonperiodic patterns of dynamical quantum phase transitions, with their appearance determined by an interplay of the band structure and the frequency of the drive. A case study of the anisotropic XY chain in a transverse magnetic field, transcribed to an effective two-band model, shows that the modes come with quantized geometric phases, allowing for the construction of an effective dynamical order parameter. Numerical studies in the limit of a strong magnetic field reveal distinct signals of precursors of dynamical quantum phase transitions also when the initial state of the XY chain is perturbed slightly away from maximal mixing, suggesting that the transitions may be accessible experimentally. A blueprint for an experiment built around laser-trapped circular Rydberg atoms is proposed.


title: Quality and Cost Trade-offs in Passage Re-ranking Task

abstract: Deep learning models named transformers achieved state-of-the-art results in a vast majority of NLP tasks at the cost of increased computational complexity and high memory consumption. Using the transformer model in real-time inference becomes a major challenge when implemented in production, because it requires expensive computational resources. The more executions of a transformer are needed the lower the overall throughput is, and switching to the smaller encoders leads to the decrease of accuracy. Our paper is devoted to the problem of how to choose the right architecture for the ranking step of the information retrieval pipeline, so that the number of required calls of transformer encoder is minimal with the maximum achievable quality of ranking. We investigated several late-interaction models such as Colbert and Poly-encoder architectures along with their modifications. Also, we took care of the memory footprint of the search index and tried to apply the learning-to-hash method to binarize the output vectors from the transformer encoders. The results of the evaluation are provided using TREC 2019-2021 and MS Marco dev datasets.


title: New physics searches with the International Large Detector at the ILC

abstract: Although the LHC experiments have searched for and excluded many proposed new particles up to masses close to 1 TeV, there are many scenarios that are difficult to address at a hadron collider. This talk will review a number of these scenarios and present the expectations for searches at an electron-positron collider such as the International Linear Collider. The cases discussed include the light Higgsino, the \stau~ slepton in the coannihilation region relevant to dark matter, as well as other BSM signatures. The studies are based on the ILD concept at the ILC.


title: Scaling violation in interaction of cosmic ray hadrons and the nature of the 3 PeV knee in the spectrum of primary cosmic rays

abstract: Energy spectra of the most energetic hadrons in the core region of extensive air showers (EAS) were studied in dependence on the EAS energy $E_0$ in the hybrid experiment HADRON at the Tien Shan mountain cosmic ray station. For the first time by direct measurement it was found in this experiment that the slope of the power energy spectrum of EAS hadrons diminishes itself locally, and the average hadron energy, correspondingly, increases between the $E_0$ borders 3PeV and 20PeV. Such behavior agrees with threshold appearance in EAS, in the same energy range, of a long-flying penetrative component which was earlier revealed at the Tien Shan station. Now we reconsider this Tien Shan result in comparison with the new data of modern collider experiments. The analysis permits to state unambiguously an astrophysical nature of the penetrative EAS component, and to discuss its probable connection with the origin of the famous knee in the energy spectrum of cosmic rays at $E_0\simeq 3$ PeV.


title: Learning To Estimate Regions Of Attraction Of Autonomous Dynamical Systems Using Physics-Informed Neural Networks

abstract: When learning to perform motor tasks in a simulated environment, neural networks must be allowed to explore their action space to discover new potentially viable solutions. However, in an online learning scenario with physical hardware, this exploration must be constrained by relevant safety considerations in order to avoid damage to the agent's hardware and environment. We aim to address this problem by training a neural network, which we will refer to as a "safety network", to estimate the region of attraction (ROA) of a controlled autonomous dynamical system. This safety network can thereby be used to quantify the relative safety of proposed control actions and prevent the selection of damaging actions. Here we present our development of the safety network by training an artificial neural network (ANN) to represent the ROA of several autonomous dynamical system benchmark problems. The training of this network is predicated upon both Lyapunov theory and neural solutions to partial differential equations (PDEs). By learning to approximate the viscosity solution to a specially chosen PDE that contains the dynamics of the system of interest, the safety network learns to approximate a particular function, similar to a Lyapunov function, whose zero level set is boundary of the ROA. We train our safety network to solve these PDEs in a semi-supervised manner following a modified version of the Physics Informed Neural Network (PINN) approach, utilizing a loss function that penalizes disagreement with the PDE's initial and boundary conditions, as well as non-zero residual and variational terms. In future work we intend to apply this technique to reinforcement learning agents during motor learning tasks.


title: DawDreamer: Bridging the Gap Between Digital Audio Workstations and Python Interfaces

abstract: Audio production techniques which previously only existed in GUI-constrained digital audio workstations, livecoding environments, or C++ APIs are now accessible with our new Python module called DawDreamer. DawDreamer therefore bridges the gap between real sound engineers and coders imitating them with offline batch-processing. Like contemporary modules in this domain, DawDreamer can create directed acyclic graphs of audio processors such as VSTs which generate or manipulate audio streams. DawDreamer can also dynamically compile and execute code from Faust, a powerful signal processing language which can be deployed to many platforms and microcontrollers. We discuss DawDreamer's unique features in detail and potential applications across music information retrieval including source separation, transcription, and audio effect parameter inference. We provide fully cross-platform PyPI installers, a Linux Dockerfile, and an example Jupyter notebook.


title: Minimum Resource Threshold Policy Under Partial Interference

abstract: When developing policies for prevention of infectious diseases, policymakers often set specific, outcome-oriented targets to achieve. For example, when developing a vaccine allocation policy, policymakers may want to distribute them so that at least a certain fraction of individuals in a census block are disease-free and spillover effects due to interference within blocks are accounted for. The paper proposes methods to estimate a block-level treatment policy that achieves a pre-defined, outcome-oriented target while accounting for spillover effects due to interference. Our policy, the minimum resource threshold policy (MRTP), suggests the minimum fraction of treated units required within a block to meet or exceed the target level of the outcome. We estimate the MRTP from empirical risk minimization using a novel, nonparametric, doubly robust loss function. We then characterize statistical properties of the estimated MRTP in terms of the excess risk bound. We apply our methodology to design a water, sanitation, and hygiene allocation policy for Senegal with the goal of increasing the proportion of households with no children experiencing diarrhea to a level exceeding a specified threshold. Our policy outperforms competing policies and offers new approaches to design allocation policies, especially in international development for communicable diseases.


title: Loss Functions for Discrete Contextual Pricing with Observational Data

abstract: We study a pricing setting where each customer is offered a contextualized price based on customer and/or product features. Often only historical sales data are available, so we observe whether a customer purchased a product at the price prescribed rather than the customer's true valuation. Such observational data are influenced by historical pricing policies, which introduce difficulties in evaluating the effectiveness of future policies. The goal of this paper is to formulate loss functions that can be used for evaluating pricing policies directly from observational data, rather than going through an intermediate demand estimation stage, which may suffer from bias. To achieve this, we adapt ideas from machine learning with corrupted labels, where we consider each observed purchase decision as a known probabilistic transformation of the customer's valuation. From this transformation, we derive a class of unbiased loss functions. Within this class, we identify minimum variance estimators and estimators robust to poor demand estimation. Furthermore, we show that for contextual pricing, estimators popular in the off-policy evaluation literature fall within this class of loss functions. We offer managerial insights into scenarios under which these estimators are effective.


title: Constraint-based Diversification of JOP Gadgets

abstract: Modern software deployment process produces software that is uniform and hence vulnerable to large-scale code-reuse attacks, such as Jump-Oriented Programming (JOP) attacks. Compiler-based diversification improves the resilience of software systems by automatically generating different assembly code versions of a given program. Existing techniques are efficient but do not have a precise control over the quality of the generated variants. This paper introduces Diversity by Construction (DivCon), a constraint-based approach to software diversification. Unlike previous approaches, DivCon allows users to control and adjust the conflicting goals of diversity and code quality. A key enabler is the use of Large Neighborhood Search (LNS) to generate highly diverse code efficiently. For larger problems, we propose a combination of LNS with a structural decomposition of the problem. To further improve the diversification efficiency of DivCon against JOP attacks, we propose an application-specific distance measure tailored to the characteristics of JOP attacks. We evaluate DivCon with 20 functions from a popular benchmark suite for embedded systems. These experiments show that the combination of LNS and our application-specific distance measure generates binary programs that are highly resilient against JOP attacks. Our results confirm that there is a trade-off between the quality of each assembly code version and the diversity of the entire pool of versions. In particular, the experiments show that DivCon generates near-optimal binary programs that share a small number of gadgets. For constraint programming researchers and practitioners, this paper demonstrates that LNS is a valuable technique for finding diverse solutions. For security researchers and software engineers, DivCon extends the scope of compiler-based diversification to performance-critical and resource-constrained applications.


title: A Conformer-based ASR Frontend for Joint Acoustic Echo Cancellation, Speech Enhancement and Speech Separation

abstract: We present a frontend for improving robustness of automatic speech recognition (ASR), that jointly implements three modules within a single model: acoustic echo cancellation, speech enhancement, and speech separation. This is achieved by using a contextual enhancement neural network that can optionally make use of different types of side inputs: (1) a reference signal of the playback audio, which is necessary for echo cancellation; (2) a noise context, which is useful for speech enhancement; and (3) an embedding vector representing the voice characteristic of the target speaker of interest, which is not only critical in speech separation, but also helpful for echo cancellation and speech enhancement. We present detailed evaluations to show that the joint model performs almost as well as the task-specific models, and significantly reduces word error rate in noisy conditions even when using a large-scale state-of-the-art ASR model. Compared to the noisy baseline, the joint model reduces the word error rate in low signal-to-noise ratio conditions by at least 71% on our echo cancellation dataset, 10% on our noisy dataset, and 26% on our multi-speaker dataset. Compared to task-specific models, the joint model performs within 10% on our echo cancellation dataset, 2% on the noisy dataset, and 3% on the multi-speaker dataset.


title: From the Harmonic Oscillator to Time-Frequency Analysis of Chirp Signals

abstract: This paper presents a novel approach to understanding the role of harmonic dynamics and gaining a deeper appreciation for its impact within and outside of quantum mechanics. This includes consequences of harmonic dynamics and the uncertainty principle for anomalous diffusion and for the time-frequency analysis of chirp signals. In this approach, we consider a contact transformation to view a system of canonical variables with coordinate $x$ and momentum $p_x$ in the context of a new system of "generalized" coordinates and momentum. This new system is first studied in the context of non-relativistic quantum mechanics. The classical analog is then explored by use of the Poisson bracket equation. From this, new implications are demonstrated in classical phenomena. One is for a new model of Anomalous and Normal Diffusion. In another, we introduce the concept of the "Mixed Fourier Transform" which explores a new Gaussian Fourier Transform kernel in terms of the generalized variables. This has the ultimate objective of "harmonizing" chirp signals or producing a harmonic signal from an otherwise non-harmonic chirp.


title: Anomalous convective flows carve pinnacles and scallops in melting ice

abstract: We report on the shape dynamics of ice suspended in cold fresh water and subject to the natural convective flows generated during melting. Experiments reveal shape motifs for increasing far-field temperature: Sharp pinnacles directed downward at low temperatures, scalloped waves for intermediate temperatures between 5 and $7^\circ$C, and upward pointing pinnacles at higher temperatures. Phase-field simulations reproduce these morphologies, which are closely linked to the anomalous density-temperature profile of liquid water. Boundary layer flows yield pinnacles that sharpen with accelerating growth of tip curvature while scallops emerge from a Kelvin-Helmholtz-like instability caused by counterflowing currents that roll up to form vortex arrays. By linking the molecular-scale effects underlying water's density anomaly to the macro-scale flows that imprint the surface, these results show that the morphology of melted ice is a sensitive indicator of ambient temperature.


title: Multiplicative summations into algebraically closed fields

abstract: In this paper, extending our earlier program, we derive maximal canonical extensions for multiplicative summations into algebraically closed fields. We show that there is a well-defined analogue to minimal polynomials for a series algebraic over a ring of series, the "scalar polynomial". When that ring is the domain of a summation $\mathfrak{S}$, we derive the related concepts of the $\mathfrak{S}$-minimal polynomial for a series, which is mapped by $\mathfrak{S}$ to a scalar polynomial. When the scalar polynomial for a series has the form $(t-a)^n$, $a$ is the unique value to which the series can be mapped by an extension of the original summation.


title: Explainable predictions of different machine learning algorithms used to predict Early Stage diabetes

abstract: Machine Learning and Artificial Intelligence can be widely used to diagnose chronic diseases so that necessary precautionary treatment can be done in critical time. Diabetes Mellitus which is one of the major diseases can be easily diagnosed by several Machine Learning algorithms. Early stage diagnosis is crucial to prevent dangerous consequences. In this paper we have made a comparative analysis of several machine learning algorithms viz. Random Forest, Decision Tree, Artificial Neural Networks, K Nearest Neighbor, Support Vector Machine, and XGBoost along with feature attribution using SHAP to identify the most important feature in predicting the diabetes on a dataset collected from Sylhet Hospital. As per the experimental results obtained, the Random Forest algorithm has outperformed all the other algorithms with an accuracy of 99 percent on this particular dataset.


title: Chiral Phase Change Nanomaterials

abstract: Chiral nanostructures offer the ability to respond to the vector nature of a light beam at the nanoscale. While naturally chiral materials offer a path towards scalability, engineered structures offer a path to wavelength tunability through geometric manipulation. Neither approach, however, allows for temporal control of chirality. Therefore, in the best of all worlds, it is crucial to realize chiral materials that possess the quality of scalability, tailored wavelength response, and dynamic control at high speeds. Here, a new class of intrinsically chiral phase change nanomaterials (PCNMs) is proposed and explored, based on a scalable bottom-up fabrication technique with a high degree of control in three dimensions. Angular resolved Mueller Matrix and spectroscopic ellipsometry are performed to characterize the optical birefringence and dichroism, and a numerical model is provided to explain the origin of optical activity. This work achieves the critical goal of demonstrating high-speed dynamic switching of chirality over 50,000 cycles via the underlying PCNM.


title: Dyson gas on a curved contour

abstract: We introduce and study a model of a logarithmic gas with inverse temperature $\beta$ on an arbitrary smooth closed contour in the plane. This model generalizes Dyson's gas (the $\beta$-ensemble) on the unit circle. We compute the non-vanishing terms of the large $N$ expansion of the free energy ($N$ is the number of particles) by iterating the "loop equation" that is the Ward identity with respect to reparametrizations and dilatation of the contour. We show that the main contribution to the free energy is expressed through the spectral determinant of the Neumann jump operator associated with the contour, or equivalently through the Fredholm determinant of the Neumann-Poincare (double layer) operator. This result connects the statistical mechanics of the Dyson gas to the spectral geometry of the interior and exterior domains of the supporting contour.


title: An incompressible Eulerian method for fluid-structure interaction with mixed soft and rigid solids

abstract: We present a general simulation approach for incompressible fluid--structure interactions in a fully Eulerian framework using the reference map technique (RMT). The approach is suitable for modeling one or more rigid or finitely-deformable objects or soft objects with rigid components interacting with the fluid and with each other. It is also extended to control the kinematics of structures in fluids. The model is based on our previous Eulerian fluid--soft solver, and generalized to rigid structures by constraining the deformation-rate tensor in a projection framework. Several numerical examples are presented to illustrate the capability of the method.


title: Real-Time Magnetometry Using Dark States of a Nitrogen Vacancy Center

abstract: We demonstrate real-time magnetometry by detecting fluorescence from a nitrogen vacancy center in the setting of coherent population trapping and by estimating magnetic field from the time series of the observed photon counts, which are correlated with the underlying field. The proof-of-principle experiment uses an external time-varying magnetic field that follows an Ornstein-Uhlenbeck (OU) process. By taking into consideration the statistical properties of the OU process, a Bayesian inference-based estimator can effectively update dynamical information of the field in real time with the detection of just a single photon.


title: 3N potentials in the Faddeev coordinate space approach to Nd scattering

abstract: In the last decade, for studying 3$N$ bound states and $Nd$ scattering the Tucson-Melbourne (TM) and Urbana 3$N$ force derived from the chiral EFT have been applied. We plan to use the TM 3$N$ force for studying the $Nd$ scattering on the basis of the Faddeev equations in configuration space. In the given paper, we present our final formulas for components of the TM 3$N$ potential obtained in the coordinate space.


title: Metallic magnetic calorimeter arrays for the first phase of the ECHo experiment

abstract: The ECHo experiment has been designed for the determination of the effective electron neutrino mass by means of the analysis of the end-point region of the Ho-163 electron capture spectrum. Metallic magnetic calorimeters enclosing Ho-163 are used for the high energy resolution calorimetric measurement of the Ho-163 spectrum. For the first phase of the experiment, ECHo-1k, a 72-pixel MMC array has been developed. The single-pixel design has been optimised to reach 100% stopping power for the radiation emitted in the Ho-163 electron capture process (besides the electron neutrino) and an energy resolution < 10 eV FWHM. We describe the design of the ECHo-1k detector chip, the fabrication steps and the characterisation at room temperature, at 4 K and at the final operation temperatures. In particular, a detailed analysis of the results from these tests allowed to define a quality check protocol based on parameters measurable at room temperature. We discuss the performance achieved with the two ECHo-1k detector chips - the first one with Ho-163 implanted in gold and the second one with Ho-163 implanted in silver - which have been used for the high statistics measurement of the ECHo-1k experiment. An average activity per pixel of 0.81 Bq and 0.71 Bq and an average energy resolution of 6.07 eV FWHM and 5.55 eV FWHM have been achieved with these two detectors, fulfilling the requirements for the first phase of the ECHo experiment.


title: More on Emergent Dark Energy from Unparticles

abstract: In a recent paper \cite{Artymowski:2020zwy} we suggested the possibility that the present acceleration of the Universe is due to thermodynamical behavior of unparticles. The model is free of scalar fields, modified gravity, a Cosmological Constant (CC), the coincidence problem, initial conditions problem and possesses interesting distinct predictions regarding the equation of state of Dark Energy, the growth rate and the number of relativistic degrees of freedom at BBN and CMB decoupling. In this work, we relate to a recent paper \cite{Abchouyeh:2021wey}, which discusses a similar setup of unparticles with and without a CC as an external source of late-time acceleration. The authors have shown how such a model is inconsistent with the data. We show that these claims are viable only in a particular part of the parameter space and that model \cite{Artymowski:2020zwy} stands tall. We further suggest a consistency condition in terms of observables. We then fit publicly available supernovae data to derive the expected Hubble parameter and constrain the parameters of the model.


title: Parallel Algorithms for Masked Sparse Matrix-Matrix Products

abstract: Computing the product of two sparse matrices (SpGEMM) is a fundamental operation in various combinatorial and graph algorithms as well as various bioinformatics and data analytics applications for computing inner-product similarities. For an important class of algorithms, only a subset of the output entries are needed, and the resulting operation is known as Masked SpGEMM since a subset of the output entries is considered to be "masked out". Existing algorithms for Masked SpGEMM usually do not consider mask as part of multiplication and either first compute a regular SpGEMM followed by masking, or perform a sparse inner product only for output elements that are not masked out. In this work, we investigate various novel algorithms and data structures for this rather challenging and important computation, and provide guidelines on how to design a fast Masked-SpGEMM for shared-memory architectures. Our evaluations show that factors such as matrix and mask density, mask structure and cache behavior play a vital role in attaining high performance for Masked SpGEMM. We evaluate our algorithms on a large number of matrices using several real-world benchmarks and show that our algorithms in most cases significantly outperform the state of the art for Masked SpGEMM implementations.


title: Algebraic Presentations of Dependent Type Theories

abstract: C-systems were defined by Cartmell as models of generalized algebraic theories. B-systems were defined by Voevodsky in his quest to formulate and prove an initiality conjecture for type theories. They play a crucial role in Voevodsky's construction of a syntactic C-system from a term monad. In this work, we construct an equivalence between the category of C-systems and the category of B-systems, thus proving a conjecture by Voevodsky. We construct this equivalence as the restriction of an equivalence between more general structures, called CE-systems and E-systems, respectively. To this end, we identify C-systems and B-systems as "stratified" CE-systems and E-systems, respectively; that is, systems whose contexts are built iteratively via context extension, starting from the empty context.


title: A fast algorithm for computing the Smith normal form with multipliers for a nonsingular integer matrix

abstract: A Las Vegas randomized algorithm is given to compute the Smith multipliers for a nonsingular integer matrix $A$, that is, unimodular matrices $U$ and $V$ such that $AV=US$, with $S$ the Smith normal form of $A$. The expected running time of the algorithm is about the same as required to multiply together two matrices of the same dimension and size of entries as $A$. Explicit bounds are given for the size of the entries in both unimodular multipliers. The main tool used by the algorithm is the Smith massager, a relaxed version of $V$, the unimodular matrix specifying the column operations of the Smith computation. From the perspective of efficiency, the main tools used are fast linear solving and partial linearization of integer matrices. As an application of the Smith with multipliers algorithm, a fast algorithm is given to find the fractional part of the inverse of the input matrix.


title: Correcting Face Distortion in Wide-Angle Videos

abstract: Video blogs and selfies are popular social media formats, which are often captured by wide-angle cameras to show human subjects and expanded background. Unfortunately, due to perspective projection, faces near corners and edges exhibit apparent distortions that stretch and squish the facial features, resulting in poor video quality. In this work, we present a video warping algorithm to correct these distortions. Our key idea is to apply stereographic projection locally on the facial regions. We formulate a mesh warp problem using spatial-temporal energy minimization and minimize background deformation using a line-preservation term to maintain the straight edges in the background. To address temporal coherency, we constrain the temporal smoothness on the warping meshes and facial trajectories through the latent variables. For performance evaluation, we develop a wide-angle video dataset with a wide range of focal lengths. The user study shows that 83.9% of users prefer our algorithm over other alternatives based on perspective projection.


title: Time-Optimal Paths for Simple Cars with Moving Obstacles in the Hamilton-Jacobi Formulation

abstract: We consider the problem of time-optimal path planning for simple nonholonomic vehicles. In previous similar work, the vehicle has been simplified to a point mass and the obstacles have been stationary. Our formulation accounts for a rectangular vehicle, and involves the dynamic programming principle and a time-dependent Hamilton-Jacobi-Bellman (HJB) formulation which allows for moving obstacles. To our knowledge, this is the first HJB formulation of the problem which allows for moving obstacles. We design an upwind finite difference scheme to approximate the equation and demonstrate the efficacy of our model with a few synthetic examples.


title: Dispersion chain of Vlasov equations

abstract: On the basis of the Vlasov chain of equations, a new infinite dispersion chain of equations is obtained for the distribution functions of mixed higher order kinematical values. In contrast to the Vlasov chain, the dispersion chain contains distribution functions with an arbitrary set of kinematical values and has a tensor form of writing. For the dispersion chain, new equations for mixed Boltzmann functions and the corresponding chain of conservation laws for fluid dynamics are obtained. The probability is proved to be a constant value for a particle to belong the region where the quasi-probability density is negative (Wigner function).


title: Golden rhythms as a theoretical framework for cross-frequency organization

abstract: While brain rhythms appear fundamental to brain function, why brain rhythms consistently organize into the small set of discrete frequency bands observed remains unknown. Here we propose that rhythms separated by factors of the golden ratio ($\phi=(1+ \sqrt{5})/2$) optimally support segregation and cross-frequency integration of information transmission in the brain. Organized by the golden ratio, pairs of transient rhythms support multiplexing by reducing interference between separate communication channels, and triplets of transient rhythms support integration of signals to establish a hierarchy of cross-frequency interactions. We illustrate this framework in simulation and apply this framework to propose four hypotheses.


title: MS-nowcasting: Operational Precipitation Nowcasting with Convolutional LSTMs at Microsoft Weather

abstract: We present the encoder-forecaster convolutional long short-term memory (LSTM) deep-learning model that powers Microsoft Weather's operational precipitation nowcasting product. This model takes as input a sequence of weather radar mosaics and deterministically predicts future radar reflectivity at lead times up to 6 hours. By stacking a large input receptive field along the feature dimension and conditioning the model's forecaster with predictions from the physics-based High Resolution Rapid Refresh (HRRR) model, we are able to outperform optical flow and HRRR baselines by 20-25% on multiple metrics averaged over all lead times.


title: SmartSlice: Dynamic, self-optimization of applications QoS requests to 5G networks

abstract: Applications can tailor a network slice by specifying a variety of QoS attributes related to application-specific performance, function or operation. However, some QoS attributes like guaranteed bandwidth required by the application do vary over time. For example, network bandwidth needs of video streams from surveillance cameras can vary a lot depending on the environmental conditions and the content in the video streams. In this paper, we propose a novel, dynamic QoS attribute prediction technique that assists any application to make optimal resource reservation requests at all times. Standard forecasting using traditional cost functions like MAE, MSE, RMSE, MDA, etc. don't work well because they do not take into account the direction (whether the forecasting of resources is more or less than needed), magnitude (by how much the forecast deviates, and in which direction), or frequency (how many times the forecast deviates from actual needs, and in which direction). The direction, magnitude and frequency have a direct impact on the application's accuracy of insights, and the operational costs. We propose a new, parameterized cost function that takes into account all three of them, and guides the design of a new prediction technique. To the best of our knowledge, this is the first work that considers time-varying application requirements and dynamically adjusts slice QoS requests to 5G networks in order to ensure a balance between application's accuracy and operational costs. In a real-world deployment of a surveillance video analytics application over 17 cameras, we show that our technique outperforms other traditional forecasting methods, and it saves 34% of network bandwidth (over a ~24 hour period) when compared to a static, one-time reservation.


title: Anomalous Lifetimes of Ultracold Complexes Decaying into a Single Channel: What's Taking So Long in There?

abstract: We investigate the lifetimes of complexes formed in ultracold molecule collisions. Employing both transition-state-theory and an optical model approach we examine processes that can extend the lifetime of complexes beyond that predicted by Rice-Ramsperger-Kassel-Marcus theory. We focus on complexes that possess only one open channel, and find that the extreme distribution of widths for this case favors low decay rates. Thus decay from a complex into a single energetically available channel can be anomalously slow, and moreover nonexponential in time. We apply the theory to two systems of current experimental interest, RbCs and NaRb, finding qualitatively that the empirical time scales can be accounted for in the theory.


title: Rethinking Dilated Convolution for Real-time Semantic Segmentation

abstract: The field-of-view is an important metric when designing a model for semantic segmentation. To obtain a large field-of-view, previous approaches generally choose to rapidly downsample the resolution, usually with average poolings or stride 2 convolutions. We take a different approach by using dilated convolutions with large dilation rates throughout the backbone, allowing the backbone to easily tune its field-of-view by adjusting its dilation rates, and show that it's competitive with existing approaches. To effectively use the dilated convolution, we show a simple upper bound on the dilation rate in order to not leave gaps in between the convolutional weights, and design an SE-ResNeXt inspired block structure that uses two parallel $3\times 3$ convolutions with different dilation rates to preserve the local details. Manually tuning the dilation rates for every block can be difficult, so we also introduce a differentiable neural architecture search method that uses gradient descent to optimize the dilation rates. In addition, we propose a lightweight decoder that restores local information better than common alternatives. To demonstrate the effectiveness of our approach, our model RegSeg achieves competitive results on real-time Cityscapes and CamVid datasets. Using a T4 GPU with mixed precision, RegSeg achieves 78.3 mIOU on Cityscapes test set at $37$ FPS, and 80.9 mIOU on CamVid test set at $112$ FPS, both without ImageNet pretraining.


title: A Nodal Immersed Finite Element-Finite Difference Method

abstract: The immersed finite element-finite difference (IFED) method is a computational approach to modeling interactions between a fluid and an immersed structure. This method uses a finite element (FE) method to approximate the stresses and forces on a structural mesh and a finite difference (FD) method to approximate the momentum of the entire fluid-structure system on a Cartesian grid. The fundamental approach used by this method follows the immersed boundary framework for modeling fluid-structure interaction (FSI), in which a force spreading operator prolongs structural forces to a Cartesian grid, and a velocity interpolation operator restricts a velocity field defined on that grid back onto the structural mesh. Force spreading and velocity interpolation both require projecting data onto the finite element space. Consequently, evaluating either coupling operator requires solving a matrix equation at every time step. Mass lumping, in which the projection matrices are replaced by diagonal approximations, has the potential to accelerate this method considerably. Constructing the coupling operators also requires determining the locations on the structure mesh where the forces and velocities are sampled. Here we show that sampling the forces and velocities at the nodes of the structural mesh is equivalent to using lumped mass matrices in the coupling operators. A key theoretical result of our analysis is that if both of these approaches are used together, the IFED method permits the use of lumped mass matrices derived from nodal quadrature rules for any standard interpolatory element. This is different from standard FE methods, which require specialized treatments to accommodate mass lumping with higher-order shape functions. Our theoretical results are confirmed by numerical benchmarks, including standard solid mechanics tests and examination of a dynamic model of a bioprosthetic heart valve.


title: Predictive Scheduling of Collaborative Mobile Robots for Improved Crop-transport Logistics of Manually Harvested Crops

abstract: Mechanizing the manual harvesting of fresh market fruits constitutes one of the biggest challenges to the sustainability of the fruit industry. During manual harvesting of some fresh-market crops like strawberries and table grapes, pickers spend significant amounts of time walking to carry full trays to a collection station at the edge of the field. A step toward increasing harvest automation for such crops is to deploy harvest-aid collaborative robots (co-bots) that transport the empty and full trays, thus increasing harvest efficiency by reducing pickers' non-productive walking times. This work presents the development of a co-robotic harvest-aid system and its evaluation during commercial strawberry harvesting. At the heart of the system lies a predictive stochastic scheduling algorithm that minimizes the expected non-picking time, thus maximizing the harvest efficiency. During the evaluation experiments, the co-robots improved the mean harvesting efficiency by around 10% and reduced the mean non-productive time by 60%, when the robot-to-picker ratio was 1:3. The concepts developed in this work can be applied to robotic harvest-aids for other manually harvested crops that involve walking for crop transportation.


title: Reining in Mobile Web Performance with Document and Permission Policies

abstract: The quality of experience with the mobile web remains poor, partially as a result of complex websites and design choices that worsen performance, particularly for users in suboptimal networks or devices. Prior proposed solutions have seen limited adoption due in part to the demand they place on developers and content providers, and the performing infrastructure needed to support them. We argue that Document and Permissions Policies -- an ongoing effort to enforce good practices on web design -- may offer the basis for a readily-available and easily-adoptable solution. In this paper, we evaluate the potential performance cost of violating well-understood policies and how common such violations are in today's web. Our analysis show, for example, that controlling for unsized-media policy, something applicable to 70% of the top-1million websites, can indeed reduce Cumulative Layout Shift metric.


title: A Review of Adversarial Attack and Defense for Classification Methods

abstract: Despite the efficiency and scalability of machine learning systems, recent studies have demonstrated that many classification methods, especially deep neural networks (DNNs), are vulnerable to adversarial examples; i.e., examples that are carefully crafted to fool a well-trained classification model while being indistinguishable from natural data to human. This makes it potentially unsafe to apply DNNs or related methods in security-critical areas. Since this issue was first identified by Biggio et al. (2013) and Szegedy et al.(2014), much work has been done in this field, including the development of attack methods to generate adversarial examples and the construction of defense techniques to guard against such examples. This paper aims to introduce this topic and its latest developments to the statistical community, primarily focusing on the generation and guarding of adversarial examples. Computing codes (in python and R) used in the numerical experiments are publicly available for readers to explore the surveyed methods. It is the hope of the authors that this paper will encourage more statisticians to work on this important and exciting field of generating and defending against adversarial examples.