title: Accurate Clinical Toxicity Prediction using Multi-task Deep Neural Nets and Contrastive Molecular Explanations

abstract: Explainable ML for molecular toxicity prediction is a promising approach for efficient drug development and chemical safety. A predictive ML model of toxicity can reduce experimental cost and time while mitigating ethical concerns by significantly reducing animal and clinical testing. Herein, we use a deep learning framework for simultaneously modeling in vitro, in vivo, and clinical toxicity data. Two different molecular input representations are used: Morgan fingerprints and pre-training SMILES embeddings. A multi-task deep learning model accurately predicts toxicity for all endpoints, including clinical, as indicated by AUROC and balanced accuracy. In particular, SMILES embeddings as input to the multi-task model improved clinical toxicity predictions compared to existing models in MoleculeNet benchmark. Additionally, our multi-task approach is comprehensive in the sense that it is comparable to state-of-the-art approaches for specific endpoints in in vitro, in vivo and clinical platforms. Through both the multi-task model and transfer learning, we were able to indicate the minimal need of in vivo data for clinical toxicity predictions. To provide confidence and explain the model's predictions, we adapt a post-hoc contrastive explanation method that returns pertinent positive and pertinent negative features, which correspond well to known mutagenic and reactive toxicophores, such as unsubstituted bonded heteroatoms, aromatic amines, and Michael receptors. Furthermore, toxicophore recovery by pertinent feature analysis captures more of the in vitro (53%) and in vivo (56%), rather than of the clinical (8%), endpoints, and indeed uncovers a preference in known toxicophore data towards in vitro and in vivo experimental data. To our knowledge, this is the first contrastive explanation, using both present and absent substructures, for predictions of clinical and in vivo molecular toxicity.


title: Polarised random k-SAT

abstract: In this paper we study a variation of the random $k$-SAT problem, called polarized random $k$-SAT. In this model there is a polarization parameter $p$, and in half of the clauses each variable occurs negated with probability $p$ and pure otherwise, while in the other half the probabilities are interchanged. For $p=1/2$ we get the classical random $k$-SAT model, and at the other extreme we have the fully polarized model where $p=0$, or $1$. Here there are only two types of clauses: clauses where all $k$ variables occur pure, and clauses where all $k$ variables occur negated. That is, for $p=0$ we get an instance of random monotone $k$-SAT. We show that the threshold of satisfiability does not decrease as $p$ moves away from $\frac{1}{2}$ and thus that the satisfiability threshold for polarized random $k$-SAT is an upper bound on the threshold for random $k$-SAT. In fact, we conjecture that asymptotically the two thresholds coincide.


title: Predicting score distribution to improve non-intrusive speech quality estimation

abstract: Deep noise suppressors (DNS) have become an attractive solution to remove background noise, reverberation, and distortions from speech and are widely used in telephony/voice applications. They are also occasionally prone to introducing artifacts and lowering the perceptual quality of the speech. Subjective listening tests that use multiple human judges to derive a mean opinion score (MOS) are a popular way to measure these models' performance. Deep neural network based non-intrusive MOS estimation models have recently emerged as a popular cost-efficient alternative to these tests. These models are trained with only the MOS labels, often discarding the secondary statistics of the opinion scores. In this paper, we investigate several ways to integrate the distribution of opinion scores (e.g. variance, histogram information) to improve the MOS estimation performance. Our model is trained on a corpus of 419K denoised samples by 320 different DNS models and model variations and evaluated on 18K test samples from DNSMOS. We show that with very minor modification of a single task MOS estimation pipeline, these freely available labels can provide up to a 0.016 RMSE and 1% SRCC improvement.


title: Model knot solutions for the twisted Bogomolny equations

abstract: In this paper we prove existence for model knot solutions to the dimensionally reduced twisted Kapustin-Witten equations on $\mathbb{R}^3_+$ for any twisting parameter $t\in(0,\infty)$. We start with the explicit solutions for $t = 1$ derived in \cite{WFivebranes} and perform a continuity argument in $t$. This corroborates a prediction of Gaiotto and Witten \cite[][p. 961]{GW}.


title: Formal Language Recognition by Hard Attention Transformers: Perspectives from Circuit Complexity

abstract: This paper analyzes three formal models of Transformer encoders that differ in the form of their self-attention mechanism: unique hard attention (UHAT); generalized unique hard attention (GUHAT), which generalizes UHAT; and averaging hard attention (AHAT). We show that UHAT and GUHAT Transformers, viewed as string acceptors, can only recognize formal languages in the complexity class AC$^0$, the class of languages recognizable by families of Boolean circuits of constant depth and polynomial size. This upper bound subsumes Hahn's (2020) results that GUHAT cannot recognize the DYCK languages or the PARITY language, since those languages are outside AC$^0$ (Furst et al., 1984). In contrast, the non-AC$^0$ languages MAJORITY and DYCK-1 are recognizable by AHAT networks, implying that AHAT can recognize languages that UHAT and GUHAT cannot.


title: The heat kernel on the diagonal for a compact metric graph

abstract: We analyze the heat kernel associated to the Laplacian on a compact metric graph, with standard Kirchoff-Neumann vertex conditions. An explicit formula for the heat kernel as a sum over loops, developed by Roth and Kostrykin, Potthoff, and Schrader, allows for a straightforward analysis of small-time asymptotics. We show that the restriction of the heat kernel to the diagonal satisfies a modified version of the heat equation. This observation leads to an "edge" heat trace formula, expressing the a sum over eigenfunction amplitudes on a single edge as a sum over closed loops containing that edge. The proof of this formula relies on a modified heat equation satisfied by the diagonal restriction of the heat kernel. Further study of this equation leads to explicit formulas for completely symmetric graphs.


title: Cartesian convexity as the key notion in the variational existence theory for nonlocal supremal functionals

abstract: Motivated by the direct method in the calculus of variations in $L^{\infty}$, our main result identifies the notion of convexity characterizing the weakly$^*$ lower semicontinuity of nonlocal supremal functionals: Cartesian level convexity. This new concept coincides with separate level convexity in the one-dimensional setting and is strictly weaker for higher dimensions. We discuss relaxation in the vectorial case, showing that the relaxed functional will not generally maintain the supremal form. Apart from illustrating this fact with examples of multi-well type, we present precise criteria for structure-preservation. When the structure is preserved, a representation formula is given in terms of the Cartesian level convex envelope of the (diagonalized) original supremand. This work does not only complete the picture of the analysis initiated in [Kreisbeck \& Zappale, Calc.~Var.~PDE, 2020], but also establishes a connection with double integrals. We relate the two classes of functionals via an $L^p$-approximation in the sense of $\Gamma$-convergence for diverging integrability exponents. The proofs exploit recent results on nonlocal inclusions and their asymptotic behavior, and use tools from Young measure theory and convex analysis.


title: Crystalline prisms: Reflections and diffractions, present and past

abstract: Let $Y/S$ be a $p$-completely smooth morphism of $p$-torsion free $p$-adic formal schemes endowed with a Frobenius lift, and let $\overline Y/\overline S$ denote its reduction modulo $p$. We show that the category of crystals on the prismatic site of $\overline Y/S$ is equivalent to the category of $O_Y$-modules with integrable and quasi-nilpotent $p$-connection, and that the cohomology of such a crystal is computed by the associated $p$-de Rham complex. More generally, if $X$ is a closed subscheme of $\overline Y$, smooth over $\overline S$, then the prismatic envelope $\Delta_X(Y)$ of $X$ in $Y$ admits such a $p$-connection, the category of prismatic crystals on $X/S$ is equivalent to the category of $O_ {\Delta_X(Y)}$-modules with compatible integrable and quasi-nilpotent $p$-connection, and the cohomology of such a crystal is again computed by its $p$-de Rham complex. We also give a geometric construction of the ``prismatic Sen operator.'' Namely, we show that a lifting of $X$ (mod $p^2$) in $Y$ defines a vector field on the reduction modulo $p$ of $\Delta_X(Y)$ and on a ``diffracted'' Higgs complex which calculates the mod $p$ prismatic and de Rham cohomologies of $X$. Surprisingly, this complex is not the reduction modulo $p$ of the afore-mentioned $p$-de Rham complexbut is rather its ``$\alpha$-transform.'' As a consequence, we get a fairly explicit description of the action of the group scheme $G^\gamma$ on $R\Gamma(X, \Omega^\bullet_{X/S})$, Drinfeld's strengthening of the Deligne-Illusie decomposition theorem. We also explain how earlier work by several authors relating Higgs fields, $p$-connections, and connections can be placed in the prismatic context.


title: PDE-constrained optimization for electroencephalographic source reconstruction

abstract: This paper introduces a novel numerical method for the inverse problem of electroencephalography(EEG). We pose the inverse EEG problem as an optimal control (OC) problem for Poisson's equation. The optimality conditions lead to a variational system of differential equations. It is discretized directly in finite-element spaces leading to a system of linear equations with a sparse Karush-Kuhn-Tucker matrix. The method uses finite-element discretization and thus can handle MRI-based meshes of almost arbitrary complexity. It extends the well-known mixed quasi-reversibility method (mQRM) in that pointwise noisy data explicitly appear in the formulation making unnecessary tedious interpolation of the noisy data from the electrodes to the scalp surface. The resulting algebraic problem differs considerably from that obtained in the mixed quasi-reversibility, but only slightly larger. An interesting feature of the algorithm is that it does not require the formation of the lead-field matrix. Our tests, both with spherical and MRI-based meshes, demonstrates that the method accurately reconstructs cortical activity.


title: Lags of the KiloHertz Quasi-Periodic Oscillations in the transient source XTE J1701$-$462

abstract: We analysed 14 observations with kilohertz quasi-periodic oscillations (kHz QPOs) of the neutron star X-ray binary XTE J1701$-$462, the first source to show a clear transition between atoll and Z-like behaviour during a single outburst. We calculated the average cross-spectrum of both atoll and Z-phase observations of XTE J1701$-$462 between a reference/hard band (6.1 - 25.7 keV) and a subject/soft band (2.1 - 5.7 keV) to obtain, using a novel technique, the average time lags of the lower and upper kHz QPOs. During the atoll phase, we found that at the frequency of the lower kHz QPO the soft photons lag behind the hard ones by $18 \pm 8$ $\mu$s, whereas during the Z phase the lags are $33\pm35$ $\mu$s, consistent with zero. This difference in the lags of both phases suggests that in XTE J1701$-$462, as observed in other sources, the lags decrease with increasing luminosity. We found that for both the atoll and Z phase observations the fractional rms amplitude increases with energy up to $\sim$10 keV and remains more or less constant at higher energies. Since these changes in the variability of XTE J1701$-$462 occur within the same outburst, properties like the mass of the neutron star or the inclination of the system cannot be responsible for the differences in the timing properties of the kHz QPOs in the atoll and Z phase. Here we suggest that these differences are driven by a Comptonizing component or corona, possibly oscillating in a coupled mode with the innermost regions of the accretion disc.


title: A Natural Language Processing Approach for Instruction Set Architecture Identification

abstract: Binary analysis of software is a critical step in cyber forensics applications such as program vulnerability assessment and malware detection. This involves interpreting instructions executed by software and often necessitates converting the software's binary file data to assembly language. The conversion process requires information about the binary file's target instruction set architecture (ISA). However, ISA information might not be included in binary files due to compilation errors, partial downloads, or adversarial corruption of file metadata. Machine learning (ML) is a promising methodology that can be used to identify the target ISA using binary data in the object code section of binary files. In this paper we propose a binary code feature extraction model to improve the accuracy and scalability of ML-based ISA identification methods. Our feature extraction model can be used in the absence of domain knowledge about the ISAs. Specifically, we adapt models from natural language processing (NLP) to i) identify successive byte patterns commonly observed in binary codes, ii) estimate the significance of each byte pattern to a binary file, and iii) estimate the relevance of each byte pattern in distinguishing between ISAs. We introduce character-level features of encoded binaries to identify fine-grained bit patterns inherent to each ISA. We use a dataset with binaries from 12 different ISAs to evaluate our approach. Empirical evaluations show that using our byte-level features in ML-based ISA identification results in an 8% higher accuracy than the state-of-the-art features based on byte-histograms and byte pattern signatures. We observe that character-level features allow reducing the size of the feature set by up to 16x while maintaining accuracy above 97%.


title: CAMERO: Consistency Regularized Ensemble of Perturbed Language Models with Weight Sharing

abstract: Model ensemble is a popular approach to produce a low-variance and well-generalized model. However, it induces large memory and inference costs, which are often not affordable for real-world deployment. Existing work has resorted to sharing weights among models. However, when increasing the proportion of the shared weights, the resulting models tend to be similar, and the benefits of using model ensemble diminish. To retain ensemble benefits while maintaining a low memory cost, we propose a consistency-regularized ensemble learning approach based on perturbed models, named CAMERO. Specifically, we share the weights of bottom layers across all models and apply different perturbations to the hidden representations for different models, which can effectively promote the model diversity. Meanwhile, we apply a prediction consistency regularizer across the perturbed models to control the variance due to the model diversity. Our experiments using large language models demonstrate that CAMERO significantly improves the generalization performance of the ensemble model. Specifically, CAMERO outperforms the standard ensemble of 8 BERT-base models on the GLUE benchmark by 0.7 with a significantly smaller model size (114.2M vs. 880.6M).


title: Adaptive Memory Management for Video Object Segmentation

abstract: Matching-based networks have achieved state-of-the-art performance for video object segmentation (VOS) tasks by storing every-k frames in an external memory bank for future inference. Storing the intermediate frames' predictions provides the network with richer cues for segmenting an object in the current frame. However, the size of the memory bank gradually increases with the length of the video, which slows down inference speed and makes it impractical to handle arbitrary length videos. This paper proposes an adaptive memory bank strategy for matching-based networks for semi-supervised video object segmentation (VOS) that can handle videos of arbitrary length by discarding obsolete features. Features are indexed based on their importance in the segmentation of the objects in previous frames. Based on the index, we discard unimportant features to accommodate new features. We present our experiments on DAVIS 2016, DAVIS 2017, and Youtube-VOS that demonstrate that our method outperforms state-of-the-art that employ first-and-latest strategy with fixed-sized memory banks and achieves comparable performance to the every-k strategy with increasing-sized memory banks. Furthermore, experiments show that our method increases inference speed by up to 80% over the every-k and 35% over first-and-latest strategies.


title: Performance Assessment of different Machine Learning Algorithm for Life-Time Prediction of Solder Joints based on Synthetic Data

abstract: This paper proposes a computationally efficient methodology to predict the damage progression in solder contacts of electronic components using temperature-time curves. For this purpose, two machine learning algorithms, a Multilayer Perceptron and a Long Short-Term Memory network, are trained and compared with respect to their prediction accuracy and the required amount of training data. The training is performed using synthetic, normally distributed data that is realistic for automotive applications. A finite element model of a simple bipolar chip resistor in surface mount technology configuration is used to numerically compute the synthetic data. As a result, both machine learning algorithms show a relevant accuracy for the prediction of accumulated creep strains. With a training data length of 350 hours (12.5% of the available training data), both models show a constantly good fitting performance of $R^2$ of 0.72 for the Multilayer Perceptron and $R^2$ of 0.87 for the Long Short-Term Memory network. The prediction errors of the accumulated creep strains are less than 10% with an amount of 350 hours training data and decreases to less than 5 % when using further data. Therefore, both approaches are promising for the lifetime prediction directly on the electronic device.


title: Spectroscopic analysis of VVV CL001 cluster with MUSE

abstract: Like most spiral galaxies, the Milky Way contains a population of blue, metal-poor globular clusters and another of red, metal-rich ones. Most of the latter belong to the bulge, and therefore they are poorly studied compared to the blue (halo) ones because they suffer higher extinction and larger contamination from field stars. These intrinsic difficulties, together with a lack of low-mass bulge globular clusters, are reasons to believe that their census is not complete yet. Indeed, a few new clusters have been confirmed in the last few years. One of them is VVV CL001, the subject of the present study. We present a new spectroscopic analysis of the recently confirmed globular cluster VVV CL001, made by means of MUSE@VLT integral field data. Individual spectra were extracted for stars in the VVV CL001 field. Radial velocities were derived by cross-correlation with synthetic templates. Coupled with PMs from the VVV survey, these data allow us to select 55 potential cluster members, for which we derive metallicities using the public code The Cannon. The mean radial velocity of the cluster is Vhelio = -324.9 +- 0.8 km/s,as estimated from 55 cluster members. This high velocity, together with a low metallicity [Fe/H] = -2.04 +- 0.02 dex suggests that VVV CL001 could be a very old cluster. The estimated distance is d = 8.23 +- 0.46 kpc, placing the cluster in the Galactic bulge. Furthermore, both its current position and the orbital parameters suggest that VVV CL001 is most probably a bulge globular cluster.


title: A fast direct solver for two dimensional quasi-periodic multilayered media scattering problems, Part II

abstract: This manuscript is the second in a series presenting fast direct solution techniques for solving two-dimensional wave scattering problems from quasi-periodic multilayered structures. The fast direct solvers presented in the series are for the linear system that results from the discretization of a robust integral formulation. The fast direct solver presented in this manuscript has a computational cost that scales linearly with respect to the number of discretization points on the interfaces and the number of layers. The latter is an improvement over the previous solver and makes the new solver more efficient especially for problems involving multiple incident angles and changes to the layered media. Numerical results illustrate the improved performance of the new solver over the previous one.


title: Colourings of path systems

abstract: A $P_m$ path in a graph is a path on $m$ vertices. A $P_m$ system of order $n>1$ is a partition of the edges of the complete graph $K_n$ into $P_m$ paths. A $P_m$ system is said to be $k$-colourable if the vertex set of $K_n$ can be partitioned into $k$ sets called colour classes such that no path in the system is monochromatic. The system is $k$-chromatic if it is $k$-colourable but is not $(k-1)$-colourable. If every $k$-colouring of a $P_m$ system can be obtained from some $k$-colouring $\phi$ by a permutation of the colours, we say that the system is uniquely $k$-colourable. In this paper, we first observe that there exists a $k$-chromatic $P_m$ system for any $k\geq 2$ and $m\geq 4$ where $m$ is even. Next, we prove that there exists an equitably 2-chromatic $P_4$ system of order $n$ for each admissible order $n$. We then show that for all $k\geq 3$, there exists a $k$-chromatic $P_4$ system of order $n$ for all sufficiently large admissible $n$. Finally, we show that there exists a uniquely 2-chromatic $P_4$ system of order $n$ for each admissible $n \geq 109$.


title: A self-consistent analytical theory for rotator networks under stochastic forcing: effects of intrinsic noise and common input

abstract: Despite the incredible complexity of our brains' neural networks, theoretical descriptions of neural dynamics have led to profound insights into possible network states and dynamics. It remains challenging to develop theories that apply to spiking networks and thus allow one to characterize the dynamic properties of biologically more realistic networks. Here, we build on recent work by van Meegen & Lindner who have shown that "rotator networks," while considerably simpler than real spiking networks and therefore more amenable to mathematical analysis, still allow to capture dynamical properties of networks of spiking neurons. This framework can be easily extended to the case where individual units receive uncorrelated stochastic input which can be interpreted as intrinsic noise. However, the assumptions of the theory do not apply anymore when the input received by the single rotators is strongly correlated among units. As we show, in this case the network fluctuations become significantly non-Gaussian, which calls for a reworking of the theory. Using a cumulant expansion, we develop a self-consistent analytical theory that accounts for the observed non-Gaussian statistics. Our theory provides a starting point for further studies of more general network setups and information transmission properties of these networks.


title: Recurrent event analysis in the presence of real-time high frequency data via random subsampling

abstract: Digital monitoring studies collect real-time high frequency data via mobile sensors in the subjects' natural environment. This data can be used to model the impact of changes in physiology on recurrent event outcomes such as smoking, drug use, alcohol use, or self-identified moments of suicide ideation. Likelihood calculations for the recurrent event analysis, however, become computationally prohibitive in this setting. Motivated by this, a random subsampling framework is proposed for computationally efficient, approximate likelihood-based estimation. A subsampling-unbiased estimator for the derivative of the cumulative hazard enters into an approximation of log-likelihood. The estimator has two sources of variation: the first due to the recurrent event model and the second due to subsampling. The latter can be reduced by increasing the sampling rate; however, this leads to increased computational costs. The approximate score equations are equivalent to logistic regression score equations, allowing for standard, "off-the-shelf" software to be used in fitting these models. Simulations demonstrate the method and efficiency-computation trade-off. We end by illustrating our approach using data from a digital monitoring study of suicidal ideation.


title: Assessing the data-analysis impact of LISA orbit approximations using a GPU-accelerated response model

abstract: The analysis of gravitational wave (GW) datasets is based on the comparison of measured time series with theoretical templates of the detector's response to a variety of source parameters. For LISA, the main scientific observables will be the so-called time-delay interferometry (TDI) combinations, which suppress the otherwise overwhelming laser noise. Computing the TDI response to GW involves projecting the GW polarizations onto the LISA constellation arms, and then combining projections delayed by a multiple of the light propagation time along the arms. Both computations are difficult to perform efficiently for generic LISA orbits and GW signals. Various approximations are currently used in practice, e.g., assuming constant and equal armlengths, which yields analytical TDI expressions. In this article, we present 'fastlisaresponse', a new efficient GPU-accelerated code that implements the generic TDI response to GWs in the time domain. We use it to characterize the parameter-estimation bias incurred by analyzing loud Galactic-binary signals using the equal-armlength approximation. We conclude that equal-armlength parameter-estimation codes should be upgraded to the generic response if they are to achieve optimal accuracy for high (but reasonable) SNR sources within the actual LISA data.


title: Seaweed algebras

abstract: The index of a Lie algebra is an important algebraic invariant, but it is notoriously difficult to compute. However, for the suggestively-named seaweed algebras, the computation of the index can be reduced to a combinatorial formula based on the connected components of a "meander": a planar graph associated with the algebra. Our index analysis on seaweed algebras requires only basic linear and abstract algebra. Indeed, the main goal of this survey-type article is to introduce a broader audience to seaweed algebras with minimal appeal to specialized language and notation from Lie theory. This said, we present several results that do not appear elsewhere and do appeal to more advanced language in the Introduction to provide added context.


title: A Novel Approach for Optimum-Path Forest Classification Using Fuzzy Logic

abstract: In the past decades, fuzzy logic has played an essential role in many research areas. Alongside, graph-based pattern recognition has shown to be of great importance due to its flexibility in partitioning the feature space using the background from graph theory. Some years ago, a new framework for both supervised, semi-supervised, and unsupervised learning named Optimum-Path Forest (OPF) was proposed with competitive results in several applications, besides comprising a low computational burden. In this paper, we propose the Fuzzy Optimum-Path Forest, an improved version of the standard OPF classifier that learns the samples' membership in an unsupervised fashion, which are further incorporated during supervised training. Such information is used to identify the most relevant training samples, thus improving the classification step. Experiments conducted over twelve public datasets highlight the robustness of the proposed approach, which behaves similarly to standard OPF in worst-case scenarios.


title: Characterization of fractional Sobolev--Poincar\'e and (localized) Hardy inequalities

abstract: In this paper, we prove capacitary versions of the fractional Sobolev--Poincar\'e inequalities. We characterize localized variant of the boundary fractional Sobolev--Poincar\'e inequalities through uniform fatness condition of the domain in $\mathbb{R}^n$. Existence type results on the fractional Hardy inequality are established in the supercritical case $sp>n$ for $s\in(0,1)$, $p>1$. Characterization of the fractional Hardy inequality through weak supersolution of the associate problem is also addressed.


title: Nonparametric Identification of Differentiated Products Demand Using Micro Data

abstract: We examine identification of differentiated products demand when one has "micro data" linking individual consumers' characteristics and choices. Our model nests standard specifications featuring rich observed and unobserved consumer heterogeneity as well as product/market-level unobservables that introduce the problem of econometric endogeneity. Previous work establishes identification of such models using market-level data and instruments for all prices and quantities. Micro data provides a panel structure that facilitates richer demand specifications and reduces requirements on both the number and types of instrumental variables. We address identification of demand in the standard case in which non-price product characteristics are assumed exogenous, but also cover identification of demand elasticities and other key features when product characteristics are endogenous. We discuss implications of these results for applied work.


title: Clifford Circuits can be Properly PAC Learned if and only if $\textsf{RP}=\textsf{NP}$

abstract: Given a dataset of input states, measurements, and probabilities, is it possible to efficiently predict the measurement probabilities associated with a quantum circuit? Recent work of Caro and Datta (2020) studied the problem of PAC learning quantum circuits in an information theoretic sense, leaving open questions of computational efficiency. In particular, one candidate class of circuits for which an efficient learner might have been possible was that of Clifford circuits, since the corresponding set of states generated by such circuits, called stabilizer states, are known to be efficiently PAC learnable (Rocchetto 2018). Here we provide a negative result, showing that proper learning of CNOT circuits is hard for classical learners unless $\textsf{RP} = \textsf{NP}$. As the classical analogue and subset of Clifford circuits, this naturally leads to a hardness result for Clifford circuits as well. Additionally, we show that if $\textsf{RP} = \textsf{NP}$ then there would exist efficient proper learning algorithms for CNOT and Clifford circuits. By similar arguments, we also find that an efficient proper quantum learner for such circuits exists if and only if $\textsf{NP} \subseteq \textsf{RQP}$.


title: Observation of bosonic stimulation in light scattering

abstract: For bosons, the transition rate into an already occupied quantum state is enhanced by its occupation number: the effect of bosonic stimulation. Bosonic enhancement of light scattering has been predicted more than 30 years ago but has not been observed before. Here we theoretically investigate and experimentally demonstrate this effect. We show that the bosonic enhancement factor for a harmonically trapped gas is bounded by a universal constant $\zeta(2)/\zeta(3)$ above the phase transition to a Bose-Einstein condensate (BEC), and depends linearly on the BEC fraction just below the phase transition. Bosonically enhanced light scattering is observed and characterized above and below the phase transition, and the effect of interactions is discussed. For a multi-level system, bosonic enhancement is reduced because bosonic stimulation occurs only for Rayleigh scattering, but not for Raman scattering.


title: Temperature-Dependent Diffuse Reflectance Measurements of Ceramic Powders in the Near- and Mid-Infrared Spectra

abstract: This study focuses on experimentally measuring temperature-dependent diffuse reflectance in the near- and mid-infrared spectra for ceramic particles with applications as heat-transfer and thermal-storage media in concentrated solar power (CSP) plants. Specifically, a commercially available sintered bauxite ceramic powder, ACCUCAST ID80, and its primary chemical constituents, alumina (Al2O3) and silica (SiO2), are measured using a Fourier transform infrared spectrometer (FTIR) coupled with a specialized diffuse reflectance accessory and a heated stage. Room-temperature diffuse reflectance measurements show increased absorption in tests with greater mass fractions of the ceramic samples. There is a strong correlation in the measured reflectance spectra of ACCUCAST with alumina and silica in the spectral range 2000-500 cm-1. For the first time, temperature-dependent diffuse reflectance measurements are reported for ACCUCAST, including a novel technique for accessing reflectance values above the limiting temperature of the background material KBr. All three materials exhibit a calculated emittance of ~0.9 at room temperature. However, this value drops to 0.68 at 1000 C for ACCUCAST and ~0.43 for alumina and silica. Thermal cycling in air from 25 C to 1000 C resulted in a visible color change from dark grey to light orange for ACCUCAST and a subsequent 5X greater increase in reflectance at 4000 cm-1 as compared to ACCUCAST thermally cycled at 1000 C in vacuum. Alumina and silica spectra proved to be largely unaffected by thermal cycling under atmospheric and evacuated conditions.


title: Multiple higher-order topological phases with even and odd pairs of zero-energy corner modes in a $C_3$ symmetry broken model

abstract: The higher-order corner modes for quantum anomalous Hall insulators in $C_3$ symmetry broken honeycomb lattice have been engineered recently. Here we consider an extended Haldane model in presence of inversion symmetry breaking sub-lattice mass, time-reversal symmetry breaking Zeeman field and spin-orbit coupling interaction where we find that only the quantum spin Hall insulator can host the second-order dipolar phase while the remaining two first-order topological phases do not morph into the latter. Remarkably, four-fold degeneracy of zero-energy dipolar states can be reduced to two-fold under the application (withdrawn) of sub-lattice mass (Zeeman field) term when the spin-orbit coupling is already present. On the other hand, the sub-lattice mass and Zeeman field terms compete with each other to pin down the two mid-gap states at zero-energy in the absence or presence of spin-orbit coupling. Interestingly, the bulk-polarization can topologically characterize the dipolar phase irrespective of the energy of the mid-gap states as long as inversion symmetry is preserved. The effective gap criterion can qualitatively mimic the extent of SOT phase originated by the interplay between finite Zeeman exchange field, sub-lattice mass, and SOC interaction.


title: Optimal resource allocation for flexible-grid entanglement distribution networks

abstract: We use a genetic algorithm (GA) as a design aid for determining the optimal provisioning of entangled photon spectrum in flex-grid quantum networks with arbitrary numbers of channels and users. After introducing a general model for entanglement distribution based on frequency-polarization hyperentangled biphotons, we derive upper bounds on fidelity and entangled bit rate for networks comprising one-to-one user connections. Simple conditions based on user detector quality and link efficiencies are found that determine whether entanglement is possible. We successfully apply a GA to find optimal resource allocations in four different representative network scenarios and validate features of our model experimentally in a quantum local area network in deployed fiber. Our results show promise for the rapid design of large-scale entanglement distribution networks.


title: Fix Bugs with Transformer through a Neural-Symbolic Edit Grammar

abstract: We introduce NSEdit (neural-symbolic edit), a novel Transformer-based code repair method. Given only the source code that contains bugs, NSEdit predicts an editing sequence that can fix the bugs. The edit grammar is formulated as a regular language, and the Transformer uses it as a neural-symbolic scripting interface to generate editing programs. We modify the Transformer and add a pointer network to select the edit locations. An ensemble of rerankers are trained to re-rank the editing sequences generated by beam search. We fine-tune the rerankers on the validation set to reduce over-fitting. NSEdit is evaluated on various code repair datasets and achieved a new state-of-the-art accuracy ($24.04\%$) on the Tufano small dataset of the CodeXGLUE benchmark. NSEdit performs robustly when programs vary from packages to packages and when buggy programs are concrete. We conduct detailed analysis on our methods and demonstrate the effectiveness of each component.


title: METRO: Efficient Denoising Pretraining of Large Scale Autoencoding Language Models with Model Generated Signals

abstract: We present an efficient method of pretraining large-scale autoencoding language models using training signals generated by an auxiliary model. Originated in ELECTRA, this training strategy has demonstrated sample-efficiency to pretrain models at the scale of hundreds of millions of parameters. In this work, we conduct a comprehensive empirical study, and propose a recipe, namely "Model generated dEnoising TRaining Objective" (METRO), which incorporates some of the best modeling techniques developed recently to speed up, stabilize, and enhance pretrained language models without compromising model effectiveness. The resultant models, METRO-LM, consisting of up to 5.4 billion parameters, achieve new state-of-the-art on the GLUE, SuperGLUE, and SQuAD benchmarks. More importantly, METRO-LM are efficient in that they often outperform previous large models with significantly smaller model sizes and lower pretraining cost.


title: Wassmap: Wasserstein Isometric Mapping for Image Manifold Learning

abstract: In this paper, we propose Wasserstein Isometric Mapping (Wassmap), a nonlinear dimensionality reduction technique that provides solutions to some drawbacks in existing global nonlinear dimensionality reduction algorithms in imaging applications. Wassmap represents images via probability measures in Wasserstein space, then uses pairwise Wasserstein distances between the associated measures to produce a low-dimensional, approximately isometric embedding. We show that the algorithm is able to exactly recover parameters of some image manifolds including those generated by translations or dilations of a fixed generating measure. Additionally, we show that a discrete version of the algorithm retrieves parameters from manifolds generated from discrete measures by providing a theoretical bridge to transfer recovery results from functional data to discrete data. Testing of the proposed algorithms on various image data manifolds show that Wassmap yields good embeddings compared with other global and local techniques.


title: The Many-Faceted Light Curves of Young Disk-bearing Stars in Taurus as Seen by K2

abstract: We present a comprehensive study of the variability properties of young disk-bearing stars in the Taurus star-forming region, paralleling our previous investigation in rho Oph and Upper Sco. A sample of 99 confirmed Taurus association members is placed in the diagnostic Q-M plane of flux asymmetry (M) and quasi-periodicity (Q), which guides our assignment of variability classes. We find a similar proportion of flux-symmetric variables in Taurus, but more bursters and fewer dippers relative to Upper Sco. The regions also differ in that the amplitudes for periodic and quasi-periodic sources are larger in Taurus relative to the more evolved Upper Sco star/disk systems. The relationship between photometric variability patterns at optical wavelengths, which arise in the inner disk and at the stellar surface, are assessed relative to available disk inclination measurements.


title: Copiloting Autonomous Multi-Robot Missions: A Game-inspired Supervisory Control Interface

abstract: Real-world deployment of new technology and capabilities can be daunting. The recent DARPA Subterranean (SubT) Challenge, for instance, aimed at the advancement of robotic platforms and autonomy capabilities in three one-year development pushes. While multi-agent systems are traditionally deployed in controlled and structured environments that allow for controlled testing (e.g., warehouses), the SubT challenge targeted various types of unknown underground environments that imposed the risk of robot loss in the case of failure. In this work, we introduce a video game-inspired interface, an autonomous mission assistant, and test and deploy these using a heterogeneous multi-agent system in challenging environments. This work leads to improved human-supervisory control for a multi-agent system reducing overhead from application switching, task planning, execution, and verification while increasing available exploration time with this human-autonomy teaming platform.


title: Simplicial quantum contextuality

abstract: We introduce a new framework for contextuality based on simplicial sets, combinatorial models of topological spaces that play a prominent role in modern homotopy theory. Our approach extends measurement scenarios to consist of spaces (rather than sets) of measurements and outcomes, and thereby generalizes nonsignaling distributions to simplicial distributions, which are distributions on spaces modeled by simplicial sets. Using this formalism we present a topologically inspired new proof of Fine's theorem for characterizing noncontextuality in Bell scenarios. Strong contextuality is generalized suitably for simplicial distributions, allowing us to define cohomological witnesses that extend the earlier topological constructions restricted to algebraic relations among quantum observables to the level of probability distributions. Foundational theorems of quantum theory such as the Gleason's theorem and Kochen-Specker theorem can be expressed naturally within this new language.


title: Pointed Gromov-Hausdorff Topological Stability for non-compact metric spaces

abstract: We combine the pointed Gromov-Hausdorff metric [Ron10] with the locally $C^0$ distance to obtain the pointed $C^0$-Gromov-Hausdorff distance between maps of possibly different non-compact pointed metric spaces. The latter is then combined with Walters's locally topological stability [LNY18] and $GH$-stability from [AMR17] to obtain the notion of topologically $GH$-stable pointed homeomorphism. We give one example to show the difference between the distance when take different base point in a pointed metric space.


title: $\mathsf{SOCA}$ and $\mathsf{OGA}$ for $\mathsf{HL}$ spaces with strong properties

abstract: We study open colorings in certain classes of hereditary Lindel\"{o}f ($\mathsf{HL}$) spaces and submetrizable spaces. In particular, we show that the definible version for the Open Graph Axiom ($\mathsf{OGA}$) holds for the class of $\mathsf{HL}$ strong Choquet submetrizable spaces extending a well-known result of Feng. Furthermore, we show the consistency of the Open Graph Axiom for regular spaces that have countable spread and it's square also has it, reaching closer to a well known conjecture of Todor\v{c}evi\'{c}: "It is consistent that all regular spaces with countable spread satisfy $\mathsf{OGA}$".


title: Observability Analysis of Visual-Inertial Odometry with Online Calibration of Velocity-Control Based Kinematic Motion Models

abstract: In this paper, we analyze the observability of the visual-inertial odometry (VIO) using stereo cameras with a velocity-control based kinematic motion model. Previous work shows that in general case the global position and yaw are unobservable in VIO system, additionally the roll and pitch become also unobservable if there is no rotation. We prove that by integrating a planar motion constraint roll and pitch become observable. We also show that the parameters of the motion model are observable.


title: Joint Coreset Construction and Quantization for Distributed Machine Learning

abstract: Coresets are small, weighted summaries of larger datasets, aiming at providing provable error bounds for machine learning (ML) tasks while significantly reducing the communication and computation costs. To achieve a better trade-off between ML error bounds and costs, we propose the first framework to incorporate quantization techniques into the process of coreset construction. Specifically, we theoretically analyze the ML error bounds caused by a combination of coreset construction and quantization. Based on that, we formulate an optimization problem to minimize the ML error under a fixed budget of communication cost. To improve the scalability for large datasets, we identify two proxies of the original objective function, for which efficient algorithms are developed. For the case of data on multiple nodes, we further design a novel algorithm to allocate the communication budget to the nodes while minimizing the overall ML error. Through extensive experiments on multiple real-world datasets, we demonstrate the effectiveness and efficiency of our proposed algorithms for a variety of ML tasks. In particular, our algorithms have achieved more than 90% data reduction with less than 10% degradation in ML performance in most cases.


title: Sketching Algorithms and Lower Bounds for Ridge Regression

abstract: We give a sketching-based iterative algorithm that computes a $1+\varepsilon$ approximate solution for the ridge regression problem $\min_x \|Ax-b\|_2^2 +\lambda\|x\|_2^2$ where $A \in R^{n \times d}$ with $d \ge n$. Our algorithm, for a constant number of iterations (requiring a constant number of passes over the input), improves upon earlier work (Chowdhury et al.) by requiring that the sketching matrix only has a weaker Approximate Matrix Multiplication (AMM) guarantee that depends on $\varepsilon$, along with a constant subspace embedding guarantee. The earlier work instead requires that the sketching matrix has a subspace embedding guarantee that depends on $\varepsilon$. For example, to produce a $1+\varepsilon$ approximate solution in $1$ iteration, which requires $2$ passes over the input, our algorithm requires the OSNAP embedding to have $m= O(n\sigma^2/\lambda\varepsilon)$ rows with a sparsity parameter $s = O(\log(n))$, whereas the earlier algorithm of Chowdhury et al. with the same number of rows of OSNAP requires a sparsity $s = O(\sqrt{\sigma^2/\lambda\varepsilon} \cdot \log(n))$, where $\sigma = \opnorm{A}$ is the spectral norm of the matrix $A$. We also show that this algorithm can be used to give faster algorithms for kernel ridge regression. Finally, we show that the sketch size required for our algorithm is essentially optimal for a natural framework of algorithms for ridge regression by proving lower bounds on oblivious sketching matrices for AMM. The sketch size lower bounds for AMM may be of independent interest.


title: Exploring the Impacts of Power Grid Signals on Data Center Operations using a Receding-Horizon Scheduling Model

abstract: Data centers (DCs) can help decarbonize the power grid by helping absorb renewable power (e.g., wind and solar) due to their ability to shift power loads across space and time. However, to harness such load-shifting flexibility, it is necessary to understand how grid signals (carbon signals and market price/load allocations) affect DC operations. An obstacle that arises here is the lack of computationally-tractable DC operation models that can capture objectives, constraints, and information flows that arise at the interface of DCs and the power grid. To address this gap, we present a receding-horizon resource management model (a mixed-integer programming model) that captures the resource management layer between the DC scheduler and the grid while accounting for logical constraints, different types of objectives, and forecasts of incoming job profiles and of available computing capacity. We use our model to conduct extensive case studies based on public data from Microsoft Azure and MISO. Our studies show that DCs can provide significant temporal load-shifting flexibility that results in reduced carbon emissions and peak demand charges. Models and case studies are shared as easy-to-use Julia code.


title: Performance of a convolutional autoencoder designed to remove electronic noise from p-type point contact germanium detector signals

abstract: We present a convolutional autoencoder to denoise pulses from a p-type point contact high-purity germanium detector similar to those used in several rare event searches. While we focus on training procedures that rely on detailed detector physics simulations, we also present implementations requiring only noisy detector pulses to train the model. We validate our autoencoder on both simulated data and calibration data from an $^{241}$Am source, the latter of which is used to show that the denoised pulses are statistically compatible with data pulses. We demonstrate that our denoising method is able to preserve the underlying shapes of the pulses well, offering improvement over traditional denoising methods. We also show that the shaping time used to calculate energy with a trapezoidal filter can be significantly reduced while maintaining a comparable energy resolution. Under certain circumstances, our denoising method can improve the overall energy resolution. The methods we developed to remove electronic noise are straightforward to extend to other detector technologies. Furthermore, the latent representation from the encoder is also of use in quantifying shape-based characteristics of the signals. Our work has great potential to be used in particle physics experiments and beyond.


title: Photodynamical analysis of the nearly resonant planetary system WASP-148: Accurate transit-timing variations and mutual orbital inclination

abstract: WASP-148 is a recently announced extra-solar system harbouring at least two giant planets. The inner planet transits its host star. The planets travel on eccentric orbits and are near the 4:1 mean-motion resonance, which implies significant mutual gravitational interactions. In particular, this causes transit-timing variations of a few minutes, which were detected based on ground-based photometry. This made WASP-148 one of the few cases where such a phenomenon was detected without space-based photometry. Here, we present a self-consistent model of WASP-148 that takes into account the gravitational interactions between all known bodies in the system. Our analysis simultaneously fits the available radial velocities and transit light curves. In particular, we used the photometry secured by the TESS space telescope and made public after the WASP-148 discovery announcement. The TESS data confirm the transit-timing variations, but only in combination with previously measured transit times. The system parameters we derived agree with those previously reported and have a significantly improved precision, including the mass of the non-transiting planet. We found a significant mutual inclination between the orbital planes of the two planets: I=41.0 +6.2 -7.6 deg based on the modelling of the observations, although we found I=20.8 +/- 4.6 deg when we imposed a constraint on the model enforcing long-term dynamical stability. When a third planet was added to the model - based on a candidate signal in the radial velocity - the mutual inclination between planets b and c changed significantly allowing solutions closer to coplanar. We conclude that more data are needed to establish the true architecture of the system. If the significant mutual inclination is confirmed, WASP-148 would become one of the only few candidate non-coplanar planetary systems. We discuss possible origins for this misalignment.


title: A Bayesian Machine Learning Approach for Estimating Heterogeneous Survivor Causal Effects: Applications to a Critical Care Trial

abstract: Motivated by the Acute Respiratory Distress Syndrome Network (ARDSNetwork) ARDS respiratory management (ARMA) trial, we developed a flexible Bayesian machine learning approach to estimate the average causal effect and heterogeneous causal effects among the always-survivors stratum when clinical outcomes are subject to truncation. We adopted Bayesian additive regression trees (BART) to flexibly specify separate models for the potential outcomes and latent strata membership. In the analysis of the ARMA trial, we found that the low tidal volume treatment had an overall benefit for participants sustaining acute lung injuries on the outcome of time to returning home, but substantial heterogeneity in treatment effects among the always-survivors, driven most strongly by sex and the alveolar-arterial oxygen gradient at baseline (a physiologic measure of lung function and source of hypoxemia). These findings illustrate how the proposed methodology could guide the prognostic enrichment of future trials in the field. We also demonstrated through a simulation study that our proposed Bayesian machine learning approach outperforms other parametric methods in reducing the estimation bias in both the average causal effect and heterogeneous causal effects for always-survivors.


title: On Periodic Decompositions and Nonexpansive Lines

abstract: In his Ph.D. thesis, Michal Szabados conjectured that for a not fully periodic configuration with a minimal periodic decomposition the nonexpansive lines are exactly the lines that contain a period for some periodic configuration in such decomposition. In this paper, we study Szabados's conjecture. First, we show that we may consider a minimal periodic decomposition where each periodic configuration is defined on a finite alphabet. Then we prove that Szabados's conjecture holds for a wide class of configurations, which includes all not fully periodic low convex pattern complexity configurations.


title: The Transferability Limits of Static Benchmarks

abstract: Every practical method to solve the Schr\"odinger equation for interacting many-particle systems introduces approximations. Such methods are therefore plagued by systematic errors. For computational chemistry, it is decisive to quantify the specific error for some system under consideration. Traditionally, the primary resource for such an error assessment have been benchmarking results, usually taken from the literature. However, their transferability to a specific molecular system, and hence, the reliability of the traditional approach always remains uncertain to some degree. In this communication, we elaborate on the shortcomings of this traditional way of static benchmarking by exploiting statistical analyses at the example of one of the largest quantum chemical benchmark sets available. We demonstrate the uncertainty of error estimates in the light of the choice of reference data selected for a benchmark study. To alleviate the issues with static benchmarks, we advocate to rely instead on a rolling and system-focused approach for rigorously quantifying the uncertainty of a quantum chemical result.


title: Second Order Regret Bounds Against Generalized Expert Sequences under Partial Bandit Feedback

abstract: We study the problem of expert advice under partial bandit feedback setting and create a sequential minimax optimal algorithm. Our algorithm works with a more general partial monitoring setting, where, in contrast to the classical bandit feedback, the losses can be revealed in an adversarial manner. Our algorithm adopts a universal prediction perspective, whose performance is analyzed with regret against a general expert selection sequence. The regret we study is against a general competition class that covers many settings (such as the switching or contextual experts settings) and the expert selection sequences in the competition class are determined by the application at hand. Our regret bounds are second order bounds in terms of the sum of squared losses and the normalized regret of our algorithm is invariant under arbitrary affine transforms of the loss sequence. Our algorithm is truly online and does not use any preliminary information about the loss sequences.

categories: cs.LG cs.IT eess.SP math.IT math.OC stat.ML


title: ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

abstract: Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet.


title: Cosmic ray spectrum of protons plus helium nuclei between 6 TeV and 158 TeV from HAWC data

abstract: A measurement with high statistics of the differential energy spectrum of light elements in cosmic rays, in particular, of primary H plus He nuclei, is reported. The spectrum is presented in the energy range from $6$ to $158$ TeV per nucleus. Data was collected with the High Altitude Water Cherenkov (HAWC) Observatory between June 2015 and June 2019. The analysis was based on a Bayesian unfolding procedure, which was applied on a subsample of vertical HAWC data that was enriched to $82\%$ of events induced by light nuclei. To achieve the mass separation, a cut on the lateral age of air shower data was set guided by predictions of CORSIKA/QGSJET-II-04 simulations. The measured spectrum is consistent with a broken power-law spectrum and shows a kneelike feature at around $E = 24.0^{+3.6}_{-3.1} $ TeV, with a spectral index $\gamma = -2.51 \pm 0.02$ before the break and with $\gamma = -2.83 \pm 0.02$ above it. The feature has a statistical significance of $4.1 \, \sigma$. Within systematic uncertainties, the significance of the spectral break is $0.8 \, \sigma$.


title: Renormalized Area for Minimal Hypersurfaces of 5D Poincar\'e-Einstein Spaces

abstract: In this paper we derive a Gauss-Bonnet formula for the renormalized area of Graham-Witten minimal hypersurfaces of 5-dimensional Poincar\'e-Einstein spaces. The formula we derive expresses the renormalized area in terms of integrals of pointwise scalar Riemannian invariants. We also prove a result which gives a characterization of minimal hypersurfaces with $L^2$ second fundamental form in terms of conformal geometry at infinity.


title: Achieving Representative Data via Convex Hull Feasibility Sampling Algorithms

abstract: Sampling biases in training data are a major source of algorithmic biases in machine learning systems. Although there are many methods that attempt to mitigate such algorithmic biases during training, the most direct and obvious way is simply collecting more representative training data. In this paper, we consider the task of assembling a training dataset in which minority groups are adequately represented from a given set of data sources. In essence, this is an adaptive sampling problem to determine if a given point lies in the convex hull of the means from a set of unknown distributions. We present adaptive sampling methods to determine, with high confidence, whether it is possible to assemble a representative dataset from the given data sources. We also demonstrate the efficacy of our policies in simulations in the Bernoulli and a multinomial setting.


title: On a System of Weakly Null Semilinear Equations

abstract: We develop a new method for addressing certain weakly null systems of wave equations. This approach does not rely on Lorentz invariance nor on the use of null foliations, both of which restrict applications to, e.g., multiple speed systems. The proof uses a class of space-time Klainerman-Sobolev estimates of the first author, Tataru, and Tohaneanu, which pair nicely with local energy estimates that combine the $r^{p}$-weighted method of Dafermos and Rodnianski with the ghost weight method of Alinhac. We further refine the standard local energy estimate with a modification of the $\partial_{t} - \partial_{r}$ portion of the multiplier.


title: Explicit caching HYB: a new high-performance SpMV framework on GPGPU

abstract: Sparse Matrix-Vector Multiplication (SpMV) is a critical operation for the iterative solver of Finite Element Methods on computer simulation. Since the SpMV operation is a memory-bound algorithm, the efficiency of data movements heavily influenced the performance of the SpMV on GPU. In recent years, many research is conducted in accelerating the performance of SpMV on the graphic processing units (GPU). The performance optimization methods used in existing studies focus on the following areas: improve the load balancing between GPU processors, and reduce the execution divergence between GPU threads. Although some studies have made preliminary optimization on the input vector fetching, the effect of explicitly caching the input vector on GPU base SpMV has not been studied in depth yet. In this study, we are trying to minimize the data movements cost for GPU-based SpMV using a new framework named "explicit caching Hybrid (EHYB)". The EHYB framework achieved significant performance improvement by using the following methods: 1. Improve the speed of data movements by partitioning and explicitly caching the input vector to the shared memory of the CUDA kernel. 2. Reduce the volume of data movements by storing the major part of the column index with a compact format. We tested our implementation with sparse matrices derived from FEM applications in different areas. The experiment results show that our implementation can overperform the state-of-the-arts implementation with significant speedup, and leads to higher FLOPs than the theoryperformance up-boundary of the existing GPU-based SpMV implementations.


title: Review of CMS contribution to Hadron Spectroscopy (Snowmass 2021 White Paper)

abstract: In this brief White Paper we glance over the studies carried out so far by the CMS Collaboration in the field of conventional and exotic hadron spectroscopy. We highlight the most relevant scientific achievements and discuss the future perspectives of this engagement.


title: High-precision half-life determination of $^{14}$O via direct $\beta$ counting

abstract: The half-life of the superallowed Fermi $\beta^+$ emitter $^{14}$O was determined to high precision via a direct $\beta$ counting experiment performed at the Isotope Separator and Accelerator (ISAC) facility at TRIUMF. The result, $T_{1/2}$($^{14}$O) = 70619.2(76) ms, is consistent with, but is more precise than, the world average obtained from 11 previous measurements. Combining the $^{14}$O half-life deduced in the present work with the previous most precise measurements of this quantity leads to a reduction in the overall uncertainty, by nearly a factor of 2. The new world average is $T_{1/2}$($^{14}$O) = 70619.6(63) ms with a reduced $\chi^2$ value of 0.87 obtained from 8 degrees of freedom.


title: Percolation in a triangle on a square lattice

abstract: Percolation on a plane is usually associated with clusters spanning two opposite sides of a rectangular system. Here we investigate three-leg clusters generated on a square lattice and spanning the three sides of equilateral triangles. If the position and orientation of the triangles relative to the lattice are uniformly randomized, one obtains an efficient method of determining the percolation threshold, on par with the most advanced Monte Carlo methods developed for the rectangular geometry. The universal crossing probability for three-leg clusters is geometry-independent, which opens a way for further improvements of the method.


title: SkiQL: A Unified Schema Query Language

abstract: Most NoSQL systems are schema-on-read: data can be stored without first having to declare a Schema that imposes a structure. This schemaless feature offers flexibility to evolve data-intensive applications when data frequently change. However, freeing from declaring schemas does not mean their absence, but rather that they are implicit in data and code. Therefore, diagramming tools similar to those available for relational systems are also needed to help developers and administrators understanding NoSQL schemas. Visualizing diagrams is not practical if schemas contain hundreds of database entities, and exploration or query facilities are then needed. In schemaless NoSQL stores, data of the same entity can be stored with different structure which can increase the difficulty of having readable diagrams. NoSQL schema management tools should therefore have three main components: schema extraction, schema visualization, and schema query. Since that there exist four main NoSQL data models, it is convenient that such tools can be built on a generic data model that provide platform-independence to query and visualize schemas. With the aim of favoring the creation of generic database tools, the authors of this paper defined the U-Schema unified data model that integrates the four main NoSQL data models and the relational model. This paper is focused on querying NoSQL and relational schemas which are represented as U-Schema models. We present the SkiQL language designed on U-Schema to achieve a platform-independent schema query service. SkiQL provides two constructs: schema-query and relationship-query. The former allows to obtain information of entity or relationship types, and the latter that of the aggregations or references (relations among types). We will show how SkiQL was evaluated by calculating well-known metrics for languages and using a survey.


title: A deep learning algorithm for reducing false positives in screening mammography

abstract: Screening mammography improves breast cancer outcomes by enabling early detection and treatment. However, false positive callbacks for additional imaging from screening exams cause unnecessary procedures, patient anxiety, and financial burden. This work demonstrates an AI algorithm that reduces false positives by identifying mammograms not suspicious for breast cancer. We trained the algorithm to determine the absence of cancer using 123,248 2D digital mammograms (6,161 cancers) and performed a retrospective study on 14,831 screening exams (1,026 cancers) from 15 US and 3 UK sites. Retrospective evaluation of the algorithm on the largest of the US sites (11,592 mammograms, 101 cancers) a) left the cancer detection rate unaffected (p=0.02, non-inferiority margin 0.25 cancers per 1000 exams), b) reduced callbacks for diagnostic exams by 31.1% compared to standard clinical readings, c) reduced benign needle biopsies by 7.4%, and d) reduced screening exams requiring radiologist interpretation by 41.6% in the simulated clinical workflow. This work lays the foundation for semi-autonomous breast cancer screening systems that could benefit patients and healthcare systems by reducing false positives, unnecessary procedures, patient anxiety, and expenses.


title: Effective Field Theory of Gravitating Continuum: Solids, Fluids, and Aether Unified

abstract: We investigate the relativistic effective field theory (EFT) describing a non-dissipative gravitating continuum. In addition to ordinary continua, namely solids and fluids, we find an extraordinary more symmetric continuum, aether. In particular, the symmetry of the aether concludes that a homogeneous and isotropic state behaves like a cosmological constant. We formulate the EFT in the unitary/comoving gauge in which the dynamical degrees of freedom of the continuum (phonons) are eaten by the spacetime metric. This gauge choice, which is interpreted as the Lagrangian description in hydrodynamics, offers a neat geometrical understanding of continua. We examine a thread-based spacetime decomposition with respect to the four-velocity of the continuum which is different from the foliation-based Arnowitt-Deser-Misner one. Our thread-based decomposition respects the symmetries of the continua and, therefore, makes it possible to systematically find invariant building blocks of the EFT for each continuum even at higher orders in the derivative expansion. We also discuss the linear dynamics of the system and show that both gravitons and phonons acquire "masses" in a gravitating background.


title: Adaptive Modulation for Wobbling UAV Air-to-Ground Links in Millimeter-wave Bands

abstract: The emerging millimeter-wave (mm-wave) unmanned aerial vehicle (UAV) air-to-ground (A2G) communications are facing the Doppler effect problem that arises from the inevitable wobbling of the UAV. The fast time-varying channel for UAV A2G communications may lead to the outdated channel state information (CSI) from the channel estimation. In this paper, we introduce two detectors to demodulate the received signal and get the instantaneous bit error probability (BEP) of a mm-wave UAV A2G link under imperfect CSI. Based on the designed detectors, we propose an adaptive modulation scheme to maximize the average transmission rate under imperfect CSI by optimizing the data transmission time subject to the maximum tolerable BEP. A power control policy is in conjunction with adaptive modulation to minimize the transmission power while maintaining both the BEP under the threshold and the maximized average transmission rate. Numerical results show that the proposed adaptive modulation scheme in conjunction with the power control policy could maximize the temporally averaged transmission rate, while saves as much as 50\% energy.


title: GAP: A Graph-aware Language Model Framework for Knowledge Graph-to-Text Generation

abstract: Recent improvements in KG-to-text generation are due to additional auxiliary pre-training tasks designed to give the fine-tune task a boost in performance. These tasks require extensive computational resources while only suggesting marginal improvements. Here, we demonstrate that by fusing graph-aware elements into existing pre-trained language models, we are able to outperform state-of-the-art models and close the gap imposed by additional pre-training tasks. We do so by proposing a mask structure to capture neighborhood information and a novel type encoder that adds a bias to the graph-attention weights depending on the connection type. Experiments on two KG-to-text benchmark datasets show our models are competitive while involving fewer parameters and no additional pre-training tasks. By formulating the problem as a framework, we can interchange the various proposed components and begin interpreting KG-to-text generative models based on the topological and type information found in a graph.


title: Geometric Understanding of Sketches

abstract: Sketching is used as a ubiquitous tool of expression by novices and experts alike. In this thesis I explore two methods that help a system provide a geometric machine-understanding of sketches, and in-turn help a user accomplish a downstream task. The first work deals with interpretation of a 2D-line drawing as a graph structure, and also illustrates its effectiveness through its physical reconstruction by a robot. We setup a two-step pipeline to solve the problem. Formerly, we estimate the vertices of the graph with sub-pixel level accuracy. We achieve this using a combination of deep convolutional neural networks learned under a supervised setting for pixel-level estimation followed by the connected component analysis for clustering. Later we follow it up with a feedback-loop-based edge estimation method. To complement the graph-interpretation, we further perform data-interchange to a robot legible ASCII format, and thus teach a robot to replicate a line drawing. In the second work, we test the 3D-geometric understanding of a sketch-based system without explicit access to the information about 3D-geometry. The objective is to complete a contour-like sketch of a 3D-object, with illumination and texture information. We propose a data-driven approach to learn a conditional distribution modelled as deep convolutional neural networks to be trained under an adversarial setting; and we validate it against a human-in-the-loop. The method itself is further supported by synthetic data generation using constructive solid geometry following a standard graphics pipeline. In order to validate the efficacy of our method, we design a user-interface plugged into a popular sketch-based workflow, and setup a simple task-based exercise, for an artist. Thereafter, we also discover that form-exploration is an additional utility of our application.


title: DRAGON (Differentiable Graph Execution) : A suite of Hardware Simulation and Optimization tools for Modern AI/Non-AI Workloads

abstract: We introduce DRAGON, an open-source, fast and explainable hardware simulation and optimization toolchain that enables hardware architects to simulate hardware designs, and to optimize hardware designs to efficiently execute workloads. The DRAGON toolchain provides the following tools: Hardware Model Generator (DGen), Hardware Simulator (DSim) and Hardware Optimizer (DOpt). DSim provides the simulation of running algorithms (represented as data-flow graphs) on hardware described. DGen describes the hardware in detail, with user input architectures/technology (represented in a custom description language). A novel methodology of gradient descent from the simulation allows us optimize the hardware model (giving the directions for improvements in technology parameters and design parameters), provided by Dopt. DRAGON framework (DSim) is much faster than previously avaible works for simulation, which is possible through performance-first code writing practices, mathematical formulas for common computing operations to avoid cycle-accurate simulation steps, efficient algorithms for mapping, and data-structure representations for hardware state. DRAGON framework (Dopt) generates performance optimized architectures for both AI and Non-AI Workloads, and provides technology improvement directions for 100x-1000x better future computing systems.


title: Dynamic Schema Graph Fusion Network for Multi-Domain Dialogue State Tracking

abstract: Dialogue State Tracking (DST) aims to keep track of users' intentions during the course of a conversation. In DST, modelling the relations among domains and slots is still an under-studied problem. Existing approaches that have considered such relations generally fall short in: (1) fusing prior slot-domain membership relations and dialogue-aware dynamic slot relations explicitly, and (2) generalizing to unseen domains. To address these issues, we propose a novel \textbf{D}ynamic \textbf{S}chema \textbf{G}raph \textbf{F}usion \textbf{Net}work (\textbf{DSGFNet}), which generates a dynamic schema graph to explicitly fuse the prior slot-domain membership relations and dialogue-aware dynamic slot relations. It also uses the schemata to facilitate knowledge transfer to new domains. DSGFNet consists of a dialogue utterance encoder, a schema graph encoder, a dialogue-aware schema graph evolving network, and a schema graph enhanced dialogue state decoder. Empirical results on benchmark datasets (i.e., SGD, MultiWOZ2.1, and MultiWOZ2.2), show that DSGFNet outperforms existing methods.


title: Rotational Solitons for the Curve Shortening Flow on Revolution Surfaces

abstract: We present a characterization for the rotational soliton for the curve shortening flow (CSF) on the revolution surfaces of $\mathbb{R}^3$. Furthermore, we describe the behavior of such curves by showing that the two ends of each open curve are asymptotic to a parallel geodesic.


title: Weighted homological regularities

abstract: Let $A$ be a noetherian connected graded algebra. We introduce and study homological invariants that are weighted sums of the homological and internal degrees of cochain complexes of graded $A$-modules, providing weighted versions of Castelnuovo--Mumford regularity, Tor-regularity, Artin--Schelter regularity, and concavity. In some cases an invariant (such as Tor-regularity) that is infinite can be replaced with a weighted invariant that is finite, and several homological invariants of complexes can be expressed as weighted homological regularities. We prove a few weighted homological identities some of which unify different classical homological identities and produce interesting new ones.


title: Double version of the Rashba and Dresselhaus spin-orbit coupling

abstract: The Rashba and Dresselhaus types of spin-orbit coupling are two typical linear coupling forms. We establish the fundamental physics of a model which can be viewed as the double version of the Rashba and Dresselhaus spin-orbit coupling. This model describes the low energy physics of a class of massless Dirac fermions in spin-orbit systems. The physical properties of the massless Dirac fermions are determined by the mathematical relations of spin-orbit coefficients. For equal Rashba and Dresselhauss coupling constants, k-independent eigenspinors and a persistent spin helix combined with massless birefringent Dirac fermions emerge in this model. The spin-orbit coupled systems described by this model have potential technological applications from spintronics to quantum computation.


title: MINSU (Mobile Inventory And Scanning Unit):Computer Vision and AI

abstract: The MINSU(Mobile Inventory and Scanning Unit) algorithm uses the computational vision analysis method to record the residual quantity/fullness of the cabinet. To do so, it goes through a five-step method: object detection, foreground subtraction, K-means clustering, percentage estimation, and counting. The input image goes through the object detection method to analyze the specific position of the cabinets in terms of coordinates. After doing so, it goes through the foreground subtraction method to make the image more focus-able to the cabinet itself by removing the background (some manual work may have to be done such as selecting the parts that were not grab cut by the algorithm). In the K-means clustering method, the multi-colored image turns into a 3 colored monotonous image for quicker and more accurate analysis. At last, the image goes through percentage estimation and counting. In these two methods, the proportion that the material inside the cabinet is found in percentage which then is used to approximate the number of materials inside. Had this project been successful, the residual quantity management could solve the problem addressed earlier in the introduction.


title: GM-TOuNN: Graded Multiscale Topology Optimization using Neural Networks

abstract: Multiscale topology optimization (M-TO) entails generating an optimal global topology, and an optimal set of microstructures at a smaller scale, for a physics-constrained problem. With the advent of additive manufacturing, M-TO has gained significant prominence. However, generating optimal microstructures at various locations can be computationally very expensive. As an alternate, graded multiscale topology optimization (GM-TO) has been proposed where one or more pre-selected and graded (parameterized) microstructural topologies are used to fill the domain optimally. This leads to a significant reduction in computation while retaining many of the benefits of M-TO. A successful GM-TO framework must: (1) be capable of efficiently handling numerous pre-selected microstructures, (2) be able to continuously switch between these microstructures during optimization, (3) ensure that the partition of unity is satisfied, and (4) discourage microstructure mixing at termination. In this paper, we propose to meet these requirements by exploiting the unique classification capacity of neural networks. Specifically, we propose a graded multiscale topology optimization using neural-network (GM-TOuNN) framework with the following features: (1) the number of design variables is only weakly dependent on the number of pre-selected microstructures, (2) it guarantees partition of unity while discouraging microstructure mixing, and (3) it supports automatic differentiation, thereby eliminating manual sensitivity analysis. The proposed framework is illustrated through several examples.


title: Revisiting Transformer-based Models for Long Document Classification

abstract: The recent literature in text classification is biased towards short text sequences (e.g., sentences or paragraphs). In real-world applications, multi-page multi-paragraph documents are common and they cannot be efficiently encoded by vanilla Transformer-based models. We compare different Transformer-based Long Document Classification (TrLDC) approaches that aim to mitigate the computational overhead of vanilla transformers to encode much longer text, namely sparse attention and hierarchical encoding methods. We examine several aspects of sparse attention (e.g., size of local attention window, use of global attention) and hierarchical (e.g., document splitting strategy) transformers on four document classification datasets covering different domains. We observe a clear benefit from being able to process longer text, and, based on our results, we derive practical advice of applying Transformer-based models on long document classification tasks.


title: Multifidelity deep neural operators for efficient learning of partial differential equations with application to fast inverse design of nanoscale heat transport

abstract: Deep neural operators can learn operators mapping between infinite-dimensional function spaces via deep neural networks and have become an emerging paradigm of scientific machine learning. However, training neural operators usually requires a large amount of high-fidelity data, which is often difficult to obtain in real engineering problems. Here, we address this challenge by using multifidelity learning, i.e., learning from multifidelity datasets. We develop a multifidelity neural operator based on a deep operator network (DeepONet). A multifidelity DeepONet includes two standard DeepONets coupled by residual learning and input augmentation. Multifidelity DeepONet significantly reduces the required amount of high-fidelity data and achieves one order of magnitude smaller error when using the same amount of high-fidelity data. We apply a multifidelity DeepONet to learn the phonon Boltzmann transport equation (BTE), a framework to compute nanoscale heat transport. By combining a trained multifidelity DeepONet with genetic algorithm or topology optimization, we demonstrate a fast solver for the inverse design of BTE problems.


title: Hard X-ray helical dichroism of disordered molecular media

abstract: Chirality is a structural property of molecules lacking mirror symmetry that has strong implications in diverse fields, ranging from life to materials sciences. Established spectroscopic methods that are sensitive to chirality, such as circular dichroism (CD), exhibit weak signal contributions on an achiral background. Helical dichroism (HD), which is based on the orbital angular momentum (OAM) of light, offers a new approach to probe molecular chirality, but it has never been demonstrated on disordered samples. Furthermore, in the optical domain the challenge lies in the need to transfer the OAM of the photon to an electron that is localized on an {\AA}-size orbital. Here, we overcome this challenge using hard X-rays with spiral Fresnel zone, which can induce an OAM. We present the first HD spectra of a disordered powder sample of enantiopure molecular complexes of [Fe(4,4'-diMebpy)3]2+ at the iron K-edge (7.1 keV) with OAM-carrying beams. The HD spectra exhibit the expected inversions of signs switching from a left to a right helical wave front or from an enantiomer to the other. The asymmetry ratios for the HD spectra are within one to five percent for OAM beams with topological charges of one and three. These results open a new window into the studies of molecular chirality and its interaction with the orbital angular momentum of light.


title: Isoperimetric Inequalities Made Simpler

abstract: We give an alternative, simple method to prove isoperimetric inequalities over the hypercube. In particular, we show: 1. An elementary proof of classical isoperimetric inequalities of Talagrand, as well as a stronger isoperimetric result conjectured by Talagrand and recently proved by Eldan and Gross. 2. A strengthening of the Friedgut junta theorem, asserting that if the $p$-moment of the sensitivity of a function is constant for some $p\geq 1/2 + \varepsilon$, then the function is close to a junta. In this language, Friedgut's theorem is the special case that $p=1$.


title: Designing Experiments Toward Shrinkage Estimation

abstract: We consider how increasingly available observational data can be used to improve the design of randomized controlled trials (RCTs). We seek to design a prospective RCT, with the intent of using an Empirical Bayes estimator to shrink the causal estimates from our trial toward causal estimates obtained from an observational study. We ask: how might we design the experiment to better complement the observational study in this setting? We propose using an estimator that shrinks each component of the RCT causal estimator toward its observational counterpart by a factor proportional to its variance. First, we show that the risk of this estimator can be computed efficiently via numerical integration. We then propose algorithms for determining the best allocation of units to strata (the best "design"). We consider three options: Neyman allocation; a "naive" design assuming no unmeasured confounding in the observational study; and a "defensive" design accounting for the imperfect parameter estimates we would obtain from the observational study with unmeasured confounding. We also incorporate results from sensitivity analysis to establish guardrails on the designs, so that our experiment could be reasonably analyzed with and without shrinkage. We demonstrate the superiority of these experimental designs with a simulation study involving causal inference on a rare, binary outcome.


title: Time Series of Non-Additive Metrics: Identification and Interpretation of Contributing Factors of Variance by Linear Decomposition

abstract: The research paper addresses linear decomposition of time series of non-additive metrics that allows for the identification and interpretation of contributing factors (input features) of variance. Non-additive metrics, such as ratios, are widely used in a variety of domains. It commonly requires preceding aggregations of underlying variables that are used to calculate the metric of interest. The latest poses a dimensionality challenge when the input features and underlying variables are formed as two-dimensional arrays along elements, such as account or customer identifications, and time points. It rules out direct modeling of the time series of a non-additive metric as a function of input features. The article discusses a five-step approach: (1) segmentations of input features and the underlying variables of the metric that are supported by unsupervised autoencoders, (2) univariate or joint fittings of the metric by the aggregated input features on the segmented domains, (3) transformations of pre-screened input features according to the fitted models, (4) aggregation of the transformed features as time series, and (5) modelling of the metric time series as a sum of constrained linear effects of the aggregated features. Alternatively, approximation by numerical differentiation has been considered to linearize the metric. It allows for element level univariate or joint modeling of step (2). The process of these analytical steps allows for a backward-looking explanatory decomposition of the metric as a sum of time series of the survived input features. The paper includes a synthetic example that studies loss-to-balance monthly rates of a hypothetical retail credit portfolio. To validate that no latent factors other than the survived input features have significant impacts on the metric, Statistical Process Control has been introduced for the residual time series.


title: Experimental quantum state discrimination using the optimal fixed rate of inconclusive outcomes strategy

abstract: The problem of non-orthogonal state discrimination underlies crucial quantum information tasks, such as cryptography and computing protocols. Therefore, it is decisive to find optimal scenarios for discrimination among quantum states. We experimentally investigate the strategy for the optimal discrimination of two non-orthogonal states considering a fixed rate of inconclusive outcomes (FRIO). The main advantage of the FRIO strategy is to interpolate between unambiguous and minimum error discrimination by solely adjusting the rate of inconclusive outcomes. We present a versatile experimental scheme that performs the optimal FRIO measurement for any pair of generated non-orthogonal states with arbitrary a priori probabilities and for any fixed rate of inconclusive outcomes. Considering different values of the free parameters in the FRIO protocol, we implement it upon qubit states encoded in the polarization mode of single photons generated in the spontaneous parametric down-conversion process. Moreover, we resort to a newfangled double-path Sagnac interferometer to perform a three-outcome non-projective measurement required for the discrimination task, showing excellent agreement with the theoretical prediction. This experiment provides a practical toolbox for a wide range of quantum state discrimination strategies using the FRIO scheme, which can greatly benefit quantum information applications and fundamental studies in quantum theory.


title: Creep tide theory. Equations for differentiated bodies with aligned layers

abstract: The creep tide theory is used to establish the basic equations of the tidal evolution of differentiated bodies formed by aligned homogeneous layers in co-rotation. The mass concentration of the body is given by the fluid Love number $k_f$. The formulas are given by series expansions valid for high eccentricity systems. They are equivalent to Darwin's equations, but formally more compact. An application to the case of Enceladus, with $k_f=0.942$, is discussed.


title: Gas and Dust Shadows in the TW Hydrae Disk

abstract: We present new observations of CO J=2-1 emission from the protoplanetary disk around TW Hya. Emission is detected out to 240 au (4") and found to exhibit azimuthal variations up to 20% beyond 180 au (3"), with the west side of the disk brighter than the east. This asymmetry is interpreted as tracing the shadow previously seen in scattered light. A reanalysis of the multi-epoch observations of the dust shadow in scattered light from Debes et al. (2017) suggests that an oscillatory motion would provide a better model of the temporal evolution of the dust shadow rather than orbital motion. Both models predict an angular offset between the dust shadow and the gas shadow of up to ~100 deg. We attribute this offset to the finite rate at which dust grains and gas molecules can exchange heat, dominated by the collisional rate between gas molecules and dust grains, $t_{\rm coll}$. The angular offsets derived are equivalent to collisional timescales that range from the near instantaneous up to $t_{\rm coll}$ ~ 10 years, depending on whether a straight or a curved dust shadow, as suggested by HST observations reported by Debes et al. (2017), is adopted. The inferred range of $t_{\rm coll}$ are consistent with those predictions based on representative gas densities, temperatures, gas-to-dust ratios and grain sizes. These results represent the first time empirical constraints can be placed on $t_{\rm coll}$.


title: Stability of China's Stock Market: Measure and Forecast by Ricci Curvature on Network

abstract: The systemic stability of a stock market is one of the core issues in the financial field. The market can be regarded as a complex network whose nodes are stocks connected by edges that signify their correlation strength. Since the market is a strongly nonlinear system, it is difficult to measure the macroscopic stability and depict market fluctuations in time. In this paper, we use a geometric measure derived from discrete Ricci curvature to capture the higher-order nonlinear architecture of financial networks. In order to confirm the effectiveness of our method, we use it to analyze the CSI 300 constituents of China's stock market from 2005--2020 and the systemic stability of the market is quantified through the network's Ricci type curvatures. Furthermore, we use a hybrid model to analyze the curvature time series and predict the future trends of the market accurately. As far as we know, this is the first paper to apply Ricci curvature to forecast the systemic stability of domestic stock market, and our results show that Ricci curvature has good explanatory power for the market stability and can be a good indicator to judge the future risk and volatility of the domestic market.


title: Learning fixed-complexity polyhedral Lyapunov functions from counterexamples

abstract: We study the problem of synthesizing polyhedral Lyapunov functions for hybrid linear systems. Such functions are defined as convex piecewise linear functions, with a finite number of pieces. We first prove that deciding whether there exists an $m$-piece polyhedral Lyapunov function for a given hybrid linear system is NP-hard. We then present a counterexample-guided algorithm for solving this problem. The algorithm alternates between choosing a candidate polyhedral function based on a finite set of counterexamples and verifying whether the candidate satisfies the Lyapunov conditions. If the verification fails, we find a new counterexample that is added to our set. We prove that if the algorithm terminates, it discovers a valid Lyapunov function or concludes that no such Lyapunov function exists. However, our initial algorithm can be non-terminating. We modify our algorithm to provide a terminating version based on the so-called cutting-plane argument from nonsmooth optimization. We demonstrate our algorithm on numerical examples.


title: Feasibility of the experimental study of $D_{s}^{\ast}$ ${\to}$ ${\phi}{\pi}$ decay

abstract: The current knowledge on the $D_{s}^{\ast}$ meson are very limited. Besides the dominant electromagnetic decays, the $D_{s}^{\ast}$ weak decays are legal and offer the valuable opportunities to explore the wanted $D_{s}^{\ast}$ meson. In this paper, the $D_{s}^{\ast}$ ${\to}$ ${\phi}{\pi}$ decay was studied with the factorization approach. It is found that the branching ratio ${\cal B}(D_{s}^{\ast}{\to}{\phi}{\pi})$ ${\sim}$ ${\cal O}(10^{-7})$, which corresponds to several thousands of events at the $e^{+}e^{-}$ collider experiments including STCF, SuperKEKB, CEPC and FCC-ee, and several millions of events at the hadron collider experiments, such as LHCb@HL-LHC. It is feasible to experimentally study the $D_{s}^{\ast}$ ${\to}$ ${\phi}{\pi}$ weak decay in the future, even considering the identification efficiency.


title: Breakdown voltage and linear temperature drift in a single-molecule junction

abstract: Using first-principles calculations based on density functional theory combined with the non-equilibrium Green's function approach, the transport behaviors of a single-molecule junction formed by benzenedithiol connected to gold electrodes are investigated. The breakdown voltage for the model of benzenedithiol plus gold electrodes is 0.7 V, which is close to the experimental value. A linear response between the conductance and temperature (known as linear temperature drift) is found in the molecular device, which indicates that it could be used to maintain the stability of molecular circuits. Meanwhile, input and output with the same accuracies would be useful for designing multi-level circuits, which would be used to improve the resolution ratio in analog-to-digital converters. The present findings indicate that benzenedithiol-based single-molecule junctions would be promising functional units for molecular sensors.


title: Dendrite formation in rechargeable lithium-metal batteries: Phase-field modeling using open-source finite element library

abstract: We describe a phase-field model for the electrodeposition process that forms dendrites within metal-anode batteries. We derive the free energy functional model, arriving at a system of partial differential equations that describe the evolution of a phase field, the lithium-ion concentration, and an electric potential. We formulate, discretize, and solve the set of partial differential equations describing the coupled electrochemical interactions during a battery charge cycle using an open-source finite element library. The open-source library allows us to use parallel solvers and time-marching adaptivity. We describe two- and three-dimensional simulations; these simulations agree with experimentally-observed dendrite growth rates and morphologies reported in the literature.