title: Optimal Collusion-Free Teaching

abstract: Formal models of learning from teachers need to respect certain criteria to avoid collusion. The most commonly accepted notion of collusion-freeness was proposed by Goldman and Mathias (1996), and various teaching models obeying their criterion have been studied. For each model $M$ and each concept class $\mathcal{C}$, a parameter $M$-$\mathrm{TD}(\mathcal{C})$ refers to the teaching dimension of concept class $\mathcal{C}$ in model $M$---defined to be the number of examples required for teaching a concept, in the worst case over all concepts in $\mathcal{C}$. This paper introduces a new model of teaching, called no-clash teaching, together with the corresponding parameter $\mathrm{NCTD}(\mathcal{C})$. No-clash teaching is provably optimal in the strong sense that, given any concept class $\mathcal{C}$ and any model $M$ obeying Goldman and Mathias's collusion-freeness criterion, one obtains $\mathrm{NCTD}(\mathcal{C})\le M$-$\mathrm{TD}(\mathcal{C})$. We also study a corresponding notion $\mathrm{NCTD}^+$ for the case of learning from positive data only, establish useful bounds on $\mathrm{NCTD}$ and $\mathrm{NCTD}^+$, and discuss relations of these parameters to the VC-dimension and to sample compression. In addition to formulating an optimal model of collusion-free teaching, our main results are on the computational complexity of deciding whether $\mathrm{NCTD}^+(\mathcal{C})=k$ (or $\mathrm{NCTD}(\mathcal{C})=k$) for given $\mathcal{C}$ and $k$. We show some such decision problems to be equivalent to the existence question for certain constrained matchings in bipartite graphs. Our NP-hardness results for the latter are of independent interest in the study of constrained graph matchings.


title: Optimal Path Planning for Connected and Automated Vehicles at Urban Intersections

abstract: In earlier work, a decentralized optimal control framework was established for coordinating online connected and automated vehicles (CAVs) at urban intersections. The policy designating the sequence that each CAV crosses the intersection, however, was based on a first-in-first-out queue, imposing limitations on the optimal solution. Moreover, no lane changing, or left and right turns were considered. In this paper, we formulate an upper-level optimization problem, the solution of which yields, for each CAV, the optimal sequence and lane to cross the intersection. The effectiveness of the proposed approach is illustrated through simulation.


title: High-temperature domain walls of QCD with imaginary chemical potentials

abstract: We study QCD with massless quarks on $\mathbb{R}^3\times S^1$ under symmetry-twisted boundary conditions with small compactification radius, i.e. at high temperatures. Under suitable boundary conditions, the theory acquires a part of the center symmetry and it is spontaneously broken at high temperatures. We show that these vacua at high temperatures can be regarded as different symmetry-protected topological orders, and the domain walls between them support nontrivial massless gauge theories as a consequence of anomaly-inflow mechanism. At sufficiently high temperatures, we can perform the semiclassical analysis to obtain the domain-wall theory, and $2$d $U(N_\mathrm{c}-1)$ gauge theories with massless fermions match the 't~Hooft anomaly. We perform these analysis for the high-temperature domain wall of $\mathbb{Z}_{N_\mathrm{c}}$-QCD and also of Roberge-Weiss phase transitions.


title: NormalNet: Learning-based Normal Filtering for Mesh Denoising

abstract: Mesh denoising is a critical technology in geometry processing that aims to recover high-fidelity 3D mesh models of objects from their noise-corrupted versions. In this work, we propose a learning-based normal filtering scheme for mesh denoising called NormalNet, which maps the guided normal filtering (GNF) into a deep network. The scheme follows the iterative framework of filtering-based mesh denoising. During each iteration, first, the voxelization strategy is applied on each face in a mesh to transform the irregular local structure into the regular volumetric representation, therefore, both the structure and face normal information are preserved and the convolution operations in CNN(Convolutional Neural Network) can be easily performed. Second, instead of the guidance normal generation and the guided filtering in GNF, a deep CNN is designed, which takes the volumetric representation as input, and outputs the learned filtered normals. At last, the vertex positions are updated according to the filtered normals. Specifically, the iterative training framework is proposed, in which the generation of training data and the network training are alternately performed, whereas the ground truth normals are taken as the guidance normals in GNF to get the target normals. Compared to state-of-the-art works, NormalNet can effectively remove noise while preserving the original features and avoiding pseudo-features.


title: $\beta^3$-IRT: A New Item Response Model and its Applications

abstract: Item Response Theory (IRT) aims to assess latent abilities of respondents based on the correctness of their answers in aptitude test items with different difficulty levels. In this paper, we propose the $\beta^3$-IRT model, which models continuous responses and can generate a much enriched family of Item Characteristic Curve (ICC). In experiments we applied the proposed model to data from an online exam platform, and show our model outperforms a more standard 2PL-ND model on all datasets. Furthermore, we show how to apply $\beta^3$-IRT to assess the ability of machine learning classifiers. This novel application results in a new metric for evaluating the quality of the classifier's probability estimates, based on the inferred difficulty and discrimination of data instances.


title: A Superconvergent Ensemble HDG Method for Parameterized Convection Diffusion Equations

abstract: In this paper, we first devise an ensemble hybridizable discontinuous Galerkin (HDG) method to efficiently simulate a group of parameterized convection diffusion PDEs. These PDEs have different coefficients, initial conditions, source terms and boundary conditions. The ensemble HDG discrete system shares a common coefficient matrix with multiple right hand side (RHS) vectors; it reduces both computational cost and storage. We have two contributions in this paper. First, we derive an optimal $L^2$ convergence rate for the ensemble solutions on a general polygonal domain, which is the first such result in the literature. Second, we obtain a superconvergent rate for the ensemble solutions after an element-by-element postprocessing under some assumptions on the domain and the coefficients of the PDEs. We present numerical experiments to confirm our theoretical results.


title: Limit theorems for some time dependent expanding dynamical systems

abstract: In this paper we will prove various probabilistic limit theorems for some classes of distance expanding sequential dynamical systems (SDS). Our starting point here is certain sequential complex Ruelle-Perron-Frobenius (RPF) theorems which were proved in \cite{book} and \cite{ASIP me} using contraction properties of a complex version of the projective Hilbert metric developed in \cite{Rug}. We will start from the growth rate of the variances of the underlying partial sums. This is well understood in the random dynamics setup, when the maps are stationary, but not in the SDS setup, where various growth rates can occur. Some of our results in this direction rely on certain type of stability in these RPF theorems, which is one of the novelties of this paper. Then we will provide general conditions for several classical limit theorems to hold true in the sequential setup. Some of our general results mostly have applications for composition of random non-stationary map, while the conditions of the other results hold true for general type of SDS. In the latter setup, results such as the Berry-Esse\'en theorem and the local central limit theorem were not obtained so far even for independent but not identically distributed maps, which is a particular case of the setup considered in this paper.


title: Deep Reinforcement Learning of Volume-guided Progressive View Inpainting for 3D Point Scene Completion from a Single Depth Image

abstract: We present a deep reinforcement learning method of progressive view inpainting for 3D point scene completion under volume guidance, achieving high-quality scene reconstruction from only a single depth image with severe occlusion. Our approach is end-to-end, consisting of three modules: 3D scene volume reconstruction, 2D depth map inpainting, and multi-view selection for completion. Given a single depth image, our method first goes through the 3D volume branch to obtain a volumetric scene reconstruction as a guide to the next view inpainting step, which attempts to make up the missing information; the third step involves projecting the volume under the same view of the input, concatenating them to complete the current view depth, and integrating all depth into the point cloud. Since the occluded areas are unavailable, we resort to a deep Q-Network to glance around and pick the next best view for large hole completion progressively until a scene is adequately reconstructed while guaranteeing validity. All steps are learned jointly to achieve robust and consistent results. We perform qualitative and quantitative evaluations with extensive experiments on the SUNCG data, obtaining better results than the state of the art.


title: Simulated Refraction-Enhanced X-Ray Radiography of Laser-Driven Shocks

abstract: Refraction-enhanced x-ray radiography (REXR) is used to infer shock-wave positions of more than one shock wave, launched by a multiple-picket pulse in a planar plastic foil. This includes locating shock waves before the shocks merge, during the early time and the main drive of the laser pulse that is not possible with the velocity interferometer system for any reflector. Simulations presented in this paper of REXR show that it is necessary to incorporate refraction and attenuation of x rays along with the appropriate opacity and refractive-index tables to interpret experimental images. Simulated REXR shows good agreement with an experiment done on the OMEGA laser facility to image a shock wave. REXR can be applied to design multiple-picket pulses with a better understanding of the shock locations. This will be beneficial to obtain the required adiabats for inertial confinement fusion implosions.


title: Quantitative unique continuation for Schr\"odinger operators

abstract: We investigate the quantitative unique continuation properties of solutions to second order elliptic equations with singular lower order terms. The main theorem presents a quantification of the strong unique continuation property for $\Delta + V$. That is, for any non-trivial $u$ that solves $\Delta u + V u = 0$ in some open, connected subset of $\mathbb{R}^n$, we estimate the vanishing order of solutions in terms of the $L^t$-norm of $V$. Our results apply to all $t > \frac n 2$ and $n \ge 3$. With these maximal order of vanishing estimates, we employ a scaling argument to produce quantitative unique continuation at infinity estimates for global solutions to $\Delta u + V u = 0$. To handle $V \in L^t$ for every $t \in (\frac n 2, \infty]$, we prove a novel $L^p - L^q$ Carleman estimate by interpolating a known $L^p - L^2$ estimate with a new endpoint Carleman estimate. This new Carleman estimate may also be used to establish improved order of vanishing estimates for equations with a first order term, those of the form $\Delta u + W \cdot \nabla u + V u = 0$.


title: Coherent long-distance spin-qubit-transmon coupling

abstract: Spin qubits and superconducting qubits are among the promising candidates for a solid state quantum computer. For the implementation of a hybrid architecture which can profit from the advantages of either world, a coherent long-distance link is necessary that integrates and couples both qubit types on the same chip. We realize such a link with a frequency-tunable high impedance SQUID array resonator. The spin qubit is a resonant exchange qubit hosted in a GaAs triple quantum dot. It can be operated at zero magnetic field, allowing it to coexist with superconducting qubits on the same chip. We find a working point for the spin qubit, where the ratio between its coupling strength and decoherence rate is optimized. We observe coherent interaction between the resonant exchange qubit and a transmon qubit in both resonant and dispersive regimes, where the interaction is mediated either by real or virtual resonator photons.


title: Constitutive law of nonlocal gravity

abstract: We analyze the structure of a recent nonlocal generalization of Einstein's theory of gravitation by Mashhoon et al. By means of a covariant technique, we derive an expanded version of the nonlocality tensor which constitutes the theory. At the lowest orders of approximation, this leads to a simplification which sheds light on the fundamental structure of the theory and may prove useful in the search for exact solutions of nonlocal gravity.


title: Semiclassical probability of radiation of twisted photons in the ultrarelativistic limit

abstract: The semiclassical general formula for the probability of radiation of twisted photons by ultrarelativistic scalar and Dirac particles moving in the electromagnetic field of a general form is derived. This formula is the analog of the Baier-Katkov formula for the probability of radiation of one plane wave photon with the quantum recoil taken into account. The derived formula is used to describe the radiation of twisted photons by charged particles in undulators and laser waves. Thus, the general theory of undulator radiation of twisted photons and radiation of twisted photons in the nonlinear Compton process is developed with account for the quantum recoil. The explicit formulas for the probability to record a twisted photon are obtained in these cases. In particular, we found that the quantum recoil and spin degrees of freedom increase the radiation probability of twisted photons in comparison with the formula for scalar particles without recoil. In the range of applicability of the semiclassical formula, the selection rules for undulator radiation established in the purely classical framework are not violated. The manifestation of the blossoming out rose effect in the nonlinear Compton process in a strong laser wave with circular polarization and in the wiggler radiation is revealed. Several examples are studied: the radiation of MeV twisted photons by $180$ GeV electrons in the wiggler; the radiation of twisted photons by $256$ MeV electrons in strong electromagnetic waves produced by the CO$_2$ and Ti:Sa lasers; and the radiation of MeV twisted photons by $51.1$ MeV electrons in the electromagnetic wave generated by the FEL with photon energy $1$ keV.


title: Group-wise Correlation Stereo Network

abstract: Stereo matching estimates the disparity between a rectified image pair, which is of great importance to depth sensing, autonomous driving, and other related tasks. Previous works built cost volumes with cross-correlation or concatenation of left and right features across all disparity levels, and then a 2D or 3D convolutional neural network is utilized to regress the disparity maps. In this paper, we propose to construct the cost volume by group-wise correlation. The left features and the right features are divided into groups along the channel dimension, and correlation maps are computed among each group to obtain multiple matching cost proposals, which are then packed into a cost volume. Group-wise correlation provides efficient representations for measuring feature similarities and will not lose too much information like full correlation. It also preserves better performance when reducing parameters compared with previous methods. The 3D stacked hourglass network proposed in previous works is improved to boost the performance and decrease the inference computational cost. Experiment results show that our method outperforms previous methods on Scene Flow, KITTI 2012, and KITTI 2015 datasets. The code is available at https://github.com/xy-guo/GwcNet


title: A note on a recent attempt to prove Sendov's conjecture

abstract: Recently GM Sofi & SA Shabir [arXive: 1903.01850v2 [math.GM] 6 Mar 2019] made an attempt to prove the Sendov's conjecture. But unfortunately the proof is not correct. In this note, we discuss the fallacy in the proof.


title: Just-Enough Interaction Approach to Knee MRI Segmentation: Data from the Osteoarthritis Initiative

abstract: State-of-the-art automated segmentation algorithms are not 100\% accurate especially when segmenting difficult to interpret datasets like those with severe osteoarthritis (OA). We present a novel interactive method called just-enough interaction (JEI), which adds a fast correction step to the automated layered optimal graph segmentation of multiple objects and surfaces (LOGISMOS). After LOGISMOS segmentation in knee MRI, the JEI user interaction does not modify boundary surfaces of the bones and cartilages directly. Local costs of underlying graph nodes are modified instead and the graph is re-optimized, providing globally optimal corrected results. Significant performance improvement ($p \ll 0.001$) was observed when comparing JEI-corrected results to the automated. The algorithm was extended from 3D JEI to longitudinal multi-3D (4D) JEI allowing simultaneous visualization and interaction of multiple-time points of the same patient.


title: Tools for designing atom interferometers in a microgravity environment

abstract: We present a variational model suitable for rapid preliminary design of atom interferometers in a microgravity environment. The model approximates the solution of the 3D rotating--frame Gross--Pitaevskii equation (GPE) as the sum of Nc Gaussian clouds. Each Gaussian cloud is assumed to have time--dependent center positions, widths, and linear and quadratic phase parameters. We applied the Lagrangian Variational Method (LVM) with this trial wave function to derive equations of motion for these parameters that can be adapted to any external potential. We also present a 1D version of this variational model. As an example we apply the model to a 1D atom interferometry scheme for measuring Newton's gravitational constant, G, in a microgravity environment. We show how the LVM model can (1) constrain the experimental parameter space size, (2) show how the value of G can be obtained from the experimental conditions and interference pattern characteristics, and (3) show how to improve the sensitivity of the measurement and construct a preliminary error budget.


title: Continuous Data Assimilation Reduced Order Models of Fluid Flow

abstract: We propose, analyze, and test a novel continuous data assimilation reduced order model (DA-ROM) for simulating incompressible flows. While ROMs have a long history of success on certain problems with recurring dominant structures, they tend to lose accuracy on more complicated problems and over longer time intervals. Meanwhile, continuous data assimilation (DA) has recently been used to improve accuracy and, in particular, long time accuracy in fluid simulations by incorporating measurement data into the simulation. This paper synthesizes these two ideas, in an attempt to address inaccuracies in ROM by applying DA, especially over long time intervals and when only inaccurate snapshots are available. We prove that with a properly chosen nudging parameter, the proposed DA-ROM algorithm converges exponentially fast in time to the true solution, up to discretization and ROM truncation errors. Finally, we propose a strategy for nudging adaptively in time, by adjusting dissipation arising from the nudging term to better match true solution energy. Numerical tests confirm all results, and show that the DA-ROM strategy with adaptive nudging can be highly effective at providing long time accuracy in ROMs.


title: High- and Low-$\alpha$ Disk Stars Separate Dynamically at all Ages

abstract: There is a dichotomy in the Milky Way in the $[\alpha/$Fe]-[Fe/H] plane, in which stars fall into high-$\alpha$, and low-$\alpha$ sequences. The high-$\alpha$ sequence comprises mostly old stars, and the low-$\alpha$ sequence comprises primarily young stars. The origin of this dichotomy is uncertain. To better understand how the high- and low-$\alpha$ stars are affiliated, we examine if the high- and low-$\alpha$ sequences have distinct orbits at all ages, or if age sets the orbital properties of stars irrespective of their $\alpha$-enhancement. Orbital actions $J_R$, $J_z$, and $J_\phi$ (or $L_z$) are our labels of stellar dynamics. We use ages for 58,278 LAMOST stars (measured to a precision of 40\%) within $\leq$2kpc of the Sun and we calculate orbital actions from proper motions and parallaxes given by Gaia's DR2. We find that \emph{at all ages}, the high- and low-$\alpha$ sequences are dynamically distinct. This implies separate formation and evolutionary histories for the two sequences; a star's membership in the high- or low-$\alpha$ sequence indicates its dynamical properties at a given time. We use action space to make an efficient selection of halo stars and subsequently report a group of old, low-$\alpha$ stars in the halo, which may be a discrete population from an infall event.


title: High-frequency EPR study on Cu$_4$Cu- and Co$_4$Co-metallacrown complexes

abstract: High-frequency/high-field electron paramagnetic resonance studies on two homo-nuclear 12-MC-4 metallacrown complexes Cu$_4$Cu and Co$_4$Co are presented. For Cu$_4$Cu, our data imply axial-type $g$-anisotropy with $g_{x}$ = 2.03 $\pm$ 0.01, $g_{y}$ = 2.04 $\pm$ 0.01, and $g_{z}$ = 2.23 $\pm$ 0.01, yielding $g=2.10 \pm 0.02$. No significant zero field splitting (ZFS) of the ground state mode is observed. In Co$_4$Co, we find a m$_S$ = $\pm$3/2 ground state with $g=2.66$. The data suggest large anisotropy $D$ of negative sign.


title: Electron-impact fine-structure excitation of Fe II at low temperature

abstract: Fe II emission lines are observed from nearly all classes of astronomical objects over a wide spectral range from the infrared to the ultraviolet. To meaningfully interpret these lines, reliable atomic data are necessary. In work presented here we focused on low-lying fine-structure transitions, within the ground term, due to electron impact. We provide effective collision strengths together with estimated uncertainties as functions of temperature of astrophysical importance ($10 - 100,000$ K). Due to the importance of fine-structure transitions within the ground term, the focus of this work is on obtaining accurate rate coefficients at the lower end of this temperature range, for applications in low temperature environments such as the interstellar medium. We performed three different flavours of scattering calculations: i) a intermediate coupling frame transformation (ICFT) $R$-matrix method, ii) a Breit-Pauli (BP) $R$-matrix method, and iii) a Dirac $R$-matrix method. The ICFT and BP $R$-matrix calculations involved three different AUTOSTRUCTURE target models each. The Dirac $R$-matrix calculation was based on a reliable 20 configuration, 6069 level atomic structure model. Good agreement was found with our BP and Dirac $R$-matrix collision results compared to previous $R$-matrix calculations. We present a set of recommended effective collision strengths for the low-lying forbidden transitions together with associated uncertainty estimates.


title: Oriented Local Moves and Divisibility of the Jones Polynomial

abstract: For any virtual link $L = S \cup T$ that may be decomposed into a pair of oriented $n$-tangles $S$ and $T$, an oriented local move of type $T \mapsto T'$ is a replacement of $T$ with the $n$-tangle $T'$ in a way that preserves the orientation of $L$. After developing a general decomposition for the Jones polynomial of the virtual link $L = S \cup T$ in terms of various (modified) closures of $T$, we analyze the Jones polynomials of virtual links $L_1,L_2$ that differ via a local move of type $T \mapsto T'$. Succinct divisibility conditions on $V(L_1)-V(L_2)$ are derived for broad classes of local moves that include the $\Delta$-move and the double-$\Delta$-move as special cases. As a consequence of our divisibility result for the double-$\Delta$-move, we introduce a necessary condition for any pair of classical knots to be $S$-equivalent.


title: Reversible Quaternionic Hyperbolic Isometries

abstract: Let $G$ be a group. An element $g$ in $G$ is called reversible if it is conjugate to $g^{-1}$ within $G$, and called strongly reversible if it is conjugate to its inverse by an order two element of $G$. Let $\textbf{H}_{\mathbb H}^n$ be the $n$-dimensional quaternionic hyperbolic space. Let $\mathrm{PSp}(n,1)$ be the isometry group of $\textbf{H}_{\mathbb H}^n$. In this paper, we classify reversible and strongly reversible elements in $\mathrm{Sp}(n)$ and $\mathrm{Sp}(n,1)$. Also, we prove that all the elements of $\mathrm{PSp}(n,1)$ are strongly reversible.


title: Retailer response to wholesale stockouts

abstract: The purpose of this paper is to identify the immediate and future retailer response to wholesale stockouts. We perform a statistical analysis of historical customer order and delivery data of a local tool wholesaler and distributor, whose customers are retailers, over a period of four years. We investigate the effect of customer service on the order fill rate and the rate of future demand, where the customer service is defined in terms of timely delivery and the fill rate is defined as the fraction of the order that is eventually materialized, i.e., is not cancelled following a stockout.We find that for customers who order frequently, stockouts have an adverse effect on the fill rate of their orders and on the frequency but not the value of their future demand; however, this latter effect seems to be more short-term than long-term. Practically all studies on the effects of stockouts measure immediate reported/intended consumer purchase incidence and choice decision behavior in response to stockouts in retail environments, mostly based on surveys. This study looks at how stockouts affect future demand in a wholesale environment, based on historical behavioral data analysis.


title: Reminiscences

abstract: I would like to thank Solar Physics colleagues for asking me to write this chapter on my professional life. My main interest has always been focused on the Sun, our star, from the heating of the corona, to the dynamics of prominences and their eruptions, ares and coronal mass ejections until their impact on the Earth. I built a new group in solar physics and gave to them my enthusiasm. They brought to me a lot of satisfaction. We have made important advances in solar physics with a step forward to understand the triggers of solar activity and their terrestrial effects. Our avant-garde research and discovery has opened new topics for the solar community. Mixing observations obtained on the ground and in space with theory and numerical simulations brings a new perspective in research.


title: Self-Localized Solutions of the Nonlinear Quantum Harmonic Oscillator

abstract: We analyze the existences, properties, and stabilities of the self-localized solutions of the nonlinear quantum harmonic oscillator (NQHO) using spectral renormalization method (SRM). We show that self-localized single, dual and triple soliton solutions of the NQHO do exists, however, only single and dual soliton solutions satisfy the necessary Vakhitov and Kolokolov slope condition, therefore triple soliton solution is found to be unstable, at least for the parameter ranges considered. Additionally, we investigate the stability characteristics of the single and dual soliton solutions using a split-step Fourier scheme. We show that single and dual soliton solutions are pulsating during time stepping. We discuss our findings and comment on our results.


title: Crossover of Charge Fluctuations across the Strange Metal Phase Diagram

abstract: A normal metal exhibits a valence plasmon, which is a sound wave in its conduction electron density. The mysterious strange metal is characterized by non-Boltzmann transport and violates most fundamental Fermi liquid scaling laws. A fundamental question is: Do strange metals have plasmons? Using momentum-resolved inelastic electron scattering (M-EELS) we recently showed that, rather than a plasmon, optimally-doped Bi$_{2.1}$Sr$_{1.9}$Ca$_{1.0}$Cu$_{2.0}$O$_{8+x}$ (Bi-2212) exhibits a featureless, temperature-independent continuum with a power-law form over most energy and momentum scales [M. Mitrano, PNAS 115, 5392-5396 (2018)]. Here, we show that this continuum is present throughout the fan-shaped, strange metal region of the phase diagram. Outside this region, dramatic changes in spectral weight are observed: In underdoped samples, spectral weight up to 0.5 eV is enhanced at low temperature, biasing the system towards a charge order instability. The situation is reversed in the overdoped case, where spectral weight is strongly suppressed at low temperature, increasing quasiparticle coherence in this regime. Optimal doping corresponds to the boundary between these two opposite behaviors at which the response is temperature-independent. Our study suggests that plasmons do not exist as well-defined excitations in Bi-2212, and that a featureless continuum is a defining property of the strange metal, which is connected to a peculiar crossover where the spectral weight change undergoes a sign reversal.


title: Knowledge compilation languages as proof systems

abstract: In this paper, we study proof systems in the sense of Cook-Reckhow for problems that are higher in the polynomial hierarchy than coNP, in particular, #SAT and maxSAT. We start by explaining how the notion of Cook-Reckhow proof systems can be apply to these problems and show how one can twist existing languages in knowledge compilation such as decision DNNF so that they can be seen as proof systems for problems such as #SAT and maxSAT.


title: Measurement of electroweak production of a W boson in association with two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A measurement is presented of electroweak (EW) production of a W boson in association with two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV. The data sample was recorded by the CMS Collaboration at the LHC and corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The measurement is performed for the $\ell\nu$jj final state (with $\ell\nu$ indicating a lepton-neutrino pair, and j representing the quarks produced in the hard interaction) in a kinematic region defined by invariant mass $m_\mathrm{jj}$ $>$ 120 GeV and transverse momenta $p_\mathrm{T j}$ $>$ 25 GeV. The cross section of the process is measured in the electron and muon channels yielding $\sigma_\mathrm{EW}$(Wjj) = 6.23 $\pm$ 0.12 (stat) $\pm$ 0.61 (syst) pb per channel, in agreement with leading-order standard model predictions. The additional hadronic activity of events in a signal-enriched region is studied, and the measurements are compared with predictions. The final state is also used to perform a search for anomalous trilinear gauge couplings. Limits on anomalous trilinear gauge couplings associated with dimension-six operators are given in the framework of an effective field theory. The corresponding 95% confidence level intervals are $-$2.3 $<$ $c_{\mathrm{WWW}}/\Lambda^2$ $<$ 2.5 TeV$^{-2}$, $-$8.8 $<$ $c_{\mathrm{W}}/\Lambda^2$ $<$ 16 TeV$^{-2}$, and $-$45 $<$ $c_{\mathrm{B}}/\Lambda^2$ $<$ 46 TeV$^{-2}$. These results are combined with the CMS EW Zjj analysis, yielding the constraint on the $c_{\mathrm{WWW}}$ coupling: $-$1.8 $<$ $c_{\mathrm{WWW}}/\Lambda^2$ $<$ 2.0 TeV$^{-2}$.


abstract: The decisive role of magnetic anisotropy even in systems with small anisotropy is illustrated for the honeycomb-layered antiferromagnets The decisive role of magnetic anisotropy in honeycomb layered Li$_3$Ni$_2$SbO$_6$ and Na$_3$Ni$_2$SbO$_6$ with $A$ = Li and Na. Both systems evolve long range magnetic order below $T_{\rm N}$ = 14 and 16.5~K, respectively. The magnetic phase diagrams obtained from static magnetisation studies up to 15~T imply competing antiferromagnetic phases and a tricritical point at $T_{\rm N}$. The phase boundaries are visible in the dynamic response of the antiferromagnetic resonance modes, too, which investigation by means of high frequency/high field electron spin resonance enables precise determination of magnetic anisotropy. The anisotropy gap amounts to $\Delta = 360 \pm 2$~GHz in Na$_3$Ni$_2$SbO$_6$ while in Li$_3$Ni$_2$SbO$_6$ orthorhombicity is associated with $\Delta = 198 \pm 4$ and $218 \pm 4$~GHz. Above $T_{\rm N}$, the data imply short-range antiferromagnetic order up to at least 80~K. The data suggest a crucial role of anisotropy for selecting the actual spin structure at $B=0$~T.


title: Algorithms for an Efficient Tensor Biclustering

abstract: Consider a data set collected by (individuals-features) pairs in different times. It can be represented as a tensor of three dimensions (Individuals, features and times). The tensor biclustering problem computes a subset of individuals and a subset of features whose signal trajectories over time lie in a low-dimensional subspace, modeling similarity among the signal trajectories while allowing different scalings across different individuals or different features. This approach are based on spectral decomposition in order to build the desired biclusters. We evaluate the quality of the results from each algorithms with both synthetic and real data set.


title: Streamlined Variational Inference for Higher Level Group-Specific Curve Models

abstract: A two-level group-specific curve model is such that the mean response of each member of a group is a separate smooth function of a predictor of interest. The three-level extension is such that one grouping variable is nested within another one, and higher level extensions are analogous. Streamlined variational inference for higher level group-specific curve models is a challenging problem. We confront it by systematically working through two-level and then three-level cases and making use of the higher level sparse matrix infrastructure laid down in Nolan and Wand (2018). A motivation is analysis of data from ultrasound technology for which three-level group-specific curve models are appropriate. Whilst extension to the number of levels exceeding three is not covered explicitly, the pattern established by our systematic approach sheds light on what is required for even higher level group-specific curve models.


title: Fluid statics of a self-gravitating perfect-gas isothermal sphere

abstract: We open the paper with introductory considerations describing the motivations of our long-term research plan targeting gravitomagnetism, illustrating the fluid-dynamics numerical test case selected for that purpose, that is, a perfect-gas sphere contained in a solid shell located in empty space sufficiently away from other masses, and defining the main objective of this study: the determination of the gravitofluid-static field required as initial field ($t=0$) in forthcoming fluid-dynamics calculations. The determination of the gravitofluid-static field requires the solution of the isothermal-sphere Lane-Emden equation. We do not follow the habitual approach of the literature based on the prescription of the central density as boundary condition; we impose the gravitational field at the solid-shell internal wall. As the discourse develops, we point out differences and similarities between the literature's and our approach. We show that the nondimensional formulation of the problem hinges on a unique physical characteristic number that we call gravitational number because it gauges the self-gravity effects on the gas' fluid statics. We illustrate and discuss numerical results; some peculiarities, such as gravitational-number upper bound and multiple solutions, lead us to investigate the thermodynamics of the physical system, particularly entropy and energy, and preliminarily explore whether or not thermodynamic-stability reasons could provide justification for either selection or exclusion of multiple solutions. We close the paper with a summary of the present study in which we draw conclusions and describe future work.


title: Exceptional points of two-dimensional random walks at multiples of the cover time

abstract: We study exceptional sets of the local time of the continuous-time simple random walk in scaled-up (by $N$) versions $D_N\subseteq \mathbb Z^2$ of bounded open domains $D\subseteq \mathbb R^2$. Upon exit from $D_N$, the walk lands on a "boundary vertex" and then reenters $D_N$ through a random boundary edge in the next step. In the parametrization by the local time at the "boundary vertex" we prove that, at times corresponding to a $\theta$-multiple of the cover time of $D_N$, the sets of suitably defined $\lambda$-thick (i.e., heavily visited) and $\lambda$-thin (i.e., lightly visited) points are, as $N\to\infty$, distributed according to the Liouville Quantum Gravity $Z^D_\lambda$ with parameter $\lambda$-times the critical value. For $\theta<1$, also the set of avoided vertices (a.k.a. late points) and the set where the local time is of order unity are distributed according to $Z^D_{\sqrt\theta}$. The local structure of the exceptional sets is described as well, and is that of a pinned Discrete Gaussian Free Field for the thick and thin points and that of random-interlacement occupation-time field for the avoided points. The results demonstrate universality of the Gaussian Free Field for these extremal problems.


title: The Local Density Approximation in Density Functional Theory

abstract: We give the first mathematically rigorous justification of the Local Density Approximation in Density Functional Theory. We provide a quantitative estimate on the difference between the grand-canonical Levy-Lieb energy of a given density (the lowest possible energy of all quantum states having this density) and the integral over the Uniform Electron Gas energy of this density. The error involves gradient terms and justifies the use of the Local Density Approximation in the situation where the density is very flat on sufficiently large regions in space.

categories: math-ph cond-mat.stat-mech math.MP physics.chem-ph quant-ph


title: SacCalib: Reducing Calibration Distortion for Stationary Eye Trackers Using Saccadic Eye Movements

abstract: Recent methods to automatically calibrate stationary eye trackers were shown to effectively reduce inherent calibration distortion. However, these methods require additional information, such as mouse clicks or on-screen content. We propose the first method that only requires users' eye movements to reduce calibration distortion in the background while users naturally look at an interface. Our method exploits that calibration distortion makes straight saccade trajectories appear curved between the saccadic start and end points. We show that this curving effect is systematic and the result of distorted gaze projection plane. To mitigate calibration distortion, our method undistorts this plane by straightening saccade trajectories using image warping. We show that this approach improves over the common six-point calibration and is promising for reducing distortion. As such, it provides a non-intrusive solution to alleviating accuracy decrease of eye tracker during long-term use.


title: Geometric and numerical methods for a state constrained minimum time control problem of an electric vehicle

abstract: In this article, the minimum time control problem of an electric vehicle is modeled as a Mayer problem in optimal control, with affine dynamics with respect to the control and with state constraints. The candidates as minimizers are selected among a set of extremals, solutions of a Hamiltonian system given by the maximum principle. An analysis, with the techniques of geometric control, is used first to reduce the set of candidates and then to construct the numerical methods. This leads to a numerical investigation based on indirect methods using the HamPath software. Multiple shooting and homotopy techniques are used to build a synthesis with respect to the bounds of the boundary sets.


title: Exploration of Interesting Dense Regions in Spatial Data

abstract: Nowadays, spatial data are ubiquitous in various fields of science, such as transportation and the social Web. A recent research direction in analyzing spatial data is to provide means for "exploratory analysis" of such data where analysts are guided towards interesting options in consecutive analysis iterations. Typically, the guidance component learns analyst's preferences using her explicit feedback, e.g., picking a spatial point or selecting a region of interest. However, it is often the case that analysts forget or don't feel necessary to explicitly express their feedback in what they find interesting. Our approach captures implicit feedback on spatial data. The approach consists of observing mouse moves (as a means of analyst's interaction) and also the explicit analyst's interaction with data points in order to discover interesting spatial regions with dense mouse hovers. In this paper, we define, formalize and explore Interesting Dense Regions (IDRs) which capture preferences of analysts, in order to automatically find interesting spatial highlights. Our approach involves a polygon-based abstraction layer for capturing preferences. Using these IDRs, we highlight points to guide analysts in the analysis process. We discuss the efficiency and effectiveness of our approach through realistic examples and experiments on Airbnb and Yelp datasets.


abstract: The chapter "Solar Active Region Electric Currents Before and During Eruptive Flares" is a discussion on electric currents in the pre-eruption state and in the course of eruptions of solar magnetic structures, using information from solar observations, nonlinear force-free field extrapolations relying on these observations, and three-dimensional magnetohydrodynamic (MHD) models. The discussion addresses the issue of neutralized vs. non-neutralized currents in active regions and concludes that MHD models are able to explain non-neutralized currents in active regions by the existence of strong magnetic shear along the polarity inversion lines, thus confirming previous observations that already contained this result. The models have also captured the essence of the behavior of electric currents in active regions during solar eruptions, predicting current-density increases and decreases inside flare ribbons and in the interior of expanding flux ropes respectively. The observed photospheric current density maps, inferred from vector magnetic field observations, exhibit similar whirling ribbon patterns to the MHD model results, that are interpreted as the signatures of flux ropes and of quasi-separatrix layers (QSLs) between the magnetic systems in active regions. Enhancement of the total current in these QSLs during the eruptions and decreasing current densities at the footpoint of erupting flux ropes, has been confirmed in the observations.


title: Demand Prediction for Electric Vehicle Sharing

abstract: Electric Vehicle (EV) sharing systems have recently experienced unprecedented growth across the globe. Many car sharing service providers as well as automobile manufacturers are entering this competition by expanding both their EV fleets and renting/returning station networks, aiming to seize a share of the market and bring car sharing to the zero emissions level. During their fast expansion, one fundamental determinant for success is the capability of dynamically predicting the demand of stations. In this paper we propose a novel demand prediction approach, which is able to model the dynamics of the system and predict demand accordingly. We use a local temporal encoding process to handle the available historical data at individual stations, and a spatial encoding process to take correlations between stations into account with graph convolutional neural networks. The encoded features are fed to a prediction network, which forecasts both the long-term expected demand of the stations. We evaluate the proposed approach on real-world data collected from a major EV sharing platform. Experimental results demonstrate that our approach significantly outperforms the state of the art.


title: Space-time coupled evolution equations and their stochastic solutions

abstract: We consider a class of space-time coupled evolution equations (CEEs), obtained by a subordination of the heat operator. Our CEEs reformulate and extend known governing equations of non-Markovian processes arising as scaling limits of continuous time random walks, with widespread applications. In particular we allow for initial conditions imposed on the past, general spatial operators on Euclidean domains and a forcing term. We prove existence, uniqueness and stochastic representation for solutions.


title: Affordance Learning for End-to-End Visuomotor Robot Control

abstract: Training end-to-end deep robot policies requires a lot of domain-, task-, and hardware-specific data, which is often costly to provide. In this work, we propose to tackle this issue by employing a deep neural network with a modular architecture, consisting of separate perception, policy, and trajectory parts. Each part of the system is trained fully on synthetic data or in simulation. The data is exchanged between parts of the system as low-dimensional latent representations of affordances and trajectories. The performance is then evaluated in a zero-shot transfer scenario using Franka Panda robot arm. Results demonstrate that a low-dimensional representation of scene affordances extracted from an RGB image is sufficient to successfully train manipulator policies. We also introduce a method for affordance dataset generation, which is easily generalizable to new tasks, objects and environments, and requires no manual pixel labeling.


title: Asymptotically faster algorithm for counting self-avoiding walks and self-avoiding polygons

abstract: We give an algorithm for counting self-avoiding walks or self-avoiding polygons that runs in time $\exp(C\sqrt{n\log n})$ on 2-dimensional lattices and time $\exp(C_dn^{(d-1)/d}\log n)$ on $d$-dimensional lattices for $d>2$.


title: Does Unit-Tested Code Crash? A Case Study of Eclipse

abstract: Context: Software development projects increasingly adopt unit testing as a way to identify and correct program faults early in the construction process. Code that is unit tested should therefore have fewer failures associated with it. Objective: Compare the number of field failures arising in unit tested code against those arising in code that has not been unit tested. Method: We retrieved 2,083,979 crash incident reports associated with the Eclipse integrated development environment project, and processed them to obtain a set of 126,026 unique program failure stack traces associated with a specific popular release. We then run the JaCoCo code test coverage analysis on the same release, obtaining results on the line, instruction, and branch-level coverage of 216,392 methods. We also extracted from the source code the classes that are linked to a corresponding test class so as to limit test code coverage results to 1,267 classes with actual tests. Finally, we correlated unit tests with failures at the level of 9,523 failing tested methods. Results: Unit-tested code does not appear to be associated with fewer failures. Conclusion: Unit testing on its own may not be a sufficient method for preventing program failures.


title: One-Pass Sparsified Gaussian Mixtures

abstract: We present a one-pass sparsified Gaussian mixture model (SGMM). Given $N$ data points in $P$ dimensions, $X$, the model fits $K$ Gaussian distributions to $X$ and (softly) classifies each point to these clusters. After paying an up-front cost of $\mathcal{O}(NP\log P)$ to precondition the data, we subsample $Q$ entries of each data point and discard the full $P$-dimensional data. SGMM operates in $\mathcal{O}(KNQ)$ time per iteration for diagonal or spherical covariances, independent of $P$, while estimating the model parameters in the full $P$-dimensional space, making it one-pass and hence suitable for streaming data. We derive the maximum likelihood estimators for the parameters in the sparsified regime, demonstrate clustering on synthetic and real data, and show that SGMM is faster than GMM while preserving accuracy.


title: Likelihood-free MCMC with Amortized Approximate Ratio Estimators

abstract: Posterior inference with an intractable likelihood is becoming an increasingly common task in scientific domains which rely on sophisticated computer simulations. Typically, these forward models do not admit tractable densities forcing practitioners to make use of approximations. This work introduces a novel approach to address the intractability of the likelihood and the marginal model. We achieve this by learning a flexible amortized estimator which approximates the likelihood-to-evidence ratio. We demonstrate that the learned ratio estimator can be embedded in MCMC samplers to approximate likelihood-ratios between consecutive states in the Markov chain, allowing us to draw samples from the intractable posterior. Techniques are presented to improve the numerical stability and to measure the quality of an approximation. The accuracy of our approach is demonstrated on a variety of benchmarks against well-established techniques. Scientific applications in physics show its applicability.


title: Classifying the unknown: discovering novel gravitational-wave detector glitches using similarity learning

abstract: The observation of gravitational waves from compact binary coalescences by LIGO and Virgo has begun a new era in astronomy. A critical challenge in making detections is determining whether loud transient features in the data are caused by gravitational waves or by instrumental or environmental sources. The citizen-science project \emph{Gravity Spy} has been demonstrated as an efficient infrastructure for classifying known types of noise transients (glitches) through a combination of data analysis performed by both citizen volunteers and machine learning. We present the next iteration of this project, using similarity indices to empower citizen scientists to create large data sets of unknown transients, which can then be used to facilitate supervised machine-learning characterization. This new evolution aims to alleviate a persistent challenge that plagues both citizen-science and instrumental detector work: the ability to build large samples of relatively rare events. Using two families of transient noise that appeared unexpectedly during LIGO's second observing run (O2), we demonstrate the impact that the similarity indices could have had on finding these new glitch types in the Gravity Spy program.


title: Hidden tail chains and recurrence equations for dependence parameters associated with extremes of higher-order Markov chains

abstract: We derive some key extremal features for $k$th order Markov chains that can be used to understand how the process moves between an extreme state and the body of the process. The chains are studied given that there is an exceedance of a threshold, as the threshold tends to the upper endpoint of the distribution. Unlike previous studies with $k>1$, we consider processes where standard limit theory describes each extreme event as a single observation without any information about the transition to and from the body of the distribution. Our work uses different asymptotic theory which results in non-degenerate limit laws for such processes. We study the extremal properties of the initial distribution and the transition probability kernel of the Markov chain under weak assumptions for broad classes of extremal dependence structures that cover both asymptotically dependent and asymptotically independent Markov chains. For chains with $k>1$, the transition of the chain away from the exceedance involves novel functions of the $k$ previous states, in comparison to just the single value, when $k=1$. This leads to an increase in the complexity of determining the form of this class of functions, their properties and the method of their derivation in applications. We find that it is possible to derive an affine normalization, dependent on the threshold excess, such that non-degenerate limiting behaviour of the process is assured for all lags. These normalization functions have an attractive structure that has parallels to the Yule-Walker equations. Furthermore, the limiting process is always linear in the innovations. We illustrate the results with the study of $k$th order stationary Markov chains with exponential margins based on widely studied families of copula dependence structures.


title: Stackelberg Independence

abstract: The standard model of sequential capacity choices is the Stackelberg quantity leadership model with linear demand. I show that under the standard assumptions, leaders' actions are informative about market conditions and independent of leaders' beliefs about the arrivals of followers. However, this Stackelberg independence property relies on all standard assumptions being satisfied. It fails to hold whenever the demand function is non-linear, marginal cost is not constant, goods are differentiated, firms are non-identical, or there are any externalities. I show that small deviations from the linear demand assumption may make the leaders' choices completely uninformative.


title: Stern-Gerlach splitting of low-energy ion beams

abstract: We present a feasibility study with several magnetic field configurations for creating spin-dependent forces that can split a low-energy ion beam by the Stern-Gerlach effect. To the best of our knowledge, coherent spin-splittings of charged particles have yet to be realised. Our proposal is based on ion source parameters taken from a recent experiment that demonstrated single-ion implantation from a high-brightness ion source combined with a radio-frequency Paul trap. The inhomogeneous magnetic fields can be created by permanently magnetised microstructures or from current-carrying wires with sizes in the micron range, such as those recently used in a successful implementation of the Stern-Gerlach effect with neutral atoms. All relevant forces (Lorentz force and image charges) are taken into account, and measurable splittings are found by analytical and numerical calculations.


title: Moser Polynomials and Eulerian Numbers

abstract: Article presents a short investigation into some properties of the Moser polynomials which appear in various problems from algebraic combinatorics. For instance, these polynomials can be used to solve the Generalized Moser's Problem on multiset recovery: Can a collection (multiset) of $n$ numbers can be uniquely restored given the collection of its $s$-sums? We prove some explicit formulas showing relationships between Moser polynomials and such popular algebraic combinatorial sequences as Eulerian and Stirling numbers.


title: What is the nature and origin of the highest-energy particles in the universe?

abstract: This white paper was submitted to the US Astronomy and Astrophysics Decadal Survey (Astro2020) and defines the science questions to be answered in the next decade in the field of Ultra-High Energy Cosmic-Rays. Following a review of the recent experimental and theoretical advances in the field, the paper outlines strategies and requirements desirable for the design of future experiments.


title: Sliced Wasserstein Discrepancy for Unsupervised Domain Adaptation

abstract: In this work, we connect two distinct concepts for unsupervised domain adaptation: feature distribution alignment between domains by utilizing the task-specific decision boundary and the Wasserstein metric. Our proposed sliced Wasserstein discrepancy (SWD) is designed to capture the natural notion of dissimilarity between the outputs of task-specific classifiers. It provides a geometrically meaningful guidance to detect target samples that are far from the support of the source and enables efficient distribution alignment in an end-to-end trainable fashion. In the experiments, we validate the effectiveness and genericness of our method on digit and sign recognition, image classification, semantic segmentation, and object detection.


title: Accretion in Stellar-Mass Black Holes at High X-ray Spectral Resolution

abstract: Accretion disks around stellar-mass black holes offer unique opportunities to study the fundamental physics of standard thin disks, super-Eddington disks, and structure that may be connected to flux variability. These local analogues of active galactic nuclei (AGN) are particularly attractive for their proximity, high flux, and peak emissivity in the X-ray band. X-ray calorimeter spectrometers, with energy resolutions of 2-5 eV, are ideally suited to study accretion in stellar-mass black holes. The results will make strong tests of seminal disk theory that applies in a broad range of circumstances, help to drive new numerical simulations, and will inform our understanding of AGN fueling, evolution, and feedback.


title: Strong pinning in the hole-doped pnictide superconductor La$_{0.34}$Na$_{0.66}$Fe$_2$As$_2$

abstract: We present magnetization studies as a function of time, temperature and magnetic field for $H$ $\parallel$ c-axis, in a hole-doped pnictide superconductor, La$_{0.34}$Na$_{0.66}$Fe$_2$As$_2$, with, $T_c$ $\approx$ 27 K. The obtained vortex phase-diagram shows that the magnetic irreversibility line is very close to the mean-field superconducting transition line, similar to the low $T_c$ superconductors, evidencing a strong pinning behavior. The irreversibility line does not follow a power law behavior with ($T_c$-$T$), however, it is well described using an expression developed in the literature considering the effect of disorder in the system. The critical current density estimated using the Bean's critical-state model is found to be of the order of 10$^5$ A/cm$^2$ below 12 K in the limit of zero magnetic field. A plot of the normalized pinning force density as a function of the reduced magnetic field at different temperatures shows a good scaling and the analysis suggests that the vortex pinning is due to normal point like pinning centers. The temperature dependence of the critical current density suggests that the pinning due to the variation in charge carrier mean free path alone is not sufficient to explain the experimental data. Magnetic relaxation rate as a function of temperature and magnetic field is also studied.


title: Hybrid photonic-crystal fiber

abstract: This article offers an extensive survey of results obtained using hybrid photonic crystal fibers (PCFs) which constitute one of the most active research fields in contemporary fiber optics. The ability to integrate novel and functional materials in solid- and hollow-core PCFs through various post-processing methods has enabled new directions towards understanding fundamental linear and nonlinear phenomena as well as novel application aspects, within the fields of optoelectronics, material and laser science, remote sensing and spectroscopy. Here the recent progress in the field of hybrid PCFs is reviewed from scientific and technological perspectives, focusing on how different fluids, solids and gases can significantly extend the functionality of PCFs. In the first part of this review we discuss the most important efforts by research groups around the globe to develop tunable linear and nonlinear fiber-optic devices using PCFs infiltrated with various liquids, glasses, semiconductors and metals. The second part is concentrated on the most recent and state-of-the-art advances in the field of gas-filled hollow-core PCFs. Extreme ultrafast gas-based nonlinear optics towards light generation in the extreme wavelength regions of vacuum ultraviolet (VUV), pulse propagation and compression dynamics in both atomic and molecular gases, and novel soliton - plasma interactions are reviewed. A discussion of future prospects and directions is also included.


title: Electrical generation and propagation of spin waves in antiferromagnetic thin films

abstract: Electrical generation of THz spin waves is theoretically explored in an antiferromangetic nanostrip via the current-induced spin-orbit torque. The analysis based on micromagnetic simulations clearly illustrates that the Neel-vector oscillations excited at one end of the magnetic strip can propagate in the form of a traveling wave when the nanostrip axis aligns with the magnetic easy-axis. A sizable threshold is observed in the driving current density or the torque to overcome the unfavorable anisotropy as expected. The generated spin waves are found to travel over a long distance while the angle of rotation undergoes continuous decay in the presence of non-zero damping. The oscillation frequency is tunable via the strength of the spin-orbit torque, reaching the THz regime. Other key characteristics of the spin waves such as the phase and the chirality can also be modulated actively. The simulation results further indicate the possibility of wave-like superposition between the excited spin oscillations, illustrating its application as an efficient source of spin-wave signals for information processing.


title: What we can learn from multi-band observations of black hole binaries

abstract: The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to ground-based GW detectors. This white paper discusses some of the synergies that result when the same binaries are observed by instruments in space and on the ground. We consider intermediate-mass black hole binaries (with total mass $M \sim 10^2 -10^4 M_{\odot}$) as well as stellar-mass black hole binaries. We illustrate how combining space-based and ground-based data sets can break degeneracies and thereby improve our understanding of the binary's physical parameters. While early work focused on how space-based observatories can forecast precisely when some mergers will be observed on the ground, the reverse is also important: ground-based detections will allow us to "dig deeper" into archived, space-based data to confidently identify black hole inspirals whose signal-to-noise ratios were originally sub-threshold, increasing the number of binaries observed in both bands by a factor of $\sim 4 - 7$.


title: Orbital stabilization of nonlinear systems via Mexican sombrero energy shaping and pumping-and-damping injection

abstract: In this paper we show that a slight modification to the widely popular interconnection and damping assignment passivity-based control method---originally proposed for stabilization of equilibria of nonlinear systems---allows us to provide a solution to the more challenging orbital stabilization problem. Two different, though related, ways how this procedure can be applied are proposed. First, the assignment of an energy function that has a minimum in a closed curve, i.e., with the shape of a Mexican sombrero. Second, the use of a damping matrix that changes "sign" according to the position of the state trajectory relative to the desired orbit, that is, pumping or dissipating energy. The proposed methodologies are illustrated with the example of the induction motor and prove that it yields the industry standard field oriented control.


title: Constraining the Johannsen deformation parameter $\epsilon_3$ with black hole X-ray data

abstract: We extend our previous work on tests of the Kerr black hole hypothesis using X-ray reflection spectroscopy and we infer observational constraints on the Johannsen deformation parameter $\epsilon_3$. We analyze a NuSTAR observation of the stellar-mass black hole in GS 1354-645, Suzaku data of the supermassive black holes in Ark 564, 1H0419-577, and PKS 0558-504, and some simultaneous XMM-Newton and NuSTAR observations of the supermassive black hole in MCG-6-30-15. All our measurements of $\epsilon_3$ are consistent with the Kerr metric and we find quite strong constraints from Ark 564 and MCG-6-30-15. We discuss the implications of our results and the next steps to improve our tests.


title: 30-micron sources in galaxies with different metallicities

abstract: We present an analysis and comparison of the 30 micron dust features seen in the Spitzer Space Telescope spectra of 207 carbon-rich asymptotic giant branch (AGB) stars, post-AGB objects, and planetary nebulae located in the Milky Way, the Magellanic Clouds (MCs), or the Sagittarius dwarf spheroidal galaxy (Sgr dSph), which are characterised by different average metallicities. We investigated whether the formation of the 30 micron feature carrier may be a function of the metallicity. Through this study we expect to better understand the late stages of stellar evolution of carbon-rich stars in these galaxies. Our analysis uses the `Manchester method' as a basis for estimating the temperature of dust for the carbon-rich AGB stars and the planetary nebulae in our sample. We used a black-body function with a single temperature deduced from the Manchester method or its modification to approximate the continuum under the 30 micron feature. The most important conclusion of our work is the fact that the formation of the 30 micron feature is affected by metallicity. Specifically that, as opposed to more metal-poor samples of AGB stars in the MCs, the feature is seen at lower mass-loss rates, higher temperatures, and has seen to be more prominent in Galactic carbon stars. The averaged feature (profile) in the AGB, post-AGB objects, and PNe seems unaffected by metallicity at least between a fifth and solar metallicity, but in the case of PNe it is shifted to significantly longer wavelengths.


title: An Adaptive Observer Design for Charge-State and Crossover Estimation in Disproportionation Redox Flow Batteries undergoing Self-Discharge

abstract: This article considers a model formulation and an adaptive observer design for the simultaneous estimation of the state of charge and crossover flux in disproportionation redox flow batteries. This novel nonaqueous battery chemistry allows a simple isothermal lumped parameter model to be formulated. The transport of vanadium through the porous separator is a key unknown function of battery variables and it is approximated in the space of continuous functions. The state and parameter observer adaptation laws are derived using Lyapunov analysis applied to the estimation error, the stability and convergence of which are proved. Numerical values of observer gains are calculated by solving a polytopic linear matrix inequality and equality problem via convex optimization. The performance of this design is evaluated on a laboratory flow battery prototype, and it is shown that the crossover flux can be considered a linear function of state of charge for this battery configuration during self-discharge.


title: Dilogarithm identities for solutions to Pell's equation in terms of continued fraction convergents

abstract: In this paper we give describe a new connection between the dilogarithm function and solutions to Pell's equation $x^2-ny^2 = \pm 1$. For each solution $x,y$ to Pell's equation we obtain a dilogarithm identity whose terms are given by the continued fraction expansion of the associated unit $x+y\sqrt{n} \in \Z[\sqrt{n}]$. We further show that Ramanujan's dilogarithm value-identities correspond to an identity for the regular ideal hyperbolic hexagon.


title: Nature of System Calls in CPU-centric Computing Paradigm

abstract: Modern operating systems are typically POSIX-compliant with major system calls specified decades ago. The next generation of non-volatile memory (NVM) technologies raise concerns about the efficiency of the traditional POSIX-based systems. As one step toward building high performance NVM systems, we explore the potential dependencies between system call performance and major hardware components (e.g., CPU, memory, storage) under typical user cases (e.g., software compilation, installation, web browser, office suite) in this paper. We build histograms for the most frequent and time-consuming system calls with the goal to understand the nature of distribution on different platforms. We find that there is a strong dependency between the system call performance and the CPU architecture. On the other hand, the type of persistent storage plays a less important role in affecting the performance.


title: Cluster structure of light nuclei

abstract: We review recent studies of the cluster structure of light nuclei within the framework of the algebraic cluster model (ACM) for nuclei composed of k alpha-particles and within the framework of the cluster shell model (CSM) for nuclei composed of k alpha-particles plus x additional nucleons. The calculations, based on symmetry considerations and thus for the most part given in analytic form, are compared with experiments in light cluster nuclei. The comparison shows evidence for Z_2, D_{3h} and T_d symmetry in the even-even nuclei 8Be (k=2), 12C (k=3) and 16O (k=4), respectively, and for the associated double groups Z'_2 and D'_{3h} in the odd nuclei 9Be, 9B (k=2, x=1) and 13C (k=3, x=1), respectively.


title: Sidechains and interoperability

abstract: There appears to be an insatiable desire for spawning new bespoke blockchains to harness the functionality provided by blockchain technologies, resulting in a constant stream of blockchain start-up companies entering the market with their own unique vision and mission. Some target a particular niche market such as supply chain and financial services, while others strive to differentiate themselves from the increasingly saturated market by offering new functionality. This dynamic and constantly changing blockchain ecosystem makes it very challenging to keep abreast of all the latest breakthroughs and research. It is evident that there is also a growing desire to collaborate with others developing blockchain solutions, which brings new impetus to blockchain interoperability research. We review the strategies that some key players in the blockchain ecosystem have implemented, or are proposing to develop, to satisfy this increasing demand for cross-chain communication and transactions between sidechains. Interoperability presents a complex and challenging stumbling block to the wider uptake of blockchain technology. We find that although there is a plethora of blockchains and interoperability implementations, or proposals, at a higher level of abstraction there is only a handful of approaches. However, the way they are implemented can differ quite substantially. We present a summary of the reviews we conducted in a table for ease of comparing and contrasting.


title: The First Luminous Quasars and Their Host Galaxies

abstract: The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.


title: Nano-scale imaging of the full strain tensor of specific dislocations extracted from a bulk sample

abstract: Lattice defects play a key role in determining the properties of crystalline materials. Probing the 3D lattice strains that govern their interactions remains a challenge. Bragg Coherent Diffraction Imaging (BCDI) allows strain to be measured with nano-scale 3D resolution. However, it is currently limited to materials that form micro-crystals. Here we introduce a new technique that allows the manufacture of BCDI samples from bulk materials. Using tungsten as an example, we show that focussed ion beam (FIB) machining can be used to extract, from macroscopic crystals, micron-sized BCDI samples containing specific pre-selected defects. To interpret the experimental data, we develop a new displacement-gradient-based analysis for multi-reflection BCDI. This allows accurate recovery of the full lattice strain tensor from samples containing multiple dislocations. These new capabilities open the door to BCDI as a microscopy tool for studying complex real-world materials.


title: [Plasma 2020 Decadal] The Material Properties of Weakly Collisional, High-Beta Plasmas

abstract: This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.


title: Redditors in Recovery: Text Mining Reddit to Investigate Transitions into Drug Addiction

abstract: Increasing rates of opioid drug abuse and heightened prevalence of online support communities underscore the necessity of employing data mining techniques to better understand drug addiction using these rapidly developing online resources. In this work, we obtain data from Reddit, an online collection of forums, to gather insight into drug use/misuse using text data from users themselves. Specifically, using user posts, we trained 1) a binary classifier which predicts transitions from casual drug discussion forums to drug recovery forums and 2) a Cox regression model that outputs likelihoods of such transitions. In doing so, we found that utterances of select drugs and certain linguistic features contained in one's posts can help predict these transitions. Using unfiltered drug-related posts, our research delineates drugs that are associated with higher rates of transitions from recreational drug discussion to support/recovery discussion, offers insight into modern drug culture, and provides tools with potential applications in combating the opioid crisis.


title: The Paracrystalline Nature of Lattice Distortion in a High Entropy Alloy

abstract: Severe lattice distortion is suggested for high entropy alloys (HEAs), however, evidence for such effect so far is lacking, and the nature of distortion is yet to be understood. Here, we reveal the distortion in an fcc HEA, Al0.1CrFeCoNi, by direct imaging using electron nanodiffraction. Information about crystal symmetry, lattice strain and atomic distortion are data-mined and mapped from many (~10^4) diffraction patterns. Application to the HEA reveals two embodiments of distortion, nm-sized mosaic blocks of paracrystals and strained nano-clusters. Their interaction gives rise to fractal strain field across nanoscopic to mesoscopic scales. As lattice distortion impedes dislocation motion and contributes to strengthening, results here thus provide critical insights about the complex nature of distortion in a HEA.


title: The Advanced X-ray Imaging Satellite

abstract: Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.


title: Exploring OpenStreetMap Availability for Driving Environment Understanding

abstract: With the great achievement of artificial intelligence, vehicle technologies have advanced significantly from human centric driving towards fully automated driving. An intelligent vehicle should be able to understand the driver's perception of the environment as well as controlling behavior of the vehicle. Since high digital map information has been available to provide rich environmental context about static roads, buildings and traffic infrastructures, it would be worthwhile to explore map data capability for driving task understanding. Alternative to commercial used maps, the OpenStreetMap (OSM) data is a free open dataset, which makes it unique for the exploration research. This study is focused on two tasks that leverage OSM for driving environment understanding. First, driving scenario attributes are retrieved from OSM elements, which are combined with vehicle dynamic signals for the driving event recognition. Utilizing steering angle changes and based on a Bi-directional Recurrent Neural Network (Bi-RNN), a driving sequence is segmented and classified as lane-keeping, lane-change-left, lane-change-right, turn-left, and turn-right events. Second, for autonomous driving perception, OSM data can be used to render virtual street views, represented as prior knowledge to fuse with vision/laser systems for road semantic segmentation. Five different types of road masks are generated from OSM, images, and Lidar points, and fused to characterize the drivable space at the driver's perspective. An alternative data-driven approach is based on a Fully Convolutional Network (FCN), OSM availability for deep learning methods are discussed to reveal potential usage on compensating street view images and automatic road semantic annotation.


title: Real Polynomial Gram Matrices Without Real Spectral Factors

abstract: It is well known that a non-negative definite polynomial matrix (a polynomial Gramian) $G(t)$ can be written as a product of its polynomial spectral factors, $G(t) = X(t)^H X(t)$. In this paper, we give a new algebraic characterization of spectral factors when $G(t)$ is real-valued. The key idea is to construct a representation set that is in bijection with the set of real polynomial Gramians. We use the derived characterization to identify the set of all complex polynomial matrices that generate real-valued Gramians, and we formulate a conjecture that typical rank-deficient real polynomial Gramians have real spectral factors.


title: Quantitative acoustic models for superfluid circuits

abstract: We experimentally realize a highly tunable superfluid oscillator circuit in a quantum gas of ultracold atoms and develop and verify a simple lumped-element description of this circuit. At low oscillator currents, we demonstrate that the circuit is accurately described as a Helmholtz resonator, a fundamental element of acoustic circuits. At larger currents, the breakdown of the Helmholtz regime is heralded by a turbulent shedding of vortices and density waves. Although a simple phase-slip model offers qualitative insights into the circuit's resistive behavior, our results indicate deviations from the phase-slip model. A full understanding of the dissipation in superfluid circuits will thus require the development of empirical models of the turbulent dynamics in this system, as have been developed for classical acoustic systems.


title: Uniting the wave and the particle in quantum mechanics

abstract: We present a unified field theory of wave and particle in quantum mechanics. This emerges from an investigation of three weaknesses in the de Broglie-Bohm (deBB) theory: its reliance on the quantum probability formula to justify the particle law; its insouciance regarding the absence of reciprocal action of the particle on the guiding wave; and its lack of a unified model to represent its inseparable components. These problems are resolved within an analytical framework by requiring that the wave-particle composite exhibits no observable differences with a quantum system. This scheme is implemented by appealing to symmetries (global gauge and spacetime translations) and imposing equality of the corresponding conserved Noether densities (matter, momentum and energy) with their Schrodinger counterparts. In conjunction with the condition of time reversal covariance this implies the deBB law for the particle where the quantum potential mediates the wave-particle interaction (we also show how the time reversal assumption may be replaced by a statistical condition). The method clarifies the nature of the mass of the composite, and its energy-momentum conservation laws. Our principal result is the unification of the Schrodinger equation and the deBB law in an inhomogeneous equation whose solution amalgamates the wavefunction and a singular soliton model of the particle in a unified spacetime field. The wavefunction suffers no reaction from the particle since it is the homogeneous part of the unified field to whose source the particle contributes. The theory applies to many-body systems. We review the objections of de Broglie to the pilot-wave theory and suggest that our field-theoretic model provides a realization of his hitherto unfulfilled double-solution programme. A revised set of postulates for the deBB theory is given in which the unified field is taken as the basic descriptive element.


title: Probing the Jet Turnover Frequency Dependence on Black Hole Mass and Mass Accretion Rate

abstract: We have examined a sample of 13 sub-Eddington supermassive black holes hosted by galaxies spanning a variety of morphological classifications to further understand the empirical fundamental plane of black hole activity. This plane describes black holes from stellar-mass to supermassive and relates the mass of an accreting black hole and its radio and X-ray luminosities. A key factor in studying the fundamental plane is the turnover frequency, the frequency at which the radio continuum emission becomes optically thin. We measured this turnover frequency using new VLA observations combined, when necessary, with archival Chandra observations. Radio observations are in the range of 5--40 GHz across four frequency bands in B-configuration, giving high spatial resolution to focus on the core emission. We use Markov Chain Monte Carlo methods to fit the continuum emission in order to find the turnover frequency. After testing for correlations, the turnover frequency does not display a significant dependence on either mass or mass accretion rate, indicating that more complicated physics than simple scaling and optical depth effects are at play, as has been suggested by recent theoretical work.


title: From conjugacy classes in the Weyl group to semisimple conjugacy classes

abstract: Suppose $G$ is a connected complex semisimple group and $W$ is its Weyl group. The lifting of an element of $W$ to $G$ is semisimple. This induces a well-defined map from the set of elliptic conjugacy classes of $W$ to the set of semisimple conjugacy classes of $G$. In this paper, we give a uniform algorithm to compute this map. We also consider the twisted case.


title: A Hybrid Framework for Action Recognition in Low-Quality Video Sequences

abstract: Vision-based activity recognition is essential for security, monitoring and surveillance applications. Further, real-time analysis having low-quality video and contain less information about surrounding due to poor illumination, and occlusions. Therefore, it needs a more robust and integrated model for low quality and night security operations. In this context, we proposed a hybrid model for illumination invariant human activity recognition based on sub-image histogram equalization enhancement and k-key pose human silhouettes. This feature vector gives good average recognition accuracy on three low exposure video sequences subset of original actions video datasets. Finally, the performance of the proposed approach is tested over three manually downgraded low qualities Weizmann action, KTH, and Ballet Movement dataset. This model outperformed on low exposure videos over existing technique and achieved comparable classification accuracy to similar state-of-the-art methods.


title: A study of Tate homology via the approximation theory with applications to the depth formula

abstract: In this paper we are concerned with absolute, relative and Tate Tor modules. In the first part of the paper we generalize a result of Avramov and Martsinkovsky by using the Auslander-Buchweitz approximation theory, and obtain a new exact sequence connecting absolute Tor modules with relative and Tate Tor modules. In the second part of the paper we consider a depth equality, called the depth formula, which has been initially introduced by Auslander and developed further by Huneke and Wiegand. As an application of our main result, we generalize a result of Yassemi and give a new sufficient condition implying the depth formula to hold for modules of finite Gorenstein and finite injective dimension.


title: MTRNet: A Generic Scene Text Eraser

abstract: Text removal algorithms have been proposed for uni-lingual scripts with regular shapes and layouts. However, to the best of our knowledge, a generic text removal method which is able to remove all or user-specified text regions regardless of font, script, language or shape is not available. Developing such a generic text eraser for real scenes is a challenging task, since it inherits all the challenges of multi-lingual and curved text detection and inpainting. To fill this gap, we propose a mask-based text removal network (MTRNet). MTRNet is a conditional adversarial generative network (cGAN) with an auxiliary mask. The introduced auxiliary mask not only makes the cGAN a generic text eraser, but also enables stable training and early convergence on a challenging large-scale synthetic dataset, initially proposed for text detection in real scenes. What's more, MTRNet achieves state-of-the-art results on several real-world datasets including ICDAR 2013, ICDAR 2017 MLT, and CTW1500, without being explicitly trained on this data, outperforming previous state-of-the-art methods trained directly on these datasets.


title: Restriction estimates to complex hypersurfaces

abstract: The restriction problem is better understood for hypersurfaces and recent progresses have been made by bilinear and multilinear approaches and most recently polynomial partitioning method which is combined with those estimates. However, for surfaces with codimension bigger than 1, bilinear and multilinear generalization of restriction estimates are more involved and effectiveness of these multilinear estimates is not so well understood yet. Regarding the restriction problem for the surfaces with codimensions bigger than 1, the current state of the art is still at the level of $TT^*$ method which is known to be useful for obtaining $L^q$--$L^2$ restriction estimates. In this paper, we consider a special type of codimension 2 surfaces which are given by graphs of complex analytic functions and attempt to make progress beyond the $L^2$ restriction estimates.


title: Light and strange quark masses from $N_f=2+1$ simulations with Wilson fermions

abstract: We present a nearly final analysis of the $u/d$ and $s$ quark masses, extracted using the PCAC quark masses reported in [PRD 95 (2017) 074504]. The data is based on the CLS $N_f = 2 + 1$ simulations with Wilson/Clover quarks and L\"uscher-Weisz gauge action, at four $\beta$ values (i.e. lattice spacings) and a range of quark masses. We use the ALPHA results of [EPJC 78 (2018) 387] for non-perturbative quark mass renormalisation and RG-running from hadronic to electroweak scales in the Schr\"odinger Functional scheme. Quark masses are quoted both in the $\overline{\rm MS}$ scheme and as RGI quantities.


title: Earth wind as a possible source of lunar surface hydration

abstract: Understanding the sources of lunar water is crucial for studying the history of lunar evolution, and also the solar wind interaction with the Moon and other airless bodies. Recent observations revealed lunar hydration is very likely a surficial dynamic process driven by solar wind. Solar wind is shielded over a period of 3-5 days as the Moon passes through the Earth's magnetosphere, during which a significant loss of hydration is expected from previous works.Here we study lunar hydration inside the magnetosphere using orbital spectral data, which unexpectedly found that the polar surficial OH/H2O abundance remains at the same level when in the solar wind and in the magnetosphere. We suggest that particles from the magnetosphere (Earth wind, naturally different from solar wind) contribute to lunar hydration. From lunar orbital plasma observations, we find the existence of optimal energy ranges, other than 1 keV as previously thought, for surface hydration formation. These optimal energy ranges deduced from space observations may provide strong implications for laboratory experiments simulating lunar hydration processes.


title: On Separable Quadratic Lyapunov Functions for Convex Design of Distributed Controllers

abstract: We consider the problem of designing a stabilizing and optimal static controller with a pre-specified sparsity pattern. Since this problem is NP-hard in general, it is necessary to resort to approximation approaches. In this paper, we characterize a class of convex restrictions of this problem that are based on designing a separable quadratic Lyapunov function for the closed-loop system. This approach generalizes previous results based on optimizing over diagonal Lyapunov functions, thus allowing for improved feasibility and performance. Moreover, we suggest a simple procedure to compute favourable structures for the Lyapunov function yielding high-performance distributed controllers. Numerical examples validate our results.


title: Bus Manufacturing Workshop Scheduling Method with Routing Buffer

abstract: Aiming at solving the problem that the moving route is complicated and the scheduling is difficult in the routing buffer of the bus in the manufacturing workshop, a routing buffer mathematical programming model for bus manufacturing workshop is proposed. We design a moving approach for minimizing the total setup cost for moving in routing buffer. The framework and the solution ofthe optimization problem of such a bus manufacturing workshop scheduling with routing buffer arepresented. The evaluation results show that, comparing with the irregularly guided moving method, the proposed method can better guide the bus movement in routing buffer by reducing the total setup time of all buses processed at the next stage, and obtaining a better scheduling optimization solution with minimize maximum total completion time.


title: Zariski $N$-ples for a smooth cubic and its tangent lines

abstract: In this paper, we study the geometry of two-torsion points of elliptic curves in order to distinguish the embedded topology of reducible plane curves consisting of a smooth cubic and its tangent lines. As a result, we obtain a new family of Zariski N-ples consisting of such curves.


title: The fast scalar auxiliary variable approach with unconditional energy stability for nonlocal Cahn-Hilliard equation

abstract: Comparing with the classical local gradient flow and phase field models, the nonlocal models such as nonlocal Cahn-Hilliard equations equipped with nonlocal diffusion operator can describe more practical phenomena for modeling phase transitions. In this paper, we construct an accurate and efficient scalar auxiliary variable approach for the nonlocal Cahn-Hilliard equation with general nonlinear potential. The first contribution is that we have proved the unconditional energy stability for nonlocal Cahn-Hilliard model and its semi-discrete schemes carefully and rigorously. Secondly, what we need to focus on is that the non-locality of the nonlocal diffusion term will lead the stiffness matrix to be almost full matrix which generates huge computational work and memory requirement. For spatial discretizaion by finite difference method, we find that the discretizaition for nonlocal operator will lead to a block-Toeplitz-Toeplitz-block (BTTB) matrix by applying four transformation operators. Based on this special structure, we present a fast procedure to reduce the computational work and memory requirement. Finally, several numerical simulations are demonstrated to verify the accuracy and efficiency of our proposed schemes.


title: Conformal Symplectic and Relativistic Optimization

abstract: Arguably, the two most popular accelerated or momentum-based optimization methods in machine learning are Nesterov's accelerated gradient and Polyaks's heavy ball, both corresponding to different discretizations of a particular second order differential equation with friction. Such connections with continuous-time dynamical systems have been instrumental in demystifying acceleration phenomena in optimization. Here we study structure-preserving discretizations for a certain class of dissipative (conformal) Hamiltonian systems, allowing us to analyze the symplectic structure of both Nesterov and heavy ball, besides providing several new insights into these methods. Moreover, we propose a new algorithm based on a dissipative relativistic system that normalizes the momentum and may result in more stable/faster optimization. Importantly, such a method generalizes both Nesterov and heavy ball, each being recovered as distinct limiting cases, and has potential advantages at no additional cost.


title: Large Scale Learning of Agent Rationality in Two-Player Zero-Sum Games

abstract: With the recent advances in solving large, zero-sum extensive form games, there is a growing interest in the inverse problem of inferring underlying game parameters given only access to agent actions. Although a recent work provides a powerful differentiable end-to-end learning frameworks which embed a game solver within a deep-learning framework, allowing unknown game parameters to be learned via backpropagation, this framework faces significant limitations when applied to boundedly rational human agents and large scale problems, leading to poor practicality. In this paper, we address these limitations and propose a framework that is applicable for more practical settings. First, seeking to learn the rationality of human agents in complex two-player zero-sum games, we draw upon well-known ideas in decision theory to obtain a concise and interpretable agent behavior model, and derive solvers and gradients for end-to-end learning. Second, to scale up to large, real-world scenarios, we propose an efficient first-order primal-dual method which exploits the structure of extensive-form games, yielding significantly faster computation for both game solving and gradient computation. When tested on randomly generated games, we report speedups of orders of magnitude over previous approaches. We also demonstrate the effectiveness of our model on both real-world one-player settings and synthetic data.


title: Blameworthiness in Multi-Agent Settings

abstract: We provide a formal definition of blameworthiness in settings where multiple agents can collaborate to avoid a negative outcome. We first provide a method for ascribing blameworthiness to groups relative to an epistemic state (a distribution over causal models that describe how the outcome might arise). We then show how we can go from an ascription of blameworthiness for groups to an ascription of blameworthiness for individuals using a standard notion from cooperative game theory, the Shapley value. We believe that getting a good notion of blameworthiness in a group setting will be critical for designing autonomous agents that behave in a moral manner.