title: Non-local Corrections to Collisional Transport in Magnetised Plasmas

abstract: In modern inertial fusion experiments there is a complex interplay between non-locality and magnetisation that can greatly influence transport. In this work we use a matrix recursion method to include higher-order corrections beyond the diffusion approximation usually used for magnetised plasmas. Working in the linear regime, we show this can account for arbitrary orders of the distribution function expansion in Knudsen non-locality parameter $k\lambda_{ei}$. Transport coefficients, such as thermal conductivity, deviate significantly from the magnetised diffusive approximation. In particular we show how higher orders of the expansion contribute to transport asynchronously parallel, perpendicular and cross-perpendicular to a uniform magnetic field perpendicular to plasma perturbation.


title: Reducing approximate Longest Common Subsequence to approximate Edit Distance

abstract: Given a pair of strings, the problems of computing their Longest Common Subsequence and Edit Distance have been extensively studied for decades. For exact algorithms, LCS and Edit Distance (with character insertions and deletions) are equivalent; the state of the art running time is (almost) quadratic and this is tight under plausible fine-grained complexity assumptions. But for approximation algorithms the picture is different: there is a long line of works with improved approximation factors for Edit Distance, but for LCS (with binary strings) only a trivial $1/2$-approximation was known. In this work we give a reduction from approximate LCS to approximate Edit Distance, yielding the first efficient $(1/2+\epsilon)$-approximation algorithm for LCS for some constant $\epsilon>0$.


title: What Storage Access Privacy is Achievable with Small Overhead?

abstract: Oblivious RAM (ORAM) and private information retrieval (PIR) are classic cryptographic primitives used to hide the access pattern to data whose storage has been outsourced to an untrusted server. Unfortunately, both primitives require considerable overhead compared to plaintext access. For large-scale storage infrastructure with highly frequent access requests, the degradation in response time and the exorbitant increase in resource costs incurred by either ORAM or PIR prevent their usage. In an ideal scenario, a privacy-preserving storage protocols with small overhead would be implemented for these heavily trafficked storage systems to avoid negatively impacting either performance and/or costs. In this work, we study the problem of the best $\mathit{storage\ access\ privacy}$ that is achievable with only $\mathit{small\ overhead}$ over plaintext access. To answer this question, we consider $\mathit{differential\ privacy\ access}$ which is a generalization of the $\mathit{oblivious\ access}$ security notion that are considered by ORAM and PIR. Quite surprisingly, we present strong evidence that constant overhead storage schemes may only be achieved with privacy budgets of $\epsilon = \Omega(\log n)$. We present asymptotically optimal constructions for differentially private variants of both ORAM and PIR with privacy budgets $\epsilon = \Theta(\log n)$ with only $O(1)$ overhead. In addition, we consider a more complex storage primitive called key-value storage in which data is indexed by keys from a large universe (as opposed to consecutive integers in ORAM and PIR). We present a differentially private key-value storage scheme with $\epsilon = \Theta(\log n)$ and $O(\log\log n)$ overhead. This construction uses a new oblivious, two-choice hashing scheme that may be of independent interest.


title: Energy-Based Continuous Inverse Optimal Control

abstract: The problem of continuous inverse optimal control (over finite time horizon) is to learn the unknown cost function over the sequence of continuous control variables from expert demonstrations. In this article, we study this fundamental problem in the framework of energy-based model, where the observed expert trajectories are assumed to be random samples from a probability density function defined as the exponential of the negative cost function up to a normalizing constant. The parameters of the cost function are learned by maximum likelihood via an "analysis by synthesis" scheme, which iterates (1) synthesis step: sample the synthesized trajectories from the current probability density using the Langevin dynamics via back-propagation through time, and (2) analysis step: update the model parameters based on the statistical difference between the synthesized trajectories and the observed trajectories. Given the fact that an efficient optimization algorithm is usually available for an optimal control problem, we also consider a convenient approximation of the above learning method, where we replace the sampling in the synthesis step by optimization. Moreover, to make the sampling or optimization more efficient, we propose to train the energy-based model simultaneously with a top-down trajectory generator via cooperative learning, where the trajectory generator is used to fast initialize the synthesis step of the energy-based model. We demonstrate the proposed methods on autonomous driving tasks, and show that they can learn suitable cost functions for optimal control.


title: Robust Two-Step phase estimation using the Simplified Lissajous Ellipse Fitting method with Gabor Filter Banks preprocessing

abstract: We present the Simplified Lissajous Ellipse Fitting (SLEF) method for the calculation of the random phase step and the phase distribution from two phase-shifted interferograms. We consider interferograms with spatial and temporal dependency of background intensities, amplitude modulations and noise. Given these problems, the use of the Gabor Filters Bank (GFB) allows us to filter--out the noise, normalize the amplitude and eliminate the background. The normalized patterns permit to implement the SLEF algorithm, which is based on reducing the number of estimated coefficients of the ellipse equation, from five terms to only two. Our method consists of three stages. First, we preprocess the interferograms with GFB methodology in order to normalize the fringe patterns. Second, we calculate the phase step by using the proposed SLEF technique and third, we estimate the phase distribution using a two--steps formula. For the calculation of the phase step, we present two alternatives: the use of the Least Squares (LS) method to approximate the values of the coefficients and, in order to improve the LS estimation, a robust estimation based on the Leclerc's potential. The SLEF method's performance is evaluated through synthetic and experimental data to demonstrate its feasibility.


title: Enhancement of spin transparency by interfacial alloying

abstract: We report that atomic-layer alloying (intermixing) at a Pt/Co interface can increase, by approximately 30%, rather than degrade the interfacial spin transparency, and thereby strengthen the efficiency of the dampinglike spin-orbit torque arising from the spin Hall effect in the Pt. At the same time, this interfacial alloying substantially reduces fieldlike spin-orbit torque. Insertion of an ultrathin magnetic alloy layer at heavy-metal/ferromagnet interfaces represents an effective approach for improving interfacial spin transparency that may enhance not only spin-orbit torques but also the spin detection efficiency in inverse spin Hall experiments.


title: R-Storm: Resource-Aware Scheduling in Storm

abstract: The era of big data has led to the emergence of new systems for real-time distributed stream processing, e.g., Apache Storm is one of the most popular stream processing systems in industry today. However, Storm, like many other stream processing systems lacks an intelligent scheduling mechanism. The default round-robin scheduling currently deployed in Storm disregards resource demands and availability, and can therefore be inefficient at times. We present R-Storm (Resource-Aware Storm), a system that implements resource-aware scheduling within Storm. R-Storm is designed to increase overall throughput by maximizing resource utilization while minimizing network latency. When scheduling tasks, R-Storm can satisfy both soft and hard resource constraints as well as minimizing network distance between components that communicate with each other. We evaluate R-Storm on set of micro-benchmark Storm applications as well as Storm applications used in production at Yahoo! Inc. From our experimental results we conclude that R-Storm achieves 30-47% higher throughput and 69-350% better CPU utilization than default Storm for the micro-benchmarks. For the Yahoo! Storm applications, R-Storm outperforms default Storm by around 50% based on overall throughput. We also demonstrate that R-Storm performs much better when scheduling multiple Storm applications than default Storm.


title: Instance Segmentation based Semantic Matting for Compositing Applications

abstract: Image compositing is a key step in film making and image editing that aims to segment a foreground object and combine it with a new background. Automatic image compositing can be done easily in a studio using chroma-keying when the background is pure blue or green. However, image compositing in natural scenes with complex backgrounds remains a tedious task, requiring experienced artists to hand-segment. In order to achieve automatic compositing in natural scenes, we propose a fully automated method that integrates instance segmentation and image matting processes to generate high-quality semantic mattes that can be used for image editing task. Our approach can be seen both as a refinement of existing instance segmentation algorithms and as a fully automated semantic image matting method. It extends automatic image compositing techniques such as chroma-keying to scenes with complex natural backgrounds without the need for any kind of user interaction. The output of our approach can be considered as both refined instance segmentations and alpha mattes with semantic meanings. We provide experimental results which show improved performance results as compared to existing approaches.


title: Using Multiwavelength Variability to Explore the Connection between X-ray Emission, the Far-Ultraviolet H2 Bump, and Accretion in T Tauri Stars

abstract: The high-energy radiation fields of T Tauri stars (TTS) should affect the surrounding circumstellar disk, having implications for disk transport and heating. Yet, observational evidence of the effect of high-energy fields on disks is scarce. Here we investigate the connection between X-ray emission and the innermost gas disk by leveraging the variability of TTS. We obtained multiple epochs of coordinated data (taken either simultaneously or within a few hours) of accreting TTS with the Hubble Space Telescope, the Neil Gehrels Swift Observatory, and the Chandra X-ray Observatory. We measured the far-ultraviolet (FUV) H2 bump feature at 1600 A, which traces gas <1 AU from the star; the near-ultraviolet (NUV) emission, from which we extract the accretion luminosity; and also the X-ray luminosity. We do not find a correlation between the FUV H2 bump and X-ray luminosity. Therefore, an observable tracer of the effect of X-ray ionization in the innermost disk remains elusive. We report a correlation between the FUV H2 bump and accretion luminosity, linking this feature to the disk surface density. We also see a correlation between the X-ray luminosity and the accretion column density, implying that flaring activity may influence accretion. These results stress the importance of coordinated multiwavelength work to understand TTS.


title: Constant factor approximations to edit distance on far input pairs in nearly linear time

abstract: For any $T \geq 1$, there are constants $R=R(T) \geq 1$ and $\zeta=\zeta(T)>0$ and a randomized algorithm that takes as input an integer $n$ and two strings $x,y$ of length at most $n$, and runs in time $O(n^{1+\frac{1}{T}})$ and outputs an upper bound $U$ on the edit distance $ED(x,y)$ that with high probability, satisfies $U \leq R(ED(x,y)+n^{1-\zeta})$. In particular, on any input with $ED(x,y) \geq n^{1-\zeta}$ the algorithm outputs a constant factor approximation with high probability. A similar result has been proven independently by Brakensiek and Rubinstein (2019).


title: Least Squares Auto-Tuning

abstract: Least squares is by far the simplest and most commonly applied computational method in many fields. In almost all applications, the least squares objective is rarely the true objective. We account for this discrepancy by parametrizing the least squares problem and automatically adjusting these parameters using an optimization algorithm. We apply our method, which we call least squares auto-tuning, to data fitting.


title: Less is More: Real-time Failure Localization in Power Systems

abstract: Cascading failures in power systems exhibit non-local propagation patterns which make the analysis and mitigation of failures difficult. In this work, we propose a distributed control framework inspired by the recently proposed concepts of unified controller and network tree-partition that offers strong guarantees in both the mitigation and localization of cascading failures in power systems. In this framework, the transmission network is partitioned into several control areas which are connected in a tree structure, and the unified controller is adopted by generators or controllable loads for fast timescale disturbance response. After an initial failure, the proposed strategy always prevents successive failures from happening, and regulates the system to the desired steady state where the impact of initial failures are localized as much as possible. For extreme failures that cannot be localized, the proposed framework has a configurable design, that progressively involves and coordinates more control areas for failure mitigation and, as a last resort, imposes minimal load shedding. We compare the proposed control framework with Automatic Generation Control (AGC) on the IEEE 118-bus test system. Simulation results show that our novel framework greatly improves the system robustness in terms of the N-1 security standard, and localizes the impact of initial failures in majority of the load profiles that are examined. Moreover, the proposed framework incurs significantly less load loss, if any, compared to AGC, in all of our case studies.


title: The Spinful Large Charge Sector of Non-Relativistic CFTs: From Phonons to Vortex Crystals

abstract: We study operators in Schr\"odinger invariant field theories (non-relativistic conformal field theories or NRCFTs) with large charge (particle number) and spin. Via the state-operator correspondence for NRCFTs, such operators correspond to states of a superfluid in a harmonic trap with phonons or vortices. Using the effective field theory of the Goldstone mode, we compute the dimensions of operators to leading order in the angular momentum $L$ and charge $Q$. We find a diverse set of scaling behaviors for NRCFTs in both $d=2$ and $d=3$.


title: Optomechanically Induced Birefringence and Faraday Effect

abstract: We demonstrate an optomechanical platform where optical mode conversion mediated by mechanical motion enables arbitrary tailoring of polarization states of propagating light fields. Optomechanical interactions are realized in a Fabry-Per\'ot resonator, which naturally supports two polarization-degenerate states while an optical control field induces rotational symmetry breaking. Applying such principles, the entire Poincar\'e sphere is spanned by just optical control of the driving field, realizing reciprocal and non-reciprocal optomechanically-induced birefringence for linearly polarized and circularly polarized control driving. A straightforward extension of this setup also enables all-optical tunable isolation and circulation. Our findings open new avenues to exploit optomechanics for arbitrary manipulation of light polarization.


title: Vehicle Routing for the Last-Mile Logistics Problem

abstract: Energy consumption is the major contributor associated with large and growing transportation cost in logistics. Optimal vehicle routing approaches can provide solutions to reduce their operating costs and address implications on energy. This paper outlines a solution to the single-depot capacitated vehicle routing problem with the objective of minimizing daily operation cost with a homogeneous fleet of delivery vehicles. The problem is solved using Simulated Annealing, to provide optimal routes for the vehicles traveling between the depot and destinations. Simulation results demonstrate that the proposed approach is effective to recommend an optimal route and reduce operation cost. Supplementary information and video of our proposed approach can be found at: https://sites.google.com/view/ud-ids-lab/last-mile


title: Reconstruction of tunnel exit time and exit momentum in strong field ionization, based on phase space methods

abstract: We analyze tunnel ionization of a single atom based on the Wigner function over the classical phase space which inspires improved classical electron trajectories: these start with exit momenta based on the quantum momentum function and correspond very well to the subsequent quantum evolution. We derive an approximate analytic formula to reconstruct the tunnel exit time and exit momentum from electron momentum data that can be measured e.g. with a usual time-of-flight electron detector.


title: Exact sequences on Powell-Sabin splits

abstract: We construct smooth finite elements spaces on Powell-Sabin triangulations that form an exact sequence. The first space of the sequence coincides with the classical $C^1$ Powell-Sabin space, while the others form stable and divergence-free yielding pairs for the Stokes problem. We develop degrees of freedom for these spaces that induce projections that commute with the differential operators.


title: Improved visual function in a case of ultra-low vision following ischemic encephalopathy following transcranial electrical stimulation; A case study

abstract: Cortical visual impairment is amongst the key pathological causes of pediatric visual abnormalities predominantly resulting from hypoxic-ischemic brain injury. Such an injury results in profound visual impairments which severely impairs the patient's quality of life. Given the nature of the pathology, treatments are mostly limited to rehabilitation strategies such as transcranial electrical stimulation and visual rehabilitation therapy. Here, we discussed an 11 year-old girl case with cortical visual impairment who underwent concurrent visual rehabilitation therapy and transcranial electrical stimulation resulting in her improved visual function. This novel and noninvasive therapeutic intervention has shown potential for application in neuro-visual rehabilitation therapy (nVRT).


title: Incidence estimates for well spaced tubes

abstract: We prove analogues of the Szemer\'edi-Trotter theorem and other incidence theorems using $\delta$-tubes in place of straight lines, assuming that the $\delta$-tubes are well-spaced in a strong sense.


title: Exponential cosmological solutions with two factor spaces in EGB model with $\Lambda = 0$ revisited

abstract: We study exact cosmological solutions in $D$-dimensional Einstein-Gauss-Bonnet model (with zero cosmological term) governed by two non-zero constants: $\alpha_1$ and $\alpha_2$. We deal with exponential dependence (in time) of two scale factors governed by Hubble-like parameters $H > 0$ and $h$, which correspond to factor spaces of dimensions $m > 2$ and $l > 2$, respectively, and $D = 1 + m + l$. We put $h \neq H$ and $m H + l h \neq 0$. We show that for $\alpha = \alpha_2/\alpha_1 > 0$ there are two (real) solutions with two sets of Hubble-like parameters: $(H_1, h_1)$ and $(H_2, h_2)$, which obey: $ h_1/ H_1 < - m/l < h_2/ H_2 < 0$, while for $\alpha < 0$ the (real) solutions are absent. We prove that the cosmological solution corresponding to $(H_2, h_2)$ is stable in a class of cosmological solutions with diagonal metrics, while the solution corresponding to $(H_1, h_1)$ is unstable. We present several examples of analytical solutions, e.g. stable ones with small enough variation of the effective gravitational constant $G$, for $(m, l) = (9, l > 2), (12, 11), (11,16), (15, 6)$.


title: Crystalline Nodal Topological Superconductivity and Bogolyubov Fermi Surfaces in Monolayer NbSe$_2$

abstract: We present a microscopic calculation of the phase diagram of the Ising superconductor NbSe$_{2}$ in presence of both in-plane magnetic field and Rashba spin-orbit coupling (SOC). Repulsive interactions lead to two distinct instabilities, in singlet- and triplet- interaction channels. While we recover the previously predicted nodal topological superconducting state in the absence of Rashba SOC at large magnetic field with six pairs of nodes along \(\Gamma\)-\(M\) lines, a finite Rashba SOC breaks the symmetry that protects these nodes and therefore generally lifts them, resulting in a topologically trivial phase. There is an exception when the field is applied along one of the three $\Gamma$-$K$ lines, however. In that case, a single mirror symmetry remains that can protect two pairs of nodes out of the original six, resulting in a \emph{crystalline} topological superconducting phase. Depending on the Cooper pairs' center-of-mass momentum, this superconducting state displays either Bogolyubov Fermi surfaces or point nodes. Moreover, a chiral topological superconducting phase with Chern number of 6 is realized in the regime of large Rashba SOC and dominant triplet interactions, spontaneously breaking time-reversal symmetry.


title: Signatures for quark matter from multi-messenger observations

abstract: We review the prospects for detecting quark matter in neutron star cores. We survey the proposed signatures and emphasize the importance of data from neutron star mergers, which provide access to dynamical properties that operate on short timescales that are not probed by other neutron star observables.


title: Theory of Cryptocurrency Interest Rates

abstract: A term structure model in which the short rate is zero is developed as a candidate for a theory of cryptocurrency interest rates. The price processes of crypto discount bonds are worked out, along with expressions for the instantaneous forward rates and the prices of interest-rate derivatives. The model admits functional degrees of freedom that can be calibrated to the initial yield curve and other market data. Our analysis suggests that strict local martingales can be used for modelling the pricing kernels associated with virtual currencies based on distributed ledger technologies.


title: Dimension transcendence and anomalous charge transport in magnets with moving multiple-$Q$ spin textures

abstract: Multiple-$Q$ spin textures, such as magnetic bubble and skyrmion lattices, can be driven into motion by external stimuli. The motion of spin textures affects the electronic states. Here we show that to describe correctly the electronic dynamics, the momentum space needs to be transcended to higher dimensions by including the ancillary dimensions associated with phason modes of the translational motion of the spin textures. The electronic states have non-trivial topology characterized by the first and second Chern numbers in the high dimensional hybrid momentum space. This gives rise to an anomalous electric charge transport due to the motion of spin textures. By deforming the spin textures, a nonlinear response associated with the second Chern number can be induced. The charge transport is derived from the semi-classical equation of motion of electrons that depends on the Berry curvature in the hybrid momentum space. Our results suggest that the motion of multiple-$Q$ spin textures has significant effects on the electronic dynamics and provides a new platform to explore high dimensional topological physics.


title: Efficient Distributed Workload (Re-)Embedding

abstract: Modern networked systems are increasingly reconfigurable, enabling demand-aware infrastructures whose resources can be adjusted according to the workload they currently serve. Such dynamic adjustments can be exploited to improve network utilization and hence performance, by moving frequently interacting communication partners closer, e.g., collocating them in the same server or datacenter. However, dynamically changing the embedding of workloads is algorithmically challenging: communication patterns are often not known ahead of time, but must be learned. During the learning process, overheads related to unnecessary moves (i.e., re-embeddings) should be minimized. This paper studies a fundamental model which captures the tradeoff between the benefits and costs of dynamically collocating communication partners on $\ell$ servers, in an online manner. Our main contribution is a distributed online algorithm which is asymptotically almost optimal, i.e., almost matches the lower bound (also derived in this paper) on the competitive ratio of any (distributed or centralized) online algorithm. As an application, we show that our algorithm can be used to solve a distributed union find problem in which the sets are stored across multiple servers.


title: Learning to Generate Synthetic Data via Compositing

abstract: We present a task-aware approach to synthetic data generation. Our framework employs a trainable synthesizer network that is optimized to produce meaningful training samples by assessing the strengths and weaknesses of a `target' network. The synthesizer and target networks are trained in an adversarial manner wherein each network is updated with a goal to outdo the other. Additionally, we ensure the synthesizer generates realistic data by pairing it with a discriminator trained on real-world images. Further, to make the target classifier invariant to blending artefacts, we introduce these artefacts to background regions of the training images so the target does not over-fit to them. We demonstrate the efficacy of our approach by applying it to different target networks including a classification network on AffNIST, and two object detection networks (SSD, Faster-RCNN) on different datasets. On the AffNIST benchmark, our approach is able to surpass the baseline results with just half the training examples. On the VOC person detection benchmark, we show improvements of up to 2.7% as a result of our data augmentation. Similarly on the GMU detection benchmark, we report a performance boost of 3.5% in mAP over the baseline method, outperforming the previous state of the art approaches by up to 7.5% on specific categories.


title: Smart Home Survey on Security and Privacy

abstract: Smart homes are a special use-case of the Internet-of-Things (IoT) paradigm. Security and privacy are two prime concern in smart home networks. A threat-prone smart home can reveal lifestyle and behavior of the occupants, which may be a significant concern. This article shows security requirements and threats to a smart home and focuses on a privacy-preserving security model. We classify smart home services based on the spatial and temporal properties of the underlying device-to-device and owner-to-cloud interaction. We present ways to adapt existing security solutions such as distance-bounding protocols, ISO-KE, SIGMA, TLS, Schnorr, Okamoto Identification Scheme (IS), Pedersen commitment scheme for achieving security and privacy in a cloud-assisted home area network.


title: On Entanglement Entropy of Maxwell fields in 3+1 dimensions

abstract: We consider entanglement entropy between two halves of space separated by a plane, in the theory of free photon in 3+1 dimensions. We show how to separate local gauge invariant quantities that belong to the two spatial regions. We calculate the entanglement entropy by integrating over the degrees of freedom in one half space using an approximation that assumes slow variation of the magnetic fields in longitudinal direction. We find that the entropy is proportional to the transverse area as expected. Interestingly the entanglement properties of the 2D transverse and longitudinal modes of magnetic field are quite different. While the transverse fields are entangled mostly in the neighborhood of the separation surface as expected, the longitudinal fields are entangled through an infrared mode which extends to large distances from the entanglement surface. This long range entanglement arises due to necessity to solve the no-monopole constraint condition for magnetic field.


title: Predicting Progression of Age-related Macular Degeneration from Fundus Images using Deep Learning

abstract: Background: Patients with neovascular age-related macular degeneration (AMD) can avoid vision loss via certain therapy. However, methods to predict the progression to neovascular age-related macular degeneration (nvAMD) are lacking. Purpose: To develop and validate a deep learning (DL) algorithm to predict 1-year progression of eyes with no, early, or intermediate AMD to nvAMD, using color fundus photographs (CFP). Design: Development and validation of a DL algorithm. Methods: We trained a DL algorithm to predict 1-year progression to nvAMD, and used 10-fold cross-validation to evaluate this approach on two groups of eyes in the Age-Related Eye Disease Study (AREDS): none/early/intermediate AMD, and intermediate AMD (iAMD) only. We compared the DL algorithm to the manually graded 4-category and 9-step scales in the AREDS dataset. Main outcome measures: Performance of the DL algorithm was evaluated using the sensitivity at 80% specificity for progression to nvAMD. Results: The DL algorithm's sensitivity for predicting progression to nvAMD from none/early/iAMD (78+/-6%) was higher than manual grades from the 9-step scale (67+/-8%) or the 4-category scale (48+/-3%). For predicting progression specifically from iAMD, the DL algorithm's sensitivity (57+/-6%) was also higher compared to the 9-step grades (36+/-8%) and the 4-category grades (20+/-0%). Conclusions: Our DL algorithm performed better in predicting progression to nvAMD than manual grades. Future investigations are required to test the application of this DL algorithm in a real-world clinical setting.


title: On a global supersonic-sonic patch characterized by 2-D steady full Euler equations

abstract: Supersonic-sonic patches are ubiquitous in regions of transonic flows and they boil down to a family of degenerate hyperbolic problems in regions surrounded by a streamline, a characteristic curve and a possible sonic curve. This paper establishes the global existence of solutions in a whole supersonic-sonic patch characterized by the two-dimensional full system of steady Euler equations and studies solution behaviors near sonic curves, depending on the proper choice of boundary data extracted from the airfoil problem and related contexts. New characteristic decompositions are developed for the full system and a delicate local partial hodograph transformation is introduced for the solution estimates. It is shown that the solution is uniformly $C^{1,\frac{1}{6}}$ continuous up to the sonic curve and the sonic curve is also $C^{1,\frac{1}{6}}$ continuous.


title: On the properties of a newborn magnetar powering the X-ray transient CDF-S XT2

abstract: Very recently \citet{XueYQ2019} reported an important detection of the X-ray transient, CDF-S XT2, whose light curve is analogous to X-ray plateau features of gamma-ray burst afterglows. They suggested that this transient is powered by a remnant stable magnetar from a binary neutron star merger since several pieces of evidence (host galaxy, location, and event rate) all point toward such an assumption. In this paper, we revisit this scenario and confirm that this X-ray emission can be well explained by the internal gradual magnetic dissipation process in an ultra-relativistic wind of the newborn magnetar. We show that both the light curve and spectral evolution of CDF-S XT2 can be well fitted by such a model. Furthermore, we can probe some key properties of the central magnetar, such as its initial spin period, surface magnetic field strength and wind saturation Lorentz factor.


title: Stochastic Comparative Statics in Markov Decision Processes

abstract: In multi-period stochastic optimization problems, the future optimal decision is a random variable whose distribution depends on the parameters of the optimization problem. We analyze how the expected value of this random variable changes as a function of the dynamic optimization parameters in the context of Markov decision processes. We call this analysis \emph{stochastic comparative statics}. We derive both \emph{comparative statics} results and \emph{stochastic comparative statics} results showing how the current and future optimal decisions change in response to changes in the single-period payoff function, the discount factor, the initial state of the system, and the transition probability function. We apply our results to various models from the economics and operations research literature, including investment theory, dynamic pricing models, controlled random walks, and comparisons of stationary distributions.


title: PSR J1926-0652: A Pulsar with Interesting Emission Properties Discovered at FAST

abstract: We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.


title: Parallels Between Phase Transitions and Circuit Complexity?

abstract: In many natural average-case problems, there are or there are believed to be critical values in the parameter space where the structure of the space of solutions changes in a fundamental way. These phase transitions are often believed to coincide with drastic changes in the computational complexity of the associated problem. In this work, we study the circuit complexity of inference in the broadcast tree model, which has important applications in phylogenetic reconstruction and close connections to community detection. We establish a number of qualitative connections between phase transitions and circuit complexity in this model. Specifically, we show that there is a $\mathbf{TC}^0$ circuit that competes with the Bayes optimal predictor in some range of parameters above the Kesten-Stigum bound. We also show that there is a $16$ label broadcast tree model beneath the Kesten-Stigum bound in which it is possible to accurately guess the label of the root, but beating random guessing is $\mathbf{NC}^1$-hard on average. The key to locating phase transitions is often to study some intrinsic notions of complexity associated with belief propagation \--- e.g. where do linear statistics fail, or when is the posterior sensitive to noise? Ours is the first work to study the complexity of belief propagation in a way that is grounded in circuit complexity.


title: On Code Design for Wireless Channels with Additive Radar Interference

abstract: This paper considers the problem of code design for a channel where communications and radar systems coexist, modeled as having both Additive White Gaussian Noise (AWGN) and Additive Radar Interference (ARI). The issue of how to adapt or re-design convolutional codes (decoded by the Viterbi algorithm) and LDPC codes (decoded by the sum-product algorithm and optimized by using the EXIT chart method) to effectively handle the overall non-Gaussian ARI noise is investigated. A decoding metric is derived from the non-Gaussian ARI channel transition probability as a function of the Signal-to-Noise Ratio (SNR) and Interference-to-Noise Ratio (INR). Two design methodologies are benchmarked against a baseline "unaltered legacy system", where a code designed for AWGN-only noise, but used on the non-Gaussian ARI channel, is decoded by using the AWGN-only metric (i.e., as if INR is zero). The methodologies are: M1) codes designed for AWGN-only noise, but decoded with the new metric that accounts for both SNR and INR; and M2) codes optimized for the overall non-Gaussian ARI channel. Both methodologies give better average Bit Error Rate (BER) in the high INR regime compared to the baseline. In the low INR regime, both methodologies perform as the baseline since in this case the radar interference is weak. Interestingly, the performance improvement of M2 over M1 is minimal. In practice, this implies that specifications in terms of channel error correcting codes for commercially available wireless systems need not be changed, and that it suffices to use an appropriate INR-based decoding metric in order to effectively cope with the ARI.


title: O-minimal de Rham cohomology

abstract: O-minimal geometry generalizes both semialgebraic and subanalytic geometries, and has been very successful in solving special cases of some problems in arithmetic geometry, such as Andr\'e-Oort conjecture. Among the many tools developed in an o-minimal setting are cohomology theories for abstract-definable continuous manifolds such as singular cohomology, sheaf cohomology and \v Cech cohomology, which have been used for instance to prove Pillay's conjecture concerning definably compact groups. In the present paper we elaborate an o-minimal de Rham cohomology theory for abstract-definable $\mathcal{C}^\infty$ manifolds in an o-minimal expansion of the real field which admits smooth cell decomposition and defines the exponential function. We can specify the o-minimal cohomology groups and attain some properties such as the existence of Mayer-Vietoris sequence and the invariance under abstract-definable $\mathcal{C}^\infty$ diffeomorphisms. However, in order to obtain the invariance of our o-minimal cohomology under abstract-definable homotopy we must, working in a tame context that defines sufficiently many primitives, assume the validity of a statement related to Br\"ocker's question.


title: Spin-exchange-induced exotic superfluids in a Bose-Fermi spinor mixture

abstract: We consider a mixture of spin-1/2 bosons and fermions, where only the bosons are subjected to the spin-orbit coupling induced by Raman beams. The fermions, although not directly coupled to the Raman lasers, acquire an effective spin-orbit coupling through the spin-exchange interaction between the two species. Our calculation shows that this is a promising way of obtaining spin-orbit coupled Fermi gas without Raman-induced heating, where the long-sought topological Fermi superfluids and topological bands can be realized. Conversely, we find that the presence of fermions not only provides a new way to create the supersolid stripe phase of the bosons, but more strikingly it can also greatly increase the spatial period of the bosonic density stripes, and hence makes this phase directly observable in the experiment. This system provides a new and practical platform to explore the physics of spin-orbit coupling, which possesses a dynamic nature through the interaction between the two species.


title: Centralized Coded Caching with User Cooperation

abstract: In this paper, we consider the coded-caching broadcast network with user cooperation, where a server connects with multiple users and the users can cooperate with each other through a cooperation network. We propose a centralized coded caching scheme based on a new deterministic placement strategy and a parallel delivery strategy. It is shown that the new scheme optimally allocate the communication loads on the server and users, obtaining cooperation gain and parallel gain that greatly reduces the transmission delay. Furthermore, we show that the number of users who parallelly send information should decrease when the users' caching size increases. In other words, letting more users parallelly send information could be harmful. Finally, we derive a constant multiplicative gap between the lower bound and upper bound on the transmission delay, which proves that our scheme is order optimal.


title: Deep Neural Network Ensembles

abstract: Current deep neural networks suffer from two problems; first, they are hard to interpret, and second, they suffer from overfitting. There have been many attempts to define interpretability in neural networks, but they typically lack causality or generality. A myriad of regularization techniques have been developed to prevent overfitting, and this has driven deep learning to become the hot topic it is today; however, while most regularization techniques are justified empirically and even intuitively, there is not much underlying theory. This paper argues that to extract the features used in neural networks to make decisions, it's important to look at the paths between clusters existing in the hidden spaces of neural networks. These features are of particular interest because they reflect the true decision making process of the neural network. This analysis is then furthered to present an ensemble algorithm for arbitrary neural networks which has guarantees for test accuracy. Finally, a discussion detailing the aforementioned guarantees is introduced and the implications to neural networks, including an intuitive explanation for all current regularization methods, are presented. The ensemble algorithm has generated state-of-the-art results for Wide-ResNets on CIFAR-10 (top 5 for all models) and has improved test accuracy for all models it has been applied to.


title: Tracking excited states in wave function optimization using density matrices and variational principles

abstract: We present a method for finding individual excited states' energy stationary points in complete active space self-consistent field theory that is compatible with standard optimization methods and highly effective at overcoming difficulties due to root flipping and near-degeneracies. Inspired by both the maximum overlap method and recent progress in excited state variational principles, our approach combines these ideas in order to track individual excited states throughout the orbital optimization process. In a series of tests involving root flipping, near-degeneracies, charge transfers, and double excitations, we show that this approach is more effective for state-specific optimization than either the naive selection of roots based on energy ordering or a more direct generalization of the maximum overlap method. Furthermore, we provide evidence that this state-specific approach improves the performance of complete active space perturbation theory. With a simple implementation, a low cost, and compatibility with large active space methods, the approach is designed to be useful in a wide range of excited state investigations.


title: A centrally concentrated sub-solar mass starless core in the Taurus L1495 filamentary complex

abstract: The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka Morita array) and the IRAM 30 m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H$_2$ volume density of $\sim$10$^6$ cm$^{-3}$, MC5-N, connected to a narrow (width $\sim$0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim$0.2-0.4 $M_{\odot}$, which is an order of magnitude smaller than typical low-mass prestellar cores. Taking into account a typical core to star formation efficiency for prestellar cores ($\sim$20%-40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N$_2$H$^+$ and N$_2$D$^+$ observations using the Nobeyama 45 m telescope have confirmed the high-deuterium fractionation ($\sim$30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.


title: Absence of energy currents in an equilibrium state and chiral anomalies

abstract: A long time ago F. Bloch showed that in a system of interacting non-relativistic particles the net particle-number current must vanish in any equilibrium state. Bloch's argument does not generalize easily to the energy current. We devise an alternative argument which proves the vanishing of the net energy currents in equilibrium states of lattice systems as well as systems of non-relativistic particles with finite-range potential interactions. We discuss some applications of these results. In particular, we show that neither a 1d lattice system nor a 1d system of non-relativistic particles with finite-range potential interactions can flow to a Conformal Field Theory with unequal left-moving and right-moving central charges.


title: New Identities for Padovan Numbers

abstract: In \cite{Choi-Jo}, the $am+b$ ($0\leq b<a$) subscripted Tribonacci numbers are studied. This work is devoted to study a new generalization of Fibonacci numbers called Padovan numbers. In particular, the $am+b$ subscripted Padovan numbers will be expressed by three $a$ step apart Padovan numbers for any $0\leq b<a$, where $a\in \mathbb{Z}$.


title: Non-locally Encoder-Decoder Convolutional Network for Whole Brain QSM Inversion

abstract: Quantitative Susceptibility Mapping (QSM) reconstruction is a challenging inverse problem driven by ill conditioning of its field-to -susceptibility transformation. State-of-art QSM reconstruction methods either suffer from image artifacts or long computation times, which limits QSM clinical translation efforts. To overcome these limitations, a non-locally encoder-decoder gated convolutional neural network is trained to infer whole brain susceptibility map, using the local field and brain mask as the inputs. The performance of the proposed method is evaluated relative to synthetic data, a publicly available challenge dataset, and clinical datasets. The proposed approach can outperform existing methods on quantitative metrics and visual assessment of image sharpness and streaking artifacts. The estimated susceptibility maps can preserve conspicuity of fine features and suppress streaking artifacts. The demonstrated methods have potential value in advancing QSM clinical research and aiding in the translation of QSM to clinical operations.


title: Classical simulation of boson sampling with sparse output

abstract: Boson sampling can simulate physical problems for which classical simulations are inefficient. However, not all problems simulated by boson sampling are classically intractable. We consider a situation in which it is known that the outcome from boson sampling is sparse. It can be determined from a few marginal distributions which are classically calculable. Still, recovering of the joint distribution can be of high complexity. We show classically efficient methods of the recovery assuming high sparsity of the joint distribution. Various extensions are discussed including a version involving quantum annealing.


title: Optimal Nonergodic Sublinear Convergence Rate of Proximal Point Algorithm for Maximal Monotone Inclusion Problems

abstract: We establish the optimal nonergodic sublinear convergence rate of the proximal point algorithm for maximal monotone inclusion problems. First, the optimal bound is formulated by the performance estimation framework, resulting in an infinite dimensional nonconvex optimization problem, which is then equivalently reformulated as a finite dimensional semidefinite programming (SDP) problem. By constructing a feasible solution to the dual SDP, we obtain an upper bound on the optimal nonergodic sublinear rate. Finally, an example in two dimensional space is constructed to provide a lower bound on the optimal nonergodic sublinear rate. Since the lower bound provided by the example matches exactly the upper bound obtained by the dual SDP, we have thus established the worst case nonergodic sublinear convergence rate which is optimal in terms of both the order as well as the constants involved. Our result sharpens the understanding of the fundamental proximal point algorithm.


title: A new small satellite sunspot triggering recurrent standard- and blowout-coronal jets

abstract: In this paper,we report a detailed analysis of recurrent jets originated from a location with emerging, canceling and converging negative magnetic field at the east edge of NOAA active region AR11166 from 2011 March 09 to 10. The event presented several interesting features. First, a satellite sunspot appeared and collided with a pre-existing opposite polarity magnetic field and caused a recurrent solar jet event. Second, the evolution of the jets showed blowout-like nature and standard characteristics. Third, the satellite sunspot exhibited a motion toward southeast of AR11166 and merged with the emerging flux near the opposite polarity sunspot penumbra, which afterward, due to flux convergence and cancellation episodes, caused recurrent jets. Fourth, three of the blowout jets associated with coronal mass ejections (CMEs), were observed from field of view of the Solar Terrestrial Relations Observatory. Fifth, almost all the blowout jet eruptions were accompanied with flares or with more intense brightening in the jet base region, while almost standard jets did not manifest such obvious feature during eruptions. The most important, the blowout jets were inclined to faster and larger scale than the standard jets. The standard jets instead were inclined to relative longer-lasting. The obvious shearing and twisting motions of the magnetic field may be interpreted as due to the shearing and twisting motions for a blowout jet eruption. From the statistical results, about 30% blowout jets directly developed into CMEs. It suggests that the blowout jets and CMEs should have a tight relationship.


title: Generation of vortices and stabilization of vortex lattices in holographic superfluids

abstract: Within the simplest holographic superfluid model and without any ingredients put by hand, it is shown that vortices can be generated when the angular velocity of rotating superfluids exceeds certain critical values, which can be precisely determined by linear perturbation analyses (quasi-normal modes of the bulk AdS black brane). These vortices appear at the edge of the superfluid system first, and then automatically move into the bulk of the system, where they are eventually stabilized into certain vortex lattices. For the case of 18 vortices generated, we find (at least) five different patterns of the final lattices formed due to different initial perturbations, which can be compared to the known result for such lattices in weakly coupled Bose-Einstein condensates from free energy analyses.


title: Synthesizing Database Programs for Schema Refactoring

abstract: Many programs that interact with a database need to undergo schema refactoring several times during their life cycle. Since this process typically requires making significant changes to the program's implementation, schema refactoring is often non-trivial and error-prone. Motivated by this problem, we propose a new technique for automatically synthesizing a new version of a database program given its original version and the source and target schemas. Our method does not require manual user guidance and ensures that the synthesized program is equivalent to the original one. Furthermore, our method is quite efficient and can synthesize new versions of database programs (containing up to 263 functions) that are extracted from real-world web applications with an average synthesis time of 69.4 seconds.


title: On the 2-Adic Complexity of the Ding-Helleseth-Martinsen Binary Sequences

abstract: We determine the 2-adic complexity of the Ding-Helleseth-Martinsen (DHM) binary sequences by using cyclotomic numbers of order four, "Gauss periods" and "quadratic Gauss sum" on finite field $\mathbb{F}_q$ and valued in $\mathbb{Z}_{2^N-1}$ where $q \equiv 5\pmod 8$ is a prime number and $N=2q$ is the period of the DHM sequences.


title: Uniquely-Wilf classes

abstract: Two permutations in a class are Wilf-equivalent if, for every size, $n$, the number of permutations in the class of size $n$ containing each of them is the same. Those infinite classes that have only one equivalence class in each size for this relation are characterised provided either that they avoid at least one permutation of size 3, or at least three permutations of size 4.


title: Inverse Problems of Determining Sources of the Fractional Partial Differential Equations

abstract: In this chapter, we mainly review theoretical results on inverse source problems for diffusion equations with the Caputo time-fractional derivatives of order $\alpha\in(0,1)$. Our survey covers the following types of inverse problems: 1. determination of time-dependent functions in interior source terms 2. determination of space-dependent functions in interior source terms 3. determination of time-dependent functions appearing in boundary conditions


title: Inverse Problems of Determining Parameters of the Fractional Partial Differential Equations

abstract: When considering fractional diffusion equation as model equation in analyzing anomalous diffusion processes, some important parameters in the model related to orders of the fractional derivatives, are often unknown and difficult to be directly measured, which requires one to discuss inverse problems of identifying these physical quantities from some indirectly observed information of solutions. Inverse problems in determining these unknown parameters of the model are not only theoretically interesting, but also necessary for finding solutions to initial-boundary value problems and studying properties of solutions. This chapter surveys works on such inverse problems for fractional diffusion equations.


title: Alcock-Paczynski Test with the Evolution of Redshift-Space Galaxy Clustering Anisotropy

abstract: We develop an improved Alcock-Paczynski (AP) test method that uses the redshift-space two-point correlation function (2pCF) of galaxies. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF, which should not change apart from the expected small non-linear evolution. An incorrect choice of cosmology used to convert redshift to comoving distance will manifest itself as redshift-dependent 2pCF. Our method decomposes the redshift difference of the two-dimensional correlation function into the Legendre polynomials whose amplitudes are modeled by radial fitting functions. Our likelihood analysis with this 2-D fitting scheme tightens the constraints on $\Omega_m$ and ${w}$ by $\sim 40\%$ compared to the method of Li et al. (2016, 2017, 2018) that uses one dimensional angular dependence only. We also find that the correction for the non-linear evolution in the 2pCF has a non-negligible cosmology dependence, which has been neglected in previous similar studies by Li et al.. With an accurate accounting for the non-linear systematics and use of full two-dimensional shape information of the 2pCF down to scales as small as $5~h^{-1}{\rm Mpc}$ it is expected that the AP test with redshift-space galaxy clustering anisotropy can be a powerful method to constrain the expansion history of the universe.


title: Intermediate Phase in Interacting Dirac Fermions with Staggered Potential

abstract: By performing exact quantum Monte Carlo simulations of a model of interacting Dirac Fermions with staggered potential, we reveal a novel intermediate phase where the electronic correlations drive a band insulator metallic, and at a larger interaction, drive the metal to Mott insulator. We also show that the Mott insulating phase is antiferromagnetic. A complete phase diagram is achieved by studying the phase transitions at large staggered potential and interaction strengths, which shows that the intermediate state is robust and occupies a large part of the phase diagram and that it should be more feasible to be detected experimentally.


title: Inverse Problems of Determining Coefficients of the Fractional Partial Differential Equations

abstract: When considering fractional diffusion equation as model equation in analyzing anomalous diffusion processes, some important parameters in the model, for example, the orders of the fractional derivative or the source term, are often unknown, which requires one to discuss inverse problems to identify these physical quantities from some additional information that can be observed or measured practically. This chapter investigates several kinds of inverse coefficient problems for the fractional diffusion equation.


title: Membership Inference Attacks on Sequence-to-Sequence Models: Is My Data In Your Machine Translation System?

abstract: Data privacy is an important issue for "machine learning as a service" providers. We focus on the problem of membership inference attacks: given a data sample and black-box access to a model's API, determine whether the sample existed in the model's training data. Our contribution is an investigation of this problem in the context of sequence-to-sequence models, which are important in applications such as machine translation and video captioning. We define the membership inference problem for sequence generation, provide an open dataset based on state-of-the-art machine translation models, and report initial results on whether these models leak private information against several kinds of membership inference attacks.


title: Topological defects and geometric memory across the nematic-smectic A liquid crystal phase transition

abstract: We study transformations of self-organized defect arrays at the nematic-smectic A liquid crystal phase transition, and show that these defect configurations are correlated, or "remembered", across the phase transition. A thin film of thermotropic liquid crystal is subjected to hybrid anchoring by an air interface and a water substrate, and viewed under polarized optical microscopy. Upon heating from smectic-A to nematic, a packing of focal conic domains melts into a dense array of boojums---nematic surface defects---which then coarsens by pair-annihilation. With the aid of Landau-de Gennes numerical modeling, we elucidate the topological and geometrical rules underlying this transformation. In the transition from nematic to smectic-A, we show that focal conic domain packings are organized over large scales in patterns that retain a geometric memory of the nematic boojum configuration, which can be recovered with remarkable fidelity.


title: Capacity and Energy-Efficiency of Delayed Access Scheme for Small Cell Networks

abstract: Data applications may typically tolerate a moderate delay before packet transmission between user equipment (UE) and cell begins. This delay can be taken advantage to reduce the communication distance, improve coverage probability, and increase overall energy-efficiency of the small cell network. To demonstrate such merits, we suggest a simple access scheme and analyze the distribution of coverage probability and throughput as a function of delay and transmit distance. Sufficient number of small base stations (SBSs) handle the peak traffic load. To improve energy-efficiency of the network, a number of SBSs are switched off at low traffic periods. Energy-efficiency can be further improved by turning all of the SBSs on and off, rather than selecting a subset and leaving them off. By doing so, coverage probability and bit-rate can be improved by delaying their transmissions and waiting for a closer SBS to become available. Results show that by turning SBSs on and off continuously and taking advantage of initial delay to connect a SBS yield an order of magnitude improvement in energy-efficiency, improves the coverage probability significantly at low signal to interference and noise (SINR) regime.


title: CNN-Based Deep Architecture for Reinforced Concrete Delamination Segmentation Through Thermography

abstract: Delamination assessment of the bridge deck plays a vital role for bridge health monitoring. Thermography as one of the nondestructive technologies for delamination detection has the advantage of efficient data acquisition. But there are challenges on the interpretation of data for accurate delamination shape profiling. Due to the environmental variation and the irregular presence of delamination size and depth, conventional processing methods based on temperature contrast fall short in accurate segmentation of delamination. Inspired by the recent development of deep learning architecture for image segmentation, the Convolutional Neural Network (CNN) based framework was investigated for the applicability of delamination segmentation under variations in temperature contrast and shape diffusion. The models were developed based on Dense Convolutional Network (DenseNet) and trained on thermal images collected for mimicked delamination in concrete slabs with different depths under experimental setup. The results suggested satisfactory performance of accurate profiling the delamination shapes.


title: Restricted Isometry Property under High Correlations

abstract: Matrices satisfying the Restricted Isometry Property (RIP) play an important role in the areas of compressed sensing and statistical learning. RIP matrices with optimal parameters are mainly obtained via probabilistic arguments, as explicit constructions seem hard. It is therefore interesting to ask whether a fixed matrix can be incorporated into a construction of restricted isometries. In this paper, we construct a new broad ensemble of random matrices with dependent entries that satisfy the restricted isometry property. Our construction starts with a fixed (deterministic) matrix $X$ satisfying some simple stable rank condition, and we show that the matrix $XR$, where $R$ is a random matrix drawn from various popular probabilistic models (including, subgaussian, sparse, low-randomness, satisfying convex concentration property), satisfies the RIP with high probability. These theorems have various applications in signal recovery, random matrix theory, dimensionality reduction, etc. Additionally, motivated by an application for understanding the effectiveness of word vector embeddings popular in natural language processing and machine learning applications, we investigate the RIP of the matrix $XR^{(l)}$ where $R^{(l)}$ is formed by taking all possible (disregarding order) $l$-way entrywise products of the columns of a random matrix $R$.


title: Designing a broadband and linear polarization metasurface carpet cloak in the visible

abstract: In the past few years, carpet cloaking attracted interests because of its feasibility at optical frequencies and potential in stealth technologies. Metasurfaces have been proposed as a method to engineer ultra-thin carpet cloaking surfaces due to their abilities to manipulate wavefronts, polarization, and phase at subwavelength scale. However, achieving broadband carpet cloaking with a significant bandwidth is one of the key remaining challenges for metasurface designs. To date, broadband carpet cloaking based on metasurfaces has not been achieved and cloaking is limited to discrete wavelengths. Here, we propose and numerically demonstrate a novel metasurface design for broadband carpet cloaking with linear polarization at visible wavelengths from 650 nm to 800 nm. Our proposed method is a promising approach for broadband structured interfaces.


title: Generalizing Monocular 3D Human Pose Estimation in the Wild

abstract: The availability of the large-scale labeled 3D poses in the Human3.6M dataset plays an important role in advancing the algorithms for 3D human pose estimation from a still image. We observe that recent innovation in this area mainly focuses on new techniques that explicitly address the generalization issue when using this dataset, because this database is constructed in a highly controlled environment with limited human subjects and background variations. Despite such efforts, we can show that the results of the current methods are still error-prone especially when tested against the images taken in-the-wild. In this paper, we aim to tackle this problem from a different perspective. We propose a principled approach to generate high quality 3D pose ground truth given any in-the-wild image with a person inside. We achieve this by first devising a novel stereo inspired neural network to directly map any 2D pose to high quality 3D counterpart. We then perform a carefully designed geometric searching scheme to further refine the joints. Based on this scheme, we build a large-scale dataset with 400,000 in-the-wild images and their corresponding 3D pose ground truth. This enables the training of a high quality neural network model, without specialized training scheme and auxiliary loss function, which performs favorably against the state-of-the-art 3D pose estimation methods. We also evaluate the generalization ability of our model both quantitatively and qualitatively. Results show that our approach convincingly outperforms the previous methods. We make our dataset and code publicly available.


title: Complete polarization control for a nanofiber waveguide using directional coupling

abstract: Optical nanofiber waveguides are widely used for near-field delivery and measurement of light. Despite their versatility and efficiency, nanofibers have a critical drawback - their inability to maintain light's polarization state on propagation. Here, we design a directional coupler consisting of two crossed nanofibers to probe the polarization state at the waist region. Directionality of coupling occurs due to asymmetric dipolar emission or spin-locking when the evanescent field pattern breaks the mirror symmetry of the crossed-nanofiber system. We demonstrate that, by monitoring the outputs from the directional coupler, two non-orthogonal polarization states can be prepared at the nanofiber waist with a fidelity higher than 99%. Based on these states, we devise a simple and reliable method for complete control of the polarization along a nanofiber waveguide.


title: Mitigating Information Leakage in Image Representations: A Maximum Entropy Approach

abstract: Image recognition systems have demonstrated tremendous progress over the past few decades thanks, in part, to our ability of learning compact and robust representations of images. As we witness the wide spread adoption of these systems, it is imperative to consider the problem of unintended leakage of information from an image representation, which might compromise the privacy of the data owner. This paper investigates the problem of learning an image representation that minimizes such leakage of user information. We formulate the problem as an adversarial non-zero sum game of finding a good embedding function with two competing goals: to retain as much task dependent discriminative image information as possible, while simultaneously minimizing the amount of information, as measured by entropy, about other sensitive attributes of the user. We analyze the stability and convergence dynamics of the proposed formulation using tools from non-linear systems theory and compare to that of the corresponding adversarial zero-sum game formulation that optimizes likelihood as a measure of information content. Numerical experiments on UCI, Extended Yale B, CIFAR-10 and CIFAR-100 datasets indicate that our proposed approach is able to learn image representations that exhibit high task performance while mitigating leakage of predefined sensitive information.


title: Triangular array of iron-oxide nanoparticles: A simulation study of intra- and inter-particle magnetism

abstract: A study of spherical maghemite nanoparticles on a two dimensional triangular array was carried out using a stochastic Landau-Lifshitz-Gilbert (sLLG) approach. The simulation method was first validated with a triangular array of simple dipoles, where results show the expected phase transition to a ferromagnetic state at a finite temperature. The ground state exhibited a continuous degeneracy that was lifted by an order-from-disorder mechanism at infinitesimal temperatures with the appearance of a six-fold planar anisotropy. The nanoparticle array consisted of 7.5 nm diameter maghemite spheres with bulk-like superexchange interactions between Fe-ions in the core, and weaker exchange between surface Fe-ions and a radial anisotropy. The triangular nanoparticle array ordered at the same reduced temperature as the simple dipole array, but exhibited different behaviour at low temperatures due to the surface anisotropy. We find that the vacancies on the octahedral sites in the nanoparticles combine with the surface anisotropy to produce an effective random temperature-dependent anisotropy for each particle. This leads to a reduction in the net magnetization of the nanoparticle array at zero temperature compared to the simple dipole array.


title: Adaptive Virtual Waveform Design for Millimeter-Wave Joint Communication-Radar

abstract: Joint communication and radar (JCR) waveforms with fully digital baseband generation and processing can now be realized at the millimeter-wave (mmWave) band. Prior work has proposed a mmWave wireless local area network (WLAN)-based JCR that exploits the WLAN preamble for radars. The performance of target velocity estimation, however, was limited. In this paper, we propose a virtual waveform design for an adaptive mmWave JCR. The proposed system transmits a few non-uniformly placed preambles to construct several receive virtual preambles for enhancing velocity estimation accuracy, at the cost of only a small reduction in the communication data rate. We evaluate JCR performance trade-offs using the Cramer-Rao Bound (CRB) metric for radar estimation and a novel distortion minimum mean square error (MMSE) metric for data communication. Additionally, we develop three different MMSE-based optimization problems for the adaptive JCR waveform design. Simulations show that an optimal virtual (non-uniform) waveform achieves a significant performance improvement as compared to a uniform waveform. For a radar CRB constrained optimization, the optimal radar range of operation and the optimal communication distortion MMSE (DMMSE) are improved. For a communication DMMSE constrained optimization with a high DMMSE constraint, the optimal radar CRB is enhanced. For a weighted MMSE average optimization, the advantage of the virtual waveform over the uniform waveform is increased with decreased communication weighting. Comparison of MMSE-based optimization with traditional virtual preamble count-based optimization indicated that the conventional solution converges to the MMSE-based one only for a small number of targets and a high signal-to-noise ratio.


title: Optical Kerr nonlinearity of disordered all-dielectric resonant high index metasurfaces with negative refraction

abstract: The optical Kerr effect of material with negative refraction is estimated for the first time. This is done via three-dimensional finite-difference time-domain (FDTD) simulations of disordered bidisperse metasurfaces consisting of high index (GaP) spheres at the wavelength of 532~nm. The metasurfaces comprise spherical particles randomly arranged on plane having two sizes close to the magnetic and electric dipole Mie resonances. The real part of the effective nonlinear refractive index of the metasurfaces is computed in the vicinity of the Mie resonances where the metasurface possesses the negative index of refraction. The optical Kerr nonlinearity has a peak under the condition for the negative refraction. Intensity-dependent refractive index of the bidisperse metasurfaces is studied through concentration transition to the negative refraction state. It is shown that the nonlinear Kerr coefficient of the monolayer metasurface has maximum when the effective linear refractive index is close to zero.


title: Classifying the formation processes of S0 galaxies using Convolutional Neural Networks

abstract: Numerous studies have demonstrated the ability of Convolutional Neural Networks (CNNs) to classify large numbers of galaxies in a manner which mimics the expertise of astronomers. Such classifications are not always physically motivated, however, such as categorising galaxies by their morphological types. In this work, we consider the use of CNNs to classify simulated S0 galaxies based on fundamental physical properties. In particular, we undertake two investigations: (1) the classification of simulated S0 galaxies into three distinct evolutionary paths (isolated, tidal interaction in a group halo, and Spiral-Spiral merger), and (2) the prediction of the mass ratio for the S0s formed via mergers. To train the CNNs, we first run several hundred N-body simulations to model the formation of S0s under idealised conditions; and then we build our training datasets by creating images of stellar density and two dimensional kinematic maps for each simulated S0. Our trained networks have remarkable accuracies exceeding 99% when classifying the S0 formation pathway. For the case of predicting merger mass ratios, the mean predictions are consistent with the true values to within roughly one standard deviation across the full range of our data. Our work demonstrates the potential of CNNs to classify galaxies by the fundamental physical properties which drive their evolution.


title: Efficient and Robust Registration on the 3D Special Euclidean Group

abstract: We present an accurate, robust and fast method for registration of 3D scans. Our motion estimation optimizes a robust cost function on the intrinsic representation of rigid motions, i.e., the Special Euclidean group $\mathbb{SE}(3)$. We exploit the geometric properties of Lie groups as well as the robustness afforded by an iteratively reweighted least squares optimization. We also generalize our approach to a joint multiview method that simultaneously solves for the registration of a set of scans. We demonstrate the efficacy of our approach by thorough experimental validation. Our approach significantly outperforms the state-of-the-art robust 3D registration method based on a line process in terms of both speed and accuracy. We also show that this line process method is a special case of our principled geometric solution. Finally, we also present scenarios where global registration based on feature correspondences fails but multiview ICP based on our robust motion estimation is successful.


title: Non-static Analysis of the Anomalous Chiral Conductivities

abstract: Given the intrinsic nonequilibrium nature of high-energy collisions the investigation of the dynamical properties of transport phenomena is important. The study of the real-time behavior of various conductivities and susceptibilities help refine the simulation tools we use to compare the theories and the experimental findings. In this contribution we take steps to give the chiral magnetic conductivity in case of magnetic field and chiral imbalance are both space-time dependent. Using linear response approximation we present the general 1-loop resummed expression for the electric current. We also suggest simple limiting cases in hope for possible implementation into a hydrodynamical framework.


title: UniVSE: Robust Visual Semantic Embeddings via Structured Semantic Representations

abstract: We propose Unified Visual-Semantic Embeddings (UniVSE) for learning a joint space of visual and textual concepts. The space unifies the concepts at different levels, including objects, attributes, relations, and full scenes. A contrastive learning approach is proposed for the fine-grained alignment from only image-caption pairs. Moreover, we present an effective approach for enforcing the coverage of semantic components that appear in the sentence. We demonstrate the robustness of Unified VSE in defending text-domain adversarial attacks on cross-modal retrieval tasks. Such robustness also empowers the use of visual cues to resolve word dependencies in novel sentences.


title: Timely-Throughput Optimal Coded Computing over Cloud Networks

abstract: In modern distributed computing systems, unpredictable and unreliable infrastructures result in high variability of computing resources. Meanwhile, there is significantly increasing demand for timely and event-driven services with deadline constraints. Motivated by measurements over Amazon EC2 clusters, we consider a two-state Markov model for variability of computing speed in cloud networks. In this model, each worker can be either in a good state or a bad state in terms of the computation speed, and the transition between these states is modeled as a Markov chain which is unknown to the scheduler. We then consider a Coded Computing framework, in which the data is possibly encoded and stored at the worker nodes in order to provide robustness against nodes that may be in a bad state. With timely computation requests submitted to the system with computation deadlines, our goal is to design the optimal computation-load allocation scheme and the optimal data encoding scheme that maximize the timely computation throughput (i.e, the average number of computation tasks that are accomplished before their deadline). Our main result is the development of a dynamic computation strategy called Lagrange Estimate-and Allocate (LEA) strategy, which achieves the optimal timely computation throughput. It is shown that compared to the static allocation strategy, LEA increases the timely computation throughput by 1.4X - 17.5X in various scenarios via simulations and by 1.27X - 6.5X in experiments over Amazon EC2 clusters


title: Relational flexibility of network elements based on inconsistent community detection

abstract: Community identification of network components enables us to understand the mesoscale clustering structure of networks. A number of algorithms have been developed to determine the most likely community structures in networks. Such a probabilistic or stochastic nature of this problem can naturally involve the ambiguity in resultant community structures. More specifically, stochastic algorithms can result in different community structures for each realization in principle. In this study, instead of trying to "solve" this community degeneracy problem, we turn the tables by taking the degeneracy as a chance to quantify how strong companionship each node has with other nodes. For that purpose, we define the concept of companionship inconsistency that indicates how inconsistently a node is identified as a member of a community regarding the other nodes. Analyzing model and real networks, we show that companionship inconsistency discloses unique characteristics of nodes, thus we suggest it as a new type of node centrality. In social networks, for example, companionship inconsistency can classify outsider nodes without firm community membership and promiscuous nodes with multiple connections to several communities. In infrastructure networks such as power grids, it can diagnose how the connection structure is evenly balanced in terms of power transmission. Companionship inconsistency, therefore, abstracts individual nodes' intrinsic property on its relationship to a higher-order organization of the network.


title: Stochastic vertex corrections: linear scaling methods for accurate quasiparticle energies

abstract: New stochastic approaches for the computation of electronic excitations are developed within the many-body perturbation theory. Three approximations to the electronic self-energy are considered: $G_0W_0$, $G_0W_0^tc$, and $G_0W_0^{tc}\Gamma_x$. All three methods are formulated in the time domain and the latter two incorporate non-local vertex corrections. In case of $G_0W_0^{tc}\Gamma_x$, the vertex corrections are included both in the screened Coulomb interaction and in the expression for the self-energy. The implementation of the three approximations is verified by comparison to deterministic results for a set of small molecules. The performance fully stochastic implementation is tested on acene molecules, C$_{60}$ and PC$_{60}$BM. The vertex correction appears crucial for the description of unoccupied states. Unlike conventional (deterministic) approaches, all three stochastic methods scale linearly with the number of electrons.


title: Subtrees of a Given size in Schroeder Trees

abstract: In this paper we find the generating function for the number of vertices which have $k$ elements in their subtree and use this generating function to calculate the probability that a vertex has a size $k$ subtree. We also show how this same technique can be applied to calculate the probabilities for other trees and specifically apply it to Motzkin trees.


title: A Selective Overview of Deep Learning

abstract: Deep learning has arguably achieved tremendous success in recent years. In simple words, deep learning uses the composition of many nonlinear functions to model the complex dependency between input features and labels. While neural networks have a long history, recent advances have greatly improved their performance in computer vision, natural language processing, etc. From the statistical and scientific perspective, it is natural to ask: What is deep learning? What are the new characteristics of deep learning, compared with classical methods? What are the theoretical foundations of deep learning? To answer these questions, we introduce common neural network models (e.g., convolutional neural nets, recurrent neural nets, generative adversarial nets) and training techniques (e.g., stochastic gradient descent, dropout, batch normalization) from a statistical point of view. Along the way, we highlight new characteristics of deep learning (including depth and over-parametrization) and explain their practical and theoretical benefits. We also sample recent results on theories of deep learning, many of which are only suggestive. While a complete understanding of deep learning remains elusive, we hope that our perspectives and discussions serve as a stimulus for new statistical research.


title: Modeling Global Syntactic Variation in English Using Dialect Classification

abstract: This paper evaluates global-scale dialect identification for 14 national varieties of English as a means for studying syntactic variation. The paper makes three main contributions: (i) introducing data-driven language mapping as a method for selecting the inventory of national varieties to include in the task; (ii) producing a large and dynamic set of syntactic features using grammar induction rather than focusing on a few hand-selected features such as function words; and (iii) comparing models across both web corpora and social media corpora in order to measure the robustness of syntactic variation across registers.


title: Review on DNA Cryptography

abstract: Cryptography is the science that secures data and communication over the network by applying mathematics and logic to design strong encryption methods. In the modern era of e-business and e-commerce the protection of confidentiality, integrity and availability (CIA triad) of stored information as well as of transmitted data is very crucial. Deoxyribonucleic acid (DNA) is a genetic molecule consisting of two linked strands that wind around each other to form a double helical structure. The backbone of each strand is made of alternating deoxyribose sugar and phosphate groups. To each sugar one of four bases are attached i.e., adenine (A), cytosine (C), guanine (G) and thymine (T). DNA molecules, having the capacity to store, process and transmit information, inspires the idea of DNA cryptography. It is the rapid emerging unconventional techniques which combines the chemical characteristics of biological DNA sequences with classical cryptography to ensure non-vulnerable transmission of data. This innovative method is based on the notion of DNA computing. The methodologies of DNA cryptography are not coded mathematically; thus, it could be too secure to be cracked easily.


title: Frequency vs. Association for Constraint Selection in Usage-Based Construction Grammar

abstract: A usage-based Construction Grammar (CxG) posits that slot-constraints generalize from common exemplar constructions. But what is the best model of constraint generalization? This paper evaluates competing frequency-based and association-based models across eight languages using a metric derived from the Minimum Description Length paradigm. The experiments show that association-based models produce better generalizations across all languages by a significant margin.


title: Recurrent Event Network: Autoregressive Structure Inference over Temporal Knowledge Graphs

abstract: Knowledge graph reasoning is a critical task in natural language processing. The task becomes more challenging on temporal knowledge graphs, where each fact is associated with a timestamp. Most existing methods focus on reasoning at past timestamps and they are not able to predict facts happening in the future. This paper proposes Recurrent Event Network (RE-NET), a novel autoregressive architecture for predicting future interactions. The occurrence of a fact (event) is modeled as a probability distribution conditioned on temporal sequences of past knowledge graphs. Specifically, our RE-NET employs a recurrent event encoder to encode past facts and uses a neighborhood aggregator to model the connection of facts at the same timestamp. Future facts can then be inferred in a sequential manner based on the two modules. We evaluate our proposed method via link prediction at future times on five public datasets. Through extensive experiments, we demonstrate the strength of RENET, especially on multi-step inference over future timestamps, and achieve state-of-the-art performance on all five datasets. Code and data can be found at https://github.com/INK-USC/RE-Net.


title: Spectroscopy of broad absorption line quasars at $3\lesssim z \lesssim 5$ -- I: evidence for quasar winds shaping broad/narrow emission line regions

abstract: We present an observational study of 22 broad absorption line quasars (BAL QSOs) at $3\lesssim z \lesssim5$ based on optical/near-IR spectroscopy, aiming to investigate quasar winds and their effects. The near-IR spectroscopy covers the \hb\ and/or \mgii\ broad emission lines (BELs) for these quasars, allowing us to estimate their central black hole (BH) masses in a robust way. We found that our BAL QSOs on average do not have a higher Eddington ratio than that from non-BAL QSOs matched in redshift and/or luminosity. In a subset consisting of seven strong BAL QSOs possessing sub-relativistic BAL outflows, we see the prevalence of large \civ-BEL blueshift ($\sim$3100 km s$^{-1}$) and weak \oiii\ emission (particularly the narrow \oiii$\lambda$5007 component), indicative of nuclear outflows affecting the narrow emission-line (NEL) regions. In another subset consisting of thirteen BAL QSOs having simultaneous observations of \mgii\ and \hb, we found a strong correlation between 3000~\AA\ and 5000~\AA\ monochromatic luminosity, consistent with that from non-BAL QSOs matched in redshift and luminosity; however, there is no correlation between \mgii\ and \hb\ in FWHM, likely due to nuclear outflows influencing the BEL regions. Our spectroscopic investigations offer strong evidence that the presence of nuclear outflows plays an important role in shaping the BEL/NEL regions of these quasars and possibly, regulating the growth of central supermassive black holes (SMBHs). We propose that BEL blueshift and BAL could be different manifestations of the same outflow system viewed at different sightlines and/or phases.


title: Beyond trace reconstruction: Population recovery from the deletion channel

abstract: \emph{Population recovery} is the problem of learning an unknown distribution over an unknown set of $n$-bit strings, given access to independent draws from the distribution that have been independently corrupted according to some noise channel. Recent work has intensively studied such problems both for the bit-flip and erasure noise channels. We initiate the study of population recovery under the \emph{deletion channel}, in which each bit is independently \emph{deleted} with some fixed probability and the surviving bits are concatenated and transmitted. This is a far more challenging noise model than bit-flip~noise or erasure noise; indeed, even the simplest case in which the population size is 1 (corresponding to a trivial distribution supported on a single string) corresponds to the \emph{trace reconstruction} problem, a challenging problem that has received much recent attention (see e.g.~\cite{DOS17,NP17,PZ17,HPP18,HHP18}). We give algorithms and lower bounds for population recovery under the deletion channel when the population size is some $\ell>1$. As our main sample complexity upper bound, we show that for any $\ell=o(\log n/\log \log n)$, a population of $\ell$ strings from $\{0,1\}^n$ can be learned under deletion channel noise using $\smash{2^{n^{1/2+o(1)}}}$ samples. On the lower bound side, we show that $n^{\Omega(\ell)}$ samples are required to perform population recovery under the deletion channel, for all $\ell \leq n^{1/2-\epsilon}$. Our upper bounds are obtained via a robust multivariate generalization of a polynomial-based analysis, due to Krasikov and Roddity \cite{KR97}, of how the $k$-deck of a bit-string uniquely identifies the string; this is a very different approach from recent algorithms for trace reconstruction (the $\ell=1$ case). Our lower bounds build on moment-matching results of Roos~\cite{Roo00} and Daskalakis and Papadimitriou~\cite{DP15}.


title: Prolific Compositions

abstract: Under what circumstances might every extension of a combinatorial structure contain more copies of another one than the original did? This property, which we call prolificity, holds universally in some cases (e.g., finite linear orders) and only trivially in others (e.g., permutations). Integer compositions, or equivalently layered permutations, provide a middle ground. In that setting, there are prolific compositions for a given pattern if and only if that pattern begins and ends with 1. For each pattern, there is an easily constructed automaton that recognises prolific compositions for that pattern. Some instances where there is a unique minimal prolific composition for a pattern are classified.


title: A Block Alternating Optimization Method for Direction-of-Arrival Estimation with Nested Array

abstract: In this paper, direction-of-arrival estimation using nested array is studied in the framework of sparse signal representation. With the vectorization operator, a new real-valued nonnegative sparse signal recovery model which has a wider virtual array aperture is built. To leverage celebrated compressive sensing algorithms, the continuous parameter space has to be discretized to a number of fixed grid points, which inevitably incurs modeling error caused by off-grid gap. To remedy this issue, a block alternating optimization method is put forth that jointly estimates the sparse signal and refines the locations of grid points. Specifically, inspired by the majorization minimization, the proposed method iteratively minimizes a surrogate function majorizing the given objective function, where only a single block of variables are updated per iteration while the remaining ones are kept fixed. The proposed method features affordable computational complexity, and numerical tests corroborate its superior performance relative to existing alternatives in both overdetermined and underdetermined scenarios.


title: An Introduction to MMPDElab

abstract: This article presents an introduction to MMPDElab, a package written in MATLAB for adaptive mesh movement and adaptive moving mesh P1 finite element solution of second-order partial different equations having continuous solutions in one, two, and three spatial dimensions. MMPDElab uses simplicial meshes.


title: Amplification of mid-infrared lasers via magnetized-plasma coupling

abstract: Plasmas may be used as gain media for amplifying intense lasers, and external magnetic fields may be applied to improve the performance. For mid-infrared lasers, the requisite magnetic field is on megagauss scale, which can already be provided by current technologies. Designing the laser amplifier requires knowing the magnetized three-wave coupling coefficient, which is mapped out systematically in this paper. By numerically evaluating its formula, we demonstrate how the coupling depends on angle of wave propagation, laser polarization, magnetic field strength, plasma temperature, and plasma density in the backscattering geometry. Since the mediation is now provided by magnetized plasma waves, the coupling can differ significantly from unmagnetized Raman and Brillouin scatterings.


title: Direct Fitting of Gaussian Mixture Models

abstract: When fitting Gaussian Mixture Models to 3D geometry, the model is typically fit to point clouds, even when the shapes were obtained as 3D meshes. Here we present a formulation for fitting Gaussian Mixture Models (GMMs) directly to a triangular mesh instead of using points sampled from its surface. Part of this work analyzes a general formulation for evaluating likelihood of geometric objects. This modification enables fitting higher-quality GMMs under a wider range of initialization conditions. Additionally, models obtained from this fitting method are shown to produce an improvement in 3D registration for both meshes and RGB-D frames. This result is general and applicable to arbitrary geometric objects, including representing uncertainty from sensor measurements.


title: Improvisation through Physical Understanding: Using Novel Objects as Tools with Visual Foresight

abstract: Machine learning techniques have enabled robots to learn narrow, yet complex tasks and also perform broad, yet simple skills with a wide variety of objects. However, learning a model that can both perform complex tasks and generalize to previously unseen objects and goals remains a significant challenge. We study this challenge in the context of "improvisational" tool use: a robot is presented with novel objects and a user-specified goal (e.g., sweep some clutter into the dustpan), and must figure out, using only raw image observations, how to accomplish the goal using the available objects as tools. We approach this problem by training a model with both a visual and physical understanding of multi-object interactions, and develop a sampling-based optimizer that can leverage these interactions to accomplish tasks. We do so by combining diverse demonstration data with self-supervised interaction data, aiming to leverage the interaction data to build generalizable models and the demonstration data to guide the model-based RL planner to solve complex tasks. Our experiments show that our approach can solve a variety of complex tool use tasks from raw pixel inputs, outperforming both imitation learning and self-supervised learning individually. Furthermore, we show that the robot can perceive and use novel objects as tools, including objects that are not conventional tools, while also choosing dynamically to use or not use tools depending on whether or not they are required.


title: On the Polarized Nonthermal Emission from AR Scorpii

abstract: We study linear polarization of optical emission from white dwarf (WD) binary system AR~Scorpii. The optical emission from this binary is modulating with the beat frequency of the system, and it is highly polarized, with the degree of the polarization reaching $\sim 40$\%. The angle of the polarization monotonically increases with the spin phase, and the total swing angle can reach $360^{\circ}$ over one spin phase. It is also observed that the morphology of the pulse profile and the degree of linear polarization evolve with the orbital phase. These polarization properties can constrain the scenario for nonthermal emission from AR Scorpii. In this paper, we study the polarization properties predicted by the emission model, in which (i) the pulsed optical emission is produced by the synchrotron emission from relativistic electrons trapped by magnetic field lines of the WD and (ii) the emission is mainly produced at magnetic mirror points of the electron motion. We find that this model can reproduce the large swing of the polarization angle, provided that the distribution of the initial pitch angle of the electrons that are leaving the M-type star is biased to a smaller angle rather than a uniform distribution. The observed direction of the swing suggests that the Earth viewing angle is less than $90^{\circ}$ measured from the WD spin axis. The current model prefers an Earth viewing angle of $50^{\circ}-60^{\circ}$ and a magnetic inclination angle of $50^{\circ}-60^{\circ}$ (or $120^{\circ}-130^{\circ}$). We discuss that the different contribution of the emission from M-type star to total emission causes a large variation in the pulsed fraction and the degree of the linear polarization along the orbital phase.