title: The climate in climate economics

abstract: To analyze climate change mitigation strategies, economists rely on simplified climate models - climate emulators. We propose a generic and transparent calibration and evaluation strategy for these climate emulators that is based on Coupled Model Intercomparison Project, Phase 5 (CMIP5). We demonstrate that the appropriate choice of the free model parameters can be of key relevance for the predicted social cost of carbon. We propose to use four different test cases: two tests to separately calibrate and evaluate the carbon cycle and temperature response, a test to quantify the transient climate response, and a final test to evaluate the performance for scenarios close to those arising from economic models. We re-calibrate the climate part of the widely used DICE-2016: the multi-model mean as well as extreme, but still permissible climate sensitivities and carbon cycle responses. We demonstrate that the functional form of the climate emulator of the DICE-2016 model is fit for purpose, despite its simplicity, but its carbon cycle and temperature equations are miscalibrated. We examine the importance of the calibration for the social cost of carbon in the context of a partial equilibrium setting where interest rates are exogenous, as well as the simple general equilibrium setting from DICE-2016. We find that the model uncertainty from different consistent calibrations of the climate system can change the social cost of carbon by a factor of four if one assumes a quadratic damage function. When calibrated to the multi-model mean, our model predicts similar values for the social cost of carbon as the original DICE-2016, but with a strongly reduced sensitivity to the discount rate and about one degree less long-term warming. The social cost of carbon in DICE-2016 is oversensitive to the discount rate, leading to extreme comparative statics responses to changes in preferences.


title: Hunt's hypothesis and symmetrization for general 1-dimensional diffusions

abstract: In this paper, we will consider the problem that how far from Hunt's hypothesis (H) to symmetrization for a general 1-dimensional diffusion. A characterization of (H) involving the classification of points for this diffusion will be first obtained. Then the main result shows that such a process is symmetrizable, if and only if (H) holds and a certain family of asymmetric shunt points is empty. Furthermore, we will also derive the representation of associated Dirichlet forms of general 1-dimensional diffusions under symmetry.


title: MaxwellNet: Physics-driven deep neural network training based on Maxwell's equations

abstract: Maxwell's equations govern light propagation and its interaction with matter. Therefore, the solution of Maxwell's equations using computational electromagnetic simulations plays a critical role in understanding light-matter interaction and designing optical elements. Such simulations are often time-consuming and recent activities have been described to replace or supplement them with trained deep neural networks (DNNs). Such DNNs typically require extensive, computationally demanding simulations using conventional electromagnetic solvers to compose the training dataset. In this paper, we present a novel scheme to train a DNN that solves Maxwell's equations speedily and accurately without relying on other computational electromagnetic solvers. Our approach is to train a DNN using the residual of Maxwell's equations as the physics-driven loss function for a network that finds the electric field given the spatial distribution of the material property. We demonstrate it by training a single network that simultaneously finds multiple solutions of various aspheric micro-lenses. Furthermore, we exploit the speed of this network in a novel inverse design scheme to design a micro-lens that maximizes a desired merit function. We believe that our approach opens up a novel way for light simulation and optical design of photonics devices.


title: 3D Parametric Wireframe Extraction Based on Distance Fields

abstract: We present a pipeline for parametric wireframe extraction from densely sampled point clouds. Our approach processes a scalar distance field that represents proximity to the nearest sharp feature curve. In intermediate stages, it detects corners, constructs curve segmentation, and builds a topological graph fitted to the wireframe. As an output, we produce parametric spline curves that can be edited and sampled arbitrarily. We evaluate our method on 50 complex 3D shapes and compare it to the novel deep learning-based technique, demonstrating superior quality.


title: Hybrid Beamforming Design for Wideband mmWave Full-Duplex Systems

abstract: Recently, full duplex (FD) has been studied in 5G LTE millimeter wave (mmWave) cellular communications for New Radio in 3GPP releases 15-17. FD allows bidirectional transmission over the same resources and has the potential to reduce latency and double spectral efficiency. Self-interference (SI) is the primary drawback. SI can be several orders of magnitude greater than the received signal power, saturate the analog-to-digital converters (ADCs) and degrade communication performance severely. Massive mmWave antenna arrays may provide enough degrees of freedom for spatial multiplexing and SI suppression. In this paper, we design spatial beamformers for the phased arrays already built into the FD basestation/relay to extend mmWave coverage to a single user. We propose alternating projections to design the precoder and combiner to maximize the sum of the uplink and downlink spectral efficiencies while bringing SI below the noise floor. Our contributions include (1) hybrid analog/digital beamformer design algorithm to cancel SI in the analog domain to avoid ADC saturation and in the digital domain on each subcarrier; (2) full-digital beamformer design algorithm; and (3) analysis of spectral efficiency, energy efficiency and outage probability. In simulation, the proposed algorithms outperform beamsteering, singular value decomposition, angle search, and half-duplex techniques.


title: Deep Learning Assisted Compact Modeling of Nanoscale Transistor

abstract: Transistors are the basic building blocks for all electronics. Accurate prediction of their current-voltage (IV) characteristics enables circuit simulations before the expensive silicon tape-out. In this work, we propose using deep neural network to improve the accuracy for the conventional, physics-based compact model for nanoscale transistors. Physics-driven requirements on the neural network are discussed. Using finite element simulation as the input dataset, together with a neural network with roughly 30 neurons, the final IV model can well-predict the IV to within 1%. The trained model can readily be implemented by the hardware description language (HDL) such as VerilogA for circuit simulation.


title: Reactivity of the Si(100)-2$\times$1-Cl surface with respect to PH$_3$, PCl$_3$, and BCl$_3$: Comparison with PH$_3$ on Si(100)-2$\times$1-H

abstract: Despite the interest in a chlorine monolayer on Si(100) as an alternative to hydrogen resist for atomic-precision doping, little is known about its interaction with dopant-containing molecules. We used the density functional theory to evaluate whether a chlorine monolayer on Si(100) is suitable as a resist for \ce{PH3}, \ce{PCl3}, and \ce{BCl3} molecules. We calculated reaction pathways for \ce{PH3}, \ce{PCl3}, and \ce{BCl3} adsorption on a bare and Cl-terminated Si(100)-2$\times$1 surface, as well as for \ce{PH3} adsorption on H-terminated Si(100)-2$\times$1, which is widely used in current technologies for atomically precise doping of Si(100) with phosphorus. It was found that the Si(100)-2$\times$1-Cl surface has a higher reactivity towards phosphine than Si(100)-2$\times$1-H, and, therefore, unpatterned areas are less protected from undesirable incorporation of \ce{PH3} fragments. On the contrary, the resistance of the Si(100)-2$\times$1-Cl surface against the chlorine-containing molecules turned out to be very high. Several factors influencing reactivity are discussed. The results reveal that phosphorus and boron trichlorides are well-suited for doping a patterned Cl-resist by donors and acceptors, respectively.


title: Gap probability for products of random matrices in the critical regime

abstract: The singular values of a product of $M$ independent Ginibre matrices of size $N\times N$ form a determinantal point process. Near the soft edge, as both $M$ and $N$ go to infinity in such a way that $M/N\to \alpha$, $\alpha>0$, a scaling limit emerges. We consider a gap probability for the corresponding limiting determinantal process, namely, the probability that there are no particles in the interval $(a,+\infty)$. We derive a Tracy-Widom-like formula in terms of the unique solution of a certain matrix Riemann-Hilbert problem of size $2 \times 2$. The right-tail asymptotics for this solution is obtained by the Deift-Zhou non-linear steepest descent analysis.


title: Robust Blind Source Separation by Soft Decision-Directed Non-Unitary Joint Diagonalization

abstract: Approximate joint diagonalization of a set of matrices provides a powerful framework for numerous statistical signal processing applications. For non-unitary joint diagonalization (NUJD) based on the least-squares (LS) criterion, outliers, also referred to as anomaly or discordant observations, have a negative influence on the performance, since squaring the residuals magnifies the effects of them. To solve this problem, we propose a novel cost function that incorporates the soft decision-directed scheme into the least-squares algorithm and develops an efficient algorithm. The influence of the outliers is mitigated by applying decision-directed weights which are associated with the residual error at each iterative step. Specifically, the mixing matrix is estimated by a modified stationary point method, in which the updating direction is determined based on the linear approximation to the gradient function. Simulation results demonstrate that the proposed algorithm outperforms conventional non-unitary diagonalization algorithms in terms of both convergence performance and robustness to outliers.


title: Suppression of vortex shedding using a slit through the circular cylinder at low Reynolds number

abstract: The present article aims to study the suppression of vortex shedding using a passive flow control technique (slit through a circular cylinder) in the laminar regime (Re=100-500). The slit width ratio S/D (slit width/diameter) on the modified cylinder plays an essential role to control the vortex shedding. The additional flow through the slit leads to the suppression of the global instability and vortex shedding, whereas a large amount of flow through the slit drastically alters the behavior of vortex shedding. The nature of vortex shedding remains periodic for all S/D, and the root mean square (rms) value of the lift coefficient decreases (in turn, vortex shedding suppression) with S/D up to Re = 300. For the range Re > 300, the root mean square (rms) value of the lift coefficient decreases up to S/D < 0.15, and the flow exhibits periodic vortex shedding, while the root mean square (rms) increases beyond S/D > 0.15 due to irregular vortex shedding downstream of cylinder. The variation of the Re for the S/D=0.20 shows bifurcation points where the flow changes its behavior from symmetric to asymmetric solution at Re=232 and again becomes symmetric at Re=304. The unsteady flow analysis over the modified cylinder also indicates the suppression in the vortex shedding; however, the analysis provides the qualitative property of suppression. The reduced order modeling, i.e. Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD), is utilized to quantify suppression and investigate the dominant vortical structure for slit through the cylinder.


title: Arrhenius.jl: A Differentiable Combustion SimulationPackage

abstract: Combustion kinetic modeling is an integral part of combustion simulation, and extensive studies have been devoted to developing both high fidelity and computationally affordable models. Despite these efforts, modeling combustion kinetics is still challenging due to the demand for expert knowledge and optimization against experiments, as well as the lack of understanding of the associated uncertainties. Therefore, data-driven approaches that enable efficient discovery and calibration of kinetic models have received much attention in recent years, the core of which is the optimization based on big data. Differentiable programming is a promising approach for learning kinetic models from data by efficiently computing the gradient of objective functions to model parameters. However, it is often challenging to implement differentiable programming in practice. Therefore, it is still not available in widely utilized combustion simulation packages such as CHEMKIN and Cantera. Here, we present a differentiable combustion simulation package leveraging the eco-system in Julia, including DifferentialEquations.jl for solving differential equations, ForwardDiff.jl for auto-differentiation, and Flux.jl for incorporating neural network models into combustion simulations and optimizing neural network models using the state-of-the-art deep learning optimizers. We demonstrate the benefits of differentiable programming in efficient and accurate gradient computations, with applications in uncertainty quantification, kinetic model reduction, data assimilation, and model discovery.


title: Orthogonal and Non-Orthogonal Signal Representations Using New Transformation Matrices Having NPM Structure

abstract: In this paper, we introduce two types of real-valued sums known as Complex Conjugate Pair Sums (CCPSs) denoted as CCPS$^{(1)}$ and CCPS$^{(2)}$, and discuss a few of their properties. Using each type of CCPSs and their circular shifts, we construct two non-orthogonal Nested Periodic Matrices (NPMs). As NPMs are non-singular, this introduces two non-orthogonal transforms known as Complex Conjugate Periodic Transforms (CCPTs) denoted as CCPT$^{(1)}$ and CCPT$^{(2)}$. We propose another NPM, which uses both types of CCPSs such that its columns are mutually orthogonal, this transform is known as Orthogonal CCPT (OCCPT). After a brief study of a few OCCPT properties like periodicity, circular shift, etc., we present two different interpretations of it. Further, we propose a Decimation-In-Time (DIT) based fast computation algorithm for OCCPT (termed as FOCCPT), whenever the length of the signal is equal to $2^v,\ v{\in} \mathbb{N}$. The proposed sums and transforms are inspired by Ramanujan sums and Ramanujan Period Transform (RPT). Finally, we show that the period (both divisor and non-divisor) and frequency information of a signal can be estimated using the proposed transforms with a significant reduction in the computational complexity over Discrete Fourier Transform (DFT).


title: National-scale electricity peak load forecasting: Traditional, machine learning, or hybrid model?

abstract: As the volatility of electricity demand increases owing to climate change and electrification, the importance of accurate peak load forecasting is increasing. Traditional peak load forecasting has been conducted through time series-based models; however, recently, new models based on machine or deep learning are being introduced. This study performs a comparative analysis to determine the most accurate peak load-forecasting model for Korea, by comparing the performance of time series, machine learning, and hybrid models. Seasonal autoregressive integrated moving average with exogenous variables (SARIMAX) is used for the time series model. Artificial neural network (ANN), support vector regression (SVR), and long short-term memory (LSTM) are used for the machine learning models. SARIMAX-ANN, SARIMAX-SVR, and SARIMAX-LSTM are used for the hybrid models. The results indicate that the hybrid models exhibit significant improvement over the SARIMAX model. The LSTM-based models outperformed the others; the single and hybrid LSTM models did not exhibit a significant performance difference. In the case of Korea's highest peak load in 2019, the predictive power of the LSTM model proved to be greater than that of the SARIMAX-LSTM model. The LSTM, SARIMAX-SVR, and SARIMAX-LSTM models outperformed the current time series-based forecasting model used in Korea. Thus, Korea's peak load-forecasting performance can be improved by including machine learning or hybrid models.


title: Triple Coding Empowered FDMA-CDMA Mode High Security CAOS Camera

abstract: For the first time, the hybrid triple coding empowered Frequency Division Multiple Access (FDMA) Code Division Multiple Access (CDMA) mode of the CAOS (i.e., Coded Access Optical Sensor) camera is demonstrated. Compared to the independent FDMA and CDMA modes, the FDMA-CDMA mode has a novel high security space-time-frequency triple signal encoding design for robust, faster, linear irradiance extraction at a moderately High Dynamic Range (HDR). Specifically, this hybrid mode simultaneously combines the linear HDR strength of the FDMA mode Fast Fourier Transform (FFT) Digital Signal Processing (DSP)-based spectrum analysis with the high Signal to Noise Ratio (SNR) provided by the many simultaneous CAOS pixels photodetection of the CDMA mode. In particular, the demonstrated FDMA CDMA mode with P FDMA channels provides a P times faster camera operation versus the equivalent linear HDR Frequency Modulation (FM)CDMA mode. The active FDMA CDMA mode CAOS camera operation is also demonstrated using P equal to 3 LED light sources, each with its unique optical spectral content driven by its independent FDMA frequency. This illuminated target spectral signature matched active CAOS mode allows simultaneous capture of P images without the use of P time multiplexed slots operation tunable optical filter.


title: An Efficient Shock Advice Algorithm based on K-nearest Neighbors for Automated External Defibrillators

abstract: Shockable rhythms, namely ventricular fibrillation and ventricular tachycardia, are the main cause of sudden cardiac arrests, which can be detected quickly by the automated external defibrillator (AED) devices. In this paper, a simple but effective algorithm is proposed as the shock advice algorithm applied in AED. The proposed algorithm consists of K-nearest neighbor classifier and an optimal set of 36 features, which are extracted from original ECG and shockable, non-shockable signals using modified variational mode decomposition technique. Cross-validation procedure and sequential forward feature selection are carefully applied to select an optimal set from entire feature space. The performance results show that the MVMD is the key element for SCA detection performance, and the proposed algorithm is simpler while remaining relatively high detection performance compared to previous publications.


title: Impedance-based Capacity Estimation for Lithium-Ion Batteries Using Generative Adversarial Network

abstract: This paper proposes a fully unsupervised methodology for the reliable extraction of latent variables representing the characteristics of lithium-ion batteries (LIBs) from electrochemical impedance spectroscopy (EIS) data using information maximizing generative adversarial networks. Meaningful representations can be obtained from EIS data even when measured with direct current and without relaxation, which are difficult to express when using circuit models. The extracted latent variables were investigated as capacity degradation progressed and were used to estimate the discharge capacity of the batteries by employing Gaussian process regression. The proposed method was validated under various conditions of EIS data during charging and discharging. The results indicate that the proposed model provides more robust capacity estimations than the direct capacity estimations obtained from EIS. We demonstrate that the latent variables extracted from the EIS data measured with direct current and without relaxation reliably represent the degradation characteristics of LIBs.


title: An Ecological Robustness-Oriented Approach for Power System Network Expansion

abstract: Power systems are critical infrastructure for reliable and secure electric energy delivery. Incidents are increasing, as unexpected multiple hazards ranging from natural disasters to cyberattacks threaten the security and functionality of society. Inspired by resilient ecosystems, this paper presents a resilient network design approach with an ecological robustness (R ECO )-oriented optimization to improve power systems' ability to maintain a secure operating state throughout unknown hazards. The approach uses R ECO , a surprisal -based metric that captures key features of an ecosystem's resilient structure, as an objective to strategically design the electrical network. The approach enables solvability and practicality by introducing a stochastic-based candidate branch creation algorithm and a Taylor series expansion for relaxation of the R ECO formulation. Finally, studies are conducted on the R ECO -oriented approach using the IEEE 24 Bus RTS and the ACTIVSg200 systems. Results demonstrate improvement of the system's reliability under multiple hazards, network properties of robust structure and equally distributed power flows, and survivability against cascading failures. From the analysis, we observe that a more redundant network structure with equally distributed power flows benefits its resilience.


title: Application of artificial intelligence techniques for automated detection of myocardial infarction: A review

abstract: Myocardial infarction (MI) results in heart muscle injury due to receiving insufficient blood flow. MI is the most common cause of mortality in middle-aged and elderly individuals around the world. To diagnose MI, clinicians need to interpret electrocardiography (ECG) signals, which requires expertise and is subject to observer bias. Artificial intelligence-based methods can be utilized to screen for or diagnose MI automatically using ECG signals. In this work, we conducted a comprehensive assessment of artificial intelligence-based approaches for MI detection based on ECG as well as other biophysical signals, including machine learning (ML) and deep learning (DL) models. The performance of traditional ML methods relies on handcrafted features and manual selection of ECG signals, whereas DL models can automate these tasks. The review observed that deep convolutional neural networks (DCNNs) yielded excellent classification performance for MI diagnosis, which explains why they have become prevalent in recent years. To our knowledge, this is the first comprehensive survey of artificial intelligence techniques employed for MI diagnosis using ECG and other biophysical signals.


title: Raspberry PI for compact autonomous home farm control

abstract: This manuscript presented an autonomous home farm for predicting metrological characteristics that can not only automate the process of growing crops but also, due to a neural network, significantly increase the productivity of the farm. The developed farm monitors and manages the following indicators: illumination, soil PH, air temperature, ground temperature, air humidity, CO2 concentration, and soil moisture. The presented farm can also be considered as a device for testing various weather conditions to determine the optimal temperature characteristics for different crops. This farm as a result is completely autonomous grows tomatoes at home.


title: Intermittent Jamming against Telemetry and Telecommand of Satellite Systems and A Learning-driven Detection Strategy

abstract: Towards sixth-generation networks (6G), satellite communication systems, especially based on Low Earth Orbit (LEO) networks, become promising due to their unique and comprehensive capabilities. These advantages are accompanied by a variety of challenges such as security vulnerabilities, management of hybrid systems, and high mobility. In this paper, firstly, a security deficiency in the physical layer is addressed with a conceptual framework, considering the cyber-physical nature of the satellite systems, highlighting the potential attacks. Secondly, a learning-driven detection scheme is proposed, and the lightweight convolutional neural network (CNN) is designed. The performance of the designed CNN architecture is compared with a prevalent machine learning algorithm, support vector machine (SVM). The results show that deficiency attacks against the satellite systems can be detected by employing the proposed scheme.


title: Predictive models for wind speed using artificial intelligence and copula

abstract: Electricity generation from burning fossil fuels is one of the major contributors to global warming. Renewable energy sources are a viable alternative to produce electrical energy and to reduce the emission from the power industry. These energy sources are the building blocks of green energy, which all have different characteristics. Their availabilities are also diverse, depending on geographical locations and other parameters. Low implementation cost and distributed availability all over the world uplifts their popularity exponentially. Therefore, it has unlocked opportunities for consumers to produce electricity locally and use it on-site, which reduces dependency on centralized utility companies. The research considers two main objectives: the prediction of wind speed that simplifies wind farm planning and feasibility study. Secondly, the need to understand the dependency structure of the wind speeds of multiple distant locations. To address the first objective, twelve artificial intelligence algorithms were used for wind speed prediction from collected meteorological parameters. The model performances were compared to determine the wind speed prediction accuracy. The results show a deep learning approach, long short-term memory (LSTM) outperforms other models with the highest accuracy of 97.8%. For dependency, a multivariate cumulative distribution function, Copula, was used to find the joint distribution of two or more distant location wind speeds, followed by a case study. We found that the appropriate copula family and the parameters vary based on the distance in between. For the case study, Joe-Frank (BB8) copula shows an efficient joint distribution fit for a wind speed pair with a standard error of 0.0094. Finally, some insights about the uncertainty aspects of wind speed dependency were addressed.


title: A Self-Regulated and Reconfigurable CMOS Physically Unclonable Function Featuring Zero-Overhead Stabilization

abstract: This article presents a reconfigurable physically unclonable function (PUF) design fabricated using 65-nm CMOS technology. A subthreshold-inverter-based static PUF cell achieves 0.3% native bit error rate (BER) at 0.062-fJ per bit core energy efficiency. A flexible, native transistor-based voltage regulation scheme achieves low-overhead supply regulation with 6-mV/V line sensitivity, making the PUF resistant against voltage variations. Additionally, the PUF cell is designed to be reconfigurable with no area overhead, which enables stabilization without redundancy on chip. Thanks to the highly stable and self-regulated PUF cell and the zero-overhead stabilization scheme, a 0.00182% native BER is obtained after reconfiguration. The proposed design shows 0.12%/10 {\deg}C and 0.057%/0.1-V bit error across the military-grade temperature range from -55 {\deg}C to 125 {\deg}C and supply voltage variation from 0.7 to 1.4 V. The total energy per bit is 15.3 fJ. Furthermore, the unstable bits can be detected by sweeping the body bias instead of temperature during enrollment, thereby significantly reducing the testing costs. Last but not least, the prototype exhibits almost ideal uniqueness and randomness, with a mean inter-die Hamming distance (HD) of 0.4998 and a 1020x inter-/intra-die HD separation. It also passes both NIST 800-22 and 800-90B randomness tests.


title: AI Algorithm for Mode Classification of PCF SPR Sensor Design

abstract: Photonic Crystal Fiber design based on surface plasmon resonance phenomenon (PCF SPR) is optimized before it is fabricated for a particular application. An artificial intelligence algorithm is evaluated here to increase the ease of the simulation process for common users. COMSOL MultiPhysics is used. The algorithm suggests best among eight standard machine learning and one deep learning model to automatically select the desired mode, chosen visually by the experts otherwise. Total seven performance indices: namely Precision, Recall, Accuracy, F1-Score, Specificity, Matthew correlation coefficient, are utilized to make the optimal decision. Robustness towards variations in sensor geometry design is also considered as an optimal parameter. Several PCF-SPR based Photonic sensor designs are tested, and a large range optimal (based on phase matching) design is proposed. For this design algorithm has selected Support Vector Machine (SVM) as the best option with an accuracy of 96%, F1-Score is 95.83%, and MCC of 92.30%. The average sensitivity of the proposed sensor design with respect to change in refractive index (1.37-1.41) is 5500 nm/RIU. Resolution is 2.0498x10^(-5) RIU^(-1). The algorithm can be integrated into commercial software as an add-on or as a module in academic codes. The proposed novel step has saved approximately 75 minutes in the overall design process. The present work is equally applicable for mode selection of sensor other than PCF-SPR sensing geometries.


title: A new method for vehicle system safety design based on data mining with uncertainty modeling

abstract: In this research, a new data mining-based design approach has been developed for designing complex mechanical systems such as a crashworthy passenger car with uncertainty modeling. The method allows exploring the big crash simulation dataset to design the vehicle at multi-levels in a top-down manner (main energy absorbing system, components, and geometric features) and derive design rules based on the whole vehicle body safety requirements to make decisions towards the component and sub-component level design. Full vehicle and component simulation datasets are mined to build decision trees, where the interrelationship among parameters can be revealed and the design rules are derived to produce designs with good performance. This method has been extended by accounting for the uncertainty in the design variables. A new decision tree algorithm for uncertain data (DTUD) is developed to produce the desired designs and evaluate the design performance variations due to the uncertainty in design variables. The framework of this method is implemented by combining the design of experiments (DOE) and crash finite element analysis (FEA) and then demonstrated by designing a passenger car subject to front impact. The results show that the new methodology could achieve the design objectives efficiently and effectively. By applying the new method, the reliability of the final designs is also increased greatly. This approach has the potential to be applied as a general design methodology for a wide range of complex structures and mechanical systems.


title: Network structure and disease risk for an endemic infectious disease

abstract: The structure of contact networks affects the likelihood of disease spread at the population scale and the risk of infection at any given node. Though this has been well characterized for both theoretical and empirical networks for the spread of epidemics on completely susceptible networks, the long-term impact of network structure on risk of infection with an endemic pathogen has been less well characterized. Here, we analyze detailed records of the transportation of cattle between farms in Turkey to characterize the global and local attributes of the shipments network between 2007-2012, building an aggregated static directed - weighted network. We study the correlation between network properties and the likelihood of infection with, or exposure to, foot-and-mouth (FMD) disease over the same time period using recorded outbreaks. The shipments network shows properties of small-worldness and scale-freeness (intermediate degrees), with an exponential cut-off for high degrees. The shipments network illustrates strong modularity and lack of assortativity. The shipments network shows signs of spatial constraints, with few long-distance connections, and a strong similarity to other spatially constrained networks. We find that farms that were either infected or at high risk of infection with FMD (within one link from an infected farm) had higher values of centrality; farms that were never less than 2 links from an infected farm had disproportionately low centrality. However, the correlation of the rankings of farms shows that central farms (high eigenvector centrality) are not necessarily those with more connections to/from it (in/out degree). Several central farms serve as bridges of densely connected farms (high betweenness centrality). These results suggest that to detect FMD spread, surveillance efforts could be focused preferentially on farms with centralities greater than the mean.


title: Deep Ranking with Adaptive Margin Triplet Loss

abstract: We propose a simple modification from a fixed margin triplet loss to an adaptive margin triplet loss. While the original triplet loss is used widely in classification problems such as face recognition, face re-identification and fine-grained similarity, our proposed loss is well suited for rating datasets in which the ratings are continuous values. In contrast to original triplet loss where we have to sample data carefully, in out method, we can generate triplets using the whole dataset, and the optimization can still converge without frequently running into a model collapsing issue. The adaptive margins only need to be computed once before the training, which is much less expensive than generating triplets after every epoch as in the fixed margin case. Besides substantially improved training stability (the proposed model never collapsed in our experiments compared to a couple of times that the training collapsed on existing triplet loss), we achieved slightly better performance than the original triplet loss on various rating datasets and network architectures.


title: The degree of asymmetry of sequences

abstract: We explore the notion of degree of asymmetry for integer sequences and related combinatorial objects. The degree of asymmetry is a new combinatorial statistic that measures how far an object is from being symmetric. We define this notion for compositions, words, matchings, binary trees and permutations, we find generating functions enumerating these objects with respect to their degree of asymmetry, and we describe the limiting distribution of this statistic in each case.


title: Bosonic fractional Chern insulating state at integer fillings in multi-band system

abstract: The integer quantum Hall state occurs when the Landau levels are fully occupied by the fermions, while the fractional quantum Hall state usually emerges when the Landau level is partially filled by the strongly correlated fermions or bosons. Here, we report two fractional Chern insulating states of the hard-core bosons in a multi-band lattice model hosting topological flat bands with high Chern number. The previously proposed $\nu=1/3$ fractional Chern insulating state inherited from the high Chern number $C=2$ of the lowest topological flat band is revisited by the infinite density matrix renormalization group algorithm. In particular, we numerically identify a bosonic $1/2$-Laughlin-like fractional Chern insulating state at the integer fillings. We show two lower topological flat bands jointly generate an effective $C=1$ Chern band with half-filling. Furthermore, we find a strictly particle-hole-like symmetry between the $\nu$ and $3-\nu$ filling in our model. These findings extend our understanding of quantum Hall states and offer a new route to realize the novel fractional states in the system with multi-bands and high-Chern numbers.


title: Towards Machine Learning-Enabled Context Adaption for Reliable Aerial Mesh Routing

abstract: In this paper, we present Context-Adaptive PARRoT (CA-PARRoT) as an extension of our previous work Predictive Ad-hoc Routing fueled by Reinforcement learning and Trajectory knowledge (PARRoT). Short-term effects, as occurring in urban surroundings, have shown to have a negative impact on the Reinforcement Learning (RL)-based routing process. Therefore, we add a timer-based compensation mechanism to the update process and introduce a hybrid Machine Learning (ML) approach to classify Radio Environment Prototypes (REPs) with a dedicated ML component and enable the protocol for autonomous context adaption. The performance of the novel protocol is evaluated in comprehensive network simulations considering different REPs and is compared to well-known established routing protocols for Mobile Ad-hoc Networks (MANETs). The results show, that CA-PARRoT is capable to compensate the challenges confronted with in different REPs and to improve its Key Performance Indicators (KPIs) up to 23% compared to PARRoT, and outperform established routing protocols by up to 50 %.


title: Weakly nonlinear analysis of the viscoelastic instability in channel flow for finite and vanishing Reynolds numbers

abstract: The recently-discovered centre-mode instability of rectilinear viscoelastic shear flow (Garg et al. Phy. Rev. Lett. 121, 024502, 2018) has offered an explanation for the origin of elasto-inertial turbulence (EIT) which occurs at lower Weissenberg ($Wi$) numbers. In support of this, we show using weakly nonlinear analysis that the subcriticality found in Page et al. (Phys. Rev. Lett. 125, 154501, 2020) is generic across the neutral curve with the instability only becoming supercritical at low Reynolds ($Re$) numbers and high $Wi$. We demonstrate that the instability can be viewed as purely elastic in origin even for $Re=O(10^3)$, rather than `elasto-inertial', as the underlying shear does not energise the instability. It is also found that the introduction of a realistic maximum polymer extension length, $L_{max}$, in the FENE-P model moves the neutral curve closer to the inertialess $Re=0$ limit at a fixed ratio of solvent-to-solution viscosities, $\beta$. In the dilute limit ($\beta \rightarrow 1$) with $L_{max} =O(100)$, the linear instability can brought down to more physically-relevant $Wi\gtrsim 110$ at $\beta=0.98$, compared with the threshold $Wi=O(10^3)$ at $\beta=0.994$ reported recently by Khalid et al. (arXiv: 2103.06794) for an Oldroyd-B fluid. Again the instability is subcritical implying that inertialess rectilinear viscoelastic shear flow is nonlinearly unstable - i.e. unstable to finite amplitude disturbances - for even lower $Wi$.


title: Correlative image learning of chemo-mechanics in phase-transforming solids

abstract: Constitutive laws underlie most physical processes in nature. However, learning such equations in heterogeneous solids (e.g., due to phase separation) is challenging. One such relationship is between composition and eigenstrain, which governs the chemo-mechanical expansion in solids. In this work, we developed a generalizable, physically-constrained image-learning framework to algorithmically learn the chemo-mechanical constitutive law at the nanoscale from correlative four-dimensional scanning transmission electron microscopy and X-ray spectro-ptychography images. We demonstrated this approach on Li$_X$FePO$_4$, a technologically-relevant battery positive electrode material. We uncovered the functional form of composition-eigenstrain relation in this two-phase binary solid across the entire composition range (0 $\leq$ X $\leq$ 1), including inside the thermodynamically-unstable miscibility gap. The learned relation directly validates Vegard's law of linear response at the nanoscale. Our physics-constrained data-driven approach directly visualizes the residual strain field (by removing the compositional and coherency strain), which is otherwise impossible to quantify. Heterogeneities in the residual strain arise from misfit dislocations and were independently verified by X-ray diffraction line profile analysis. Our work provides the means to simultaneously quantify chemical expansion, coherency strain and dislocations in battery electrodes, which has implications on rate capabilities and lifetime. Broadly, this work also highlights the potential of integrating correlative microscopy and image learning for extracting material properties and physics.


title: Large-eddy simulations of marine boundary-layer clouds associated with cold air outbreaks during the ACTIVATE campaign-part 1: Case setup and sensitivities to large-scale forcings

abstract: Large-eddy simulation (LES) is able to capture key boundary-layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LES simulations. We derive these quantities from measurements and reanalysis obtained for two cold air outbreak (CAO) events during Phase I of the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) in February-March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from the reanalysis data ERA5 agree reasonably well with those derived from ACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time evolving heat fluxes, wind and advective tendencies profiles from ERA5 reanalysis data to drive the LES. We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 reanalysis data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LES simulations. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path simulated in LES agrees reasonably well with the ACTIVATE measurements.This study paves the path to investigate aerosol-cloud-meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.


title: Geometric insights into robust portfolio construction

abstract: We investigate and extend the results of Golts and Jones (2009) that an $\alpha$-weight angle resulting from unconstrained quadratic portfolio optimisations has an upper bound dependent on the condition number of the covariance matrix. This implies that better conditioned covariance matrices produce weights from unconstrained mean-variance optimisations that are better aligned with each assets expected return. We provide further clarity on the mathematical insights that relate the inequality between the $\alpha$-weight angle and the condition number and extend the result to include portfolio optimisations with gearing constraints. We provide an extended family of robust optimisations that include the gearing constraints, and discuss their interpretation.


title: Transfer Learning in Multi-Agent Reinforcement Learning with Double Q-Networks for Distributed Resource Sharing in V2X Communication

abstract: This paper addresses the problem of decentralized spectrum sharing in vehicle-to-everything (V2X) communication networks. The aim is to provide resource-efficient coexistence of vehicle-to-infrastructure(V2I) and vehicle-to-vehicle(V2V) links. A recent work on the topic proposes a multi-agent reinforcement learning (MARL) approach based on deep Q-learning, which leverages a fingerprint-based deep Q-network (DQN) architecture. This work considers an extension of this framework by combining Double Q-learning (via Double DQN) and transfer learning. The motivation behind is that Double Q-learning can alleviate the problem of overestimation of the action values present in conventional Q-learning, while transfer learning can leverage knowledge acquired by an expert model to accelerate learning in the MARL setting. The proposed algorithm is evaluated in a realistic V2X setting, with synthetic data generated based on a geometry-based propagation model that incorporates location-specific geographical descriptors of the simulated environment(outlines of buildings, foliage, and vehicles). The advantages of the proposed approach are demonstrated via numerical simulations.


title: No Regrets for Learning the Prior in Bandits

abstract: We propose ${\tt AdaTS}$, a Thompson sampling algorithm that adapts sequentially to bandit tasks that it interacts with. The key idea in ${\tt AdaTS}$ is to adapt to an unknown task prior distribution by maintaining a distribution over its parameters. When solving a bandit task, that uncertainty is marginalized out and properly accounted for. ${\tt AdaTS}$ is a fully-Bayesian algorithm that can be implemented efficiently in several classes of bandit problems. We derive upper bounds on its Bayes regret that quantify the loss due to not knowing the task prior, and show that it is small. Our theory is supported by experiments, where ${\tt AdaTS}$ outperforms prior algorithms and works well even in challenging real-world problems.


title: Generative Adversarial Learning via Kernel Density Discrimination

abstract: We introduce Kernel Density Discrimination GAN (KDD GAN), a novel method for generative adversarial learning. KDD GAN formulates the training as a likelihood ratio optimization problem where the data distributions are written explicitly via (local) Kernel Density Estimates (KDE). This is inspired by the recent progress in contrastive learning and its relation to KDE. We define the KDEs directly in feature space and forgo the requirement of invertibility of the kernel feature mappings. In our approach, features are no longer optimized for linear separability, as in the original GAN formulation, but for the more general discrimination of distributions in the feature space. We analyze the gradient of our loss with respect to the feature representation and show that it is better behaved than that of the original hinge loss. We perform experiments with the proposed KDE-based loss, used either as a training loss or a regularization term, on both CIFAR10 and scaled versions of ImageNet. We use BigGAN/SA-GAN as a backbone and baseline, since our focus is not to design the architecture of the networks. We show a boost in the quality of generated samples with respect to FID from 10% to 40% compared to the baseline. Code will be made available.


title: On the invariant E(G) for groups of odd order

abstract: Let $G$ be a multiplicatively written finite group. We denote by $\mathsf E(G)$ the smallest integer $t$ such that every sequence of $t$ elements in $G$ contains a product-one subsequence of length $|G|$. In 1961, Erd\H{o}s, Ginzburg and Ziv proved that $\mathsf E(G)\leq 2|G|-1$ for every finite solvable group $G$ and this result is well known as the Erd\H{o}s-Ginzburg-Ziv Theorem. In 2010, Gao and Li improved this result to $\mathsf E(G)\leq\frac{7|G|}{4}-1$ and they conjectured that $\mathsf E(G)\leq \frac{3|G|}{2}$ holds for any finite non-cyclic group. In this paper, we confirm the conjecture for all finite non-cyclic groups of odd order.


title: Many-body theory for quasiparticle states in superfluid fermionic systems

abstract: We present a formalism for the fermionic quasiparticle propagator in a superfluid fermionic system. Starting from a general many-body hamiltonian confined by the two-body instantaneous interaction, the equation of motion for the fermionic propagator is obtained in the Dyson form. Before making any approximation, the interaction kernel is found to be decomposed into the static and dynamical (time-dependent) contributions, while the latter translates to the energy-dependent and the former maps to the energy-independent terms in the energy domain. The three-fermion correlation function being the heart of the dynamical part of the kernel is factorized into the two-fermion and one-fermion ones. With the relaxed particle number constraint, the normal propagator is coupled to the abnormal one via both the static and dynamical kernels, that is formalized by introducing the generalized quasiparticle propagator of the Gor'kov type. The dynamical kernel in the factorized form is associated with the quasiparticle-vibration coupling (QVC) with the vibrations unifying both the normal and pairing phonons. The QVC vertices are related to the variations of the Hamiltonian of the Bogoliubov quasiparticles, which can be obtained by the finite amplitude method.


title: Multistage Mixed Precision Iterative Refinement

abstract: Low precision arithmetic, in particular half precision floating point arithmetic, is now available in commercial hardware. Using lower precision can offer significant savings in computation and communication costs with proportional savings in energy. Motivated by this, there has been a renewed interest in mixed precision iterative refinement for solving linear systems $Ax=b$, and new variants of GMRES-based iterative refinement have been developed. Each particular variant with a given combination of precisions leads to different condition number-based constraints for convergence of the backward and forward errors, and each has different performance costs. The constraints for convergence given in the literature are, as an artifact of the analyses, often overly strict in practice, and thus could lead a user to select a more expensive variant when a less expensive one would have sufficed. In this work, we develop a multistage mixed precision iterative refinement solver which aims to combine existing mixed precision approaches to balance performance and accuracy and improve usability. For an initial combination of precisions, the algorithm begins with the least expensive approach and convergence is monitored via inexpensive computations with quantities produced during the iteration. If slow convergence or divergence is detected using particular stopping criteria, the algorithm switches to use a more expensive, but more reliable variant. A novel aspect of our approach is that, unlike existing implementations, our algorithm first attempts to use ``stronger'' solvers for the solution update before resorting to increasing the precision(s). In some scenarios, this can avoid the need to refactorize the matrix in higher precision. We perform extensive numerical experiments on random dense problems and problems from real applications which confirm the benefits of the multistage approach.


title: Hamiltonian Complexity in the Thermodynamic Limit

abstract: Despite progress in quantum Hamiltonian complexity, little is known about the computational complexity of quantum physics at the thermodynamic limit. Even defining the problem is not straight forward. We study the complexity of estimating the ground energy of a fixed, translationally invariant Hamiltonian in the thermodynamic limit, to within a given precision; the number of bits $n$ for the precision is the sole input to the problem. The complexity of this problem captures how difficult it is for the physicist to measure or compute another digit in the approximation of a physical quantity in the thermodynamic limit. We show that this problem is contained in $\mbox{FEXP}^{\mbox{QMA-EXP}}$ and is hard for $\mbox{FEXP}^{\mbox{NEXP}}$. This means that the problem is doubly exponentially hard in the size of the input. As an ingredient in our construction, we study the problem of computing the ground energy of translationally invariant finite 1D chains. A single Hamiltonian term, which is a fixed parameter of the problem, is applied to every pair of particles in a finite chain. The length of the chain is the sole input to the problem and the task is to compute an approximation of the ground energy. No thresholds are provided as in the standard formulation of the local Hamiltonian problem. We show that this problem is contained in $\mbox{FP}^{\mbox{QMA-EXP}}$ and is hard for $\mbox{FP}^{\mbox{NEXP}}$. Our techniques employ a circular clock in which the ground energy is calibrated by the length of the cycle. This requires more precise expressions for the ground states of the resulting matrices than was required for previous QMA-completeness constructions and even exact analytical bounds for the infinite case which we derive using techniques from spectral graph theory. To our knowledge, this is the first use of the circuit-to-Hamiltonian construction which shows hardness for a function class.


title: Morse theory for loop-free categories

abstract: We extend discrete Morse-Bott theory to the setting of loop-free (or acyclic) categories. First of all, we state a homological version of Quillen's Theorem A in this context and introduce the notion of cellular categories. Second, we present a notion of vector field for loop-free categories. Third, we prove a homological collapsing theorem in the absence of critical objects in order to obtain the Morse inequalities. Examples are provided through the exposition. This answers partially a question by T. John: whether there is a Morse theory for loop-free (or acyclic) categories? [14].


title: Transverse Momentum Dependent Gluon Distribution within High Energy Factorization at Next-to-Leading Order

abstract: We discuss Transverse Momentum Dependent (TMD) gluon distributions within high energy factorization at next-to-leading order in the strong coupling within the framework of Lipatov's high energy effective action. We provide a detailed discussion of both rapidity divergences related to the TMD definition and its soft factor on the one hand, and rapidity divergences due to high energy factorization on the other hand, and discuss common features and differences between Collins-Soper (CS) and Balitsky-Fadin-Kuarev-Lipatov (BFKL) evolution. While we confirm earlier results which state that the unpolarized and linearly polarized gluon TMD agree in the BFKL limit at leading order, we find that both distributions differ, once next-to-leading order corrections are being included. Unlike previous results, our framework allows to recover the complete anomalous dimension associated with Collins-Soper-Sterman (CSS) evolution of the TMD distribution, including also single-logarithmic terms in the CSS evolution. As an additional result we provide a definition of $k_T$-factorization, i.e. matching of off-shell coefficients to collinear factorization at next-to-leading order within high energy factorization and the effective action framework. We furthermore establish a link between the QCD operator definition of the TMD gluon distribution and a previously derived off-shell TMD gluon-to-gluon splitting function, which is within the present framework obtained as the real 1-loop correction.


title: Non-homogeneous $(p_1,p_2)$-fractional Laplacian systems with lack of compactness

abstract: The present paper studies the existence of weak solutions for the following type of non-homogeneous system of equations \begin{equation*} (S) \left\{\begin{aligned} (-\Delta)^{s_1}_{p_1} u &=u|u|^{\alpha-1}|v|^{\beta+1}+f_1(x) \,\mbox{ in }\, \Omega, \\ (-\Delta)^{s_2}_{p_2} v &=|u|^{\alpha+1}v|v|^{\beta-1}+f_2(x) \,\mbox{ in }\, \Omega, \\ u=v &= 0 \,\mbox{ in }\, \mathbb{R}^N \setminus \Omega, \\ \end{aligned} \right. \end{equation*} where $\Omega \subset \mathbb{R}^N$ is smooth bounded domain, $s_1,s_2 \in (0,1)$, $1<p_1,p_2<\infty$, $N>\max\{p_1s_1,p_2s_2\}$, $\alpha>-1$ and $\beta>-1$. We employ the variational techniques where the associated energy functional is minimized over Nehari manifold set while imposing appropriate bound on dual norms of $f_1,f_2$.


title: Time-Division Multiplexing Light Field Display with Learned Coded Aperture

abstract: Conventional stereoscopic displays suffer from vergence-accommodation conflict and cause visual fatigue. Integral-imaging-based displays resolve the problem by directly projecting the sub-aperture views of a light field into the eyes using a microlens array or a similar structure. However, such displays have an inherent trade-off between angular and spatial resolutions. In this paper, we propose a novel coded time-division multiplexing technique that projects encoded sub-aperture views to the eyes of a viewer with correct cues for vergence-accommodation reflex. Given sparse light field sub-aperture views, our pipeline can provide a perception of high-resolution refocused images with minimal aliasing by jointly optimizing the sub-aperture views for display and the coded aperture pattern. This is achieved via deep learning in an end-to-end fashion by simulating light transport and image formation with Fourier optics. To our knowledge, this work is among the first that optimize the light field display pipeline with deep learning. We verify our idea with objective image quality metrics (PSNR, SSIM, and LPIPS) and perform an extensive study on various customizable design variables in our display pipeline. Experimental results show that light fields displayed using the proposed technique indeed have higher quality than that of baseline display designs.


title: ML-Quest: A Game for Introducing Machine Learning Concepts to K-12 Students

abstract: Today, Machine Learning (ML) is of a great importance to society due to the availability of huge data and high computational resources. This ultimately led to the introduction of ML concepts at multiple levels of education including K-12 students to promote computational thinking. However, teaching these concepts to K-12 through traditional methodologies such as video lectures and books is challenging. Many studies in the literature have reported that using interactive environments such as games to teach computational thinking and programming improves retention capacity and motivation among students. Therefore, introducing ML concepts using a game might enhance students' understanding of the subject and motivate them to learn further. However, we are not aware of any existing game which explicitly focuses on introducing ML concepts to students using game play. Hence, in this paper, we propose ML-Quest, a 3D video game to provide conceptual overview of three ML concepts: Supervised Learning, Gradient Descent and K-Nearest Neighbor (KNN) Classification. The crux of the game is to introduce the definition and working of these concepts, which we call conceptual overview, in a simulated scenario without overwhelming students with the intricacies of ML. The game has been predominantly evaluated for its usefulness and player experience using the Technology Acceptance Model (TAM) model with the help of 23 higher-secondary school students. The survey result shows that around 70% of the participants either agree or strongly agree that the ML-Quest is quite interactive and useful in introducing them to ML concepts.


title: Adaptive Machine Learning for Time-Varying Systems: Low Dimensional Latent Space Tuning

abstract: Machine learning (ML) tools such as encoder-decoder convolutional neural networks (CNN) can represent incredibly complex nonlinear functions which map between combinations of images and scalars. For example, CNNs can be used to map combinations of accelerator parameters and images which are 2D projections of the 6D phase space distributions of charged particle beams as they are transported between various particle accelerator locations. Despite their strengths, applying ML to time-varying systems, or systems with shifting distributions, is an open problem, especially for large systems for which collecting new data for re-training is impractical or interrupts operations. Particle accelerators are one example of large time-varying systems for which collecting detailed training data requires lengthy dedicated beam measurements which may no longer be available during regular operations. We present a recently developed method of adaptive ML for time-varying systems. Our approach is to map very high (N>100k) dimensional inputs (a combination of scalar parameters and images) into the low dimensional (N~2) latent space at the output of the encoder section of an encoder-decoder CNN. We then actively tune the low dimensional latent space-based representation of complex system dynamics by the addition of an adaptively tuned feedback vector directly before the decoder sections builds back up to our image-based high-dimensional phase space density representations. This method allows us to learn correlations within and to quickly tune the characteristics of incredibly high parameter systems and to track their evolution in real time based on feedback without massive new data sets for re-training.


title: Learning intermolecular forces at liquid-vapor interfaces

abstract: By adopting a perspective informed by contemporary liquid state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. We find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces, but typically fail for properties that depend on unbalanced long-ranged interactions which build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly-varying long-ranged interactions and training neural networks only on the short ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid where the exact long-ranged interaction is known, and in an ab initio water model where it is approximated.


title: Learning a Discriminant Latent Space with Neural Discriminant Analysis

abstract: Discriminative features play an important role in image and object classification and also in other fields of research such as semi-supervised learning, fine-grained classification, out of distribution detection. Inspired by Linear Discriminant Analysis (LDA), we propose an optimization called Neural Discriminant Analysis (NDA) for Deep Convolutional Neural Networks (DCNNs). NDA transforms deep features to become more discriminative and, therefore, improves the performances in various tasks. Our proposed optimization has two primary goals for inter- and intra-class variances. The first one is to minimize variances within each individual class. The second goal is to maximize pairwise distances between features coming from different classes. We evaluate our NDA optimization in different research fields: general supervised classification, fine-grained classification, semi-supervised learning, and out of distribution detection. We achieve performance improvements in all the fields compared to baseline methods that do not use NDA. Besides, using NDA, we also surpass the state of the art on the four tasks on various testing datasets.


title: Sissy That Walk: Transportation to Work by Sexual Orientation

abstract: We analyze differences in mode of transportation to work by sexual orientation, using the American Community Survey 2008-2019. Individuals in same-sex couples are significantly less likely to drive to work than men and women in different-sex couples. This gap is particularly stark among men: on average, almost 12 percentage point (or 13%) lower likelihood of driving to work for men in same-sex couples. Individuals in same-sex couples are also more likely to use public transport, walk, or bike to work: on average, men and women are 7 and 3 percentage points more likely, respectively, to take public transportation to work than those in different-sex couples. These differences persist after controlling for demographic characteristics, partner's characteristics, location, fertility, and marital status. Additional evidence from the General Social Survey 2008-2018 suggests that these disparities by sexual orientation may be due to lesbian, gay, and bisexual individuals caring more for the environment than straight individuals.


title: Attention-Guided Progressive Neural Texture Fusion for High Dynamic Range Image Restoration

abstract: High Dynamic Range (HDR) imaging via multi-exposure fusion is an important task for most modern imaging platforms. In spite of recent developments in both hardware and algorithm innovations, challenges remain over content association ambiguities caused by saturation, motion, and various artifacts introduced during multi-exposure fusion such as ghosting, noise, and blur. In this work, we propose an Attention-guided Progressive Neural Texture Fusion (APNT-Fusion) HDR restoration model which aims to address these issues within one framework. An efficient two-stream structure is proposed which separately focuses on texture feature transfer over saturated regions and multi-exposure tonal and texture feature fusion. A neural feature transfer mechanism is proposed which establishes spatial correspondence between different exposures based on multi-scale VGG features in the masked saturated HDR domain for discriminative contextual clues over the ambiguous image areas. A progressive texture blending module is designed to blend the encoded two-stream features in a multi-scale and progressive manner. In addition, we introduce several novel attention mechanisms, i.e., the motion attention module detects and suppresses the content discrepancies among the reference images; the saturation attention module facilitates differentiating the misalignment caused by saturation from those caused by motion; and the scale attention module ensures texture blending consistency between different coder/decoder scales. We carry out comprehensive qualitative and quantitative evaluations and ablation studies, which validate that these novel modules work coherently under the same framework and outperform state-of-the-art methods.


title: 'CADSketchNet' -- An Annotated Sketch dataset for 3D CAD Model Retrieval with Deep Neural Networks

abstract: Ongoing advancements in the fields of 3D modelling and digital archiving have led to an outburst in the amount of data stored digitally. Consequently, several retrieval systems have been developed depending on the type of data stored in these databases. However, unlike text data or images, performing a search for 3D models is non-trivial. Among 3D models, retrieving 3D Engineering/CAD models or mechanical components is even more challenging due to the presence of holes, volumetric features, presence of sharp edges etc., which make CAD a domain unto itself. The research work presented in this paper aims at developing a dataset suitable for building a retrieval system for 3D CAD models based on deep learning. 3D CAD models from the available CAD databases are collected, and a dataset of computer-generated sketch data, termed 'CADSketchNet', has been prepared. Additionally, hand-drawn sketches of the components are also added to CADSketchNet. Using the sketch images from this dataset, the paper also aims at evaluating the performance of various retrieval system or a search engine for 3D CAD models that accepts a sketch image as the input query. Many experimental models are constructed and tested on CADSketchNet. These experiments, along with the model architecture, choice of similarity metrics are reported along with the search results.


title: Modeling the Impact of Social Distancing and Targeted Vaccination on the Spread of COVID-19 through a Real City-Scale Contact Network

abstract: We use mobile device data to construct empirical interpersonal physical contact networks in the city of Portland, Oregon, both before and after social distancing measures were enacted during the COVID-19 pandemic. These networks reveal how social distancing measures and the public's reaction to the incipient pandemic affected the connectivity patterns within the city. We find that as the pandemic developed there was a substantial decrease in the number of individuals with many contacts. We further study the impact of these different network topologies on the spread of COVID-19 by simulating an SEIR epidemic model over these networks, and find that the reduced connectivity greatly suppressed the epidemic. We then investigate how the epidemic responds when part of the population is vaccinated, and we compare two vaccination distribution strategies, both with and without social distancing. Our main result is that the heavy-tailed degree distribution of the contact networks causes a targeted vaccination strategy that prioritizes high-contact individuals to reduce the number of cases far more effectively than a strategy that vaccinates individuals at random. Combining both targeted vaccination and social distancing leads to the greatest reduction in cases, and we also find that the marginal benefit of a targeted strategy as compared to a random strategy exceeds the marginal benefit of social distancing for reducing the number of cases. These results have important implications for ongoing vaccine distribution efforts worldwide.


title: Measuring kinetic parameters using quantum plasmonic sensing

abstract: The measurement of parameters that describe kinetic processes is important in the study of molecular interactions. It enables a deeper understanding of the physical mechanisms underlying how different biological entities interact with each other, such as viruses with cells, vaccines with antibodies, or new drugs with specific diseases. In this work, we study theoretically the use of quantum sensing techniques for measuring the kinetic parameters of molecular interactions. The sensor we consider is a plasmonic resonance sensor -- a label-free photonic sensor that is one of the most widely used in research and industry. The first type of interaction we study is the antigen BSA interacting with antibody IgG1, which provides a large sensor response. The second type is the enzyme carbonic anhydrase interacting with the tumor growth inhibitor benzenesulfonamide, which produces a small sensor response. For both types of interaction we consider the use of two-mode Fock states, squeezed vacuum states and squeezed displaced states. We find that these quantum states offer an enhancement in the measurement precision of kinetic parameters when compared to that obtained with classical light. The results may help in the design of more precise quantum-based sensors for studying kinetics in the life sciences.


title: Scoring from Pairwise Winning Indices

abstract: The pairwise winning indices, computed in the Stochastic Multicriteria Acceptability Analysis, give the probability with which an alternative is preferred to another taking into account all the instances of the assumed preference model compatible with the information provided by the Decision Maker in terms of pairwise preference comparisons of reference alternatives. In this paper we present a new scoring method assigning a value to each alternative summarizing the results of the pairwise winning indices. Several procedures assigning a score to each alternative on the basis of the pairwise winning indices have been provided in literature. However, while all of them compute this score just to rank the alternatives under consideration, our method, expressing the score in terms of an additive value function, permits to disaggregate the overall evaluation of each alternative in the sum of contributions of considered criteria. This will permit not only to rank the alternatives but also to explain the reasons for which an alternative obtains its evaluation and, consequently, fills a certain ranking position. To prove the efficiency of the method in representing the preferences of the Decision Maker, we performed an extensive set of simulations varying the number of alternatives and criteria. The results of the simulations, analyzed from a statistical point of view, show the goodness of our procedure. The applicability of the method to decision making problems is explained by means of a case study related to the evaluation of financial funds.


title: Bag-of-Tasks Scheduling on Related Machines

abstract: We consider online scheduling to minimize weighted completion time on related machines, where each job consists of several tasks that can be concurrently executed. A job gets completed when all its component tasks finish. We obtain an $O(K^3 \log^2 K)$-competitive algorithm in the non-clairvoyant setting, where $K$ denotes the number of distinct machine speeds. The analysis is based on dual-fitting on a precedence-constrained LP relaxation that may be of independent interest.


title: What classifiers know what they don't?

abstract: Being uncertain when facing the unknown is key to intelligent decision making. However, machine learning algorithms lack reliable estimates about their predictive uncertainty. This leads to wrong and overly-confident decisions when encountering classes unseen during training. Despite the importance of equipping classifiers with uncertainty estimates ready for the real world, prior work has focused on small datasets and little or no class discrepancy between training and testing data. To close this gap, we introduce UIMNET: a realistic, ImageNet-scale test-bed to evaluate predictive uncertainty estimates for deep image classifiers. Our benchmark provides implementations of eight state-of-the-art algorithms, six uncertainty measures, four in-domain metrics, three out-domain metrics, and a fully automated pipeline to train, calibrate, ensemble, select, and evaluate models. Our test-bed is open-source and all of our results are reproducible from a fixed commit in our repository. Adding new datasets, algorithms, measures, or metrics is a matter of a few lines of code-in so hoping that UIMNET becomes a stepping stone towards realistic, rigorous, and reproducible research in uncertainty estimation. Our results show that ensembles of ERM classifiers as well as single MIMO classifiers are the two best alternatives currently available to measure uncertainty about both in-domain and out-domain classes.


title: Tempering stochastic density functional theory

abstract: We introduce a tempering approach with stochastic density functional theory (sDFT), labeled t-sDFT, which reduces the statistical errors in the estimates of observable expectation values. This is achieved by rewriting the electronic density as a sum of a "warm" component complemented by "colder" correction(s). Since the "warm" component is larger in magnitude but faster to evaluate, we use many more stochastic orbitals for its evaluation than for the smaller-sized colder correction(s). This results in a significant reduction of the statistical fluctuations and the bias compared to sDFT for the same computational effort. We the method's performance on large hydrogen-passivated silicon nanocrystals (NCs), finding a reduction in the systematic error in the energy by more than an order of magnitude, while the systematic errors in the forces are also quenched. Similarly, the statistical fluctuations are reduced by factors of around 4-5 for the total energy and around 1.5-2 for the forces on the atoms. Since the embedding in t-sDFT is fully stochastic, it is possible to combine t-sDFT with other variants of sDFT such as energy-window sDFT and embedded-fragmented sDFT.


title: Towards Unsupervised Domain Generalization

abstract: Domain generalization (DG) aims to help models trained on a set of source domains generalize better on unseen target domains. The performances of current DG methods largely rely on sufficient labeled data, which are usually costly or unavailable, however. Since unlabeled data are far more accessible, we seek to explore how unsupervised learning can help deep models generalize across domains. Specifically, we study a novel generalization problem called unsupervised domain generalization (UDG), which aims to learn generalizable models with unlabeled data and analyze the effects of pre-training on DG. In UDG, models are pretrained with unlabeled data from various source domains before being trained on labeled source data and eventually tested on unseen target domains. Then we propose a method named Domain-Aware Representation LearnING (DARLING) to cope with the significant and misleading heterogeneity within unlabeled pretraining data and severe distribution shifts between source and target data. Surprisingly we observe that DARLING can not only counterbalance the scarcity of labeled data but also further strengthen the generalization ability of models when the labeled data are insufficient. As a pretraining approach, DARLING shows superior or comparable performance compared with ImageNet pretraining protocol even when the available data are unlabeled and of a vastly smaller amount compared to ImageNet, which may shed light on improving generalization with large-scale unlabeled data.


title: Lattice associated to a Shi variety

abstract: Let $W$ be a irreducible Weyl group and $W_a$ its affine Weyl group. In a previous article the author defined an affine variety $\widehat{X}_{W_a}$, called the Shi variety of $W_a$, whose integral points are in bijection with $W_a$. The set of irreducible components of $\widehat{X}_{W_a}$, denoted $H^0(\widehat{X}_{W_a})$, is of some interest and we show in this article that $H^0(\widehat{X}_{W_a})$ has a structure of semidistributive lattice.


title: Direct mapping of electronic orbitals in graphene using electron energy-loss spectroscopy

abstract: The spatial distributions of anti-bonding $\pi^\ast$ and $\sigma^\ast$ states in epitaxial graphene multilayers are mapped using electron energy-loss spectroscopy in a scanning transmission electron microscope. Inelastic channeling simulations validate the interpretation of the spatially-resolved signals in terms of electronic orbitals, and demonstrate the crucial effect of the material thickness on the experimental capability to resolve the distribution of unoccupied states. This work illustrates the current potential of core-level electron energy-loss spectroscopy towards the direct visualization of electronic orbitals in a wide range of materials, of huge interest to better understand chemical bonding among many other properties at interfaces and defects in solids.


title: Bottomonium production in heavy-ion collisions using quantum trajectories: Differential observables and momentum anisotropy

abstract: We report predictions for the suppression and elliptic flow of the $\Upsilon(1S)$, $\Upsilon(2S)$, and $\Upsilon(3S)$ as a function of centrality and transverse momentum in ultra-relativistic heavy-ion collisions. We obtain our predictions by numerically solving a Lindblad equation for the evolution of the heavy-quarkonium reduced density matrix derived using potential nonrelativistic QCD and the formalism of open quantum systems. To numerically solve the Lindblad equation, we make use of a stochastic unraveling called the quantum trajectories algorithm. This unraveling allows us to solve the Lindblad evolution equation efficiently on large lattices with no angular momentum cutoff. The resulting evolution describes the full 3D quantum and non-abelian evolution of the reduced density matrix for bottomonium states. We expand upon our previous work by treating differential observables and elliptic flow; this is made possible by a newly implemented Monte-Carlo sampling of physical trajectories. Our final results are compared to experimental data collected in $\sqrt{s_{NN}} = 5.02$ TeV Pb-Pb collisions by the ALICE, ATLAS, and CMS collaborations.


title: Impact of heat waves and cold spells on cause-specific mortality in the city of Sao Paulo, Brazil

abstract: The impact of heat waves and cold spells on mortality has become a major public health problem worldwide, especially among older adults living in low- to middle-income countries. This study aimed to investigate the effects of heat waves and cold spells under different definitions on cause-specific mortality among people aged 65 years and over in Sao Paulo from 2006 to 2015. A quasi-Poisson generalized linear model with a distributed lag model was used to investigate the association between cause-specific mortality and extreme air temperature events. To evaluate the effects of the intensity under different durations, we considered 12 heat wave and nine cold spell definitions. Our results showed an increase in cause-specific deaths related to heat waves and cold spells under several definitions. The highest risk of death related to heat waves was identified mostly at higher temperature thresholds with longer events. We verified that men were more vulnerable to die from an ischemic stroke on heat waves and cold spells days than women, while women presented a higher risk of dying from ischemic heart diseases during cold spells and tended to have a higher risk of chronic obstructive pulmonary disease than men. Identification of heat wave- and cold spell-related mortality is important for the development and promotion of public health measures.


title: T product Tensors Part II: Tail Bounds for Sums of Random T product Tensors

abstract: This paper is the Part II of a serious work about T product tensors focusing at establishing new probability bounds for sums of random, independent, T product tensors. These probability bounds characterize large deviation behavior of the extreme eigenvalue of the sums of random T product tensors. We apply Lapalace transform method and Lieb concavity theorem for T product tensors obtained from our Part I paper, and apply these tools to generalize the classical bounds associated with the names Chernoff, and Bernstein from the scalar to the T product tensor setting. Tail bounds for the norm of a sum of random rectangular T product tensors are also derived from corollaries of random Hermitian T product tensors cases. The proof mechanism is also applied to T product tensor valued martingales and T product tensor based Azuma, Hoeffding and McDiarmid inequalities are derived.


title: On string functions and double-sum formulas

abstract: String functions are important building blocks of characters of integrable highest modules over affine Kac--Moody algebras. Kac and Peterson computed string functions for affine Lie algebras of type $A_{1}^{(1)}$ in terms of Dedekind eta functions. We produce new relations between string functions by writing them as double-sums and then using certain symmetry relations. We evaluate the series using special double-sum formulas that express Hecke-type double-sums in terms of Appell--Lerch functions and theta functions, where we point out that Appell--Lerch functions are the building blocks of Ramanujan's classical mock theta functions.


title: Pessimistic Model-based Offline Reinforcement Learning under Partial Coverage

abstract: We study model-based offline Reinforcement Learning with general function approximation without a full coverage assumption on the offline data distribution. We present an algorithm named Constrained Pessimistic Policy Optimization (CPPO)which leverages a general function class and uses a constraint over the model class to encode pessimism. Under the assumption that the ground truth model belongs to our function class (i.e., realizability in the function class), CPPO has a PAC guarantee with offline data only providing partial coverage, i.e., it can learn a policy that competes against any policy that is covered by the offline data. We then demonstrate that this algorithmic framework can be applied to many specialized Markov Decision Processes where additional structural assumptions can further refine the concept of partial coverage. Two notable examples are: (1) low-rank MDP with representation learning where the partial coverage condition is defined using a relative condition number measured by the unknown ground truth feature representation; (2) factored MDP where the partial coverage condition is defined using density ratio based concentrability coefficients associated with individual factors.


title: Microlensing of the broad emission lines in 27 gravitationally lensed quasars. Broad line region structure and kinematics

abstract: We aim to study the structure and kinematics of the broad line region (BLR) of a sample of 27 gravitationally lensed quasars with up to five different epochs of observation. This sample is composed of ~100 spectra from the literature plus 22 unpublished spectra of 11 systems. We measure the magnitude differences in the broad emission line (BEL) wings and statistically model the distribution of microlensing magnifications to determine a maximum likelihood estimate for the sizes of the C IV, C III], and Mg II emitting regions. The BELs in lensed quasars are expected to be magnified differently owing to the different sizes of the regions from which they originate. Focusing on the most common BELs in our spectra (C IV, C III], and Mg II), we find that the low-ionization line Mg II is only weakly affected by microlensing. In contrast, the high-ionization line C IV shows strong microlensing in some cases, indicating that its emission region is more compact. Thus, the BEL profiles are deformed differently depending on the geometry and kinematics of the corresponding emitting region. We detect microlensing in either the blue or the red wing (or in both wings with different amplitudes) of C IV in more than 50% of the systems and find outstanding asymmetries in the wings of QSO 0957+561, SDSS J1004+4112, SDSS J1206+4332, and SDSS J1339+1310. This observation indicates that the BLR is, in general, not spherically symmetric and supports the existence of two regions in the BLR, one insensitive to microlensing and another that only shows up when it is magnified by microlensing.


title: Dominating ideals and closed neighborhood ideals of graphs

abstract: We study the closed neighborhood ideals and the dominating ideals of graphs, in particular, of trees and cycles. We prove that the closed neighborhood ideals and the dominating ideals of trees are normally torsion-free. The closed neighborhood ideals and the dominating ideals of cycles fail to be normally torsion-free. However, we prove that the closed neighborhood ideals of cycles admit the (strong) persistence property and the dominating ideals of cycles are nearly normally torsion-free.


title: Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions

abstract: Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not only probed in astrophysical observations, but also in terrestrial heavy-ion collision experiments. In this work, we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars and from heavy-ion collisions of gold nuclei at relativistic energies with microscopic nuclear theory calculations to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent NICER observations. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment, and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

categories: nucl-th astro-ph.HE astro-ph.SR gr-qc nucl-ex


title: Spiral Density Enhancements in Be Binary Systems

abstract: We use a smoothed particle hydrodynamics (SPH) code to examine the effects of a binary companion on a Be star disk for a range of disk viscosities and misalignment angles, i.e. the angle between the orbital plane and the primary's spin axis. The density structures in the disk due to the tidal interaction with the binary companion are investigated. Expanding on our previous work, the shape and density structure of density enhancements due to the binary companion are analyzed and the changes in observed interferometric features due to these orbiting enhancements are also predicted. We find that larger misalignment angles and viscosity values result in more tightly wound spiral arms with densities that fall-off more slowly with radial distance from the central star. We show that the orbital phase has very little effect on the structure of the spiral density enhancements. We demonstrate that these spiral features can be detected with an interferometer in H$\alpha$ and K-band emission. We also show that the spiral features affect the axis ratios determined by interferometry depending on the orientation of these features and the observer. For example, our simulations show that the axis ratios can vary by 20% for our co-planar binary disk system depending on the location of the disk density enhancements.


title: Timbre Classification of Musical Instruments with a Deep Learning Multi-Head Attention-Based Model

abstract: The aim of this work is to define a model based on deep learning that is able to identify different instrument timbres with as few parameters as possible. For this purpose, we have worked with classical orchestral instruments played with different dynamics, which are part of a few instrument families and which play notes in the same pitch range. It has been possible to assess the ability to classify instruments by timbre even if the instruments are playing the same note with the same intensity. The network employed uses a multi-head attention mechanism, with 8 heads and a dense network at the output taking as input the log-mel magnitude spectrograms of the sound samples. This network allows the identification of 20 instrument classes of the classical orchestra, achieving an overall F$_1$ value of 0.62. An analysis of the weights of the attention layer has been performed and the confusion matrix of the model is presented, allowing us to assess the ability of the proposed architecture to distinguish timbre and to establish the aspects on which future work should focus.


title: Maintaining $\mathsf{CMSO}_2$ properties on dynamic structures with bounded feedback vertex number

abstract: Let $\varphi$ be a sentence of $\mathsf{CMSO}_2$ (monadic second-order logic with quantification over edge subsets and counting modular predicates) over the signature of graphs. We present a dynamic data structure that for a given graph $G$ that is updated by edge insertions and edge deletions, maintains whether $\varphi$ is satisfied in $G$. The data structure is required to correctly report the outcome only when the feedback vertex number of $G$ does not exceed a fixed constant $k$, otherwise it reports that the feedback vertex number is too large. With this assumption, we guarantee amortized update time ${\cal O}_{\varphi,k}(\log n)$. By combining this result with a classic theorem of Erd\H{o}s and P\'osa, we give a fully dynamic data structure that maintains whether a graph contains a packing of $k$ vertex-disjoint cycles with amortized update time ${\cal O}_{k}(\log n)$. Our data structure also works in a larger generality of relational structures over binary signatures.


title: Generalizations of the Yao-Yao partition theorem and the central transversal theorem

abstract: We generalize the Yao-Yao partition theorem by showing that for any smooth measure in $R^d$ there exist equipartitions using $(t+1)2^{d-1}$ convex regions such that every hyperplane misses the interior of at least $t$ regions. In addition, we present tight bounds on the smallest number of hyperplanes whose union contains the boundary of an equipartition of a measure into $n$ regions. We also present a simple proof of a Borsuk-Ulam type theorem for Stiefel manifolds that allows us to generalize the central transversal theorem and prove results bridging the Yao--Yao partition theorem and the central transversal theorem.


title: Thermal variational quantum simulation on a superconducting quantum processor

abstract: Solving finite-temperature properties of quantum many-body systems is generally challenging to classical computers due to their high computational complexities. In this article, we present experiments to demonstrate a hybrid quantum-classical simulation of thermal quantum states. By combining a classical probabilistic model and a 5-qubit programmable superconducting quantum processor, we prepare Gibbs states and excited states of Heisenberg XY and XXZ models with high fidelity and compute thermal properties including the variational free energy, energy, and entropy with a small statistical error. Our approach combines the advantage of classical probabilistic models for sampling and quantum co-processors for unitary transformations. We show that the approach is scalable in the number of qubits, and has a self-verifiable feature, revealing its potentials in solving large-scale quantum statistical mechanics problems on near-term intermediate-scale quantum computers.


title: Exploiting Image Translations via Ensemble Self-Supervised Learning for Unsupervised Domain Adaptation

abstract: We introduce an unsupervised domain adaption (UDA) strategy that combines multiple image translations, ensemble learning and self-supervised learning in one coherent approach. We focus on one of the standard tasks of UDA in which a semantic segmentation model is trained on labeled synthetic data together with unlabeled real-world data, aiming to perform well on the latter. To exploit the advantage of using multiple image translations, we propose an ensemble learning approach, where three classifiers calculate their prediction by taking as input features of different image translations, making each classifier learn independently, with the purpose of combining their outputs by sparse Multinomial Logistic Regression. This regression layer known as meta-learner helps to reduce the bias during pseudo label generation when performing self-supervised learning and improves the generalizability of the model by taking into consideration the contribution of each classifier. We evaluate our method on the standard UDA benchmarks, i.e. adapting GTA V and Synthia to Cityscapes, and achieve state-of-the-art results in the mean intersection over union metric. Extensive ablation experiments are reported to highlight the advantageous properties of our proposed UDA strategy.


title: An FPT algorithm for the embeddability of graphs into two-dimensional simplicial complexes

abstract: We consider the embeddability problem of a graph G into a two-dimensional simplicial complex C: Given G and C, decide whether G admits a topological embedding into C. The problem is NP-hard, even in the restricted case where C is homeomorphic to a surface. It is known that the problem admits an algorithm with running time f(c).n^{O(c)}, where n is the size of the graph G and c is the size of the two-dimensional complex C. In other words, that algorithm is polynomial when C is fixed, but the degree of the polynomial depends on C. We prove that the problem is fixed-parameter tractable in the size of the two-dimensional complex, by providing a deterministic f(c).n^3-time algorithm. We also provide a randomized algorithm with expected running time 2^{c^{O(1)}}.n^{O(1)}. Our approach is to reduce to the case where G has bounded branchwidth via an irrelevant vertex method, and to apply dynamic programming. We do not rely on any component of the existing linear-time algorithms for embedding graphs on a fixed surface; the only elaborated tool that we use is an algorithm to compute grid minors.


title: Boson dark matter halos with a dominant noncondensed component

abstract: We consider galaxy halos formed by dark matter bosons with mass in the range of about a few tens or hundreds eV. A major part of the particles is in a noncondensed state and described under the Thomas-Fermi approach. Derived equations are solved numerically to find the halo density profile. The noncondensed state is supported in the entire halo except compact gravitationally bounded Bose-Einstein condensates. Although the size of these compact objects, also known as Bose stars, depends on interactions between the particles, its upper limit is only about 100 astronomical units. The Bose stars collect the condensed bosons providing a density cusp avoidance in the halo as well as a natural mechanism to prevent overproduction of small halos. Clusters of the Bose stars can also contribute to the halo density profile. The model is analyzed by confronting its predictions with observations of galaxy rotation curves. We employ 22 low surface brightness galaxies and obtain that the model is consistent with the observational data when the particle mass is in the range above about 50 eV and the best fit corresponds to the mass $m=86$ eV. This mass is appropriate for relic dark matter bosons, which decouple just after QCD phase transition.


title: GENIUS-MAWII: For Robust Mendelian Randomization with Many Weak Invalid Instruments

abstract: Mendelian randomization (MR) has become a popular approach to study causal effects by using genetic variants as instrumental variables. We propose a new MR method, GENIUS-MAWII, which simultaneously addresses the two salient phenomena that adversely affect MR analyses: many weak instruments and widespread horizontal pleiotropy. Similar to MR GENIUS (Tchetgen Tchetgen et al., 2021), we achieve identification of the treatment effect by leveraging heteroscedasticity of the exposure. We then derive the class of influence functions of the treatment effect, based on which, we construct a continuous updating estimator and establish its consistency and asymptotic normality under a many weak invalid instruments asymptotic regime by developing novel semiparametric theory. We also provide a measure of weak identification, an overidentification test, and a graphical diagnostic tool. We demonstrate in simulations that GENIUS-MAWII has clear advantages in the presence of directional or correlated horizontal pleiotropy compared to other methods. We apply our method to study the effect of body mass index on systolic blood pressure using UK Biobank.