title: Review of approximations for the exchange-correlation energy in density-functional theory

abstract: In this chapter, we provide a review of ground-state Kohn-Sham density-functional theory of electronic systems and some of its extensions, we present exact expressions and constraints for the exchange and correlation density functionals, and we discuss the main families of approximations for the exchange-correlation energy: semilocal approximations, single-determinant hybrid approximations, multideterminant hybrid approximations, dispersion-corrected approximations, as well as orbital-dependent exchange-correlation density functionals. The chapter aims at providing both a consistent bird's-eye view of the field and a detailed description of some of the most used approximations. It is intended to be readable by chemists/physicists and applied mathematicians.


title: Critical Slowing Down Near Topological Transitions in Rate-Distortion Problems

abstract: In rate-distortion (RD) problems one seeks reduced representations of a source that meet a target distortion constraint. Such optimal representations undergo topological transitions at some critical rate values, when their cardinality or dimensionality change. We study the convergence time of the Arimoto-Blahut alternating projection algorithms, used to solve such problems, near those critical points, both for the rate-distortion and information bottleneck settings. We argue that they suffer from critical slowing down -- a diverging number of iterations for convergence -- near the critical points. This phenomenon can have theoretical and practical implications for both machine learning and data compression problems.


title: Nonlinearly Stable Flux Reconstruction High-Order Methods in Split Form

abstract: The flux reconstruction (FR) method has gained popularity in the research community as it recovers promising high-order methods through modally filtered correction fields, such as the discontinuous Galerkin method, amongst others, on unstructured grids over complex geometries. Moreover, FR schemes, specifically energy stable FR (ESFR) schemes also known as Vincent-Castonguay-Jameson-Huynh schemes, have proven attractive as they allow for design flexibility as well as stability proofs for the linear advection problem on affine elements. Additionally, split forms have recently seen a resurgence in research activity due to their resultant nonlinear (entropy) stability proofs. This paper derives for the first time nonlinearly stable ESFR schemes in split form that enable nonlinear stability proofs for, uncollocated, modal, ESFR split forms with different volume and surface cubature nodes. The critical enabling technology is applying the splitting to the discrete stiffness operator. This naturally leads to appropriate surface and numerical fluxes, enabling both entropy stability and conservation proofs. When these schemes are recast in strong form, they differ from schemes found in the ESFR literature as the ESFR correction functions are incorporated on the volume integral. Furthermore, numerical experiments are conducted verifying that the new class of proposed ESFR split forms is nonlinearly stable in contrast to the standard split form ESFR approach. Lastly, the new ESFR split form is shown to obtain the correct orders of accuracy.


title: The Effect of Super-spreader Events in Epidemics

abstract: The spread of infectious epidemics is often accelerated by super-spreader events. Understanding their effect is important, particularly in the context of standard epidemiological models, which require estimates for parameters such as $R_0$. In this letter, we show that the effective value of $R_0$ in super-spreader situations is significantly large, of the order of hundreds, suggesting a delta-function-like behavior during the event. Use of a well-mixed room model supports these findings. They elucidate infection kinetic modeling in enclosed environments, which differ from the standard SIR model, and provide expressions for $R_0$ in terms of physical and operational parameters. The overall impact of super-spreader events can be significant, depending on the state of the epidemic and how the infections generated by the event subsequently spread in the community.


title: Self-play Learning Strategies for Resource Assignment in Open-RAN Networks

abstract: Open Radio Access Network (ORAN) is being developed with an aim to democratise access and lower the cost of future mobile data networks, supporting network services with various QoS requirements, such as massive IoT and URLLC. In ORAN, network functionality is dis-aggregated into remote units (RUs), distributed units (DUs) and central units (CUs), which allows flexible software on Commercial-Off-The-Shelf (COTS) deployments. Furthermore, the mapping of variable RU requirements to local mobile edge computing centres for future centralized processing would significantly reduce the power consumption in cellular networks. In this paper, we study the RU-DU resource assignment problem in an ORAN system, modelled as a 2D bin packing problem. A deep reinforcement learning-based self-play approach is proposed to achieve efficient RU-DU resource management, with AlphaGo Zero inspired neural Monte-Carlo Tree Search (MCTS). Experiments on representative 2D bin packing environment and real sites data show that the self-play learning strategy achieves intelligent RU-DU resource assignment for different network conditions.


title: Successor Feature Sets: Generalizing Successor Representations Across Policies

abstract: Successor-style representations have many advantages for reinforcement learning: for example, they can help an agent generalize from past experience to new goals, and they have been proposed as explanations of behavioral and neural data from human and animal learners. They also form a natural bridge between model-based and model-free RL methods: like the former they make predictions about future experiences, and like the latter they allow efficient prediction of total discounted rewards. However, successor-style representations are not optimized to generalize across policies: typically, we maintain a limited-length list of policies, and share information among them by representation learning or GPI. Successor-style representations also typically make no provision for gathering information or reasoning about latent variables. To address these limitations, we bring together ideas from predictive state representations, belief space value iteration, successor features, and convex analysis: we develop a new, general successor-style representation, together with a Bellman equation that connects multiple sources of information within this representation, including different latent states, policies, and reward functions. The new representation is highly expressive: for example, it lets us efficiently read off an optimal policy for a new reward function, or a policy that imitates a new demonstration. For this paper, we focus on exact computation of the new representation in small, known environments, since even this restricted setting offers plenty of interesting questions. Our implementation does not scale to large, unknown environments -- nor would we expect it to, since it generalizes POMDP value iteration, which is difficult to scale. However, we believe that future work will allow us to extend our ideas to approximate reasoning in large, unknown environments.


title: Experimental Body-input Three-stage DC offset Calibration Scheme for Memristive Crossbar

abstract: Reading several ReRAMs simultaneously in a neuromorphic circuit increases power consumption and limits scalability. Applying small inference read pulses is a vain attempt when offset voltages of the read-out circuit are decisively more. This paper presents an experimental validation of a three-stage calibration scheme to calibrate the DC offset voltage across the rows of the memristive crossbar. The proposed method is based on biasing the body terminal of one of the differential pair MOSFETs of the buffer through a series of cascaded resistor banks arranged in three stages: coarse, fine and finer stages. The circuit is designed in a 130 nm CMOS technology, where the OxRAM-based binary memristors are built on top of it. A dedicated PCB and other auxiliary boards have been designed for testing the chip. Experimental results validate the presented approach, which is only limited by mismatch and electrical noise.


title: A derivation of the Doppler factor in the Lienard-Wiechert potentials

abstract: We present an elegant derivation of the Doppler factor in the Lienard-Wiechert potentials of an extended particle, based on a theorem authored by J. L. Synge and dealing with cross-sections through worldtubes.


title: On the optimal controllability time for linear hyperbolic systems with time-dependent coefficients

abstract: The optimal time for the controllability of linear hyperbolic systems in one dimensional space with one-side controls has been obtained recently for time-independent coefficients in our previous works. In this paper, we consider linear hyperbolic systems with time-varying zero-order terms. We show the possibility that the optimal time for the null-controllability becomes significantly larger than the one of the time-invariant setting even when the zero-order term is indefinitely differentiable. When the analyticity with respect to time is imposed for the zero-order term, we also establish that the optimal time is the same as in the time-independent setting.


title: A Robust Adversarial Network-Based End-to-End Communications System With Strong Generalization Ability Against Adversarial Attacks

abstract: We propose a novel defensive mechanism based on a generative adversarial network (GAN) framework to defend against adversarial attacks in end-to-end communications systems. Specifically, we utilize a generative network to model a powerful adversary and enable the end-to-end communications system to combat the generative attack network via a minimax game. We show that the proposed system not only works well against white-box and black-box adversarial attacks but also possesses excellent generalization capabilities to maintain good performance under no attacks. We also show that our GAN-based end-to-end system outperforms the conventional communications system and the end-to-end communications system with/without adversarial training.


title: Ray-Wave Duality of Electromagnetic Fields: A Feynman Path integral Approach to Classical Vectorial Imaging

abstract: We consider how vectorial aspects (polarization) of light propagation can be implemented, and its origin, within a Feynman path integral approach. A key part of this scheme is in generalising the standard optical path length integral from a scalar to a matrix quantity. Reparametrization invariance along the rays allows a covariant formulation where propagation can take place along a general curve. A general gradient index background is used to demonstrate the scheme. This affords a description of classical imaging optics when the polarization aspects may be varying rapidly and cannot be neglected.


title: Global well-posedness for the one-phase Muskat problem

abstract: The free boundary problem for a two-dimensional fluid filtered in porous media is studied. This is known as the one-phase Muskat problem and is mathematically equivalent to the vertical Hele-Shaw problem driven by gravity force. We prove that if the initial free boundary is the graph of a periodic Lipschitz function, then there exists a global-in-time Lipschitz solution in the strong $L^\infty_t L^2_x$ sense and it is the unique viscosity solution. The proof requires quantitative estimates for layer potentials and pointwise elliptic regularity in Lipschitz domains. This is the first construction of unique global strong solutions for the Muskat problem with initial data of arbitrary size.


title: Numerical investigation of some reductions for the Gatenby-Gawlinski model

abstract: The Gatenby-Gawlinski model for cancer invasion is object of analysis in order to investigate the mathematical framework behind the model working by means of suitable reductions. We perform numerical simulations to study the sharpness/smoothness of the traveling fronts starting from a brief overview about the full model and proceed by examining the case of a two-equations-based and one-equation-based reduction. We exploit a numerical strategy depending on a finite volume approximation and employ a space-averaged wave speed estimate to quantitatively approach the traveling waves phenomenon. Concerning the one equation-based model, we propose a reduction framed within the degenerate reaction-diffusion equations field, which proves to be effective in order to qualitatively recover the typical trends arising from the Gatenby-Gawlinski model. Finally, we carry out some numerical tests in a specific case where the analytical solution is available.


title: The Importance of Funding Space-Based Research

abstract: When it comes to conversations about funding, the questions of whether the United States should be spending its resources on space-based research often rears its head. Opponents of the idea tend to share the opinion that the resources would be better spent helping citizens on the ground. With an estimated homeless population around 562,000 throughout the country, roughly 39.4 million Americans (12.3% of the population) living below the poverty level, and 63.1 million tons of food waste per year, it's hard to argue that the United States does not have its share of problems that need to be addressed. However, a history of space-based research has proven time and time again to bring forth advances in technology and scientific understanding that benefit humans across the globe and provide crucial protection for life on Earth.


title: Stochastic optimization for numerical evaluation of imprecise probabilities

abstract: In applications of imprecise probability, analysts must compute lower (or upper) expectations, defined as the infimum of an expectation over a set of parameter values. Monte Carlo methods consistently approximate expectations at fixed parameter values, but can be costly to implement in grid search to locate minima over large subsets of the parameter space. We investigate the use of stochastic iterative root-finding methods for efficiently computing lower expectations. In two examples we illustrate the use of various stochastic approximation methods, and demonstrate their superior performance in comparison to grid search.


title: Elastocaloric effect in amorphous polymer networks undergoing mechanotropic phase transitions

abstract: Deformations of amorphous polymer networks prepared with significant concentrations of liquid crystalline mesogens have been recently reported to undergo mechanotropic phase transitions. Here, we report that these mechanotropic phase transitions are accompanied by an elastocaloric response ($\Delta T = 2.9 \text{ K}$). Applied uniaxial strain to the elastomeric polymer network transitions the organization of the material from a disordered, amorphous state (order parameter $Q=0$) to the nematic phase ($Q=0.47$). Both the magnitude of the elastocaloric temperature change and mechanically induced order parameter are dependent on the concentration of liquid crystal mesogens in the material. While the observed temperature changes in these materials are smaller than those observed in shape memory alloys, the responsivity, defined as the temperature change divided by the input stress, is larger by an order of magnitude.


title: Large-amplitude prominence oscillations following the impact by a coronal jet

abstract: Observational evidence shows that coronal jets can hit prominences and set them in motion. The impact leads to large-amplitude oscillations (LAOs) of the prominence. In this paper we attempt to understand this process via 2.5D MHD numerical experiments. In our model, the jets are generated in a sheared magnetic arcade above a parasitic bipolar region located in one of the footpoints of the filament channel (FC) supporting the prominence. The shear is imposed with velocities not far above observed photospheric values; it leads to a multiple reconnection process, as in previous jet models. Both a fast Alfv\'enic perturbation and a slower supersonic front preceding a plasma jet are issued from the reconnection site; in the later phase, a more violent (eruptive) jet is produced. The perturbation and jets run along the FC; they are partially reflected at the prominence and partially transmitted through it. There results a pattern of counter-streaming flows along the FC and oscillations of the prominence. The oscillations are LAOs (with amplitude above $10~\mathrm{km\, s^{-1}}$) in parts of the prominence both in the longitudinal and transverse directions. In some field lines, the impact is so strong that the prominence mass is brought out of the dip and down to the chromosphere along the FC. Two cases are studied with different heights of the arcade above the parasitic bipolar region, leading to different heights for the region of the prominence perturbed by the jets. The obtained oscillation amplitudes and periods are in general agreement with the observations.


title: Deep Clustering by Semantic Contrastive Learning

abstract: Whilst contrastive learning has recently brought notable benefits to deep clustering of unlabelled images by learning sample-specific discriminative visual features, its potential for explicitly inferring class decision boundaries is less well understood. This is because its instance discrimination strategy is not class sensitive, therefore, the clusters derived on the resulting sample-specific feature space are not optimised for corresponding to meaningful class decision boundaries. In this work, we solve this problem by introducing Semantic Contrastive Learning (SCL). SCL imposes explicitly distance-based cluster structures on unlabelled training data by formulating a semantic (cluster-aware) contrastive learning objective. Moreover, we introduce a clustering consistency condition to be satisfied jointly by both instance visual similarities and cluster decision boundaries, and concurrently optimising both to reason about the hypotheses of semantic ground-truth classes (unknown/unlabelled) on-the-fly by their consensus. This semantic contrastive learning approach to discovering unknown class decision boundaries has considerable advantages to unsupervised learning of object recognition tasks. Extensive experiments show that SCL outperforms state-of-the-art contrastive learning and deep clustering methods on six object recognition benchmarks, especially on the more challenging finer-grained and larger datasets.


title: Stability of Neural Networks on Riemannian Manifolds

abstract: Convolutional Neural Networks (CNNs) have been applied to data with underlying non-Euclidean structures and have achieved impressive successes. This brings the stability analysis of CNNs on non-Euclidean domains into notice because CNNs have been proved stable on Euclidean domains. This paper focuses on the stability of CNNs on Riemannian manifolds. By taking the Laplace-Beltrami operators into consideration, we construct an $\alpha$-frequency difference threshold filter to help separate the spectrum of the operator with an infinite dimensionality. We further construct a manifold neural network architecture with these filters. We prove that both the manifold filters and neural networks are stable under absolute perturbations to the operators. The results also implicate a trade-off between the stability and discriminability of manifold neural networks. Finally we verify our conclusions with numerical experiments in a wireless adhoc network scenario.


title: Waveguide lattice based architecture for multichannel optical transformations

abstract: We consider coupled waveguide lattices as an architecture that implement a wide range of multiport transformations. In this architecture, a particular transfer matrix is obtained through setting the step-wise profiles of the propagation constants seen by the field evolving in the lattice. To investigate the transformation capabilities, the implementation of a set of transfer matrices taken at random and particular cases of discrete Fourier transform, Hadamard and permutation matrices have been described. Because the waveguide lattices schemes are more compact than their traditional lumped-parameter counterparts, our architecture may be beneficial for using in photonic information processing systems of the future.


title: The Anatomy of an Unusual Edge-on Protoplanetary Disk I. Dust Settling in a Cold Disk

abstract: As the earliest stage of planet formation, massive, optically thick, and gas rich protoplanetary disks provide key insights into the physics of star and planet formation. When viewed edge-on, high resolution images offer a unique opportunity to study both the radial and vertical structures of these disks and relate this to vertical settling, radial drift, grain growth, and changes in the midplane temperatures. In this work, we present multi-epoch HST and Keck scattered light images, and an ALMA 1.3 mm continuum map for the remarkably flat edge-on protoplanetary disk SSTC2DJ163131.2-242627, a young solar-type star in $\rho$ Ophiuchus. We model the 0.8 $\mu$m and 1.3 mm images in separate MCMC runs to investigate the geometry and dust properties of the disk using the MCFOST radiative transfer code. In scattered light, we are sensitive to the smaller dust grains in the surface layers of the disk, while the sub-millimeter dust continuum observations probe larger grains closer to the disk midplane. An MCMC run combining both datasets using a covariance-based log-likelihood estimation was marginally successful, implying insufficient complexity in our disk model. The disk is well characterized by a flared disk model with an exponentially tapered outer edge viewed nearly edge-on, though some degree of dust settling is required to reproduce the vertically thin profile and lack of apparent flaring. A colder than expected disk midplane, evidence for dust settling, and residual radial substructures all point to a more complex radial density profile to be probed with future, higher resolution observations.


title: The Anatomy of an Unusual Edge-on Protoplanetary Disk. II. Gas temperature and a warm outer region

abstract: We present high-resolution $^{12}$CO and $^{13}$CO 2-1 ALMA observations, as well as optical and near-infrared spectroscopy, of the highly-inclined protoplanetary disk around SSTC2D J163131.2-242627. The spectral type we derive for the source is consistent with a $\rm 1.2 \, M_{\odot}$ star inferred from the ALMA observations. Despite its massive circumstellar disk, we find little to no evidence for ongoing accretion on the star. The CO maps reveal a disk that is unusually compact along the vertical direction, consistent with its appearance in scattered light images. The gas disk extends about twice as far away as both the submillimeter continuum and the optical scattered light. CO is detected from two surface layers separated by a midplane region in which CO emission is suppressed, as expected from freeze-out in the cold midplane. We apply a modified version of the Topographically Reconstructed Distribution method presented by Dutrey et al. 2017 to derive the temperature structure of the disk. We find a temperature in the CO-emitting layers and the midplane of $\sim$33 K and $\sim$20 K at $\rm R<200$ au, respectively. Outside of $\rm R>200$ au, the disk's midplane temperature increases to $\sim$30 K, with a nearly vertically isothermal profile. The transition in CO temperature coincides with a dramatic reduction in the sub-micron and sub-millimeter emission from the disk. We interpret this as interstellar UV radiation providing an additional source of heating to the outer part of the disk.


title: Out of Distribution Generalization in Machine Learning

abstract: Machine learning has achieved tremendous success in a variety of domains in recent years. However, a lot of these success stories have been in places where the training and the testing distributions are extremely similar to each other. In everyday situations when models are tested in slightly different data than they were trained on, ML algorithms can fail spectacularly. This research attempts to formally define this problem, what sets of assumptions are reasonable to make in our data and what kind of guarantees we hope to obtain from them. Then, we focus on a certain class of out of distribution problems, their assumptions, and introduce simple algorithms that follow from these assumptions that are able to provide more reliable generalization. A central topic in the thesis is the strong link between discovering the causal structure of the data, finding features that are reliable (when using them to predict) regardless of their context, and out of distribution generalization.


title: Revisiting the Concrete Security of Goldreich's Pseudorandom Generator

abstract: Local pseudorandom generators are a class of fundamental cryptographic primitives having very broad applications in theoretical cryptography. Following Couteau et al.'s work in ASIACRYPT 2018, this paper further studies the concrete security of one important class of local pseudorandom generators, i.e., Goldreich's pseudorandom generators. Our first attack is of the guess-and-determine type. Our result significantly improves the state-of-the-art algorithm proposed by Couteau et al., in terms of both asymptotic and concrete complexity, and breaks all the challenge parameters they proposed. For instance, for a parameter set suggested for 128 bits of security, we could solve the instance faster by a factor of about $2^{61}$, thereby destroying the claimed security completely. Our second attack further exploits the extremely sparse structure of the predicate $P_5$ and combines ideas from iterative decoding. This novel attack, named guess-and-decode, substantially improves the guess-and-determine approaches for cryptographic-relevant parameters. All the challenge parameter sets proposed in Couteau et al.'s work in ASIACRYPT 2018 aiming for 80-bit (128-bit) security levels can be solved in about $2^{58}$ ($2^{78}$) operations. We suggest new parameters for achieving 80-bit (128-bit) security with respect to our attacks. We also extend the attack to other promising predicates and investigate their resistance.


title: An Empirical Analysis of UI-based Flaky Tests

abstract: Flaky tests have gained attention from the research community in recent years and with good reason. These tests lead to wasted time and resources, and they reduce the reliability of the test suites and build systems they affect. However, most of the existing work on flaky tests focus exclusively on traditional unit tests. This work ignores UI tests that have larger input spaces and more diverse running conditions than traditional unit tests. In addition, UI tests tend to be more complex and resource-heavy, making them unsuited for detection techniques involving rerunning test suites multiple times. In this paper, we perform a study on flaky UI tests. We analyze 235 flaky UI test samples found in 62 projects from both web and Android environments. We identify the common underlying root causes of flakiness in the UI tests, the strategies used to manifest the flaky behavior, and the fixing strategies used to remedy flaky UI tests. The findings made in this work can provide a foundation for the development of detection and prevention techniques for flakiness arising in UI tests.


title: Dynamical inverse seesaw mechanism as a simple benchmark for electroweak breaking and Higgs boson studies

abstract: The Standard Model(SM) vacuum is unstable for the measured values of the top Yukawa coupling and Higgs mass. Here we study the issue of vacuum stability when neutrino masses are generated through spontaneous low-scale lepton number violation. In the simplest dynamical inverse seesaw, the SM Higgs has two siblings: a massive CP-even scalar plus a massless Nambu-Goldstone boson, called majoron. For TeV scale breaking of lepton number, Higgs bosons can have a sizeable decay into the invisible majorons. We examine the interplay and complementarity of vacuum stability and perturbativity restrictions, with collider constraints on visible and invisible Higgs boson decay channels. This simple framework may help guiding further studies, for example, at the proposed FCC facility.


title: Estimating the impact of mixed QCD-electroweak corrections on the $W$-mass determination at the LHC

abstract: We study the impact of the recently computed mixed QCD-electroweak corrections to the production of $W$ and $Z$ bosons at the LHC on the value of the $W$ mass extracted from the transverse momentum distribution of charged leptons from $W$ decays. Using the average lepton transverse momenta in $W$ and $Z$ decays as simplified observables for the determination of the $W$ mass, we estimate that mixed QCD-electroweak corrections can shift the extracted value of the $W$ mass by up to ${\cal O}(20)~{\rm MeV}$, depending on the kinematic cuts employed to define fiducial cross sections for $Z$ and $W$ production. Since the target precision of the $W$-mass measurement at the LHC is ${\cal O}(10)~{\rm MeV}$, our results emphasize the need for fully-differential computations of mixed QCD-electroweak corrections and a careful analysis of their potential impact on the determination of the $W$ mass.


title: Processing Dynamics of 3D-Printed Carbon Nanotubes-Epoxy Composites

abstract: Carbon Nanotubes (CNTs)-polymer composites are promising candidates for a myriad of applications. Ad-hoc CNTs-polymer composite fabrication techniques inherently pose roadblock to optimized processing resulting in microstructural defects i.e., void formation, poor interfacial adhesion, wettability, and agglomeration of CNTs inside the polymer matrix. Although improvement in the microstructures can be achieved via additional processing steps such as-mechanical methods and/or chemical functionalization, the resulting composites are somewhat limited in structural and functional performances. Here, we demonstrate that 3D printing technique like-direct ink writing offers improved processing of CNTs-polymer composites. The shear-induced flow of an engineered nanocomposite ink through the micronozzle offers some benefits including reducing the number of voids within the epoxy, improving CNTs dispersion and adhesion with epoxy, and partially aligns the CNTs. Such microstructural changes result in superior mechanical performance and heat transfer in the composites compared to their mold-casted counterparts. This work demonstrates the advantages of 3D printing over traditional fabrication methods, beyond the ability to rapidly fabricate complex architectures, to achieve improved processing dynamics for fabricating CNT-polymer nanocomposites with better structural and functional properties.


title: Closed-formed ab initio solutions of geometric albedos and reflected light phase curves of exoplanets

abstract: Studying the albedos of the planets and moons of the Solar System dates back at least a century. Of particular interest is the relationship between the albedo measured at superior conjunction, known as the ``geometric albedo", and the albedo considered over all orbital phase angles, known as the ``spherical albedo". Determining the relationship between the geometric and spherical albedos usually involves complex numerical calculations and closed-form solutions are restricted to simple reflection laws. Here we report the discovery of closed-form solutions for the geometric albedo and integral phase function, which apply to any law of reflection that only depends on the scattering angle. The shape of a reflected light phase curve, quantified by the integral phase function, and the secondary eclipse depth, quantified by the geometric albedo, may now be self-consistently inverted to retrieve globally averaged physical parameters. Fully Bayesian phase curve inversions for reflectance maps and simultaneous light curve detrending may now be performed due to the efficiency of computation. Demonstrating these innovations for the hot Jupiter Kepler-7b, we infer a geometric albedo of $0.25^{+0.01}_{-0.02}$, a phase integral of $1.77 \pm0.07$, a spherical albedo of $0.44^{+0.02}_{-0.03}$ and a scattering asymmetry factor of $0.07^{+0.12}_{-0.11}$.


title: Sorkin parameter for type-I spontaneous parametric down-conversion biphotons and matter waves

abstract: We propose experimental measurements of the logarithmic negativity, which quantifies quantum correlations using Gouy phase measurements in an asymmetric double-slit interference experiment for twin photons. This is possible because both quantities have analogous dependence with the spatial confinement by the slits and enables one to manipulate the portion of entanglement by the Gouy phase. In order to obtain those measurements, we need to work in a regime where the position correlations between particles are strong, therefore we investigate such correlations for biphotons. Since we would like to handle entanglement quantifiers through the Gouy phase, we analyze the Gouy phase difference for two entangled photons in an asymmetric double-slit interference experiment.


title: Solitary waves in a Whitham equation with small surface tension

abstract: Using a nonlocal version of the center manifold theorem and a normal form reduction, we prove the existence of small-amplitude generalized solitary-wave solutions and modulated solitary-wave solutions to the steady gravity-capillary Whitham equation with weak surface tension. Through the application of the center manifold theorem, the nonlocal equation for the solitary wave profiles is reduced to a four-dimensional system of ODEs inheriting reversibility. Along particular parameter curves, relating directly to the classical gravity-capillary water wave problem, the associated linear operator is seen to undergo either a reversible $0^{2+}(i k_0)$ bifurcation or a reversible $(i s)^2$ bifurcation. Through a normal form transformation, the reduced system of ODEs along each relevant parameter curve is seen to be well approximated by a truncated system retaining only second-order or third-order terms. These truncated systems relate directly to systems obtained in the study of the full gravity-capillary water wave equation and, as such, the existence of generalized and modulated solitary waves for the truncated systems is guaranteed by classical works, and they are readily seen to persist as solutions of the gravity-capillary Whitham equation due to reversibility. Consequently, this work illuminates further connections between the gravity-capillary Whitham equation and the full two-dimensional gravity-capillary water wave problem.


title: Efficient UAV Trajectory-Planning using Economic Reinforcement Learning

abstract: Advances in unmanned aerial vehicle (UAV) design have opened up applications as varied as surveillance, firefighting, cellular networks, and delivery applications. Additionally, due to decreases in cost, systems employing fleets of UAVs have become popular. The uniqueness of UAVs in systems creates a novel set of trajectory or path planning and coordination problems. Environments include many more points of interest (POIs) than UAVs, with obstacles and no-fly zones. We introduce REPlanner, a novel multi-agent reinforcement learning algorithm inspired by economic transactions to distribute tasks between UAVs. This system revolves around an economic theory, in particular an auction mechanism where UAVs trade assigned POIs. We formulate the path planning problem as a multi-agent economic game, where agents can cooperate and compete for resources. We then translate the problem into a Partially Observable Markov decision process (POMDP), which is solved using a reinforcement learning (RL) model deployed on each agent. As the system computes task distributions via UAV cooperation, it is highly resilient to any change in the swarm size. Our proposed network and economic game architecture can effectively coordinate the swarm as an emergent phenomenon while maintaining the swarm's operation. Evaluation results prove that REPlanner efficiently outperforms conventional RL-based trajectory search.


title: Online adaptive algorithm for Constraint Energy Minimizing Generalized Multiscale Discontinuous Galerkin Method

abstract: In this research, we propose an online basis enrichment strategy within the framework of a recently developed constraint energy minimizing generalized multiscale discontinuous Galerkin method (CEM-GMsDGM). Combining the technique of oversampling, one makes use of the information of the current residuals to adaptively construct basis functions in the online stage to reduce the error of multiscale approximation. A complete analysis of the method is presented, which shows the proposed online enrichment leads to a fast convergence from multiscale approximation to the fine-scale solution. The error reduction can be made sufficiently large by suitably selecting oversampling regions and the number of oversampling layers. Further, the convergence rate of the enrichment algorithm depends on a factor of exponential decay regarding the number of oversampling layers and a user-defined parameter. Numerical results are provided to demonstrate the effectiveness and efficiency of the proposed online adaptive algorithm.


title: Analysis of a dynamical system modeling lasers and applications for optical neural networks

abstract: An analytical study of dynamical properties of a semiconductor laser with optical injection of arbitrary polarization is presented. It is shown that if the injected field is sufficiently weak, then the laser has nine equilibrium points, however, only one of them is stable. Even if the injected field is linearly polarized, six of the equilibrium points have a state of polarization that is elliptical. Dependence of the equilibrium points on the injected field is described, and it is shown that as the intensity of the injected field increases, the number of equilibrium points decreases, with only a single equilibrium point remaining for strong enough injected fields. As an application, a complex-valued optical neural network with working principle based on injection locking is proposed.


title: New plasma regimes with small ELMs and high confinement at the Joint European Torus

abstract: New plasma regimes with high confinement, low core impurity accumulation and small Edge localized mode (ELMs) perturbations have been obtained close to ITER conditions in magnetically confined plasmas from the Joint European torus (JET) tokamak. Such regimes are achieved by means of optimized particle fuelling conditions which trigger a self-organize state with a strong increase in rotation and ion temperature and a decrease of the edge density. An interplay between core and edge plasma regions leads to reduced turbulence levels and outward impurity convection. These results pave the way to an attractive alternative to the standard plasmas considered for fusion energy generation in a tokamak with metallic wall environment such as the ones expected in ITER


title: Weighted-Graph-Based Change Point Detection

abstract: We consider the detection and localization of change points in the distribution of an offline sequence of observations. Based on a nonparametric framework that uses a similarity graph among observations, we propose new test statistics when at most one change point occurs and generalize them to multiple change points settings. The proposed statistics leverage edge weight information in the graphs, exhibiting substantial improvements in testing power and localization accuracy in simulations. We derive the null limiting distribution, provide accurate analytic approximations to control type I error, and establish theoretical guarantees on the power consistency under contiguous alternatives for the one change point setting, as well as the minimax localization rate. In the multiple change points setting, the asymptotic correctness of the number and location of change points are also guaranteed. The methods are illustrated on the MIT proximity network data.


title: The Proximity Operator of the Log-Sum Penalty

abstract: The log-sum penalty is often adopted as a replacement for the $\ell_0$ pseudo-norm in compressive sensing and low-rank optimization. The hard-thresholding operator, i.e., the proximity operator of the $\ell_0$ penalty, plays an essential role in applications; similarly, we require an efficient method for evaluating the proximity operator of the log-sum penalty. Due to the nonconvexity of this function, its proximity operator is commonly computed through the iteratively reweighted $\ell_1$ method, which replaces the log-sum term with its first-order approximation. This paper reports that the proximity operator of the log-sum penalty actually has an explicit expression. With it, we show that the iteratively reweighted $\ell_1$ solution disagrees with the true proximity operator of the log-sum penalty in certain regions. As a by-product, the iteratively reweighted $\ell_1$ solution is precisely characterized in terms of the chosen initialization. We also give the explicit form of the proximity operator for the composition of the log-sum penalty with the singular value function, as seen in low-rank applications. These results should be useful in the development of efficient and accurate algorithms for optimization problems involving the log-sum penalty.


title: $W^+W^-H$ production through bottom quarks fusion at hadron colliders

abstract: With the standard model working well in describing the collider data, the focus is now on determining the standard model parameters as well as for any hint of deviation. In particular, the determination of the couplings of the Higgs boson with itself and with other particles of the model is important to better understand the electroweak symmetry breaking sector of the model. In this letter, we look at the process $pp \to WWH$, in particular through the fusion of bottom quarks. Due to the non-negligible coupling of the Higgs boson with the bottom quarks, there is a dependence on the $WWHH$ coupling in this process. This sub-process receives largest contribution when the $W$ bosons are longitudinally polarized. We compute one-loop QCD corrections to various final states with polarized $W$ bosons. We find that the corrections to the final state with the longitudinally polarized $W$ bosons are large. It is shown that the measurement of the polarization of the $W$ bosons can be used as a tool to probe the $WWHH$ coupling in this process. We also examine the effect of varying $WWHH$ coupling in the $\kappa$-framework.


title: Preventing Unauthorized Use of Proprietary Data: Poisoning for Secure Dataset Release

abstract: Large organizations such as social media companies continually release data, for example user images. At the same time, these organizations leverage their massive corpora of released data to train proprietary models that give them an edge over their competitors. These two behaviors can be in conflict as an organization wants to prevent competitors from using their own data to replicate the performance of their proprietary models. We solve this problem by developing a data poisoning method by which publicly released data can be minimally modified to prevent others from train-ing models on it. Moreover, our method can be used in an online fashion so that companies can protect their data in real time as they release it.We demonstrate the success of our approach onImageNet classification and on facial recognition.


title: Gauge-underdetermination and shades of locality in the Aharonov-Bohm effect

abstract: I address the view that the classical electromagnetic potentials are shown by the Aharonov-Bohm effect to be physically real (which I dub: 'the potentials view'). I give a historico-philosophical presentation of this view and assess its prospects, more precisely than has so far been done in the literature. Taking the potential as physically real runs prima facie into 'gauge-underdetermination': different gauge choices represent different physical states of affairs and hence different theories. I then illustrate this theme by what I take to be the basic insight of the AB effect for the potentials view, namely that the gauge equivalence class that directly corresponds to the electric and magnetic fields (which I call the Wide Equivalence Class) is too wide, i.e., the Narrow Equivalence Class encodes additional physical degrees of freedom: these only play a distinct role in a multiply-connected space. There is a trade-off between explanatory power and gauge symmetries. Although this narrower equivalence class gives a explanation in terms of `local interactions', locality is not satisfied in the sense of signal locality. It is therefore intellectually mandatory to seek desiderata that will distinguish even within these narrower equivalence classes, i.e. will prefer some elements of such an equivalence class over others. I consider various formulations of locality, such as Bell locality, local interaction Hamiltonians, and signal locality. I show that Bell locality can only be evaluated if one fixes the gauge freedom completely. Yet, an explanation in terms of signal locality can be accommodated by the Lorenz gauge: the potentials propagate in waves at finite speed. I therefore suggest the Lorenz gauge potentials theory -- an even narrower gauge equivalence relation -- as the ontology of electrodynamics.


title: A transiting warm giant planet around the young active star TOI-201

abstract: We present the confirmation of the eccentric warm giant planet TOI-201 b, first identified as a candidate in \textit{TESS} photometry (Sectors 1-8, 10-13, and 27-28) and confirmed using ground-based photometry from NGTS and radial velocities from FEROS, HARPS, CORALIE, and \textsc{Minerva}-Australis. TOI-201 b orbits a young ($\mathrm{0.87^{+0.46}_{-0.49} \, Gyr}$) and bright(V=9.07 mag) F-type star with a $\mathrm{52.9781 \, d}$ period. The planet has a mass of $\mathrm{0.42^{+0.05}_{-0.03}\, M_J}$, a radius of $\mathrm{1.008^{+0.012}_{-0.015}\, R_J}$, and an orbital eccentricity of $0.28^{+0.06}_{-0.09}$; it appears to still be undergoing fairly rapid cooling, as expected given the youth of the host star. The star also shows long-term variability in both the radial velocities and several activity indicators, which we attribute to stellar activity. The discovery and characterization of warm giant planets such as TOI-201 b is important for constraining formation and evolution theories for giant planets.


title: The studies on phonons and magnons in (Ni80Fe20/Au/Co/Au) multilayers of different number of repetitions

abstract: In Ni80Fe20/Au/Co/Au deposited on silicon substrate the interaction between spin waves and surface acoustic waves is observed by Brillion light scattering spectroscopy. We show that by changing the number of receptions in the multilayer we can tune the dispersion of surface acoustic waves (SAWs) not affecting significantly fundamentals spin wave (SW) mode. As a result, we can control (activate or deactivate) the magnetoelastic interaction between fundamental SW mode and SAWs by the thickness of the magnetostrictive multilayer.


title: Lazy FSCA for Unsupervised Variable Selection

abstract: Various unsupervised greedy selection methods have been proposed as computationally tractable approximations to the NP-hard subset selection problem. These methods rely on sequentially selecting the variables that best improve performance with respect to a selection criterion. Theoretical results exist that provide performance bounds and enable "lazy greedy" efficient implementations for selection criteria that satisfy a diminishing returns property known as submodularity. This has motivated the development of variable selection algorithms based on mutual information and frame potential. Recently, the authors introduced Forward Selection Component Analysis (FSCA) which uses variance explained as its selection criterion. While this criterion is not submodular, FSCA has been shown to be highly effective for applications such as measurement plan optimisation. In this paper a "lazy" implementation of the FSCA algorithm (L-FSCA) is proposed, which, although not equivalent to FSCA due to the absence of submodularity, has the potential to yield comparable performance while being up to an order of magnitude faster to compute. The efficacy of L-FSCA is demonstrated by performing a systematic comparison with FSCA and five other unsupervised variable selection methods from the literature using simulated and real-world case studies. Experimental results confirm that L-FSCA yields almost identical performance to FSCA while reducing computation time by between 22% and 94% for the case studies considered.


title: Influence of gravitational waves upon light. Part II. Electric field propagation and interference pattern in a gravitational wave detector

abstract: In this second article of the series, we apply our recently derived equation for the electric field propagation along light rays [arXiv:2004.03496], valid on the electromagnetic geometrical optics limit, to the special case of a toy interferometer used to detect gravitational waves in a flat background. Such an equation shows that, assuming the detector is in the transverse-traceless frame, which has a local shearing relative motion due to the gravitational wave perturbations, the electric field does not propagate as in an inertial reference frame in Minkowski spacetime. We present the electric field at the end of the interferometric process, for arbitrary arm configurations with respect to the plane gravitational wave packet propagation direction. Then, for normal incidence, we compute the interference pattern and, in addition to the usual term associated with the difference in path traveled by light in the arms, we deduce two new contributions to the final intensity, arising from: (i) the round-trip electromagnetic frequency shift and (ii) the divergence of the light beam. Their quantitative relevance is compared to the traditional contribution and shown to be typically negligible due to the geometrical optics regime of light. Moreover, a non-parallel transport of the polarization vector takes place, in general, because of the gravitational wave, a feature which could generate further contributions. However, we conclude that for the normal incidence case such vector is parallel transported, preventing this kind of correction.


title: Boosting energy-time entanglement using coherent time-delayed feedback

abstract: The visibility of the two-photon interference in the Franson interferometer serves as a measure of the energy-time entanglement of the photons. We propose to control the visibility of the interference in the second-order coherence function by implementing a coherent time-delayed feedback mechanism. Simulating the non-Markovian dynamics within the matrix product state framework, we find that the visibility for two photons emitted from a three-level system (3LS) in ladder configuration can be enhanced significantly for a wide range of parameters by slowing down the decay of the upper level of the 3LS.


title: A comprehensive survey on point cloud registration

abstract: Registration is a transformation estimation problem between two point clouds, which has a unique and critical role in numerous computer vision applications. The developments of optimization-based methods and deep learning methods have improved registration robustness and efficiency. Recently, the combinations of optimization-based and deep learning methods have further improved performance. However, the connections between optimization-based and deep learning methods are still unclear. Moreover, with the recent development of 3D sensors and 3D reconstruction techniques, a new research direction emerges to align cross-source point clouds. This survey conducts a comprehensive survey, including both same-source and cross-source registration methods, and summarize the connections between optimization-based and deep learning methods, to provide further research insight. This survey also builds a new benchmark to evaluate the state-of-the-art registration algorithms in solving cross-source challenges. Besides, this survey summarizes the benchmark data sets and discusses point cloud registration applications across various domains. Finally, this survey proposes potential research directions in this rapidly growing field.


title: Natural Language Understanding for Argumentative Dialogue Systems in the Opinion Building Domain

abstract: This paper introduces a natural language understanding (NLU) framework for argumentative dialogue systems in the information-seeking and opinion building domain. The proposed framework consists of two sub-models, namely intent classifier and argument similarity. Intent classifier model stacks BiLSTM with attention mechanism on top of the pre-trained BERT model and fine-tune the model for recognizing the user intent, whereas the argument similarity model employs BERT+BiLSTM for identifying system arguments the user refers to in his or her natural language utterances. Our model is evaluated in an argumentative dialogue system that engages the user to inform him-/herself about a controversial topic by exploring pro and con arguments and build his/her opinion towards the topic. In order to evaluate the proposed approach, we collect user utterances for the interaction with the respective system labeling intent and referenced argument in an extensive online study. The data collection includes multiple topics and two different user types (native English speakers from the UK and non-native English speakers from China). Additionally, we evaluate the proposed intent classifier and argument similarity models separately on the publicly available Banking77 and STS benchmark datasets. The evaluation indicates a clear advantage of the utilized techniques over baseline approaches on several datasets, as well as the robustness of the proposed approach against new topics and different language proficiency as well as the cultural background of the user. Furthermore, results show that our intent classifier model outperforms DIET, DistillBERT, and BERT fine-tuned models in few-shot setups (i.e., with 10, 20, or 30 labeled examples per intent) and full data setup.


title: Micrometer-thick, atomically random Si0.06Ge0.90Sn0.04 for silicon-integrated infrared optoelectronics

abstract: A true monolithic infrared photonics platform is within reach if strain and bandgap energy can be independently engineered in SiGeSn semiconductors. Herein, we investigate the structural and optoelectronic properties of a 1.5 {\mu}m-thick Si0.06Ge0.90Sn0.04 layer that is nearly lattice-matched to a Ge on Si substrate. Atomic-level studies demonstrate high crystalline quality and uniform composition and show no sign of short-range ordering and clusters. Room temperature spectroscopic ellipsometry and transmission measurements show direct bandgap absorption at 0.83 eV and a reduced indirect bandgap absorption at lower energies. Si0.06Ge0.90Sn0.04 photoconductive devices operating at room temperature exhibit dark current and spectral responsivity (1 A/W below 1.5 {\mu}m wavelengths) similar to Ge on Si devices, with the advantage of a near-infrared band gap tunable by alloy composition. These results underline the relevance of SiGeSn semiconductors in implementing a group IV material platform for silicon-integrated infrared optoelectronics.


title: Eulerian numerical modeling of contaminant transport in Lower Manhattan, New York City, from a point-source release under the dominant wind condition: Insights gained via LES

abstract: Released pollutants and poisonous chemicals in the highly-populated urban areas can spread via wind flow and affect the health, safety, and wellbeing of those close to the source and others downstream. Being able to predict the transport of airborne contaminants in such urban areas could help stakeholders to plan emergency responses. Airflow due to the wind in urban areas with closely packed buildings results in complicated aerodynamics around buildings that, in turn, can lead to intricate contaminant transport phenomena. We conduct high-resolution large-eddy simulation (LES) of airflow in Southern Manhattan, a highly-populated area in New York City, due to a dominant wind blowing from south to north. The LES results show complex flow dynamics around the buildings of various heights and, consequently, complicated contaminant transport processes from a source-point which is located in front of the New York Stock Exchange building at 30 m above ground. The flow tends to align with the streets that are fairly parallel to the wind direction creating high-velocity core regions, while the airflow loses its momentum markedly farther downwind in the urban area. For as long as the source of contamination exists, the pollution plume rapidly travels through the streets reaching downwind areas. When the source is removed, contaminants that are trapped in the low-momentum urban area linger to exit the flow domain.


title: The future lifespan of Earth's oxygenated Atmosphere

abstract: Earth's modern atmosphere is highly oxygenated and is a remotely detectable signal of its surface biosphere. However, the lifespan of oxygen-based biosignatures in Earth's atmosphere remains uncertain, particularly for the distant future. Here we use a combined biogeochemistry and climate model to examine the likely timescale of oxygen-rich atmospheric conditions on Earth. Using a stochastic approach, we find that the mean future lifespan of Earth's atmosphere with oxygen levels more than 1% of the present atmospheric level is 1.08+-0.14 billion years. The model projects that a deoxygenation of the atmosphere, with atmospheric oxygen dropping sharply to levels reminiscnet of the Archaean Earth, will most probably be triggered before the inception of moist greenhouse conditions in Earth's climate system and before the extensive loss of surface water from the atmosphere. We find that future deoxygenation is an inevitable consequence of increasing solar fluxes, whereas its precise timing is modulated by the exchange flux of reducing power between the mantle and the ocean-atmosphere-crust system. Our results suggest that the planetary carbonate-silicate cycle will tend to lead to terminally CO2-limited biospheres and rapid atmospheric deoxygenation, emphasizing the need for robust atmospheric biosignatures applicable to weakly oxygenated and anoxic exoplanet atmospheres and highlighting the potential importance of atmospheric organic haze during the terminal stages of planetary habitability.


title: Shift Invariance Can Reduce Adversarial Robustness

abstract: Shift invariance is a critical property of CNNs that improves performance on classification. However, we show that invariance to circular shifts can also lead to greater sensitivity to adversarial attacks. We first characterize the margin between classes when a shift-invariant linear classifier is used. We show that the margin can only depend on the DC component of the signals. Then, using results about infinitely wide networks, we show that in some simple cases, fully connected and shift-invariant neural networks produce linear decision boundaries. Using this, we prove that shift invariance in neural networks produces adversarial examples for the simple case of two classes, each consisting of a single image with a black or white dot on a gray background. This is more than a curiosity; we show empirically that with real datasets and realistic architectures, shift invariance reduces adversarial robustness. Finally, we describe initial experiments using synthetic data to probe the source of this connection.


title: On the Importance of Sampling in Training GCNs: Tighter Analysis and Variance Reduction

abstract: Graph Convolutional Networks (GCNs) have achieved impressive empirical advancement across a wide variety of semi-supervised node classification tasks. Despite their great success, training GCNs on large graphs suffers from computational and memory issues. A potential path to circumvent these obstacles is sampling-based methods, where at each layer a subset of nodes is sampled. Although recent studies have empirically demonstrated the effectiveness of sampling-based methods, these works lack theoretical convergence guarantees under realistic settings and cannot fully leverage the information of evolving parameters during optimization. In this paper, we describe and analyze a general doubly variance reduction schema that can accelerate any sampling method under the memory budget. The motivating impetus for the proposed schema is a careful analysis of the variance of sampling methods where it is shown that the induced variance can be decomposed into node embedding approximation variance (zeroth-order variance) during forward propagation and layerwise-gradient variance (first-order variance) during backward propagation. We theoretically analyze the convergence of the proposed schema and show that it enjoys an $\mathcal{O}(1/T)$ convergence rate. We complement our theoretical results by integrating the proposed schema in different sampling methods and applying them to different large real-world graphs.


title: New Coalescences for the Painlev\'e Equations

abstract: The Painlev\'e equations are here connected to other classes of equations with the Painlev\'e Property (Ince's equations) by the same degeneracy procedure that connects the Painlev\'e equations (coalescence). These Ince's equations here are also connected among themselves like in the traditional Painlev\'e's coalescence cascade. Such degeneracy is considered also for the special equations, symmetric equations and B\"acklund transformations.


title: A Radiation-Balanced Silica Fiber Amplifier

abstract: We report what we believe to be the first radiation-balanced fiber amplifier - a device that provides optical gain while experiencing no temperature rise. The gain medium is a silica fiber with a 21 um-diameter core highly doped with Yb3+ (2.52 wt.%) and co doped with 2.00 wt.% Al to reduce concentration quenching. The amplifier was core-pumped with 1040 nm light to create anti-Stokes fluorescence (ASF) cooling and gain in the core at 1064 nm. Using a custom slow-light FBG sensor with mK resolution, temperature measurements were performed at multiple locations along the amplifier fiber. A 4.35-m fiber pumped with 2.62 W produced 17 dB of gain while the average fiber temperature remained slightly below room temperature. This advancement is a fundamental step toward the creation of ultra-stable lasers necessary to many applications, especially low-noise sensing and high-precision metrology.


title: Compactness statistics for spanning tree recombination

abstract: Ensemble analysis has become an important tool for quantifying gerrymandering; the main idea is to generate a large, random sample of districting plans (an "ensemble") to which any proposed plan may be compared. If a proposed plan is an extreme outlier compared to the ensemble with regard to various redistricting criteria, this may indicate that the plan was deliberately engineered to produce a specific outcome. Many methods have been used to construct ensembles, and a fundamental question that arises is: Given a method for constructing plans, can we identify a probability distribution on the space of plans that describes the probability of constructing any particular plan by that method? Recently, MCMC methods have become a predominant tool for constructing ensembles. Here we focus on the MCMC method known as "ReCom," which was introduced in 2018 by the MGGG Redistricting Lab. ReCom tends to produce plans with more compact districts than some other methods, and we sought to better understand this phenomenon. We adopted a discrete analog of district perimeter called "cut edges" as a quantitative measure for district compactness; this measure was proposed by Duchin and Tenner, and it avoids some of the difficulties associated with compactness measures based on geographic perimeter, such as the Polsby-Popper score. To model the basic ReCom step, we constructed ensembles of 2-district plans for two grid graphs and for the precinct graph of Boulder County, CO. We found that the probability of sampling any particular plan -- which is roughly proportional to the product of the numbers of spanning trees for each of the two districts -- is also approximately proportional to an exponentially decaying function of the number of cut edges in the plan. This is an important step towards understanding compactness properties for districting plans produced by the ReCom method.


title: Decoding supercodes of Gabidulin codes and applications to cryptanalysis

abstract: This article discusses the decoding of Gabidulin codes and shows how to extend the usual decoder to any supercode of a Gabidulin code at the cost of a significant decrease of the decoding radius. Using this decoder, we provide polynomial time attacks on the rank-metric encryption schemes RAMESSES and LIGA.


title: Exploring the roles of local mobility patterns, socioeconomic conditions, and lockdown policies in shaping the patterns of COVID-19 spread

abstract: The COVID-19 crisis has shown that we can only prevent the risk of mass contagion through timely, large-scale, coordinated, and decisive actions. However, frequently the models used by experts [from whom decision-makers get their main advice] focus on a single perspective [for example, the epidemiological one] and do not consider many of the multiple forces that affect the COVID-19 outbreak patterns. The epidemiological, socioeconomic, and human mobility context of COVID-19 can be considered as a complex adaptive system. So, these interventions (for example, lock-downs) could have many and/or unexpected ramifications. This situation makes it difficult to understand the overall effect produced by any public policy measure and, therefore, to assess its real effectiveness and convenience. By using mobile phone data, socioeconomic data, and COVID-19 cases data recorded throughout the pandemic development, we aim to understand and explain [make sense of] the observed heterogeneous regional patterns of contagion across time and space. We will also consider the causal effects produced by confinement policies by developing data-based models to explore, simulate, and estimate these policies' effectiveness. We intend to develop a methodology to assess and improve public policies' effectiveness associated with the fight against the pandemic, emphasizing its convenience, the precise time of its application, and extension. The contributions of this work can be used regardless of the region. The only likely impediment is the availability of the appropriate data.


title: Robust PDF Files Forensics Using Coding Style

abstract: Identifying how a file has been created is often interesting in security. It can be used by both attackers and defenders. Attackers can exploit this information to tune their attacks and defenders can understand how a malicious file has been created after an incident. In this work, we want to identify how a PDF file has been created. This problem is important because PDF files are extremely popular: many organizations publish PDF files online and malicious PDF files are commonly used by attackers. Our approach to detect which software has been used to produce a PDF file is based on coding style: given patterns that are only created by certain PDF producers. We have analyzed the coding style of 900 PDF files produced using 11 PDF producers on 3 different Operating Systems. We have obtained a set of 192 rules which can be used to identify 11 PDF producers. We have tested our detection tool on 508836 PDF files published on scientific preprints servers. Our tool is able to detect certain producers with an accuracy of 100%. Its overall detection is still high (74%). We were able to apply our tool to identify how online PDF services work and to spot inconsistency.


title: The Spatial Selective Auditory Attention of Cochlear Implant Users in Different Conversational Sound Levels

abstract: In multi speakers environments, cochlear implant (CI) users may attend to a target sound source in a different manner from the normal hearing (NH) individuals during a conversation. This study attempted to investigate the effect of conversational sound levels on the mechanisms adopted by CI and NH listeners in selective auditory attention and how it affects their daily conversation. Nine CI users (five bilateral, three unilateral, and one bimodal) and eight NH listeners participated in this study. The behavioral speech recognition scores were collected using a matrix sentences test and neural tracking to speech envelope was recorded using electroencephalography (EEG). Speech stimuli were presented at three different levels (75, 65, and 55 dB SPL) in the presence of two maskers from three spatially separated speakers. Different combinations of assisted/impaired hearing modes were evaluated for CI users and the outcomes were analyzed in three categories: electric hearing only, acoustic hearing only, and electric+acoustic hearing. Our results showed that increasing the conversational sound level degraded the selective auditory attention in electrical hearing. On the other hand, increasing the sound level improved the selective auditory attention for the acoustic hearing group. In NH listeners, however, increasing the sound level did not cause a significant change in the auditory attention. Our result implies that the effect of the sound level on the selective auditory attention varies depending on hearing modes and the loudness control is necessary for the ease of attending to the conversation by CI users.


title: Tribracket Polynomials

abstract: We introduce a six-variable polynomial invariant of Niebrzydowski tribrackets analogous to quandle,rack and biquandle polynomials. Using the subtribrackets of a tribracket, we additionally define subtribracket polynomials and establish a sufficient condition for isomorphic subtribrackets to have the same polynomial regardless of their embedding in the ambient tribracket. As an application, we enhance the tribracket counting invariant of knots and links using subtribracket polynomials and provide examples to demonstrate that this enhancement is proper.


title: Radion dynamics, heavy Kaluza-Klein resonances and gravitational waves

abstract: We study the confinement/deconfinement phase transition of the radion field in a warped model with a polynomial bulk potential. The backreaction of the radion on the metric is taken into account by using the superpotential formalism, while the radion effective potential is obtained from a novel formulation which can incorporate the backreaction. The phase transition leads to a stochastic gravitational wave background that depends on the energy scale of the first Kaluza-Klein resonance, $m_{\textrm{KK}}$. This work completes previous studies in the following aspects: i) we detail the evaluation of the radion spectrum; ii) we report on the mismatches between the thick wall approximation and the numerical bounce solution; iii) we include a suppression factor in the spectrum of sound waves accounting for their finite lifetime; and, iv) we update the bound on $m_{\textrm{KK}}$ in view of the O3 LIGO and Virgo data. We find that the forthcoming gravitational wave interferometers can probe scenarios where $m_{\textrm{KK}} \lesssim 10^9$ TeV, while the O3-run bounds rule out warped models with $10^4 \textrm{TeV} \lesssim m_{\textrm{KK}} \lesssim 10^7$ TeV exhibiting an extremely strong confinement/deconfinement phase transition.


title: Quantifying the spin mixing conductance of EuO/W heterostructures by spin Hall magnetoresistance experiments

abstract: The spin Hall magnetoresistance (SMR) allows to investigate the magnetic textures of magnetically ordered insulators in heterostructures with normal metals by magnetotransport experiments. We here report the observation of the SMR in in-situ prepared ferromagnetic EuO/W thin film bilayers with magnetically and chemically well-defined interfaces. We characterize the magnetoresistance effects utilizing angle-dependent and field-dependent magnetotransport measurements as a function of temperature. Applying the established SMR model, we derive and quantify the real and imaginary parts of the complex spin mixing interface conductance. We find that the imaginary part is by one order of magnitude larger than the real part. Both decrease with increasing temperature. This reduction is in agreement with thermal fluctuations in the ferromagnet.


title: Exploitation and Sanitization of Hidden Data in PDF Files

abstract: Organizations publish and share more and more electronic documents like PDF files. Unfortunately, most organizations are unaware that these documents can compromise sensitive information like authors names, details on the information system and architecture. All these information can be exploited easily by attackers to footprint and later attack an organization. In this paper, we analyze hidden data found in the PDF files published by an organization. We gathered a corpus of 39664 PDF files published by 75 security agencies from 47 countries. We have been able to measure the quality and quantity of information exposed in these PDF files. It can be effectively used to find weak links in an organization: employees who are running outdated software. We have also measured the adoption of PDF files sanitization by security agencies. We identified only 7 security agencies which sanitize few of their PDF files before publishing. Unfortunately, we were still able to find sensitive information within 65% of these sanitized PDF files. Some agencies are using weak sanitization techniques: it requires to remove all the hidden sensitive information from the file and not just to remove the data at the surface. Security agencies need to change their sanitization methods.


title: Search for resonant and nonresonant new phenomena in high-mass dilepton final states at $\sqrt{s} = $ 13 TeV

abstract: A search is presented for physics beyond the standard model (SM) using electron or muon pairs with high invariant mass. A data set of proton-proton collisions collected by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV from 2016 to 2018 corresponding to a total integrated luminosity of up to 140 fb$^{-1}$ is analyzed. No significant deviation is observed with respect to the SM background expectations. Upper limits are presented on the ratio of the product of the production cross section and the branching fraction to dileptons of a new narrow resonance to that of the Z boson. These provide the most stringent lower limits to date on the masses for various spin-1 particles, spin-2 gravitons in the Randall--Sundrum model, as well as spin-1 mediators between the SM and dark matter particles. Lower limits on the ultraviolet cutoff parameter are set both for four-fermion contact interactions and for the Arkani-Hamed, Dimopoulos, and Dvali model with large extra dimensions. Lepton flavor universality is tested at the TeV scale for the first time by comparing the dimuon and dielectron mass spectra. No significant deviation from the SM expectation of unity is observed.


title: A Gaussian Process Regression Reveals No Evidence for Planets Orbiting Kapteyn's Star

abstract: Radial-velocity (RV) planet searches are often polluted by signals caused by gas motion at the star's surface. Stellar activity can mimic or mask changes in the RVs caused by orbiting planets, resulting in false positives or missed detections. Here we use Gaussian process (GP) regression to disentangle the contradictory reports of planets vs. rotation artifacts in Kapteyn's star (Anglada-Escude et al. 2014, Robertson et al. 2015, Anglada-Escude et al. 2016). To model rotation, we use joint quasi-periodic kernels for the RV and H-alpha signals, requiring that their periods and correlation timescales be the same. We find that the rotation period of Kapteyn's star is 125 days, while the characteristic active-region lifetime is 694 days. Adding a planet to the RV model produces a best-fit orbital period of 100~years, or 10 times the observing time baseline, indicating that the observed RVs are best explained by star rotation only. We also find no significant periodic signals in residual RV data sets constructed by subtracting off realizations of the best-fit rotation model and conclude that both previously reported "planets" are artifacts of the star's rotation and activity. Our results highlight the pitfalls of using sinusoids to model quasi-periodic rotation signals.


title: G-Family Polynomials

abstract: We introduce two notions of quandle polynomials for G-families of quandles: the quandle polynomial of the associated quandle and a G-family polynomial with coefficients in the group ring of G. As an application we define image subquandle polynomial enhancements of the G-family counting invariant for trivalent spatial graphs and handlebody-links. We provide examples to show that the new enhancements are proper.


title: Malware Classification with Word Embedding Features

abstract: Malware classification is an important and challenging problem in information security. Modern malware classification techniques rely on machine learning models that can be trained on features such as opcode sequences, API calls, and byte $n$-grams, among many others. In this research, we consider opcode features. We implement hybrid machine learning techniques, where we engineer feature vectors by training hidden Markov models -- a technique that we refer to as HMM2Vec -- and Word2Vec embeddings on these opcode sequences. The resulting HMM2Vec and Word2Vec embedding vectors are then used as features for classification algorithms. Specifically, we consider support vector machine (SVM), $k$-nearest neighbor ($k$-NN), random forest (RF), and convolutional neural network (CNN) classifiers. We conduct substantial experiments over a variety of malware families. Our experiments extend well beyond any previous work in this field.


title: A classification of ideals in Steinberg and Leavitt path algebras over arbitrary rings

abstract: We give a one-to-one correspondence between ideals in the Steinberg algebra of a Hausdorff ample groupoid $G$, and certain families of ideals in the group algebras of isotropy groups in $G$. This generalises a known ideal correspondence theorem for Steinberg algebras of strongly effective groupoids. We use this to give a complete graph-theoretic description of the ideal lattice of Leavitt path algebras over arbitrary commutative rings, generalising the classification of ideals in Leavitt path algebras over fields.


title: Field Theories for type-II fractons

abstract: We derive an effective field theory for a type-II fracton starting from the Haah code on the lattice. The effective topological theory is not given exclusively in terms of an action; it must be supplemented with a condition that selects physical states. Without the constraint, the action only describes a type-I fracton. The constraint emerges from a condition that cube operators multiply to the identity, and it cannot be consistently implemented in the continuum theory at the operator level, but only in a weaker form, in terms of matrix elements of physical states. Informed by these studies and starting from the opposite end, i.e., the continuum, we discuss a Chern-Simons-like theory that does not need a constraint or projector, and yet has no mobile excitations. Whether this continuum theory admits a lattice counterpart remains unanswered.


title: The effect of chaos on the simulation of quantum critical phenomena in analog quantum simulators

abstract: We study how chaos, introduced by a weak perturbation, affects the reliability of the output of analog quantum simulation. As a toy model, we consider the Lipkin-Meshkov-Glick (LMG) model. Inspired by the semiclassical behavior of the order parameter in the thermodynamic limit, we propose a protocol to measure the quantum phase transition in the ground state and the dynamical quantum phase transition associated with quench dynamics. We show that the presence of a small time-dependent perturbation can render the dynamics of the system chaotic. We then show that the estimates of the critical points of these quantum phase transitions, obtained from the quantum simulation of its dynamics, are robust to the presence of this chaotic perturbation, while other aspects of the system, such as the mean magnetization are fragile, and therefore cannot be reliably extracted from this simulator. This can be understood in terms of the simulated quantities that depend on the global structure of phase space vs. those that depend on local trajectories.


title: Dynamical System Analysis of a Dirac-Born-Infeld Model : A Center Manifold Perspective

abstract: In this paper we present the cosmological dynamics of a perfect fluid and the Dark Energy (DE) component of the Universe, where our model of the dark energy is the string-theoritic Dirac-Born-Infeld (DBI) model. We assume that the potential of the scalar field and the warp factor of the warped throat region of the compact space in the extra dimension for the DBI model are both exponential in nature. In the background of spatially flat Friedman-Robertson-Walker-Lemaitre Universe, the Einstein field equations for the DBI dark energy reduce to a system of autonomous dynamical system. We then perform a dynamical system analysis for this system. Our analysis is motivated by the invariant manifold approach of the mathematical dynamics. In this method, it is possible to reach a definite conclusion even when the critical points of a dynamical system are non-hyperbolic in nature. Since we find the complete set of critical points for this system, the center manifold analysis ensures that our investigation of this model leaves no stone unturned. We find some interesting results such as that for some critical points there are situations where scaling solutions exist. Finally we present various topologically different phase planes and stability diagrams and discuss the corresponding cosmological scenario.


title: Policy Decomposition: Approximate Optimal Control with Suboptimality Estimates

abstract: Numerically computing global policies to optimal control problems for complex dynamical systems is mostly intractable. In consequence, a number of approximation methods have been developed. However, none of the current methods can quantify by how much the resulting control underperforms the elusive globally optimal solution. Here we propose policy decomposition, an approximation method with explicit suboptimality estimates. Our method decomposes the optimal control problem into lower-dimensional subproblems, whose optimal solutions are recombined to build a control policy for the entire system. Many such combinations exist, and we introduce the value error and its LQR and DDP estimates to predict the suboptimality of possible combinations and prioritize the ones that minimize it. Using a cart-pole, a 3-link balancing biped and N-link planar manipulators as example systems, we find that the estimates correctly identify the best combinations, yielding control policies in a fraction of the time it takes to compute the optimal control without a notable sacrifice in closed-loop performance. While more research will be needed to find ways of dealing with the combinatorics of policy decomposition, the results suggest this method could be an effective alternative for approximating optimal control in intractable systems.


title: Uncertainty in Self-Adaptive Systems: A Research Community Perspective

abstract: One of the primary drivers for self-adaptation is ensuring that systems achieve their goals regardless of the uncertainties they face during operation. Nevertheless, the concept of uncertainty in self-adaptive systems is still insufficiently understood. Several taxonomies of uncertainty have been proposed, and a substantial body of work exists on methods to tame uncertainty. Yet, these taxonomies and methods do not fully convey the research community's perception on what constitutes uncertainty in self-adaptive systems and how to tackle it. To understand this perception and learn from it, we conducted a survey comprising two complementary stages. In the first stage, we focused on current research and development. In the second stage, we focused on directions for future research. The key findings of the first stage are: a) an overview of uncertainty sources considered in self-adaptive systems, b) an overview of existing methods used to tackle uncertainty in concrete applications, c) insights into the impact of uncertainty on non-functional requirements, d) insights into different opinions in the perception of uncertainty within the community, and the need for standardised uncertainty-handling processes to facilitate uncertainty management in self-adaptive systems. The key findings of the second stage are: a) the insight that over 70% of the participants believe that self-adaptive systems can be engineered to cope with unanticipated change, b) a set of potential approaches for dealing with unanticipated change, c) a set of open challenges in mitigating uncertainty in self-adaptive systems, in particular in those with safety-critical requirements. From these findings, we outline an initial reference process to manage uncertainty in self-adaptive systems.


title: A Modified Drake Equation for Assessing Adversarial Risk to Machine Learning Models

abstract: Machine learning models present a risk of adversarial attack when deployed in production. Quantifying the contributing factors and uncertainties using empirical measures could assist the industry with assessing the risk of downloading and deploying common model types. This work proposes modifying the traditional Drake Equation's formalism to estimate the number of potentially successful adversarial attacks on a deployed model. The Drake Equation is famously used for parameterizing uncertainties and it has been used in many research fields outside of its original intentions to estimate the number of radio-capable extra-terrestrial civilizations. While previous work has outlined methods for discovering vulnerabilities in public model architectures, the proposed equation seeks to provide a semi-quantitative benchmark for evaluating and estimating the potential risk factors for adversarial attacks.


title: Delay-dependent and delay-independent stability of Cournot duopoly model with tax evasion and time-delay

abstract: In this paper a stability analysis for a Cournot duopoly model with tax evasion and time-delay in a continuous-time framework is presented. The mathematical model under consideration follows a gradient dynamics approach, is nonlinear and four-dimensional with state variables given by the production and declared revenue of each competitor. We prove that both the marginal cost rate and time delay play roles as bifurcation parameters. More precisely, if the marginal cost rate lies in certain closed interval then the equilibrium point is delay-independent stable, otherwise it is delay-dependent stable and a Hopf bifurcation necessarily occurs. Some numerical simulations are presented in order to confirm the proposed theoretical results and illustrate the effect of the bifurcation parameters on model stability.


title: Pressureless stationary solutions in a Newton-Yukawa gravity model

abstract: Non-minimally coupled curvature-matter gravity models are an interesting alternative to the Theory of General Relativity and to address the dark energy and dark matter cosmological problems. These models have complex field equations that prevent a full analytical study. Nonetheless, in a particular limit, the behavior of a matter distribution can, in these models, be described by a Schr\"odinger-Newton system. In nonlinear optics, the Schr\"odinger-Newton system can be used to tackle a wide variety of relevant situations and several numerical tools have been developed for this purpose. Interestingly, these methods can be adapted to study General Relativity problems as well as its extensions. In this work, we report the use of these numerical tools to study a particular non-minimal coupling model that introduces two new potentials, an attractive Yukawa potential and a repulsive potential proportional to the energy density. Using the imaginary-time propagation method we have shown that stationary solutions arise even at low energy density regimes.


title: Importance Sampling with the Integrated Nested Laplace Approximation

abstract: The Integrated Nested Laplace Approximation (INLA) is a deterministic approach to Bayesian inference on latent Gaussian models (LGMs) and focuses on fast and accurate approximation of posterior marginals for the parameters in the models. Recently, methods have been developed to extend this class of models to those that can be expressed as conditional LGMs by fixing some of the parameters in the models to descriptive values. These methods differ in the manner descriptive values are chosen. This paper proposes to combine importance sampling with INLA (IS-INLA), and extends this approach with the more robust adaptive multiple importance sampling algorithm combined with INLA (AMIS-INLA). This paper gives a comparison between these approaches and existing methods on a series of applications with simulated and observed datasets and evaluates their performance based on accuracy, efficiency, and robustness. The approaches are validated by exact posteriors in a simple bivariate linear model; then, they are applied to a Bayesian lasso model, a Bayesian imputation of missing covariate values, and lastly, in parametric Bayesian quantile regression. The applications show that the AMIS-INLA approach, in general, outperforms the other methods, but the IS-INLA algorithm could be considered for faster inference when good proposals are available.


title: A Markov jump process associated with the matrix-exponential distribution

abstract: Let $f$ be the density function associated to a matrix-exponential distribution of parameters $(\alpha, T,s)$. By exponentially tilting $f$, we find a probabilistic interpretation which generalises the one associated to phase-type distributions. More specifically, we show that for any sufficiently large $\lambda\ge 0$, the function $x\mapsto \left(\int_0^\infty e^{-\lambda r}f(r) dr\right)^{-1}e^{-\lambda x}f(x)$ can be described in terms of a Markov jump process whose generator is tied to $T$. Finally, we show how to revert the exponential tilting in order to assign a probabilistic interpretation to $f$ itself.


title: Identification and Development of Therapeutics for COVID-19

abstract: After emerging in China in late 2019, the novel Severe acute respiratory syndrome-like coronavirus 2 (SARS-CoV-2) spread worldwide and as of early 2021, continues to significantly impact most countries. Only a small number of coronaviruses are known to infect humans, and only two are associated with the severe outcomes associated with SARS-CoV-2: Severe acute respiratory syndrome-related coronavirus, a closely related species of SARS-CoV-2 that emerged in 2002, and Middle East respiratory syndrome-related coronavirus, which emerged in 2012. Both of these previous epidemics were controlled fairly rapidly through public health measures, and no vaccines or robust therapeutic interventions were identified. However, previous insights into the immune response to coronaviruses gained during the outbreaks of severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) have proved beneficial to identifying approaches to the treatment and prophylaxis of novel coronavirus disease 2019 (COVID-19). A number of potential therapeutics against SARS-CoV-2 and the resultant COVID-19 illness were rapidly identified, leading to a large number of clinical trials investigating a variety of possible therapeutic approaches being initiated early on in the pandemic. As a result, a small number of therapeutics have already been authorized by regulatory agencies such as the Food and Drug Administration (FDA) in the United States, and many other therapeutics remain under investigation. Here, we describe a range of approaches for the treatment of COVID-19, along with their proposed mechanisms of action and the current status of clinical investigation into each candidate. The status of these investigations will continue to evolve, and this review will be updated as progress is made.


title: Skyrmion Logic Clocked via Voltage Controlled Magnetic Anisotropy

abstract: Magnetic skyrmions are exciting candidates for energy-efficient computing due to their non-volatility, detectability,and mobility. A recent proposal within the paradigm of reversible computing enables large-scale circuits composed ofdirectly-cascaded skyrmion logic gates, but it is limited by the manufacturing difficulty and energy costs associated withthe use of notches for skyrmion synchronization. To overcome these challenges, we therefore propose a skyrmion logicsynchronized via modulation of voltage-controlled magnetic anisotropy (VCMA). In addition to demonstrating theprinciple of VCMA synchronization through micromagnetic simulations, we also quantify the impacts of current den-sity, skyrmion velocity, and anisotropy barrier height on skyrmion motion. Further micromagnetic results demonstratethe feasibility of cascaded logic circuits in which VCMA synchronizers enable clocking and pipelining, illustrating afeasible pathway toward energy-efficient large-scale computing systems based on magnetic skyrmions.


title: Laser soliton microcombs on silicon

abstract: Silicon photonics enables wafer-scale integration of optical functionalities on chip. A silicon-based laser frequency combs could significantly expand the applications of silicon photonics, by providing integrated sources of mutually coherent laser lines for terabit-per-second transceivers, parallel coherent LiDAR, or photonics-assisted signal processing. Here, we report on heterogeneously integrated laser soliton microcombs combining both InP/Si semiconductor lasers and ultralow-loss silicon nitride microresonators on monolithic silicon substrate. Thousands of devices are produced from a single wafer using standard CMOS techniques. Using on-chip electrical control of the microcomb-laser relative optical phase, these devices can output single-soliton microcombs with 100 GHz repetition rate. Our approach paves the way for large-volume, low-cost manufacturing of chip-based frequency combs for next-generation high-capacity transceivers, datacenters, space and mobile platforms.


title: Implicit Methods with Reduced Memory for Thermal Radiative Transfer

abstract: This paper presents approximation methods for time-dependent thermal radiative transfer problems in high energy density physics. It is based on the multilevel quasidiffusion method defined by the high-order radiative transfer equation (RTE) and the low-order quasidiffusion (aka VEF) equations for the moments of the specific intensity. A large part of data storage in TRT problems between time steps is determined by the dimensionality of grid functions of the radiation intensity. The approximate implicit methods with reduced memory for the time-dependent Boltzmann equation are applied to the high-order RTE, discretized in time with the backward Euler (BE) scheme. The high-dimensional intensity from the previous time level in the BE scheme is approximated by means of the low-rank proper orthogonal decomposition (POD). Another version of the presented method applies the POD to the remainder term of P2 expansion of the intensity. The accuracy of the solution of the approximate implicit methods depends of the rank of the POD. The proposed methods enable one to reduce storage requirements in time dependent problems. Numerical results of a Fleck-Cummings TRT test problem are presented.


title: Preference-based Learning of Reward Function Features

abstract: Preference-based learning of reward functions, where the reward function is learned using comparison data, has been well studied for complex robotic tasks such as autonomous driving. Existing algorithms have focused on learning reward functions that are linear in a set of trajectory features. The features are typically hand-coded, and preference-based learning is used to determine a particular user's relative weighting for each feature. Designing a representative set of features to encode reward is challenging and can result in inaccurate models that fail to model the users' preferences or perform the task properly. In this paper, we present a method to learn both the relative weighting among features as well as additional features that help encode a user's reward function. The additional features are modeled as a neural network that is trained on the data from pairwise comparison queries. We apply our methods to a driving scenario used in previous work and compare the predictive power of our method to that of only hand-coded features. We perform additional analysis to interpret the learned features and examine the optimal trajectories. Our results show that adding an additional learned feature to the reward model enhances both its predictive power and expressiveness, producing unique results for each user.