title: Causal Sufficiency and Actual Causation

abstract: Pearl opened the door to formally defining actual causation using causal models. His approach rests on two strategies: first, capturing the widespread intuition that X=x causes Y=y iff X=x is a Necessary Element of a Sufficient Set for Y=y, and second, showing that his definition gives intuitive answers on a wide set of problem cases. This inspired dozens of variations of his definition of actual causation, the most prominent of which are due to Halpern & Pearl. Yet all of them ignore Pearl's first strategy, and the second strategy taken by itself is unable to deliver a consensus. This paper offers a way out by going back to the first strategy: it offers six formal definitions of causal sufficiency and two interpretations of necessity. Combining the two gives twelve new definitions of actual causation. Several interesting results about these definitions and their relation to the various Halpern & Pearl definitions are presented. Afterwards the second strategy is evaluated as well. In order to maximize neutrality, the paper relies mostly on the examples and intuitions of Halpern & Pearl. One definition comes out as being superior to all others, and is therefore suggested as a new definition of actual causation.


title: Eigenstates of Quasi-Keplerian Self-Gravitating Particle Discs

abstract: Although quasi-Keplerian discs are among the most common astrophysical structures, computation of secular angular momentum transport within them routinely presents a considerable practical challenge. In this work, we investigate the secular small-inclination dynamics of a razor-thin particle disc as the continuum limit of a discrete Lagrange-Laplace secular perturbative theory and explore the analogy between the ensuing secular evolution -- including non-local couplings of self-gravitating discs -- and quantum mechanics. We find the 'quantum' Hamiltonian that describes the time evolution of the system and demonstrate the existence of a conserved inner product. The lowest-frequency normal modes are numerically approximated by performing a Wick rotation on the equations of motion. These modes are used to quantify the accuracy of a much simpler local-coupling model, revealing that it predicts the shape of the normal modes to a high degree of accuracy, especially in narrow annuli, even though it fails to predict their eigenfrequencies.


title: Decoupling of Asymmetric Dark Matter During an Early Matter Dominated Era

abstract: In models of Asymmetric Dark Matter (ADM) the relic density is set by a particle asymmetry in an analogous manner to the baryons. Here we explore the scenario in which ADM decouples from the Standard Model thermal bath during an early period of matter domination. We first present a model independent analysis for a generic ADM candidate with s-wave annihilation cross section with fairly general assumptions regarding the origin of the early matter dominated period. We contrast our results to those from conventional ADM models which assume radiation domination during decoupling. Subsequently, we examine an explicit example of this scenario in the context of an elegant SO(10) implementation of ADM in which the matter dominated era is due to a long lived heavy right-handed neutrino. In the concluding remarks we discuss the prospects for superheavy ADM in this setting.


title: Transient heat transfer of superfluid $^4$He in nonhomogeneous geometries -- Part I: Second sound, rarefaction, and thermal layer

abstract: Transient heat transfer in superfluid $^4$He (He II) is a complex process that involves the interplay of the unique counterflow heat-transfer mode, the emission of second-sound waves, and the creation of quantized vortices. Many past researches focused on homogeneous heat transfer of He II in a uniform channel driven by a planar heater. In this paper, we report our systematic study of He II transient heat transfer in nonhomogeneous geometries that are pertinent to emergent applications. By solving the He II two-fluid equation of motion coupled with the Vinen's equation for vortex-density evolution, we examine and compare the characteristics of transient heat transfer from planar, cylindrical, and spherical heaters in He II. Our results show that as the heater turns on, an outgoing second-sound pulse emerges, in which the vortex density grows rapidly. These vortices attenuate the second sound and result in a heated He II layer in front of the heater, i.e., the thermal layer. In the planar case where the vortices are created throughout the space, the second-sound pulse is continuously attenuated, leading to a strong thermal layer that diffusively spreads following the heat pulse. On the contrary, in the cylindrical and the spherical heater cases, vortices are created mainly in a thin thermal layer near the heater surface. As the heat pulse ends, a rarefaction tail develops following the second-sound pulse, in which the temperature drops. This rarefaction tail can promptly suppress the thermal layer and take away all the thermal energy deposited in it. The effects of the heater size, heat flux, pulse duration, and temperature on the thermal-layer dynamics are discussed. We also show how the peak heat flux for the onset of boiling in He II can be studied in our model.


title: Real-Time Optimal Trajectory Planning for Autonomous Vehicles and Lap Time Simulation Using Machine Learning

abstract: Widespread development of driverless vehicles has led to the formation of autonomous racing, where technological development is accelerated by the high speeds and competitive environment of motorsport. A particular challenge for an autonomous vehicle is that of identifying a target trajectory - or, in the case of a competition vehicle, the racing line. Many existing approaches to finding the racing line are either not time-optimal solutions, or are computationally expensive - rendering them unsuitable for real-time application using on-board processing hardware. This study describes a machine learning approach to generating an accurate prediction of the racing line in real-time on desktop processing hardware. The proposed algorithm is a feed-forward neural network, trained using a dataset comprising racing lines for a large number of circuits calculated via traditional optimal control lap time simulation. The network predicts the racing line with a mean absolute error of +/-0.27m, and just +/-0.11m at corner apex - comparable to human drivers, and autonomous vehicle control subsystems. The approach generates predictions within 33ms, making it over 9,000 times faster than traditional methods of finding the optimal trajectory. Results suggest that for certain applications data-driven approaches to find near-optimal racing lines may be favourable to traditional computational methods.


title: Dust Traps and the Formation of Cavities in Transition Discs: A millimetre to sub-millimetre comparison survey

abstract: The origin of the inner dust cavities observed in transition discs remains unknown. The segregation of dust and size of the cavity is expected to vary depending on which clearing mechanism dominates grain evolution. We present the results from the Discs Down Under program, an 8.8 mm continuum Australia Telescope Compact Array (ATCA) survey targeting 15 transition discs with large (> 20 au) cavities, and compare the resulting dust emission to Atacama Large millimetre/sub-millimetre Array (ALMA) observations. Our ATCA observations resolve the inner cavity for 8 of the 14 detected discs. We fit the visibilities and reconstruct 1D radial brightness models for 10 sources with a S/N > 5sigma. We find that, for sources with a resolved cavity in both wavebands, the 8.8 mm and sub-mm brightness distributions peak at the same radius from the star. We suggest that a similar cavity size for 8.8 mm and sub-mm dust grains is due to a dust trap induced by the presence of a companion.


title: Nearly hyperuniform, nonhyperuniform, and antihyperuniform density fluctuations in two-dimensional transition metal dichalcogenides with defects

abstract: Hyperuniform many-body systems in $d$-dimensional Euclidean space are characterized by completely suppressed (normalized) infinite-wavelength density fluctuations, and appear to be endowed with novel exotic physical properties. In this work, we analyze the effects of localized defects on the density fluctuations across length scales and on the hyperuniformity property of experimental samples of two-dimensional transition metal dichalcogenides. In particular, we extract atomic coordinates from time series annular dark field-scanning transmission electron microscopy imaging data of 2D tungsten chalcogenides with the 2H structure showing continuous development and evolution of electron-beam induced defects, and construct the corresponding chemical-bonding informed coordination networks between the atoms. We then compute a variety of pair statistics and bond-orientational statistics to characterize the samples. At low defect concentrations, the corresponding materials are nearly hyperuniform. As more defects are introduced, the (approximate) hyperuniformity of the materials is gradually destroyed, and the system becomes non-hyperuniform even when the material still contains a significant amount of crystalline regions. At high defect concentrations, the structures become antihyperuniform with diverging (normalized) large-scale density fluctuations. Overall, the defected materials possess varying degrees of orientation order, and there is apparently no intermediate hexatic phase emerging. We also construct a minimalist structural model and demonstrate that the experimental samples can be essentially viewed as perturbed honeycomb crystals with small correlated displacements and double chalcogen vacancies.


title: System Intelligence for UAV-Based Mission Critical with Challenging 5G/B5G Connectivity

abstract: Unmanned aerial vehicles (UAVs) and communication systems are fundamental elements in Mission Critical services, such as search and rescue. In this article, we introduce an architecture for managing and orchestrating 5G and beyond networks that operate over a heterogeneous infrastructure with UAVs' aid. UAVs are used for collecting and processing data, as well as improving communications. The proposed System Intelligence (SI) architecture was designed to comply with recent standardization works, especially the ETSI Experiential Networked Intelligence specifications. Another contribution of this article is an evaluation using a testbed based on a virtualized non-standalone 5G core and a 4G Radio Access Network (RAN) implemented with open-source software. The experimental results indicate, for instance, that SI can substantially improve the latency of UAV-based services by splitting deep neural networks between UAV and edge or cloud equipment. Other experiments explore the slicing of RAN resources and efficient placement of virtual network functions to assess the benefits of incorporating intelligence in UAV-based mission-critical services.


title: Spitzer IRAC observations of JWST calibration stars

abstract: We present infrared photometry of all 36 potential JWST calibrators for which there is archival Spitzer IRAC data. This photometry can then be used to inform stellar models necessary to provide absolute calibration for all JWST instruments. We describe in detail the steps necessary to measure IRAC photometry from archive retrieval to photometric corrections. To validate our photometry we examine the distribution of uncertainties from all detections in all four IRAC channels as well as compare the photometry and its uncertainties to those from models, ALLWISE, and the literature. 75% of our detections have standard deviations per star of all observations within each channel of less than three percent. The median standard deviations are 1.2, 1.3, 1.1, and 1.9% in [3.6] - [8.0] respectively. We find less than 8% standard deviations in differences of our photometry with ALLWISE, and excellent agreement with literature values (less than 3% difference) lending credence to our measured fluxes. JWST is poised to do ground-breaking science, and accurate calibration and cross-calibration with other missions will be part of the underpinnings of that science.


title: One Label, One Billion Faces: Usage and Consistency of Racial Categories in Computer Vision

abstract: Computer vision is widely deployed, has highly visible, society altering applications, and documented problems with bias and representation. Datasets are critical for benchmarking progress in fair computer vision, and often employ broad racial categories as population groups for measuring group fairness. Similarly, diversity is often measured in computer vision datasets by ascribing and counting categorical race labels. However, racial categories are ill-defined, unstable temporally and geographically, and have a problematic history of scientific use. Although the racial categories used across datasets are superficially similar, the complexity of human race perception suggests the racial system encoded by one dataset may be substantially inconsistent with another. Using the insight that a classifier can learn the racial system encoded by a dataset, we conduct an empirical study of computer vision datasets supplying categorical race labels for face images to determine the cross-dataset consistency and generalization of racial categories. We find that each dataset encodes a substantially unique racial system, despite nominally equivalent racial categories, and some racial categories are systemically less consistent than others across datasets. We find evidence that racial categories encode stereotypes, and exclude ethnic groups from categories on the basis of nonconformity to stereotypes. Representing a billion humans under one racial category may obscure disparities and create new ones by encoding stereotypes of racial systems. The difficulty of adequately converting the abstract concept of race into a tool for measuring fairness underscores the need for a method more flexible and culturally aware than racial categories.


title: Hamiltonicity of graphs perturbed by a random geometric graph

abstract: We study Hamiltonicity in graphs obtained as the union of a deterministic $n$-vertex graph $H$ with linear degrees and a $d$-dimensional random geometric graph $G^d(n,r)$, for any $d\geq1$. We obtain an asymptotically optimal bound on the minimum $r$ for which a.a.s. $H\cup G^d(n,r)$ is Hamiltonian. Our proof provides a linear time algorithm to find a Hamilton cycle in such graphs.


title: Query Complexity of Least Absolute Deviation Regression via Robust Uniform Convergence

abstract: Consider a regression problem where the learner is given a large collection of $d$-dimensional data points, but can only query a small subset of the real-valued labels. How many queries are needed to obtain a $1+\epsilon$ relative error approximation of the optimum? While this problem has been extensively studied for least squares regression, little is known for other losses. An important example is least absolute deviation regression ($\ell_1$ regression) which enjoys superior robustness to outliers compared to least squares. We develop a new framework for analyzing importance sampling methods in regression problems, which enables us to show that the query complexity of least absolute deviation regression is $\Theta(d/\epsilon^2)$ up to logarithmic factors. We further extend our techniques to show the first bounds on the query complexity for any $\ell_p$ loss with $p\in(1,2)$. As a key novelty in our analysis, we introduce the notion of robust uniform convergence, which is a new approximation guarantee for the empirical loss. While it is inspired by uniform convergence in statistical learning, our approach additionally incorporates a correction term to avoid unnecessary variance due to outliers. This can be viewed as a new connection between statistical learning theory and variance reduction techniques in stochastic optimization, which should be of independent interest.


title: Volumes spanned by $k$ point configurations in $\mathbb{R}^d$

abstract: Given a $k$-point configuration $x\in (\mathbb{R}^d)^k$, we consider the $\binom{k}{d}$-vector of volumes determined by choosing any $d$ points of $x$. We prove that a compact set $E\subset \R^d$ determines a positive measure of such volume types if the Hausdorff dimension of $E$ is greater than $d-\frac{d-1}{2k-d}$. This generalizes results of Greenleaf, Iosevich, and Mourgoglou, Greenleaf, Iosevich, and Taylor, and the second listed author.


title: The SPIRou wavelength calibration for precise radial velocities in the near infrared

abstract: SPIRou is a near-infrared (nIR) spectropolarimeter at the CFHT, covering the YJHK nIR spectral bands ($980-2350\,\mathrm{nm}$). We describe the development and current status of the SPIRou wavelength calibration in order to obtain precise radial velocities (RVs) in the nIR. We make use of a UNe hollow-cathode lamp and a Fabry-P\'erot \'etalon to calibrate the pixel-wavelength correspondence for SPIRou. Different methods are developed for identifying the hollow-cathode lines, for calibrating the wavelength dependence of the Fabry-P\'erot cavity width, and for combining the two calibrators. The hollow-cathode spectra alone do not provide a sufficiently accurate wavelength solution to meet the design requirements of an internal error of $\mathrm{<0.45\,m\,s^{-1}}$, for an overall RV precision of $\mathrm{1\,m\,s^{-1}}$. However, the combination with the Fabry-P\'erot spectra allows for significant improvements, leading to an internal error of $\mathrm{\sim 0.15\,m\,s^{-1}}$. We examine the inter-night stability, intra-night stability, and impact on the stellar RVs of the wavelength solution.


title: Spin-flip pair-density functional theory: A practical approach to treat static and dynamical correlations in large molecules

abstract: We present a practical approach to treat static and dynamical correlation accurately in large multi-configurational systems. The static correlation is accounted for using the spin-flip approach which is well known for capturing static correlation accurately at low-computational expense. Unlike previous approaches to add dynamical correlation to spin-flip models which use perturbation theory or coupled-cluster theory, we explore the ability to use the on-top pair-density functional theory approaches recently developed by Gagliardi and co-workers (JCTC, 10, 3669, 2014). External relaxations are carried out in the spin-flip calculations though a restricted active space framework for which a truncation scheme for the orbitals used in the external excitation is presented. The performance of the approach is demonstrated by computing energy gaps between ground and excited states for diradicals, triradicals and linear polyacene chains ranging from naphthalene to dodecacene. Accurate results are obtained using the new approach for these challenging open-shell molecular systems.


title: Effects of Number of Filters of Convolutional Layers on Speech Recognition Model Accuracy

abstract: Inspired by the progress of the End-to-End approach [1], this paper systematically studies the effects of Number of Filters of convolutional layers on the model prediction accuracy of CNN+RNN (Convolutional Neural Networks adding to Recurrent Neural Networks) for ASR Models (Automatic Speech Recognition). Experimental results show that only when the CNN Number of Filters exceeds a certain threshold value is adding CNN to RNN able to improve the performance of the CNN+RNN speech recognition model, otherwise some parameter ranges of CNN can render it useless to add the CNN to the RNN model. Our results show a strong dependency of word accuracy on the Number of Filters of convolutional layers. Based on the experimental results, the paper suggests a possible hypothesis of Sound-2-Vector Embedding (Convolutional Embedding) to explain the above observations. Based on this Embedding hypothesis and the optimization of parameters, the paper develops an End-to-End speech recognition system which has a high word accuracy but also has a light model-weight. The developed LVCSR (Large Vocabulary Continuous Speech Recognition) model has achieved quite a high word accuracy of 90.2% only by its Acoustic Model alone, without any assistance from intermediate phonetic representation and any Language Model. Its acoustic model contains only 4.4 million weight parameters, compared to the 35~68 million acoustic-model weight parameters in DeepSpeech2 [2] (one of the top state-of-the-art LVCSR models) which can achieve a word accuracy of 91.5%. The light-weighted model is good for improving the transcribing computing efficiency and also useful for mobile devices, Driverless Vehicles, etc. Our model weight is reduced to ~10% the size of DeepSpeech2, but our model accuracy remains close to that of DeepSpeech2. If combined with a Language Model, our LVCSR system is able to achieve 91.5% word accuracy.


title: Two- and three-nucleon contact interactions and ground-state energies of light- and medium-mass nuclei

abstract: Classes of two-nucleon ($2N$) contact interactions are developed in configuration space at leading order (LO), next-to-leading order (NLO), and next-to-next-to-next-to-leading order (N3LO) by fitting the experimental singlet $np$ scattering length and deuteron binding energy at LO, and $np$ and $pp$ scattering data in the laboratory-energy range 0--15 MeV at NLO and 0--25 MeV at N3LO. These interactions are regularized by including two Gaussian cutoffs, one for $T\,$=$\,0$ and the other for $T\,$=$\,1$ channels. The cutoffs are taken to vary in the ranges $R_0\,$=$(1.5$--2.3) fm and $R_1\,$=$(1.5$--3.0) fm. The 780 (1,100) data points up to 15 (25) MeV energy, primarily differential cross sections, are fitted by the NLO (N3LO) models with a $\chi^2$/datum about 1.7 or less (well below 1.5), when harder cutoff values are adopted. As a first application, we report results for the binding energies of nuclei with mass numbers $A\,$=$\,3$--6 and 16 obtained with selected LO and NLO $2N$ models both by themselves as well as in combination with a LO three-nucleon ($3N$) contact interaction. The latter is characterized by a single low-energy constant that is fixed to reproduce the experimental $^3$H binding energy. The inclusion of the $3N$ interaction largely removes the sensitivity to cutoff variations in the few-nucleon systems and leads to predictions for the $^3$He and $^4$He binding energies that cluster around 7.8 MeV and 30 MeV, respectively. However, in $^{16}$O this cutoff sensitivity remains rather strong. Finally, predictions at LO only are also reported for medium-mass nuclei with $A\,$=$\,40$, 48, and 90.


title: Slowly rotating black holes modeled by Solv geometries

abstract: We present a slowly rotating generalization of a black hole modeled by a Solv horizon geometry, in five dimensional General Relativity. A separable ansatz compatible with Einstein equations is proposed, which is integrated in terms of hypergeometric and exponential functions. The solution highly simplifies after requiring to maintain the same asymptotic behaviour as the non-rotating black hole and to have finite, non-vanishing charges. We compute the charges of the solution and also comment on its importance for studying the phase space of the theory in this sector.


title: Modeling Financial Products and their Supply Chains

abstract: The objective of this paper is to explore how financial big data and machine learning methods can be applied to model and understand financial products. We focus on residential mortgage backed securities, resMBS, which were at the heart of the 2008 US financial crisis. These securities are contained within a prospectus and have a complex waterfall payoff structure. Multiple financial institutions form a supply chain to create prospectuses. To model this supply chain, we use unsupervised probabilistic methods, particularly dynamic topics models (DTM), to extract a set of features (topics) reflecting community formation and temporal evolution along the chain. We then provide insight into the performance of the resMBS securities and the impact of the supply chain through a series of increasingly comprehensive models. First, models at the security level directly identify salient features of resMBS securities that impact their performance. We then extend the model to include prospectus level features and demonstrate that the composition of the prospectus is significant. Our model also shows that communities along the supply chain that are associated with the generation of the prospectuses and securities have an impact on performance. We are the first to show that toxic communities that are closely linked to financial institutions that played a key role in the subprime crisis can increase the risk of failure of resMBS securities.


title: Function Delivery Network: Extending Serverless Computing for Heterogeneous Platforms

abstract: Serverless computing has rapidly grown following the launch of Amazon's Lambda platform. Function-as-a-Service (FaaS) a key enabler of serverless computing allows an application to be decomposed into simple, standalone functions that are executed on a FaaS platform. The FaaS platform is responsible for deploying and facilitating resources to the functions. Several of today's cloud applications spread over heterogeneous connected computing resources and are highly dynamic in their structure and resource requirements. However, FaaS platforms are limited to homogeneous clusters and homogeneous functions and do not account for the data access behavior of functions before scheduling. We introduce an extension of FaaS to heterogeneous clusters and to support heterogeneous functions through a network of distributed heterogeneous target platforms called Function Delivery Network (FDN). A target platform is a combination of a cluster of homogeneous nodes and a FaaS platform on top of it. FDN provides Function-Delivery-as-a-Service (FDaaS), delivering the function to the right target platform. We showcase the opportunities such as varied target platform's characteristics, possibility of collaborative execution between multiple target platforms, and localization of data that the FDN offers in fulfilling two objectives: Service Level Objective (SLO) requirements and energy efficiency when scheduling functions by evaluating over five distributed target platforms using the FDNInspector, a tool developed by us for benchmarking distributed target platforms. Scheduling functions on an edge target platform in our evaluation reduced the overall energy consumption by 17x without violating the SLO requirements in comparison to scheduling on a high-end target platform.


title: Low Radio Frequency Observations from the Moon Enabled by NASA Landed Payload Missions

abstract: A new era of exploration of the low radio frequency Universe from the Moon will soon be underway with landed payload missions facilitated by NASA's Commercial Lunar Payload Services (CLPS) program. CLPS landers are scheduled to deliver two radio science experiments, ROLSES to the nearside and LuSEE to the farside, beginning in 2021. These instruments would be pathfinders for a 10-km diameter interferometric array, FARSIDE, composed of 128 pairs of dipole antennas proposed to be delivered to the lunar surface later in the decade. ROLSES and LuSEE, operating at frequencies from 100 kHz to a few tens of MHz, will investigate the plasma environment above the lunar surface and measure the fidelity of radio spectra on the surface. Both use electrically-short, spiral-tube deployable antennas and radio spectrometers based upon previous flight models. ROLSES will measure the photoelectron sheath density to better understand the charging of the lunar surface via photoionization and impacts from the solar wind, charged dust, and current anthropogenic radio frequency interference. LuSEE will measure the local magnetic field and exo-ionospheric density, interplanetary radio bursts, Jovian and terrestrial natural radio emission, and the galactic synchrotron spectrum. FARSIDE, and its precursor risk-reduction six antenna-node array PRIME, would be the first radio interferometers on the Moon. FARSIDE would break new ground by imaging radio emission from Coronal Mass Ejections (CME) beyond 2 solar radii, monitor auroral radiation from the B-fields of Uranus and Neptune (not observed since Voyager), and detect radio emission from stellar CMEs and the magnetic fields of nearby potentially habitable exoplanets.


title: The Enigma of Complexity

abstract: In this paper we examine the concept of complexity as it applies to generative art and design. Complexity has many different, discipline specific definitions, such as complexity in physical systems (entropy), algorithmic measures of information complexity and the field of "complex systems". We apply a series of different complexity measures to three different generative art datasets and look at the correlations between complexity and individual aesthetic judgement by the artist (in the case of two datasets) or the physically measured complexity of 3D forms. Our results show that the degree of correlation is different for each set and measure, indicating that there is no overall "better" measure. However, specific measures do perform well on individual datasets, indicating that careful choice can increase the value of using such measures. We conclude by discussing the value of direct measures in generative and evolutionary art, reinforcing recent findings from neuroimaging and psychology which suggest human aesthetic judgement is informed by many extrinsic factors beyond the measurable properties of the object being judged.


title: Entanglement Cost for Steering Assemblages

abstract: In this paper we pose the question of how one can quantify the amount of entanglement necessary to create a given steering assemblage, when the original state and measurements are unknown. To do this, we extend the concepts of entanglement cost and entanglement of formation to steering assemblages, and show how to upper bound them with semidefinite programming. We prove that the entanglement of formation for assemblages is not generally continuous and not a flat roof extension; and use numerical analysis to illustrate these properties. Finally, we discuss the consequences of these results for assemblage-to-assemblage conversion and directions for further research.


title: Zero-state Coupled Markov Switching Count Models for Spatio-temporal Infectious Disease Spread

abstract: Spatio-temporal counts of infectious disease cases often contain an excess of zeros. With existing zero inflated count models applied to such data it is difficult to quantify space-time heterogeneity in the effects of disease spread between areas. Also, existing methods do not allow for separate dynamics to affect the reemergence and persistence of the disease. As an alternative, we develop a new zero-state coupled Markov switching negative binomial model, under which the disease switches between periods of presence and absence in each area through a series of partially hidden nonhomogeneous Markov chains coupled between neighboring locations. When the disease is present, an autoregressive negative binomial model generates the cases with a possible 0 representing the disease being undetected. Bayesian inference and prediction is illustrated using spatio-temporal counts of dengue fever cases in Rio de Janeiro, Brazil.


title: Self-Supervised Claim Identification for Automated Fact Checking

abstract: We propose a novel, attention-based self-supervised approach to identify "claim-worthy" sentences in a fake news article, an important first step in automated fact-checking. We leverage "aboutness" of headline and content using attention mechanism for this task. The identified claims can be used for downstream task of claim verification for which we are releasing a benchmark dataset of manually selected compelling articles with veracity labels and associated evidence. This work goes beyond stylistic analysis to identifying content that influences reader belief. Experiments with three datasets show the strength of our model. Data and code available at https://github.com/architapathak/Self-Supervised-ClaimIdentification


title: On the Approximation Power of Two-Layer Networks of Random ReLUs

abstract: This paper considers the following question: how well can depth-two ReLU networks with randomly initialized bottom-level weights represent smooth functions? We give near-matching upper- and lower-bounds for $L_2$-approximation in terms of the Lipschitz constant, the desired accuracy, and the dimension of the problem, as well as similar results in terms of Sobolev norms. Our positive results employ tools from harmonic analysis and ridgelet representation theory, while our lower-bounds are based on (robust versions of) dimensionality arguments.


title: Environment Predictive Coding for Embodied Agents

abstract: We introduce environment predictive coding, a self-supervised approach to learn environment-level representations for embodied agents. In contrast to prior work on self-supervised learning for images, we aim to jointly encode a series of images gathered by an agent as it moves about in 3D environments. We learn these representations via a zone prediction task, where we intelligently mask out portions of an agent's trajectory and predict them from the unmasked portions, conditioned on the agent's camera poses. By learning such representations on a collection of videos, we demonstrate successful transfer to multiple downstream navigation-oriented tasks. Our experiments on the photorealistic 3D environments of Gibson and Matterport3D show that our method outperforms the state-of-the-art on challenging tasks with only a limited budget of experience.


title: Microscopic Patterns in the 2D Phase-Field-Crystal Model

abstract: Using the recently developed theory of rigorously validated numerics, we address the Phase-Field-Crystal (PFC) model at the microscopic (atomistic) level. We show the existence of critical points and local minimizers associated with "classical" candidates, grain boundaries, and localized patterns. We further address the dynamical relationships between the observed patterns for fixed parameters and across parameter space, then formulate several conjectures on the dynamical connections (or orbits) between steady states.


title: Simulated annealing from continuum to discretization: a convergence analysis via the Eyring--Kramers law

abstract: We study the convergence rate of continuous-time simulated annealing $(X_t; \, t \ge 0)$ and its discretization $(x_k; \, k =0,1, \ldots)$ for approximating the global optimum of a given function $f$. We prove that the tail probability $\mathbb{P}(f(X_t) > \min f +\delta)$ (resp. $\mathbb{P}(f(x_k) > \min f +\delta)$) decays polynomial in time (resp. in cumulative step size), and provide an explicit rate as a function of the model parameters. Our argument applies the recent development on functional inequalities for the Gibbs measure at low temperatures -- the Eyring-Kramers law. In the discrete setting, we obtain a condition on the step size to ensure the convergence.


title: MUFASA: Multimodal Fusion Architecture Search for Electronic Health Records

abstract: One important challenge of applying deep learning to electronic health records (EHR) is the complexity of their multimodal structure. EHR usually contains a mixture of structured (codes) and unstructured (free-text) data with sparse and irregular longitudinal features -- all of which doctors utilize when making decisions. In the deep learning regime, determining how different modality representations should be fused together is a difficult problem, which is often addressed by handcrafted modeling and intuition. In this work, we extend state-of-the-art neural architecture search (NAS) methods and propose MUltimodal Fusion Architecture SeArch (MUFASA) to simultaneously search across multimodal fusion strategies and modality-specific architectures for the first time. We demonstrate empirically that our MUFASA method outperforms established unimodal NAS on public EHR data with comparable computation costs. In addition, MUFASA produces architectures that outperform Transformer and Evolved Transformer. Compared with these baselines on CCS diagnosis code prediction, our discovered models improve top-5 recall from 0.88 to 0.91 and demonstrate the ability to generalize to other EHR tasks. Studying our top architecture in depth, we provide empirical evidence that MUFASA's improvements are derived from its ability to both customize modeling for each data modality and find effective fusion strategies.


title: A Poisson map from kinetic theory to hydrodnamics with non-constant entropy

abstract: Kinetic theory describes a dilute monatomic gas using a distribution function $f(q,p,t)$, the expected phase-space density of particles. The distribution function evolves according to the collisionless Boltzmann equation in the high Knudsen number limit. Fluid dynamics is an alternative description of the gas using hydrodynamic variables that are functions of position and time only. These hydrodynamic variables evolve according to the compressible Euler equations in the inviscid limit. Both systems are noncanonical Hamiltonian systems. Each configuration space is an infinite-dimensional Poisson manifold, and the dynamics is the flow generated by a Hamiltonian functional via a Poisson bracket. We construct a map $\mathcal{J}_1$ from the space of distribution functions to the space of hydrodynamic variables that respects the Poisson brackets on the two spaces i.e. a Poisson map. It maps the $p$-integral of the Boltzmann entropy $f\log f$ to the hydrodynamic entropy density. This map belongs to a family of Poisson maps to spaces that include generalised entropy densities as additional hydrodynamic variables. The whole family can be generated from the Taylor expansion of a further Poisson map that depends on a formal parameter. If the kinetic-theory Hamiltonian factors through the Poisson map $\mathcal{J}_1$, an exact reduction of kinetic theory to fluid dynamics is possible. This is not the case, but by ignoring the relative entropy of $f$ to its local Maxwellian, we can construct an approximate Hamiltonian that factors through the map. The resulting reduced Hamiltonian generates the compressible Euler equations. We can thus derive the compressible Euler equations as a Hamiltonian approximation to kinetic theory. We also give an analogous Hamiltonian derivation of the compressible Euler--Poisson equations with non-constant entropy, starting from the Vlasov--Poisson equation.


title: Rate-dependent effects of lidocaine on cardiac dynamics: Development and analysis of a low-dimensional drug-channel interaction model

abstract: State-dependent Na+ channel blockers are often prescribed to treat cardiac arrhythmias, but many Na+ channel blockers are known to have pro-arrhythmic side effects. While the anti and proarrhythmic potential of a Na+ channel blocker is thought to depend on the characteristics of its rate-dependent block, the mechanisms linking these two attributes are unclear. Furthermore, how specific properties of rate-dependent block arise from the binding kinetics of a particular drug is poorly understood. Here, we examine the rate-dependent effects of the Na+ channel blocker lidocaine by constructing and analyzing a novel drug-channel interaction model. First, we identify the predominant mode of lidocaine binding in a 24 variable Markov model for lidocaine-Na+ channel interaction by Moreno et al. We then develop a novel 3-variable lidocaine-Na+ channel interaction model that incorporates only the predominant mode of drug binding. Our low-dimensional model replicates the extensive voltage-clamp data used to parameterize the Moreno et al. model. Furthermore, the effects of lidocaine on action potential upstroke velocity and conduction velocity in our model are similar to those predicted by the Moreno et al. model. By exploiting the low-dimensionality of our model, we derive an algebraic expression for level of rate-dependent block as a function of pacing frequency, restitution properties, diastolic and plateau potentials, and drug binding rate constants. Our model predicts that the level of rate-dependent block is sensitive to alterations in restitution properties and increases in diastolic potential, but it is insensitive to variations in the shape of the action potential waveform and lidocaine binding rates.


title: Bianchi IX and VIII Quantum Cosmology with a Cosmological Constant, Aligned Electromagnetic Field, and Scalar Field

abstract: We investigate the quantum cosmologies of the Bianchi IX and VIII models when a cosmological constant, aligned electromagnetic field and free scalar field are present. The conserved quantity $p_{\phi}$ associated with our free scalar field results in $\phi$ classically being a quantity which monotonically increases with respect to time, thus allowing it to fulfill the role of an 'emergent' internal clock for our constrained quantum systems. We embark on this investigation to better understand how matter sources can affect general anisotropic quantum cosmologies. To aid us we use the Euclidean-signature semi classical method to obtain our wave functions and analyze them. In addition we study briefly the quantum Taub models when an $e^{ \phi}$ potential and aligned electromagnetic field are present. One of the interesting things we found was that our aligned electromagnetic field, depending on how strong it is, can create or destroy geometric states in our 'excited' state wave functions that quantum Bianchi IX universes can tunnel in and out of. This creation of a state is somewhat similar to how non-commutativity in the minisuperspace variables can cause new quantum states to emerge in the quantum Kantowski-Sachs and Bianchi I models. Our results further show the utility of the Euclidean-signature semi classical method towards tackling Lorentzian signature problems without having to invoke a Wick rotation. This feature of not needing to apply a Wick rotation makes this method potentially very useful for tackling a variety of problems in bosonic relativistic field theory and quantum gravity.


title: Horizontally Fused Training Array: An Effective Hardware Utilization Squeezer for Training Novel Deep Learning Models

abstract: Driven by the tremendous effort in researching novel deep learning (DL) algorithms, the training cost of developing new models increases staggeringly in recent years. We analyze GPU cluster usage statistics from a top research institute for more insights into the hardware efficiency achieved by typical DL training jobs. Our study reveals that single-accelerator training jobs can dominate the cluster-wide resource consumption when launched repetitively (e.g., for hyper-parameter tuning) while severely under-utilizing the hardware. Fortunately, we observe that such workloads have the following unique characteristics: (i) the models among jobs often have the same types of operators with the same shapes, and (ii) the inter-model horizontal fusion of such operators is mathematically equivalent to other already well-optimized operators. Thus, to help DL researchers and practitioners effectively improve the hardware utilization of their novel DL training workloads, we propose Horizontally Fused Training Array (HFTA). HFTA is a new DL framework extension library that horizontally fuses the models from different repetitive jobs deeply down to operators and then trains them simultaneously on a shared accelerator. To show the generality of our solution, we apply HFTA to six DL models training on state-of-the-art accelerators (GPUs and TPUs). Our results indicate that HFTA is highly effective in improving hardware utilization and achieves up to $15.1 \times$ higher training throughput vs. the standard practice of running each job on a separate accelerator.


title: LAMMPS Implementation of Rapid Artificial Neural Network Derived Interatomic Potentials

abstract: While machine learning approaches have been successfully used to represent interatomic potentials, their speed has typically lagged behind conventional formalisms. This is often due to the complexity of the structural fingerprints used to describe the local atomic environment and the large cutoff radii and neighbor lists used in the calculation of these fingerprints. Even recent machine learned methods are at least 10 times slower than traditional formalisms. An implementation of a rapid artificial neural network (RANN) style potential in the LAMMPS molecular dynamics package is presented here which utilizes angular screening to reduce computational complexity without reducing accuracy. For the smallest neural network architectures, this formalism rivals the modified embedded atom method (MEAM) for speed and accuracy, while the networks approximately one third as fast as MEAM were capable of reproducing the training database with chemical accuracy. The numerical accuracy of the LAMMPS implementation is assessed by verifying conservation of energy and agreement between calculated forces and pressures and the observed derivatives of the energy as well as by assessing the stability of the potential in dynamic simulation. The potential style is tested using a force field for magnesium and the computational efficiency for a variety of architectures is compared to a traditional potential models as well as alternative ANN formalisms. The predictive accuracy is found to rival that of slower methods.


title: Microscopic Calculation of Fission Product Yields with Particle Number Projection

abstract: Fission fragments' charge and mass distribution is an important input to applications ranging from basic science to energy production or nuclear non-proliferation. In simulations of nucleosynthesis or calculations of superheavy elements, these quantities must be computed from models, as they are needed in nuclei where no experimental information is available. Until now, standard techniques to estimate these distributions were not capable of accounting for fine-structure effects, such as the odd-even staggering of the charge distributions. In this work, we combine a fully-microscopic collective model of fission dynamics with a recent extension of the particle number projection formalism to provide the highest-fidelity prediction of the primary fission fragment distributions for the neutron-induced fission of $^{235}$U and $^{239}$Pu. We show that particle number projection is an essential ingredient to reproduce odd-even staggering in the charge yields and benchmark the performance of various empirical probability laws that could simulate its effect. This new approach also enables for the first time the realistic determination of two-dimensional isotopic yields within nuclear density functional theory.


title: Control of Plasmons in Topological Insulators via Local Perturbations

abstract: We use a fully quantum mechanical approach to demonstrate control of plasmonic excitations in prototype models of topological insulators by molecule-scale perturbations. Strongly localized surface plasmons are present in the host systems, arising from the topologically non-trivial single-particle edge states. A numerical evaluation of the RPA equations for the perturbed systems reveals how the positions and the internal electronic structure of the added molecules affect the degeneracy of the locally confined collective excitations, i.e., shifting the plasmonic energies of the host system and changing their spatial charge density profile. In particular, we identify conditions under which significant charge transfer from the host system to the added molecules occurs. Furthermore, the induced field energy density in the perturbed topological systems due to external electric fields is determined.


title: Enhanced Coalbed Methane Extraction by Geothermal Stimulation in Deep Coal Mines: An Appraisal

abstract: Coalbed methane embedded in coal seams, is an unconventional energy resource as well as a hazardous gas existing in mining industries, which attracts lots of global attention. As the largest coal producer, the mining industry in China had to deal with many hazards induced by methane for decades. To solve this issue, underground methane extraction is commonly used in underground coal mines. However, underground methane extraction is hampered by low production rate and low efficiency because of slow gas emission from coal primarily controlled by gas desorption and permeability. It is well known that temperature has a great impact on gas sorption. The higher the temperature the larger the desorption rate. As the depth of coal mines increases beyond 1000m coal mines suffer elevated air temperatures caused by the natural geothermal gradient. The elevated temperature in such mines provides a potential economical way for geothermal energy extraction and utilization in deep coal mines which can largely cut the expenses of installation and operation maintenance. Therefore, a novel method is proposed to enhance underground methane extraction by deep heat stimulation. This paper mainly presents an assessment of previous and ongoing research in the related field and provides a first feasibility analysis of this method applied in the underground environment. The technique proposed in this early appraisal is deemed significant for coalbed methane drainage enhancing the productivity of deep coal mines by geothermal technology and can also be extended for many applications in relevant areas such as shale gas, and tight oil.


title: Levelness versus almost Gorensteinness of edge rings of complete multipartite graphs

abstract: Levelness and almost Gorensteinness are well-studied properties on graded rings as a generalized notion of Gorensteinness. In the present paper, we study those properties for the edge rings of the complete multipartite graphs, denoted by $\Bbbk[K_{r_1,\ldots,r_n}]$ with $1 \leq r_1 \leq \cdots \leq r_n$. We give the complete characterization of which $\Bbbk[K_{r_1,\ldots,r_n}]$ is level in terms of $n$ and $r_1,\ldots,r_n$. Similarly, we also give the complete characterization of which $\Bbbk[K_{r_1,\ldots,r_n}]$ is almost Gorenstein in terms of $n$ and $r_1,\ldots,r_n$.


title: Tournaments with maximal decomposability

abstract: Given a tournament $T$, a module of $T$ is a subset $M$ of $V(T)$ such that for $x, y\in M$ and $v\in V(T)\setminus M$, $(x,v)\in A(T)$ if and only if $(y,v)\in A(T)$. The trivial modules of $T$ are $\emptyset$, $\{u\}$ $(u\in V(T))$ and $V(T)$. The tournament $T$ is indecomposable if all its modules are trivial; otherwise it is decomposable. The decomposability index of $T$, denoted by $\delta(T)$, is the smallest number of arcs of $T$ that must be reversed to make $T$ indecomposable. In a previous paper, we proved that for $n \geq 5$, we have $\delta(n) = \left\lceil \frac{n+1}{4} \right\rceil$, where $\delta(n)$ is the maximum of $\delta(T)$ over the tournaments $T$ with $n$ vertices. In this paper, we characterize the tournaments $T$ with $\delta$-maximal decomposability, i.e., such that $\delta(T)=\delta(\vert T\vert)$.


title: On Multi-Human Multi-Robot Remote Interaction: A Study of Transparency, Inter-Human Communication, and Information Loss in Remote Interaction

abstract: In this paper, we investigate how to design an effective interface for remote multi-human multi-robot interaction. While significant research exists on interfaces for individual human operators, little research exists for the multi-human case. Yet, this is a critical problem to solve to make complex, large-scale missions achievable in which direct human involvement is impossible or undesirable, and robot swarms act as a semi-autonomous agents. This paper's contribution is twofold. The first contribution is an exploration of the design space of computer-based interfaces for multi-human multi-robot operations. In particular, we focus on information transparency and on the factors that affect inter-human communication in ideal conditions, i.e., without communication issues. Our second contribution concerns the same problem, but considering increasing degrees of information loss, defined as intermittent reception of data with noticeable gaps between individual receipts. We derived a set of design recommendations based on two user studies involving 48 participants.


title: Binary-driven stellar rotation evolution at the main-sequence turn-off in star clusters

abstract: The impact of stellar rotation on the morphology of star cluster colour-magnitude diagrams is widely acknowledged. However, the physics driving the distribution of the equatorial rotation velocities of main-sequence turn-off (MSTO) stars is as yet poorly understood. Using Gaia Data Release 2 photometry and new Southern African Large Telescope medium-resolution spectroscopy, we analyse the intermediate-age ($\sim1\,$Gyr-old) Galactic open clusters NGC 3960, NGC 6134 and IC 4756 and develop a novel method to derive their stellar rotation distributions based on SYCLIST stellar rotation models. Combined with literature data for the open clusters NGC 5822 and NGC 2818, we find a tight correlation between the number ratio of slow rotators and the clusters' binary fractions. The blue-main-sequence stars in at least two of our clusters are more centrally concentrated than their red-main-sequence counterparts. The origin of the equatorial stellar rotation distribution and its evolution remains as yet unidentified. However, the observed correlation in our open cluster sample suggests a binary-driven formation mechanism.


title: Disorder-robust high-field superconducting phase of FeSe single crystals

abstract: When exposed to high magnetic fields, certain materials manifest an exotic superconducting (SC) phase that has attracted considerable attention. A proposed explanation for the origin of the high-field SC phase is the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. This state is characterized by inhomogeneous superconductivity, where the Cooper pairs have finite center-of-mass momenta. Recently, the high-field SC phase was observed in FeSe, and it was deemed to originate from the FFLO state. Here, we synthesize FeSe single crystals with different levels of disorder. The level of disorder is expressed by the ratio of the mean free path to the coherence length and ranges between 35 and 1.2. The upper critical field \textit{B}$_{\rm{c}2}$ was obtained by both resistivity and magnetic torque measurements over a wide range of temperatures, which went as low as $\sim$0.5 K, and magnetic fields, which went up to $\sim$38 T along the \textit{c} axis and in the \textit{ab} plane. In the high-field region parallel to the \textit{ab} plane, an unusual SC phase was confirmed in all the crystals, and the phase was found to be robust against disorder. This result suggests that the high-field SC phase in FeSe is not a conventional FFLO state.


title: NENCI-2021 Part I: A Large Benchmark Database of Non-Equilibrium Non-Covalent Interactions Emphasizing Close Intermolecular Contacts

abstract: In this work, we present NENCI-2021, a benchmark database of approximately $8,000$ non-equilibrium non-covalent interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the $101$ molecular dimers in the widely used S66 and S101 databases, and extends the scope of these works by: (i) including $40$ cation- and anion-$\pi$ complexes, a fundamentally important class of non-covalent interactions (NCIs) that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all $141$ intermolecular potential energy surfaces (PES) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from $0.7\times\mathrm{-}1.1\times$ the equilibrium separation) as well as nine intermolecular angles per distance (five for each ion-$\pi$ complex), yielding an extensive database of $7,763$ benchmark intermolecular interaction energies ($E_{\rm int}$) obtained at the CCSD(T)/CBS level of theory. In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent -- these interactions in particular pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed SAPT-based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021.


title: Entropic Dynamics of Networks

abstract: Here we present the entropic dynamics formalism for networks. That is, a framework for the dynamics of graphs meant to represent a network derived from the principle of maximum entropy and the rate of transition is obtained taking into account the natural information geometry of probability distributions. We apply this framework to the Gibbs distribution of random graphs obtained with constraints on the node connectivity. The information geometry for this graph ensemble is calculated and the dynamical process is obtained as a diffusion equation. We compare the steady state of this dynamics to degree distributions found on real-world networks.


title: Convergence of Combinatorial Gravity

abstract: We present a new regularisation of Euclidean Einstein gravity in terms of (sequences of) graphs. In particular, we define a discrete Einstein-Hilbert action that converges to its manifold counterpart on sufficiently dense random geometric graphs (more generally on any sequence of graphs that converges to the manifold in the sense of Gromov-Hausdorff). Our construction relies crucially on the Ollivier curvature of optimal transport theory. Our methods also allow us to define an analogous discrete action for Klein-Gordon fields. These results may be taken as the basis for a combinatorial approach to quantum gravity where we seek to generate graphs that approximate manifolds as metric-measure structures.


title: Speed of sound in dense matter and two families of compact stars

abstract: The existence of massive compact stars $(M\gtrsim 2.1 M_{\odot})$ implies that the conformal limit of the speed of sound $c_s^2=1/3$ is violated if those stars have a crust of ordinary nuclear matter. Here we show that, if the most massive objects are strange quark stars, i.e. stars entirely composed of quarks, the conformal limit can be respected while observational limits on those objects are also satisfied. By using astrophysical data associated with those massive stars, derived from electromagnetic and gravitational wave signals, we show, within a Bayesian analysis framework and by adopting a constant speed of sound equation of state, that the posterior distribution of $c_s^2$ is peaked around 0.3, and the maximum mass of the most probable equation of state is $\sim 2.13 M_{\odot}$. We discuss which new data would require a violation of the conformal limit even when considering strange quark stars, in particular we analyze the possibility that the maximum mass of compact stars is larger than $2.5M_{\odot}$, as it would be if the secondary component of GW190814 is a compact star and not a black hole. Finally, we discuss how the new data for PSR J0740+6620 obtained by the NICER collaboration compare with our analysis (not based on them) and with other possible interpretations.


title: A Practical Coding Scheme for the BSC with Feedback

abstract: We provide a practical implementation of the rubber method of Ahlswede et al. for binary alphabets. The idea is to create the "skeleton" sequence therein via an arithmetic decoder designed for a particular $k$-th order Markov chain. For the stochastic binary symmetric channel, we show that the scheme is nearly optimal in a strong sense for certain parameters.


title: Wave-function engineering via conditional quantum teleportation with non-Gaussian entanglement resource

abstract: We propose and analyze a setup to tailor the wave functions of the quantum states. Our setup is based on the quantum teleportation circuit, but instead of the usual two-mode squeezed state, two-mode non-Gaussian entangled state is used. Using this setup, we can generate various classes of quantum states such as Schr\"odinger cat states, four-component cat states, superpositions of Fock states, and cubic phase states. These results demonstrate the versatility of our system as a state generator and suggest that conditioning using homodyne measurements is an important tool in the generations of the non-Gaussian states in complementary to the photon number detection.


title: Proof of two multivariate $q$-binomial sums arising in Gromov-Witten theory

abstract: We prove two multivariate $q$-binomial identities conjectured by Bousseau, Brini and van Garrel [Stable maps to Looijenga pairs, arXiv:2011.08830] which give generating series for Gromov-Witten invariants of two specific log Calabi-Yau surfaces. The key identity in all the proofs is Jackson's $q$-analogue of the Pfaff-Saalsch\"utz summation formula from the theory of basic hypergeometric series.


title: Quantum Phases of Self-Bound Droplets of Bose-Bose Mixtures

abstract: We systematically investigate the ground-state properties of self-bound droplets of quasi-two-dimensional binary Bose gases by using the Gaussian state theory. We find that quantum droplets consists two macroscopic squeezed phases and a macroscopic coherent phase. We map out the phase diagram and determine all phase boundaries via both numerical and nearly analytical methods. In particular, we find three easily accessible signatures for the quantum phases and the stablization mechanism of the self-bound droplets by precisely measuring their radial size. Our studies indicate that binary droplets represent an ideal platform for in-depth investigations of the quantum nature of the droplet state.


title: Resolving the History of Life on Earth by Seeking Life As We Know It on Mars

abstract: An origin of Earth life on Mars would resolve significant inconsistencies between the inferred history of life and Earth's geologic history. Life as we know it utilizes amino acids, nucleic acids, and lipids for the metabolic, informational, and compartment-forming subsystems of a cell. Such building blocks may have formed simultaneously from cyanosulfidic chemical precursors in a planetary surface scenario involving ultraviolet light, wet-dry cycling, and volcanism. However, early Earth was a water world, and the timing of the rise of oxygen on Earth is inconsistent with final fixation of the genetic code in response to oxidative stress. A cyanosulfidic origin of life could have taken place on Mars via photoredox chemistry, facilitated by orders of magnitude more sub-aerial crust than early Earth, and an earlier transition to oxidative conditions. Meteoritic bombardment may have generated transient habitable environments and ejected and transferred life to Earth. The Mars 2020 Perseverance Rover offers an unprecedented opportunity to confirm or refute evidence consistent with a cyanosulfidic origin of life on Mars, search for evidence of ancient life, and constrain the evolution of Mars' oxidation state over time. We should seek to prove or refute a Martian origin for life on Earth alongside other possibilities.


title: Operational Semantics with Hierarchical Abstract Syntax Graphs

abstract: This is a motivating tutorial introduction to a semantic analysis of programming languages using a graphical language as the representation of terms, and graph rewriting as a representation of reduction rules. We show how the graphical language automatically incorporates desirable features, such as alpha-equivalence and how it can describe pure computation, imperative store, and control features in a uniform framework. The graph semantics combines some of the best features of structural operational semantics and abstract machines, while offering powerful new methods for reasoning about contextual equivalence. All technical details are available in an extended technical report by Muroya and the author and in Muroya's doctoral dissertation.


title: On Stochastic Rewriting and Combinatorics via Rule-Algebraic Methods

abstract: Building upon the rule-algebraic stochastic mechanics framework, we present new results on the relationship of stochastic rewriting systems described in terms of continuous-time Markov chains, their embedded discrete-time Markov chains and certain types of generating function expressions in combinatorics. We introduce a number of generating function techniques that permit a novel form of static analysis for rewriting systems based upon marginalizing distributions over the states of the rewriting systems via pattern-counting observables.


title: Wind Field Reconstruction with Adaptive Random Fourier Features

abstract: We investigate the use of spatial interpolation methods for reconstructing the horizontal near-surface wind field given a sparse set of measurements. In particular, random Fourier features is compared to a set of benchmark methods including Kriging and Inverse distance weighting. Random Fourier features is a linear model $\beta(\pmb x) = \sum_{k=1}^K \beta_k e^{i\omega_k \pmb x}$ approximating the velocity field, with frequencies $\omega_k$ randomly sampled and amplitudes $\beta_k$ trained to minimize a loss function. We include a physically motivated divergence penalty term $|\nabla \cdot \beta(\pmb x)|^2$, as well as a penalty on the Sobolev norm. We derive a bound on the generalization error and derive a sampling density that minimizes the bound. Following (arXiv:2007.10683 [math.NA]), we devise an adaptive Metropolis-Hastings algorithm for sampling the frequencies of the optimal distribution. In our experiments, our random Fourier features model outperforms the benchmark models.


title: Parallel Independence in Attributed Graph Rewriting

abstract: In order to define graph transformations by the simultaneous application of concurrent rules, we have adopted in previous work a structure of attributed graphs stable by unions. We analyze the consequences on parallel independence, a property that characterizes the possibility to resort to sequential rewriting. This property turns out to depend not only on the left-hand side of rules, as in algebraic approaches to graph rewriting, but also on their right-hand side. It is then shown that, of three possible definitions of parallel rewriting, only one is convenient in the light of parallel independence.


title: Finite Difference Weerakoon-Fernando Method to solve nonlinear equations without using derivatives

abstract: This research was mainly conducted to explore the possibility of formulating an efficient algorithm to find roots of nonlinear equations without using the derivative of the function. The Weerakoon-Fernando method had been taken as the base in this project to find a new method without the derivative since Weerakoon-Fernando method gives 3rd order convergence. After several unsuccessful attempts we were able to formulate the Finite Difference Weerakoon-Fernando Method (FDWFM) presented here. We noticed that the FDWFM approaches the root faster than any other existing method in the absence of the derivatives as an example, the popular nonlinear equation solver such as secant method (order of convergence is 1.618) in the absence of the derivative. And the FDWFM had three function evaluations and secant method had two function evaluations. By implementing FDWFM on nonlinear equations with complex roots and also on systems of nonlinear equations, we received very encouraging results. When applying the FDWFM to systems of nonlinear equations, we resolved the involvement of the Jacobian problem by following the procedure in the Broyden's method. The computational order of convergence of the FDWFM was close to 2.5 for all these cases. This will undoubtedly provide scientists the efficient numerical algorithm, that doesn't need the derivative of the function to solve nonlinear equations, that they were searching for over centuries.


title: Verifying Security Vulnerabilities in Large Software Systems using Multi-Core k-Induction

abstract: Computer-based systems have been used to solve several domain problems, such as industrial, military, education, and wearable. Those systems need high-quality software to guarantee security and safety. We advocate that Bounded Model Checking (BMC) techniques can detect security vulnerabilities in the early stages of development processes. However, this technique struggles to scale up and verify large software commonly found on computer-based systems. Here, we develop and evaluate a pragmatic approach to verify large software systems using a state-of-the-art bounded model checker. In particular, we pre-process the input source-code files and then guide the model checker to explore the code systematically. We also present a multi-core implementation of the k-induction proof algorithm to verify and falsify large software systems iteratively. Our experimental results using the Efficient SMT-based Model Checker (ESBMC) show that our approach can guide ESBMC to efficiently verify large software systems. We evaluate our approach using the PuTTY application to verify 136 files and 2803 functions in less than 86 minutes, and the SlimGuard allocator, where we have found real security vulnerabilities confirmed by the developers. We conclude that our approach can successfully guide a bounded model checker to verify large software systems systematically.


title: Direct Fidelity Estimation of Quantum States using Machine Learning

abstract: In almost all quantum applications, one of the key steps is to verify that the fidelity of the prepared quantum state meets expectations. In this Letter, we propose a new approach solving this problem using machine-learning techniques. Compared to other fidelity estimation methods, our method is applicable to arbitrary quantum states, the number of required measurement settings is small, and this number does not increase with the size of the system. For example, for a general five-qubit quantum state, only four measurement settings are required to predict its fidelity with $\pm1\%$ precision in a nonadversarial scenario. This machine-learning-based approach for estimating quantum state fidelity has the potential to be widely used in the field of quantum information.


title: Analytic Design of Accelerated Adiabatic Gates in Realistic Qubits: General Theory and Applications to Superconducting Circuits

abstract: Shortcuts to adiabaticity is a general method for speeding up adiabatic quantum protocols, and has many potential applications in quantum information processing. Unfortunately, analytically constructing shortcuts to adiabaticity for systems having complex interactions and more than a few levels is a challenging task. This is usually overcome by assuming an idealized Hamiltonian [e.g., only a limited subset of energy levels are retained, and the rotating-wave approximation (RWA) is made]. Here we develop an $analytic$ approach that allows one to go beyond these limitations. Our method is general and results in analytically derived pulse shapes that correct both nonadiabatic errors as well as non-RWA errors. We also show that our approach can yield pulses requiring a smaller driving power than conventional nonadiabatic protocols. We show in detail how our ideas can be used to analytically design high-fidelity single-qubit "tripod" gates in a realistic superconducting fluxonium qubit.


title: MeInGame: Create a Game Character Face from a Single Portrait

abstract: Many deep learning based 3D face reconstruction methods have been proposed recently, however, few of them have applications in games. Current game character customization systems either require players to manually adjust considerable face attributes to obtain the desired face, or have limited freedom of facial shape and texture. In this paper, we propose an automatic character face creation method that predicts both facial shape and texture from a single portrait, and it can be integrated into most existing 3D games. Although 3D Morphable Face Model (3DMM) based methods can restore accurate 3D faces from single images, the topology of 3DMM mesh is different from the meshes used in most games. To acquire fidelity texture, existing methods require a large amount of face texture data for training, while building such datasets is time-consuming and laborious. Besides, such a dataset collected under laboratory conditions may not generalized well to in-the-wild situations. To tackle these problems, we propose 1) a low-cost facial texture acquisition method, 2) a shape transfer algorithm that can transform the shape of a 3DMM mesh to games, and 3) a new pipeline for training 3D game face reconstruction networks. The proposed method not only can produce detailed and vivid game characters similar to the input portrait, but can also eliminate the influence of lighting and occlusions. Experiments show that our method outperforms state-of-the-art methods used in games.


title: Coevolution of theoretical and applied research: a case study of graphene research by temporal and geographic analysis

abstract: As a part of science of science (SciSci) research, the evolution of scientific disciplines has been attracting a great deal of attention recently. This kind of discipline level analysis not only give insights of one particular field but also shed light on general principles of scientific enterprise. In this paper we focus on graphene research, a fast growing field covers both theoretical and applied study. Using co-clustering method, we split graphene literature into two groups and confirm that one group is about theoretical research (T) and another corresponds to applied research (A). We analyze the proportion of T/A and found applied research becomes more and more popular after 2007. Geographical analysis demonstrated that countries have different preference in terms of T/A and they reacted differently to research trend. The interaction between two groups has been analyzed and shows that T extremely relies on T and A heavily relies on A, however the situation is very stable for T but changed markedly for A. No geographic difference is found for the interaction dynamics. Our results give a comprehensive picture of graphene research evolution and also provide a general framework which is able to analyze other disciplines.


title: Moments of central values of quartic Dirichlet $L$-functions

abstract: In this paper, we study moments of the central values of quartic Dirichlet $L$-functions and establish quantitative non-vanishing result for these $L$-values.


title: Sampling in Combinatorial Spaces with SurVAE Flow Augmented MCMC

abstract: Hybrid Monte Carlo is a powerful Markov Chain Monte Carlo method for sampling from complex continuous distributions. However, a major limitation of HMC is its inability to be applied to discrete domains due to the lack of gradient signal. In this work, we introduce a new approach based on augmenting Monte Carlo methods with SurVAE Flows to sample from discrete distributions using a combination of neural transport methods like normalizing flows and variational dequantization, and the Metropolis-Hastings rule. Our method first learns a continuous embedding of the discrete space using a surjective map and subsequently learns a bijective transformation from the continuous space to an approximately Gaussian distributed latent variable. Sampling proceeds by simulating MCMC chains in the latent space and mapping these samples to the target discrete space via the learned transformations. We demonstrate the efficacy of our algorithm on a range of examples from statistics, computational physics and machine learning, and observe improvements compared to alternative algorithms.


title: Self-Interacting Dark Matter and Small-Scale Gravitational Lenses in Galaxy Clusters

abstract: Recently, Meneghetti et al. reported an excess of small-scale gravitational lenses in galaxy clusters. We study its implications for self-interacting dark matter (SIDM), compared with standard cold dark matter (CDM). We design controlled N-body simulations that incorporate observational constraints. The presence of early-type galaxies in cluster substructures can deepen gravitational potential and reduce tidal mass loss. Both scenarios require a relatively high baryon concentration in the substructure to accommodate the lensing measurements, and their tangential caustics are similar. The SIDM substructure can experience gravothermal collapse and produce a steeper density profile than its CDM counterpart, leading to a larger radial galaxy-galaxy strong lensing cross section, although this effect is hard to observe. Our results indicate SIDM can provide a unified explanation to small-scale lenses in galaxy clusters and stellar motions in dwarf galaxies.


title: A review of motion planning algorithms for intelligent robotics

abstract: We investigate and analyze principles of typical motion planning algorithms. These include traditional planning algorithms, supervised learning, optimal value reinforcement learning, policy gradient reinforcement learning. Traditional planning algorithms we investigated include graph search algorithms, sampling-based algorithms, and interpolating curve algorithms. Supervised learning algorithms include MSVM, LSTM, MCTS and CNN. Optimal value reinforcement learning algorithms include Q learning, DQN, double DQN, dueling DQN. Policy gradient algorithms include policy gradient method, actor-critic algorithm, A3C, A2C, DPG, DDPG, TRPO and PPO. New general criteria are also introduced to evaluate performance and application of motion planning algorithms by analytical comparisons. Convergence speed and stability of optimal value and policy gradient algorithms are specially analyzed. Future directions are presented analytically according to principles and analytical comparisons of motion planning algorithms. This paper provides researchers with a clear and comprehensive understanding about advantages, disadvantages, relationships, and future of motion planning algorithms in robotics, and paves ways for better motion planning algorithms.


title: Noise and error analysis and optimization in particle-based kinetic plasma simulations

abstract: We analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We describe kernel density estimation for continuous values of the spatial variable, $x$, and compute the covariance matrix for uniform true density. The band width of the covariance matrix is related to the width of the kernel. We find the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step. The effect of these negative correlations on the electric field, computed by Gauss's law, is that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge. For non-constant density still with continuous $x$, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length $l$ and kernel width $h$, having too few particles within $h$ leads to too much variance; for $h$ large relative to $l$, there is too much smoothing of the density. The optimum between these two limits is found by BVO. We repeat the analysis for $x$ discretized on a grid, connecting particle shapes and estimation kernels. If the particle shape satisfies a sum rule, the charge deposited on the grid is conserved exactly. Particle shapes expressed as the convolution of a kernel with another kernel that satisfies the sum rule also obey the sum rule. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. We show good agreement with numerical computations.


title: Histogram Specification by Assignment of Optimal Unique Values

abstract: In this paper, we propose two novel algorithms for histogram specification and quantile transformation of data without local information. These are core techniques that can serve as building blocks for applications that require specifying the sample distribution of a given set of data. Histogram specification is best known for its image enhancement applications, whereas quantile transformation is typically employed in data preprocessing for data normalization. In signal processing, methods often require temporal or spatial information; in data preprocessing, methods work by interpolation or by approximation, drawing from results in computational statistics, and have a trade-off between speed and quality. It is nontrivial to accommodate for cases that do not have local information (e.g., tabular data) while also providing a fast, exact solution. For that, we take up a concept in image processing called group mapping law and propose an extension. The proposed extension allows us to formulate a convex functional where we look for the best approximation between the output unique values and the reference histogram. Then, we apply the ordered assignment solution, a result in optimal transport, to reconstruct the output from the optimal unique values. Two sets of results show the effectiveness of the proposed algorithms when compared to traditional and state-of-the-art methods. The proposed algorithms are fast, exact, and least $p$-norm optimal. Further, we define the algorithms as generic data processing methods. Thus, contributions from this paper can be easily incorporated in applications spanning many disciplines, especially in applied data science.


title: Airport Capacity and Performance in Europe -- A study of transport economics, service quality and sustainability

abstract: The purpose of this dissertation is to present an overview of the operational and financial performance of airports in Europe. In benchmarking studies, airports are assessed and compared with other airports based on key indicators from a technical and an economic point of view. The interest lies primarily in the question, which key figures can best measure the perception of quality of service from the point of view of the passenger for the services at an airport.


title: Reinforcing increase of {\Delta}TC in MgB2 smart meta-superconductors by adjusting the concentration of inhomogeneous phases

abstract: Incorporating with inhomogeneous phases with high electroluminescence (EL) intensity to prepare smart meta-superconductors (SMSCs) is an effective method of increasing the superconducting transition temperature (Tc) and has been confirmed in both MgB2 and Bi(Pb)SrCaCuO systems. However, the increase of {\Delta}Tc has been quite small because of the low optimal concentrations of inhomogeneous phases. In this work, three kinds of MgB2 raw materials, namely, aMgB2, bMgB2, and cMgB2, were prepared with particle sizes decreasing in order. Inhomogeneous phases, Y2O3:Eu3+ and Y2O3:Eu3+/Ag, were also prepared and doped into MgB2 to study the influence of doping concentration on the {\Delta}Tc of MgB2 with different particle sizes. Results show that reducing the MgB2 particle size increases the optimal doping concentration of inhomogeneous phases, thereby increasing {\Delta}Tc. The optimal doping concentrations for aMgB2, bMgB2, and cMgB2 are 0.5%, 0.8%, and 1.2%, respectively. The corresponding {\Delta}Tc values are 0.4, 0.9, and 1.2 K, respectively. This work open a new approach to reinforcing increase of {\Delta}Tc in MgB2 SMSCs.


title: Tilted Nonparametric Regression Function Estimation

abstract: This paper provides the theory about the convergence rate of the tilted version of linear smoother. We study tilted linear smoother, a nonparametric regression function estimator, which is obtained by minimizing the distance to an infinite order flat-top trapezoidal kernel estimator. We prove that the proposed estimator achieves a high level of accuracy. Moreover, it preserves the attractive properties of the infinite order flat-top kernel estimator. We also present an extensive numerical study for analysing the performance of two members of the tilted linear smoother class named tilted Nadaraya-Watson and tilted local linear in the finite sample. The simulation study shows that tilted Nadaraya-Watson and tilted local linear perform better than their classical analogs in some conditions in terms of Mean Integrated Squared Error (MISE). Finally, the performance of these estimators as well as the conventional estimators were illustrated by curve fitting to COVID-19 data for 12 countries and a dose-response data set.


title: Mainstreaming of conspiracy theories and misinformation

abstract: Parents - particularly moms - increasingly consult social media for support when taking decisions about their young children, and likely also when advising other family members such as elderly relatives. Minimizing malignant online influences is therefore crucial to securing their assent for policies ranging from vaccinations, masks and social distancing against the pandemic, to household best practices against climate change, to acceptance of future 5G towers nearby. Here we show how a strengthening of bonds across online communities during the pandemic, has led to non-Covid-19 conspiracy theories (e.g. fluoride, chemtrails, 5G) attaining heightened access to mainstream parent communities. Alternative health communities act as the critical conduits between conspiracy theorists and parents, and make the narratives more palatable to the latter. We demonstrate experimentally that these inter-community bonds can perpetually generate new misinformation, irrespective of any changes in factual information. Our findings show explicitly why Facebook's current policies have failed to stop the mainstreaming of non-Covid-19 and Covid-19 conspiracy theories and misinformation, and why targeting the largest communities will not work. A simple yet exactly solvable and empirically grounded mathematical model, shows how modest tailoring of mainstream communities' couplings could prevent them from tipping against establishment guidance. Our conclusions should also apply to other social media platforms and topics.


title: Interfacial Tension Modulation of Liquid Metal via Electrochemical Oxidation

abstract: This progress report summarizes recent studies of electrochemical oxidation to modulate the interfacial tension of gallium-based alloys. These alloys, which are liquid at ambient conditions, have the largest interfacial tension of any liquid at room temperature. The ability to modulate the tension offers the possibility to create forces that change the shape and position of the metal. It has been known since the late 1800s that electrocapillarity-the use of potential to modulate the electric double layer on the surface of metals in electrolyte-lowers the interfacial tension of liquid metal. Yet, this phenomenon can only achieve modest changes in interfacial tension since it is limited to potential windows that avoid reactions. A recent discovery suggests that reactions driven by the electrochemical oxidation of gallium alloys cause the interfacial tension to decrease from ~500 mN/m at 0 V to ~0 mN/m at ~0.8 V, a change in tension that goes well beyond what is possible via conventional electrocapillarity or surfactants. The changes in tension are reversible; reductive potentials return the metal back to a state of high interfacial tension. This report aims to summarize key work and introduce beginners to this field by including electrochemistry basics while addressing misconceptions. We discuss applications that utilize modulations in interfacial tension of liquid metal and conclude with remaining opportunities and challenges that need further investigation.


title: A Possible Artificial Intelligence Ecosystem Avatar: the Moorea case (IDEA)

abstract: High-throughput data collection techniques and largescale (cloud) computing are transforming our understanding of ecosystems at all scales by allowing the integration of multimodal data such as physics, chemistry, biology, ecology, fishing, economics and other social sciences in a common computational framework. We focus in this paper on a large scale data assimilation and prediction backbone based on Deep Stacking Networks (DSN) in the frame of the IDEA (Island Digital Ecosystem Avatars) project (Moorea Island), based on the subdivision of the island in watersheds and lagoon units. We also describe several kinds of raw data that can train and constrain such an ecosystem avatar model, as well as second level data such as ecological or physical indexes / indicators.


title: The critical behavior of Hegselmann-Krause opinion model with smart agents

abstract: The Hegselmann-Krause (HK) model allows one to characterize the continuous change of agents' opinions with the bounded confidence threshold $\varepsilon$. To consider the heterogeneity of agents in characteristics, we study the HK model on homogeneous and heterogeneous networks by introducing a kind of smart agent. Different from the averaging rule in opinion update of HK model, smart agents will consider, in updating their opinions, the environmental influence following the fact that an agent's behavior is often coupled with environmental changes. The environment is characterized by a parameter that represents the biased resource allocation between different cliques. We focus on the critical behavior of the underlying system. A phase transition point separating a complete consensus from the coexistence of different opinions is identified, which occurs at a critical value $\varepsilon_c$ for the bounded confidence threshold. We state analytically that $\varepsilon_c$ can take only one of two possible values, depending on the behavior of the average degree $k_a$ of a social graph, when agents are homogeneous in characteristics. Results also suggest that the phase transition point weakly depends on the network structure but is strongly correlated with the fraction of smart agents and the environmental parameter. We finally establish the finite size scaling law that stresses the role that the system size has in the underlying opinion dynamics. Meanwhile, introducing smart agents does not change the functional dependence between the time to reach a complete consensus and the system size. However, it can drive a complete consensus to be reached faster, for homogeneous networks that are far from the mean field limit.


title: Analysis of Temperature-to-Polarization Leakage in BICEP3 and Keck CMB Data from 2016 to 2018

abstract: The BICEP/Keck Array experiment is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background polarization from the South Pole in search of a primordial $B$-mode signature. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We use high-fidelity, in-situ measurements of the beam response to estimate the temperature-to-polarization (T $\rightarrow$ P) leakage in our latest data including observations from 2016 through 2018. This includes three years of BICEP3 observing at 95 GHz, and multifrequency data from Keck Array. Here we present band-averaged far-field beam maps, differential beam mismatch, and residual beam power (after filtering out the leading difference modes via deprojection) for these receivers. We show preliminary results of "beam map simulations," which use these beam maps to observe a simulated temperature (no $Q/U$) sky to estimate T $\rightarrow$ P leakage in our real data.


title: Self-duality of One-dimensional Quasicrystals with Spin-Orbit Interaction

abstract: Non-interacting spinless electrons in one-dimensional quasicrystals, described by the Aubry-Andr\'{e}-Harper (AAH) Hamiltonian with nearest neighbour hopping, undergoes metal to insulator transition (MIT) at a critical strength of the quasi-periodic potential. This transition is related to the self-duality of the AAH Hamiltonian. Interestingly, at the critical point, which is also known as the self-dual point, all the single particle wave functions are multifractal or non-ergodic in nature, while they are ergodic and delocalized (localized) below (above) the critical point. In this work, we have studied the one dimensional quasi-periodic AAH Hamiltonian in the presence of spin-orbit (SO) coupling of Rashba type, which introduces an additional spin conserving complex hopping and a spin-flip hopping. We have found that, although the self-dual nature of AAH Hamiltonian remains unaltered, the self-dual point gets affected significantly. Moreover, the effect of the complex and spin-flip hoppings are identical in nature. We have extended the idea of Kohn's localization tensor calculations for quasi-particles and detected the critical point very accurately. These calculations are followed by detailed multifractal analysis along with the computation of inverse participation ratio and von Neumann entropy, which clearly demonstrate that the quasi-particle eigenstates are indeed multifractal and non-ergodic at the critical point. Finally, we mapped out the phase diagram in the parameter space of quasi-periodic potential and SO coupling strength.


title: First- and Second-Moment Constrained Gaussian Channels

abstract: This paper studies the channel capacity of intensity-modulation direct-detection (IM/DD) visible light communication (VLC) systems under both optical and electrical power constraints. Specifically, it derives the asymptotic capacities in the high and low signal-to-noise ratio (SNR) regimes under peak, first-moment, and second-moment constraints. The results show that first- and second-moment constraints are never simultaneously active in the asymptotic low-SNR regime, and only in few cases in the asymptotic high-SNR regime. Moreover, the second-moment constraint is more stringent in the asymptotic low-SNR regime than in the high-SNR regime.


title: Experimental determination of the valence band offsets of $ZnGeN_2$ and $ZnGe_{0.94}Ga_{0.12}N_2$ with $GaN$

abstract: A predicted type-II staggered band alignment with an approximately $1.4 eV$ valence band offset at the $ZnGeN_2/GaN$ heterointerface has inspired novel band-engineered $III-N/ZnGeN_2$ heterostructure-based device designs for applications in high performance optoelectronics. We report on the determination of the valence band offset between metalorganic chemical vapor deposition grown $(ZnGe)_{1-x}Ga_{2x}N_2$, for $x = 0$ and $0.06$, and $GaN$ using X-ray photoemission spectroscopy. The valence band of $ZnGeN_2$ was found to lie $1.45-1.65 eV$ above that of $GaN$. This result agrees well with the value predicted by first-principles density functional theory calculations using the local density approximation for the potential profile and quasiparticle self-consistent GW calculations of the band edge states relative to the potential. For $(ZnGe)_{0.94}Ga_{0.12}N_2$ the value was determined to be $1.29 eV$, $~10-20\%$ lower than that of $ZnGeN_2$. The experimental determination of the large band offset between $ZnGeN_2$ and $GaN$ provides promising alternative solutions to address challenges faced with pure III-nitride-based structures and devices.


title: A Universal Framework for Featurization of Atomistic Systems

abstract: Molecular dynamics simulations are an invaluable tool in numerous scientific fields. However, the ubiquitous classical force fields cannot describe reactive systems, and quantum molecular dynamics are too computationally demanding to treat large systems or long timescales. Reactive force fields based on physics or machine learning can be used to bridge the gap in time and length scales, but these force fields require substantial effort to construct and are highly specific to a given chemical composition and application. A significant limitation of machine learning models is the use of element-specific features, leading to models that scale poorly with the number of elements. This work introduces the Gaussian multipole (GMP) featurization scheme that utilizes physically-relevant multipole expansions of the electron density around atoms to yield feature vectors that interpolate between element types and have a fixed dimension regardless of the number of elements present. We combine GMP with neural networks to directly compare it to the widely used Behler-Parinello symmetry functions for the MD17 dataset, revealing that it exhibits improved accuracy and computational efficiency. Further, we demonstrate that GMP-based models can achieve chemical accuracy for the QM9 dataset, and their accuracy remains reasonable even when extrapolating to new elements. Finally, we test GMP-based models for the Open Catalysis Project (OCP) dataset, revealing comparable performance to graph convolutional deep learning models. The results indicate that this featurization scheme fills a critical gap in the construction of efficient and transferable machine-learned force fields.


title: Investigating the environmental dependence of ultralight scalar dark matter with atom interferometers

abstract: We study the environmental dependence of ultralight scalar dark matter (DM) with linear interactions to the standard model particles. The solution to the DM field turns out to be a sum of the cosmic harmonic oscillation term and the local exponential fluctuation term. The amplitude of the first term depends on the local DM density and the mass of the DM field. The second term is induced by the local distribution of matter, such as the Earth. And it depends not only on the mass of the Earth, but also the density of the Earth. Then, we compute the phase shift induced by the DM field in atom interferometers (AIs), through solving the trajectories of atoms. Especially, the AI signal for the violation of weak equivalence principle (WEP) caused by the DM field is calculated. Depending on the values of the DM coupling parameters, contributions to the WEP violation from the first and second terms of the DM field can be either comparable or one larger than the other. Finally, we give some constraints to DM coupling parameters using results from the terrestrial atomic WEP tests.