title: Maximizing the Smallest Eigenvalue of Grounded Laplacian Matrix

abstract: For a connected graph $\mathcal{G}=(V,E)$ with $n$ nodes, $m$ edges, and Laplacian matrix $\boldsymbol{{\mathit{L}}}$, a grounded Laplacian matrix $\boldsymbol{{\mathit{L}}}(S)$ of $\mathcal{G}$ is a $(n-k) \times (n-k)$ principal submatrix of $\boldsymbol{{\mathit{L}}}$, obtained from $\boldsymbol{{\mathit{L}}}$ by deleting $k$ rows and columns corresponding to $k$ selected nodes forming a set $S\subseteq V$. The smallest eigenvalue $\lambda(S)$ of $\boldsymbol{{\mathit{L}}}(S)$ plays a pivotal role in various dynamics defined on $\mathcal{G}$. For example, $\lambda(S)$ characterizes the convergence rate of leader-follower consensus, as well as the effectiveness of a pinning scheme for the pinning control problem, with larger $\lambda(S)$ corresponding to smaller convergence time or better effectiveness of a pinning scheme. In this paper, we focus on the problem of optimally selecting a subset $S$ of fixed $k \ll n$ nodes, in order to maximize the smallest eigenvalue $\lambda(S)$ of the grounded Laplacian matrix $\boldsymbol{{\mathit{L}}}(S)$. We show that this optimization problem is NP-hard and that the objective function is non-submodular but monotone. Due to the difficulty to obtain the optimal solution, we first propose a na\"{\i}ve heuristic algorithm selecting one optimal node at each time for $k$ iterations. Then we propose a fast heuristic scalable algorithm to approximately solve this problem, using derivative matrix, matrix perturbations, and Laplacian solvers as tools. Our na\"{\i}ve heuristic algorithm takes $\tilde{O}(knm)$ time, while the fast greedy heuristic has a nearly linear time complexity of $\tilde{O}(km)$. We also conduct numerous experiments on different networks sized up to one million nodes, demonstrating the superiority of our algorithm in terms of efficiency and effectiveness.


title: Simulation and Model Checking for Close to Realtime Overtaking Planning

abstract: Fast and reliable trajectory planning is a key requirement of autonomous vehicles. In this paper we introduce a novel technique for planning the route of an autonomous vehicle on a straight rural road using the Spin model checker. We show how we can combine Spins ability to identify paths violating temporal properties with sensor information from a 3D Unity simulation of an autonomous vehicle, to plan and perform consecutive overtaking manoeuvres on a traffic heavy road. This involves discretising the sensory information and combining multiple sequential Spin models with a Linear Time Temporal Logic specification to generate an error path. This path provides the autonomous vehicle with an action plan. The entire process takes place in close to realtime using no precomputed data and the action plan is specifically tailored for individual scenarios. Our experiments demonstrate that the simulated autonomous vehicle implementing our approach can drive on average at least 40km and overtake 214 vehicles before experiencing a collision, which is usually caused by inaccuracies in the sensory system. While the proposed system has some drawbacks, we believe that our novel approach demonstrates a potentially powerful future tool for efficient trajectory planning for autonomous vehicles.


title: Improving Online Railway Deadlock Detection using a Partial Order Reduction

abstract: Although railway dispatching on large national networks is gradually becoming more computerized, there are still major obstacles to retrofitting (semi-)autonomous control systems. In addition to requiring extensive and detailed digitalization of infrastructure models and information systems, exact optimization for railway dispatching is computationally hard. Heuristic algorithms and manual overrides are likely to be required for semi-autonomous railway operations for the foreseeable future. In this context, being able to detect problems such as deadlocks can be a valuable part of a runtime verification system. If bound-for-deadlock situations are correctly recognized as early as possible, human operators will have more time to better plan for recovery operations. Deadlock detection may also be useful for verification in a feedback loop with a heuristic or semi-autonomous dispatching algorithm if the dispatching algorithm cannot itself guarantee a deadlock-free plan. We describe a SAT-based planning algorithm for online detection of bound-for-deadlock situations. The algorithm exploits parallel updates of train positions and a partial order reduction technique to significantly reduce the number of state transitions (and correspondingly, the sizes of the formulas) in the SAT instances needed to prove whether a deadlock situation is bound to happen in the future. Implementation source code and benchmark instances are supplied, and a direct comparison against another recent study demonstrates significant performance gains.


title: Observable and Attention-Directing BDI Agents for Human-Autonomy Teaming

abstract: Human-autonomy teaming (HAT) scenarios feature humans and autonomous agents collaborating to meet a shared goal. For effective collaboration, the agents must be transparent and able to share important information about their operation with human teammates. We address the challenge of transparency for Belief-Desire-Intention agents defined in the Conceptual Agent Notation (CAN) language. We extend the semantics to model agents that are observable (i.e. the internal state of tasks is available), and attention-directing (i.e. specific states can be flagged to users), and provide an executable semantics via an encoding in Milner's bigraphs. Using an example of unmanned aerial vehicles, the BigraphER tool, and PRISM, we show and verify how the extensions work in practice.


title: Identifying optimal photovoltaic technologies for underwater applications

abstract: Improving solar energy collection in aquatic environments would allow for superior environmental monitoring and remote sensing, but the identification of optimal photovoltaic technologies for such applications is challenging as evaluation requires either field deployment or access to large water tanks. Here, we present a simple bench-top characterization technique that does not require direct access to water and therefore circumvents the need for field testing during initial trials of development. Employing LEDs to simulate underwater solar spectra at various depths, we compare Si and CdTe solar cells, two commercially available technologies, with GaInP cells, a technology with a wide band gap close to ideal for underwater solar harvesting. We use this method to show that while Si cells outperform both CdTe and GaInP under terrestrial AM1.5G solar irradiance, both CdTe and GaInP outperform Si at depths > 2 m, with GaInP cells operating with underwater efficiencies exceeding 51%.


title: Towards a Formalisation of Justification and Justifiability

abstract: We introduce the logic QKSD which is a normal multi-modal logic over finitely many modalities that additionally supports bounded quantification of modalities. An important feature of this logic is that it allows to quantify over the information components of systems and, hence, can be used to derive justifications. We compare the proposed logic with Artemov's justification logic and also report on a prototypical implementation of a satisfiability solver of this logic and show some examples.


title: Weighted Mean Difference Statistics for Paired Data in Presence of Missing Values

abstract: Missing data is a common issue in many biomedical studies. Under a paired design, some subjects may have missing values in either one or both of the conditions due to loss of follow-up, insufficient biological samples, etc. Such partially paired data complicate statistical comparison of the distribution of the variable of interest between the two conditions. In this paper, we propose a general class of test statistics based on the difference in weighted sample means without imposing any distributional or model assumption. An optimal weight is derived for this class of tests. Simulation studies show that our proposed test with the optimal weight performs well and outperforms existing methods in practical situations. Two cancer biomarker studies are provided for illustration.


title: Extending Urban Multi-Lane Spatial Logic to Formalise Road Junction Rules

abstract: During the design of autonomous vehicles (AVs), several stages should include a verification process to guarantee that the AV is driving safely on the roads. One of these stages is to assure the AVs abide by the road traffic rules. To include road traffic rules in the design of an AV, a precise and unambiguous formalisation of these rules is needed. However, only recently this has been pointed out as an issue for the design of AVs and the few works on this only capture the temporal aspects of the rules, leaving behind the spatial aspects. Here, we extend the spatial traffic logic, Urban Multi-lane Spatial Logic, to formalise a subset of the UK road junction rules, where both temporal and spatial aspects of the rules are captured. Our approach has an abstraction level for urban road junctions that could easily promote the formalisation of the whole set of road junction rules and we exemplarily formalise three of the UK road junction rules. Once we have the whole set formalised, we will model, implement, and formally verify the behaviour of an AV against road traffic rules so that guidelines for the creation of a Digital Highway Code for AVs can be established.


title: Chords of 2-factors in planar cubic bridgeless graphs

abstract: We show that every edge in a 2-edge-connected planar cubic graph is either contained in a 2-edge-cut or is a chord of some cycle that is contained in a 2-factor of the graph. As a consequence, we show that every edge in a cyclically 4-edge-connected planar cubic graph, except $K_2^3$ and $K_4$, is contained in a perfect matching whose removal disconnects the graph. We obtain a complete characterization of 2-edge-connected planar cubic graphs that have an edge such that every 2-factor containing the edge is a Hamiltonian cycle, and also of those that have an edge such that the complement of every perfect matching containing the edge is a Hamiltonian cycle. Another immediate consequence of the main result is that for any two edges contained in a facial cycle of a 2-edge-connected planar cubic graph, there exists a 2-factor in the graph such that both edges are contained in the same cycle of the 2-factor. We conjecture that this property holds for any two edges in a 2-edge-connected planar cubic graph, and prove it for planar cubic bipartite graphs. The main result is proved in the dual form by showing that every plane triangulation admits a vertex 3-coloring such that no face is monochromatic and there is exactly one specified edge between a specified pair of color classes.


title: Towards Partial Monitoring: It is Always too Soon to Give Up

abstract: Runtime Verification is a lightweight formal verification technique. It is used to verify at runtime whether the system under analysis behaves as expected. The expected behaviour is usually formally specified by means of properties, which are used to automatically synthesise monitors. A monitor is a device that, given a sequence of events representing a system execution, returns a verdict symbolising the satisfaction or violation of the formal property. Properties that can (resp. cannot) be verified at runtime by a monitor are called monitorable and non-monitorable, respectively. In this paper, we revise the notion of monitorability from a practical perspective, where we show how non-monitorable properties can still be used to generate partial monitors, which can partially check the properties. Finally, we present the implications both from a theoretical and practical perspectives.


title: Complete Agent-driven Model-based System Testing for Autonomous Systems

abstract: In this position paper, a novel approach to testing complex autonomous transportation systems (ATS) in the automotive, avionic, and railway domains is described. It is intended to mitigate some of the most critical problems regarding verification and validation (V&V) effort for ATS. V&V is known to become infeasible for complex ATS, when using conventional methods only. The approach advocated here uses complete testing methods on the module level, because these establish formal proofs for the logical correctness of the software. Having established logical correctness, system-level tests are performed in simulated cloud environments and on the target system. To give evidence that 'sufficiently many' system tests have been performed with the target system, a formally justified coverage criterion is introduced. To optimise the execution of very large system test suites, we advocate an online testing approach where multiple tests are executed in parallel, and test steps are identified on-the-fly. The coordination and optimisation of these executions is achieved by an agent-based approach. Each aspect of the testing approach advocated here is shown to either be consistent with existing standards for development and V&V of safety-critical transportation systems, or it is justified why it should become acceptable in future revisions of the applicable standards.


title: Formal Guarantees of Timely Progress for Distributed Knowledge Propagation

abstract: Autonomous air traffic management (ATM) operations for urban air mobility (UAM) will necessitate the use of distributed protocols for decentralized coordination between aircraft. As UAM operations are time-critical, it will be imperative to have formal guarantees of progress for the distributed protocols used in ATM. Under asynchronous settings, message transmission and processing delays are unbounded, making it impossible to provide deterministic bounds on the time required to make progress. We present an approach for formally guaranteeing timely progress in a Two-Phase Acknowledge distributed knowledge propagation protocol by probabilistically modeling the delays using theories of the Multicopy Two-Hop Relay protocol and the M/M/1 queue system. The guarantee states a probabilistic upper bound to the time for progress as a function of the probabilities of the total transmission and processing delays being less than two given values. We also showcase the development of a library of formal theories, that is tailored towards reasoning about timely progress in distributed protocols deployed in airborne networks, in the Athena proof assistant.


title: QuantifyML: How Good is my Machine Learning Model?

abstract: The efficacy of machine learning models is typically determined by computing their accuracy on test data sets. However, this may often be misleading, since the test data may not be representative of the problem that is being studied. With QuantifyML we aim to precisely quantify the extent to which machine learning models have learned and generalized from the given data. Given a trained model, QuantifyML translates it into a C program and feeds it to the CBMC model checker to produce a formula in Conjunctive Normal Form (CNF). The formula is analyzed with off-the-shelf model counters to obtain precise counts with respect to different model behavior. QuantifyML enables i) evaluating learnability by comparing the counts for the outputs to ground truth, expressed as logical predicates, ii) comparing the performance of models built with different machine learning algorithms (decision-trees vs. neural networks), and iii) quantifying the safety and robustness of models.


title: Complete Test of Synthesised Safety Supervisors for Robots and Autonomous Systems

abstract: Verified controller synthesis uses world models that comprise all potential behaviours of humans, robots, further equipment, and the controller to be synthesised. A world model enables quantitative risk assessment, for example, by stochastic model checking. Such a model describes a range of controller behaviours some of which -- when implemented correctly -- guarantee that the overall risk in the actual world is acceptable, provided that the stochastic assumptions have been made to the safe side. Synthesis then selects an acceptable-risk controller behaviour. However, because of crossing abstraction, formalism, and tool boundaries, verified synthesis for robots and autonomous systems has to be accompanied by rigorous testing. In general, standards and regulations for safety-critical systems require testing as a key element to obtain certification credit before entry into service. This work-in-progress paper presents an approach to the complete testing of synthesised supervisory controllers that enforce safety properties in domains such as human-robot collaboration and autonomous driving. Controller code is generated from the selected controller behaviour. The code generator, however, is hard, if not infeasible, to verify in a formal and comprehensive way. Instead, utilising testing, an abstract test reference is generated, a symbolic finite state machine with simpler semantics than code semantics. From this reference, a complete test suite is derived and applied to demonstrate the observational equivalence between the synthesised abstract test reference and the generated concrete controller code running on a control system platform.


title: Online Strategy Synthesis for Safe and Optimized Control of Steerable Needles

abstract: Autonomous systems are often applied in uncertain environments, which require prospective action planning and retrospective data evaluation for future planning to ensure safe operation. Formal approaches may support these systems with safety guarantees, but are usually expensive and do not scale well with growing system complexity. In this paper, we introduce online strategy synthesis based on classical strategy synthesis to derive formal safety guarantees while reacting and adapting to environment changes. To guarantee safety online, we split the environment into region types which determine the acceptance of action plans and trigger local correcting actions. Using model checking on a frequently updated model, we can then derive locally safe action plans (prospectively), and match the current model against new observations via reachability checks (retrospectively). As use case, we successfully apply online strategy synthesis to medical needle steering, i.e., navigating a (flexible and beveled) needle through tissue towards a target without damaging its surroundings.


title: Assuring Increasingly Autonomous Systems in Human-Machine Teams: An Urban Air Mobility Case Study

abstract: As aircraft systems become increasingly autonomous, the human-machine role allocation changes and opportunities for new failure modes arise. This necessitates an approach to identify the safety requirements for the increasingly autonomous system (IAS) as well as a framework and techniques to verify and validate that an IAS meets its safety requirements. We use Crew Resource Management techniques to identify requirements and behaviors for safe human-machine teaming behaviors. We provide a methodology to verify that an IAS meets its requirements. We apply the methodology to a case study in Urban Air Mobility, which includes two contingency scenarios: unreliable sensor and aborted landing. For this case study, we implement an IAS agent in the Soar language that acts as a copilot for the selected contingency scenarios and performs takeoff and landing preparation, while the pilot maintains final decision authority. We develop a formal human-machine team architecture model in the Architectural Analysis and Design Language (AADL), with operator and IAS requirements formalized in the Assume Guarantee REasoning Environment (AGREE) Annex to AADL. We formally verify safety requirements for the human-machine team given the requirements on the IAS and operator. We develop an automated translator from Soar to the nuXmv model checking language and formally verify that the IAS agent satisfies its requirements using nuXmv. We share the design and requirements errors found in the process as well as our lessons learned.


title: Collective plasma effects of electron-positron pairs in beam-driven QED cascades

abstract: Understanding the interplay of strong-field QED and collective plasma effects is important for explaining extreme astrophysical environments like magnetars. It has been shown that QED pair plasmas is possible to be produced and observed by passing a relativistic electron beam through an intense laser field. This paper presents in detail multiple sets of 3D QED-PIC simulations to show the creation of pair plasmas in the QED cascade. The beam driven method enables a high pair particle density and also a low particle gamma factor, which both play equal rolls on exhibiting large collective plasma effects. Finite laser frequency upshift is observed with both ideal parameters (24 PW laser laser colliding with 300 GeV electron beam) and with existing technologies (3 PW laser laser colliding with 30 GeV electron beam).


title: First-principles study of defect formation energies in LaO$X$S$_2$ ($X=$ Sb, Bi)

abstract: We theoretically investigate defect formation energies in LaO$X$S$_2$ ($X=$Sb, Bi) using first-principles calculation. We find that the oxygen vacancy is relatively stable, where its formation energy is higher in $X=$ Sb than in $X=$ Bi. An interesting feature of $X=$ Sb is that the vacancy of the in-plane sulfur atom becomes more stable than in $X=$ Bi, caused by the formation of an Sb$_2$ dimer and the electron occupation of the impurity energy levels. The formation energies of cation defects and anion-cation antisite defects are positive for the chemical equilibrium condition used in this study. Fluorine likely replaces oxygen, and its defect formation energy is negative for both $X=$ Sb and Bi, while that for $X=$ Sb is much higher than $X=$ Bi. Our study clarifies the stability of several point defects and suggests that the in-plane structural instability is enhanced in $X=$ Sb, which seems to affect a structural change caused by some in-plane point defects.


title: Developing a Meta-suggestion Engine for Search Queries

abstract: Typically, search engines provide query suggestions to assist users in the search process. Query suggestions are very important for improving users search experience. However, most query suggestions are based on the user's search logs, and they can be influenced by infrequently searched queries. Depending on the user's query, query suggestions can be ineffective in global search engines but effective in a domestic search engine. Conversely, it can be effective in global engines and weak in domestic engines. In addition, log-based query suggestions require many search logs, which makes them difficult to construct outside of a large search engine. Some search engines do not provide query suggestions, making searches difficult for users. These query suggestion vulnerabilities degrade the user's search experience. In this study, we develop a meta-suggestion, a new query suggestion scheme. Similar to meta-searches, meta-suggestions retrieves candidate queries of suggestions from other search engines. Meta-suggestions generate suggestions by reranking the aggregated candidate queries. We develop a meta-suggestion engine (MSE) browser extension that generates meta-suggestions. It can provide query suggestions for any webpage and does not require a search log. Comparing our meta-suggestions to major search engines such as Google, showed a 17% performance improvement on normalized discounted cumulative gain (NDCG) and a 31% improvement on precision. If more detailed factors, such as user preferences are discovered through continued research, it is expected that user searches will greatly improve. An enhanced user search experience is possible if factors, such as user preference, are examined in future work.


title: Fast Rank-1 NMF for Missing Data with KL Divergence

abstract: We propose a fast non-gradient-based method of rank-1 non-negative matrix factorization (NMF) for missing data, called A1GM, that minimizes the KL divergence from an input matrix to the reconstructed rank-1 matrix. Our method is based on our new finding of an analytical closed-formula of the best rank-1 non-negative multiple matrix factorization (NMMF), a variety of NMF. NMMF is known to exactly solve NMF for missing data if positions of missing values satisfy a certain condition, and A1GM transforms a given matrix so that the analytical solution to NMMF can be applied. We empirically show that A1GM is more efficient than a gradient method with competitive reconstruction errors.


title: GeoSneakPique: Visual Autocompletion for Geospatial Queries

abstract: How many crimes occurred in the city center? And exactly which part of town is the 'city center'? While location is at the heart of many data questions, geographic location can be difficult to specify in natural language (NL) queries. This is especially true when working with fuzzy cognitive regions or regions that may be defined based on data distributions instead of absolute administrative location (e.g., state, country). GeoSneakPique presents a novel method for using a mapping widget to support the NL query process, allowing users to specify location via direct manipulation with data-driven guidance on spatial distributions to help select the area of interest. Users receive feedback to help them evaluate and refine their spatial selection interactively and can save spatial definitions for re-use in subsequent queries. We conduct a qualitative evaluation of the GeoSneakPique that indicates the usefulness of the interface as well as opportunities for better supporting geospatial workflows in visual analysis tasks employing cognitive regions.


title: On Gromov-Yomdin type theorems and a categorical interpretation of holomorphicity

abstract: In topological dynamics, the Gromov--Yomdin theorem states that the topological entropy of a holomorphic automorphism $f$ of a smooth projective variety is equal to the logarithm of the spectral radius of the induced map $f^*$. In order to establish a categorical analogue of the Gromov--Yomdin theorem, one first needs to find a categorical analogue of a holomorphic automorphism. In this paper, we propose a categorical analogue of a holomorphic automorphism and prove that the Gromov--Yomdin type theorem holds for them.


title: Development of an Undergraduate Quantum Engineering Degree

abstract: Quantum technology is exploding. Computing, communication, and sensing are just a few areas likely to see breakthroughs in the next few years. Worldwide, national governments, industries, and universities are moving to create a new class of workforce - the Quantum Engineers. Demand for such engineers is predicted to be in the tens of thousands within a five-year timescale. However, how best to train this next generation of engineers is far from obvious. Quantum mechanics - long a pillar of traditional physics undergraduate degrees - must now be merged with traditional engineering offerings. This paper discusses the history, development, and first year of operation of the world's first undergraduate degree in quantum engineering. The main purpose of the paper is to inform the wider debate, now being held by many institutions worldwide, on how best to formally educate the Quantum Engineer.


title: Sparse varying-coefficient functional linear model

abstract: We consider the problem of variable selection in varying-coefficient functional linear models, where multiple predictors are functions and a response is a scalar and depends on an exogenous variable. The varying-coefficient functional linear model is estimated by the penalized maximum likelihood method with the sparsity-inducing penalty. Tuning parameters that controls the degree of the penalization are determined by a model selection criterion. The proposed method can reveal which combination of functional predictors relates to the response, and furthermore how each predictor relates to the response by investigating coefficient surfaces. Simulation studies are provided to investigate the effectiveness of the proposed method. We also apply it to the analysis of crop yield data to investigate which combination of environmental factors relates to the amount of a crop yield.


title: Local $\varepsilon$-regularity criteria for the five dimensional stationary Navier-Stokes equations

abstract: In this paper, we establish $\varepsilon$-regularity criteria at one scale for suitable weak solutions to the five dimensional stationary incompressible Navier-Stokes equations in both the unit ball $B_1$ and the unit half ball $B_1^+$, respectively, which improve previous results in [J. Math. Fluid Mech., 6(2004), pp 78--101] and [Comm. Pure Appl. Math., 41(1988), pp 437--458]. The proofs of the main results are based on an iteration argument, Campanato's method, and interpolation techniques.


title: Semantic Resizing of Charts Through Generalization:A Case Study with Line Charts

abstract: Inspired by cartographic generalization principles, we present a generalization technique for rendering line charts at different sizes, preserving the important semantics of the data at that display size. The algorithm automatically determines the generalization operators to be applied at that size based on spatial density, distance, and the semantic importance of the various visualization elements in the line chart. A qualitative evaluation of the prototype that implemented the algorithm indicates that the generalized line charts pre-served the general data shape, while minimizing visual clutter. We identify future opportunities where generalization can be extended and applied to other chart types and visual analysis authoring tools.


title: Dynamic Influence Maximization

abstract: We initiate a systematic study on $\mathit{dynamic}$ $\mathit{influence}$ $\mathit{maximization}$ (DIM). In the DIM problem, one maintains a seed set $S$ of at most $k$ nodes in a dynamically involving social network, with the goal of maximizing the expected influence spread while minimizing the amortized updating cost. We consider two evolution models. In the $\mathit{incremental}$ model, the social network gets enlarged over time and one only introduces new users and establishes new social links, we design an algorithm that achieves $(1-1/e-\epsilon)$-approximation to the optimal solution and has $k \cdot\mathsf{poly}(\log n, \epsilon^{-1})$ amortized running time, which matches the state-of-art offline algorithm with only poly-logarithmic overhead. In the $\mathit{fully}$ $\mathit{dynamic}$ model, users join in and leave, influence propagation gets strengthened or weakened in real time, we prove that under the Strong Exponential Time Hypothesis (SETH), no algorithm can achieve $2^{-(\log n)^{1-o(1)}}$-approximation unless the amortized running time is $n^{1-o(1)}$. On the technical side, we exploit novel adaptive sampling approaches that reduce DIM to the dynamic MAX-k coverage problem, and design an efficient $(1-1/e-\epsilon)$-approximation algorithm for it. Our lower bound leverages the recent developed distributed PCP framework.


title: Common Information based Approximate State Representations in Multi-Agent Reinforcement Learning

abstract: Due to information asymmetry, finding optimal policies for Decentralized Partially Observable Markov Decision Processes (Dec-POMDPs) is hard with the complexity growing doubly exponentially in the horizon length. The challenge increases greatly in the multi-agent reinforcement learning (MARL) setting where the transition probabilities, observation kernel, and reward function are unknown. Here, we develop a general compression framework with approximate common and private state representations, based on which decentralized policies can be constructed. We derive the optimality gap of executing dynamic programming (DP) with the approximate states in terms of the approximation error parameters and the remaining time steps. When the compression is exact (no error), the resulting DP is equivalent to the one in existing work. Our general framework generalizes a number of methods proposed in the literature. The results shed light on designing practically useful deep-MARL network structures under the "centralized learning distributed execution" scheme.


title: Capillary gravity water waves linearized at monotone shear flows: eigenvalues and inviscid damping

abstract: This paper is concerned with the eigenvalues and linear inviscid damping of the 2D capillary gravity water waves of finite depth $x_2\in(-h,0)$ linearized at a monotone shear flow $U(x_2)$. Unlike the linearized Euler equation in a fixed channel where eigenvalues exist only in low horizontal wave number $k$, we first prove the linearized capillary gravity wave has two branches of eigenvalues $-ikc^\pm(k)$, where the wave speeds $c^\pm(k)=O(\sqrt{|k|})$ for $|k|\gg1$ have the same asymptotics as the those of the linear irrotational capillary gravity waves. Under the additional assumption of $U"\ne0$, we obtain the complete continuation of these two branches, which are all the eigenvalues in this (and some other) case(s). Particularly $-ikc^-(k)$ could bifurcate into unstable eigenvalues at $c^-(k)=U(-h)$. The bifurcation of unstable eigenvalues from inflection values of $U$ is also proved. Assuming no singular modes, i.e. no embedded eigenvalues for any wave number $k$, linear solutions $(v(t,x),\eta(t,x_1))$ are studieded in both periodic-in-$x_1$ and $x_1\in R$ cases, where $v$ is the velocity and $\eta$ the surface profile. Solutions can be split into $(v^p,\eta^p)$ and $(v^c,\eta^c)$ whose $k$-th Fourier mode in $x_1$ correspond to the eigenvalues and the continuous spectra of wave number $k$, respectively. The component $(v^p,\eta^p)$ is governed by a (possibly unstable) dispersion relation given by the eigenvalues, which are simply $k\to-ikc^\pm(k)$ in the case of $x_1\in R$. The other component $(v^c,\eta^c)$ satisfies the inviscid damping as fast as $|v_1^c|_{L_x^2},|\eta^c|_{L_x^2}=O(|t|^{-1})$ and $|v_2^c|_{L_x^2}=O(t^{-2})$ as $|t|\gg1$. Additional decay of $tv_1^c,t^2v_2^c$ in $L_x^2L_t^q$, $q\in(2,\infty]$, is obtained after leading asymptotic terms are removed, which are in the forms of $t$-dependent translations in $x_1$ of certain functions of $x$.


title: Investigations on the optical properties and X-ray emission in compact radio AGN

abstract: Compact radio active galactic nuclei (compact radio AGN) are compact ($\leq$ 20 kpc), powerful radio sources. Currently, the preferred scenario is that they are at the early stage of AGN evolution. At present, the research of compact radio AGN mainly focuses on the radio band, other bands have not been extensively studied. We present the systemic optical properties and X-ray emission studies for compact radio AGN, to investigate the accretion properties, AGN evolution and their X-ray origin. We find that compact radio AGN have various accretion modes indicated by the accretion rate analysis. In the radio power-linear size diagram they generally follow the evolutionary trend towards large-scale radio galaxies with increasing linear size and decreasing accretion rate. Their hard X-ray emission may be from jet based on the radio/X-ray relation and fundamental plane of black hole activity.


title: MUSE: Feature Self-Distillation with Mutual Information and Self-Information

abstract: We present a novel information-theoretic approach to introduce dependency among features of a deep convolutional neural network (CNN). The core idea of our proposed method, called MUSE, is to combine MUtual information and SElf-information to jointly improve the expressivity of all features extracted from different layers in a CNN. We present two variants of the realization of MUSE -- Additive Information and Multiplicative Information. Importantly, we argue and empirically demonstrate that MUSE, compared to other feature discrepancy functions, is a more functional proxy to introduce dependency and effectively improve the expressivity of all features in the knowledge distillation framework. MUSE achieves superior performance over a variety of popular architectures and feature discrepancy functions for self-distillation and online distillation, and performs competitively with the state-of-the-art methods for offline distillation. MUSE is also demonstrably versatile that enables it to be easily extended to CNN-based models on tasks other than image classification such as object detection.


title: Least Square Calibration for Peer Review

abstract: Peer review systems such as conference paper review often suffer from the issue of miscalibration. Previous works on peer review calibration usually only use the ordinal information or assume simplistic reviewer scoring functions such as linear functions. In practice, applications like academic conferences often rely on manual methods, such as open discussions, to mitigate miscalibration. It remains an important question to develop algorithms that can handle different types of miscalibrations based on available prior knowledge. In this paper, we propose a flexible framework, namely least square calibration (LSC), for selecting top candidates from peer ratings. Our framework provably performs perfect calibration from noiseless linear scoring functions under mild assumptions, yet also provides competitive calibration results when the scoring function is from broader classes beyond linear functions and with arbitrary noise. On our synthetic dataset, we empirically demonstrate that our algorithm consistently outperforms the baseline which select top papers based on the highest average ratings.


title: Effective-Medium Theory with No Approximation for One Dimensional Metamaterials

abstract: Transfer matrix theory (TMT) is used to study the effective-medium theory (EMT) of one dimensional metamaterials (1D MMs). 1D MMs with equal diagonal elements of periodic transfer matrix are defined as 1D perfect MMs (1D PMMs), which can be suitable for EMT with no approximation.


title: No News is Good News: A Critique of the One Billion Word Benchmark

abstract: The One Billion Word Benchmark is a dataset derived from the WMT 2011 News Crawl, commonly used to measure language modeling ability in natural language processing. We train models solely on Common Crawl web scrapes partitioned by year, and demonstrate that they perform worse on this task over time due to distributional shift. Analysis of this corpus reveals that it contains several examples of harmful text, as well as outdated references to current events. We suggest that the temporal nature of news and its distribution shift over time makes it poorly suited for measuring language modeling ability, and discuss potential impact and considerations for researchers building language models and evaluation datasets.


title: Antenna Array Enabled Space/Air/Ground Communications and Networking for 6G

abstract: Antenna arrays have a long history of more than 100 years and have evolved closely with the development of electronic and information technologies, playing an indispensable role in wireless communications and radar. With the rapid development of electronic and information technologies, the demand for all-time, all-domain, and full-space network services has exploded, and new communication requirements have been put forward on various space/air/ground platforms. To meet the ever increasing requirements of the future sixth generation (6G) wireless communications, such as high capacity, wide coverage, low latency, and strong robustness, it is promising to employ different types of antenna arrays with various beamforming technologies in space/air/ground communication networks, bringing in advantages such as considerable antenna gains, multiplexing gains, and diversity gains. However, enabling antenna array for space/air/ground communication networks poses specific, distinctive and tricky challenges, which has aroused extensive research attention. This paper aims to overview the field of antenna array enabled space/air/ground communications and networking. The technical potentials and challenges of antenna array enabled space/air/ground communications and networking are presented first. Subsequently, the antenna array structures and designs are discussed. We then discuss various emerging technologies facilitated by antenna arrays to meet the new communication requirements of space/air/ground communication systems. Enabled by these emerging technologies, the distinct characteristics, challenges, and solutions for space communications, airborne communications, and ground communications are reviewed. Finally, we present promising directions for future research in antenna array enabled space/air/ground communications and networking.


title: Spontaneously Broken Subsystem Symmetries

abstract: We investigate the spontaneous breaking of subsystem symmetries directly in the context of continuum field theories by calculating the correlation function of charged operators. Our methods confirm the lack of spontaneous symmetry breaking in some of the existing continuum field theories with subsystem symmetries, as had previously been established based on a careful analysis of the spectrum. We present some novel continuum field theory constructions that do exhibit spontaneous symmetry breaking whenever allowed by general principles. These interesting patterns of symmetry breaking occur despite the fact that all the theories we study are non-interacting.


title: DelightfulTTS: The Microsoft Speech Synthesis System for Blizzard Challenge 2021

abstract: This paper describes the Microsoft end-to-end neural text to speech (TTS) system: DelightfulTTS for Blizzard Challenge 2021. The goal of this challenge is to synthesize natural and high-quality speech from text, and we approach this goal in two perspectives: The first is to directly model and generate waveform in 48 kHz sampling rate, which brings higher perception quality than previous systems with 16 kHz or 24 kHz sampling rate; The second is to model the variation information in speech through a systematic design, which improves the prosody and naturalness. Specifically, for 48 kHz modeling, we predict 16 kHz mel-spectrogram in acoustic model, and propose a vocoder called HiFiNet to directly generate 48 kHz waveform from predicted 16 kHz mel-spectrogram, which can better trade off training efficiency, modelling stability and voice quality. We model variation information systematically from both explicit (speaker ID, language ID, pitch and duration) and implicit (utterance-level and phoneme-level prosody) perspectives: 1) For speaker and language ID, we use lookup embedding in training and inference; 2) For pitch and duration, we extract the values from paired text-speech data in training and use two predictors to predict the values in inference; 3) For utterance-level and phoneme-level prosody, we use two reference encoders to extract the values in training, and use two separate predictors to predict the values in inference. Additionally, we introduce an improved Conformer block to better model the local and global dependency in acoustic model. For task SH1, DelightfulTTS achieves 4.17 mean score in MOS test and 4.35 in SMOS test, which indicates the effectiveness of our proposed system


abstract: Face presentation attack detection (fPAD) plays a critical role in the modern face recognition pipeline. The generalization ability of face presentation attack detection models to unseen attacks has become a key issue for real-world deployment, which can be improved when models are trained with face images from different input distributions and different types of spoof attacks. In reality, due to legal and privacy issues, training data (both real face images and spoof images) are not allowed to be directly shared between different data sources. In this paper, to circumvent this challenge, we propose a Federated Test-Time Adaptive Face Presentation Attack Detection with Dual-Phase Privacy Preservation framework, with the aim of enhancing the generalization ability of fPAD models in both training and testing phase while preserving data privacy. In the training phase, the proposed framework exploits the federated learning technique, which simultaneously takes advantage of rich fPAD information available at different data sources by aggregating model updates from them without accessing their private data. To further boost the generalization ability, in the testing phase, we explore test-time adaptation by minimizing the entropy of fPAD model prediction on the testing data, which alleviates the domain gap between training and testing data and thus reduces the generalization error of a fPAD model. We introduce the experimental setting to evaluate the proposed framework and carry out extensive experiments to provide various insights about the proposed method for fPAD.


title: On the average hitting times of the squares of cycles

abstract: The exact formula for the average hitting time (HT, as an abbreviation) of simple random walks from one vertex to any other vertex on the square $C^2_N$ of an $N$-vertex cycle graph $C_N$ was given by N. Chair [\textit{Journal of Statistical Physics}, \textbf{154} (2014) 1177-1190]. In that paper, the author gives the expression for the even $N$ case and the expression for the odd $N$ case separately. In this paper, by using an elementary method different from Chair (2014), we give a much simpler single formula for the HT's of simple random walks on $C^2_N$. Our proof is considerably short and fully combinatorial, in particular, has no-need of any spectral graph theoretical arguments. Not only the formula itself but also intermediate results through the process of our proof describe clear relations between the HT's of simple random walks on $C^2_N$ and the Fibonacci numbers.


title: Linear Contextual Bandits with Adversarial Corruptions

abstract: We study the linear contextual bandit problem in the presence of adversarial corruption, where the interaction between the player and a possibly infinite decision set is contaminated by an adversary that can corrupt the reward up to a corruption level $C$ measured by the sum of the largest alteration on rewards in each round. We present a variance-aware algorithm that is adaptive to the level of adversarial contamination $C$. The key algorithmic design includes (1) a multi-level partition scheme of the observed data, (2) a cascade of confidence sets that are adaptive to the level of the corruption, and (3) a variance-aware confidence set construction that can take advantage of low-variance reward. We further prove that the regret of the proposed algorithm is $\tilde{O}(C^2d\sqrt{\sum_{t = 1}^T \sigma_t^2} + C^2R\sqrt{dT})$, where $d$ is the dimension of context vectors, $T$ is the number of rounds, $R$ is the range of noise and $\sigma_t^2,t=1\ldots,T$ are the variances of instantaneous reward. We also prove a gap-dependent regret bound for the proposed algorithm, which is instance-dependent and thus leads to better performance on good practical instances. To the best of our knowledge, this is the first variance-aware corruption-robust algorithm for contextual bandits. Experiments on synthetic data corroborate our theory.


title: On query complexity measures and their relations for symmetric functions

abstract: The main reason for query model's prominence in complexity theory and quantum computing is the presence of concrete lower bounding techniques: polynomial and adversary method. There have been considerable efforts to give lower bounds using these methods, and to compare/relate them with other measures based on the decision tree. We explore the value of these lower bounds on quantum query complexity and their relation with other decision tree based complexity measures for the class of symmetric functions, arguably one of the most natural and basic sets of Boolean functions. We show an explicit construction for the dual of the positive adversary method and also of the square root of private coin certificate game complexity for any total symmetric function. This shows that the two values can't be distinguished for any symmetric function. Additionally, we show that the recently introduced measure of spectral sensitivity gives the same value as both positive adversary and approximate degree for every total symmetric Boolean function. Further, we look at the quantum query complexity of Gap Majority, a partial symmetric function. It has gained importance recently in regard to understanding the composition of randomized query complexity. We characterize the quantum query complexity of Gap Majority and show a lower bound on noisy randomized query complexity (Ben-David and Blais, FOCS 2020) in terms of quantum query complexity. Finally, we study how large certificate complexity and block sensitivity can be as compared to sensitivity for symmetric functions (even up to constant factors). We show tight separations, i.e., give upper bounds on possible separations and construct functions achieving the same.


title: Black Hole-String Correspondence

abstract: I was asked to give a brief review of the black hole-string correspondence as a warm-up for a longer SITP-group discussion of a recent paper by Chen, Maldacena, and Witten. Here are my notes in written form.


title: Learning Insertion Primitives with Discrete-Continuous Hybrid Action Space for Robotic Assembly Tasks

abstract: This paper introduces a discrete-continuous action space to learn insertion primitives for robotic assembly tasks. Primitive is a sequence of elementary actions with certain exit conditions, such as "pushing down the peg until contact". Since the primitive is an abstraction of robot control commands and encodes human prior knowledge, it reduces the exploration difficulty and yields better learning efficiency. In this paper, we learn robot assembly skills via primitives. Specifically, we formulate insertion primitives as parameterized actions: hybrid actions consisting of discrete primitive types and continuous primitive parameters. Compared with the previous work using a set of discretized parameters for each primitive, the agent in our method can freely choose primitive parameters from a continuous space, which is more flexible and efficient. To learn these insertion primitives, we propose Twin-Smoothed Multi-pass Deep Q-Network (TS-MP-DQN), an advanced version of MP-DQN with twin Q-network to reduce the Q-value over-estimation. Extensive experiments are conducted in the simulation and real world for validation. From experiment results, our approach achieves higher success rates than three baselines: MP-DQN with parameterized actions, primitives with discrete parameters, and continuous velocity control. Furthermore, learned primitives are robust to sim-to-real transfer and can generalize to challenging assembly tasks such as tight round peg-hole and complex shaped electric connectors with promising success rates. Experiment videos are available at https://msc.berkeley.edu/research/insertion-primitives.html.


title: Resonance order-dependent plasmon-induced transparency in orthogonally-arranged nanoscale cavities

abstract: In this study, we investigate plasmon-induced transparency (PIT) in a resonator structure consisting of two orthogonally-arranged metal-insulator-metal (MIM) nanocavities with the aim of spectral modulation of a specific resonant order of the resonator. Our FDTD simulations demonstrate that when both cavities in this structure resonate at the same frequency, the PIT effect can be used to induce spectral modulation. This spectral modulation depends on the resonance order of the cavity coupled directly to the external field, occurring when first-order resonance is exhibited, but not with second-order resonance. We confirmed that this behavior is caused by the discrepancies between odd-order and even-order resonances using classical mechanical models analogous to the nanocavities. By tuning the resonance frequency and resonance order of the cavities, one can modulate the spectrum of the resonator structure in an order-selective manner.


title: Constraints on small-scale primordial density fluctuation from cosmic microwave background through dark matter annihilation

abstract: The cosmic microwave background (CMB) observation by the Planck satellite precisely determines primordial curvature fluctuations on larger scales than $\mathcal O(1)\,\mathrm{Mpc}$, while the small-scale curvature fluctuation is still less constrained. The constraint on small-scale fluctuations is highly improved if we assume the standard thermal relic dark matter scenario. When small-scale fluctuations are large enough, dense regions collapse to form small halos even in a redshift $z\gtrsim10^3$, which is called ultracompact minihalos. These minihalos enhance the annihilation of the dark matter and it is constrained by observations such as extragalactic gamma rays and the CMB. We revisit the effect of minihalos formed by the small-scale density fluctuations and calculate the ionization history modified by the dark matter annihilation. We perform the Markov Chain Monte Carlo method to constrain the size of small-scale curvature fluctuations by the CMB power spectrum. It is found that the constraint from the CMB power spectrum is comparable to that from the extragalactic gamma rays. We confirm that our constraint mainly comes from the energy injection in early time ($z\gtrsim 100$) and hence it is independent of the uncertainty of minihalo properties in the late time.


title: Accelerate 3D Object Processing via Spectral Layout

abstract: 3D image processing is an important problem in computer vision and pattern recognition fields. Compared with 2D image processing, its computation difficulty and cost are much higher due to the extra dimension. To fundamentally address this problem, we propose to embed the essential information in a 3D object into 2D space via spectral layout. Specifically, we construct a 3D adjacency graph to capture spatial structure of the 3D voxel grid. Then we calculate the eigenvectors corresponding to the second and third smallest eigenvalues of its graph Laplacian and perform spectral layout to map each voxel into a pixel in 2D Cartesian coordinate plane. The proposed method can achieve high quality 2D representations for 3D objects, which enables to use 2D-based methods to process 3D objects. The experimental results demonstrate the effectiveness and efficiency of our method.


title: Radical pairs may explain reactive oxygen species-mediated effects of hypomagnetic field on neurogenesis

abstract: Exposures to a hypomagnetic field can affect biological processes. Recently, it has been observed that hypomagnetic field exposure can adversely affect adult hippocampal neurogenesis and hippocampus-dependent cognition in mice. In the same study, the role of reactive oxygen species (ROS) in hypomagnetic field effects has been demonstrated. However, the mechanistic reasons behind this effect are not clear. This study proposes a radical pair mechanism based on a flavin-superoxide radical pair to explain the modulation of ROS production and the attenuation of adult hippocampal neurogenesis in a hypomagnetic field. The results of our calculations favor a singlet-born radical pair over a triplet-born radical pair. Our model predicts hypomagnetic field effects on the triplet/singlet yield of comparable strength as the effects observed in experimental studies on adult hippocampal neurogenesis. Our predictions are also in qualitative agreement with experimental results on superoxide concentration and other observed ROS effects. We also predict the effects of applied magnetic fields and oxygen isotopic substitution on adult hippocampal neurogenesis. Our findings strengthen the idea that nature might harness quantum resources in the context of the brain.


title: Ultra Light OCR Competition Technical Report

abstract: Ultra Light OCR Competition is a Chinese scene text recognition competition jointly organized by CSIG (China Society of Image and Graphics) and Baidu, Inc. In addition to focusing on common problems in Chinese scene text recognition, such as long text length and massive characters, we need to balance the trade-off of model scale and accuracy since the model size limitation in the competition is 10M. From experiments in aspects of data, model, training, etc, we proposed a general and effective method for Chinese scene text recognition, which got us second place among over 100 teams with accuracy 0.817 in TestB dataset. The code is available at https://aistudio.baidu.com/aistudio/projectdetail/2159102.


title: Adapting Quantum Approximation Optimization Algorithm (QAOA) for Unit Commitment

abstract: In the present Noisy Intermediate-Scale Quantum (NISQ), hybrid algorithms that leverage classical resources to reduce quantum costs are particularly appealing. We formulate and apply such a hybrid quantum-classical algorithm to a power system optimization problem called Unit Commitment, which aims to satisfy a target power load at minimal cost. Our algorithm extends the Quantum Approximation Optimization Algorithm (QAOA) with a classical minimizer in order to support mixed binary optimization. Using Qiskit, we simulate results for sample systems to validate the effectiveness of our approach. We also compare to purely classical methods. Our results indicate that classical solvers are effective for our simulated Unit Commitment instances with fewer than 400 power generation units. However, for larger problem instances, the classical solvers either scale exponentially in runtime or must resort to coarse approximations. Potential quantum advantage would require problem instances at this scale, with several hundred units.


title: A radiation tolerant clock generator for the CMS Endcap Timing Layer readout chip

abstract: We present the test results of a low jitter Phase Locked Loop (PLL) prototype chip for the CMS Endcap Timing Layer readout chip (ETROC). This chip is based on the improved version of a clock synthesis circuit named ljCDR from the Low-Power Gigabit Transceiver (lpGBT) project. The ljCDR is tested in its PLL mode. An automatic frequency calibration (AFC) block with the Triple Modular Redundancy (TMR) register is developed for the LC-oscillator calibration. The chip was manufactured in a 65 nm CMOS process with 10 metal layers. The chip has been extensively tested, including Total Ionizing Dose (TID) testing up to 300 Mrad and Single Event Upset (SEU) testing with heavy ions possessing a Linear energy transfer (LET) from 1.3 to 62.5 MeV*cm^2/mg.


title: QTIA, a 2.5 or 10 Gbps 4-Channel Array Optical Receiver ASIC in a 65 nm CMOS Technology

abstract: The Quad transimpedance and limiting amplifier (QTIA) is a 4-channel array optical receiver ASIC, developed using a 65 nm CMOS process. It is configurable between the bit rate of 2.56 Gbps and 10 Gbps per channel. QTIA offers careful matching to both GaAs and InGaAs photodiodes. At this R&D stage, each channel has a different biasing scheme to the photodiode for optimal coupling. A charge pump is implemented in one channel to provide a higher reverse bias voltage, which is especially important to mitigate radiation effects on the photodiodes. The circuit functions of QTIA successfully passed the lab tests with GaAs photodiodes.


title: Orthogonal variance-based feature selection for intrusion detection systems

abstract: In this paper, we apply a fusion machine learning method to construct an automatic intrusion detection system. Concretely, we employ the orthogonal variance decomposition technique to identify the relevant features in network traffic data. The selected features are used to build a deep neural network for intrusion detection. The proposed algorithm achieves 100% detection accuracy in identifying DDoS attacks. The test results indicate a great potential of the proposed method.


title: Recurrent Off-policy Baselines for Memory-based Continuous Control

abstract: When the environment is partially observable (PO), a deep reinforcement learning (RL) agent must learn a suitable temporal representation of the entire history in addition to a strategy to control. This problem is not novel, and there have been model-free and model-based algorithms proposed for this problem. However, inspired by recent success in model-free image-based RL, we noticed the absence of a model-free baseline for history-based RL that (1) uses full history and (2) incorporates recent advances in off-policy continuous control. Therefore, we implement recurrent versions of DDPG, TD3, and SAC (RDPG, RTD3, and RSAC) in this work, evaluate them on short-term and long-term PO domains, and investigate key design choices. Our experiments show that RDPG and RTD3 can surprisingly fail on some domains and that RSAC is the most reliable, reaching near-optimal performance on nearly all domains. However, one task that requires systematic exploration still proved to be difficult, even for RSAC. These results show that model-free RL can learn good temporal representation using only reward signals; the primary difficulty seems to be computational cost and exploration. To facilitate future research, we have made our PyTorch implementation publicly available at https://github.com/zhihanyang2022/off-policy-continuous-control.


title: Cylindric plane partitions, Lambda determinants, Commutants in semicircular systems

abstract: This thesis is divided into three parts. The first part deals with cylindric plane partitions. The second with lambda-determinants and the third with commutators in semi-circular systems. For more detailed abstract please see inside. Cylindric plane partitions may be thought of as a natural generalization of reverse plane partitions. A generating series for the enumeration of cylindric plane partitions was recently given by Borodin. The first result of section one is a new bijective proof of Borodin's identity which makes use of Fomin's growth diagram framework for generalized RSK correspondences. The second result is a $(q,t)$-analog of Borodin's identity which extends previous work by Okada in the reverse plane partition case. The third result is an explicit combinatorial interpretation of the Macdonald weight occurring in the $(q,t)$-analog using the non-intersecting lattice path model for cylindric plane partitions. Alternating sign matrices were discovered by Robbins and Rumsey whilst studying \emph{$\lambda$-determinants}. In the second part of this thesis we prove a multi-parameter generalization of the $\lambda$-determinant, generalizing a recent result by di Francesco. Like the original $\lambda$-determinant, our formula exhibits the \emph{Laurent phenomenon}. Semicircular systems were first introduced by Voiculescu as a part of his study of von Neumann algebras. In the third part of this thesis we study certain commutator subalgebras of the semicircular system. We find a projection matrix with an interesting self-similar structure. Making use of our projection formula we given an alternative, elementary proof that the semicircular system is a factor.


title: Unified Smoothing Approach for Best Hyperparameter Selection Problem Using a Bilevel Optimization Strategy

abstract: Strongly motivated from use in various fields including machine learning, the methodology of sparse optimization has been developed intensively so far. Especially, the recent advance of algorithms for solving problems with nonsmooth regularizers is remarkable. However, those algorithms suppose that weight parameters of regularizers, called hyperparameters hereafter, are pre-fixed, and it is a crucial matter how the best hyperparameter should be selected. In this paper, we focus on the hyperparameter selection of regularizers related to the $\ell_p$ function with $0<p\le 1$ and apply a bilevel programming strategy, wherein we need to solve a bilevel problem, whose lower-level problem is nonsmooth, possibly nonconvex and non-Lipschitz. Recently, for solving a bilevel problem for hyperparameter selection of the pure $\ell_p\ (0<p \le 1)$ regularizer Okuno et al. discovered new necessary optimality conditions, called SB(scaled bilevel)-KKT conditions, and further proposed a smoothing-type algorithm using a certain smoothing function. However, this optimality measure is loose in the sense that there could be many points that satisfy the SB-KKT conditions. In this work, we propose new bilevel KKT conditions, which are new necessary optimality conditions tighter than the ones proposed by Okuno et al. Moreover, we propose a unified smoothing approach using smoothing functions that belong to the Chen-Mangasarian class, and then prove that generated iteration points accumulate at \alert{bilevel KKT points under milder constraint qualifications. Another contribution is that our approach and analysis are applicable to a wider class of regularizers. Numerical comparisons demonstrate which smoothing functions work well for hyperparameter optimization via bilevel optimization approach.


title: Time series signal recovery methods: comparative study

abstract: Signal data often contains missing values. Effective replacement (imputation) of the missing values can have significant positive effects on processing the signal. In this paper, we compare three commonly employed methods for estimating missing values in time series data: forward fill, backward fill, and mean fill. We carry out a large scale experimental analysis using 3,600 AR(1)-based simulated time series to determine the optimal method for estimating missing values. The results of the numerical experiments show that the forward and backward fill methods are better suited for times series with large positive correlations, while the mean fill method is better suited for times series with low or negative correlations. The extensive and exhaustive nature of the numerical experiments provides a definitive answer to the comparison of the three imputation methods.


title: Effects of shallow carbon and deep N++ layer on the radiation hardness of IHEP-IME LGAD sensors

abstract: Low Gain Avalanche Diode (LGAD) is applied for the High-Granularity Timing Detector (HGTD), and it will be used to upgrade the ATLAS experiment. The first batch IHEP-IME LGAD sensors were designed by the Institute of High Energy Physics (IHEP) and fabricated by the Institute of Microelectronics (IME). Three IHEP-IME sensors (W1, W7 and W8) were irradiated by the neutrons up to the fluence of 2.5 x 10^15 n_eq/cm^2 to study the effect of the shallow carbon and deep N++ layer on the irradiation hardness. Taking W7 as a reference, W1 has an extra shallow carbon applied, and W8 has a deeper N++ layer.


title: Age and Gender Prediction using Deep CNNs and Transfer Learning

abstract: The last decade or two has witnessed a boom of images. With the increasing ubiquity of cameras and with the advent of selfies, the number of facial images available in the world has skyrocketed. Consequently, there has been a growing interest in automatic age and gender prediction of a person using facial images. We in this paper focus on this challenging problem. Specifically, this paper focuses on age estimation, age classification and gender classification from still facial images of an individual. We train different models for each problem and we also draw comparisons between building a custom CNN (Convolutional Neural Network) architecture and using various CNN architectures as feature extractors, namely VGG16 pre-trained on VGGFace, Res-Net50 and SE-ResNet50 pre-trained on VGGFace2 dataset and training over those extracted features. We also provide baseline performance of various machine learning algorithms on the feature extraction which gave us the best results. It was observed that even simple linear regression trained on such extracted features outperformed training CNN, ResNet50 and ResNeXt50 from scratch for age estimation.


title: Accelerated Almost-Sure Convergence Rates for Nonconvex Stochastic Gradient Descent using Stochastic Learning Rates

abstract: Large-scale optimization problems require algorithms both effective and efficient. One such popular and proven algorithm is Stochastic Gradient Descent which uses first-order gradient information to solve these problems. This paper studies almost-sure convergence rates of the Stochastic Gradient Descent method when instead of deterministic, its learning rate becomes stochastic. In particular, its learning rate is equipped with a multiplicative stochasticity, producing a stochastic learning rate scheme. Theoretical results show accelerated almost-sure convergence rates of Stochastic Gradient Descent in a nonconvex setting when using an appropriate stochastic learning rate, compared to a deterministic-learning-rate scheme. The theoretical results are verified empirically.


title: Progressively Select and Reject Pseudo-labelled Samples for Open-Set Domain Adaptation

abstract: Domain adaptation solves image classification problems in the target domain by taking advantage of the labelled source data and unlabelled target data. Usually, the source and target domains share the same set of classes. As a special case, Open-Set Domain Adaptation (OSDA) assumes there exist additional classes in the target domain but not present in the source domain. To solve such a domain adaptation problem, our proposed method learns discriminative common subspaces for the source and target domains using a novel Open-Set Locality Preserving Projection (OSLPP) algorithm. The source and target domain data are aligned in the learned common spaces class-wisely. To handle the open-set classification problem, our method progressively selects target samples to be pseudo-labelled as known classes and rejects the outliers if they are detected as from unknown classes. The common subspace learning algorithm OSLPP simultaneously aligns the labelled source data and pseudo-labelled target data from known classes and pushes the rejected target data away from the known classes. The common subspace learning and the pseudo-labelled sample selection/rejection facilitate each other in an iterative learning framework and achieves state-of-the-art performance on benchmark datasets Office-31 and Office-Home with the average HOS of 87.4% and 67.0% respectively.


title: MOVER confidence intervals for a difference or ratio effect parameter under stratified sampling

abstract: Stratification is commonly employed in clinical trials to reduce the chance covariate imbalances and increase the precision of the treatment effect estimate. We propose a general framework for constructing the confidence interval (CI) for a difference or ratio effect parameter under stratified sampling by the method of variance estimates recovery (MOVER). We consider the additive variance and additive CI approaches for the difference, in which either the CI for the weighted difference, or the CI for the weighted effect in each group, or the variance for the weighted difference is calculated as the weighted sum of the corresponding stratum-specific statistics. The CI for the ratio is derived by the Fieller and log-ratio methods. The weights can be random quantities under the assumption of a constant effect across strata, but this assumption is not needed for fixed weights. These methods can be easily applied to different endpoints in that they require only the point estimate, CI, and variance estimate for the measure of interest in each group across strata. The methods are illustrated with two real examples. In one example, we derive the MOVER CIs for the risk difference and risk ratio for binary outcomes. In the other example, we compare the restricted mean survival time and milestone survival in stratified analysis of time-to-event outcomes. Simulations show that the proposed MOVER CIs generally outperform the standard large sample CIs, and that the additive CI approach performs slightly better than the additive variance approach.


title: Thermodynamics of spin crossover in ferropericlase: an improved LDA+$U_{sc}$ calculation

abstract: We present LDA+$U_{sc}$ calculations of high-spin (HS) and low-spin (LS) states in ferropericlase (fp) with an iron concentration of 18.75$\%$. The Hubbard parameter $U$ is determined self-consistently with structures optimized at arbitrary pressures. We confirm a strong dependence of $U$ on the pressure and spin state. Static calculations confirm that the antiferromagnetic configuration is more stable than the ferromagnetic one in the HS state, consistent with low-temperature measurements. Phonon calculations guarantee the dynamical stability of HS and LS states throughout the pressure range of the Earth mantle. Compression curves for HS and LS states agree well with experiments. Using a non-ideal mixing model for the HS to LS states solid solution, we obtain a crossover starting at $\sim$45 GPa at room temperature and considerably broader than previous results. The spin-crossover phase diagram is calculated, including vibrational, magnetic, electronic, and non-ideal HS-LS entropic contributions. Our results suggest the mixed-spin state predominates in fp in most of the lower mantle.


title: Deep Learning for UAV-based Object Detection and Tracking: A Survey

abstract: Owing to effective and flexible data acquisition, unmanned aerial vehicle (UAV) has recently become a hotspot across the fields of computer vision (CV) and remote sensing (RS). Inspired by recent success of deep learning (DL), many advanced object detection and tracking approaches have been widely applied to various UAV-related tasks, such as environmental monitoring, precision agriculture, traffic management. This paper provides a comprehensive survey on the research progress and prospects of DL-based UAV object detection and tracking methods. More specifically, we first outline the challenges, statistics of existing methods, and provide solutions from the perspectives of DL-based models in three research topics: object detection from the image, object detection from the video, and object tracking from the video. Open datasets related to UAV-dominated object detection and tracking are exhausted, and four benchmark datasets are employed for performance evaluation using some state-of-the-art methods. Finally, prospects and considerations for the future work are discussed and summarized. It is expected that this survey can facilitate those researchers who come from remote sensing field with an overview of DL-based UAV object detection and tracking methods, along with some thoughts on their further developments.


title: Spatial profile of accelerated electrons from ponderomotive scattering in hydrogen cluster targets

abstract: We study the laser-driven acceleration of electrons from overdense hydrogen clusters to energies of up to 13 MeV in laser forward direction and several hundreds of keV in an outer ring-like structure. The use of cryogenic hydrogen allows for high repetition-rate operation and examination of the influence of source parameters like temperature and gas flow. The outer ring-like structure of accelerated electrons, originating from the interaction, that is robust against the change of laser and target parameters can be observed for low electron densities of ca. 3$\times$10$^{16}$ cm$^{-3}$. For higher electron densities, an additional central spot of electrons in the laser forward direction can be observed. Utilizing 3D-PIC simulations, it is revealed that both electron populations mainly stem from ponderomotive scattering.


title: A sufficient condition for the quasipotential to be the rate function of the invariant measure of countable-state mean-field interacting particle systems

abstract: This paper considers the family of invariant measures of Markovian mean-field interacting particle systems on a countably infinite state space and studies its large deviation asymptotics. The Freidlin-Wentzell quasipotential is the usual candidate rate function for the sequence of invariant measures indexed by the number of particles. The paper provides two counterexamples where the quasipotential is not the rate function. The quasipotential arises from finite horizon considerations. However there are certain barriers that cannot be surmounted easily in any finite time horizon, but these barriers can be crossed in the stationary regime. Consequently, the quasipotential is infinite at some points where the rate function is finite. After highlighting this phenomenon, the paper studies some sufficient conditions on a class of interacting particle systems under which one can continue to assert that the Freidlin-Wentzell quasipotential is indeed the rate function.


title: Freeze-frame approach for robust single-molecule tip-enhnaced Raman spectroscopy at room temperature

abstract: A quantitative single-molecule tip-enhanced Raman spectroscopy (TERS) study at room temperature remained a challenge due to the rapid structural dynamics of molecules exposed to air. Here, we demonstrate the single-molecule level hyperspectral TERS imaging of brilliant cresyl blue (BCB) at room temperature for the first time, along with quantitative spectral analyses. Freeze-frame approach using a thin Al2O3 capping layer, which suppresses spectral diffusions and inhibits chemical reactions and contaminations in air, enabled reliable and robust chemical imaging. For the molecules resolved spatially in the TERS image, a clear Raman peak variation up to 7.5 cm-1 is observed, which cannot be found in molecular ensembles. From density functional theory-based quantitative analyses of the varied TERS peaks, we reveal the conformational heterogeneity at the single-molecule level. This work provides a facile way to investigate the single-molecule properties in interacting media, expanding the scope of single-molecule vibrational spectroscopy.


title: Exploration of a new reconstructed structure on GaN(0001) surface by Bayesian optimization

abstract: GaN(0001) surfaces with Ga- and H-adsorbates are fundamental stages for epitaxial growth of semiconductor thin films. We explore stable surface structures with nanometer scale by the density-functional calculations combined with Bayesian optimization, and succeed to reach a single structure with satisfactorily low mixing enthalpy among hundreds of thousand possible candidate structures. We find that the obtained structure is free from any postulated high symmetry previously introduced by human intuition, satisfies electron counting rule locally, and shows new adsorbate arrangement, reflecting characteristics of nitride semiconductors.


title: A diophantine problem concerning third order matrices

abstract: In this paper we find a third order unimodular matrix, none of whose entries is $1$ or $-1$, such that when each entry of the matrix is replaced by its cube, the resulting matrix is also unimodular. Further, we find third order square integer matrices $(a_{ij})$, none of the integers $a_{ij}$ being $1$ or $-1$, such that $\det{(a_{ij})}=k$ and $\det{(a_{ij}^3)}=k^3$, where $k$ is a nonzero integer.


title: Kernel density estimation-based sampling for neural network classification

abstract: Imbalanced data occurs in a wide range of scenarios. The skewed distribution of the target variable elicits bias in machine learning algorithms. One of the popular methods to combat imbalanced data is to artificially balance the data through resampling. In this paper, we compare the efficacy of a recently proposed kernel density estimation (KDE) sampling technique in the context of artificial neural networks. We benchmark the KDE sampling method against two base sampling techniques and perform comparative experiments using 8 datasets and 3 neural networks architectures. The results show that KDE sampling produces the best performance on 6 out of 8 datasets. However, it must be used with caution on image datasets. We conclude that KDE sampling is capable of significantly improving the performance of neural networks.


title: SgSum: Transforming Multi-document Summarization into Sub-graph Selection

abstract: Most of existing extractive multi-document summarization (MDS) methods score each sentence individually and extract salient sentences one by one to compose a summary, which have two main drawbacks: (1) neglecting both the intra and cross-document relations between sentences; (2) neglecting the coherence and conciseness of the whole summary. In this paper, we propose a novel MDS framework (SgSum) to formulate the MDS task as a sub-graph selection problem, in which source documents are regarded as a relation graph of sentences (e.g., similarity graph or discourse graph) and the candidate summaries are its sub-graphs. Instead of selecting salient sentences, SgSum selects a salient sub-graph from the relation graph as the summary. Comparing with traditional methods, our method has two main advantages: (1) the relations between sentences are captured by modeling both the graph structure of the whole document set and the candidate sub-graphs; (2) directly outputs an integrate summary in the form of sub-graph which is more informative and coherent. Extensive experiments on MultiNews and DUC datasets show that our proposed method brings substantial improvements over several strong baselines. Human evaluation results also demonstrate that our model can produce significantly more coherent and informative summaries compared with traditional MDS methods. Moreover, the proposed architecture has strong transfer ability from single to multi-document input, which can reduce the resource bottleneck in MDS tasks. Our code and results are available at: \url{https://github.com/PaddlePaddle/Research/tree/master/NLP/EMNLP2021-SgSum}.


title: Zero-Shot Dialogue Disentanglement by Self-Supervised Entangled Response Selection

abstract: Dialogue disentanglement aims to group utterances in a long and multi-participant dialogue into threads. This is useful for discourse analysis and downstream applications such as dialogue response selection, where it can be the first step to construct a clean context/response set. Unfortunately, labeling all~\emph{reply-to} links takes quadratic effort w.r.t the number of utterances: an annotator must check all preceding utterances to identify the one to which the current utterance is a reply. In this paper, we are the first to propose a~\textbf{zero-shot} dialogue disentanglement solution. Firstly, we train a model on a multi-participant response selection dataset harvested from the web which is not annotated; we then apply the trained model to perform zero-shot dialogue disentanglement. Without any labeled data, our model can achieve a cluster F1 score of 25. We also fine-tune the model using various amounts of labeled data. Experiments show that with only 10\% of the data, we achieve nearly the same performance of using the full dataset\footnote{Code is released at \url{https://github.com/chijames/zero_shot_dialogue_disentanglement}}.


title: Bone Marrow Cell Recognition: Training Deep Object Detection with A New Loss Function

abstract: For a long time, bone marrow cell morphology examination has been an essential tool for diagnosing blood diseases. However, it is still mainly dependent on the subjective diagnosis of experienced doctors, and there is no objective quantitative standard. Therefore, it is crucial to study a robust bone marrow cell detection algorithm for a quantitative automatic analysis system. Currently, due to the dense distribution of cells in the bone marrow smear and the diverse cell classes, the detection of bone marrow cells is difficult. The existing bone marrow cell detection algorithms are still insufficient for the automatic analysis system of bone marrow smears. This paper proposes a bone marrow cell detection algorithm based on the YOLOv5 network, trained by minimizing a novel loss function. The classification method of bone marrow cell detection tasks is the basis of the proposed novel loss function. Since bone marrow cells are classified according to series and stages, part of the classes in adjacent stages are similar. The proposed novel loss function considers the similarity between bone marrow cell classes, increases the penalty for prediction errors between dissimilar classes, and reduces the penalty for prediction errors between similar classes. The results show that the proposed loss function effectively improves the algorithm's performance, and the proposed bone marrow cell detection algorithm has achieved better performance than other cell detection algorithms.


title: Review-Based Domain Disentanglement without Duplicate Users or Contexts for Cross-Domain Recommendation

abstract: A cross-domain recommendation has shown promising results in solving data-sparsity and cold-start problems. Despite such progress, existing methods focus on domain-shareable information (overlapped users or same contexts) for a knowledge transfer, and they fail to generalize well without such requirements. To deal with these problems, we suggest utilizing review texts that are general to most e-commerce systems. Our model (named SER) uses three text analysis modules, guided by a single domain discriminator for disentangled representation learning. Here, we suggest a novel optimization strategy that can enhance the quality of domain disentanglement, and also debilitates detrimental information of a source domain. Also, we extend the encoding network from a single to multiple domains, which has proven to be powerful for review-based recommender systems. Extensive experiments and ablation studies demonstrate that our method is efficient, robust, and scalable compared to the state-of-the-art single and cross-domain recommendation methods.


title: BFKL phenomenology: resummation of high-energy logs in inclusive processes

abstract: We present recent phenomenological studies, tailored on kinematic configurations typical of current and forthcoming analyses at the LHC, for two novel probe channels of the BFKL resummation of energy logarithms. Particular attention is drawn to the behavior of distributions differential in azimuthal angle and rapidity, where significant high-energy effects are expected.


title: Sparser Kernel Herding with Pairwise Conditional Gradients without Swap Steps

abstract: The Pairwise Conditional Gradients (PCG) algorithm is a powerful extension of the Frank-Wolfe algorithm leading to particularly sparse solutions, which makes PCG very appealing for problems such as sparse signal recovery, sparse regression, and kernel herding. Unfortunately, PCG exhibits so-called swap steps that might not provide sufficient primal progress. The number of these bad steps is bounded by a function in the dimension and as such known guarantees do not generalize to the infinite-dimensional case, which would be needed for kernel herding. We propose a new variant of PCG, the so-called Blended Pairwise Conditional Gradients (BPCG). This new algorithm does not exhibit any swap steps, is very easy to implement, and does not require any internal gradient alignment procedures. The convergence rate of BPCG is basically that of PCG if no drop steps would occur and as such is no worse than PCG but improves and provides new rates in many cases. Moreover, we observe in the numerical experiments that BPCG's solutions are much sparser than those of PCG. We apply BPCG to the kernel herding setting, where we derive nice quadrature rules and provide numerical results demonstrating the performance of our method.


title: Doping evolution of superconductivity, charge order and band topology in hole-doped topological kagome superconductors Cs(V$_{1-x}$Ti$_x$)$_3$Sb$_5$

abstract: The newly discovered Kagome superconductors $A$V$_3$Sb$_5$ ($A$ = K, Rb, Cs) exhibit superconductivity, charge order, and band topology simultaneously. To explore the intricate interplay between the superconducting and charge orders, we investigate the doping evolution of superconductivity, charge-density-wave (CDW) order, and band topology in doped topological kagome superconductors Cs(V$_{1-x}$Ti$_x$)$_3$Sb$_5$ where the Ti-dopant introduces hole-like charge carriers. Despite the absence of the CDW phase transition in doped compounds even for the lowest doping level of $x=0.047$, the superconductivity survives in all doped samples with enhanced critical temperatures. The high-resolution angle-resolved photoemission spectroscopy (ARPES) measurements reveal that the Ti-dopant in the kagome plane lowers the chemical potential, pushing the van Hove singularity (VHS) at $M$ point above the Fermi level. First-principle simulations corroborate the doping evolution of the band structure observed in ARPES, and affirm that the CDW instability does not occur once the VHS is pushed above the Fermi level, explaining the absence of the CDW ordering in our doped samples Cs(V$_{1-x}$Ti$_x$)$_3$Sb$_5$. Our results demonstrate a competition between the CDW and superconducting orders in the kagome-metal superconductor CsV$_3$Sb$_5$, although the superconductivity is likely inconsequential of the CDW order.


title: Extractors for Sum of Two Sources

abstract: We consider the problem of extracting randomness from \textit{sumset sources}, a general class of weak sources introduced by Chattopadhyay and Li (STOC, 2016). An $(n,k,C)$-sumset source $\mathbf{X}$ is a distribution on $\{0,1\}^n$ of the form $\mathbf{X}_1 + \mathbf{X}_2 + \ldots + \mathbf{X}_C$, where $\mathbf{X}_i$'s are independent sources on $n$ bits with min-entropy at least $k$. Prior extractors either required the number of sources $C$ to be a large constant or the min-entropy $k$ to be at least $0.51 n$. As our main result, we construct an explicit extractor for sumset sources in the setting of $C=2$ for min-entropy $\mathrm{poly}(\log n)$ and polynomially small error. We can further improve the min-entropy requirement to $(\log n) \cdot (\log \log n)^{1 + o(1)}$ at the expense of worse error parameter of our extractor. We find applications of our sumset extractor for extracting randomness from other well-studied models of weak sources such as affine sources, small-space sources, and interleaved sources. Interestingly, it is unknown if a random function is an extractor for sumset sources. We use techniques from additive combinatorics to show that it is a disperser, and further prove that an affine extractor works for an interesting subclass of sumset sources which informally corresponds to the "low doubling" case (i.e., the support of $\mathbf{X_1} + \mathbf{X_2}$ is not much larger than $2^k$).


title: Index of Embedded Networks in the Sphere

abstract: In this paper, we will compute the Morse index and nullity for the stationary embedded networks in spheres. The key theorem in the computation is that the index (and nullity) for the whole network is related to the index (and nullity) of small networks and the Dirichlet-to-Neumann map defined in this paper. Finally, we will show that for all stationary triple junction networks in $\mathbb{S}^2$, there is only one eigenvalue (without multiplicity) $-1$, which is less than 0, and the corresponding eigenfunctions are locally constant. Besides, the multiplicity of eigenvalues 0 is 3 for these networks, and their eigenfunctions are generated by the rotations on the sphere.


title: Facilitating Database Tuning with Hyper-Parameter Optimization: A Comprehensive Experimental Evaluation

abstract: Recently, using automatic configuration tuning to improve the performance of modern database management systems (DBMSs) has attracted increasing interest from the database community. This is embodied with a number of systems featuring advanced tuning capabilities being developed. However, it remains a challenge to select the best solution for database configuration tuning, considering the large body of algorithm choices. In addition, beyond the applications on database systems, we could find more potential algorithms designed for configuration tuning. To this end, this paper provides a comprehensive evaluation of configuration tuning techniques from a broader perspective, hoping to better benefit the database community. In particular, we summarize three key modules of database configuration tuning systems and conduct extensive ablation studies using various challenging cases. Our evaluation demonstrates that the hyper-parameter optimization algorithms can be borrowed to further enhance the database configuration tuning. Moreover, we identify the best algorithm choices for different modules. Beyond the comprehensive evaluations, we offer an efficient and unified database configuration tuning benchmark via surrogates that reduces the evaluation cost to a minimum, allowing for extensive runs and analysis of new techniques.


title: Search for sterile neutrinos by shower events at a future neutrino telescope

abstract: It is pointed out that searching for sterile neutrinos at high energy regions at a future IceCube-like facility has advantages compared with that of reactor or short baseline accelerator neutrino experiments, in which the size of the detector and energy resolution make it difficult to gain good sensitivity for a large mass-squared difference. In this study we show that it is possible to improve constraints on sterile neutrino mixing $\theta_{14}$ for 1 eV$^2~\lesssim \Delta m_{41}^2\lesssim$ 100 eV$^2$ by looking for a dip in the shower events at an IceCube-like neutrino telescope whose volume is at least 10 times as large as that of IceCube and duration is 10 years. We also give an analytic expression for the oscillation probabilities in two cases where the condition of one mass scale dominance is satisfied.


title: Uniformization of surfaces with boundary and the application to the triple junction surfaces with negative Euler characteristic

abstract: The conformal structure on minimal surfaces plays a key role in studying the properties of minimal surfaces. Here we extend the results of uniformization of surfaces with boundary to get the (weak) uniformization results for triple junction surfaces.


title: Nonlinear electronic density response of the warm dense electron gas: multiple perturbations and mode coupling

abstract: We present extensive new ab initio path integral Monte Carlo (PIMC) results for an electron gas at warm dense matter conditions that is subject to multiple harmonic perturbations. In addition to the previously investigated nonlinear effects at the original wave number [Dornheim \emph{et al.}, PRL \textbf{125}, 085001 (2020)] and the excitation of higher harmonics [Dornheim \emph{et al.}, PRR \textbf{3}, 033231 (2021)], the presence of multiple external potentials leads to mode-coupling effects, which constitute the dominant nonlinear effect and lead to a substantially more complicated density response compared to linear response theory. One possibility to estimate mode-coupling effects from a PIMC simulation of the unperturbed system is given in terms of generalized imaginary-time correlation functions that have been recently introduced by Dornheim \emph{et al.}~[JCP \textbf{155}, 054110 (2021)]. In addition, we extend our previous analytical theory of the nonlinear density response of the electron gas in terms of the static local field correction [Dornheim \emph{et al.}, PRL \textbf{125}, 235001 (2020)], which allows for a highly accurate description of the PIMC results with negligible computational cost.


title: Operator Shifting for Model-based Policy Evaluation

abstract: In model-based reinforcement learning, the transition matrix and reward vector are often estimated from random samples subject to noise. Even if the estimated model is an unbiased estimate of the true underlying model, the value function computed from the estimated model is biased. We introduce an operator shifting method for reducing the error introduced by the estimated model. When the error is in the residual norm, we prove that the shifting factor is always positive and upper bounded by $1+O\left(1/n\right)$, where $n$ is the number of samples used in learning each row of the transition matrix. We also propose a practical numerical algorithm for implementing the operator shifting.

categories: cs.LG cs.NA math.NA math.ST stat.ML stat.TH


title: Interactions between irregular wave fields and sea ice: A physical model for wave attenuation and ice breakup in an ice tank

abstract: Irregular, unidirectional surface water waves incident on model ice in an ice tank are used as a physical model of ocean surface wave interactions with sea ice. Results are given for an experiment consisting of three tests, starting with a continuous ice cover and in which the incident wave steepness increases between tests. The incident waves range from causing no breakup of the ice cover to breakup of the full length of ice cover. Temporal evolution of the ice edge, breaking front and mean floe sizes are reported. Floe size distributions in the different tests are analysed. The evolution of the wave spectrum with distance into the ice-covered water is analysed in terms of changes of energy content, mean wave period and spectral bandwidth relative to their incident counterparts, and pronounced differences are found between the tests. Further, an empirical attenuation coefficient is derived from the measurements and shown to have a power-law dependence on frequency comparable to that found in field measurements. Links between wave properties and ice breakup are discussed.


title: Metal pollution of the solar white dwarf by solar system small bodies

abstract: White dwarfs (WDs) often show metal lines in their spectra, indicating accretion of asteroidal material. Our Sun is to become a WD in several Gyr. Here, we examine how the solar WD accretes from the three major small body populations: the main belt asteroids (MBAs), Jovian trojan asteroids (JTAs), and trans-Neptunian objects (TNOs). Owing to the solar mass loss during the giant branch, 40\% of the JTAs are lost but the vast majority of MBAs and TNOs survive. During the WD phase, objects from all three populations are sporadically scattered onto the WD, implying ongoing accretion. For young cooling ages $\lesssim 100$ Myr, accretion of MBAs predominates; our predicted accretion rate $\sim10^6$ g/s falls short of observations by two orders of magnitude. On Gyr timescales, thanks to the consumption of the TNOs that kicks in $\gtrsim 100$ Myr, the rate oscillates around $10^6-10^7$ g/s until several Gyr and drops to $\sim10^5$ g/s at 10 Gyr. Our solar WD accretion rate from 1 Gyr and beyond agrees well with those of the extrasolar WDs. We show that for the solar WD, the accretion source region evolves in an inside-out pattern. Moreover, in a realistic small body population with individual sizes covering a wide range as WD pollutants, the accretion is dictated by the largest objects. As a consequence, the accretion rate is lower by an order of magnitude than that from a population of bodies of a uniform size and the same total mass and shows greater scatter.