title: Bayesian non-parametric ordinal regression under a monotonicity constraint

abstract: Compared to the nominal scale, the ordinal scale for a categorical outcome variable has the property of making a monotonicity assumption for the covariate effects meaningful. This assumption is encoded in the commonly used proportional odds model, but there it is combined with other parametric assumptions such as linearity and additivity. Herein, the considered models are non-parametric and the only condition imposed is that the effects of the covariates on the outcome categories are stochastically monotone according to the ordinal scale. We are not aware of the existence of other comparable multivariable models that would be suitable for inference purposes. We generalize our previously proposed Bayesian monotonic multivariable regression model to ordinal outcomes, and propose an estimation procedure based on reversible jump Markov chain Monte Carlo. The model is based on a marked point process construction, which allows it to approximate arbitrary monotonic regression function shapes, and has a built-in covariate selection property. We study the performance of the proposed approach through extensive simulation studies, and demonstrate its practical application in two real data examples.


title: Secure Beamforming and Ergodic Secrecy Rate Analysis for Amplify-and-Forward Relay Networks with Wireless Powered Jammer

abstract: In this correspondence, we consider an amplify-and-forward relay network in which relayed information is overheard by an eavesdropper. In order to confound the eavesdropper, a wireless-powered jammer is also considered which harvests energy from a multiple-antenna source. We proposed a new secure beamforming scheme in which beamforming vector is a linear combination of the energy beamforming (EB) and information beamforming (IB) vectors. We also present a new closed-form solution for the proposed beamforming vector which is shown to achieve a higher secrecy rate as compared to the trivial EB and IB vectors. Moreover, a tight closed-form approximation for the ergodic secrecy rate is also derived for the asymptotic regime of a large number of antennas at the source. Finally, numerical examples and simulations are provided which validate our analytical results.


title: Tubes of finite $II$-type Gauss map

abstract: In this paper, we continue the classification of finite type Gauss map surfaces in the 3-dimensional Euclidean space $\mathbb{E}^{3}$. We present an important family of surfaces, namely, tubes in $\mathbb{E}^{3}$. We show that the Gauss map of a tube is of an infinite type corresponding to the second fundamental form.


title: Observation of tunable mobility edges in generalized Aubry-Andr\'{e} lattices

abstract: Using synthetic lattices of laser-coupled atomic momentum modes, we experimentally realize a recently proposed family of nearest-neighbor tight-binding models having quasiperiodic site energy modulation that host an exact mobility edge protected by a duality symmetry. These one-dimensional tight-binding models can be viewed as a generalization of the well-known Aubry-Andr\'{e} (AA) model, with an energy-dependent self duality condition that constitutes an analytical mobility edge relation. By adiabatically preparing the lowest and highest energy eigenstates of this model system and performing microscopic measurements of their participation ratio, we track the evolution of the mobility edge as the energy-dependent density of states is modified by the model's tuning parameter. Our results show strong deviations from single-particle predictions, consistent with attractive interactions causing both enhanced localization of the lowest energy state due to self-trapping and inhibited localization of the highest energy state due to screening. This study paves the way for quantitative studies of interaction effects on self duality induced mobility edges.


title: Robust Linear Regression: Optimal Rates in Polynomial Time

abstract: We obtain robust and computationally efficient estimators for learning several linear models that achieve statistically optimal convergence rate under minimal distributional assumptions. Concretely, we assume our data is drawn from a $k$-hypercontractive distribution and an $\epsilon$-fraction is adversarially corrupted. We then describe an estimator that converges to the optimal least-squares minimizer for the true distribution at a rate proportional to $\epsilon^{2-2/k}$, when the noise is independent of the covariates. We note that no such estimator was known prior to our work, even with access to unbounded computation. The rate we achieve is information-theoretically optimal and thus we resolve the main open question in Klivans, Kothari and Meka [COLT'18]. Our key insight is to identify an analytic condition that serves as a polynomial relaxation of independence of random variables. In particular, we show that when the moments of the noise and covariates are negatively-correlated, we obtain the same rate as independent noise. Further, when the condition is not satisfied, we obtain a rate proportional to $\epsilon^{2-4/k}$, and again match the information-theoretic lower bound. Our central technical contribution is to algorithmically exploit independence of random variables in the "sum-of-squares" framework by formulating it as the aforementioned polynomial inequality.


title: Scalable Comparative Visualization of Ensembles of Call Graphs

abstract: Optimizing the performance of large-scale parallel codes is critical for efficient utilization of computing resources. Code developers often explore various execution parameters, such as hardware configurations, system software choices, and application parameters, and are interested in detecting and understanding bottlenecks in different executions. They often collect hierarchical performance profiles represented as call graphs, which combine performance metrics with their execution contexts. The crucial task of exploring multiple call graphs together is tedious and challenging because of the many structural differences in the execution contexts and significant variability in the collected performance metrics (e.g., execution runtime). In this paper, we present an enhanced version of CallFlow to support the exploration of ensembles of call graphs using new types of visualizations, analysis, graph operations, and features. We introduce ensemble-Sankey, a new visual design that combines the strengths of resource-flow (Sankey) and box-plot visualization techniques. Whereas the resource-flow visualization can easily and intuitively describe the graphical nature of the call graph, the box plots overlaid on the nodes of Sankey convey the performance variability within the ensemble. Our interactive visual interface provides linked views to help explore ensembles of call graphs, e.g., by facilitating the analysis of structural differences, and identifying similar or distinct call graphs. We demonstrate the effectiveness and usefulness of our design through case studies on large-scale parallel codes.


title: Conduction-radiation coupling between two closely-separated solids

abstract: In the theory of radiative heat exchanges between two closely-spaced bodies introduced by Polder and van Hove, no interplay between the heat carriers inside the materials and the photons crossing the separation gap is assumed. Here we release this constraint by developing a general theory to describe the conduction-radiation coupling between two solids of arbitrary size separated by a subwavelength separation gap. We show that, as a result of the temperature profile induced by the coupling with conduction, the radiative heat flux exchanged between two parallel slabs at nanometric distances can be several orders of magnitude smaller than the one predicted by the conventional theory. These results could have important implications in the fields of nanoscale thermal management, near-field solid-state cooling and nanoscale energy conversion.


title: Adaptive Braking for Mitigating Gradient Delay

abstract: Neural network training is commonly accelerated by using multiple synchronized workers to compute gradient updates in parallel. Asynchronous methods remove synchronization overheads and improve hardware utilization at the cost of introducing gradient delay, which impedes optimization and can lead to lower final model performance. We introduce Adaptive Braking (AB), a modification for momentum-based optimizers that mitigates the effects of gradient delay. AB dynamically scales the gradient based on the alignment of the gradient and the velocity. This can dampen oscillations along high curvature directions of the loss surface, stabilizing and accelerating asynchronous training. We show that applying AB on top of SGD with momentum enables training ResNets on CIFAR-10 and ImageNet-1k with delays $D \geq$ 32 update steps with minimal drop in final test accuracy.


title: Estimation of pure quantum states in high dimension at the limit of quantum accuracy through complex optimization and statistical inference

abstract: Quantum tomography has become a key tool for the assessment of quantum states, processes, and devices. This drives the search for tomographic methods that achieve greater accuracy. In the case of mixed states of a single 2-dimensional quantum system adaptive methods have been recently introduced that achieve the theoretical accuracy limit deduced by Hayashi and Gill and Massar. However, accurate estimation of higher-dimensional quantum states remains poorly understood. This is mainly due to the existence of incompatible observables, which makes multiparameter estimation difficult. Here we present an adaptive tomographic method and show through numerical simulations that, after a few iterations, it is asymptotically approaching the fundamental Gill-Massar lower bound for the estimation accuracy of pure quantum states in high dimension. The method is based on a combination of stochastic optimization on the field of the complex numbers and statistical inference, exceeds the accuracy of any mixed-state tomographic method, and can be demonstrated with current experimental capabilities. The proposed method may lead to new developments in quantum metrology.


title: Sheaves of Twisted Cherednik Algebras as Universal Filtered Formal Deformations

abstract: According to a statement by Pavel Etingof, in the special case of an affine variety $X$ with a faithful action by a finite group $G$, the sheaf of (twisted) Cherednik algebras $\mathcal{H}_{1, c, \psi, X, G}$ with formal parameters $c, \psi$ is a universal formal deformation of $\mathcal{D}_X\rtimes G$ where $\mathcal{D}_X$ is the sheaf of differential operators on $X$. In the current note, we generalize Etingof's result to the non-affine case. We prove that for a generic smooth analytic or algebraic variety $X$, the sheaf $\mathcal{H}_{1, c, \psi, X, G}$ with formal $c$ and $\psi$ is a universal filtered formal deformation of $\mathcal{D}_X\rtimes G$. To that aim, we first construct quasi-isomorphisms between the Hochschild (co)chain complex of $\mathcal{D}_X\rtimes G$ and the $G$-invariant part of the direct sum over all elements $g$ in $G$ of sheaves of holomorphic differential forms on the cotangent bundles of the $g$-fixed point submanifolds in $X$. Finally, we combine these quasi-isomorphisms with results from the theory of algebraic extensions for sheaves of filtered associative algebras to establish a bijective correspondence between the space of isomorphism classes of filtered infinitesimal deformations of $\mathcal{D}_X\rtimes G$ and the parameter space of $\mathcal{H}_{1, c, \psi, X, G}$.


title: Necessary condition on the weight for maximal and integral operators with rough kernels

abstract: Let $0\leq \alpha<n$, $m\in \mathbb{N}$ and let consider $T_{\alpha,m}$ be a of integral operator, given by kernel of the form $$K(x,y)=k_1(x-A_1y)k_2(x-A_2y)\dots k_m(x-A_my),$$ where $A_i$ are invertible matrices and each $k_i$ satisfies a fractional size and generalized fractional H\"ormander condition. In [Iba\~nez-Firnkorn, G. H., and Riveros, M. S. (2018). Certain fractional type operators with H\"ormander conditions. To appear in Ann. Acad. Sci. Fenn. Math.] it was proved that $T_{\alpha,m}$ is controlled in $L^p(w)$-norms, $w\in A_{\infty}$, by the sum of maximal operators $M_{A_i^{-1},\alpha}$. In this paper we present the class of weights $\mathcal{A}_{A,p,q}$, where $A$ is an invertible matrix. This class are the good weights for the weak-type estimate of $M_{A^{-1},\alpha}$. For certain kernels $k_i$ we can characterize the weights for the strong-type estimate of $T_{\alpha,m}$. Also, we give a the strong-type estimate using testing conditions.


title: Hybrid deep learning architecture for general disruption prediction across tokamaks

abstract: In this paper, we present a new deep learning disruption prediction algorithm based on important findings from explorative data analysis which effectively allows knowledge transfer from existing devices to new ones, thereby predicting disruptions using very limited disruptive data from the new devices. The explorative data analysis conducted via unsupervised clustering techniques confirms that time-sequence data are much better separators of disruptive and non-disruptive behavior than the instantaneous plasma state data with further advantageous implications for a sequence-based predictor. Based on such important findings, we have designed a new algorithm for multi-machine disruption prediction that achieves high predictive accuracy on the C-Mod (AUC=0.801), DIII-D (AUC=0.947) and EAST (AUC=0.973) tokamaks with limited hyperparameter tuning. Through numerical experiments, we show that boosted accuracy (AUC=0.959) is achieved on EAST predictions by including in the training only 20 disruptive discharges, thousands of non-disruptive discharges from EAST, and combining this with more than a thousand discharges from DIII-D and C-Mod. The improvement of predictive ability obtained by combining disruptive data from other devices is found to be true for all permutations of the three devices. Furthermore, by comparing the predictive performance of each individual numerical experiment, we find that non-disruptive data are machine-specific while disruptive data from multiple devices contain device-independent knowledge that can be used to inform predictions for disruptions occurring on a new device.


title: Schr\"odinger Operators with Thin Spectra

abstract: The determination of the spectrum of a Schr\"odinger operator is a fundamental problem in mathematical quantum mechanics. We discuss a series of results showing that Schr\"odinger operators can exhibit spectra that are remarkably thin in the sense of Lebesgue measure and fractal dimensions. We begin with a brief discussion of results in the periodic theory, and then move to a discussion of aperiodic models with thin spectra.


title: Finite Temperature Off-Diagonal Long-Range Order for Interacting Bosons

abstract: Characterizing the scaling with the total particle number ($N$) of the largest eigenvalue of the one--body density matrix ($\lambda_0$), provides informations on the occurrence of the off-diagonal long-range order (ODLRO) according to the Penrose-Onsager criterion. Setting $\lambda_0\sim N^{\mathcal{C}_0}$, then $\mathcal{C}_0=1$ corresponds to ODLRO. The intermediate case, $0<\mathcal{C}_0<1$, corresponds for translational invariant systems to the power-law decaying of (non-connected) correlation functions and it can be seen as identifying quasi-long-range order. The goal of the present paper is to characterize the ODLRO properties encoded in $\mathcal{C}_0$ [and in the corresponding quantities $\mathcal{C}_{k \neq 0}$ for excited natural orbitals] exhibited by homogeneous interacting bosonic systems at finite temperature for different dimensions. We show that $\mathcal{C}_{k \neq 0}=0$ in the thermodynamic limit. In $1D$ it is $\mathcal{C}_0=0$ for non-vanishing temperature, while in $3D$ $\mathcal{C}_0=1$ ($\mathcal{C}_0=0$) for temperatures smaller (larger) than the Bose-Einstein critical temperature. We then focus our attention to $D=2$, studying the $XY$ and the Villain models, and the weakly interacting Bose gas. The universal value of $\mathcal{C}_0$ near the Berezinskii--Kosterlitz--Thouless temperature $T_{BKT}$ is $7/8$. The dependence of $\mathcal{C}_0$ on temperatures between $T=0$ (at which $\mathcal{C}_0=1$) and $T_{BKT}$ is studied in the different models. An estimate for the (non-perturbative) parameter $\xi$ entering the equation of state of the $2D$ Bose gases, is obtained using low temperature expansions and compared with the Monte Carlo result. We finally discuss a double jump behaviour for $\mathcal{C}_0$, and correspondingly of the anomalous dimension $\eta$, right below $T_{BKT}$ in the limit of vanishing interactions.


title: Crowdfunding for Design Innovation: Prediction Model with Critical Factors

abstract: Online reward-based crowdfunding campaigns have emerged as an innovative approach for validating demands, discovering early adopters, and seeking learning and feedback in the design processes of innovative products. However, crowdfunding campaigns for innovative products are faced with a high degree of uncertainty and suffer meager rates of success to fulfill their values for design. To guide designers and innovators for crowdfunding campaigns, this paper presents a data-driven methodology to build a prediction model with critical factors for crowdfunding success, based on public online crowdfunding campaign data. Specifically, the methodology filters 26 candidate factors in the Real-Win-Worth framework and identifies the critical ones via step-wise regression to predict the amount of crowdfunding. We demonstrate the methodology via deriving prediction models and identifying essential factors from 3D printer and smartwatch campaign data on Kickstarter and Indiegogo. The critical factors can guide campaign developments, and the prediction model may evaluate crowdfunding potential of innovations in contexts, to increase the chance of crowdfunding success of innovative products.


title: Rigidity results for automorphisms of Hardy-Toeplitz $C^*$-algebras

abstract: We prove a number of results on the automorphisms of and isomorphisms between Hardy-Toeplitz algebras $\mathcal{T}(D)$ associated to bounded symmetric domains $D$: that the stable isomorphism class of $\mathcal{T}(D)$ determines $D$ (even when it is reducible), that for reducible domains $D=D_1\times\cdots \times D_s$ the automorphisms of the Shilov boundary $\check{S}(D)$ induced by those of $\mathcal{T}(D)$ permute the Shilov boundaries $\check{S}(D_i)$, and that by contrast to arbitrary solvable algebras, automorphisms of $\mathcal{T}(D)$ that are trivial on their character spaces $\check{S}(D)$ are trivial on the entire spectrum $\widehat{\mathcal{T}(D)}$.


title: Radial regular and rupture solutions for a MEMS model with fringing field

abstract: We investigate radial solutions for the problem \[ \begin{cases} \displaystyle -\Delta U=\frac{\lambda+\delta|\nabla U|^2}{1-U},\; U>0 & \textrm{in}\ B,\\ U=0 & \textrm{on}\ \partial B, \end{cases} \] which is related to the study of Micro-Electromechanical Systems (MEMS). Here, $B\subset \mathbb{R}^N$ $(N\geq 2)$ denotes the open unit ball and $\lambda, \delta>0$ are real numbers. Two classes of solutions are considered in this work: (i) {\it regular solutions}, which satisfy $0<U<1$ in $B$ and (ii) {\it rupture solutions} which satisfy $U(0)=1$, and thus make the equation singular at the origin. Bifurcation with respect to parameter $\lambda>0$ is also discussed.


title: Theoretical approaches for nanoscale thermoelectric phenomena

abstract: Focus of the chapter is on the theoretical approaches aimed to analyze thermoelectric properties at the nanoscale. We discuss several relevant theoretical approaches for different set-ups of nano-devices providing estimations of the thermoelectric parameters in the linear and non-linear regime, in particular the thermoelectric figure of merit and the power-efficiency trade-off. Moreover, we analyze the role of not only electronic, but also of vibrational degrees of freedom. First, nanoscale thermoelectric phenomena are considered in the quantum coherent regime using the Landauer-B\"uttiker method and focusing on effects of energy filtering. Then, we analyze the effects of many-body couplings between nano-structure degrees of freedom, such as electron-electron and electron-vibration interactions, which can strongly affect the thermoelectric conversion. In particular, we discuss the enhancement of the thermoelectric figure of merit in the Coulomb blockade regime for a quantum dot model starting from the master equation for charge state probabilities and the tunneling rates through the electrodes. Finally, within the non-equilibrium Green function formalism, we quantify the reduction of the thermoelectric performance in simple models of molecular junctions due to the effects of the electron-vibration coupling and phonon transport at room temperature.


title: Creating a content delivery network for general science on the internet backbone using XCaches

abstract: A general problem faced by computing on the grid for opportunistic users is that delivering cycles is simpler than delivering data to those cycles. In this project we show how we integrated XRootD caches placed on the internet backbone to implement a content delivery network for general science workflows. We will show that for some workflows on different science domains like high energy physics, gravitational waves, and others the combination of data reuse from the workflows together with the use of caches increases CPU efficiency while decreasing network bandwidth use.


title: A (Slightly) Improved Approximation Algorithm for Metric TSP

abstract: For some $\epsilon > 10^{-36}$ we give a randomized $3/2-\epsilon$ approximation algorithm for metric TSP.


title: Weighted estimates of the Cayley transform method for boundary value problems in a Banach space

abstract: We consider the boundary value problems (BVPs) for linear secondorder ODEs with a strongly positive operator coefficient in a Banach space. The solutions are given in the form of the infinite series by means of the Cayley transform of the operator, the Meixner type polynomials of the independent variable, the operator Green function and the Fourier series representation for the right-hand side of the equation. The approximate solution of each problem is a partial sum of N (or expressed through N) summands. We prove the weighted error estimates depending on the discretization parameter N, the distance of the independent variable to the boundary points of the interval and some smoothness properties of the input data.


title: Hospitalization dynamics during the first COVID-19 pandemic wave: SIR modelling compared to Belgium, France, Italy, Switzerland and New York City data

abstract: Using the classical Susceptible-Infected-Recovered epidemiological model, an analytical formula is derived for the number of beds occupied by Covid-19 patients. The analytical curve is fitted to data in Belgium, France, New York City and Switzerland, with a correlation coefficient exceeding 98.8%, suggesting that finer models are unnecessary with such macroscopic data. The fitting is used to extract estimates of the doubling time in the ascending phase of the epidemic, the mean recovery time and, for those who require medical intervention, the mean hospitalization time. Large variations can be observed among different outbreaks.


title: Asymptotics and estimates for spectral minimal partitions of metric graphs

abstract: We study properties of spectral minimal partitions of metric graphs within the framework recently introduced in [Kennedy et al, Calc. Var. 60 (2021), 61]. We provide sharp lower and upper estimates for minimal partition energies in different classes of partitions; while the lower bounds are reminiscent of the classic isoperimetric inequalities for metric graphs, the upper bounds are more involved and mirror the combinatorial structure of the metric graph as well. Combining them, we deduce that these spectral minimal energies also satisfy a Weyl-type asymptotic law similar to the well-known one for eigenvalues of quantum graph Laplacians with various vertex conditions. Drawing on two examples we show that in general no second term in the asymptotic expansion for minimal partition energies can exist, but show that various kinds of behaviour are possible. We also study certain aspects of the asymptotic behaviour of the minimal partitions themselves.


title: Wearable Respiration Monitoring: Interpretable Inference with Context and Sensor Biomarkers

abstract: Breathing rate (BR), minute ventilation (VE), and other respiratory parameters are essential for real-time patient monitoring in many acute health conditions, such as asthma. The clinical standard for measuring respiration, namely Spirometry, is hardly suitable for continuous use. Wearables can track many physiological signals, like ECG and motion, yet not respiration. Deriving respiration from other modalities has become an area of active research. In this work, we infer respiratory parameters from wearable ECG and wrist motion signals. We propose a modular and generalizable classification-regression pipeline to utilize available context information, such as physical activity, in learning context-conditioned inference models. Morphological and power domain novel features from the wearable ECG are extracted to use with these models. Exploratory feature selection methods are incorporated in this pipeline to discover application-specific interpretable biomarkers. Using data from 15 subjects, we evaluate two implementations of the proposed pipeline: for inferring BR and VE. Each implementation compares generalized linear model, random forest, support vector machine, Gaussian process regression, and neighborhood component analysis as contextual regression models. Permutation, regularization, and relevance determination methods are used to rank the ECG features to identify robust ECG biomarkers across models and activities. This work demonstrates the potential of wearable sensors not only in continuous monitoring, but also in designing biomarker-driven preventive measures.


title: Minkowski gauges and deviation measures

abstract: We propose to derive deviation measures through the Minkowski gauge of a given set of acceptable positions. We show that, given a suitable acceptance set, any positive homogeneous deviation measure can be accommodated in our framework. In doing so, we provide a new interpretation for such measures, namely, that they quantify how much one must shrink or deleverage a position for it to become acceptable. In particular, the Minkowski Deviation of a set which is convex, stable under scalar addition, and radially bounded at non-constants, is a generalized deviation measure. Furthermore, we explore the relations existing between mathematical and financial properties attributable to an acceptance set, and the corresponding properties of the induced measure. Hence, we fill the gap that is the lack of an acceptance set for deviation measures. Dual characterizations in terms of polar sets and support functionals are provided.


title: The Supernova Remnant Population of NGC6946 as Observed in [Fe II] 1.644 $\mu$m with HST

abstract: NGC6946 is a high star formation rate face-on spiral galaxy that has hosted ten supernovae since 1917. Not surprisingly, a large number of supernova remnants and candidates have been identified either as optical nebulae with high [S II]:H$\alpha$ line ratios (147) or as compact non-thermal radio sources (35). However, there are only seven overlaps between these two samples. Here, we apply [Fe II] 1.644 $\mu$m emission as a new diagnostic to search for supernova remnants in an attempt to resolve this discrepancy. [Fe II] is expected to be relatively strong in the radiative shocks of supernova remnants and almost absent in HII regions. It is less susceptible to the effects of absorption along the line of sight than the optical lines normally used to identify remnants. Using data from the WFC3 camera on HST}, we identify 132 [Fe II] emission nebulae in NGC6946 as likely supernova remnants. Of these, 54 align with previously known optical supernova remnants. The remaining 78 objects are new; of these 44 are visible in new HST imagery in H$\alpha$ and [S II]. This brings the total number of supernova remnant candidates (from optical and/or IR data) in NGC6946 to 225. A total of 14 coincidences with radio supernova remnant candidates (out of 30 in our search area) are found in this expanded list. The identification of so many new remnant candidates validates the use of [Fe II] imagery for finding remnants, and suggests that previous remnant searches in other galaxies may be far from complete.


title: An order-adaptive compact approximation Taylor method for systems of conservation laws

abstract: We present a new family of high-order shock-capturing finite difference numerical methods for systems of conservation laws. These methods, called Adaptive Compact Approximation Taylor (ACAT) schemes, use centered $(2p + 1)$-point stencils, where $p$ may take values in $\{1, 2, \dots, P\}$ according to a new family of smoothness indicators in the stencils. The methods are based on a combination of a robust first order scheme and the Compact Approximate Taylor (CAT) methods of order $2p$-order, $p=1,2,\dots, P$ so that they are first order accurate near discontinuities and have order $2p$ in smooth regions, where $(2p +1)$ is the size of the biggest stencil in which large gradients are not detected. CAT methods, introduced in \cite{CP2019}, are an extension to nonlinear problems of the Lax-Wendroff methods in which the Cauchy-Kovalesky (CK) procedure is circumvented following the strategy introduced in \cite{ZBM2017} that allows one to compute time derivatives in a recursive way using high-order centered differentiation formulas combined with Taylor expansions in time. The expression of ACAT methods for 1D and 2D systems of balance laws are given and the performance is tested in a number of test cases for several linear and nonlinear systems of conservation laws, including Euler equations for gas dynamics.


title: Approximate Quantum Circuit Synthesis using Block-Encodings

abstract: One of the challenges in quantum computing is the synthesis of unitary operators into quantum circuits with polylogarithmic gate complexity. Exact synthesis of generic unitaries requires an exponential number of gates in general. We propose a novel approximate quantum circuit synthesis technique by relaxing the unitary constraints and interchanging them for ancilla qubits via block-encodings. This approach combines smaller block-encodings, which are easier to synthesize, into quantum circuits for larger operators. Due to the use of block-encodings, our technique is not limited to unitary operators and can also be applied for the synthesis of arbitrary operators. We show that operators which can be approximated by a canonical polyadic expression with a polylogarithmic number of terms can be synthesized with polylogarithmic gate complexity with respect to the matrix dimension.


title: Learning Orientation Distributions for Object Pose Estimation

abstract: For robots to operate robustly in the real world, they should be aware of their uncertainty. However, most methods for object pose estimation return a single point estimate of the object's pose. In this work, we propose two learned methods for estimating a distribution over an object's orientation. Our methods take into account both the inaccuracies in the pose estimation as well as the object symmetries. Our first method, which regresses from deep learned features to an isotropic Bingham distribution, gives the best performance for orientation distribution estimation for non-symmetric objects. Our second method learns to compare deep features and generates a non-parameteric histogram distribution. This method gives the best performance on objects with unknown symmetries, accurately modeling both symmetric and non-symmetric objects, without any requirement of symmetry annotation. We show that both of these methods can be used to augment an existing pose estimator. Our evaluation compares our methods to a large number of baseline approaches for uncertainty estimation across a variety of different types of objects.


title: Persistent Neurons

abstract: Neural networks (NN)-based learning algorithms are strongly affected by the choices of initialization and data distribution. Different optimization strategies have been proposed for improving the learning trajectory and finding a better optima. However, designing improved optimization strategies is a difficult task under the conventional landscape view. Here, we propose persistent neurons, a trajectory-based strategy that optimizes the learning task using information from previous converged solutions. More precisely, we utilize the end of trajectories and let the parameters explore new landscapes by penalizing the model from converging to the previous solutions under the same initialization. Persistent neurons can be regarded as a stochastic gradient method with informed bias where individual updates are corrupted by deterministic error terms. Specifically, we show that persistent neurons, under certain data distribution, is able to converge to more optimal solutions while initializations under popular framework find bad local minima. We further demonstrate that persistent neurons helps improve the model's performance under both good and poor initializations. We evaluate the full and partial persistent model and show it can be used to boost the performance on a range of NN structures, such as AlexNet and residual neural network (ResNet).


title: CoPhy-PGNN: Learning Physics-guided Neural Networks with Competing Loss Functions for Solving Eigenvalue Problems

abstract: Physics-guided Neural Networks (PGNNs) represent an emerging class of neural networks that are trained using physics-guided (PG) loss functions (capturing violations in network outputs with known physics), along with the supervision contained in data. Existing work in PGNNs has demonstrated the efficacy of adding single PG loss functions in the neural network objectives, using constant trade-off parameters, to ensure better generalizability. However, in the presence of multiple PG functions with competing gradient directions, there is a need to adaptively tune the contribution of different PG loss functions during the course of training to arrive at generalizable solutions. We demonstrate the presence of competing PG losses in the generic neural network problem of solving for the lowest (or highest) eigenvector of a physics-based eigenvalue equation, which is commonly encountered in many scientific problems. We present a novel approach to handle competing PG losses and demonstrate its efficacy in learning generalizable solutions in two motivating applications of quantum mechanics and electromagnetic propagation. All the code and data used in this work is available at https://github.com/jayroxis/Cophy-PGNN.


title: Semi-Supervised Training of Optical Flow Convolutional Neural Networks in Ultrasound Elastography

abstract: Convolutional Neural Networks (CNN) have been found to have great potential in optical flow problems thanks to an abundance of data available for training a deep network. The displacement estimation step in UltraSound Elastography (USE) can be viewed as an optical flow problem. Despite the high performance of CNNs in optical flow, they have been rarely used for USE due to unique challenges that both input and output of USE networks impose. Ultrasound data has much higher high-frequency content compared to natural images. The outputs are also drastically different, where displacement values in USE are often smooth without sharp motions or discontinuities. The general trend is currently to use pre-trained networks and fine-tune them on a small simulation ultrasound database. However, realistic ultrasound simulation is computationally expensive. Also, the simulation techniques do not model complex motions, nonlinear and frequency-dependent acoustics, and many sources of artifact in ultrasound imaging. Herein, we propose an unsupervised fine-tuning technique which enables us to employ a large unlabeled dataset for fine-tuning of a CNN optical flow network. We show that the proposed unsupervised fine-tuning method substantially improves the performance of the network and reduces the artifacts generated by networks trained on computer vision databases.


title: Driven-dissipative phase transition in a Kerr oscillator: From semiclassical $\mathcal{PT}$ symmetry to quantum fluctuations

abstract: We study a minimal model that has a driven-dissipative quantum phase transition, namely a Kerr non-linear oscillator subject to driving and dissipation. Using mean-field theory, exact diagonalization, and the Keldysh formalism, we analyze the critical phenomena in this system, showing which aspects can be captured by each approach and how the approaches complement each other. Then critical scaling and finite-size scaling are calculated analytically using the quantum Langevin equation. The physics contained in this simple model is surprisingly rich: it includes a continuous phase transition, $Z_{2}$ symmetry breaking, $\mathcal{PT}$ symmetry, state squeezing, and critical fluctuations. Due to its simplicity and solvability, this model can serve as a paradigm for exploration of open quantum many-body physics.


title: Maximizing Cohesion and Separation in Graph Representation Learning: A Distance-aware Negative Sampling Approach

abstract: The objective of unsupervised graph representation learning (GRL) is to learn a low-dimensional space of node embeddings that reflect the structure of a given unlabeled graph. Existing algorithms for this task rely on negative sampling objectives that maximize the similarity in node embeddings at nearby nodes (referred to as "cohesion") by maintaining positive and negative corpus of node pairs. While positive samples are drawn from node pairs that co-occur in short random walks, conventional approaches construct negative corpus by uniformly sampling random pairs, thus ignoring valuable information about structural dissimilarity among distant node pairs (referred to as "separation"). In this paper, we present a novel Distance-aware Negative Sampling (DNS) which maximizes the separation of distant node-pairs while maximizing cohesion at nearby node-pairs by setting the negative sampling probability proportional to the pair-wise shortest distances. Our approach can be used in conjunction with any GRL algorithm and we demonstrate the efficacy of our approach over baseline negative sampling methods over downstream node classification tasks on a number of benchmark datasets and GRL algorithms. All our codes and datasets are available at https://github.com/Distance-awareNS/DNS/.


title: Active Control and Sustained Oscillations in actSIS Epidemic Dynamics

abstract: An actively controlled Susceptible-Infected-Susceptible (actSIS) contagion model is presented for studying epidemic dynamics with continuous-time feedback control of infection rates. Our work is inspired by the observation that epidemics can be controlled through decentralized disease-control strategies such as quarantining, sheltering in place, social distancing, etc., where individuals actively modify their contact rates with others in response to observations of infection levels in the population. Accounting for a time lag in observations and categorizing individuals into distinct sub-populations based on their risk profiles, we show that the actSIS model manifests qualitatively different features as compared with the SIS model. In a homogeneous population of risk-averters, the endemic equilibrium is always reduced, although the transient infection level can exhibit overshoot or undershoot. In a homogeneous population of risk-tolerating individuals, the system exhibits bistability, which can also lead to reduced infection. For a heterogeneous population comprised of risk-tolerators and risk-averters, we prove conditions on model parameters for the existence of a Hopf bifurcation and sustained oscillations in the infected population.


title: Continuous Time Limit of the DTQW in 2D+1 and Plasticity

abstract: A Plastic Quantum Walk admits both continuous time and continuous spacetime. The model has been recently proposed by one of the authors in \cite{molfetta2019quantum}, leading to a general quantum simulation scheme for simulating fermions in the relativistic and non relativistic regimes. The extension to two physical dimensions is still missing and here, as a novel result, we demonstrate necessary and sufficient conditions concerning which discrete time quantum walks can admit plasticity, showing the resulting Hamiltonians. We consider coin operators as general $4$ parameter unitary matrices, with parameters which are function of the lattice step size $\varepsilon$. This dependence on $\varepsilon$ encapsulates all functions of $\varepsilon$ for which a Taylor series expansion in $\varepsilon$ is well defined, making our results very general.


title: Risk Modelling on Liquidations with L\'{e}vy Processes

abstract: It has been decades since the academic world of ruin theory defined the insolvency of an insurance company as the time when its surplus falls below zero. This simplification, however, needs careful adaptions to imitate the real-world liquidation process. Inspired by Broadie et al. (2007) and Li et al. (2020), this paper uses a three-barrier model to describe the financial stress towards bankruptcy of an insurance company. The financial status of the insurer is divided into solvent, insolvent and liquidated three states, where the insurer's surplus process at the state of solvent and insolvent is modelled by two spectrally negative L\'{e}vy processes, which have been taken as good candidates to model insurance risks. We provide a rigorous definition of the time of liquidation ruin in this three-barrier model. By adopting the techniques of excursions in the fluctuation theory, we study the joint distribution of the time of liquidation, the surplus at liquidation and the historical high of the surplus until liquidation, which generalizes the known results on the classical expected discounted penalty function in Gerber and Shiu (1998). The results have semi-explicit expressions in terms of the scale functions and the L\'{e}vy triplets associated with the two underlying L\'{e}vy processes. The special case when the two underlying L\'{e}vy processes coincide with each other is also studied, where our results are expressed compactly via only the scale functions. The corresponding results have good consistency with the existing literatures on Parisian ruin with (or without) a lower barrier in Landriault et al. (2014), Baurdoux et al. (2016) and Frostig and Keren-Pinhasik (2019). Besides, numerical examples are provided to illustrate the underlying features of liquidation ruin.


title: Atomic-Scale Defect Detection by Nonlinear Light Scattering and Localization

abstract: Hetero-epitaxial crystalline films underlie many electronic and optical technologies but are prone to forming defects at their hetero-interfaces. Atomic-scale defects such as threading dislocations that propagate into a film impede the flow of charge carriers and light degrading electrical-optical performance of devices. Diagnosis of subsurface defects traditionally requires time consuming invasive techniques such as cross sectional transmission electron microscopy. Using III-V films grown on Si, we have demonstrated noninvasive, bench-top diagnosis of sub-surface defects by optical second-harmonic scanning probe microscope. We observed a high-contrast pattern of sub-wavelength hot spots caused by scattering and localization of fundamental light by defect scattering sites. Size of these observed hotspots are strongly correlated to the density of dislocation defects. Our results not only demonstrate a global and versatile method for diagnosing sub-surface scattering sites but uniquely elucidate optical properties of disordered media. An extension to third harmonics would enable irregularities detection in non-X(2) materials making the technique universally applicable.


title: Medial Axis Isoperimetric Profiles

abstract: Recently proposed as a stable means of evaluating geometric compactness, the isoperimetric profile of a planar domain measures the minimum perimeter needed to inscribe a shape with prescribed area varying from 0 to the area of the domain. While this profile has proven valuable for evaluating properties of geographic partitions, existing algorithms for its computation rely on aggressive approximations and are still computationally expensive. In this paper, we propose a practical means of approximating the isoperimetric profile and show that for domains satisfying a "thick neck" condition, our approximation is exact. For more general domains, we show that our bound is still exact within a conservative regime and is otherwise an upper bound. Our method is based on a traversal of the medial axis which produces efficient and robust results. We compare our technique with the state-of-the-art approximation to the isoperimetric profile on a variety of domains and show significantly tighter bounds than were previously achievable.


title: The Global Landscape of Neural Networks: An Overview

abstract: One of the major concerns for neural network training is that the non-convexity of the associated loss functions may cause bad landscape. The recent success of neural networks suggests that their loss landscape is not too bad, but what specific results do we know about the landscape? In this article, we review recent findings and results on the global landscape of neural networks. First, we point out that wide neural nets may have sub-optimal local minima under certain assumptions. Second, we discuss a few rigorous results on the geometric properties of wide networks such as "no bad basin", and some modifications that eliminate sub-optimal local minima and/or decreasing paths to infinity. Third, we discuss visualization and empirical explorations of the landscape for practical neural nets. Finally, we briefly discuss some convergence results and their relation to landscape results.


title: Portfolio Optimization of 40 Stocks Using the DWave Quantum Annealer

abstract: We investigate the use of quantum computers for building a portfolio out of a universe of U.S. listed, liquid equities that contains an optimal set of stocks. Starting from historical market data, we look at various problem formulations on the D-Wave Systems Inc. D-Wave 2000Q(TM) System (hereafter called DWave) to find the optimal risk vs return portfolio; an optimized portfolio based on the Markowitz formulation and the Sharpe ratio, a simplified Chicago Quantum Ratio (CQR), then a new Chicago Quantum Net Score (CQNS). We approach this first classically, then by our new method on DWave. Our results show that practitioners can use a DWave to select attractive portfolios out of 40 U.S. liquid equities.


title: Graphs without gap-vertex-labellings: families and bounds

abstract: A proper labelling of a graph $G$ is a pair $({\pi},c_{\pi})$ in which ${\pi}$ is an assignment of numeric labels to some elements of $G$, and $c_{\pi}$ is a colouring induced by ${\pi}$ through some mathematical function over the set of labelled elements. In this work, we consider gap-vertex-labellings, in which the colour of a vertex is determined by a function considering the largest difference between the labels assigned to its neighbours. We present the first upper-bound for the vertex-gap number of arbitrary graphs, which is the least number of labels required to properly label a graph. We investigate families of graphs which do not admit any gap-vertex-labelling, regardless of the number of labels. Furthermore, we introduce a novel parameter associated with this labelling and provide bounds for it for complete graphs ${K_n}$.


title: Probabilistic Inversions for Time-Distance Helioseismology

abstract: Time-distance helioseismology is a set of powerful tools to study features below the Sun's surface. Inverse methods are needed to interpret time-distance measurements, with many examples in the literature. However, techniques that utilize a more statistical approach to inferences, and broadly used in the astronomical community, are less-commonly found in helioseismology. This article aims to introduce a potentially powerful inversion scheme based on Bayesian probability theory and Monte Carlo sampling that is suitable for local helioseismology. We describe the probabilistic method and how it is conceptually different from standard inversions used in local helioseismology. Several example calculations are carried out to compare and contrast the setup of the problems and the results that are obtained. The examples focus on two important phenomena studied with helioseismology: meridional circulation and supergranulation. Numerical models are used to compute synthetic observations, providing the added benefit of knowing the solution against which the results can be tested. For demonstration purposes, the problems are formulated in two and three dimensions, using both ray- and Born-theoretical approaches. The results seem to indicate that the probabilistic inversions not only find a better solution with much more realistic estimation of the uncertainties, but they also provide a broader view of the range of solutions possible for any given model, making the interpretation of the inversion more quantitative in nature. Unlike the progress being made in fundamental measurement schemes in local helioseismology that image the far side of the Sun, or have detected signatures of global Rossby waves, among many others, inversions of those measurements have had significantly less success. Such statistical methods may help overcome some of these barriers to move the field forward.


title: Two-dimensional thermal finite element model of directed energy deposition: matching melt pool temperature profile with pyrometer measurement

abstract: An open source two-dimensional (2D) thermal finite element (FE) model of the Directed Energy Deposition (DED) process is developed using the Python-based FEniCS framework. The model incrementally deposits material ahead of the laser focus point according to the geometry of the part. The laser heat energy is supplied by a Gaussian-distributed heat source while the phase change is represented by increased heat capacity around the solidus-liquidus temperature range. Experimental validation of the numerical model is performed by matching with the melt pool temperature measurements taken by a dual wavelength pyrometer during the build process of a box-shaped Ti--6Al--4V part with large geometrical voids. Effects of large geometrical voids on the melt pool shape and maximum melt pool temperature are examined. Both the numerical and experimental data show an increase in the melt pool size and temperature during deposition above large voids. The trailing edge of the melt pool's temperature profile obtained using the developed numerical model closely matches pyrometer measurements.


title: In Search of Lost Domain Generalization

abstract: The goal of domain generalization algorithms is to predict well on distributions different from those seen during training. While a myriad of domain generalization algorithms exist, inconsistencies in experimental conditions -- datasets, architectures, and model selection criteria -- render fair and realistic comparisons difficult. In this paper, we are interested in understanding how useful domain generalization algorithms are in realistic settings. As a first step, we realize that model selection is non-trivial for domain generalization tasks. Contrary to prior work, we argue that domain generalization algorithms without a model selection strategy should be regarded as incomplete. Next, we implement DomainBed, a testbed for domain generalization including seven multi-domain datasets, nine baseline algorithms, and three model selection criteria. We conduct extensive experiments using DomainBed and find that, when carefully implemented, empirical risk minimization shows state-of-the-art performance across all datasets. Looking forward, we hope that the release of DomainBed, along with contributions from fellow researchers, will streamline reproducible and rigorous research in domain generalization.


title: A variational formulation for motion design of adaptive compliant structures

abstract: Adaptive structures are characterized by their ability to adjust their geometrical and other properties to changing loads or requirements during service. This contribution deals with a method for the design of quasi-static motions of structures between two prescribed geometrical configurations that are optimal with regard to a specified quality function while taking large deformations into account. It is based on a variational formulation and the solution by two finite element discretizations, the spatial discretization (the standard finite element mesh) and an additional discretization of the deformation path or trajectory. For the investigations, an exemplary objective function, the minimization of the internal energy, integrated along the deformation path, is used. The method for motion design presented herein uses the Newton-Raphson method as a second order optimization algorithm and allows for analytical sensitivity analysis. The proposed method is verified and its properties are investigated by benchmark examples including rigid body motions, instability phenomena and determination of inextensible deformations of shells.


title: Attribution Methods Reveal Flaws in Fingerprint-Based Virtual Screening

abstract: Fingerprint-based models for protein-ligand binding have demonstrated outstanding success on benchmark datasets; however, these models may not learn the correct binding rules. To assess this concern, we use in silico datasets with known binding rules to develop a general framework for evaluating model attribution. This framework identifies fragments that a model considers necessary to achieve a particular score, sidestepping the need for a model to be differentiable. Our results confirm that high-performing models may not learn the correct binding rule, and suggest concrete steps that can remedy this situation. We show that adding fragment-matched inactive molecules (decoys) to the data reduces attribution false negatives, while attribution false positives largely arise from the background correlation structure of molecular data. Normalizing for these background correlations helps to reveal the true binding logic. Our work highlights the danger of trusting attributions from high-performing models and suggests that a closer examination of fingerprint correlation structure and better decoy selection may help reduce misattributions.


title: Genuine Parrondo's paradox in quantum walks with time-dependent coin operators

abstract: We show that a genuine Parrondo paradox can emerge in two-state quantum walks without resorting to experimentally intricate high-dimensional coins. To achieve such goal we employ a time-dependent coin operator without breaking the translation spatial invariance of the system.


title: Hyper-massive Black Holes have Faint Broad and Narrow Emission Lines

abstract: The EUV provides most of the ionization that creates the high equivalent width (EW) broad and narrow emission lines (BELs, NELs) of quasars. Spectra of Hypermassive Schwarzschild black holes (HMBHs, $M_{BH} \geq 10^{10} M_{\odot}$) with $\alpha$-discs, decline rapidly in the EUV suggesting much lower EWs. Model spectra for black holes of mass $10^{6}-10^{12} M_{\odot}$ and accretion rates $0.03 \leq L_{bol}/L_{edd} \leq 1.0$ were input to the CLOUDY photoionization code. BELs become $\sim$100 times weaker in EW from $M_{BH} \sim 10^8 M_{\odot}$ to $M_{BH} \sim 10^{10} M_{\odot}$. The high ionization BELs (O VI 1034 $\overset{\circ}{\mathrm {A}}$, C IV 1549 $\overset{\circ}{\mathrm {A}}$, He II 1640 $\overset{\circ}{\mathrm {A}}$) decline in EW from ($M_{BH} \geq 10^6 M_{\odot}$, reproducing the Baldwin effect, but regain EW for $M_{BH} \geq 10^{10} M_{\odot}$). The low ionization lines (MgII 2798 $\overset{\circ}{\mathrm {A}}$, H$\beta$ 4861 $\overset{\circ}{\mathrm {A}}$ and H$\alpha$ 6563 $\overset{\circ}{\mathrm {A}}$) remain weak. Lines for maximally spinning HMBHs behave similarly. Line ratio diagrams for the BELs show that high OVI/H$\beta$ and low CIV/H$\alpha$ may pick out HMBH, although OVI is often hard to observe. In NEL BPT diagrams HMBHs lie among star-forming regions, except for highly spinning, high accretion rate HMBHs. In summary, the BELs expected from HMBHs would be hard to detect using the current optical facilities. From 100 to $10^{12} M_{\odot}$, the emission lines used to detect AGN only have high EW in the $10^6 - 10^9 M_{\odot}$ window, where most AGN are found. This selection effect may be distorting reported distributions of $M_{BH}$.


title: Computational modeling of extreme wildland fire events: a synthesis of scientific understanding with applications to forecasting, land management, and firefighter safety

abstract: The understanding and prediction of large wildland fire events around the world is a growing interdisciplinary research area advanced rapidly by development and use of computational models. Recent models bidirectionally couple computational fluid dynamics models including weather prediction models with modules containing algorithms representing fire spread and heat release, simulating fire-atmosphere interactions across scales spanning three orders of magnitude. Integrated with weather data and airborne and satellite remote sensing data on wildland fuels and active fire detection, modern coupled weather-fire modeling systems are being used to solve current science problems. Compared to legacy tools, these dynamic computational modeling systems increase cost and complexity but have produced breakthrough insights notably into the mechanisms underlying extreme wildfire events such as fine-scale extreme winds associated with interruptions of the electricity grid and have been configured to forecast a fire's growth, expanding our ability to anticipate how they will unfold. We synthesize case studies of recent extreme events, expanding applications, and the challenges and limitations in our remote sensing systems, fire prediction tools, and meteorological models that add to wildfires' mystery and apparent unpredictability.


title: Searching for axion-like particles under strong gravitational lenses

abstract: We establish strong gravitational lens systems as robust probes of axion-like particles (ALPs) -- a candidate for dark matter. A tiny interaction of photons with ALPs induces birefringence. Multiple images of gravitationally lensed polarised objects allow differential birefringence measurement, alleviating systematics and astrophysical dependencies. We apply this novel method to the lens system CLASS B1152+199 and constrain ALP-photon coupling $\le 9.2\times 10^{-11}\, {\rm GeV}^{-1} \textrm{ to } 7.7\times 10^{-8}\, {\rm GeV}^{-1}$ ($95\%$ C.L.) for ALP mass between $3.6\times 10^{-21} \,{\rm eV}$ and $4.6 \times 10^{-18} \,{\rm eV}$. A larger sample will improve the constraints.


title: Joint Frequency and Image Space Learning for MRI Reconstruction and Analysis

abstract: We propose neural network layers that explicitly combine frequency and image feature representations and show that they can be used as a versatile building block for reconstruction from frequency space data. Our work is motivated by the challenges arising in MRI acquisition where the signal is a corrupted Fourier transform of the desired image. The proposed joint learning schemes enable both correction of artifacts native to the frequency space and manipulation of image space representations to reconstruct coherent image structures at every layer of the network. This is in contrast to most current deep learning approaches for image reconstruction that treat frequency and image space features separately and often operate exclusively in one of the two spaces. We demonstrate the advantages of joint convolutional learning for a variety of tasks, including motion correction, denoising, reconstruction from undersampled acquisitions, and combined undersampling and motion correction on simulated and real world multicoil MRI data. The joint models produce consistently high quality output images across all tasks and datasets. When integrated into a state of the art unrolled optimization network with physics-inspired data consistency constraints for undersampled reconstruction, the proposed architectures significantly improve the optimization landscape, which yields an order of magnitude reduction of training time. This result suggests that joint representations are particularly well suited for MRI signals in deep learning networks. Our code and pretrained models are publicly available at https://github.com/nalinimsingh/interlacer.


title: Multi-Agent Low-Dimensional Linear Bandits

abstract: We study a multi-agent stochastic linear bandit with side information, parameterized by an unknown vector $\theta^* \in \mathbb{R}^d$. The side information consists of a finite collection of low-dimensional subspaces, one of which contains $\theta^*$. In our setting, agents can collaborate to reduce regret by sending recommendations across a communication graph connecting them. We present a novel decentralized algorithm, where agents communicate subspace indices with each other and each agent plays a projected variant of LinUCB on the corresponding (low-dimensional) subspace. By distributing the search for the optimal subspace across users and learning of the unknown vector by each agent in the corresponding low-dimensional subspace, we show that the per-agent finite-time regret is much smaller than the case when agents do not communicate. We finally complement these results through simulations.


title: Machine Learning Approach for Transforming Scattering Parameters to Complex Permittivity

abstract: This study investigates the application of an artificial neural network to predict the complex dielectric properties of granular catalysts commonly used in microwave reaction chemistry. The study utilizes finite element electromagnetic simulations and two-dimensional convolutional neural networks to solve for a large solution space of varying dielectrics. This convolutional neural network was trained using a supervised learning approach and a common backpropagation. The frequency range of interest was between 0.1 to 13.5 GHz with the real part of the dielectric constants ranging from 1 to 100 and the imaginary part ranging from 0.0 to 0.2. The network was double validated using experimental data collected from a coaxial airline. The model was demonstrated to convert either experimental or computational derived scattering parameter to complex permittivities. Moreover, the model eliminates the need for iterative solutions that often have difficulty with the piecewise continuous nature of frequency dependent scattering parameters.


title: Generative Modeling for Atmospheric Convection

abstract: While cloud-resolving models can explicitly simulate the details of small-scale storm formation and morphology, these details are often ignored by climate models for lack of computational resources. Here, we explore the potential of generative modeling to cheaply recreate small-scale storms by designing and implementing a Variational Autoencoder (VAE) that performs structural replication, dimensionality reduction, and clustering of high-resolution vertical velocity fields. Trained on ~6*10^6 samples spanning the globe, the VAE successfully reconstructs the spatial structure of convection, performs unsupervised clustering of convective organization regimes, and identifies anomalous storm activity, confirming the potential of generative modeling to power stochastic parameterizations of convection in climate models.


title: Minimizing Convex Functions with Rational Minimizers

abstract: Given a separation oracle $\mathsf{SO}$ for a convex function $f$ defined on $\mathbb{R}^n$ that has an integral minimizer inside a box with radius $R$, we show how to find an exact minimizer of $f$ using at most (a) $O(n (n \log \log (n)/\log (n) + \log(R)))$ calls to $\mathsf{SO}$ and $\mathsf{poly}(n, \log(R))$ arithmetic operations, or (b) $O(n \log(nR))$ calls to $\mathsf{SO}$ and $\exp(O(n)) \cdot \mathsf{poly}(\log(R))$ arithmetic operations. When the set of minimizers of $f$ has integral extreme points, our algorithm outputs an integral minimizer of $f$. This improves upon the previously best oracle complexity of $O(n^2 (n + \log(R)))$ for polynomial time algorithms and $O(n^2\log(nR))$ for exponential time algorithms obtained by [Gr\"otschel, Lov\'asz and Schrijver, Prog. Comb. Opt. 1984, Springer 1988] over thirty years ago. Our improvement on Gr\"otschel, Lov\'asz and Schrijver's result generalizes to the setting where the set of minimizers of $f$ is a rational polyhedron with bounded vertex complexity. For the Submodular Function Minimization problem, our result immediately implies a strongly polynomial algorithm that makes at most $O(n^3 \log \log (n)/\log (n))$ calls to an evaluation oracle, and an exponential time algorithm that makes at most $O(n^2 \log(n))$ calls to an evaluation oracle. These improve upon the previously best $O(n^3 \log^2(n))$ oracle complexity for strongly polynomial algorithms given in [Lee, Sidford and Wong, FOCS 2015] and [Dadush, V\'egh and Zambelli, SODA 2018], and an exponential time algorithm with oracle complexity $O(n^3 \log(n))$ given in the former work. Our result is achieved via a reduction to the Shortest Vector Problem in lattices. We analyze its oracle complexity using a potential function that simultaneously captures the size of the search set and the density of the lattice.


title: Losing Oceans: The Effects of Composition on the Thermal Component of Impact-driven Atmospheric Loss

abstract: The formation of the solar system's terrestrial planets concluded with a period of giant impacts. Previous works examining the volatile loss caused by the impact shock in the moon-forming impact find atmospheric losses of at most 20-30 per cent and essentially no loss of oceans. However, giant impacts also result in thermal heating, which can lead to significant atmospheric escape via a Parker-type wind. Here we show that H2O and other high-mean molecular weight outgassed species can be efficiently lost through this thermal wind if present in a hydrogen-dominated atmosphere, substantially altering the final volatile inventory of terrestrial planets. Specifically, we demonstrate that a giant impact of a Mars-sized embryo with a proto-Earth can remove several Earth oceans' worth of H2O, and other heavier volatile species, together with a primordial hydrogen-dominated atmosphere. These results may offer an explanation for the observed depletion in Earth's light noble gas budget and for its depleted xenon inventory, which suggest that Earth underwent significant atmospheric loss by the end of its accretion. Because planetary embryos are massive enough to accrete primordial hydrogen envelopes and because giant impacts are stochastic and occur concurrently with other early atmospheric evolutionary processes, our results suggest a wide diversity in terrestrial planet volatile budgets.


title: Improved Preterm Prediction Based on Optimized Synthetic Sampling of EHG Signal

abstract: Preterm labor is the leading cause of neonatal morbidity and mortality and has attracted research efforts from many scientific areas. The inter-relationship between uterine contraction and the underlying electrical activities makes uterine electrohysterogram (EHG) a promising direction for preterm detection and prediction. Due the scarcity of EHG signals, especially those of preterm patients, synthetic algorithms are applied to create artificial samples of preterm type in order to remove prediction bias towards term, at the expense of a reduction of the feature effectiveness in machine-learning based automatic preterm detecting. To address such problem, we quantify the effect of synthetic samples (balance coefficient) on features' effectiveness, and form a general performance metric by utilizing multiple feature scores with relevant weights that describe their contributions to class separation. Combined with the activation/inactivation functions that characterizes the effect of the abundance of training samples in term and preterm prediction precision, we obtain an optimal sample balance coefficient that compromise the effect of synthetic samples in removing bias towards the majority and the side-effect of reducing features' importance. Substantial improvement in prediction precision has been achieved through a set of numerical tests on public available TPEHG database, and it verifies the effectiveness of the proposed method.


title: From Fear to Hate: How the Covid-19 Pandemic Sparks Racial Animus in the United States

abstract: We estimate the effect of the Coronavirus (Covid-19) pandemic on racial animus, as measured by Google searches and Twitter posts including a commonly used anti-Asian racial slur. Our empirical strategy exploits the plausibly exogenous variation in the timing of the first Covid-19 diagnosis across regions in the United States. We find that the first local diagnosis leads to an immediate increase in racist Google searches and Twitter posts, with the latter mainly coming from existing Twitter users posting the slur for the first time. This increase could indicate a rise in future hate crimes, as we document a strong correlation between the use of the slur and anti-Asian hate crimes using historic data. Moreover, we find that the rise in the animosity is directed at Asians rather than other minority groups and is stronger on days when the connection between the disease and Asians is more salient, as proxied by President Trump's tweets mentioning China and Covid-19 at the same time. In contrast, the negative economic impact of the pandemic plays little role in the initial increase in racial animus. Our results suggest that de-emphasizing the connection between the disease and a particular racial group can be effective in curbing current and future racial animus.


title: The Two-Dimensional Fractional Discrete Nonlinear Schrodinger Equation

abstract: We study a fractional version of the two-dimensional discrete nonlinear Schr\"{o}dinger (DNLS) equation, where the usual discrete Laplacian is replaced by its fractional form that depends on a fractional exponent $s$ that interpolates between the case of an identity operator ($s=0$) and that of the usual discrete 2D Laplacian ($s=1$). This replacement leads to a long-range coupling among sites that, at low values of $s$, decreases the bandwidth and leads to quasi-degenerate states. The mean square displacement of an initially-localized excitation is shown to be ballistic at all times with a `speed' that increases monotonically with the fractional exponent $s$. We also compute the nonlinear modes and their stability for both, bulk and surface modes. The modulational stability is seen to increase with an increase in the fractional exponent. The trapping of an initially localized excitation shows a selftrapping transition as a function of nonlinearity strength, whose threshold increases with the value of $s$. In the linear limit, there persists a linear trapping at small $s$ values. This behavior is connected with the decrease of the bandwidth and its associated increase in quasi-degeneracy.


title: Perfect submonoids of dominant weights

abstract: Let G be a connected semisimple group. Vinberg introduced the notion of perfect submonoids of dominant weights of G in the study of Vinberg monoids. In this paper, we give explicit descriptions of the perfect submonoids.


title: Improved RIP-Based Bounds for Guaranteed Performance of two Compressed Sensing Algorithms

abstract: Iterative hard thresholding (IHT) and compressive sampling matching pursuit (CoSaMP) are two types of mainstream compressed sensing algorithms using hard thresholding operators for signal recovery and approximation. The guaranteed performance for signal recovery via these algorithms has mainly been analyzed under the condition that the restricted isometry constant of a sensing matrix, denoted by $ \delta_K$ (where $K$ is an integer number), is smaller than a certain threshold value in the interval $(0,1).$ The condition $ \delta_{K}< \delta^*$ for some constant $ \delta^* \leq 1 $ ensuring the success of signal recovery with a specific algorithm is called the restricted-isometry-property-based (RIP-based) bound for guaranteed performance of the algorithm. At the moment, the best known RIP-based bound for the guaranteed recovery of $k$-sparse signals via IHT is $\delta_{3k}< 1/\sqrt{3}\approx 0.5774,$ and the bound for guaranteed recovery via CoSaMP is $\delta_{4k} < 0.4782. $ A fundamental question in this area is whether such theoretical results can be further improved. The purpose of this paper is to affirmatively answer this question and rigorously show that the RIP-based bounds for guaranteed performance of IHT can be significantly improved to $ \delta_{3k} < (\sqrt{5}-1)/2 \approx 0.618, $ and the bound for CoSaMP can be improved and pushed to $ \delta_{4k}< 0.5102. $ These improvements are achieved through a deep property of the hard thresholding operator.


title: Modeling from Features: a Mean-field Framework for Over-parameterized Deep Neural Networks

abstract: This paper proposes a new mean-field framework for over-parameterized deep neural networks (DNNs), which can be used to analyze neural network training. In this framework, a DNN is represented by probability measures and functions over its features (that is, the function values of the hidden units over the training data) in the continuous limit, instead of the neural network parameters as most existing studies have done. This new representation overcomes the degenerate situation where all the hidden units essentially have only one meaningful hidden unit in each middle layer, and further leads to a simpler representation of DNNs, for which the training objective can be reformulated as a convex optimization problem via suitable re-parameterization. Moreover, we construct a non-linear dynamics called neural feature flow, which captures the evolution of an over-parameterized DNN trained by Gradient Descent. We illustrate the framework via the standard DNN and the Residual Network (Res-Net) architectures. Furthermore, we show, for Res-Net, when the neural feature flow process converges, it reaches a global minimal solution under suitable conditions. Our analysis leads to the first global convergence proof for over-parameterized neural network training with more than $3$ layers in the mean-field regime.


title: Intrinsic bulk quantum oscillations in a bulk unconventional insulator SmB$_6$

abstract: The finding of bulk quantum oscillations in the bulk Kondo insulator SmB$_6$, which has been proposed to be a correlated topological insulator, proved a considerable surprise. The subsequent measurement of bulk quantum oscillations in other correlated insulators including YbB$_{12}$ have lent support to our discovery of a class of unconventional insulators that are host to bulk quantum oscillations, of which SmB$_6$ was the first example. Here we perform a series of experiments to examine evidence for the intrinsic character of bulk quantum oscillations in floating zone-grown single crystals of SmB$_6$ that have been the subject of our quantum oscillation studies thus far. We present results of experiments including chemical composition analysis, magnetisation, thermal conductivity, electrical transport, and heat capacity on floating zone-grown single crystals of SmB$_6$, and a series of quantum oscillation experiments as a function of magnetic field, temperature, and magnetic field-orientation on single crystals of floating-zone grown SmB$_6$, LaB$_6$, and elemental Aluminium. Results of these experimental studies establish the intrinsic origin of quantum oscillations from the bulk of pristine floating zone-grown single crystals of SmB$_6$. The origin of the underlying bulk Fermi surface that bears close similarity with the unhybridised Fermi surface in metallic hexaborides despite the bulk insulating character of SmB$_6$ is thus at the heart of a theoretical mystery.


title: A new fixed point approach to hyperstability of radical-type functional equations in quasi-$(2,\beta)$-Banach spaces

abstract: The main focus of this paper is to define the notion of quasi-$(2,\beta)$-Banach space and show some properties in this new space, by help of it and under some natural assumptions, we prove that the fixed point theorem [16, Theorem 2.1] is still valid in the setting of quasi-$(2,\beta)$-Banach spaces, this is also an extension of the fixed point result of Brzd\k{e}k et al. [12, Theorem 1] in $2$-Banach spaces to quasi-$(2,\beta)$-Banach spaces. In the next part, we give a general solution of the radical-type functional equation (1.2). In addition, we study the hyperstability results for these functional equation by applying the aforementioned fixed point theorem, and at the end of this paper we will derive some consequences.


title: Text-based Emotion Aware Recommender

abstract: We apply the concept of users' emotion vectors (UVECs) and movies' emotion vectors (MVECs) as building components of Emotion Aware Recommender System. We built a comparative platform that consists of five recommenders based on content-based and collaborative filtering algorithms. We employed a Tweets Affective Classifier to classify movies' emotion profiles through movie overviews. We construct MVECs from the movie emotion profiles. We track users' movie watching history to formulate UVECs by taking the average of all the MVECs from all the movies a user has watched. With the MVECs, we built an Emotion Aware Recommender as one of the comparative platforms' algorithms. We evaluated the top-N recommendation lists generated by these Recommenders and found the top-N list of Emotion Aware Recommender showed serendipity recommendations.


title: Ferri-chiral compounds with potentially switchable Dresselhaus spin split-ting

abstract: Spin splitting of energy bands can be induced by relativistic spin-orbit interactions in materials without inversion symmetry. Whereas polar space group symmetries permit Rashba (R-1) spin splitting with helical spin textures in momentum space, which could be reversed upon switching a ferroelectric polarization via applied electric fields, the ordinary Dresselhaus effect (D-1A) is ac-tive only in materials exhibiting nonpolar noncentrosymmetric crystal classes with atoms occupy-ing exclusively non-polar lattice sites. Consequently, the spin-momentum locking induced by D-1A is not electric field-switchable. An alternative type of Dresselhaus symmetry, referred to as D-1B, exhibits crystal class constraints similar to D-1A (all dipoles add up to zero), but unlike D-1A, at least one polar site is occupied. We find that this behavior exists in a class of ferri-chiral crys-tals, which in principle could be reversed upon a change in chirality. Focusing on alkali metal chalcogenides, we identify NaCu5S3 in the nonentiamorphic ferri-chiral structure, which exhibits CuS3 chiral units differing in the magnitude of their Cu displacements. We then synthesize NaCu5S3 (space group P6322) and confirm its ferri-chiral structure with power x-ray diffraction. Our electronic structure calculations demonstrate it exhibits D-1B spin splitting as well as a ferri-chiral phase transition, revealing spin splitting interdependent on chirality. Our electronic struc-ture calculations show that a few percent biaxial tensile strain can reduce (or nearly quench) the switching barrier separating the monodomain ferri-chiral P6322 states. We discuss what type of external stimuli might switch the chirality so as to reverse the (non-helical) Dresselhaus D-1B spin texture, offering an alternative to the traditional reversal of the (helical) Rashba spin texture.


title: A generalized stochastic control problem of bounded noise process under ambiguity arising in biological management

abstract: The objectives and contributions of this paper are mathematical and numerical analyses of a stochastic control problem of bounded population dynamics under ambiguity, an important but not well-studied problem, focusing on the optimality equation as a nonlinear degenerate parabolic partial integro-differential equation (PIDE). The ambiguity comes from lack of knowledge on the continuous and jump noises in the dynamics, and its optimization appears as nonlinear and nonlocal terms in the PIDE. Assuming a strong dynamic programming principle for continuous value functions, we characterize its solutions from both viscosity and distribution viewpoints. Numerical computation focusing on an ergodic case are presented as well to complement the mathematical analysis.


title: Confidence-Aware Learning for Deep Neural Networks

abstract: Despite the power of deep neural networks for a wide range of tasks, an overconfident prediction issue has limited their practical use in many safety-critical applications. Many recent works have been proposed to mitigate this issue, but most of them require either additional computational costs in training and/or inference phases or customized architectures to output confidence estimates separately. In this paper, we propose a method of training deep neural networks with a novel loss function, named Correctness Ranking Loss, which regularizes class probabilities explicitly to be better confidence estimates in terms of ordinal ranking according to confidence. The proposed method is easy to implement and can be applied to the existing architectures without any modification. Also, it has almost the same computational costs for training as conventional deep classifiers and outputs reliable predictions by a single inference. Extensive experimental results on classification benchmark datasets indicate that the proposed method helps networks to produce well-ranked confidence estimates. We also demonstrate that it is effective for the tasks closely related to confidence estimation, out-of-distribution detection and active learning.


title: A New Theoretical Framework of Pyramid Markov Processes for Blockchain Selfish Mining

abstract: In this paper, we provide a new theoretical framework of pyramid Markov processes to solve some open and fundamental problems of blockchain selfish mining under a rigorous mathematical setting. We first describe a more general model of blockchain selfish mining with both a two-block leading competitive criterion and a new economic incentive mechanism. Then we establish a pyramid Markov process and show that it is irreducible and positive recurrent, and its stationary probability vector is matrix-geometric with an explicitly representable rate matrix. Also, we use the stationary probability vector to study the influence of many orphan blocks on the waste of computing resource. Next, we set up a pyramid Markov reward process to investigate the long-run average profits of the honest and dishonest mining pools, respectively. As a by-product, we build three approximative Markov processes and provide some new interesting interpretation on the Markov chain and the revenue analysis reported in the seminal work by Eyal and Sirer (2014). Note that the pyramid Markov (reward) processes can open up a new avenue in the study of blockchain selfish mining. Thus we hope that the methodology and results developed in this paper shed light on the blockchain selfish mining such that a series of promising research can be developed potentially.


title: Least Squares Estimator for Vasicek Model Driven by Sub-fractional Brownian Processes from Discrete Observations

abstract: We study the parameter estimation problem of Vasicek Model driven by sub-fractional Brownian processes from discrete observations, and let {S_t^H,t>=0} denote a sub-fractional Brownian motion whose Hurst parameter 1/2<H<1 . The studies are as follows: firstly, two unknown parameters in the model are estimated by the least squares method. Secondly, the strong consistency and the asymptotic distribution of the estimators are studied respectively. Finally, our estimators are validated by numerical simulation.


title: Diffusion Limit of the Vlasov-Poisson-Boltzmann System

abstract: In the present paper, we study the diffusion limit of the classical solution to the unipolar Vlasov-Poisson-Boltzmann (VPB) system with initial data near a global Maxwellian. We prove the convergence and establish the convergence rate of the global strong solution to the unipolar VPB system towards the solution to an incompressible Navier-Stokes-Poisson-Fourier system based on the spectral analysis with precise estimation on the initial layer.


title: Rigorous Quantum Formulation of Parity-Time Symmetric Coupled Resonators

abstract: Rigorous quantum formulation of the Parity-Time (PT) symmetry phenomenon in the RF/microwave regime for a coupled coil resonators with lump elements has been presented. The coil resonator is described by the lump-element model that consists of an inductor (L), a resistor (R) and a capacitor (C). Rigorous quantum Hamiltonian for the coupled LRC coil resonators system has been derived through twice basis transforms of the original basis. The first basis transform rotates the original basis such that off-diagonal terms of the governing matrix of the equation system of the coupled coil resonators reduces to constants. Then a second basis transform obtains the quantum Hamiltonian, including the diagonal effective complex frequencies and the off-diagonal coupling terms, together with the transformed basis. With the obtain quantum Hamiltonian, the eigenvalues and eigenvectors of the coupled coil resonators can be obtained as usual as the quantum Hamiltonian. Finally, numerical simulation verifies the correctness of the theory. The quantum formulation of the coupled coil resonators can provide better guideline to design a better PT-symmetric system.


title: Flexibility in an asymmetric system with prolonged service time at non-dedicated servers

abstract: The prolonged service time at non-dedicated servers has been observed in [1]. Motivated by such real problems, we propose a stylized model which characterizes the feature of the prolonged service time at non-dedicated servers in an asymmetric system. We study the independent system, the full flexibility system and the partial flexibility system when the occupation rate of the system, the degree of the prolonged service time and the degree of the asymmetry are allowed to change. We show that under certain circumstances, the partial flexibility scheme outperforms the full flexibility system and the independent system in such a model. Our results also provide instructions on how to introduce flexibility when the service time at non-dedicated servers is prolonged in an asymmetric system.


title: Anatomy-Aware Siamese Network: Exploiting Semantic Asymmetry for Accurate Pelvic Fracture Detection in X-ray Images

abstract: Visual cues of enforcing bilaterally symmetric anatomies as normal findings are widely used in clinical practice to disambiguate subtle abnormalities from medical images. So far, inadequate research attention has been received on effectively emulating this practice in CAD methods. In this work, we exploit semantic anatomical symmetry or asymmetry analysis in a complex CAD scenario, i.e., anterior pelvic fracture detection in trauma PXRs, where semantically pathological (refer to as fracture) and non-pathological (e.g., pose) asymmetries both occur. Visually subtle yet pathologically critical fracture sites can be missed even by experienced clinicians, when limited diagnosis time is permitted in emergency care. We propose a novel fracture detection framework that builds upon a Siamese network enhanced with a spatial transformer layer to holistically analyze symmetric image features. Image features are spatially formatted to encode bilaterally symmetric anatomies. A new contrastive feature learning component in our Siamese network is designed to optimize the deep image features being more salient corresponding to the underlying semantic asymmetries (caused by pelvic fracture occurrences). Our proposed method have been extensively evaluated on 2,359 PXRs from unique patients (the largest study to-date), and report an area under ROC curve score of 0.9771. This is the highest among state-of-the-art fracture detection methods, with improved clinical indications.


title: Deep-PowerX: A Deep Learning-Based Framework for Low-Power Approximate Logic Synthesis

abstract: This paper aims at integrating three powerful techniques namely Deep Learning, Approximate Computing, and Low Power Design into a strategy to optimize logic at the synthesis level. We utilize advances in deep learning to guide an approximate logic synthesis engine to minimize the dynamic power consumption of a given digital CMOS circuit, subject to a predetermined error rate at the primary outputs. Our framework, Deep-PowerX, focuses on replacing or removing gates on a technology-mapped network and uses a Deep Neural Network (DNN) to predict error rates at primary outputs of the circuit when a specific part of the netlist is approximated. The primary goal of Deep-PowerX is to reduce the dynamic power whereas area reduction serves as a secondary objective. Using the said DNN, Deep-PowerX is able to reduce the exponential time complexity of standard approximate logic synthesis to linear time. Experiments are done on numerous open source benchmark circuits. Results show significant reduction in power and area by up to 1.47 times and 1.43 times compared to exact solutions and by up to 22% and 27% compared to state-of-the-art approximate logic synthesis tools while having orders of magnitudes lower run-time.


title: Task-agnostic Temporally Consistent Facial Video Editing

abstract: Recent research has witnessed the advances in facial image editing tasks. For video editing, however, previous methods either simply apply transformations frame by frame or utilize multiple frames in a concatenated or iterative fashion, which leads to noticeable visual flickers. In addition, these methods are confined to dealing with one specific task at a time without any extensibility. In this paper, we propose a task-agnostic temporally consistent facial video editing framework. Based on a 3D reconstruction model, our framework is designed to handle several editing tasks in a more unified and disentangled manner. The core design includes a dynamic training sample selection mechanism and a novel 3D temporal loss constraint that fully exploits both image and video datasets and enforces temporal consistency. Compared with the state-of-the-art facial image editing methods, our framework generates video portraits that are more photo-realistic and temporally smooth.


title: Quantum Pricing with a Smile: Implementation of Local Volatility Model on Quantum Computer

abstract: Applications of the quantum algorithm for Monte Carlo simulation to pricing of financial derivatives have been discussed in previous papers. However, up to now, the pricing model discussed in such papers is Black-Scholes model, which is important but simple. Therefore, it is motivating to consider how to implement more complex models used in practice in financial institutions. In this paper, we then consider the local volatility (LV) model, in which the volatility of the underlying asset price depends on the price and time. We present two types of implementation. One is the register-per-RN way, which is adopted in most of previous papers. In this way, each of random numbers (RNs) required to generate a path of the asset price is generated on a separated register, so the required qubit number increases in proportion to the number of RNs. The other is the PRN-on-a-register way, which is proposed in the author's previous work. In this way, a sequence of pseudo-random numbers (PRNs) generated on a register is used to generate paths of the asset price, so the required qubit number is reduced with a trade-off against circuit depth. We present circuit diagrams for these two implementations in detail and estimate required resources: qubit number and T-count.


title: Weak type $(p,p)$ bounds for Schr\"odinger groups via generalized Gaussian estimates

abstract: Let $L$ be a non-negative self-adjoint operator acting on $L^2(X)$, where $X$ is a space of homogeneous type with a dimension $n$. Suppose that the heat operator $e^{-tL}$ satisfies the generalized Gaussian $(p_0, p'_0)$-estimates of order $m$ for some $1\leq p_0 < 2$. It is known that the operator $(I+L)^{-s } e^{itL}$ is bounded on $L^p(X)$ for $s\geq n|{1/ 2}-{1/p}| $ and $ p\in (p_0, p_0')$ (see for example, \cite{Blunck2, BDN, CCO, CDLY, DN, Mi1}). In this paper we study the endpoint case $p=p_0$ and show that for $s_0= n\big|{1\over 2}-{1\over p_0}\big|$, the operator $(I+L)^{-{s_0}}e^{itL} $ is of weak type $(p_{0},p_{0})$, that is, there is a constant $C>0$, independent of $t$ and $f$ so that \begin{eqnarray*} \mu\left(\left\{x: \big|(I+L)^{-s_0}e^{itL} f(x)\big|>\alpha \right\} \right)\leq C (1+|t|)^{n(1 - {p_0\over 2}) } \left( {\|f\|_{p_0} \over \alpha} \right)^{p_0} , \ \ \ t\in{\mathbb R} \end{eqnarray*} for $\alpha>0$ when $\mu(X)=\infty$, and $\alpha>\big(\|f\|_{p_{0}}/\mu(X) \big)^{p_{0}}$ when $\mu(X)<\infty$. Our results can be applied to Schr\"odinger operators with rough potentials and %second order elliptic operators with rough lower order terms, or higher order elliptic operators with bounded measurable coefficients although in general, their semigroups fail to satisfy Gaussian upper bounds.


title: The Schr\"odinger equation in $L^p$ spaces for operators with heat kernel satisfying Poisson type bounds

abstract: Let $L$ be a non-negative self-adjoint operator acting on $L^2(X)$ where $X$ is a space of homogeneous type with a dimension $n$. In this paper, we study sharp endpoint $L^p$-Sobolev estimates for the solution of the initial value problem for the Schr\"odinger equation, $i \partial_t u + L u=0 $ and show that for all $f\in L^p(X), 1<p<\infty,$ \begin{eqnarray*} \left\| e^{itL} (I+L)^{-{\sigma n}} f\right\|_{p} \leq C(1+|t|)^{\sigma n} \|f\|_{p}, \ \ \ t\in{\mathbb R}, \ \ \ \sigma\geq \big|{1\over 2}-{1\over p}\big|, \end{eqnarray*} where the semigroup $e^{-tL}$ generated by $L$ satisfies a Poisson type upper bound. This extends the previous result in \cite{CDLY1} in which the semigroup $e^{-tL}$ generated by $L$ satisfies the exponential decay.


title: Operational, gauge-free quantum tomography

abstract: As increasingly impressive quantum information processors are realized in laboratories around the world, robust and reliable characterization of these devices is now more urgent than ever. These diagnostics can take many forms, but one of the most popular categories is tomography, where an underlying parameterized model is proposed for a device and inferred by experiments. Here, we introduce and implement efficient operational tomography, which uses experimental observables as these model parameters. This addresses a problem of ambiguity in representation that arises in current tomographic approaches (the gauge problem). Solving the gauge problem enables us to efficiently implement operational tomography in a Bayesian framework computationally, and hence gives us a natural way to include prior information and discuss uncertainty in fit parameters. We demonstrate this new tomography in a variety of different experimentally-relevant scenarios, including standard process tomography, Ramsey interferometry, randomized benchmarking, and gate set tomography.


title: An adaptive multiresolution ultra-weak discontinuous Galerkin method for nonlinear Schrodinger equations

abstract: This paper develops a high order adaptive scheme for solving nonlinear Schrodinger equations. The solutions to such equations often exhibit solitary wave and local structures, which makes adaptivity essential in improving the simulation efficiency. Our scheme uses the ultra-weak discontinuous Galerkin (DG) formulation and belongs to the framework of adaptive multiresolution schemes. Various numerical experiments are presented to demonstrate the excellent capability of capturing the soliton waves and the blow-up phenomenon.


title: Increasing Trustworthiness of Deep Neural Networks via Accuracy Monitoring

abstract: Inference accuracy of deep neural networks (DNNs) is a crucial performance metric, but can vary greatly in practice subject to actual test datasets and is typically unknown due to the lack of ground truth labels. This has raised significant concerns with trustworthiness of DNNs, especially in safety-critical applications. In this paper, we address trustworthiness of DNNs by using post-hoc processing to monitor the true inference accuracy on a user's dataset. Concretely, we propose a neural network-based accuracy monitor model, which only takes the deployed DNN's softmax probability output as its input and directly predicts if the DNN's prediction result is correct or not, thus leading to an estimate of the true inference accuracy. The accuracy monitor model can be pre-trained on a dataset relevant to the target application of interest, and only needs to actively label a small portion (1% in our experiments) of the user's dataset for model transfer. For estimation robustness, we further employ an ensemble of monitor models based on the Monte-Carlo dropout method. We evaluate our approach on different deployed DNN models for image classification and traffic sign detection over multiple datasets (including adversarial samples). The result shows that our accuracy monitor model provides a close-to-true accuracy estimation and outperforms the existing baseline methods.


title: Fujita decomposition and Massey product for fibered varieties

abstract: Let $f\colon X\to B$ be a semistable fibration where $X$ is a smooth variety of dimension $n\geq 2$ and $B$ is a smooth curve. We give the structure theorem for the local system of the relative $1$-forms and of the relative top forms. This gives a neat interpretation of the second Fujita decomposition of $f_*\omega_{X/B}$. We apply our interpretation to show the existence, up to base change, of higher irrational pencils and on the finiteness of the associated monodromy representations under natural Castelnuovo-type hypothesis on local subsystems. Finally we give a criterion to have that $X$ is not of Albanese general type if $B=\mathbb{P}^1$.