title: Regolith behavior under asteroid-level gravity conditions: Low-velocity impacts into mm- and cm-sized grain targets

abstract: In situ observations of small asteroids show that surfaces covered by boulders and coarse terrain are frequent on such bodies. Regolith grain sizes have distributions on approximately mm and cm scales, and the behavior of such large grains in the very low-gravity environments of small body surfaces dictates their morphology and evolution. In order to support the understanding of natural processes (e.g., the recapturing of impact ejecta) or spacecraft-induced interactions (e.g., the fate of a small lander), we aim to experimentally investigate the response of coarse-grained target surfaces to very-low-speed impacts (below 2 m/s). We present the outcome of 86 low-speed impacts of a cm-sized spherical projectile into a bed of simulated regolith, composed of irregular mm- and cm-sized grains. These impacts were performed under vacuum and microgravity conditions. Our results include measurements for the projectile coefficient of restitution and penetration depth, as well as ejecta production, speed, and mass estimation. We find that impact outcomes include the frequent occurrence of projectile bouncing and tangential rolling on the target surface upon impact. Ejecta is produced for impact speeds higher than about 12 cm/s, and ejecta speeds scale with the projectile to target the grain size ratio and the impact speed. Ejected mass estimations indicate that ejecta is increasingly difficult to produce for increasing grain sizes. Coefficients of restitution of rebounding projectiles do not display a dependency on the target grain size, unlike their maximum penetration depth, which can be scaled with the projectile to target grain size ratio. Finally, we compare our experimental measurements to spacecraft data and numerical work on Hayabusa 2's MASCOT landing on the surface of the asteroid Ryugu.


title: The Quantum Supremacy Tsirelson Inequality

abstract: A leading proposal for verifying near-term quantum supremacy experiments on noisy random quantum circuits is linear cross-entropy benchmarking. For a quantum circuit $C$ on $n$ qubits and a sample $z \in \{0,1\}^n$, the benchmark involves computing $|\langle z|C|0^n \rangle|^2$, i.e. the probability of measuring $z$ from the output distribution of $C$ on the all zeros input. Under a strong conjecture about the classical hardness of estimating output probabilities of quantum circuits, no polynomial-time classical algorithm given $C$ can output a string $z$ such that $|\langle z|C|0^n\rangle|^2$ is substantially larger than $\frac{1}{2^n}$ (Aaronson and Gunn, 2019). On the other hand, for a random quantum circuit $C$, sampling $z$ from the output distribution of $C$ achieves $|\langle z|C|0^n\rangle|^2 \approx \frac{2}{2^n}$ on average (Arute et al., 2019). In analogy with the Tsirelson inequality from quantum nonlocal correlations, we ask: can a polynomial-time quantum algorithm do substantially better than $\frac{2}{2^n}$? We study this question in the query (or black box) model, where the quantum algorithm is given oracle access to $C$. We show that, for any $\varepsilon \ge \frac{1}{\mathrm{poly}(n)}$, outputting a sample $z$ such that $|\langle z|C|0^n\rangle|^2 \ge \frac{2 + \varepsilon}{2^n}$ on average requires at least $\Omega\left(\frac{2^{n/4}}{\mathrm{poly}(n)}\right)$ queries to $C$, but not more than $O\left(2^{n/3}\right)$ queries to $C$, if $C$ is either a Haar-random $n$-qubit unitary, or a canonical state preparation oracle for a Haar-random $n$-qubit state. We also show that when $C$ samples from the Fourier distribution of a random Boolean function, the naive algorithm that samples from $C$ is the optimal 1-query algorithm for maximizing $|\langle z|C|0^n\rangle|^2$ on average.


title: Electron-phonon drag enhancement of transport properties from fully coupled \textit{ab initio} Boltzmann formalism

abstract: We present a combined treatment of the non-equilibrium dynamics and transport of electrons and phonons by carrying out \textit{ab initio} calculations of the fully coupled electron and phonon Boltzmann transport equations. We find that the presence of mutual drag between the two carriers causes the thermopower to be enhanced and dominated by the transport of phonons, rather than electrons as in the traditional semiconductor picture. Drag also strongly boosts the intrinsic electron mobility, thermal conductivity and the Lorenz number. Impurity scattering is seen to suppress the drag-enhancement of the thermal and electrical conductivities, while having weak effects on the enhancement of the Lorenz number and thermopower. We demonstrate these effects in \textit{n}-doped 3C-SiC at room temperature, and explain their origins. This work establishes the roles of microscopic scattering mechanisms in the emergence of strong drag effects in the transport of the interacting electron-phonon gas.


title: Molecular Gas Outflow in the Starburst Galaxy NGC 1482

abstract: Galactic winds are essential to regulation of star formation in galaxies. To study the distribution and dynamics of molecular gas in a wind, we imaged the nearby starburst galaxy NGC 1482 in CO ($J=1\rightarrow0$) at a resolution of 1'' ($\approx100$ pc) using the Atacama Large Millimeter/submillimeter Array. Molecular gas is detected in a nearly edge-on disk with a radius of 3 kpc and a biconical outflow emerging from the central 1 kpc starburst and extending to at least 1.5 kpc perpendicular to the disk. In the outflow, CO gas is distributed approximately as a cylindrically symmetrical envelope surrounding the warm and hot ionized gas traced by H$\alpha$ and soft X-rays. The velocity, mass outflow rate, and kinetic energy of the molecular outflow are $v_\mathrm{w}\sim100~\mathrm{km~s^{-1}}$, $\dot{M}_\mathrm{w}\sim7~M_\odot~\mathrm{yr}^{-1}$, and $E_\mathrm{w}\sim7\times10^{54}~\mathrm{erg}$, respectively. $\dot{M}_\mathrm{w}$ is comparable to the star formation rate ($\dot{M}_\mathrm{w}/\mathrm{SFR}\sim2$) and $E_\mathrm{w}$ is $\sim1\%$ of the total energy released by stellar feedback in the past $1\times10^7~\mathrm{yr}$, which is the dynamical timescale of the outflow. The results indicate that the wind is starburst driven.


title: Global Phase Portrait and Large Degree Asymptotics for the Kissing Polynomials

abstract: We study a family of monic orthogonal polynomials which are orthogonal with respect to the varying, complex valued weight function, $\exp(nsz)$, over the interval $[-1,1]$, where $s\in\mathbb{C}$ is arbitrary. This family of polynomials originally appeared in the literature when the parameter was purely imaginary, that is $s\in i \mathbb{R}$, due to its connection with complex Gaussian quadrature rules for highly oscillatory integrals. The asymptotics for these polynomials as $n\to\infty$ have been recently studied for $s\in i\mathbb{R}$, and our main goal is to extend these results to all $s$ in the complex plane. We first use the technique of continuation in parameter space, developed in the context of the theory of integrable systems, to extend previous results on the so-called modified external field from the imaginary axis to the complex plane minus a set of critical curves, called breaking curves. We then apply the powerful method of nonlinear steepest descent for oscillatory Riemann-Hilbert problems developed by Deift and Zhou in the 1990s to obtain asymptotics of the recurrence coefficients of these polynomials when the parameter $s$ is away from the breaking curves. We then provide the analysis of the recurrence coefficients when the parameter $s$ approaches a breaking curve, by considering double scaling limits as $s$ approaches these points. We shall see a qualitative difference in the behavior of the recurrence coefficients, depending on whether or not we are approaching the points $s=\pm 2$ or some other points on the breaking curve.


title: Unraveling the role of V-V dimer on the vibrational properties of VO$_2$ by first-principles simulations and Raman spectroscopic analysis

abstract: We investigate the vibrational properties of VO2, particularly the low temperature M1 phase by first-principles calculations using the density functional theory as well as Raman spectroscopy. We perform the structural optimization using SCAN meta-GGA functional and obtain the optimized crystal structures for metallic rutile and insulating M1 phases satisfying all expected features of the experimentally derived structures. Based on the harmonic approximation around the optimized structures at zero temperature, we calculate the phonon properties and compare our results with experiments. We show that our calculated phonon density of states is in excellent agreement with the previous neutron scattering experiment. Moreover, we reproduce the phonon softening in the rutile phase as well as the phonon stiffening in the M1 phase. By comparing with the Raman experiments, we find that the Raman-active vibration modes of the M1 phase is strongly correlated with the V-V dimer distance of the crystal structure. Our combined theoretical and experimental framework demonstrates that Raman spectroscopy could serve as a reliable way to detect the subtle change of V-V dimer in the strained VO$_2$.


title: Expansions in the local and the central limit theorems for dynamical systems

abstract: We study higher order expansions both in the Berry-Ess\'een estimate (Edgeworth expansions) and in the local limit theorems for Birkhoff sums of chaotic probability preserving dynamical systems. We establish general results under technical assumptions, discuss the verification of these assumptions and illustrate our results by different examples (subshifts of finite type, Young towers, Sinai billiards, random matrix products), including situations of unbounded observables with integrability order arbitrarily close to the optimal moment condition required in the i.i.d. setting.


title: Assigning function to protein-protein interactions: a weakly supervised BioBERT based approach using PubMed abstracts

abstract: Motivation: Protein-protein interactions (PPI) are critical to the function of proteins in both normal and diseased cells, and many critical protein functions are mediated by interactions.Knowledge of the nature of these interactions is important for the construction of networks to analyse biological data. However, only a small percentage of PPIs captured in protein interaction databases have annotations of function available, e.g. only 4% of PPI are functionally annotated in the IntAct database. Here, we aim to label the function type of PPIs by extracting relationships described in PubMed abstracts. Method: We create a weakly supervised dataset from the IntAct PPI database containing interacting protein pairs with annotated function and associated abstracts from the PubMed database. We apply a state-of-the-art deep learning technique for biomedical natural language processing tasks, BioBERT, to build a model - dubbed PPI-BioBERT - for identifying the function of PPIs. In order to extract high quality PPI functions at large scale, we use an ensemble of PPI-BioBERT models to improve uncertainty estimation and apply an interaction type-specific threshold to counteract the effects of variations in the number of training samples per interaction type. Results: We scan 18 million PubMed abstracts to automatically identify 3253 new typed PPIs, including phosphorylation and acetylation interactions, with an overall precision of 46% (87% for acetylation) based on a human-reviewed sample. This work demonstrates that analysis of biomedical abstracts for PPI function extraction is a feasible approach to substantially increasing the number of interactions annotated with function captured in online databases.


title: Total Cost of Ownership Optimization for Direct Air-to-Ground System Design

abstract: Aircraft cabin is one of the last venues without mobile broadband. Considering future 5G applications and connectivity requirements, direct air-to-ground communications (DA2GC) is the only technique which can provide high capacity and low latency backhaul link for aircraft via a direct communication link. To this end, we propose an analytical framework to investigate the ground station deployment problem for DA2GC network employing multi-user beamforming with dual-polarized hybrid DA2GC antenna arrays. In addition, the proposed framework is utilized to analyze and optimize the total cost of ownership (TCO) of the DA2GC network to provide coverage for European airspace. We present the interplay between different network parameters: the number of ground stations, array size, transmit power and bandwidth, and TCO optimizing deployment parameters are calculated in order to satisfy capacity requirements. At the end, we show that, depending on the cost of different network resources, a terrestrial cellular network can be designed to cover the whole European airspace with limited number of ground stations with a certain array size, i.e., 900 and 361 antenna elements for ground station and air station, respectively.


title: Gluonic Hot Spot Initial Conditions in Heavy-Ion Collisions

abstract: The initial conditions in heavy-ion collisions are calculated in many different frameworks. The importance of nucleon position fluctuations within the nucleus and sub-nucleon structure has been established when modeling initial conditions for input to hydrodynamic calculations. However, there remain outstanding puzzles regarding these initial conditions, including the measurement of the near equivalence of the elliptical $v_{2}$ and triangular $v_{3}$ flow coefficients in ultra-central 0-1% Pb+Pb collisions at the LHC. Recently a calculation termed MAGMA incorporating gluonic hot spots via two-point correlators in the Color Glass Condensate framework, and no nucleons, provided a simultaneous match to these flow coefficients measured by the ATLAS experiment, including in ultra-central 0-1% collisions. Our calculations reveal that the MAGMA initial conditions do not describe the experimental data when run through full hydrodynamic SONIC simulations or when the hot spots from one nucleus resolve hot spots from the other nucleus, as predicted in the Color Glass Condensate framework. We also explore alternative initial condition calculations and discuss their implications.


title: iPhantom: a framework for automated creation of individualized computational phantoms and its application to CT organ dosimetry

abstract: Objective: This study aims to develop and validate a novel framework, iPhantom, for automated creation of patient-specific phantoms or digital-twins (DT) using patient medical images. The framework is applied to assess radiation dose to radiosensitive organs in CT imaging of individual patients. Method: From patient CT images, iPhantom segments selected anchor organs (e.g. liver, bones, pancreas) using a learning-based model developed for multi-organ CT segmentation. Organs challenging to segment (e.g. intestines) are incorporated from a matched phantom template, using a diffeomorphic registration model developed for multi-organ phantom-voxels. The resulting full-patient phantoms are used to assess organ doses during routine CT exams. Result: iPhantom was validated on both the XCAT (n=50) and an independent clinical (n=10) dataset with similar accuracy. iPhantom precisely predicted all organ locations with good accuracy of Dice Similarity Coefficients (DSC) >0.6 for anchor organs and DSC of 0.3-0.9 for all other organs. iPhantom showed less than 10% dose errors for the majority of organs, which was notably superior to the state-of-the-art baseline method (20-35% dose errors). Conclusion: iPhantom enables automated and accurate creation of patient-specific phantoms and, for the first time, provides sufficient and automated patient-specific dose estimates for CT dosimetry. Significance: The new framework brings the creation and application of CHPs to the level of individual CHPs through automation, achieving a wider and precise organ localization, paving the way for clinical monitoring, and personalized optimization, and large-scale research.


title: GPR-based Subsurface Object Detection and Reconstruction Using Random Motion and DepthNet

abstract: Ground Penetrating Radar (GPR) is one of the most important non-destructive evaluation (NDE) devices to detect the subsurface objects (i.e. rebars, utility pipes) and reveal the underground scene. One of the biggest challenges in GPR based inspection is the subsurface targets reconstruction. In order to address this issue, this paper presents a 3D GPR migration and dielectric prediction system to detect and reconstruct underground targets. This system is composed of three modules: 1) visual inertial fusion (VIF) module to generate the pose information of GPR device, 2) deep neural network module (i.e., DepthNet) which detects B-scan of GPR image, extracts hyperbola features to remove the noise in B-scan data and predicts dielectric to determine the depth of the objects, 3) 3D GPR migration module which synchronizes the pose information with GPR scan data processed by DepthNet to reconstruct and visualize the 3D underground targets. Our proposed DepthNet processes the GPR data by removing the noise in B-scan image as well as predicting depth of subsurface objects. For DepthNet model training and testing, we collect the real GPR data in the concrete test pit at Geophysical Survey System Inc. (GSSI) and create the synthetic GPR data by using gprMax3.0 simulator. The dataset we create includes 350 labeled GPR images. The DepthNet achieves an average accuracy of 92.64% for B-scan feature detection and an 0.112 average error for underground target depth prediction. In addition, the experimental results verify that our proposed method improve the migration accuracy and performance in generating 3D GPR image compared with the traditional migration methods.


title: Discovery of Light-induced Metastable Martensitic Anomaly Controlled by Single-Cycle Terahertz Pulses

abstract: We report on an ultrafast photoinduced phase transition with a strikingly long-lived Martensitic anomaly driven by above-threshold single-cycle terahertz (THz) pulses in Nb$_3$Sn. A non-thermal, THz-induced depletion of low frequency conductivity indicates increased gap splitting of high energy $\Gamma_{12}$ bands by removal of their degeneracies which enhances the Martensitic phase. In contrast, optical pumping leads to a $\Gamma_{12}$ gap melting. Such light-induced non-equilibrium Martensitic instability persists up to a critical temperature $\sim$100 K, i.e., more than twice the equilibrium temperature, and can be stabilized beyond technologically-relevant, nanosecond timescales. Together with first-principle simulations, we identify a compelling THz tuning of structural fluctuations via E$_u$ phonons to achieve a non-equilibrium ordering at high temperatures far exceeding those for equilibrium states.


title: Optimal Network Compression

abstract: This paper introduces a formulation of the optimal network compression problem for financial systems. This general formulation is presented for different levels of network compression or rerouting allowed from the initial interbank network. We prove that this problem is, generically, NP-hard. We focus on objective functions generated by systemic risk measures under shocks to the financial network. We use this framework to study the (sub)optimality of the maximally compressed network. We conclude by studying the optimal compression problem for specific networks; this permits us to study, e.g., the so-called robust fragility of certain network topologies more generally as well as the potential benefits and costs of network compression. In particular, under systematic shocks and heterogeneous financial networks the robust fragility results of Acemoglu et al. (2015) no longer hold generally.


title: Model-free optimal control of discrete-time systems with additive and multiplicative noises

abstract: This paper investigates the optimal control problem for a class of discrete-time stochastic systems subject to additive and multiplicative noises. A stochastic Lyapunov equation and a stochastic algebra Riccati equation are established for the existence of the optimal admissible control policy. A model-free reinforcement learning algorithm is proposed to learn the optimal admissible control policy using the data of the system states and inputs without requiring any knowledge of the system matrices. It is proven that the learning algorithm converges to the optimal admissible control policy. The implementation of the model-free algorithm is based on batch least squares and numerical average. The proposed algorithm is illustrated through a numerical example, which shows our algorithm outperforms other policy iteration algorithms.


title: Simultaneously-Collected Multimodal Lying Pose Dataset: Towards In-Bed Human Pose Monitoring under Adverse Vision Conditions

abstract: Computer vision (CV) has achieved great success in interpreting semantic meanings from images, yet CV algorithms can be brittle for tasks with adverse vision conditions and the ones suffering from data/label pair limitation. One of this tasks is in-bed human pose estimation, which has significant values in many healthcare applications. In-bed pose monitoring in natural settings could involve complete darkness or full occlusion. Furthermore, the lack of publicly available in-bed pose datasets hinders the use of many successful pose estimation algorithms for this task. In this paper, we introduce our Simultaneously-collected multimodal Lying Pose (SLP) dataset, which includes in-bed pose images from 109 participants captured using multiple imaging modalities including RGB, long wave infrared, depth, and pressure map. We also present a physical hyper parameter tuning strategy for ground truth pose label generation under extreme conditions such as lights off and being fully covered by a sheet/blanket. SLP design is compatible with the mainstream human pose datasets, therefore, the state-of-the-art 2D pose estimation models can be trained effectively with SLP data with promising performance as high as 95% at PCKh@0.5 on a single modality. The pose estimation performance can be further improved by including additional modalities through collaboration.


title: Unconventional four-terminal thermoelectric transport due to inelastic transport: cooling by transverse current, transverse thermoelectric effect and Maxwell demon

abstract: We show that in mesoscopic four-terminal thermoelectric devices with two electrodes (the source and the drain) and two heat baths, inelastic scattering processes can lead to unconventional thermoelectric transport. The source (or the drain) can be cooled by passing a thermal current between the two heat baths, with no net heat exchange between the heat baths and the electrodes. This effect, termed as cooling by heat current, is a mesoscopic heat drag effect. In addition, there is a transverse thermoelectric effect where electrical current and power can be generated by a transverse temperature bias (i.e., the temperature bias between the two heat baths). This transverse thermoelectric effect, originates from inelastic scattering processes, may have advantages for improved figures of merit and power factor due to spatial separation of charge and heat transport. We study the Onsager current-affinity relations, the linear-response transport properties, and the transverse thermoelectric figure of merit of the four-terminal thermoelectric devices for various system parameters. In addition, we investigate the efficiency and power of the cooling by transverse current effect in both linear and nonlinear transport regimes. We also demonstrate that by exploiting the inelastic transport in the quantum-dot four-terminal systems, a type of Maxwell's demon can be realized using nonequilibrium heat baths.


title: The Koopman Expectation: An Operator Theoretic Method for Efficient Analysis and Optimization of Uncertain Hybrid Dynamical Systems

abstract: For dynamical systems involving decision making, the success of the system greatly depends on its ability to make good decisions with incomplete and uncertain information. By leveraging the Koopman operator and its adjoint property, we introduce the Koopman Expectation, an efficient method for computing expectations as propagated through a dynamical system. Unlike other Koopman operator-based approaches in the literature, this is possible without an explicit representation of the Koopman operator. Furthermore, the efficiencies enabled by the Koopman Expectation are leveraged for optimization under uncertainty when expected losses and constraints are considered. We show how the Koopman Expectation is applicable to discrete, continuous, and hybrid non-linear systems driven by process noise with non-Gaussian initial condition and parametric uncertainties. We finish by demonstrating a 1700x acceleration for calculating probabilistic quantities of a hybrid dynamical system over the naive Monte Carlo approach with many orders of magnitudes improvement in accuracy.


title: Hidden Devil's staircase in a two-dimensional elastic model of spin crossover materials

abstract: Spin crossover (SCO) materials are reversible molecular switches found in a wide range of transition metal complexes and metal organic frameworks (MOFs). They exhibit diverse spin state orderings and transitions between them. We present an exact mapping from an elastic lattice mismatch model to a long-range Ising model, with an inverse square decay of the interaction strengths at large distances (on the square lattice). This provides a microscopic justification for an Ising model description, which has previously only been justified on phenomenological grounds. Elastic frustration is required for non-zero Ising interactions, but whether or not the short-range interactions in the Ising model are geometrically frustrated depends on the ratio of the bulk and shear moduli or equivalently Poisson's ratio. We show that, for a simple square lattice model with realistic parameters, sweeping the enthalpy difference between the two spin-states at zero temperature leads to a large (probably infinite) number of spin-state orderings and corresponding steps in the fraction of high-spin ions, consistent with a Devil's staircase. The staircase can also be climbed by varying the temperature, but then some of the steps are hidden and only a finite number remain, consistent with experiments on relevant framework materials, such as {(Fe[Hg(SCN)$_3$]$_2$(4,4'-bipy)$_2$)}$_n$. Our results are also relevant to other binary systems with lattice mismatch, e.g., heterogeneous solids.

categories: cond-mat.str-el cond-mat.other cond-mat.stat-mech physics.chem-ph


title: Non-Mergeable Sketching for Cardinality Estimation

abstract: Cardinality estimation is perhaps the simplest non-trivial statistical problem that can be solved via sketching. Industrially-deployed sketches like HyperLogLog, MinHash, and PCSA are mergeable, which means that large data sets can be sketched in a distributed environment, and then merged into a single sketch of the whole data set. In the last decade a variety of sketches have been developed that are non-mergeable, but attractive for other reasons. They are simpler, their cardinality estimates are strictly unbiased, and they have substantially lower variance. We evaluate sketching schemes on a reasonably level playing field, in terms of their memory-variance product (MVP). E.g., a sketch that occupies $5m$ bits and whose relative variance is $2/m$ (standard error $\sqrt{2/m}$) has an MVP of $10$. Our contributions are as follows. Cohen and Ting independently discovered what we call the Martingale transform for converting a mergeable sketch into a non-mergeable sketch. We present a simpler way to analyze the limiting MVP of Martingale-type sketches. We prove that the \Martingale{} transform is optimal in the non-mergeable world, and that \Martingale{} \fishmonger{} in particular is optimal among linearizable sketches, with an MVP of $H_0/2 \approx 1.63$. E.g., this is circumstantial evidence that to achieve 1\% standard error, we cannot do better than a 2 kilobyte sketch. \Martingale{} \fishmonger{} is neither simple nor practical. We develop a new mergeable sketch called \Curtain{} that strikes a nice balance between simplicity and efficiency, and prove that \Martingale{} \Curtain{} has limiting $\MVP\approx 2.31$. It can be updated with $O(1)$ memory accesses and it has lower empirical variance than \Martingale{} \LogLog, a practical non-mergeable version of HyperLogLog.


title: Energy-efficiency of Massive Random Access with Individual Codebook

abstract: The massive machine-type communication has been one of the most representative services for future wireless networks. It aims to support massive connectivity of user equipments (UEs) which sporadically transmit packets with small size. In this work, we assume the number of UEs grows linearly and unboundedly with blocklength and each UE has an individual codebook. Among all UEs, an unknown subset of UEs are active and transmit a fixed number of data bits to a base station over a shared-spectrum radio link. Under these settings, we derive the achievability and converse bounds on the minimum energy-per-bit for reliable random access over quasi-static fading channels with and without channel state information (CSI) at the receiver. These bounds provide energy-efficiency guidance for new schemes suited for massive random access. Simulation results indicate that the orthogonalization scheme TDMA is energy-inefficient for large values of UE density $\mu$. Besides, the multi-user interference can be perfectly cancelled when $\mu$ is below a critical threshold. In the case of no-CSI, the energy-per-bit for random access is only a bit more than that with the knowledge UE activity.


title: Functional Data Analysis with Causation in Observational Studies: Covariate Balancing Functional Propensity Score for Functional Treatments

abstract: Functional data analysis, which handles data arising from curves, surfaces, volumes, manifolds and beyond in a variety of scientific fields, is a rapidly developing area in modern statistics and data science in the recent decades. The effect of a functional variable on an outcome is an essential theme in functional data analysis, but a majority of related studies are restricted to correlational effects rather than causal effects. This paper makes the first attempt to study the causal effect of a functional variable as a treatment in observational studies. Despite the lack of a probability density function for the functional treatment, the propensity score is properly defined in terms of a multivariate substitute. Two covariate balancing methods are proposed to estimate the propensity score, which minimize the correlation between the treatment and covariates. The appealing performance of the proposed method in both covariate balance and causal effect estimation is demonstrated by a simulation study. The proposed method is applied to study the causal effect of body shape on human visceral adipose tissue.


title: Massive Unsourced Random Access for Massive MIMO Correlated Channels

abstract: This paper investigates the massive random access for a huge amount of user devices served by a base station (BS) equipped with a massive number of antennas. We consider a grant-free unsourced random access (U-RA) scheme where all users possess the same codebook and the BS aims at declaring a list of transmitted codewords and recovering the messages sent by active users. Most of the existing works concentrate on applying U-RA in the oversimplified independent and identically distributed (i.i.d.) channels. In this paper, we consider a fairly general joint-correlated MIMO channel model with line-of-sight components for the realistic outdoor wireless propagation environments. We conduct the activity detection for the emitted codewords by performing an improved coordinate descent approach with Bayesian learning automaton to solve a covariance-based maximum likelihood estimation problem. The proposed algorithm exhibits a faster convergence rate than traditional descent approaches. We further employ a coupled coding scheme to resolve the issue that the dimensions of the common codebook expand exponentially with user payload size in the practical massive machine-type communications scenario. Our simulations reveal that to achieve an error probability of 0.05 for reliable communications in correlated channels, one must pay a 0.9 to 1.3 dB penalty comparing to the minimum signal to noise ratio needed in i.i.d. channels on condition that a sufficient number of receiving antennas is equipped at the BS.


title: Understanding the Advisor-advisee Relationship via Scholarly Data Analysis

abstract: Advisor-advisee relationship is important in academic networks due to its universality and necessity. Despite the increasing desire to analyze the career of newcomers, however, the outcomes of different collaboration patterns between advisors and advisees remain unknown. The purpose of this paper is to find out the correlation between advisors' academic characteristics and advisees' academic performance in Computer Science. Employing both quantitative and qualitative analysis, we find that with the increase of advisors' academic age, advisees' performance experiences an initial growth, follows a sustaining stage, and finally ends up with a declining trend. We also discover the phenomenon that accomplished advisors can bring up skilled advisees. We explore the conclusion from two aspects: (1) Advisees mentored by advisors with high academic level have better academic performance than the rest; (2) Advisors with high academic level can raise their advisees' h-index ranking. This work provides new insights on promoting our understanding of the relationship between advisors' academic characteristics and advisees' performance, as well as on advisor choosing.


title: A comparison of mixed multiscale finite element methods for multiphase transport in highly heterogeneous media

abstract: In this paper, we systemically review and compare two mixed multiscale finite element methods (MMsFEM) for multiphase transport in highly heterogeneous media. In particular, we will consider the mixed multiscale finite element method using limited global information, simply denoted by MMsFEM, and the mixed generalized multiscale finite element method (MGMsFEM) with residual driven online multiscale basis functions. Both methods are under the framework of mixed multiscale finite element methods, where the pressure equation is solved in the coarse grid with carefully constructed multiscale basis functions for the velocity. The multiscale basis functions in both methods include local and global media information. In terms of MsFEM using limited global information, only one multiscale basis function is utilized in each local neighborhood while multiple basis are used in MGMsFEM. We will test and compare these two methods using the benchmark three-dimensional SPE10 model. A range of coarse grid sizes and different combinations of basis functions (offline and online) will be considered with CPU time reported for each case. In our numerical experiments, we observe good accuracy by the two above methods. Finally, we will discuss and compare the advantages and disadvantages of the two methods in terms of accuracy and computational costs.


title: High-Order Nonreciprocal Add-Drop Filter

abstract: Topological photonics have led to the robust optical behavior of the device, which has solved the problem of the influence of manufacturing defects and perturbations on the device performance. Meanwhile, temporal coupled-mode theory (t-CMT) has been developed and applied widely. However, the t-CMT of cascaded coupling cavities (CCC) system and its corresponding high-order filter has yet to be established. Here the t-CMT of CCC system is established based on the existing t-CMT. By combining the CCC with the topological waveguides, a versatile design scheme of the high-order nonreciprocal add-drop filter (HONAF) is proposed. The relationship between coupling effect of cavities and transmission and filtering performance of HONAF is analyzed quantitatively, then a method to improve the transmission efficiency and quality factor of the filter is given. Based on the combination of gyromagnetic photonic crystals and decagonal Penrose-type photonic quasicrystals, a HONAF is proposed. The transmission and filtering performance of the HONAF are numerically analyzed, which verifies the consistency between the theoretical prediction and the numerical simulation. The t-CMT of CCC system established can be widely used in the coupled resonator optical waveguides and their related systems. The designed HONAF can also be applied and compatible to microwave communication system.


title: Mie resonance engineering in meta-shell supraparticles for nanoscale nonlinear optics

abstract: Supraparticles are coordinated assemblies of discrete nanoscale building blocks into complex and hierarchical colloidal superstructures. Holistic optical responses in such assemblies are not observed in an individual building block or in their bulk counterparts. Furthermore, subwavelength dimensions of the unit building blocks enable engraving optical metamaterials within the supraparticle, which thus far has been beyond the current pool of colloidal engineering. This can lead to effective optical features in a colloidal platform with unprecedented ability to tune the electromagnetic responses of these particles. Here, we introduce and demonstrate the nanophotonics of meta-shell supraparticle (MSP), an all dielectric colloidal superstructure having an optical nonlinear metamaterial shell conformed onto a spherical core. We show that the metamaterial shell facilitates engineering the Mie resonances in the MSP that enable significant enhancement of the second harmonic generation (SHG). We show several orders of magnitude enhancement of second-harmonic generation in an MSP compared to its building blocks. Furthermore, we show an absolute conversion efficiency as high as 10^-7 far from the damage threshold, setting a new benchmark for SHG with low-index colloids. The MSP provides pragmatic solutions for instantaneous wavelength conversions with colloidal platforms that are suitable for chemical and biological applications. Their engineerability and scalability promise a fertile ground for nonlinear nanophotonics in the colloidal platforms with structural and material diversity.


title: Tuning transport coefficients of monolayer $\mathrm{MoSi_2N_4}$ with biaxial strain

abstract: Experimentally synthesized $\mathrm{MoSi_2N_4}$ (\textcolor[rgb]{0.00,0.00,1.00}{Science 369, 670-674 (2020)}) is a piezoelectric semiconductor. Here, we systematically study the large biaxial (isotropic) strain effects (0.90 to 1.10) on electronic structures and transport coefficients of monolayer $\mathrm{MoSi_2N_4}$ by density functional theory (DFT). With $a/a_0$ from 0.90 to 1.10, the energy band gap firstly increases, and then decreases, which is due to transformation of conduction band minimum (CBM). Calculated results show that the $\mathrm{MoSi_2N_4}$ monolayer is mechanically stable in considered strain range. It is found that the spin-orbital coupling (SOC) effects on Seebeck coefficient depend on the strain. In unstrained $\mathrm{MoSi_2N_4}$, the SOC has neglected influence on Seebeck coefficient. However, the SOC can produce important influence on Seebeck coefficient, when the strain is applied, for example 0.96 strain. The compressive strain can change relative position and numbers of conduction band extrema (CBE), and then the strength of conduction bands convergence can be enhanced, to the benefit of n-type $ZT_e$. Only about 0.96 strain can effectively improve n-type $ZT_e$. Our works imply that strain can effectively tune the electronic structures and transport coefficients of monolayer $\mathrm{MoSi_2N_4}$, and can motivate farther experimental exploration.


title: DPMC: Weighted Model Counting by Dynamic Programming on Project-Join Trees

abstract: We propose a unifying dynamic-programming framework to compute exact literal-weighted model counts of formulas in conjunctive normal form. At the center of our framework are project-join trees, which specify efficient project-join orders to apply additive projections (variable eliminations) and joins (clause multiplications). In this framework, model counting is performed in two phases. First, the planning phase constructs a project-join tree from a formula. Second, the execution phase computes the model count of the formula, employing dynamic programming as guided by the project-join tree. We empirically evaluate various methods for the planning phase and compare constraint-satisfaction heuristics with tree-decomposition tools. We also investigate the performance of different data structures for the execution phase and compare algebraic decision diagrams with tensors. We show that our dynamic-programming model-counting framework DPMC is competitive with the state-of-the-art exact weighted model counters cachet, c2d, d4, and miniC2D.


title: Monoclinic YFe$_{12}$ phases predicted from first principles

abstract: We searched for stable crystal structures of YFe$_{12}$ using a crystal structure prediction technique based on a genetic algorithm and first-principles calculations. We obtained two monoclinic $C2/m$ structures as metastable phases that are different from the well-known ThMn$_{12}$ structure. These two phases have advantages in their magnetism over the ThMn$_{12}$ structure: The total magnetization $M$ is increased from 25.6 $\mu_{\text{B}}$/f.u. up to 26.8 $\mu_{\text{B}}$/f.u. by the transformations. We also calculated Curie temperature $T_{\text{C}}$ for these structures within the mean-field approximation and predicted the increase of $T_{\text{C}}$ from 792 K up to 940 K, which is mainly caused by the increase of intersite magnetic couplings within the distance of 2.3--3.1\r{A}. The similar enhancements of $M$ and $T_{\text{C}}$ are also obtained in the pseudo-binary system Y(Fe$_{1-x}$Co$_{x}$)$_{12}$ with $x$ of 0--0.7.


title: Prototype-based interpretation of the functionality of neurons in winner-take-all neural networks

abstract: Prototype-based learning (PbL) using a winner-take-all (WTA) network based on minimum Euclidean distance (ED-WTA) is an intuitive approach to multiclass classification. By constructing meaningful class centers, PbL provides higher interpretability and generalization than hyperplane-based learning (HbL) methods based on maximum Inner Product (IP-WTA) and can efficiently detect and reject samples that do not belong to any classes. In this paper, we first prove the equivalence of IP-WTA and ED-WTA from a representational point of view. Then, we show that naively using this equivalence leads to unintuitive ED-WTA networks in which the centers have high distances to data that they represent. We propose $\pm$ED-WTA which models each neuron with two prototypes: one positive prototype representing samples that are modeled by this neuron and a negative prototype representing the samples that are erroneously won by that neuron during training. We propose a novel training algorithm for the $\pm$ED-WTA network, which cleverly switches between updating the positive and negative prototypes and is essential to the emergence of interpretable prototypes. Unexpectedly, we observed that the negative prototype of each neuron is indistinguishably similar to the positive one. The rationale behind this observation is that the training data that are mistaken with a prototype are indeed similar to it. The main finding of this paper is this interpretation of the functionality of neurons as computing the difference between the distances to a positive and a negative prototype, which is in agreement with the BCM theory. In our experiments, we show that the proposed $\pm$ED-WTA method constructs highly interpretable prototypes that can be successfully used for detecting outlier and adversarial examples.


title: Exact diagonalization study of the anisotropic Heisenberg model related to YbMgGaO$_4$

abstract: Employing exact diagonalization, we systematically study the anisotropic Heisenberg model which is related to rare-earth triangular-lattice materials. We probe its full 3D phase diagram afresh and identify a large region of quantum spin liquid (QSL) phase which can extend to the QSL region of the $J_1$--$J_2$ triangular Heisenberg model. Furthermore, we explore the magnetization curves of different phases and reproduce the 1/3-magnetization plateau in the quantum spin liquid phase region. More importantly, to study the possible chemical disorders in real materials, we consider the randomness of exchange interactions and find no spin glass order. And there is a large region of random-singlet phase which contains strongly random spin networks, dominated by two-spin singlets, four-spin singlets and other singlet domains. Our comprehensive ED study can give detailed insightful understanding of the microscopic Hamiltonian related to the YbMgGaO$_4$ and some other related rare-earth triangular-lattice materials.


title: Steady Rayleigh--B\'enard convection between no-slip boundaries

abstract: The central open question about Rayleigh--B\'enard convection -- buoyancy-driven flow in a fluid layer heated from below and cooled from above -- is how vertical heat flux depends on the imposed temperature gradient in the strongly nonlinear regime where the flows are typically turbulent. The quantitative challenge is to determine how the Nusselt number $Nu$ depends on the Rayleigh number $Ra$ in the $Ra\to\infty$ limit for fluids of fixed finite Prandtl number $Pr$ in fixed spatial domains. Laboratory experiments, numerical simulations, and analysis of Rayleigh's mathematical model have yet to rule out either of the proposed `classical' $Nu \sim Ra^{1/3}$ or `ultimate' $Nu \sim Ra^{1/2}$ asymptotic scaling theories. Among the many solutions of the equations of motion at high $Ra$ are steady convection rolls that are dynamically unstable but share features of the turbulent attractor. We have computed these steady solutions for $Ra$ up to $10^{14}$ with $Pr=1$ and various horizontal periods. By choosing the horizontal period of these rolls at each $Ra$ to maximize $Nu$, we find that steady convection rolls achieve classical asymptotic scaling. Moreover, they transport more heat than turbulent convection in experiments or simulations at comparable parameters. If heat transport in turbulent convection continues to be dominated by heat transport in steady rolls as $Ra\to\infty$, it cannot achieve the ultimate scaling.


title: When Homomorphic Encryption Marries Secret Sharing: Secure Large-Scale Sparse Logistic Regression and Applications in Risk Control

abstract: Logistic Regression (LR) is the most widely used machine learning model in industry for its efficiency, robustness, and interpretability. Due to the problem of data isolation and the requirement of high model performance, many applications in industry call for building a secure and efficient LR model for multiple parties. Most existing work uses either Homomorphic Encryption (HE) or Secret Sharing (SS) to build secure LR. HE based methods can deal with high-dimensional sparse features, but they incur potential security risks. SS based methods have provable security, but they have efficiency issue under high-dimensional sparse features. In this paper, we first present CAESAR, which combines HE and SS to build secure large-scale sparse logistic regression model and achieves both efficiency and security. We then present the distributed implementation of CAESAR for scalability requirement. We have deployed CAESAR in a risk control task and conducted comprehensive experiments. Our experimental results show that CAESAR improves the state-of-the-art model by around 130 times.


title: Generalizing the de Finetti--Hewitt--Savage theorem

abstract: A sequence of random variables is called exchangeable if its joint distribution is invariant under permutations. The original formulation of de Finetti's theorem says that any exchangeable sequence of $\{0,1\}$-valued random variables can be thought of as a mixture of independent and identically distributed sequences in a certain precise mathematical sense. Interpreting this statement from a convex analytic perspective, Hewitt and Savage obtained the same conclusion for more general state spaces under some topological conditions. The main contribution of this paper is in providing a new framework that explains the theorem purely as a consequence of the underlying distribution of the random variables, with no topological conditions (beyond Hausdorffness) on the state space being necessary if the distribution is Radon. We also show that it is consistent with the axioms of ZFC that de Finetti's theorem holds for all sequences of exchangeable random variables taking values in any complete metric space. The framework we use is based on nonstandard analysis. We have provided a self-contained introduction to nonstandard analysis as an appendix, thus rendering measure theoretic probability and point-set topology as the only prerequisites for this paper. Our introduction aims to develop some new ideologies that might be of interest to mathematicians, philosophers, and mathematics educators alike. Our technical tools come from nonstandard topological measure theory, in which a highlight is a new generalization of Prokhorov's theorem. Modulo such technical tools, our proof relies on properties of the empirical measures induced by hyperfinitely many identically distributed random variables -- a feature that allows us to establish de Finetti's theorem in the generality that we seek while still retaining the combinatorial intuition of proofs of simpler versions of de Finetti's theorem.


title: On $\ell_p$-norm Robustness of Ensemble Stumps and Trees

abstract: Recent papers have demonstrated that ensemble stumps and trees could be vulnerable to small input perturbations, so robustness verification and defense for those models have become an important research problem. However, due to the structure of decision trees, where each node makes decision purely based on one feature value, all the previous works only consider the $\ell_\infty$ norm perturbation. To study robustness with respect to a general $\ell_p$ norm perturbation, one has to consider the correlation between perturbations on different features, which has not been handled by previous algorithms. In this paper, we study the problem of robustness verification and certified defense with respect to general $\ell_p$ norm perturbations for ensemble decision stumps and trees. For robustness verification of ensemble stumps, we prove that complete verification is NP-complete for $p\in(0, \infty)$ while polynomial time algorithms exist for $p=0$ or $\infty$. For $p\in(0, \infty)$ we develop an efficient dynamic programming based algorithm for sound verification of ensemble stumps. For ensemble trees, we generalize the previous multi-level robustness verification algorithm to $\ell_p$ norm. We demonstrate the first certified defense method for training ensemble stumps and trees with respect to $\ell_p$ norm perturbations, and verify its effectiveness empirically on real datasets.


title: iCaps: An Interpretable Classifier via Disentangled Capsule Networks

abstract: We propose an interpretable Capsule Network, iCaps, for image classification. A capsule is a group of neurons nested inside each layer, and the one in the last layer is called a class capsule, which is a vector whose norm indicates a predicted probability for the class. Using the class capsule, existing Capsule Networks already provide some level of interpretability. However, there are two limitations which degrade its interpretability: 1) the class capsule also includes classification-irrelevant information, and 2) entities represented by the class capsule overlap. In this work, we address these two limitations using a novel class-supervised disentanglement algorithm and an additional regularizer, respectively. Through quantitative and qualitative evaluations on three datasets, we demonstrate that the resulting classifier, iCaps, provides a prediction along with clear rationales behind it with no performance degradation.


title: On Lower Bounds for Standard and Robust Gaussian Process Bandit Optimization

abstract: In this paper, we consider algorithm-independent lower bounds for the problem of black-box optimization of functions having a bounded norm is some Reproducing Kernel Hilbert Space (RKHS), which can be viewed as a non-Bayesian Gaussian process bandit problem. In the standard noisy setting, we provide a novel proof technique for deriving lower bounds on the regret, with benefits including simplicity, versatility, and an improved dependence on the error probability. In a robust setting in which every sampled point may be perturbed by a suitably-constrained adversary, we provide a novel lower bound for deterministic strategies, demonstrating an inevitable joint dependence of the cumulative regret on the corruption level and the time horizon, in contrast with existing lower bounds that only characterize the individual dependencies. Furthermore, in a distinct robust setting in which the final point is perturbed by an adversary, we strengthen an existing lower bound that only holds for target success probabilities very close to one, by allowing for arbitrary success probabilities above $\frac{2}{3}$.


title: Semi-Blind and l1 Robust System Identification for Anemia Management

abstract: Chronic diseases such as cancer, diabetes, heart diseases, chronic kidney disease (CKD) require a drug management system that ensures a stable and robust output of the patient's condition in response to drug dosage. In the case of CKD, the patients suffer from the deficiency of red blood cell count and external human recombinant erythropoietin (EPO) is required to maintain healthy levels of hemoglobin (Hb). Anemia is a common comorbidity in patients with CKD. For an efficient and robust anemia management system for CKD patients instead of traditional population-based approaches, individualized patient-specific approaches are needed. Hence, individualized system (patient) models for patient-specific drug-dose responses are required. In this research, system identification for CKD is performed for individual patients. For control-oriented system identification, two robust identification techniques are applied: (1) l1 robust identification considering zero initial conditions and (2) semi-blind robust system identification considering non-zero initial conditions. The EPO data of patients are used as the input and Hb data is used as the output of the system. For this study, individualized patient models are developed by using patient-specific data. The ARX one-step-ahead prediction technique is used for model validation at real patient data. The performance of these two techniques is compared by calculating minimum means square error (MMSE). By comparison, we show that the semi-blind robust identification technique gives better results as compared to l1 robust identification.


title: Positionality-Weighted Aggregation Methods for Cumulative Voting

abstract: Respecting minority opinions is vital in solving social problems. However, minority opinions are often ignored in general majority rules. To build consensus on pluralistic values and make social choices that consider minority opinions, we propose aggregation methods that give weighting to the minority's positionality on cardinal cumulative voting. Based on quadratic and linear voting, we formulated three weighted aggregation methods that differ in the ratio of votes to cumulative points and the weighting of the minority to all members, and assuming that the distributions of votes follow normal distributions, we calculated the frequency distributions of the aggregation results. We found that minority opinions are more likely to be reflected proportionately to the average of the distribution in two of the above three methods. This implies that Sen and Gotoh's idea of considering the social position of unfortunate people on ordinal ranking in the welfare economics, was illustrated by weighting the minority's positionality on cardinal voting. In addition, it is possible to visualize the number and positionality of the minority from the analysis of the aggregation results. These results will be useful to promote mutual understanding between the majority and minority by interactively visualizing the contents of the proposed aggregation methods in the consensus-building process. With the further development of information technology, the consensus building based on big data will be necessary. We recommend the use of our proposed aggregation methods to make social choices for pluralistic values such as social, environmental, and economic.


title: Corrections to "An Innovations Approach to Viterbi Decoding of Convolutional Codes"

abstract: We correct errors in Section III-B of the above-titled paper. The corrections are related to the joint distribution of the inputs to the main decoder which corresponds to a branch in the code trellis for the main decoder. Note that as far as the distribution of the input corresponding to a single code symbol is concerned, the results in Section III-B are correct. We have noticed that alpha in Proposition 12 (beta in Proposition 14) has another meaning, which has made it possible to derive the joint distribution of the inputs corresponding to the whole branch. Corrections are based on this important observation.


title: Reliable Traffic Monitoring Mechanisms Based on Blockchain in Vehicular Networks

abstract: The real-time traffic monitoring is a fundamental mission in a smart city to understand traffic conditions and avoid dangerous incidents. In this paper, we propose a reliable and efficient traffic monitoring system that integrates blockchain and the Internet of vehicles technologies effectively. It can crowdsource its tasks of traffic information collection to vehicles that run on the road instead of installing cameras in every corner. First, we design a lightweight blockchain-based information trading framework to model the interactions between traffic administration and vehicles. It guarantees reliability, efficiency, and security during executing trading. Second, we define the utility functions for the entities in this system and come up with a budgeted auction mechanism that motivates vehicles to undertake the collection tasks actively. In our algorithm, it not only ensures that the total payment to the selected vehicles does not exceed a given budget, but also maintains the truthfulness of auction process that avoids some vehicles to offer unreal bids for getting greater utilities. Finally, we conduct a group of numerical simulations to evaluate the reliability of our trading framework and performance of our algorithms, whose results demonstrate their correctness and efficiency perfectly.


title: Existence of partially hyperbolic motions in the N-body problem

abstract: In the context of the Newtonian N-body problem, we prove the existence of a partially hyperbolic motion with prescribed positive energy and any initial collisionless configuration. Moreover, it is a free time minimizer of the respective supercritical Newtonian action or equivalently a geodesic ray for the respective Jacobi-Maupertuis metric.


title: Scattering in the Ising Model Using Quantum Lanczos Algorithm

abstract: Time evolution and scattering simulation in phenomenological models are of great interest for testing and validating the potential for near-term quantum computers to simulate quantum field theories. Here, we simulate one-particle propagation and two-particle scattering in the one-dimensional transverse Ising model for 3 and 4 spatial sites with periodic boundary conditions on a quantum computer. We use the quantum Lanczos algorithm to obtain all energy levels and corresponding eigenstates of the system. We simplify the quantum computation by taking advantage of the symmetries of the system. These results enable us to compute one- and two-particle transition amplitudes, particle numbers for spatial sites, and the transverse magnetization as functions of time. The quantum circuits were executed on IBM 5-qubit superconducting hardware. The experimental results with readout error mitigation are in very good agreement with the values obtained using exact diagonalization.


title: Thermal Noise in Electro-Optic Devices at Cryogenic Temperatures

abstract: The quantum bits (qubits) on which superconducting quantum computers are based have energy scales corresponding to photons with GHz frequencies. The energy of photons in the gigahertz domain is too low to allow transmission through the noisy room-temperature environment, where the signal would be lost in thermal noise. Optical photons, on the other hand, have much higher energies, and signals can be detected using highly efficient single-photon detectors. Transduction from microwave to optical frequencies is therefore a potential enabling technology for quantum devices. However, in such a device the optical pump can be a source of thermal noise and thus degrade the fidelity; the similarity of input microwave state to the output optical state. In order to investigate the magnitude of this effect we model the sub-Kelvin thermal behavior of an electro-optic transducer based on a lithium niobate whispering gallery mode resonator. We find that there is an optimum power level for a continuous pump, whilst pulsed operation of the pump increases the fidelity of the conversion.


title: Kinetic Expression for Optimal Catalyst Electronic Configuration: The Case of Ammonia Decomposition

abstract: A new steady-state kinetic model of ammonia decomposition is presented and analyzed regarding the electronic properties of metal catalysts. The model is based on the classical Temkin-Ertl mechanism and modified in accordance with Wolkenstein electronic theory by implementing participation of free electrons of the catalyst to change the chemical nature of adsorbed species. Wolkenstein original theory only applied to semiconductors but by including the d-band model, the electronic theory can be extended to metals. For both simplified and full reaction mechanisms, including electronic steps, we present a steady-state rate equation where the dependence on the Fermi level of the metal creates a volcano-shaped dependence. According to the kinetic model, an increasing Fermi level of the catalyst, that approaching the antibonding state with adsorbed nitrogen molecules, will increase the fraction of neutral nitrogen molecules and enhance their the desorption. Concurrently, strong chemisorption of ammonia molecules proceeds easily through participation of additional free catalyst electrons in the adsorbate bond. As a result, the reaction rate is enhanced and reaches its maximum value. A further increasing Fermi level of the catalyst that approaches the antibonding state with ammonia molecules will result in a smaller fraction of negatively charged ammonia molecules and less dehydrogenation. Concurrently, the desorption of neutral nitrogen molecules occurs without impairment. As a result, the reaction rate decreases. The detailed kinetic model is compared to recent experimental measurements of ammonia decomposition on iron, cobalt and CoFe bimetallic catalysts.


title: Deformable PV-RCNN: Improving 3D Object Detection with Learned Deformations

abstract: We present Deformable PV-RCNN, a high-performing point-cloud based 3D object detector. Currently, the proposal refinement methods used by the state-of-the-art two-stage detectors cannot adequately accommodate differing object scales, varying point-cloud density, part-deformation and clutter. We present a proposal refinement module inspired by 2D deformable convolution networks that can adaptively gather instance-specific features from locations where informative content exists. We also propose a simple context gating mechanism which allows the keypoints to select relevant context information for the refinement stage. We show state-of-the-art results on the KITTI dataset.


title: Single Image Super-Resolution via a Holistic Attention Network

abstract: Informative features play a crucial role in the single image super-resolution task. Channel attention has been demonstrated to be effective for preserving information-rich features in each layer. However, channel attention treats each convolution layer as a separate process that misses the correlation among different layers. To address this problem, we propose a new holistic attention network (HAN), which consists of a layer attention module (LAM) and a channel-spatial attention module (CSAM), to model the holistic interdependencies among layers, channels, and positions. Specifically, the proposed LAM adaptively emphasizes hierarchical features by considering correlations among layers. Meanwhile, CSAM learns the confidence at all the positions of each channel to selectively capture more informative features. Extensive experiments demonstrate that the proposed HAN performs favorably against the state-of-the-art single image super-resolution approaches.


title: Bias of Root Numbers for Modular Newforms of Cubic Level

abstract: Let $H^{\pm}_{2k} (N^3)$ denote the set of modular newforms of cubic level $N^3$, weight $2 k$, and root number $\pm 1$. For $N > 1$ squarefree and $k>1$, we use an analytic method to establish neat and explicit formulas for the difference $|H^{+}_{2k} (N^3)| - |H^{-}_{2k} (N^3)|$ as a multiple of the product of $\varphi (N)$ and the class number of $\mathbb{Q}(\sqrt{- N})$. In particular, the formulas exhibit a strict bias towards the root number $+1$. Our main tool is a root-number weighted simple Petersson formula for such newforms.


title: Lite Training Strategies for Portuguese-English and English-Portuguese Translation

abstract: Despite the widespread adoption of deep learning for machine translation, it is still expensive to develop high-quality translation models. In this work, we investigate the use of pre-trained models, such as T5 for Portuguese-English and English-Portuguese translation tasks using low-cost hardware. We explore the use of Portuguese and English pre-trained language models and propose an adaptation of the English tokenizer to represent Portuguese characters, such as diaeresis, acute and grave accents. We compare our models to the Google Translate API and MarianMT on a subset of the ParaCrawl dataset, as well as to the winning submission to the WMT19 Biomedical Translation Shared Task. We also describe our submission to the WMT20 Biomedical Translation Shared Task. Our results show that our models have a competitive performance to state-of-the-art models while being trained on modest hardware (a single 8GB gaming GPU for nine days). Our data, models and code are available at https://github.com/unicamp-dl/Lite-T5-Translation.


title: Analysis of a nonlinear free-boundary tumor model with angiogenesis and a connection between the nonnecrotic and necrotic phases

abstract: This paper is concerned with a nonlinear free boundary problem modeling the growth of spherically symmetric tumors with angiogenesis, set with a Robin boundary condition. In which, both nonnecrotic tumors and necrotic tumors are taken into consideration. The well-posedness and asymptotic behavior of solutions are studied. It is shown that there exist two thresholds, denoted by $\tilde\sigma$ and $\sigma^*$, on the surrounding nutrient concentration $\bar\sigma$. If $\bar\sigma\leq\tilde\sigma$, then the considered problem admits no stationary solution and all evolutionary tumors will finally vanish, while if $\bar\sigma>\tilde\sigma$, then it admits a unique stationary solution and all evolutionary tumors will converge to this dormant tumor; moreover, the dormant tumor is nonnecrotic if $\tilde\sigma<\bar\sigma\leq\sigma^*$ and necrotic if $\bar\sigma>\sigma^*$. The connection and mutual transition between the nonnecrotic and necrotic phases are also given.


title: Modelling Elastically-Mediated Liquid-Liquid Phase Separation

abstract: We propose a continuum theory of the liquid-liquid phase separation in an elastic network where phase-separated microscopic droplets rich in one fluid component can form as an interplay of fluids mixing, droplet nucleation, network deformation, thermodynamic fluctuation, \emph{etc}. We find that the size of the phase separated droplets decreases with the shear modulus of the elastic network in the form of $\sim[\mathrm{modulus}]^{-1/3}$ and the number density of the droplet increases almost linearly with the shear modulus $\sim[\mathrm{modulus}]$, which are verified by the experimental observations. Phase diagrams in the space of (fluid constitution, mixture interaction, network modulus) are provided, which can help to understand similar phase separations in biological cells and also to guide fabrications of synthetic cells with desired phase properties.


title: Topological Insulator Metamaterial with Giant Circular Photogalvanic Effect

abstract: One of the most striking manifestations of electronic properties of topological insulators is the dependence of the photocurrent direction on the helicity of circularly polarized optical excitation. The helicity dependent photocurrents, underpinned by spin-momentum locking of surface Dirac electrons, are weak and easily overshadowed by bulk contributions. Here we show that the chiral response can be enhanced by nanostructuring. The tight confinement of electromagnetic fields in the resonant nanostructure enhances the photoexcitation of spin-polarized surface states of topological insulator Bi1.5Sb0.5Te1.8Se1.2, leading to an 11-fold increase of the circular photogalvanic effect and an unprecedented photocurrent dichroism (\r{ho}circ=0.87) at room temperature. The control of spin-transport in topological materials by structural design is a previously unrecognised ability of metamaterials that bridges the gap between nanophotonics and spin-electronics, providing new opportunities for developing polarization sensitive photodetectors.


title: Precision Metrology Meets Cosmology: Improved Constraints on Ultralight Dark Matter from Atom-Cavity Frequency Comparisons

abstract: We conduct frequency comparisons between a state-of-the-art strontium optical lattice clock, a cryogenic crystalline silicon cavity, and a hydrogen maser to set new bounds on the coupling of ultralight dark matter to Standard Model particles and fields in the mass range of $10^{-16}$ $-$ $10^{-21}$ eV. The key advantage of this two-part ratio comparison is the differential sensitivities to time variation of both the fine-structure constant and the electron mass, achieving a substantially improved limit on the moduli of ultralight dark matter, particularly at higher masses than typical atomic spectroscopic results. Furthermore, we demonstrate an extension of the search range to even higher masses by use of dynamical decoupling techniques. These results highlight the importance of using the best performing atomic clocks for fundamental physics applications as all-optical timescales are increasingly integrated with, and will eventually supplant, existing microwave timescales.

categories: physics.atom-ph astro-ph.IM hep-ph physics.optics


title: Super homologies associated with low dimensional Lie algebras

abstract: A Poisson structure on a manifold is characterized by the Schouten bracket. The graded algebra of the tangent bundle with the Schouten bracket is a prototype of Lie superalgebra. The Poisson condition means that a cycle in the 2-chain space. Given a graded Lie superalgebra, the 0-graded subspace is a Lie algebra. In this note, using the DGA of tangent bundle with the Schouten bracket as a model, we start from an abstract Lie algebra, construct non-trivial Lie superalgebra by Schouten-like bracket. Then it is natural to ask how the core Lie algebra control the Lie superalgebra. One trial here is to investigate the Betti numbers of the super homology groups. For abelian Lie algebras, the boundary operator is trivial, so we study super homology groups for low dimensional non-abelian Lie algebras of dimension smaller than 4.


title: Spatial--spectral FFPNet: Attention-Based Pyramid Network for Segmentation and Classification of Remote Sensing Images

abstract: We consider the problem of segmentation and classification of high-resolution and hyperspectral remote sensing images. Unlike conventional natural (RGB) images, the inherent large scale and complex structures of remote sensing images pose major challenges such as spatial object distribution diversity and spectral information extraction when existing models are directly applied for image classification. In this study, we develop an attention-based pyramid network for segmentation and classification of remote sensing datasets. Attention mechanisms are used to develop the following modules: i) a novel and robust attention-based multi-scale fusion method effectively fuses useful spatial or spectral information at different and same scales; ii) a region pyramid attention mechanism using region-based attention addresses the target geometric size diversity in large-scale remote sensing images; and iii cross-scale attention} in our adaptive atrous spatial pyramid pooling network adapts to varied contents in a feature-embedded space. Different forms of feature fusion pyramid frameworks are established by combining these attention-based modules. First, a novel segmentation framework, called the heavy-weight spatial feature fusion pyramid network (FFPNet), is proposed to address the spatial problem of high-resolution remote sensing images. Second, an end-to-end spatial--spectral FFPNet is presented for classifying hyperspectral images. Experiments conducted on ISPRS Vaihingen and ISPRS Potsdam high-resolution datasets demonstrate the competitive segmentation accuracy achieved by the proposed heavy-weight spatial FFPNet. Furthermore, experiments on the Indian Pines and the University of Pavia hyperspectral datasets indicate that the proposed spatial--spectral FFPNet outperforms the current state-of-the-art methods in hyperspectral image classification.


title: Radiation from an inertial mirror horizon

abstract: The purpose of this study is to investigate radiation from asymptotic zero acceleration motion where a horizon is formed and subsequently detected by an outside witness. A perfectly reflecting moving mirror is used to model such a system and compute the energy and spectrum. The trajectory is asymptotically inertial (zero proper acceleration)-ensuring negative energy flux (NEF), yet approaches light-speed with a null ray horizon at a finite advanced time. We compute the spectrum and energy analytically.


title: Compactness Theorem of Complete k-Curvature Manifolds with Isolated Singularities

abstract: In this paper we prove that the set of metrics conformal to the standard metric on $\mathbb{S}^{n}\backslash\{p_{1},\cdots,p_{l}\}$ is locally compact in $C^{m,\alpha}$ topology for any $m>0$, whenever the metrics have constant $\sigma_{k}$ curvature and the $k$-Dilational Pohozaev invariants have positive lower bound for $k<n/2$. Here the $k$-Dilational Pohozaev invariants come from the Kazdan-Warner type identity for the $\sigma_{k}$ curvature, which is derived by Viaclovsky \cite{Viac2000} and Han \cite{H1}. When $k=1$, Pollack \cite{Pollack} proved the compactness results for the complete metrics of constant positive scalar curvature on $\mathbb{S}^{n}\backslash\{p_{1},\cdots,p_{l}\}$.


title: Strong consistent model selection for general causal time series

abstract: We consider the strongly consistent question for model selection in a large class of causal time series models, including AR($\infty$), ARCH($\infty$), TARCH($\infty$), ARMA-GARCH and many classical others processes. We propose a penalized criterion based on the quasi likelihood of the model. We provide sufficient conditions that ensure the strong consistency of the proposed procedure. Also, the estimator of the parameter of the selected model obeys the law of iterated logarithm. It appears that, unlike the result of the weak consistency obtained by Bardet {\it et al.} \cite{Bardet2020}, a dependence between the regularization parameter and the model structure is not needed.


title: A simple 7/3-approximation algorithm for feedback vertex set in tournaments

abstract: We show that performing just one round of the Sherali-Adams hierarchy gives an easy 7/3-approximation algorithm for the Feedback Vertex Set (FVST) problem in tournaments. This matches the best deterministic approximation algorithm for FVST due to Mnich, Williams, and V\'egh, and is a significant simplification and runtime improvement of their approach.


title: Relativistic Transformation of Spherical Co ordinates(t,r,{\theta},{\phi})

abstract: With the advent of relativistic mechanics, the Lorentz transformation replaced the Galilean transformation based on classical Newtonian mechanics among inertial frames at high uniform velocities, but both transformations are based on Cartesian coordinate system, hence position of particles obtaining linear velocities in space can be obtained. In case where frames are rotating with constant angular velocity, use of Galilean rotational transformation (GRT) is replaced by Franklin transformation, proposed by Philip Franklin in 1922. The modified transformation introduced the concept of rotational motion of points in a rigid body. Both the transformations are based on cylindrical coordinate system. Here we moved a step further for making a relativistic transformation using spherical coordinate system for understanding the behaviour of rotating frames along any axis in the space passing through the center of mass of a symmetrical object (Sphere). We finally came to an understanding about how Special Theory of relativity is found to be applicable in rotational motion using different co-ordinate system.


title: Boron cage effects on Nd-Fe-B crystal structure's stability

abstract: In this study, we investigate the structure-stability relationship of hypothetical Nd-Fe-B crystal structures using descriptor-relevance analysis and the t-SNE dimensionality reduction method. 149 hypothetical Nd-Fe-B crystal structures are generated from 5967 LA-T-X host structures in Open Quantum Materials Database by using the elemental substitution method, with LA denoting lanthanides, T denoting transition metals, and X denoting light elements such as B, C, N and O. A hypothetical crystal structure is created by substituting all lanthanide sites with Nd, all transition metal sites with Fe, and all light element sites with B. High-throughput first-principle calculations are applied to evaluate the phase stability of these structures. Twenty of them are found to be potentially formable. The descriptor-relevance analysis on the orbital field matrix (OFM) materials' descriptor reveals the average atomic coordination number as the essential factor in determining the structure stability of these substituted Nd-Fe-B crystal structures. 19 among 20 hypothetical structures that are found potentially formable have an average coordination number larger than 6.5. In addition, all the local structures represented by the OFM descriptors are integrated into a visible space to study the detailed correlation between their characteristics and the stability of the crystal structure to which they belong. We discover that unstable substituted structures frequently carry Nd and Fe local structures with two prominent points: low average coordination numbers and fully occupied B neighboring atoms. Moreover, there are only three popular forms of B local structures appearing on all potentially formable substituted structures: cage networks, planar networks, and interstitial sites. The discovered relationships are promising to speed up the screening process for the new formable crystal structures.


title: The definable content of homological invariants I: $\mathrm{Ext}$ & $\mathrm{lim}^1$

abstract: This is the first installment in a series of papers in which we illustrate how classical invariants of homological algebra and algebraic topology can be enriched with additional descriptive set-theoretic information. To effect this enrichment, we show that many of these invariants can be naturally regarded as functors to the category, introduced herein, of groups with a Polish cover. The resulting definable invariants provide far stronger means of classification. In the present work we focus on the first derived functors of $\mathrm{Hom}(-,-)$ and $\mathrm{lim}(-)$. The resulting definable $\mathrm{Ext}(B,F)$ for pairs of countable abelian groups $B,F$ and definable $\mathrm{lim}^{1}(\boldsymbol{A})$ for towers $\boldsymbol{A}$ of Polish abelian groups substantially refine their classical counterparts. We show, for example, that the definable $\textrm{Ext}(-,\mathbb{Z})$ is a fully faithful contravariant functor from the category of finite rank torsion-free abelian groups $\Lambda$ with no free summands; this contrasts with the fact that there are uncountably many non-isomorphic such groups $\Lambda$ with isomorphic classical invariants $\textrm{Ext}(\Lambda,\mathbb{Z}) $. To facilitate our analysis, we introduce a general Ulam stability framework for groups with a Polish cover and we prove several rigidity results for non-Archimedean abelian groups with a Polish cover. A special case of our main result answers a question of Kanovei and Reeken regarding quotients of the $p$-adic groups. Finally, using cocycle superrigidity methods for profinite actions of property (T) groups, we obtain a hierarchy of complexity degrees for the problem $\mathcal{R}(\mathrm{Aut}(\Lambda)\curvearrowright\mathrm{Ext}(\Lambda,\mathbb{Z}))$ of classifying all group extensions of $\Lambda$ by $\mathbb{Z}$ up to base-free isomorphism, when $\Lambda =\mathbb{Z}[1/p]^{d}$ for prime numbers $p$ and $ d\geq 1$.


title: KS-pies: Kohn-Sham Inversion Toolkit

abstract: A Kohn-Sham (KS) inversion determines a KS potential and orbitals corresponding to a given electron density, a procedure that has applications in developing and evaluating functionals used in density functional theory. Despite the utility of KS inversions, application of these methods among the research community is disproportionately small. We implement the KS inversion methods of Zhao-Morrison-Parr and Wu-Yang in a framework that simplifies analysis and conversion of the resulting potential in real-space. Fully documented Python scripts integrate with PySCF, a popular electronic structure prediction software, and Fortran alternatives are provided for computational hot spots.


title: An adaptive proximal point algorithm framework and application to large-scale optimization

abstract: We investigate the proximal point algorithm (PPA) and its inexact extensions under an error bound condition, which guarantees a global linear convergence if the proximal regularization parameter is larger than the error bound condition parameter. We propose an adaptive generalized proximal point algorithm (AGPPA), which adaptively updates the proximal regularization parameters based on some implementable criteria. We show that AGPPA achieves linear convergence without any knowledge of the error bound condition parameter, and the rate only differs from the optimal one by a logarithm term. We apply AGPPA on convex minimization problem and analyze the iteration complexity bound of the resulting algorithm. Our framework and the complexity results apply to arbitrary linearly convergent inner solver and allows a hybrid with any locally fast convergent method. We illustrate the performance of AGPPA by applying it to solve large-scale linear programming (LP) problem. The resulting complexity bound has a weaker dependence on the Hoffman constant and scales with the dimension better than linearized ADMM. In numerical experiments, our algorithm demonstrates improved performance in obtaining solution of medium accuracy on large-scale LP problem.


title: Anisotropic ultrafast spin/valley dynamics in WTe2 films

abstract: WTe2 Weyl semimetal hosts the natural broken inversion symmetry and strong spin orbit coupling, making it promising for exotic spin/valley dynamics within a picosecond timescale. Here, we unveil an anisotropic ultrafast spin/valley dynamics in centimeter-scale, single-crystalline Td-WTe2 films using a femtosecond pump-probe technique at room temperature. We observe a transient (~0.8 ps) intra-valley transition and a subsequent polarization duration (~5 ps) during the whole spin/valley relaxation process. Furthermore, the relaxation exhibits the remarkable anisotropy of approximately six-fold and two-fold symmetries due to the intrinsic anisotropy along the crystalline orientation and the extrinsic matrix element effect, respectively. Our results offer a prospect for the ultrafast manipulation of spin/valleytronics in topological quantum materials for dissipationless high-speed spin/valleytronic devices.


title: Numerical solutions to Giovannini's parton branching equation up to TeV energies at the LHC

abstract: Giovannini's parton branching equation is integrated numerically using the 4th-order Runge-Kutta method. Using a simple hadronisation model, a charged-hadron multiplicity distribution is obtained. This model is then fitted to various experimental data up to the TeV scale to study how the Giovannini parameters vary with collision energy and type. The model is able to describe hadronic collisions up to the TeV scale and reveals the emergence of gluonic activity as the centre-of-mass energy increases. A prediction is made for $\sqrt{s}$ = 14 TeV.


title: Formation of magnetized spatial structures in the Beta Lyrae system II. Reflection of magnetically controlled structures in the visible spectrum

abstract: This article proposes a picture of magnetized accretion structures formed during the mass transfer in the Beta Lyrae system. It is shown that the structure of the gaseous flows between the donor and the gainer is due to the spatial configuration of the donor magnetic field. Its dipole axis is deviated substantially from the line joining the centers of the components and is inclined to the orbital plane of the binary system; the center of the magnetic dipole is displaced from the donor center toward the gainer. The surface around the donor magnetic pole, which is close to the gainer, is a region of an additional matter loss from the donor surface. The effective collision of the magnetized plasma with the accretion disk is enhanced by the fast counter-rotation of this disk, especially in the secondary quadrature phases, in which the high-temperature medium and the system of formed accretion flows are observed. This concept is demonstrated, primarily, in the obvious correlations between the phase variability of the donor magnetic field and the corresponding variability of the dynamic and energy characteristics of the various complex lines. This refers to the behavior of the radial velocity curves of the emission-absorption lines formed in the gaseous structures of type H$\alpha$, HeI $\lambda$ 7065, or the variability of their equivalent width and intensity, and the variability of conventional absorption lines of the donor atmosphere. This is true for the phase variability of the absolute flux in the H$\alpha$ emission line and the fast varying of the continuum in the H$\alpha$ region as certain parameters, which reflect the phase variability of the donor magnetic field. This approach made it possible to determine the phase boundaries of the location of the magnetic polar region on the donor surface above which the matter outflows are formed.


title: A lower bound on HMOLS with equal sized holes

abstract: It is known that $N(n)$, the maximum number of mutually orthogonal latin squares of order $n$, satisfies the lower bound $N(n) \ge n^{1/14.8}$ for large $n$. For $h\ge 2$, relatively little is known about the quantity $N(h^n)$, which denotes the maximum number of `HMOLS' or mutually orthogonal latin squares having a common equipartition into $n$ holes of a fixed size $h$. We generalize a difference matrix method that had been used previously for explicit constructions of HMOLS. An estimate of R.M. Wilson on higher cyclotomic numbers guarantees our construction succeeds in suitably large finite fields. Feeding this into a generalized product construction, we are able to establish the lower bound $N(h^n) \ge (\log n)^{1/\delta}$ for any $\delta>2$ and all $n > n_0(h,\delta)$.


title: Measurement of $\gamma$+jet and $\pi^{0}$+jet in central Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 200 GeV with the STAR experiment

abstract: We present the semi-inclusive measurement of charged jets recoiling from direct-photon and $\pi^{0}$ triggers in central Au+Au collisions at $\sqrt{s_{\mathrm {NN}}}$ = 200 GeV, using a dataset with integrated luminosity 13 $\mathrm {nb^{-1}}$ recorded by the STAR experiment in 2014. The photon and $\pi^{0}$ triggers are selected within transverse energy ($\rm E_\mathrm{T}^\mathrm{trig}$) between 9 GeV and 20 GeV. Charged jets are reconstructed with the anti-$k_{\mathrm{T}}$ algorithm with resolution parameters R = 0.2 and 0.5. A Mixed-Event technique developed previously by STAR is used to correct the recoil jet yield for uncorrelated background, enabling recoil jet measurements over a broad $\rm p_{T,jet}$ range. We report fully corrected charged-jet yields recoiling from direct-photon and $\pi^{0}$ triggers for the above two jet radii and also discuss the jet R dependence of in-medium parton energy loss at the top RHIC energy.


title: Joint Visual and Wireless Signal Feature based Approach for High-Precision Indoor Localization

abstract: The existing localization systems for indoor applications basically rely on wireless signal. With the massive deployment of low-cost cameras, the visual image based localization become attractive as well. However, in the existing literature, the hybrid visual and wireless approaches simply combine the above schemes in a straight forward manner, and fail to explore the interactions between them. In this paper, we propose a joint visual and wireless signal feature based approach for high-precision indoor localization system. In this joint scheme, WiFi signals are utilized to compute the coarse area with likelihood probability and visual images are used to fine-tune the localization result. Based on the numerical results, we show that the proposed scheme can achieve 0.62m localization accuracy with near real-time running time.


title: Facial movement synergies and Action Unit detection from distal wearable Electromyography and Computer Vision

abstract: Distal facial Electromyography (EMG) can be used to detect smiles and frowns with reasonable accuracy. It capitalizes on volume conduction to detect relevant muscle activity, even when the electrodes are not placed directly on the source muscle. The main advantage of this method is to prevent occlusion and obstruction of the facial expression production, whilst allowing EMG measurements. However, measuring EMG distally entails that the exact source of the facial movement is unknown. We propose a novel method to estimate specific Facial Action Units (AUs) from distal facial EMG and Computer Vision (CV). This method is based on Independent Component Analysis (ICA), Non-Negative Matrix Factorization (NNMF), and sorting of the resulting components to determine which is the most likely to correspond to each CV-labeled action unit (AU). Performance on the detection of AU06 (Orbicularis Oculi) and AU12 (Zygomaticus Major) was estimated by calculating the agreement with Human Coders. The results of our proposed algorithm showed an accuracy of 81% and a Cohen's Kappa of 0.49 for AU6; and accuracy of 82% and a Cohen's Kappa of 0.53 for AU12. This demonstrates the potential of distal EMG to detect individual facial movements. Using this multimodal method, several AU synergies were identified. We quantified the co-occurrence and timing of AU6 and AU12 in posed and spontaneous smiles using the human-coded labels, and for comparison, using the continuous CV-labels. The co-occurrence analysis was also performed on the EMG-based labels to uncover the relationship between muscle synergies and the kinematics of visible facial movement.


title: Families with no perfect matchings

abstract: We consider families of $k$-subsets of $\{1, \dots, n\}$, where $n$ is a multiple of $k$, which have no perfect matching. An equivalent condition for a family $\mathcal{F}$ to have no perfect matching is for there to be a blocking set, which is a set of $b$ elements of $\{1, \dots, n\}$ that cannot be covered by $b$ disjoint sets in $\mathcal{F}$. We are specifically interested in the largest possible size of a family $\mathcal{F}$ with no perfect matching and no blocking set of size less than $b$. Frankl resolved the case of families with no singleton blocking set (in other words, the $b=2$ case) for sufficiently large $n$ and conjectured an optimal construction for general $b$. Though Frankl's construction fails to be optimal for $k = 2, 3$, we show that the construction is optimal whenever $k \ge 100$ and $n$ is sufficiently large.


title: Explainable Machine Learning for Materials Discovery: Predicting the Potentially Formable Nd-Fe-B Crystal Structures and Extracting Structure-Stability Relationship

abstract: New Nd-Fe-B crystal structures can be formed via the elemental substitution of LATX host structures, including lanthanides LA, transition metals T, and light elements X as B, C, N, and O. The 5967 samples of ternary LATX materials that are collected are then used as the host structures. For each host crystal structure, a substituted crystal structure is created by substituting all lanthanide sites with Nd, all transition metal sites with Fe, and all light element sites with B. High throughput first-principles calculations are applied to evaluate the phase stability of the newly created crystal structures, and 20 of them are found to be potentially formable. A data driven approach based on supervised and unsupervised learning techniques is applied to estimate the stability and analyze the structure stability relationship of the newly created NdFeB crystal structures. For predicting the stability for the newly created NdFeB structures, three supervised learning models, kernel ridge regression, logistic classification, and decision tree model, are learned from the LATX host crystal structures; the models achieve the maximum accuracy and recall scores of 70.4 and 68.7 percent, respectively. On the other hand, our proposed unsupervised learning model based on the integration of descriptor-relevance analysis and a Gaussian mixture model achieves accuracy and recall score of 72.9 and 82.1 percent, respectively, which are significantly better than those of the supervised models. While capturing and interpreting the structure stability relationship of the NdFeB crystal structures, the unsupervised learning model indicates that the average atomic coordination number and coordination number of the Fe sites are the most important factors in determining the phase stability of the new substituted NdFeB crystal structures.


title: A Regional Bolus Tracking and Real-time B$_1$ Calibration Method for Hyperpolarized $^{13}$C MRI

abstract: Purpose: Acquisition timing and B$_1$ calibration are two key factors that affect the quality and accuracy of hyperpolarized $^{13}$C MRI. The goal of this project was to develop a new approach using regional bolus tracking to trigger Bloch-Siegert B$_1$ mapping and real-time B$_1$ calibration based on regional B$_1$ measurements, followed by dynamic imaging of hyperpolarized $^{13}C$ metabolites in vivo. Methods: The proposed approach was implemented on a system which allows real-time data processing and real-time control on the sequence. Real-time center frequency calibration upon the bolus arrival was also added. The feasibility of applying the proposed framework for in vivo hyperpolarized $^{13}$C imaging was tested on healthy rats, tumor-bearing mice and a healthy volunteer on a clinical 3T scanner following hyperpolarized [1-$^{13}$C]pyruvate injection. Multichannel receive coils were used in the human study. Results: Automatic acquisition timing based on either regional bolus peak or bolus arrival was achieved with the proposed framework. Reduced blurring artifacts in real-time reconstructed images were observed with real-time center frequency calibration. Real-time computed B$_1$ scaling factors agreed with real-time acquired B$_1$ maps. Flip angle correction using B$_1$ maps results in a more consistent quantification of metabolic activity (i.e, pyruvate-to-lactate conversion, k$_{PL}$). Experiment recordings are provided to demonstrate the real-time actions during the experiment. Conclusion: The proposed method was successfully demonstrated on animals and a human volunteer, and is anticipated to improve the efficient use of the hyperpolarized signal as well as the accuracy and robustness of hyperpolarized $^{13}$C imaging.


title: A universal frequency engineering tool for microcavity nonlinear optics: multiple selective mode splitting of whispering-gallery resonances

abstract: Whispering-gallery microcavities have been used to realize a variety of efficient parametric nonlinear optical processes through the enhanced light-matter interaction brought about by supporting multiple high quality factor and small modal volume resonances. Critical to such studies is the ability to control the relative frequencies of the cavity modes, so that frequency matching is achieved to satisfy energy conservation. Typically this is done by tailoring the resonator cross-section. Doing so modifies the frequencies of all of the cavity modes, that is, the global dispersion profile, which may be undesired, for example, in introducing competing nonlinear processes.Here, we demonstrate a frequency engineering tool, termed multiple selective mode splitting (MSMS), that is independent of the global dispersion and instead allows targeted and independent control of the frequencies of multiple cavity modes. In particular, we show controllable frequency shifts up to 0.8 nm, independent control of the splitting of up to five cavity modes with optical quality factors $\gtrsim 10^5$, and strongly suppressed frequency shifts for untargeted modes. The MSMS technique can be broadly applied to a wide variety of nonlinear optical processes across different material platforms, and can be used to both selectively enhance processes of interestand suppress competing unwanted processes.