title: Test-Time Adaptation with Shape Moments for Image Segmentation

abstract: Supervised learning is well-known to fail at generalization under distribution shifts. In typical clinical settings, the source data is inaccessible and the target distribution is represented with a handful of samples: adaptation can only happen at test time on a few or even a single subject(s). We investigate test-time single-subject adaptation for segmentation, and propose a Shape-guided Entropy Minimization objective for tackling this task. During inference for a single testing subject, our loss is minimized with respect to the batch normalization's scale and bias parameters. We show the potential of integrating various shape priors to guide adaptation to plausible solutions, and validate our method in two challenging scenarios: MRI-to-CT adaptation of cardiac segmentation and cross-site adaptation of prostate segmentation. Our approach exhibits substantially better performances than the existing test-time adaptation methods. Even more surprisingly, it fares better than state-of-the-art domain adaptation methods, although it forgoes training on additional target data during adaptation. Our results question the usefulness of training on target data in segmentation adaptation, and points to the substantial effect of shape priors on test-time inference. Our framework can be readily used for integrating various priors and for adapting any segmentation network, and our code is available.


title: The Cosmic Large-Scale Structure in X-rays (CLASSIX) cluster survey IV: Superclusters in the local Universe at z <= 0.03

abstract: It is important to map the large-scale matter distribution in the local Universe for cosmological studies, such as the tracing of the large-scale peculiar velocity flow, the characterisation of the environment for different astronomical objects, and for precision measurements of cosmological parameters. We used X-ray luminous clusters to map this matter distribution and find that about 51% of the groups and clusters are members of superclusters which occupy only a few percent of the volume. In this paper we provide a detailed description of these large-scale structures. With a friends-to-friends algorithm, we find eight superclusters with a cluster overdensity ratio of at least two with five or more galaxy group and cluster members in the cosmic volume out to z = 0.03. The four most prominent ones are the Perseus-Pisces, the Centaurus, the Coma, and the Hercules supercluster, with lengths from about 40 to over 100 Mpc and estimated masses of 0.6 - 2.2 10^16 Msun. The largest of these structures is the Perseus-Pisces supercluster. The four smaller superclusters include the Local and the Abell 400 supercluster and two superclusters in the constellations Sagittarius and Lacerta. We provide detailed maps, member catalogues, and physical descriptions of the eight superclusters. By constructing superclusters with a range of cluster sub-samples with different lower X-ray luminosity limits, we show that the main structures are always reliably recovered.


title: Expert Systems with Logic#. A Novel Modeling Framework for Logic Programming in an Object-Oriented Context of C#

abstract: We present a novel approach how logic programming for expert systems can be declared directly in an object-oriented language.


title: The twisted group algebra structure of the Cayley-Dickson algebra

abstract: The Cayley-Dickson algebra has long been a challenge due to the lack of an explicit multiplication table. Despite being constructible through inductive construction, its explicit structure has remained elusive until now. In this article, we propose a solution to this long-standing problem by revealing the Cayley-Dickson algebra as a twisted group algebra with an explicit twist function $\sigma(A,B)$. We show that this function satisfies the equation $$e_Ae_B=(-1)^{\sigma(A,B)}e_{A\oplus B}$$ and provide a formula for the relationship between the Cayley-Dickson algebra and split Cayley-Dickson algebra, thereby giving an explicit expression for the twist function of the split Cayley-Dickson algebra. Our approach not only resolves the lack of explicit structure for the Cayley-Dickson algebra and split Cayley-Dickson algebra but also sheds light on the algebraic structure underlying this fundamental mathematical object.


title: Cost Minimization in Acquisition for Gravitational Wave Surrogate Modeling

abstract: Gravitational wave science is dependent upon expensive numerical simulations, which provide the foundational understanding of binary merger radiation needed to interpret observations of massive binary black holes. The high cost of these simulations limits large-scale campaigns to explore the binary black hole parameter space. Surrogate models have been developed to efficiently interpolate between simulation results, but these models require a sufficiently comprehensive sample to train on. Acquisition functions can be used to identify points in the domain for simulation. We develop a new acquisition function which accounts for the cost of simulating new points. We show that when applied to a 3D domain of binary mass ratio and dimensionless spins, the accumulated cost of simulation is reduced by a factor of about 10.


title: Optimal method for reconstructing polychromatic maps from broadband observations with an asymmetric antenna pattern

abstract: Broadband time-ordered data obtained from telescopes with a wavelength-dependent, asymmetric beam pattern can be used to extract maps at multiple wavelengths from a single scan. This technique is especially useful when collecting data on cosmic phenomena such as the Cosmic Microwave Background (CMB) radiation, as it provides the ability to separate the CMB signal from foreground contaminants. We develop a method to determine the optimal linear combinations of wavelengths (``colors'') that can be reconstructed for a given telescope design and the number of colors that are measurable with high signal-to-noise ratio. The optimal colors are found as eigenvectors of a matrix derived from the inverse noise covariance matrix. When the telescope is able to scan the sky isotropically, it is useful to transform to a spherical harmonic basis, in which this matrix has a particularly simple form. We propose using the optimal colors determined from the isotropic case even when the actual scanning pattern is not isotropic (e.g., covers only part of the sky). We perform simulations showing that maps in multiple colors can be reconstructed accurately from both full-sky and partial-sky scans. Although the original motivation for this research comes from mapping the CMB, this method of polychromatic map-making will have broader applications throughout astrophysics.


title: Neural Eikonal Solver: improving accuracy of physics-informed neural networks for solving eikonal equation in case of caustics

abstract: The concept of physics-informed neural networks has become a useful tool for solving differential equations due to its flexibility. There are a few approaches using this concept to solve the eikonal equation which describes the first-arrival traveltimes of acoustic and elastic waves in smooth heterogeneous velocity models. However, the challenge of the eikonal is exacerbated by the velocity models producing caustics, resulting in instabilities and deterioration of accuracy due to the non-smooth solution behaviour. In this paper, we revisit the problem of solving the eikonal equation using neural networks to tackle the caustic pathologies. We introduce the novel Neural Eikonal Solver (NES) for solving the isotropic eikonal equation in two formulations: the one-point problem is for a fixed source location; the two-point problem is for an arbitrary source-receiver pair. We present several techniques which provide stability in velocity models producing caustics: improved factorization; non-symmetric loss function based on Hamiltonian; gaussian activation; symmetrization. In our tests, NES showed the relative-mean-absolute error of about 0.2-0.4% from the second-order factored Fast Marching Method, and outperformed existing neural-network solvers giving 10-60 times lower errors and 2-30 times faster training. The inference time of NES is comparable with the Fast Marching. The one-point NES provides the most accurate solution, whereas the two-point NES provides slightly lower accuracy but gives an extremely compact representation. It can be useful in various seismic applications where massive computations are required (millions of source-receiver pairs): ray modeling, traveltime tomography, hypocenter localization, and Kirchhoff migration.


title: Empowering Optimal Control with Machine Learning: A Perspective from Model Predictive Control

abstract: Solving complex optimal control problems have confronted computational challenges for a long time. Recent advances in machine learning have provided us with new opportunities to address these challenges. This paper takes model predictive control, a popular optimal control method, as the primary example to survey recent progress that leverages machine learning techniques to empower optimal control solvers. We also discuss some of the main challenges encountered when applying machine learning to develop more robust optimal control algorithms.


title: TopSort: A High-Performance Two-Phase Sorting Accelerator Optimized on HBM-based FPGAs

abstract: The emergence of high-bandwidth memory (HBM) brings new opportunities to boost the performance of sorting acceleration on FPGAs, which was conventionally bounded by the available off-chip memory bandwidth. However, it is nontrivial for designers to fully utilize this immense bandwidth. First, the existing sorter designs cannot be directly scaled at the increasing rate of available off-chip bandwidth, as the required on-chip resource usage grows at a much faster rate and would bound the sorting performance in turn. Second, designers need an in-depth understanding of HBM characteristics to effectively utilize the HBM bandwidth. To tackle these challenges, we present TopSort, a novel two-phase sorting solution optimized for HBM-based FPGAs. In the first phase, 16 merge trees work in parallel to fully utilize 32 HBM channels. In the second phase, TopSort reuses the logic from phase one to form a wider merge tree to merge the partially sorted results from phase one. TopSort also adopts HBM-specific optimizations to reduce resource overhead and improve bandwidth utilization. TopSort can sort up to 4 GB data using all 32 HBM channels, with an overall sorting performance of 15.6 GB/s. TopSort is 6.7x and 2.2x faster than state-of-the-art CPU and FPGA sorters.


title: Enhanced entanglement negativity in boundary driven monitored fermionic chains

abstract: We investigate entanglement dynamics in continuously monitored open quantum systems featuring current-carrying non-equilibrium states. We focus on a prototypical one-dimensional model of boundary-driven non-interacting fermions with monitoring of the local density, whose average Lindblad dynamics features a well-studied ballistic to diffusive crossover in transport. Here we analyze the dynamics of the fermionic negativity, mutual information, and purity along different quantum trajectories. We show that monitoring this boundary-driven system enhances its entanglement negativity at long times, which otherwise decays to zero in absence of measurements. This result is in contrast with the case of unitary evolution where monitoring suppresses entanglement production. For small values of $\gamma$, the stationary-state negativity shows a logarithmic scaling with system size, transitioning to an area-law scaling as $\gamma$ is increased beyond a critical value. Similar critical behavior is found in the mutual information, while the late-time purity shows no apparent signature of a transition, being $O(1)$ for all values of $\gamma$. Our work unveils the double role of weak monitoring in current-driven open quantum systems, simultaneously damping transport and enhancing entanglement.


title: Generalizable Task Planning through Representation Pretraining

abstract: The ability to plan for multi-step manipulation tasks in unseen situations is crucial for future home robots. But collecting sufficient experience data for end-to-end learning is often infeasible in the real world, as deploying robots in many environments can be prohibitively expensive. On the other hand, large-scale scene understanding datasets contain diverse and rich semantic and geometric information. But how to leverage such information for manipulation remains an open problem. In this paper, we propose a learning-to-plan method that can generalize to new object instances by leveraging object-level representations extracted from a synthetic scene understanding dataset. We evaluate our method with a suite of challenging multi-step manipulation tasks inspired by household activities and show that our model achieves measurably better success rate than state-of-the-art end-to-end approaches. Additional information can be found at https://sites.google.com/view/gentp


title: Joint cardiac $T_1$ mapping and cardiac function estimation using a deep manifold framework

abstract: In this work, we proposed a continuous-acquisition strategy using a gradient echo (GRE) inversion recovery sequence based on spiral trajectories to simultaneously obtain the $T_1$ mapping and CINE imaging. The acquisition is using a free-breathing and ungated fashion. An approach based on variational auto-encoder(VAE) is used for the motion estimation from the centered k-space data. The motion signal is then used to train a deep manifold reconstruction algorithm for image reconstruction. Once the network is trained, we can excite the latent vectors (the estimated motion signals and the contrast signal) in any way as we wanted to generate the image frames in the time series. We can estimate the $T_1$ mapping using the generated image frames where only contrast is varying. We can also generate the breath-hold CINE in different contrast.


title: Tapered Photonic Switching

abstract: The advent of novel nonlinear materials has stirred unprecedented interest in exploring the use of temporal inhomogeneities to achieve novel forms of wave control, amidst the greater vision of engineering metamaterials across both space and time. When the properties of an unbounded medium are abruptly switched in time, propagating waves are efficiently converted to different frequencies, and partially coupled to their back-propagating phase-conjugate partners, through a process called time-reversal. However, in realistic materials the switching time is necessarily finite, playing a central role in the resulting temporal scattering features. By identifying and leveraging the crucial role of electromagnetic momentum conservation in time-reversal processes, here we develop a very general analytical formalism to quantify time-reversal due to temporal inhomogeneities of arbitrary profile. Finally, we deploy our analytic theory to develop a formalism, analogous to spatial tapering theory, that enables the tailoring of a desired time-reversal spectral response, demonstrating its use for the realization of broadband frequency converters and filters.


title: Horizons: Nuclear Astrophysics in the 2020s and Beyond

abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.


title: Coherent light scattering from a telecom C-band quantum dot

abstract: Quantum networks have the potential to transform secure communication via quantum key distribution and enable novel concepts in distributed quantum computing and sensing. Coherent quantum light generation at telecom wavelengths is fundamental for fibre-based network implementations, but Fourier-limited emission and subnatural linewidth photons have so far only been reported from systems operating in the visible to near-infrared wavelength range. Here, we use InAs/InP quantum dots to demonstrate photons with coherence times much longer than the Fourier limit at telecom wavelength. Evidence of the responsible elastic laser scattering mechanism is observed in a distinct signature in two-photon interference measurements, and is confirmed using a direct measurement of the emission coherence. Further, we show that even the inelastically scattered photons have coherence times within the error bars of the Fourier limit. Finally, we make direct use of the minimal attenuation in fibre for these photons by measuring two-photon interference after 25 km of fibre, thereby demonstrating indistinguishability of photons emitted about 100 000 excitation cycles apart.


title: A Faber-Krahn inequality for wavelet transforms

abstract: For some special window functions $\psi_{\beta} \in H^2(\mathbb{C}^+),$ we prove that, over all sets $\Delta \subset \mathbb{C}^+$ of fixed hyperbolic measure $\nu(\Delta),$ the ones over which the Wavelet transform $W_{\overline{\psi_{\beta}}}$ with window $\overline{\psi_{\beta}}$ concentrates optimally are exactly the discs with respect to the pseudohyperbolic metric of the upper half space. This answers a question raised by Abreu and D\"orfler. Our techniques make use of a framework recently developed in a previous work by F. Nicola and the second author, but in the hyperbolic context induced by the dilation symmetry of the Wavelet transform. This leads us naturally to use a hyperbolic rearrangement function, as well as the hyperbolic isoperimetric inequality, in our analysis.


title: An Exponentially Increasing Step-size for Parameter Estimation in Statistical Models

abstract: Using gradient descent (GD) with fixed or decaying step-size is a standard practice in unconstrained optimization problems. However, when the loss function is only locally convex, such a step-size schedule artificially slows GD down as it cannot explore the flat curvature of the loss function. To overcome that issue, we propose to exponentially increase the step-size of the GD algorithm. Under homogeneous assumptions on the loss function, we demonstrate that the iterates of the proposed \emph{exponential step size gradient descent} (EGD) algorithm converge linearly to the optimal solution. Leveraging that optimization insight, we then consider using the EGD algorithm for solving parameter estimation under both regular and non-regular statistical models whose loss function becomes locally convex when the sample size goes to infinity. We demonstrate that the EGD iterates reach the final statistical radius within the true parameter after a logarithmic number of iterations, which is in stark contrast to a \emph{polynomial} number of iterations of the GD algorithm in non-regular statistical models. Therefore, the total computational complexity of the EGD algorithm is \emph{optimal} and exponentially cheaper than that of the GD for solving parameter estimation in non-regular statistical models while being comparable to that of the GD in regular statistical settings. To the best of our knowledge, it resolves a long-standing gap between statistical and algorithmic computational complexities of parameter estimation in non-regular statistical models. Finally, we provide targeted applications of the general theory to several classes of statistical models, including generalized linear models with polynomial link functions and location Gaussian mixture models.


title: Non-agency interventions for causal mediation in the presence of intermediate confounding

abstract: Recent approaches to causal inference have focused on causal effects defined as contrasts between the distribution of counterfactual outcomes under hypothetical interventions on the nodes of a graphical model. In this article we develop theory for causal effects defined with respect to a different type of intervention, one which alters the information propagated through the edges of the graph. These information transfer interventions may be more useful than node interventions in settings in which causes are non-manipulable, for example when considering race or genetics as a causal agent. Furthermore, information transfer interventions allow us to define path-specific decompositions which are identified in the presence of treatment-induced mediator-outcome confounding, a practical problem whose general solution remains elusive. We prove that the proposed effects provide valid statistical tests of mechanisms, unlike popular methods based on randomized interventions on the mediator. We propose efficient non-parametric estimators for a covariance version of the proposed effects, using data-adaptive regression coupled with semi-parametric efficiency theory to address model misspecification bias while retaining $\sqrt{n}$-consistency and asymptotic normality. We illustrate the use of our methods in two examples using publicly available data.


title: Towards Debiasing Translation Artifacts

abstract: Cross-lingual natural language processing relies on translation, either by humans or machines, at different levels, from translating training data to translating test sets. However, compared to original texts in the same language, translations possess distinct qualities referred to as translationese. Previous research has shown that these translation artifacts influence the performance of a variety of cross-lingual tasks. In this work, we propose a novel approach to reducing translationese by extending an established bias-removal technique. We use the Iterative Null-space Projection (INLP) algorithm, and show by measuring classification accuracy before and after debiasing, that translationese is reduced at both sentence and word level. We evaluate the utility of debiasing translationese on a natural language inference (NLI) task, and show that by reducing this bias, NLI accuracy improves. To the best of our knowledge, this is the first study to debias translationese as represented in latent embedding space.


title: Lost in Compression: the Impact of Lossy Image Compression on Variable Size Object Detection within Infrared Imagery

abstract: Lossy image compression strategies allow for more efficient storage and transmission of data by encoding data to a reduced form. This is essential enable training with larger datasets on less storage-equipped environments. However, such compression can cause severe decline in performance of deep Convolution Neural Network (CNN) architectures even when mild compression is applied and the resulting compressed imagery is visually identical. In this work, we apply the lossy JPEG compression method with six discrete levels of increasing compression {95, 75, 50, 15, 10, 5} to infrared band (thermal) imagery. Our study quantitatively evaluates the affect that increasing levels of lossy compression has upon the performance of characteristically diverse object detection architectures (Cascade-RCNN, FSAF and Deformable DETR) with respect to varying sizes of objects present in the dataset. When training and evaluating on uncompressed data as a baseline, we achieve maximal mean Average Precision (mAP) of 0.823 with Cascade R-CNN across the FLIR dataset, outperforming prior work. The impact of the lossy compression is more extreme at higher compression levels (15, 10, 5) across all three CNN architectures. However, re-training models on lossy compressed imagery notably ameliorated performances for all three CNN models with an average increment of ~76% (at higher compression level 5). Additionally, we demonstrate the relative sensitivity of differing object areas {tiny, small, medium, large} with respect to the compression level. We show that tiny and small objects are more sensitive to compression than medium and large objects. Overall, Cascade R-CNN attains the maximal mAP across most of the object area categories.


title: Revisit the periodicity of SGR J1935+2154 bursts with updated sample

abstract: Since FRB 200428 has been found to be associated with an X-ray burst from the Galactic magnetar SGR J1935+2154, it is interesting to explore whether the magnetar bursts also follow the similar active periodic behavior as some repeating FRBs. Previous studies show that there is possible period about 230 day in SGR J1935+2154 bursts. Here, we collected an updated burst sample from SGR J1935+2154, including all bursts reported by Fermi/GBM and GECAM till 2022 January. We also developed a targeted search pipeline to reveal more bursts from SGR J1935+2154 in the Fermi/GBM data from 2008 August to 2014 December (i.e. before the first burst detected by Swift/BAT). With this burst sample, we re-analyzed the possible periodicity of SGR J1935+2154 bursts using the Period Folding and Lomb-Scargle Periodogram methods. Our results show that the periodicity $\sim$238 day reported in literature is probably fake and the observation effects may introduce false periods (i.e. 55 day) according to simulation tests. We find that, for the current burst sample, the most probable period is 126.88$\pm$2.05 day, which could be interpreted as the precession of the magnetar. However, we note that the whole burst history is very complicated and difficult to be perfectly accommodated with any period reported thus far, therefore more monitoring observations of SGR J1935+2154 are required to test any periodicity hypothesis.


title: Multifarious skyrmion phases on a trilayer triangular lattice

abstract: The instability toward a magnetic skyrmion crystal in centrosymmetric trilayer magnets is investigated based on a spin model with layer-dependent Dzayloshinskii-Moriya interaction. We find various types of skyrmion crystal phases with different skyrmion numbers in a low-temperature phase diagram by performing the simulated annealing. In addition to the N\'eel skyrmion crystal phase that is expected to emerge in the presence of the polar-type Dzayloshinskii-Moriya interaction, we obtain the skyrmion crystal phases characteristics of the layered system: the twisted surface skyrmion crystal, anti-skyrmion crystal, and high-topological-number skyrmion crystal phases. The rich magnetic phases are brought about by the synergy among the layer-dependent Dzayloshinskii-Moriya interaction, interlayer exchange interaction, and an external magnetic field. Our results indicate that the layer degree of freedom at the surface and heterostructures provides a good platform to engineer and design the topological spin textures.


title: Controlling Viscous Fingering Instabilities of Complex Fluids

abstract: The process of one fluid pushing another is universally common while involving complex interfacial instabilities. Particularly, occurring in a myriad of natural and industrial processes, wavy fingering patterns frequently emerge when a less viscous fluid pushes another more viscous one, such as water invading oil, in a porous medium. Such finger-shaped interfaces producing partial displacement significantly affect the efficiency of numerous applications, for example, chromatography, printing devices, coating flows, oil-well cementing, as well as large-scale technologies of groundwater and enhanced oil recovery (EOR). This classical viscous fingering instability is notoriously difficult to control because the two fluids' viscosity or mobility ratio is often fixed and yet the predominant drive of the instability. Although some strategies have been recently revealed for simple fluids of constant viscosity, the feasibility of controlling the fundamental viscous fingering instability for omnipresent complex fluids has not been established. Here, we demonstrate how to control a common complex fluid (of a power-law fluid with a yield-stress) using a narrow tapered cell theoretically and experimentally.


title: A Novel Buckle-Free Large Rib Microdisk with Sub-Micron Thickness

abstract: Thin large microdisks, that are key for dense spectral microcomb generation at visible to UV wavelengths, face challenges in fabrication. One of the most difficult issues is the buckling effect that significantly reduces the cavity optical quality factor. This work introduces a novel rib disk structure that significantly mitigates the buckling effects. Using this approach, we obtained millimeter size buckle-free microdisks with sub-micron thickness and high optical quality factor exceeding $10^7$.


title: Perceptual Evaluation on Audio-visual Dataset of 360 Content

abstract: To open up new possibilities to assess the multimodal perceptual quality of omnidirectional media formats, we proposed a novel open source 360 audiovisual (AV) quality dataset. The dataset consists of high-quality 360 video clips in equirectangular (ERP) format and higher-order ambisonic (4th order) along with the subjective scores. Three subjective quality experiments were conducted for audio, video, and AV with the procedures detailed in this paper. Using the data from subjective tests, we demonstrated that this dataset can be used to quantify perceived audio, video, and audiovisual quality. The diversity and discriminability of subjective scores were also analyzed. Finally, we investigated how our dataset correlates with various objective quality metrics of audio and video. Evidence from the results of this study implies that the proposed dataset can benefit future studies on multimodal quality evaluation of 360 content.


title: Control of Viscous Fingering Instability for Complex Yield-Stress Fluids using a Tapered Cell

abstract: Being a major limiting factor for the efficiency of various technologies, such as Enhanced Oil Recovery, the viscous fingering (or Saffman--Taylor) instability has been extensively studied, especially for simple Newtonian fluids. Here, we experimentally and theoretically demonstrate a vital control of inhibiting the viscous fingering instability for complex (yield-stress) fluids to generate a complete sweep with a flat interface. Using a rectangular tapered cell, we first experimentally show the feasibility of controlling the primary fingering instability of a complex yield-stress fluid when pushed by another less viscous one. We further develop a theoretical linear stability analysis generalized for complex fluids with a yield stress and a power-law form of viscosity to provide insights. With three complex solutions yielding different viscosity contrasts, we observe stable flat and unstable wavy interfaces depending on the gap gradient ($\alpha$) and injection flow rate ($Q$). Finally, the comparison reveals an agreeable theoretical stability criterion capable of predicting stable vs. unstable displacements for yield-stress fluids under various $\alpha$.


title: GaiaHub: A method for combining data from the Gaia and Hubble space telescopes to derive improved proper motions for faint stars

abstract: We present GaiaHub, a publicly available tool that combines $Gaia$ measurements with $Hubble$ $Space$ $Telescope$ ($HST$) archival images to derive proper motions (PMs). It increases the scientific impact of both observatories beyond their individual capabilities. $Gaia$ provides PMs across the whole sky, but the limited mirror size and time baseline restrict the best PM performance to relatively bright stars. $HST$ can measure accurate PMs for much fainter stars over a small field, but this requires two epochs of observation which are not always available. GaiaHub yields considerably improved PM accuracy compared to $Gaia$-only measurements, especially for faint sources $(G \gtrsim 18)$, requiring only a single epoch of $HST$ data observed more than $\sim 7$ years ago (before 2012). This provides considerable scientific value especially for dynamical studies of stellar systems or structures in and beyond the Milky Way (MW) halo, for which the member stars are generally faint. To illustrate the capabilities and demonstrate the accuracy of GaiaHub, we apply it to samples of MW globular clusters (GCs) and classical dwarf spheroidal (dSph) satellite galaxies. This allows us, e.g., to measure the velocity dispersions in the plane of the sky for objects out to and beyond $\sim 100$ kpc. We find, on average, mild radial velocity anisotropy in GCs, consistent with existing results for more nearby samples. We observe a correlation between the internal kinematics of the clusters and their ellipticity, with more isotropic clusters being, on average, more round. Our results also support previous findings that Draco and Sculptor dSph galaxies appear to be radially anisotropic systems.


title: The e-value and the Full Bayesian Significance Test: Logical Properties and Philosophical Consequences

abstract: This article gives a conceptual review of the e-value, ev(H|X) -- the epistemic value of hypothesis H given observations X. This statistical significance measure was developed in order to allow logically coherent and consistent tests of hypotheses, including sharp or precise hypotheses, via the Full Bayesian Significance Test (FBST). Arguments of analysis allow a full characterization of this statistical test by its logical or compositional properties, showing a mutual complementarity between results of mathematical statistics and the logical desiderata lying at the foundations of this theory.


title: Level Constrained First Order Methods for Function Constrained Optimization

abstract: We present a new feasible proximal gradient method for constrained optimization where both the objective and constraint functions are given by the summation of a smooth, possibly nonconvex function and a convex simple function. The algorithm converts the original problem into a sequence of convex subproblems. Formulating those subproblems requires the evaluation of at most one gradient value of the original objective and constraint functions. Either exact or approximate subproblem solutions can be computed efficiently in many cases. An important feature of the algorithm is the constraint level parameter. By carefully increasing this level for each subproblem, we provide a simple solution to overcome the challenge of bounding the Lagrangian multipliers and show that the algorithm follows a strictly feasible solution path till convergence to the stationary point. We develop a simple, proximal gradient descent type analysis, showing that the complexity bound of this new algorithm is comparable to gradient descent for the unconstrained setting, which is new in the literature. Exploiting this new design and analysis technique, we extend our algorithms to some more challenging constrained optimization problems where 1) the objective is a stochastic or finite-sum function, and 2) structured nonsmooth functions replace smooth components of both objective and constraint functions. Complexity results for these problems also seem to be new in the literature. Finally, our method can also be applied to convex function-constrained problems where we show complexities similar to the proximal gradient method.


title: CascadER: Cross-Modal Cascading for Knowledge Graph Link Prediction

abstract: Knowledge graph (KG) link prediction is a fundamental task in artificial intelligence, with applications in natural language processing, information retrieval, and biomedicine. Recently, promising results have been achieved by leveraging cross-modal information in KGs, using ensembles that combine knowledge graph embeddings (KGEs) and contextual language models (LMs). However, existing ensembles are either (1) not consistently effective in terms of ranking accuracy gains or (2) impractically inefficient on larger datasets due to the combinatorial explosion problem of pairwise ranking with deep language models. In this paper, we propose a novel tiered ranking architecture CascadER to maintain the ranking accuracy of full ensembling while improving efficiency considerably. CascadER uses LMs to rerank the outputs of more efficient base KGEs, relying on an adaptive subset selection scheme aimed at invoking the LMs minimally while maximizing accuracy gain over the KGE. Extensive experiments demonstrate that CascadER improves MRR by up to 9 points over KGE baselines, setting new state-of-the-art performance on four benchmarks while improving efficiency by one or more orders of magnitude over competitive cross-modal baselines. Our empirical analyses reveal that diversity of models across modalities and preservation of individual models' confidence signals help explain the effectiveness of CascadER, and suggest promising directions for cross-modal cascaded architectures. Code and pretrained models are available at https://github.com/tsafavi/cascader.


title: Continual learning on 3D point clouds with random compressed rehearsal

abstract: Contemporary deep neural networks offer state-of-the-art results when applied to visual reasoning, e.g., in the context of 3D point cloud data. Point clouds are important datatype for precise modeling of three-dimensional environments, but effective processing of this type of data proves to be challenging. In the world of large, heavily-parameterized network architectures and continuously-streamed data, there is an increasing need for machine learning models that can be trained on additional data. Unfortunately, currently available models cannot fully leverage training on additional data without losing their past knowledge. Combating this phenomenon, called catastrophic forgetting, is one of the main objectives of continual learning. Continual learning for deep neural networks has been an active field of research, primarily in 2D computer vision, natural language processing, reinforcement learning, and robotics. However, in 3D computer vision, there are hardly any continual learning solutions specifically designed to take advantage of point cloud structure. This work proposes a novel neural network architecture capable of continual learning on 3D point cloud data. We utilize point cloud structure properties for preserving a heavily compressed set of past data. By using rehearsal and reconstruction as regularization methods of the learning process, our approach achieves a significant decrease of catastrophic forgetting compared to the existing solutions on several most popular point cloud datasets considering two continual learning settings: when a task is known beforehand, and in the challenging scenario of when task information is unknown to the model.


title: Accented Speech Recognition: Benchmarking, Pre-training, and Diverse Data

abstract: Building inclusive speech recognition systems is a crucial step towards developing technologies that speakers of all language varieties can use. Therefore, ASR systems must work for everybody independently of the way they speak. To accomplish this goal, there should be available data sets representing language varieties, and also an understanding of model configuration that is the most helpful in achieving robust understanding of all types of speech. However, there are not enough data sets for accented speech, and for the ones that are already available, more training approaches need to be explored to improve the quality of accented speech recognition. In this paper, we discuss recent progress towards developing more inclusive ASR systems, namely, the importance of building new data sets representing linguistic diversity, and exploring novel training approaches to improve performance for all users. We address recent directions within benchmarking ASR systems for accented speech, measure the effects of wav2vec 2.0 pre-training on accented speech recognition, and highlight corpora relevant for diverse ASR evaluations.


title: Similarity Suppresses Cyclicity: Why Similar Competitors Form Hierarchies

abstract: Competitive systems can exhibit both hierarchical (transitive) and cyclic (intransitive) structures. Despite theoretical interest in cyclic competition, which offers richer dynamics, and occupies a larger subset of the space of possible competitive systems, most real-world systems are predominantly transitive. Why? Here, we introduce a generic mechanism which promotes transitivity, even when there is ample room for cyclicity. Consider a competitive system where outcomes are mediated by competitor attributes via a performance function. We demonstrate that, if competitive outcomes depend smoothly on competitor attributes, then similar competitors compete transitively. We quantify the rate of convergence to transitivity given the similarity of the competitors and the smoothness of the performance function. Thus, we prove the adage regarding apples and oranges. Similar objects admit well ordered comparisons. Diverse objects may not. To test that theory, we run a series of evolution experiments designed to mimic genetic training algorithms. We consider a series of canonical bimatrix games and an ensemble of random performance functions that demonstrate the generality of our mechanism, even when faced with highly cyclic games. We vary the training parameters controlling the evolution process, and the shape parameters controlling the performance function, to evaluate the robustness of our results. These experiments illustrate that, if competitors evolve to optimize performance, then their traits may converge, leading to transitivity.


title: Addressable Superconductor Integrated Circuit Memory from Delay Lines

abstract: Recent advances in logic schemes and fabrication processes have renewed interest in using superconductor electronics for energy-efficient computing and quantum control processors. However, scalable superconducting memory still poses a challenge. To address this issue, we present an alternative to approaches that solely emphasize storage cell miniaturization by exploiting the minimal attenuation and dispersion properties of superconducting passive transmission lines to develop a delay-line memory system. This fully superconducting design operates at speeds between 20 GHz and 100 GHz, with $\pm$24\% and $\pm$13\% bias margins, respectively, and demonstrates data densities in the 10s of Mbit/cm$^2$ with the MIT Lincoln Laboratory SC2 fabrication process. Additionally, the circulating nature of this design allows for minimal control circuitry, eliminates the need for data splitting and merging, and enables inexpensive implementations of sequential access and content-addressable memories. Further advances in fabrication processes suggest data densities of 100s of Mbit/cm$^2$ and beyond


title: $\mathscr{H}$-Consistency Estimation Error of Surrogate Loss Minimizers

abstract: We present a detailed study of estimation errors in terms of surrogate loss estimation errors. We refer to such guarantees as $\mathscr{H}$-consistency estimation error bounds, since they account for the hypothesis set $\mathscr{H}$ adopted. These guarantees are significantly stronger than $\mathscr{H}$-calibration or $\mathscr{H}$-consistency. They are also more informative than similar excess error bounds derived in the literature, when $\mathscr{H}$ is the family of all measurable functions. We prove general theorems providing such guarantees, for both the distribution-dependent and distribution-independent settings. We show that our bounds are tight, modulo a convexity assumption. We also show that previous excess error bounds can be recovered as special cases of our general results. We then present a series of explicit bounds in the case of the zero-one loss, with multiple choices of the surrogate loss and for both the family of linear functions and neural networks with one hidden-layer. We further prove more favorable distribution-dependent guarantees in that case. We also present a series of explicit bounds in the case of the adversarial loss, with surrogate losses based on the supremum of the $\rho$-margin, hinge or sigmoid loss and for the same two general hypothesis sets. Here too, we prove several enhancements of these guarantees under natural distributional assumptions. Finally, we report the results of simulations illustrating our bounds and their tightness.


title: A Survey on Semantics in Automated Data Science

abstract: Data Scientists leverage common sense reasoning and domain knowledge to understand and enrich data for building predictive models. In recent years, we have witnessed a surge in tools and techniques for {\em automated machine learning}. While data scientists can employ various such tools to help with model building, many other aspects such as {\em feature engineering} that require semantic understanding of concepts, remain manual to a large extent. In this paper we discuss important shortcomings of current automated data science solutions and machine learning. We discuss how leveraging basic semantic reasoning on data in combination with novel tools for data science automation can help with consistent and explainable data augmentation and transformation. Moreover, semantics can assist data scientists in a new manner by helping with challenges related to {\em trust}, {\em bias}, and {\em explainability}.


title: Viscous Fingering Instability of Complex Fluids in a Tapered Geometry

abstract: Viscous fingering (VF) is an interfacial instability that occurs in a narrow confinement or porous medium when a less-viscous fluid pushes a more viscous one, producing finger-like patterns. Controlling the VF instability is essential to enhance the efficiency of various technological applications. However, the control of VF instability has been challenging and so far focused on simple Newtonian fluids of constant viscosity. Here, we extend to complex yield-stress fluids and examine the controlling feasibility by carrying out a linear stability analysis using a radial cell with a converging gap gradient. We avoid making the major assumption of a small Bingham number, Bn << 1, i.e., a negligible ratio of the yield to shear stress, and instead provide a new stability criterion predicting apparent complex VF. This criterion depends on not only the complex fluid's rheology, interfacial tension, and contact angle to the wetting wall, but also the gap gradient, the radius, gap-thickness, and velocity at the fluid-fluid interface. Finally, we compare this theoretical criterion to our experimental data with nitrogen pushing a complex yield-stress fluid in a taper and find good agreement.


title: Partial Product Aware Machine Learning on DNA-Encoded Libraries

abstract: DNA encoded libraries (DELs) are used for rapid large-scale screening of small molecules against a protein target. These combinatorial libraries are built through several cycles of chemistry and DNA ligation, producing large sets of DNA-tagged molecules. Training machine learning models on DEL data has been shown to be effective at predicting molecules of interest dissimilar from those in the original DEL. Machine learning chemical property prediction approaches rely on the assumption that the property of interest is linked to a single chemical structure. In the context of DNA-encoded libraries, this is equivalent to assuming that every chemical reaction fully yields the desired product. However, in practice, multi-step chemical synthesis sometimes generates partial molecules. Each unique DNA tag in a DEL therefore corresponds to a set of possible molecules. Here, we leverage reaction yield data to enumerate the set of possible molecules corresponding to a given DNA tag. This paper demonstrates that training a custom GNN on this richer dataset improves accuracy and generalization performance.


title: On the Homology stability range for symplectic groups

abstract: We improve, by a factor of 2, known homology stability ranges for the integral homology of symplectic groups over commutative local rings with infinite residue field and show that the obstruction to further stability is bounded below by Milnor-Witt K-theory. In particular our stability range is optimal in many cases.


title: A faster algorithm for Vertex Cover parameterized by solution size

abstract: We describe a new algorithm for vertex cover with runtime $O^*(1.25284^k)$, where $k$ is the size of the desired solution and $O^*$ hides polynomial factors in the input size. This improves over previous runtime of $O^*(1.2738^k)$ due to Chen, Kanj, & Xia (2010) standing for more than a decade. The key to our algorithm is to use a potential function which simultaneously tracks $k$ as well as the optimal value $\lambda$ of the vertex cover LP relaxation. This approach also allows us to make use of prior algorithms for Maximum Independent Set in bounded-degree graphs and Above-Guarantee Vertex Cover. The main step in the algorithm is to branch on high-degree vertices, while ensuring that both $k$ and $\mu = k - \lambda$ are decreased at each step. There can be local obstructions in the graph that prevent $\mu$ from decreasing in this process; we develop a number of novel branching steps to handle these situations.


title: A stationary model of non-intersecting directed polymers

abstract: We consider the partition function $Z_{\ell}(\vec x,0\vert \vec y,t)$ of $\ell$ non-intersecting continuous directed polymers of length $t$ in dimension $1+1$, in a white noise environment, starting from positions $\vec x$ and terminating at positions $\vec y$. When $\ell=1$, it is well known that for fixed $x$, the field $\log Z_1(x,0\vert y,t)$ solves the Kardar-Parisi-Zhang equation and admits the Brownian motion as a stationary measure. In particular, as $t$ goes to infinity, $Z_1(x,0\vert y,t)/Z_1(x,0\vert 0,t) $ converges to the exponential of a Brownian motion $B(y)$. In this article, we show an analogue of this result for any $\ell$. We show that $Z_{\ell}(\vec x,0\vert \vec y,t)/Z_{\ell}(\vec x,0\vert \vec 0,t) $ converges as $t$ goes to infinity to an explicit functional $Z_{\ell}^{\rm stat}(\vec y)$ of $\ell$ independent Brownian motions. This functional $Z_{\ell}^{\rm stat}(\vec y)$ admits a simple description as the partition sum for $\ell$ non-intersecting semi-discrete polymers on $\ell$ lines. We discuss applications to the endpoints and midpoints distribution for long non-crossing polymers and derive explicit formulas in the case of two polymers. To obtain these results, we show that the stationary measure of the O'Connell-Warren multilayer stochastic heat equation is given by a collection of independent Brownian motions. This in turn is shown via analogous results in a discrete setup for the so-called log-gamma polymer and exploit the connection between non-intersecting log-gamma polymers and the geometric RSK correspondence found in arXiv:1110.3489. .

categories: cond-mat.stat-mech math-ph math.MP math.PR


title: Harnessing Multilingual Resources to Question Answering in Arabic

abstract: The goal of the paper is to predict answers to questions given a passage of Qur'an. The answers are always found in the passage, so the task of the model is to predict where an answer starts and where it ends. As the initial data set is rather small for training, we make use of multilingual BERT so that we can augment the training data by using data available for languages other than Arabic. Furthermore, we crawl a large Arabic corpus that is domain specific to religious discourse. Our approach consists of two steps, first we train a BERT model to predict a set of possible answers in a passage. Finally, we use another BERT based model to rank the candidate answers produced by the first BERT model.


title: The Hamiltonian Path Graph is Connected for Simple $s,t$ Paths in Rectangular Grid Graphs

abstract: A \emph{simple} $s,t$ path $P$ in a rectangular grid graph $\mathbb{G}$ is a Hamiltonian path from the top-left corner $s$ to the bottom-right corner $t$ such that each \emph{internal} subpath of $P$ with both endpoints $a$ and $b$ on the boundary of $\mathbb{G}$ has the minimum number of bends needed to travel from $a$ to $b$ (i.e., $0$, $1$, or $2$ bends, depending on whether $a$ and $b$ are on opposite, adjacent, or the same side of the bounding rectangle). Here, we show that $P$ can be reconfigured to any other simple $s,t$ path of $\mathbb{G}$ by \emph{switching $2\times 2$ squares}, where at most ${5}|\mathbb{G}|/{4}$ such operations are required. Furthermore, each \emph{square-switch} is done in $O(1)$ time and keeps the resulting path in the same family of simple $s,t$ paths. Our reconfiguration result proves that the \emph{Hamiltonian path graph} $\cal{G}$ for simple $s,t$ paths is connected and has diameter at most ${5}|\mathbb{G}|/{4}$ which is asymptotically tight.


title: Shaving the Hair of Black Hole with Sagittarius A$^*$ from Event Horizon Telescope

abstract: Recently, the Event Horizon Telescope collaboration has reported the first image of the supermassive black hole Sagittarius in the Galactic Center. We attempt to test the validity of the no-hair theorem of black holes using this new shadow observation. Considering the Einstein-Maxwell-klein-Gordon theory with a minimally-coupled scalar field, we find that our numerical result is consistent with the prediction of the no-hair theorem. However, we can not rule out the possibility that black holes with scalar hair may exist in some special cases.


title: Energy transfer in quantum molecular chain -- two models of inhomogeneity

abstract: We study a linear chain of oscillators with inhomogeneity in their interactions with phonon bath. In a previous work on the Markovian master equation of the system, we investigated a model in which the difference in the site-phonon coupling between adjacent oscillators is the same throughout the chain. Here we look into another model in which the oscillators are coupled to the phonon bath with alternating strength at successive sites. Whereas in the first model all exciton modes are connected, in the second model they are coupled in pairs that are not connected to each other. Owing to this special structure in the coupling, the excitation numbers of different modes can be solved exactly in the steady state. In the first model, the minima of the excitation profile in the site basis occur at the edges of the chain, whereas in the second model the maxima occur at the edges. The energy transfer efficiency in the first model is affected by the source power whereas in the second model the efficiency is independent of it. A distinct feature in the second model is that a sink placed at the middle of the chain is able to distinguish between chains with even and odd number of sites. The energy transfer efficiency in a chain with even number of sites is higher than a chain with odd number of sites. Therefore, it reveals the discrete nature of the chain. In the limit of very long chain when the discreteness of the chain is less evident, the efficiencies approach each other.


title: Browser-based Hyperbolic Visualization of Graphs

abstract: Hyperbolic geometry offers a natural focus + context for data visualization and has been shown to underlie real-world complex networks. However, current hyperbolic network visualization approaches are limited to special types of networks and do not scale to large datasets. With this in mind, we designed, implemented, and analyzed three methods for hyperbolic visualization of networks in the browser based on inverse projections, generalized force-directed algorithms, and hyperbolic multi-dimensional scaling (H-MDS). A comparison with Euclidean MDS shows that H-MDS produces embeddings with lower distortion for several types of networks. All three methods can handle node-link representations and are available in fully functional web-based systems.


title: Automatic Error Classification and Root Cause Determination while Replaying Recorded Workload Data at SAP HANA

abstract: Capturing customer workloads of database systems to replay these workloads during internal testing can be beneficial for software quality assurance. However, we experienced that such replays can produce a large amount of false positive alerts that make the results unreliable or time consuming to analyze. Therefore, we design a machine learning based approach that attributes root causes to the alerts. This provides several benefits for quality assurance and allows for example to classify whether an alert is true positive or false positive. Our approach considerably reduces manual effort and improves the overall quality assurance for the database system SAP HANA. We discuss the problem, the design and result of our approach, and we present practical limitations that may require further research.


title: Interpretable sensitivity analysis for the Baron-Kenny approach to mediation with unmeasured confounding

abstract: Mediation analysis assesses the extent to which the exposure affects the outcome indirectly through a mediator and the extent to which it operates directly through other pathways. As the most popular method in empirical mediation analysis, the Baron-Kenny approach estimates the indirect and direct effects of the exposure on the outcome based on linear structural equation models. However, when the exposure and the mediator are not randomized, the estimates may be biased due to unmeasured confounding among the exposure, mediator, and outcome. Building on Cinelli and Hazlett (2020), we derive general omitted-variable bias formulas in linear regressions with vector responses and regressors. We then use the formulas to develop a sensitivity analysis method for the Baron-Kenny approach to mediation in the presence of unmeasured confounding. To ensure interpretability, we express the sensitivity parameters to correspond to the natural factorization of the joint distribution of the direct acyclic graph for mediation analysis. They measure the partial correlation between the unmeasured confounder and the exposure, mediator, outcome, respectively. With the sensitivity parameters, we propose a novel measure called the "robustness value for mediation" or simply the "robustness value", to assess the robustness of results based on the Baron-Kenny approach with respect to unmeasured confounding. Intuitively, the robustness value measures the minimum value of the maximum proportion of variability explained by the unmeasured confounding, for the exposure, mediator and outcome, to overturn the results of the point estimate or confidence interval for the direct and indirect effects. Importantly, we prove that all our sensitivity bounds are attainable and thus sharp.


title: Thermodynamics of Quantum Measurement and the Demon's Arrow of Time

abstract: We discuss the thermodynamic aspects of a single qubit based device, powered by weak quantum measurements, and feedback controlled by a quantum Maxwell's demon. We discuss both discrete and time-continuous operation of the measurement based device at finite temperature of the reservoir. In the discrete example where a demon acquires information via discrete weak measurements, we find that the thermodynamic variables including the heat exchanged, extractable work, and the entropy produced are completely determined by an information theoretic measure of the demon's perceived arrow of time. We also discuss a realistic time-continuous operation of the device where the feedback is applied after a sequence of weak measurements. In the time-continuous limit, we derive the exact finite-time statistics of work, heat and entropy changes along individual quantum trajectories of the quantum measurement process, and relate them to the demon's arrow of time.


title: On Algebraic Constructions of Neural Networks with Small Weights

abstract: Neural gates compute functions based on weighted sums of the input variables. The expressive power of neural gates (number of distinct functions it can compute) depends on the weight sizes and, in general, large weights (exponential in the number of inputs) are required. Studying the trade-offs among the weight sizes, circuit size and depth is a well-studied topic both in circuit complexity theory and the practice of neural computation. We propose a new approach for studying these complexity trade-offs by considering a related algebraic framework. Specifically, given a single linear equation with arbitrary coefficients, we would like to express it using a system of linear equations with smaller (even constant) coefficients. The techniques we developed are based on Siegel's Lemma for the bounds, anti-concentration inequalities for the existential results and extensions of Sylvester-type Hadamard matrices for the constructions. We explicitly construct a constant weight, optimal size matrix to compute the EQUALITY function (checking if two integers expressed in binary are equal). Computing EQUALITY with a single linear equation requires exponentially large weights. In addition, we prove the existence of the best-known weight size (linear) matrices to compute the COMPARISON function (comparing between two integers expressed in binary). In the context of the circuit complexity theory, our results improve the upper bounds on the weight sizes for the best-known circuit sizes for EQUALITY and COMPARISON.


title: Using Embeddings for Causal Estimation of Peer Influence in Social Networks

abstract: We address the problem of using observational data to estimate peer contagion effects, the influence of treatments applied to individuals in a network on the outcomes of their neighbors. A main challenge to such estimation is that homophily - the tendency of connected units to share similar latent traits - acts as an unobserved confounder for contagion effects. Informally, it's hard to tell whether your friends have similar outcomes because they were influenced by your treatment, or whether it's due to some common trait that caused you to be friends in the first place. Because these common causes are not usually directly observed, they cannot be simply adjusted for. We describe an approach to perform the required adjustment using node embeddings learned from the network itself. The main aim is to perform this adjustment nonparametrically, without functional form assumptions on either the process that generated the network or the treatment assignment and outcome processes. The key contributions are to nonparametrically formalize the causal effect in a way that accounts for homophily, and to show how embedding methods can be used to identify and estimate this effect. Code is available at https://github.com/IrinaCristali/Peer-Contagion-on-Networks.


title: DeepSim: A Reinforcement Learning Environment Build Toolkit for ROS and Gazebo

abstract: We propose DeepSim, a reinforcement learning environment build toolkit for ROS and Gazebo. It allows machine learning or reinforcement learning researchers to access the robotics domain and create complex and challenging custom tasks in ROS and Gazebo simulation environments. This toolkit provides building blocks of advanced features such as collision detection, behaviour control, domain randomization, spawner, and many more. DeepSim is designed to reduce the boundary between robotics and machine learning communities by providing Python interface. In this paper, we discuss the components and design decisions of DeepSim Toolkit.


title: Global solutions with infinitely many blowups in a mean-field neural network

abstract: We recently introduced idealized mean-field models for networks of integrate-and-fire neurons with impulse-like interactions -- the so-called delayed Poissonian mean-field models. Such models are prone to blowups: for a strong enough interaction coupling, the mean-field rate of interaction diverges in finite time with a finite fraction of neurons spiking simultaneously. Due to the reset mechanism of integrate-and-fire neurons, these blowups can happen repeatedly, at least in principle. A benefit of considering Poissonian mean-field models is that one can resolve blowups analytically by mapping the original singular dynamics onto uniformly regular dynamics via a time change. Resolving a blowup then amounts to solving the fixed-point problem that implicitly defines the time change, which can be done consistently for a single blowup and for nonzero delays. Here we extend this time-change analysis in two ways: First, we exhibit the existence and uniqueness of explosive solutions with a countable infinity of blowups in the large interaction regime. Second, we show that these delayed solutions specify "physical" explosive solutions in the limit of vanishing delays, which in turn can be explicitly constructed. The first result relies on the fact that blowups are self-sustaining but nonoverlapping in the time-changed picture. The second result follows from the continuity of blowups in the time-changed picture and incidentally implies the existence of periodic solutions. These results are useful to study the emergence of synchrony in neural network models.


title: On Semiparametric Efficiency of an Emerging Class of Regression Models for Between-subject Attributes

abstract: The semiparametric regression models have attracted increasing attention owing to their robustness compared to their parametric counterparts. This paper discusses the efficiency bound for functional response models (FRM), an emerging class of semiparametric regression that serves as a timely solution for research questions involving pairwise observations. This new paradigm is especially appealing to reduce astronomical data dimensions for those arising from wearable devices and high-throughput technology, such as microbiome Beta-diversity, viral genetic linkage, single-cell RNA sequencing, etc. Despite the growing applications, the efficiency of their estimators has not been investigated carefully due to the extreme difficulty to address the inherent correlations among pairs. Leveraging the Hilbert-space-based semiparametric efficiency theory for classical within-subject attributes, this manuscript extends such asymptotic efficiency into the broader regression involving between-subject attributes and pinpoints the most efficient estimator, which leads to a sensitive signal-detection in practice. With pairwise outcomes burgeoning immensely as effective dimension-reduction summaries, the established theory will not only fill the critical gap in identifying the most efficient semiparametric estimator but also propel wide-ranging implementations of this new paradigm for between-subject attributes.


title: A review on non-Hermitian skin effect

abstract: The past decades have witnessed the flourishing of non-Hermitian physics in non-conservative systems, leading to unprecedented phenomena of unidirectional invisibility, enhanced sensitivity and more recently the novel topological features such as bulk Fermi arcs. Among them, growing efforts have been invested to an intriguing phenomenon, known as the non-Hermitian skin effect (NHSE). Here, we review the recent progress in this emerging field. By starting from the one-dimensional (1D) case, the fundamental concepts of NHSE, its minimal model, the physical meanings and consequences are elaborated in details. In particular, we discuss the NHSE enriched by lattice symmetries, which gives rise to unique non-Hermitian topological properties with revised bulk-boundary correspondence (BBC) and new definitions of topological invariants. Then we extend the discussions to two and higher dimensions, where dimensional surprises enable even more versatile NHSE phenomena. Extensions of NHSE assisted with extra degrees of freedom such as long-range coupling, pseudospins, magnetism, non-linearity and crystal defects are also reviewed. This is followed by the contemporary experimental progress for NHSE. Finally, we provide the outlooks to possible future directions and developments.


title: Newton and interior-point methods for (constrained) nonconvex-nonconcave minmax optimization with stability and instability guarantees

abstract: We address the problem of finding a local solution to a nonconvex-nonconcave minmax optimization using Newton type methods, including interior-point ones. We modify the Hessian matrix of these methods such that, at each step, the modified Newton update direction can be seen as the solution to a quadratic program that locally approximates the minmax problem. Moreover, we show that by selecting the modification in an appropriate way, the only stable equilibrium points of the algorithm's iterations are local minmax points. As a consequence, the algorithm can only converge towards an equilibrium point if such point is a local minmax, and it will escape if the point is not a local minmax. Using numerical examples, we show that the computation time of our algorithm scales roughly linearly with the number of nonzero elements in the Hessian.


title: Duality Cascades and Parallelotopes

abstract: Duality cascades are a series of duality transformations in field theories, which can be realized as the Hanany-Witten transitions in brane configurations on a circle. In the setup of the ABJM theory and its generalizations, from the physical requirement that duality cascades always end and the final destination depends only on the initial brane configuration, we propose that the fundamental domain of supersymmetric brane configurations in duality cascades can tile the whole parameter space of relative ranks by translations, hence is a parallelotope. We provide our arguments for the proposal.


title: Zariski dense surface subgroups in $SL(n,\mathbb{Q})$ with odd $n$

abstract: For odd $n$ we construct a path $\rho_t\colon \pi_1(S) \to SL(n,\mathbb{R})$ of discrete, faithful and Zariski dense representations of a surface group such that $\rho_t(\pi_1(S)) \subset SL(n,\mathbb{Q})$ for every $t\in \mathbb{Q}$.


title: Detection and Physical Interaction with Deformable Linear Objects

abstract: Deformable linear objects (e.g., cables, ropes, and threads) commonly appear in our everyday lives. However, perception of these objects and the study of physical interaction with them is still a growing area. There have already been successful methods to model and track deformable linear objects. However, the number of methods that can automatically extract the initial conditions in non-trivial situations for these methods has been limited, and they have been introduced to the community only recently. On the other hand, while physical interaction with these objects has been done with ground manipulators, there have not been any studies on physical interaction and manipulation of the deformable linear object with aerial robots. This workshop describes our recent work on detecting deformable linear objects, which uses the segmentation output of the existing methods to provide the initialization required by the tracking methods automatically. It works with crossings and can fill the gaps and occlusions in the segmentation and output the model desirable for physical interaction and simulation. Then we present our work on using the method for tasks such as routing and manipulation with the ground and aerial robots. We discuss our feasibility analysis on extending the physical interaction with these objects to aerial manipulation applications.


title: The Impact of the Social Security Reforms on Welfare: Who benefits and Who loses across Generations, Gender, and Employment Type?

abstract: We quantitatively explore the impact of social security reforms in Japan, which is facing rapid aging and the highest government debt among developed countries, using an overlapping generations model with four types of agents distinguished by gender and employment type. We find that introducing social security reforms without extending the retirement age raises the welfare of future generations, while reforms with rising copayment rates for medical and long-term care expenditures, in particular, significantly lowers the welfare of low-income groups (females and part-timers) of the current retired and working generations. In contrast, reforms reducing the pension replacement rate lead to a greater decline in the welfare of full-timers. The combination of these reforms and the extension of the retirement age is expected to improve the welfare of the current working generations by 2--9 % over the level without reforms.


title: Explainable and Optimally Configured Artificial Neural Networks for Attack Detection in Smart Homes

abstract: In recent years cybersecurity has become a major concern in adaptation of smart applications. Specially, in smart homes where a large number of IoT devices are used having a secure and trusted mechanisms can provide peace of mind for users. Accurate detection of cyber attacks is crucial, however precise identification of the type of attacks plays a huge role if devising the countermeasure for protecting the system. Artificial Neural Networks (ANN) have provided promising results for detecting any security attacks for smart applications. However, due to complex nature of the model used for this technique it is not easy for normal users to trust ANN based security solutions. Also, selection of right hyperparameters for ANN architecture plays a crucial role in the accurate detection of security attacks, especially when it come to identifying the subcategories of attacks. In this paper, we propose a model that considers both the issues of explainability of ANN model and the hyperparameter selection for this approach to be easily trusted and adapted by users of smart home applications. Also, our approach considers a subset of the dataset for optimal selection of hyperparamters to reduce the overhead of the process of ANN architecture design. Distinctively this paper focuses on configuration, performance and evaluation of ANN architecture for identification of five categorical attacks and nine subcategorical attacks. Using a very recent IoT dataset our approach showed high performance for intrusion detection with 99.9%, 99.7%, and 97.7% accuracy for Binary, Category, and Subcategory level classification of attacks.


title: Modules at boundary points, fiberwise Bergman kernels, and log-subharmonicity II -- on Stein manifolds

abstract: In this article, we consider Bergman kernels related to modules at boundary points on Stein manifolds, and obtain a log-subharmonicity property of the Bergman kernels. As applications, we obtain a lower estimate of weighted $L^2$ integrals on Stein manifolds, and reprove an effectiveness result of strong openness property of modules at boundary points on Stein manifolds.


title: Enhanced $d$-$p$ hybridization intertwined with anomalous ground state formation in van der Waals-coupled magnetic metal Fe$_5$GeTe$_2$

abstract: Fe$_5$GeTe$_2$ is a van der Waals (vdW)-coupled unconventional ferromagnetic metal with a high Curie temperature ($T_C$) exceeding 300 K. The formation of an anomalous ground state significantly below $T_C$ has received considerable attention, resulting in increased interest in understanding the spin-polarized electronic state evolution near the Fermi energy ($E_F$) as a function of temperature. Despite recent extensive studies, a microscopic understanding of the spin-polarized electronic structure around $E_F$ has not yet been established owing to the intrinsic complexity of both the crystal and band structures. In this study, we investigate the temperature dependence of element-specific soft X-ray magnetic circular dichroism (XMCD). A systematic temperature evolution in the XMCD signal from both magnetic Fe and its ligand Te is clearly observed. More importantly, the enhancement in the hybridization between the Fe 3$d$ and Te 5$p$ states in the zero-magnetic field limit is revealed, and we discuss its implications on the possible emergence of an exotic magnetic ground state in Fe$_5$GeTe$_2$.


title: Shape complexity in cluster analysis

abstract: In cluster analysis, a common first step is to scale the data aiming to better partition them into clusters. Even though many different techniques have throughout many years been introduced to this end, it is probably fair to say that the workhorse in this preprocessing phase has been to divide the data by the standard deviation along each dimension. Like division by the standard deviation, the great majority of scaling techniques can be said to have roots in some sort of statistical take on the data. Here we explore the use of multidimensional shapes of data, aiming to obtain scaling factors for use prior to clustering by some method, like k-means, that makes explicit use of distances between samples. We borrow from the field of cosmology and related areas the recently introduced notion of shape complexity, which in the variant we use is a relatively simple, data-dependent nonlinear function that we show can be used to help with the determination of appropriate scaling factors. Focusing on what might be called "midrange" distances, we formulate a constrained nonlinear programming problem and use it to produce candidate scaling-factor sets that can be sifted on the basis of further considerations of the data, say via expert knowledge. We give results on some iconic data sets, highlighting the strengths and potential weaknesses of the new approach. These results are generally positive across all the data sets used.


title: Perfect Spectral Clustering with Discrete Covariates

abstract: Among community detection methods, spectral clustering enjoys two desirable properties: computational efficiency and theoretical guarantees of consistency. Most studies of spectral clustering consider only the edges of a network as input to the algorithm. Here we consider the problem of performing community detection in the presence of discrete node covariates, where network structure is determined by a combination of a latent block model structure and homophily on the observed covariates. We propose a spectral algorithm that we prove achieves perfect clustering with high probability on a class of large, sparse networks with discrete covariates, effectively separating latent network structure from homophily on observed covariates. To our knowledge, our method is the first to offer a guarantee of consistent latent structure recovery using spectral clustering in the setting where edge formation is dependent on both latent and observed factors.


title: A Short Introduction to the Koopman Representation of Dynamical Systems

abstract: The Koopman representation is an infinite dimensional linear representation of linear or nonlinear dynamical systems. It represents the dynamics of output maps (aka observables), which are functions on the state space whose evaluation is interpreted as an output. Conceptually simple derivations and commentary on the Koopman representation are given. We emphasize an important duality between initial conditions and output maps of the original system, and those of the Koopman representation. This duality is an important consideration when this representation is used in data-driven applications such as the Dynamic Mode Decomposition (DMD) and its variants. The adjoint relation between the Koopman representation and the transfer operator of mass transport is also shown.


title: Boundary layer problem on chemotaxis-Navier-Stokes system with Robin boundary conditions

abstract: This paper is concerned with the boundary layer problem on a chemotaxis-Navier-Stokes system modelling boundary layer formation of aerobic bacteria in fluid. Completing the system with physical Robin-type boundary conditions for oxygen, no-flux and Dirichlet boundary conditions for bacteria and fluid velocity, we show that the gradients of its radial solutions in a region between two concentric spheres possessing boundary layer effects as the oxygen diffusion rate $\varepsilon$ goes to zero and the boundary-layer thickness is of order $\mathcal{O}(\varepsilon^\alpha)$ with $0<\alpha<\frac{1}{2}$.


title: On the solidity parameter in canopy flows

abstract: We have performed high-fidelity simulations of turbulent open-channel flows over submerged rigid canopies made of cylindrical filaments of fixed length $l=0.25H$ ($H$ being the domain depth) mounted on the wall with an angle of inclination $\theta$. The inclination is the free parameter that sets the density of the canopy by varying its frontal area. The density of the canopy, based on the solidity parameter $\lambda$, is a widely accepted criterion defining the ongoing canopy flow regime, with low values ($\lambda \ll 0.1$) indicating the sparse regime, and higher values ($\lambda > 0.1$) the dense regime. All the numerical predictions have been obtained considering the same nominal bulk Reynolds number (i.e. $Re_b=U_b H/\nu = 6000$). We consider nine configurations of canopies, with $\theta$ varying symmetrically around $0{\deg}$ in the range $\theta\in [\pm 78.5{\deg}]$, where positive angles define canopies inclined in the flow direction (with the grain) and $\theta=0{\deg}$ corresponds to the wall-normally mounted canopy. The study compares canopies with identical solidity obtained inclining the filaments in opposite angles and assesses the efficacy of the solidity as a representative parameter. It is found that when the canopy is inclined, the actual flow regime differs substantially from the one of a straight canopy that shares the same solidity indicating that criteria solely based on this parameter are not robust. Finally, a new phenomenological model describing the interaction between the coherent structures populating the canopy region and the outer flow is given.


title: Gyrohydrodynamics: Relativistic spinful fluid with strong vorticity

abstract: We develop a relativistic (quasi-)hydrodynamic framework, dubbed the gyrohydrodynamics, to describe fluid dynamics of many-body systems with spin under strong vorticity based on entropy-current analysis. This framework generalizes the recently-developed spin hydrodynamics to the regime where the spin density is at the leading order in derivatives but suppressed by another small parameter, the Planck constant $\hbar$, due to its quantum nature. Our analysis shows that the complete first-order constitutive relations of gyrohydrodynamics involve seventeen transport coefficients and are highly anisotropic.


title: An Inverse Probability Weighted Regression Method that Accounts for Right-censoring for Causal Inference with Multiple Treatments and a Binary Outcome

abstract: Comparative effectiveness research often involves evaluating the differences in the risks of an event of interest between two or more treatments using observational data. Often, the post-treatment outcome of interest is whether the event happens within a pre-specified time window, which leads to a binary outcome. One source of bias for estimating the causal treatment effect is the presence of confounders, which are usually controlled using propensity score-based methods. An additional source of bias is right-censoring, which occurs when the information on the outcome of interest is not completely available due to dropout, study termination, or treatment switch before the event of interest. We propose an inverse probability weighted regression-based estimator that can simultaneously handle both confounding and right-censoring, calling the method CIPWR, with the letter C highlighting the censoring component. CIPWR estimates the average treatment effects by averaging the predicted outcomes obtained from a logistic regression model that is fitted using a weighted score function. The CIPWR estimator has a double robustness property such that estimation consistency can be achieved when either the model for the outcome or the models for both treatment and censoring are correctly specified. We establish the asymptotic properties of the CIPWR estimator for conducting inference, and compare its finite sample performance with that of several alternatives through simulation studies. The methods under comparison are applied to a cohort of prostate cancer patients from an insurance claims database for comparing the adverse effects of four candidate drugs for advanced stage prostate cancer.


title: Modeling the impact of spatial oxygen heterogeneity on radiolytic oxygen depletion during FLASH radiotherapy

abstract: It has been postulated that the delivery of radiotherapy at ultra-high dose rates ("FLASH") reduces normal tissue toxicities by depleting them of oxygen. The fraction of normal tissue and cancer cells surviving radiotherapy depends on dose and oxygen levels in an exponential manner and even a very small fraction of tissue at low oxygen levels can determine radiotherapy response. The effect of FLASH on radiation-induced normal and tumour tissue cell killing was studied by simulating oxygen diffusion, metabolism, and radiolytic oxygen depletion over domains with simulated capillary architectures. Two architectural models were used: 1.) randomly distributed capillaries and 2.) capillaries forming a regular square lattice array. The resulting oxygen partial pressure distribution histograms were used to simulate normal and tumour tissue cell survival using the linear quadratic model of cell survival, modified to incorporate oxygen-enhancement ratio (OER) effects. Tumour cell survival was found to be increased by FLASH as compared to conventional radiotherapy, with a 0-1 order of magnitude increase for expected levels of tumour hypoxia, depending on the relative magnitudes of radiolytic oxygen depletion and tissue oxygen metabolism. Interestingly, for the random capillary model, the impact of FLASH on well-oxygenated (normal) tissues was found to be much greater, with an estimated increase in cell survival by up to 10 orders of magnitude, even though reductions in mean tissue partial pressure were modest, less than 7 mmHg for the parameter values studied. The presence of very small nearly hypoxic regions in otherwise well-perfused normal tissues with high mean oxygen levels resulted in a greater proportional sparing of normal tissue than tumour cells during FLASH irradiation, possibly explaining empirical normal tissue sparing and iso-tumour control results.


title: Classification of anatomic structures in head and neck by CT-based radiomics

abstract: Background and Purpose: Radiomics features are used to identify disease types and predict therapy outcomes. However, how the radiomics features are different among different anatomical structures has never been investigated. Hence, we analyzed the radiomics features of 22 anatomical structures in the head and neck area in CT images. Furthermore, we studied whether CT radiomics can classify anatomical structures of the head and neck using unsupervised machine-learning techniques. Materials and methods: We obtained IMRT/VMAT treatment planning data from 36 patients treated for head and neck cancers in a single institution. There were 1357 contours of more than 22 anatomical structures drawn on planning CTs. We calculated 174 radiomics features using the SIBEX program. First, we tested whether the radiomics features of anatomical structures were unique enough to classify all contours into 22 groups. We then developed a two-stage clustering technique to classify 22 anatomic structures into sub-groups with similar physiological or biological characteristics. Results: The heatmap of 174 radiomics features of 22 anatomical structures showed a distinct difference among tumors and other healthy structures. Radiomics features have allowed us to identify the eyes, lens, submandibular, pituitary glands, and thyroids with over 90% accuracy. The two-stage clustering of 22 structures resulted in six subgroups, which shared common characteristics such as fatty and bony tissues. Conclusions: We have shown that anatomical structures in head and neck tumors have distinguishable radiomics features. We could observe similarities of features among subgroups of the structures. The results suggest that CT radiomics can help distinguish the biological characteristics of head and neck lesions.


title: HelixADMET: a robust and endpoint extensible ADMET system incorporating self-supervised knowledge transfer

abstract: Accurate ADMET (an abbreviation for "absorption, distribution, metabolism, excretion, and toxicity") predictions can efficiently screen out undesirable drug candidates in the early stage of drug discovery. In recent years, multiple comprehensive ADMET systems that adopt advanced machine learning models have been developed, providing services to estimate multiple endpoints. However, those ADMET systems usually suffer from weak extrapolation ability. First, due to the lack of labelled data for each endpoint, typical machine learning models perform frail for the molecules with unobserved scaffolds. Second, most systems only provide fixed built-in endpoints and cannot be customised to satisfy various research requirements. To this end, we develop a robust and endpoint extensible ADMET system, HelixADMET (H-ADMET). H-ADMET incorporates the concept of self-supervised learning to produce a robust pre-trained model. The model is then fine-tuned with a multi-task and multi-stage framework to transfer knowledge between ADMET endpoints, auxiliary tasks, and self-supervised tasks. Our results demonstrate that H-ADMET achieves an overall improvement of 4%, compared with existing ADMET systems on comparable endpoints. Additionally, the pre-trained model provided by H-ADMET can be fine-tuned to generate new and customised ADMET endpoints, meeting various demands of drug research and development requirements.


title: "What makes a question inquisitive?" A Study on Type-Controlled Inquisitive Question Generation

abstract: We propose a type-controlled framework for inquisitive question generation. We annotate an inquisitive question dataset with question types, train question type classifiers, and finetune models for type-controlled question generation. Empirical results demonstrate that we can generate a variety of questions that adhere to specific types while drawing from the source texts. We also investigate strategies for selecting a single question from a generated set, considering both an informative vs.~inquisitive question classifier and a pairwise ranker trained from a small set of expert annotations. Question selection using the pairwise ranker yields strong results in automatic and manual evaluation. Our human evaluation assesses multiple aspects of the generated questions, finding that the ranker chooses questions with the best syntax (4.59), semantics (4.37), and inquisitiveness (3.92) on a scale of 1-5, even rivaling the performance of human-written questions.


title: Machine-learning-accelerated Bose-Einstein condensation

abstract: Machine learning is emerging as a technology that can enhance physics experiment execution and data analysis. Here, we apply machine learning to accelerate the production of a Bose-Einstein condensate (BEC) of $^{87}\mathrm{Rb}$ atoms by Bayesian optimization of up to 55 control parameters. This approach enables us to prepare BECs of $2.8 \times 10^3$ optically trapped $^{87}\mathrm{Rb}$ atoms from a room-temperature gas in 575 ms. The algorithm achieves the fast BEC preparation by applying highly efficient Raman cooling to near quantum degeneracy, followed by a brief final evaporation. We anticipate that many other physics experiments with complex nonlinear system dynamics can be significantly enhanced by a similar machine-learning approach.


title: Derivation of Two-fluid Model Based on Onsager Principle

abstract: Using Onsager variational principle, we study the dynamic coupling between the stress and the composition in polymer solution. In the original derivation of the two-fluid model [Doi and Onuki, J. Phys. II France {\bf 2}, 1631 (1992)], the polymer stress was introduced \emph{a priopri}, therefore a constitutive equation is required to close the equations. Based on our previous study of viscoelastic fluids with homogeneous composition [Phys. Rev. Fluids {\bf 3}, 084004 (2018)], we start with a dumbbell model for the polymer, and derive all dynamic equations using the Onsager variational principle.


title: Evolution strategies: Application in hybrid quantum-classical neural networks

abstract: With the rapid development of quantum computers, several applications are being proposed for them. Quantum simulations, simulation of chemical reactions, solution of optimization problems and quantum neural networks (QNNs) are some examples. However, problems such as noise, limited number of qubits and circuit depth, and gradient vanishing must be resolved before we can use them to their full potential. In the field of quantum machine learning, several models have been proposed. In general, in order to train these different models, we use the gradient of a cost function with respect to the model parameters. In order to obtain this gradient, we must compute the derivative of this function with respect to the model parameters. One of the most used methods in the literature to perform this task is the parameter-shift rule method. This method consists of evaluating the cost function twice for each parameter of the QNN. A problem with this method is that the number of evaluations grows linearly with the number of parameters. In this work we study an alternative method, called Evolution Strategies (ES), which are a family of black box optimization algorithms which iteratively update the parameters using a search gradient. An advantage of the ES method is that in using it one can control the number of times the cost function will be evaluated. We apply the ES method to the binary classification task, showing that this method is a viable alternative for training QNNs. However, we observe that its performance will be strongly dependent on the hyperparameters used. Furthermore, we also observe that this method, alike the parameter shift rule method, suffers from the problem of gradient vanishing.


title: Semi-inclusive Diffractive Deep Inelastic Scattering at Small-$x$

abstract: Inspired by a recent study of Iancu, Mueller and Triantafyllopoulos [1] and earlier papers by Golec-Biernat and Wusthoff [2,3], we propose semi-inclusive diffractive deep inelastic scattering (SIDDIS) to investigate the gluon tomography in the nucleon and nuclei at small-$x$. The relevant diffractive quark and gluon parton distribution functions (DPDF) can be computed in terms of the color dipole S-matrices in the fundamental and adjoint representations, respectively. Novel correlations from the gluon tomography in the dipole S-matrix can be experimentally studied through the DPDFs in these processes at the future electron-ion collider (EIC).


title: Impact of Economic Constraints on the Projected Timeframe for Human-Crewed Deep Space Exploration

abstract: Deep space exploration offers the most profound opportunity for the expansion of humanity and our understanding of the Universe, but remains extremely challenging. Progress will continue to be paced by uncrewed missions followed up by crewed missions to ever further destinations. Major space powers continue to invest in crewed deep space exploration as an important national strategy. An improved model based on previous work is developed, which projects the earliest possible launch dates for human-crewed missions from cis-lunar space to selected destinations in the Solar System and beyond based on NASA's historic budget trend and overall development trends of deep space exploration research. The purpose of the analysis is to provide a projected timeframe for crewed missions beyond Mars. Our findings suggest the first human missions from a spacefaring nation or international collaboration to the Asteroid Belt and Jovian System could be scheduled as soon as ~2071 to ~2087 and ~2101 to ~2121, respectively, while a launch to the Saturn System may occur by the year ~2132, with an uncertainty window of ~2129 to ~2153.

categories: physics.pop-ph astro-ph.EP astro-ph.IM physics.soc-ph


title: Strong Revenue (Non-)Monotonicity of Single-parameter Auctions

abstract: Consider Myerson's optimal auction with respect to an inaccurate prior, e.g., estimated from data, which is an underestimation of the true value distribution. Can the auctioneer expect getting at least the optimal revenue w.r.t. the inaccurate prior since the true value distribution is larger? This so-called strong revenue monotonicity is known to be true for single-parameter auctions when the feasible allocations form a matroid. We find that strong revenue monotonicity fails to generalize beyond the matroid setting, and further show that auctions in the matroid setting are the only downward-closed auctions that satisfy strong revenue monotonicity. On the flip side, we recover an approximate version of strong revenue monotonicity that holds for all single-parameter auctions, even without downward-closedness. As applications, we get sample complexity upper bounds for single-parameter auctions under matroid constraints, downward-closed constraints, and general constraints. They improve the state-of-the-art upper bounds and are tight up to logarithmic factors.


title: Fast consensus of high-order multi-agent systems

abstract: In this paper, the fast consensus problem of high-order multi-agent systems under undirected topologies is considered. The direct link between the consensus convergence rate and the control gains is established. An accelerated consensus algorithm based on gradient descent is proposed to optimize the convergence rate. By applying the Routh-Hurwitz stability criterion, the lower bound on the convergence rate is derived, and explicit control gains are derived as the necessary condition to achieve the optimal convergence rate. Moreover, a protocol with time-varying control gains is designed to achieve the finite-time consensus. Explicit formulas for the time-varying control gains and the final consensus state are given. Numerical examples and simulation results are presented to illustrate the obtained theoretical results.


title: Pattern Propagation Driven by Surface Curvature

abstract: Pattern dynamics on curved surfaces are found everywhere in nature. The geometry of surfaces have been shown to influence dynamics and play a functional role, yet a comprehensive understanding is still elusive. Here, we report for the first time that a static Turing pattern on a flat surface can propagate on a curved surface, as opposed to previous studies, where the pattern is presupposed to be static irrespective of the surface geometry. To understand such significant changes on curved surfaces, we investigate reaction-diffusion systems on axisymmetric curved surfaces. Numerical and theoretical analyses reveal that both the symmetries of the surface and pattern participate in the initiation of pattern propagation. This study provides a novel and generic mechanism of pattern propagation that is caused by surface curvature, as well as insights into the general role of surface geometry.