title: Synthetic photometry of OB star clusters with stochastically sampled IMFs: analysis of models and HST observations

abstract: We present a pilot library of synthetic NUV, U, B, V, and I photometry of star clusters with stochastically sampled IMFs and ionized gas for initial masses, $M_i=10^3$, $10^4$, and $10^5$ M$_{\odot}$; $t=1$, 3, 4, and 8 Myr; $Z=0.014$ and $Z=0.002$; and log(U$_{\rm S}$) =-2 and -3. We compare the library with predictions from deterministic models and observations of isolated low-mass ($<10^4$ M$_{\odot}$) star clusters with co-spatial compact H\2 regions. The clusters are located in NGC 7793, one of the nearest galaxies observed as part of the \hst LEGUS and H$\alpha$-LEGUS surveys. 1) For model magnitudes that only account for the stars: a) the residual |deterministic mag - median stochastic mag| can be $\ge0.5$ mag, even for $M_i=10^5$ M$_{\odot}$; and b) the largest spread in stochastic magnitudes occurs when Wolf-Rayet stars are present. 2) For $M_i=10^5$ M$_{\odot}$: a) the median stochastic mag with gas can be $>$1.0 mag more luminous than the median stochastic magnitude without gas; and b) nebular emission lines can contribute with $>50\%$ and $>30\%$ to the total emission in the V and I bands, respectively. 3) Age-dating OB-star clusters via deterministic tracks in the U-B vs. V-I plane is highly uncertain at $Z=0.014$ for $M_i\sim10^3$ M$_{\odot}$ and $Z=0.002$ for $M_i\sim10^3-10^5$ M$_{\odot}$. 4) For low-mass clusters, the V-band extinction derived with stochastic models significantly depends on the value of log(U$_{\rm S}$). 5) The youngest clusters tend to have higher extinction. 6) The majority of clusters have multi-peaked age PDFs. 7) Finally, we discuss the importance of characterising the true variance in the number of stars per mass bin in nature.


title: Perspective-taking to Reduce Affective Polarization on Social Media

abstract: The intensification of affective polarization worldwide has raised new questions about how social media platforms might be further fracturing an already-divided public sphere. As opposed to ideological polarization, affective polarization is defined less by divergent policy preferences and more by strong negative emotions towards opposing political groups, and thus arguably poses a formidable threat to rational democratic discourse. We explore if prompting perspective-taking on social media platforms can help enhance empathy between opposing groups as a first step towards reducing affective polarization. Specifically, we deploy a randomized field experiment through a browser extension to 1,611 participants on Twitter, which enables participants to randomly replace their feeds with those belonging to accounts whose political views either agree with or diverge from their own. We find that simply exposing participants to "outgroup" feeds enhances engagement, but not an understanding of why others hold their political views. On the other hand, framing the experience in familiar, empathic terms by prompting participants to recall a disagreement with a friend does not affect engagement, but does increase their ability to understand opposing views. Our findings illustrate how social media platforms might take simple steps that align with business objectives to reduce affective polarization.


title: Learning to Coordinate in Multi-Agent Systems: A Coordinated Actor-Critic Algorithm and Finite-Time Guarantees

abstract: Multi-agent reinforcement learning (MARL) has attracted much research attention recently. However, unlike its single-agent counterpart, many theoretical and algorithmic aspects of MARL have not been well-understood. In this paper, we study the emergence of coordinated behavior by autonomous agents using an actor-critic (AC) algorithm. Specifically, we propose and analyze a class of coordinated actor-critic algorithms (CAC) in which individually parametrized policies have a {\it shared} part (which is jointly optimized among all agents) and a {\it personalized} part (which is only locally optimized). Such kind of {\it partially personalized} policy allows agents to learn to coordinate by leveraging peers' past experience and adapt to individual tasks. The flexibility in our design allows the proposed MARL-CAC algorithm to be used in a {\it fully decentralized} setting, where the agents can only communicate with their neighbors, as well as a {\it federated} setting, where the agents occasionally communicate with a server while optimizing their (partially personalized) local models. Theoretically, we show that under some standard regularity assumptions, the proposed MARL-CAC algorithm requires $\mathcal{O}(\epsilon^{-\frac{5}{2}})$ samples to achieve an $\epsilon$-stationary solution (defined as the solution whose squared norm of the gradient of the objective function is less than $\epsilon$). To the best of our knowledge, this work provides the first finite-sample guarantee for decentralized AC algorithm with partially personalized policies.


title: GCN-SE: Attention as Explainability for Node Classification in Dynamic Graphs

abstract: Graph Convolutional Networks (GCNs) are a popular method from graph representation learning that have proved effective for tasks like node classification tasks. Although typical GCN models focus on classifying nodes within a static graph, several recent variants propose node classification in dynamic graphs whose topologies and node attributes change over time, e.g., social networks with dynamic relationships, or literature citation networks with changing co-authorships. These works, however, do not fully address the challenge of flexibly assigning different importance to snapshots of the graph at different times, which depending on the graph dynamics may have more or less predictive power on the labels. We address this challenge by proposing a new method, GCN-SE, that attaches a set of learnable attention weights to graph snapshots at different times, inspired by Squeeze and Excitation Net (SE-Net). We show that GCN-SE outperforms previously proposed node classification methods on a variety of graph datasets. To verify the effectiveness of the attention weight in determining the importance of different graph snapshots, we adapt perturbation-based methods from the field of explainable machine learning to graphical settings and evaluate the correlation between the attention weights learned by GCN-SE and the importance of different snapshots over time. These experiments demonstrate that GCN-SE can in fact identify different snapshots' predictive power for dynamic node classification.


title: Flavored Gauge-Mediated Supersymmetry Breaking Models with Discrete Non-Abelian Symmetries

abstract: We investigate flavored gauge mediation models in which the Higgs and messenger doublets are embedded in multiplets of the discrete non-Abelian symmetry S3. In these theories, the S3 symmetry correlates the flavor structure of the quark and lepton Yukawa couplings with the structure of the messenger Yukawa couplings that contribute to the soft supersymmetry breaking mass parameters. We provide a systematic exploration of possible scenarios within this framework that can accommodate hierarchical quark and charged lepton masses, and examine the resulting phenomenological implications in each case. We find a heavier spectrum for the superpartner masses compared to flavored gauge mediation models controlled by Abelian symmetries, which can be directly traced back to the need in our scenarios for two vectorlike pairs of messenger fields for viable electroweak symmetry breaking.


title: The temperature and free energy of multi-black hole systems

abstract: In the present paper, we compute the Euclidean action of a generic system consisting of $N$ arbitrary Kerr-Newman black holes located on the symmetry axis and separated from each other by massless struts. This allows us to introduce the {\it Hawking average temperature} (HAT) $\hat T$ of the multi-black hole system via the condition of vanishing the entire set of terms involving the singular horizons due to periodic time, and the resulting formula for this temperature contains solely the surface gravities $\kappa_i$ and horizon areas $A^H_i$ of the black hole constituents. We also show that the corresponding expression for the {\it free energy} of the system defined by $\hat T$ is consistent with the first law of thermodynamics and Smarr mass relations.


title: A Time-Optimized Content Creation Workflow for Remote Teaching

abstract: We describe our workflow to create an engaging remote learning experience for a university course, while minimizing the post-production time of the educators. We make use of ubiquitous and commonly free services and platforms, so that our workflow is inclusive for all educators and provides polished experiences for students. Our learning materials provide for each lecture: 1) a recorded video, uploaded on YouTube, with exact slide timestamp indices, which enables an enhanced navigation UI; and 2) a high-quality flow-text automated transcript of the narration with proper punctuation and capitalization, improved with a student participation workflow on GitHub. All these results could be created by hand in a time consuming and costly way. However, this would generally exceed the time available for creating course materials. Our main contribution is to automate the transformation and post-production between raw narrated slides and our published materials with a custom toolchain. Furthermore, we describe our complete workflow: from content creation to transformation and distribution. Our students gave us overwhelmingly positive feedback and especially liked our use of ubiquitous platforms. The most used feature was YouTube's chapter UI enabled through our automatically generated timestamps. The majority of students, who started using the transcripts, continued to do so. Every single transcript was corrected by students, with an average word-change of 6%. We conclude with the positive feedback that our enhanced content formats are much appreciated and utilized. Important for educators is how our low overhead production workflow was sustainable throughout a busy semester.


title: PNJL model at zero temperature: three-flavor case

abstract: We propose a three-flavor version of the Polyakov-Nambu-Jona-Lasino (PNJL) model at zero temperature regime, by implementing a traced Polyakov loop ($\Phi$) dependence in the scalar, vector and 't Hooft channel strengths. We study the thermodynamics of this model, named as PNJL0, with special attention for the first order confinement/deconfinement phase transition for which $\Phi$ is the order parameter. For the symmetric quark matter case, an interesting feature observed is a strong reduction of the constituent strange quark mass ($M_s$) at the chemical potential related to point where deconfinement takes place. The emergence of $\Phi$ favors the restoration of chiral symmetry even for the strange quark. We also investigate the charge neutral system of quarks and leptons in weak equilibrium. As an application, we construct a hadron-quark phase transition with a density dependent hadronic model coupled to the SU(3) PNJL0 model. In this case, the quark side is composed by deconfined particles. This approach is used to determine mass-radius profiles compatible with recent data from the Neutron Star Interior Composition Explorer (NICER) mission.


title: Generalizing to New Domains by Mapping Natural Language to Lifted LTL

abstract: Recent work on using natural language to specify commands to robots has grounded that language to LTL. However, mapping natural language task specifications to LTL task specifications using language models require probability distributions over finite vocabulary. Existing state-of-the-art methods have extended this finite vocabulary to include unseen terms from the input sequence to improve output generalization. However, novel out-of-vocabulary atomic propositions cannot be generated using these methods. To overcome this, we introduce an intermediate contextual query representation which can be learned from single positive task specification examples, associating a contextual query with an LTL template. We demonstrate that this intermediate representation allows for generalization over unseen object references, assuming accurate groundings are available. We compare our method of mapping natural language task specifications to intermediate contextual queries against state-of-the-art CopyNet models capable of translating natural language to LTL, by evaluating whether correct LTL for manipulation and navigation task specifications can be output, and show that our method outperforms the CopyNet model on unseen object references. We demonstrate that the grounded LTL our method outputs can be used for planning in a simulated OO-MDP environment. Finally, we discuss some common failure modes encountered when translating natural language task specifications to grounded LTL.


title: A caster-wheel-aware MPC-based motion planner for mobile robotics

abstract: Differential drive mobile robots often use one or more caster wheels for balance. Caster wheels are appreciated for their ability to turn in any direction almost on the spot, allowing the robot to do the same and thereby greatly simplifying the motion planning and control. However, in aligning the caster wheels to the intended direction of motion they produce a so-called bore torque. As a result, additional motor torque is required to move the robot, which may in some cases exceed the motor capacity or compromise the motion planner's accuracy. Instead of taking a decoupled approach, where the navigation and disturbance rejection algorithms are separated, we propose to embed the caster wheel awareness into the motion planner. To do so, we present a caster-wheel-aware term that is compatible with MPC-based control methods, leveraging the existence of caster wheels in the motion planning stage. As a proof of concept, this term is combined with a a model-predictive trajectory tracking controller. Since this method requires knowledge of the caster wheel angle and rolling speed, an observer that estimates these states is also presented. The efficacy of the approach is shown in experiments on an intralogistics robot and compared against a decoupled bore-torque reduction approach and a caster-wheel agnostic controller. Moreover, the experiments show that the presented caster wheel estimator performs sufficiently well and therefore avoids the need for additional sensors.


title: Randomized Extended Kaczmarz is a Limit Point of Sketch-and-Project

abstract: The sketch-and-project (SAP) framework for solving systems of linear equations has unified the theory behind popular projective iterative methods such as randomized Kaczmarz, randomized coordinate descent, and variants thereof. The randomized extended Kaczmarz (REK) method is a popular extension of randomized Kaczmarz for solving inconsistent systems, which has not yet been shown to lie within the SAP framework. In this work we show that, in a certain sense, REK may be expressed as the limit point of a family of SAP methods, but we argue that it is unlikely that REK can be translated into a SAP method itself. We provide an extensive theoretical analysis of the family of methods comprising said limit, including convergence guarantees and further connections to REK. We follow this with an array of experiments demonstrating these methods and their connections in practice.


title: Nearest Subspace Search in The Signed Cumulative Distribution Transform Space for 1D Signal Classification

abstract: This paper presents a new method to classify 1D signals using the signed cumulative distribution transform (SCDT). The proposed method exploits certain linearization properties of the SCDT to render the problem easier to solve in the SCDT space. The method uses the nearest subspace search technique in the SCDT domain to provide a non-iterative, effective, and simple to implement classification algorithm. Experiments show that the proposed technique outperforms the state-of-the-art neural networks using a very low number of training samples and is also robust to out-of-distribution examples on simulated data. We also demonstrate the efficacy of the proposed technique in real-world applications by applying it to an ECG classification problem. The python code implementing the proposed classifier can be found in PyTransKit (https://github.com/rohdelab/PyTransKit).


title: Partial Variable Training for Efficient On-Device Federated Learning

abstract: This paper aims to address the major challenges of Federated Learning (FL) on edge devices: limited memory and expensive communication. We propose a novel method, called Partial Variable Training (PVT), that only trains a small subset of variables on edge devices to reduce memory usage and communication cost. With PVT, we show that network accuracy can be maintained by utilizing more local training steps and devices, which is favorable for FL involving a large population of devices. According to our experiments on two state-of-the-art neural networks for speech recognition and two different datasets, PVT can reduce memory usage by up to 1.9$\times$ and communication cost by up to 593$\times$ while attaining comparable accuracy when compared with full network training.


title: Tiebout sorting in online communities

abstract: This paper proposes a stylized, dynamic model to address the issue of sorting online. There are two large homogeneous groups of individuals. Everyone must choose between two online platforms, one of which has superior amenities (akin to having superior local public goods). Each individual enjoys interacting online with those from their own group but dislikes being on the same platform as those in the other group. Unlike a Tiebout model of residential sorting, both platforms have unlimited capacity so there are no constraints on cross-platform migration. It is clear how each group would like to sort themselves but, in the presence of the other type, only the two segregated outcomes are guaranteed to be equilibria. Integration on a platform can be supported in equilibrium as long as the platform is sufficiently desirable. If online integration of the two communities is a desired social outcome then the optimal policy is clear: make the preferred platform even more desirable. Revitalizing the inferior platform will never lead to integration and even increases the likelihood the segregation. Finally, integration is more elastic in response to an increase in platform amenities than to reductions in intolerance.


title: Comparison between time-independent and time-dependent quantum systems in the context of energy, Heisenberg uncertainty, average energy, force, average force and thermodynamic quantities

abstract: Exact solutions of time-dependent Schr\"odinger equation in presence of time-dependent potential is defined by point transformation and separation of variables. Energy and Heisenberg uncertainty relation are pursued for time-independent potential whereas average energy and Heisenberg uncertainty relation are defined for time-dependent potential. Forces acting on a fixed boundary wall as well as average force acting on moving boundary wall are presented along various trajectories. For high temperature, analytical forms of partition function and the corresponding thermodynamic quantities are derived following the Euler-Maclaurin summation formula over a finite as well as an infinite domain for accurate presentation. Three quantum systems are generated with the help of point transformation, separation of variables and super-symmetric quantum mechanics from one quantum system and the corresponding results are compared among all systems, where two of them are time-independent and another two are time-dependent.


title: TSK Fuzzy System Towards Few Labeled Incomplete Multi-View Data Classification

abstract: Data collected by multiple methods or from multiple sources is called multi-view data. To make full use of the multi-view data, multi-view learning plays an increasingly important role. Traditional multi-view learning methods rely on a large number of labeled and completed multi-view data. However, it is expensive and time-consuming to obtain a large number of labeled multi-view data in real-world applications. Moreover, multi-view data is often incomplete because of data collection failures, self-deficiency, or other reasons. Therefore, we may have to face the problem of fewer labeled and incomplete multi-view data in real application scenarios. In this paper, a transductive semi-supervised incomplete multi-view TSK fuzzy system modeling method (SSIMV_TSK) is proposed to address these challenges. First, in order to alleviate the dependency on labeled data and keep the model interpretable, the proposed method integrates missing view imputation, pseudo label learning of unlabeled data, and fuzzy system modeling into a single process to yield a model with interpretable fuzzy rules. Then, two new mechanisms, i.e. the bidirectional structural preservation of instance and label, as well as the adaptive multiple alignment collaborative learning, are proposed to improve the robustness of the model. The proposed method has the following distinctive characteristics: 1) it can deal with the incomplete and few labeled multi-view data simultaneously; 2) it integrates the missing view imputation and model learning as a single process, which is more efficient than the traditional two-step strategy; 3) attributed to the interpretable fuzzy inference rules, this method is more interpretable. Experimental results on real datasets show that the proposed method significantly outperforms the state-of-the-art methods.


title: Heat and Economic Preferences

abstract: The empirical evidence suggests that key accumulation decisions and risky choices associated with economic development depend, at least in part, on economic preferences such as willingness to take risk and patience. This paper studies whether temperature could be one of the potential channels that influences such economic preferences. Using data from the Indonesia Family Life Survey and NASAs Modern Era Retrospective Analysis for Research and Applications data we exploit quasi exogenous variations in outdoor temperatures caused by the random allocation of survey dates. This approach allows us to estimate the effects of temperature on elicited measures of risk aversion, rational choice violations, and impatience. We then explore three possible mechanisms behind this relationship, cognition, sleep, and mood. Our findings show that higher temperatures lead to significantly increased rational choice violations and impatience, but do not significantly increase risk aversion. These effects are mainly driven by night time temperatures on the day prior to the survey and less so by temperatures on the day of the survey. This impact is quasi linear and increasing when midnight outdoor temperatures are above 22C. The evidence shows that night time temperatures significantly deplete cognitive functioning, mathematical skills in particular. Based on these findings we posit that heat induced night time disturbances cause stress on critical parts of the brain, which then manifest in significantly lower cognitive functions that are critical for individuals to perform economically rational decision making.


title: Cluster expansion and resurgence in Polyakov model

abstract: In Polyakov model, a non-perturbative mass gap is formed at leading order semi-classics by instanton effects. By using the notions of critical points at infinity, cluster expansion and Lefschetz thimbles, we show that a third order effect in semi-classics gives an imaginary ambiguous contribution to mass gap, which is supposed to be real and unambiguous. This is troublesome for the original analysis, and it is difficult to resolve this issue directly in QFT. However, we find a new compactification of Polyakov model to quantum mechanics, by using a background 't Hooft flux (or coupling to TQFT). The compactification has the merit of remembering the monopole-instantons of the full QFT within Born-Oppenheimer (BO) approximation, while the periodic compactification does not. In QM, we prove the resurgent cancellation of the ambiguity in 3-instanton sector against ambiguity in the Borel resummation of the perturbation theory around 1-instanton. Assuming that this result holds in QFT, we provide a large-order asymptotics of perturbation theory around perturbative vacuum and instanton.


title: Virtual Knot Groups

abstract: For a knot diagram $K$, the classical knot group $\pi_1(K)$ is a free group modulo relations determined by Wirtinger-type relations on the classical crossings. The classical knot group is invariant under the Reidemeister moves. In this paper, we define a set of quotient groups associated to a knot diagram $K$. These quotient groups are invariant under the Reidemeister moves and the set includes the extended knot groups defined by Boden et al and Silver and Williams.


title: Signal Processing on Cell Complexes

abstract: The processing of signals supported on non-Euclidean domains has attracted large interest recently. Thus far, such non-Euclidean domains have been abstracted primarily as graphs with signals supported on the nodes, though the processing of signals on more general structures such as simplicial complexes has also been considered. In this paper, we give an introduction to signal processing on (abstract) regular cell complexes, which provide a unifying framework encompassing graphs, simplicial complexes, cubical complexes and various meshes as special cases. We discuss how appropriate Hodge Laplacians for these cell complexes can be derived. These Hodge Laplacians enable the construction of convolutional filters, which can be employed in linear filtering and non-linear filtering via neural networks defined on cell complexes.


title: Spirals in galaxies

abstract: Spirals in galaxies have long been thought to be caused by gravitational instability in the stellar component of the disk, but discerning the precise mechanism had proved elusive. Tidal interactions, and perhaps bars, may provoke some spiral responses, but spirals in many galaxies must be self-excited. We survey the relevant observational data and aspects of disk dynamical theory. The origin of the recurring spiral patterns in simulations of isolated disk galaxies has recently become clear and it is likely that the mechanism is the same in real galaxies, although evidence to confirm this supposition is hard to obtain. As transient spiral activity increases random motion, the patterns must fade over time unless the disk also contains a dissipative gas component. Continuing spiral activity alters the structure of the disks in other ways: reducing metallicity gradients and flattening rotation curves are two of the most significant. The overwhelming majority of spirals in galaxies have two- or three-fold rotational symmetry, indicating that the cool, thin disk component is massive. Spirals in simulations of halo-dominated disks instead manifest many arms, and consequently do not capture the expected full spiral-driven evolution. We conclude by identifying areas where further work is needed.


title: Polynomials over structured grids

abstract: We study multivariate polynomials over `structured' grids. We begin by proposing an interpretation as to what it means for a finite subset of a field to be structured; we do so by means of a numerical parameter, the nullity. We then extend several results--notably, the Combinatorial Nullstellensatz and the Coefficient Theorem--to polynomials over structured grids. The main point is that the structure of a grid allows the degree constraints on polynomials to be relaxed.


title: DFT study of tunable electronic, magnetic, thermal, and optical properties of a Ga$_2$Si$_6$ monolayer

abstract: The electrical, magnetic, thermal and optical characteristics of Gallium (Ga) doped silicene are investigated using density functional theory (DFT). The effect of doping is studied by tuning dopant concentrations as well as examining varied doping distances, and atomic dopant interactions for the same substitutional doping concentration. The results indicate that the Ga atoms alter the band structure and the band gap in the silicene monolayer at various concentrations, which can be referred back to to the repulsive interaction of Ga-Ga atoms. The band gap is determined by the interaction strength of the Ga-Ga atoms, the Coulomb repulsive force, and it does not always widen as doping concentration increases. In addition, our spin-polarized DFT calculations show that these monolayers behave like nonmagnetic semiconductors, exhibiting symmetric spin-up and spin-down channels. The repulsive interaction between the Ga atoms causes a symmetry breaking of the monolayers. As a consequence, a Ga dopant can open the band gap, leading to better thermoelectric properties such as the Seebeck coefficient and the figure of merit, as well as an increase in the optical response. As a result of our estimates, Ga doped silicene monolayers could be advantageous in thermoelectric and optoelectronic devices.


title: Synchronization for Networks of Globally Coupled Maps in the Thermodynamic Limit

abstract: We study a network of finitely many interacting clusters where each cluster is a collection of globally coupled circle maps in the thermodynamic (or mean field) limit. The state of each cluster is described by a probability measure, and its evolution is given by a self-consistent transfer operator. A cluster is synchronized if its state is a Dirac measure. We provide sufficient conditions for all clusters to synchronize and we describe setups where the conditions are met thanks to the uncoupled dynamics and/or the (diffusive) nature of the coupling. We also give sufficient conditions for partially synchronized states to arise -- i.e. states where only a subset of the clusters is synchronized -- due to the forcing of a group of cluster on the rest of the network. Lastly, we use this framework to show emergence and stability of chimera states for these systems.


title: Towards a Principled Approach for Dynamic Analysis of Android's Middleware

abstract: The Android middleware, in particular the so-called systemserver, is a crucial and central component to Android's security and robustness. To understand whether the systemserver provides the demanded security properties, it has to be thoroughly tested and analyzed. A dedicated line of research focuses exclusively on this task. While static analysis builds on established tools, dynamic testing approaches lack a common foundation, which prevents the community from comparing, reproducing, or even re-using existing results from related work. This raises questions about whether the underlying approach of any proposed solution is the only possible or optimal one, if it can be re-used as a building block for future analyses, or whether results generalize. In this work, we argue that in order to steer away from incompatible custom toolchains and towards having comparable analyses with reproducible results, a more principled approach to dynamically analyzing the Android system is required. As an important first step in this direction, we propose a unified dynamic analysis platform that provides re-usable solutions for common challenges as the building blocks for future analyses and allows to compare different approaches under the same assumptions.


title: Charge, bond, and pair density wave orders in a strongly correlated system

abstract: The coexistence of multiple quasi-degenerate orders is the hallmark of the strongly correlated materials. Experiments often reveal several spatially modulated orders in the underdoped cuprates. This has come to the forefront with the possible detection of the pair density wave states. However, microscopic calculations often struggle to stabilize such spatially modulating orders as the ground state in the strong correlation limit. This work uses the $t-t^\prime-J$-model with an additional nearest-neighbor repulsion to stabilize spatially oscillating charge, bond, and pairing orders in the underdoped regime. We employ the standard Gutzwiller approach while treating the inhomogeneity for the spatial orders using the self-consistent Hartree-Fock-Bogoliubov methodology. Our calculations reveal that unidirectional bond density states coexisting with charge and pairing modulations can have lower energy than the uniform superconducting state over an extensive doping range. These modulating states vanish monotonically as the modulation wavevector becomes shorter with increased dopings. The finite momentum orders melt upon increasing doping to a vestigial nematic state which breaks the rotational symmetry of the system. The spatial features of the ground state at each doping reveal multiple wavevectors, which potentially drives the incommensuration of charge orders. Interestingly, the spatially modulating states are absent when the strong correlations criteria are relaxed, suggesting that the removal of double occupancy aids the stabilization of density wave orders.


title: Neural Architecture Search for Efficient Uncalibrated Deep Photometric Stereo

abstract: We present an automated machine learning approach for uncalibrated photometric stereo (PS). Our work aims at discovering lightweight and computationally efficient PS neural networks with excellent surface normal accuracy. Unlike previous uncalibrated deep PS networks, which are handcrafted and carefully tuned, we leverage differentiable neural architecture search (NAS) strategy to find uncalibrated PS architecture automatically. We begin by defining a discrete search space for a light calibration network and a normal estimation network, respectively. We then perform a continuous relaxation of this search space and present a gradient-based optimization strategy to find an efficient light calibration and normal estimation network. Directly applying the NAS methodology to uncalibrated PS is not straightforward as certain task-specific constraints must be satisfied, which we impose explicitly. Moreover, we search for and train the two networks separately to account for the Generalized Bas-Relief (GBR) ambiguity. Extensive experiments on the DiLiGenT dataset show that the automatically searched neural architectures performance compares favorably with the state-of-the-art uncalibrated PS methods while having a lower memory footprint.


title: Review of Kernel Learning for Intra-Hour Solar Forecasting with Infrared Sky Images and Cloud Dynamic Feature Extraction

abstract: The uncertainty of the energy generated by photovoltaic systems incurs an additional cost for a guaranteed, reliable supply of energy (i.e., energy storage). This investigation aims to decrease the additional cost by introducing probabilistic multi-task intra-hour solar forecasting (feasible in real time applications) to increase the penetration of photovoltaic systems in power grids. The direction of moving clouds is estimated in consecutive sequences of sky images by extracting features of cloud dynamics with the objective of forecasting the global solar irradiance that reaches photovoltaic systems. The sky images are acquired using a low-cost infrared sky imager mounted on a solar tracker. The solar forecasting algorithm is based on kernel learning methods, and uses the clear sky index as predictor and features extracted from clouds as feature vectors. The proposed solar forecasting algorithm achieved 16.45\% forecasting skill 8 minutes ahead with a resolution of 15 seconds. In contrast, previous work reached 15.4\% forecasting skill with the resolution of 1 minute. Therefore, this solar forecasting algorithm increases the performances with respect to the state-of-the-art, providing grid operators with the capability of managing the inherent uncertainties of power grids with a high penetration of photovoltaic systems.


title: Photoemission spectral functions from the three-body Green's function

abstract: We present an original strategy for the calculation of direct and inverse photo-emission spectra from first principles. The main goal is to go beyond the standard Green's function approaches, such as the $GW$ method, in order to find a good description not only of the quasiparticles but also of the satellite structures, which are of particular importance in strongly correlated materials. To this end we use as a key quantity the three-body Green's function, or, more precisely, its hole-hole-electron and electron-electron-hole parts, and we show how the one-body Green's function, and hence the corresponding spectral function, can be retrieved from it. We show that, contrary to the one-body Green's function, information about satellites is already present in the non-interacting three-body Green's function. Therefore, simple approximations to the three-body self-energy, which is defined by the Dyson equation for the three-body Green's function and which contains many-body effects, can still yield accurate spectral functions. In particular, the self-energy can be chosen to be static which could simplify a self-consistent solution of the Dyson equation. We give a proof of principle of our strategy by applying it to the Hubbard dimer, for which the exact self-energy is available.


title: Magnetoelastic Coupling in Hole-doped Two-dimensional \b{eta}-PbO

abstract: The realization of intertwined ferroelasticity and ferromagnetism in two-dimensional (2D) lattices is of great interest for broad nanoscale applications, but still remains a remarkable challenge. Here, we propose an alternative approach to realize the strongly coupled ferromagnetism and ferroelasticity by carrier doping. We demonstrate that prototypical 2D \b{eta}-PbO is dynamically, thermally and mechanically stable. Under hole doping, 2D \b{eta}-PbO possesses ferromagnetism and ferroelasticity simultaneously. Moreover, the robustness of ferromagnetic and ferroelastic orders are doping tunable. In particular, 2D \b{eta}-PbO features in-plane easy magnetization axis that is coupled with lattice direction, enabling ferroelastic manipulation of spin direction. Our work highlights a new direction for 2D magnetoelastic research and enables the possibility for multifunctional devices.


title: Inferring supply networks from mobile phone data to estimate the resilience of a national economy

abstract: National economies rest on networks of millions of customer-supplier relations. Some companies -- in the case of their default -- can trigger significant cascades of shock in the supply-chain network and are thus systemically risky. Up to now, systemic risk of individual companies was practically not quantifiable, due to the unavailability of firm-level transaction data. So far, economic shocks are typically studied in the framework of input-output analysis on the industry-level that can't relate risk to individual firms. Exact firm-level supply networks based on tax or payment data exist only for very few countries. Here we explore to what extent telecommunication data can be used as an inexpensive, easily available, and real-time alternative to reconstruct national supply networks on the firm-level. We find that the conditional probability of correctly identifying a true customer-supplier link -- given a communication link exists -- is about 90%. This quality level allows us to reliably estimate a systemic risk profile of an entire country that serves as a proxy for the resilience of its economy. In particular, we are able to identify the high systemic risk companies. We find that 65 firms have the potential to trigger large cascades of disruption in production chains that could cause severe damages in the economy. We verify that the topological features of the inter-firm communication network are highly similar to national production networks with exact firm-level interactions.


title: Parameterizing Activation Functions for Adversarial Robustness

abstract: Deep neural networks are known to be vulnerable to adversarially perturbed inputs. A commonly used defense is adversarial training, whose performance is influenced by model capacity. While previous works have studied the impact of varying model width and depth on robustness, the impact of increasing capacity by using learnable parametric activation functions (PAFs) has not been studied. We study how using learnable PAFs can improve robustness in conjunction with adversarial training. We first ask the question: how should we incorporate parameters into activation functions to improve robustness? To address this, we analyze the direct impact of activation shape on robustness through PAFs and observe that activation shapes with positive outputs on negative inputs and with high finite curvature can increase robustness. We combine these properties to create a new PAF, which we call Parametric Shifted Sigmoidal Linear Unit (PSSiLU). We then combine PAFs (including PReLU, PSoftplus and PSSiLU) with adversarial training and analyze robust performance. We find that PAFs optimize towards activation shape properties found to directly affect robustness. Additionally, we find that while introducing only 1-2 learnable parameters into the network, smooth PAFs can significantly increase robustness over ReLU. For instance, when trained on CIFAR-10 with additional synthetic data, PSSiLU improves robust accuracy by 4.54% over ReLU on ResNet-18 and 2.69% over ReLU on WRN-28-10 in the $\ell_{\infty}$ threat model while adding only 2 additional parameters into the network architecture. The PSSiLU WRN-28-10 model achieves 61.96% AutoAttack accuracy, improving over the state-of-the-art robust accuracy on RobustBench (Croce et al., 2020).


title: Subnetwork Constraints for Tighter Upper Bounds and Exact Solution of the Clique Partitioning Problem

abstract: We consider a variant of the clustering problem for a complete weighted graph. The aim is to partition the nodes into clusters maximizing the sum of the edge weights within the clusters. This problem is known as the clique partitioning problem, being NP-hard in the general case of having edge weights of different signs. We propose a new method of estimating an upper bound of the objective function that we combine with the classical branch-and-bound technique to find the exact solution. We evaluate our approach on a broad range of random graphs and real-world networks. The proposed approach provided tighter upper bounds and achieved significant convergence speed improvements compared to known alternative methods.


title: The TESS-Keck Survey. VI. Two Eccentric sub-Neptunes Orbiting HIP-97166

abstract: We report the discovery of HIP-97166b (TOI-1255b), a transiting sub-Neptune on a 10.3-day orbit around a K0 dwarf 68 pc from Earth. This planet was identified in a systematic search of TESS Objects of Interest for planets with eccentric orbits, based on a mismatch between the observed transit duration and the expected duration for a circular orbit. We confirmed the planetary nature of HIP-97166b with ground-based radial velocity measurements and measured a mass of $M_{b} =$ 20 $\pm$ 2 $M_\bigoplus$ along with a radius of $R_{b} =$ 2.7 $\pm$ 0.1 $R_\bigoplus$ from photometry. We detected an additional non-transiting planetary companion with $M_{c}$ sin$i =$ 10 $\pm$ 2 $M_\bigoplus$ on a 16.8-day orbit. While the short transit duration of the inner planet initially suggested a high eccentricity, a joint RV-photometry analysis revealed a high impact parameter $b = 0.84 \pm 0.03$ and a moderate eccentricity. Modeling the dynamics with the condition that the system remain stable over $>$10$^5$ orbits yielded eccentricity constraints $e_b = 0.16 \pm 0.03$ and $e_c < 0.25$. The eccentricity we find for planet b is above average for the small population of sub-Neptunes with well-measured eccentricities. We explored the plausible formation pathways of this system, proposing an early instability and merger event to explain the high density of the inner planet at $5.3 \pm 0.9$ g/cc as well as its moderate eccentricity and proximity to a 5:3 mean-motion resonance.


title: Two-Dimensional Transition Metal Silicate Formed on Ru (0001) by Hydrogenation

abstract: Bottom-up synthesis of two-dimensional transition-metal silicates has been challenging due to strong overlayer-substrate interactions, which prevents the exfoliation of the overlayer. Here, using density functional theory calculations, we systematically investigate the hydrogenation of the overlayer as a way to decrease the substrate and overlayer interactions. Using the Fe$_2$Si$_2$O$_8\cdot$O/Ru(0001) structure as our starting point Wlodarczyk et. al. [1], we study hydrogenation levels up to Fe$_2$Si$_2$O$_9$H$_4\cdot$/Ru(0001). Structural and thermodynamical properties are studied at different hydrogenation levels to show under which conditions, the exfoliation can be feasible. Simulated core-level shifts show that Fe is primarily in 3+ state through the hydrogenation of Fe$_2$Si$_2$O$_8\cdot$O/Ru(0001). Simulated reflection adsorption infrared spectroscopy (RAIRS) yield distinctive shifts in vibrational properties with increasing hydrogenation which can guide experiments. [1] R. Wlodarczyk et al., "Atomic structure of an ultrathin Fe-silicate film grown on a metal: A monolayer of clay?", J. Am. Chem. Soc., 135, 19222, (2013).


title: Notes on kAExp(pol) problems for deterministic machines

abstract: The complexity of several logics, such as Presburger arithmetic, dependence logics and ambient logics, can only be characterised in terms of alternating Turing machines. Despite quite natural, the presence of alternation can sometimes cause neat ideas to be obfuscated inside heavy technical machinery. In these notes, we propose two problems on deterministic machines that can be used to prove lower bounds with respect to the computational class $k$AExp$_{\text{pol}}$, that is the class of all problems solvable by an alternating Turing machine running in $k$ exponential time and performing a polynomial amount of alternations, with respect to the input size. The first problem, called $k$AExp$_{\text{pol}}$-prenex TM problem, is a problem about deterministic Turing machines. The second problem, called the $k$-exp alternating multi-tiling problem, is analogous to the first one, but on tiling systems. Both problems are natural extensions of the TM alternation problem and the alternating multi-tiling problem proved AExp$_{\text{pol}}$-complete by L. Bozzelli, A. Molinari, A. Montanari and A. Peron in [GandALF, pp. 31-45, 2017]. The proofs presented in these notes follow the elegant exposition in A. Molinari's PhD thesis to extend these results from the case $k = 1$ to the case of arbitrary $k$.


title: Reeb Graph Metrics from the Ground Up

abstract: The Reeb graph has been utilized in various applications including the analysis of scalar fields. Recently, research has been focused on using topological signatures such as the Reeb graph to compare multiple scalar fields by defining distance metrics on the topological signatures themselves. Here we survey five existing metrics that have been defined on Reeb graphs: the bottleneck distance, the interleaving distance, functional distortion distance, the Reeb graph edit distance, and the universal edit distance. Our goal is to (1) provide definitions and concrete examples of these distances in order to develop the intuition of the reader, (2) visit previously proven results of stability, universality, and discriminativity, (3) identify and complete any remaining properties which have only been proven (or disproven) for a subset of these metrics, (4) expand the taxonomy of the bottleneck distance to better distinguish between variations which have been commonly miscited, and (5) reconcile the various definitions and requirements on the underlying spaces for these metrics to be defined and properties to be proven.


title: Wind-robust sound event detection and denoising for bioacoustics

abstract: Sound recordings are used in various ecological studies, including acoustic wildlife monitoring. Such surveys require automatic detection of target sound events. However, current detectors, especially those relying on band-limited energy, are severely impacted by wind. The rapid dynamics of this noise invalidate standard noise estimators, and no satisfactory method for dealing with it exists in bioacoustics, where simple training and generalization between conditions are important. We propose to estimate the transient noise level by fitting short-term spectrum models to a wavelet packet representation. This estimator is then combined with log-spectral subtraction to stabilize the background level. The resulting adjusted wavelet series can be analysed by standard energy detectors. We use real monitoring data to tune this workflow, and test it on two acoustic surveys of birds. Additionally, we show how the estimator can be incorporated in a denoising method to restore sound. The proposed noise-robust workflow greatly reduced the number of false alarms in the surveys, compared to unadjusted energy detection. As a result, the survey efficiency (precision of the estimated call density) improved for both species. Denoising was also more effective when using the short-term estimate, whereas standard wavelet shrinkage with a constant noise estimate struggled to remove the effects of wind. In contrast to existing methods, the proposed estimator can adjust for transient broadband noises without requiring additional hardware or extensive tuning to each species. It improved the detection workflow based on very little training data, making it particularly attractive for detection of rare species.


title: TCube: Domain-Agnostic Neural Time-series Narration

abstract: The task of generating rich and fluent narratives that aptly describe the characteristics, trends, and anomalies of time-series data is invaluable to the sciences (geology, meteorology, epidemiology) or finance (trades, stocks, or sales and inventory). The efforts for time-series narration hitherto are domain-specific and use predefined templates that offer consistency but lead to mechanical narratives. We present TCube (Time-series-to-text), a domain-agnostic neural framework for time-series narration, that couples the representation of essential time-series elements in the form of a dense knowledge graph and the translation of said knowledge graph into rich and fluent narratives through the transfer-learning capabilities of PLMs (Pre-trained Language Models). TCube's design primarily addresses the challenge that lies in building a neural framework in the complete paucity of annotated training data for time-series. The design incorporates knowledge graphs as an intermediary for the representation of essential time-series elements which can be linearized for textual translation. To the best of our knowledge, TCube is the first investigation of the use of neural strategies for time-series narration. Through extensive evaluations, we show that TCube can improve the lexical diversity of the generated narratives by up to 65.38% while still maintaining grammatical integrity. The practicality and deployability of TCube is further validated through an expert review (n=21) where 76.2% of participating experts wary of auto-generated narratives favored TCube as a deployable system for time-series narration due to its richer narratives. Our code-base, models, and datasets, with detailed instructions for reproducibility is publicly hosted at https://github.com/Mandar-Sharma/TCube.


title: On the AdS/CFT correspondence and quantum entanglement

abstract: String theory provides one of the most deepest insights into quantum gravity. Its single most central and profound result is the gauge/gravity duality, i.e. the emergence of gravity from gauge theory. The two examples of M(atrix)-theory and the AdS/CFT correspondence, together with the fundamental phenomena of quantum entanglement, are of paramount importance to many fundamental problems including the physics of black holes (in particular to the information loss paradox), the emergence of spacetime geometry and to the problem of the reconciliation of general relativity and quantum mechanics. In this article an account of the AdS/CFT correspondence and the role of quantum entanglement in the emergence of spacetime geometry using strictly the language of quantum field theory is put forward.


title: Real-time EEG-based Emotion Recognition using Discrete Wavelet Transforms on Full and Reduced Channel Signals

abstract: Real-time EEG-based Emotion Recognition (EEG-ER) with consumer-grade EEG devices involves classification of emotions using a reduced number of channels. These devices typically provide only four or five channels, unlike the high number of channels (32 or more) typically used in most current state-of-the-art research. In this work we propose to use Discrete Wavelet Transforms (DWT) to extract time-frequency domain features, and we use time-windows of a few seconds to perform EEG-ER classification. This technique can be used in real-time, as opposed to post-hoc on the full session data. We also apply baseline removal preprocessing, developed in prior research, to our proposed DWT Entropy and Energy features, which improves classification accuracy significantly. We consider two different classifier architectures, a 3D Convolutional Neural Network (3D CNN) and a Support Vector Machine (SVM). We evaluate both models on subject-independent and subject dependent setups to classify the Valence and Arousal dimensions of an individual's emotional state. We test them on both the full 32-channel data provided by the DEAP dataset, and also a reduced 5-channel extract of the same dataset. The SVM model performs best on all the presented scenarios, achieving an accuracy of 95.32% on Valence and 95.68% on Arousal for the full 32-channel subject-dependent case, beating prior real-time EEG-ER subject-dependent benchmarks. On the subject-independent case an accuracy of 80.70% on Valence and 81.41% on Arousal was also obtained. Reducing the input data to 5 channels only degrades the accuracy by an average of 3.54% across all scenarios, making this model appropriate for use with more accessible low-end EEG devices.


title: CAPITAL: Optimal Subgroup Identification via Constrained Policy Tree Search

abstract: Personalized medicine, a paradigm of medicine tailored to a patient's characteristics, is an increasingly attractive field in health care. An important goal of personalized medicine is to identify a subgroup of patients, based on baseline covariates, that benefits more from the targeted treatment than other comparative treatments. Most of the current subgroup identification methods only focus on obtaining a subgroup with an enhanced treatment effect without paying attention to subgroup size. Yet, a clinically meaningful subgroup learning approach should identify the maximum number of patients who can benefit from the better treatment. In this paper, we present an optimal subgroup selection rule (SSR) that maximizes the number of selected patients, and in the meantime, achieves the pre-specified clinically meaningful mean outcome, such as the average treatment effect. We derive two equivalent theoretical forms of the optimal SSR based on the contrast function that describes the treatment-covariates interaction in the outcome. We further propose a ConstrAined PolIcy Tree seArch aLgorithm (CAPITAL) to find the optimal SSR within the interpretable decision tree class. The proposed method is flexible to handle multiple constraints that penalize the inclusion of patients with negative treatment effects, and to address time to event data using the restricted mean survival time as the clinically interesting mean outcome. Extensive simulations, comparison studies, and real data applications are conducted to demonstrate the validity and utility of our method.


title: Evaporating Primordial Black Holes in Gamma Ray and Neutrino Telescopes

abstract: A primordial black hole in the last stages of evaporation and located in the local neighborhood can produce a detectable signal in gamma ray and neutrino telescopes. We re-evaluate the expected gamma ray and neutrino fluxes from these transient point events and discuss the consequences for existing constraints. For gamma rays we improve the current bounds by a factor of few, while for neutrinos we obtain significantly different results than the existing literature. The capability and advantages of neutrino telescopes in the search for primordial black holes is discussed thoroughly. The correlations of gamma ray and neutrino energy and time profiles will be promoted as a powerful tool in identifying the primordial black holes, in case of detection.


title: Searching for Replacement Classes

abstract: Software developers must often replace existing components in their systems to adapt to evolving environments or tooling. While traditional code search systems are effective at retrieving components with related functionality, it is much more challenging to retrieve components that can be used to directly replace existing functionality, as replacements must account for more fundamental program properties such as type compatibility. To address this problem, we introduce ClassFinder, a system which given a query class Q, and a search corpus S, returns a ranked subset of classes that can replace Q and its functionality. ClassFinder produces afield and method mapping between the classes that can provide useful hints to a developer and can be used to effectively refine the ranking of candidate replacement classes. Our technique leverages the complementary strengths of a distributed embeddings-based search and type-based analysis, using the former to prune down candidates for an optimization-based approach based on the latter. ClassFinder retrieves replacement classes, along with a type-aware field/method mapping between classes. We evaluate ClassFinder on a search space of ~600thousand open sourceJava classes. Querying ClassFinder with 24 known Java classes provided meaningful replacement classes and mappings, in many cases producing complete mappings with functionally identical replacement classes.


title: Steady-state accretion in magnetized protoplanetary disks

abstract: [abridged] We present a 1+1D global magnetically-driven disk accretion model that captures the essence of the MRI-driven accretion, without resorting to 3D global non-ideal MHD simulations. The gas dynamics is assumed to be solely controlled by the MRI and hydrodynamic instabilities. For given stellar and disk parameters, the Shakura-Sunyaev viscosity parameter $\alpha$ is determined self-consistently under the framework of viscously-driven accretion from detailed considerations of the MRI with non-ideal MHD effects (Ohmic resistivity and ambipolar diffusion), accounting for disk heating by stellar irradiation, non-thermal sources of ionization, and dust effects on the ionization chemistry. Additionally, the magnetic field strength is constrained and adopted to maximize the MRI activity. We demonstrate the use of our framework by investigating steady-state MRI-driven accretion in a fiducial protoplanetary disk model around a solar-type star. We find that the equilibrium solution displays no pressure maximum at the dead zone outer edge, except if a sufficient amount of dust particles have accumulated there before the disk reaches a steady-state accretion regime. Furthermore, the steady-state accretion solution describes a disk that displays a spatially extended long-lived inner disk gas reservoir (the dead zone) accreting a few $10^{-9}\, M_{\odot}.\rm{yr}^{-1}$. By conducting a detailed parameter study, we find that the extend to which the MRI can drive efficient accretion is primarily determined by the total disk gas mass, the representative grain size, the vertically-integrated dust-to-gas mass ratio, and the stellar X-ray luminosity. A self-consistent time-dependent coupling between gas, dust, stellar evolution models and our general framework on million-year timescales is required to fully understand the formation of dead zones and their potential to trap dust particles.


title: K-theory of Jones polynomials

abstract: We recover the Jones polynomials of knots and links from the K-theory of a cluster C*-algebra of the sphere with two cusps. In particular, an interplay between the Chebyshev and Jones polynomials is studied.


title: Nonparametric Group Variable Selection with Multivariate Response for Connectome-Based Modeling of Cognitive Scores

abstract: In this article, we study association between the structural connectome and cognitive profiles using a multi-response nonparametric regression model.The cognitive profiles are measured in terms of seven age-adjusted cognitive test scores. The structural connectomes are represented by undirected graphs. The connectivity properties of these graphs are available in terms of the nodal attributes. A collection of nodal centralities together can encode different patterns of connections in the brain network. In this article, we consider nine such attributes for each brain region.These nodal graph metrics may naturally be grouped together for each node, motivating us to introduce group sparsity for feature selection. We propose Gaussian RBF-nets with a novel group sparsity inducing prior to model the unknown mean functions. The covariance structure of the multivariate response is characterized in terms of a linear factor modeling framework. For posterior computation, we develop an efficient Markov chain Monte Carlo sampling algorithm. We show that the proposed method performs much better than all its competitors. Applying our proposed method to a Human Connectome Project (HCP) dataset, we identify the important brain regions and nodal attributes for cognitive functioning, as well as identify interesting low-dimensional dependency structures among the cognition related test scores. Keywords: Factor model; Group variable selection; High-dimension; Human Connectome Project (HCP); Markov chain Monte Carlo (MCMC); Neural network; Nonparametric inference; Radial basis network; Spike-and-slab prior; Variable selection.


title: Non-classical light state transfer in $su(2)$ resonator networks

abstract: We use a normal mode approach to show full and partial state transfer in a class of coupled resonator networks with underlying $su(2)$ symmetry that includes the so-called $J_{x}$ photonic lattice. Our approach defines an auxiliary Hermitian coupling matrix describing the network that yields the normal modes of the system and its time evolution in terms of orthogonal polynomials. These results provide insight on the full quantum state reconstruction time in a general $su(2)$ network of any size and the full quantum transfer time in the $J_{x}$ network of size $4 n + 1$ with $n=1,2,3,\ldots$ In the latter, our approach shows that the Fock state probability distribution of the initial state is conserved but the amplitudes suffer a phase shift proportional to $\pi/2$ that results in partial quantum state transfer for any other network size.


title: Motivating Effort with Information about Future Rewards

abstract: This paper studies the optimal mechanism to motivate effort in a dynamic principal-agent model without transfers. An agent is engaged in a task with uncertain future rewards and can shirk irreversibly at any time. The principal knows the reward of the task and provides information to the agent over time in order to motivate effort. We derive the optimal information policy in closed form and thus identify two conditions, each of which guarantees that delayed disclosure is valuable. First, if the principal is impatient compared to the agent, she prefers the front-loaded effort schedule induced by delayed disclosure. In a stationary environment, delayed disclosure is beneficial if and only if the principal is less patient than the agent. Second, if the environment makes the agent become pessimistic over time in absence of any information disclosure, then providing delayed news can counteract this downward trend in the agent's belief and encourage the agent to work longer. Notably, the level of patience remains a crucial determinant of the optimal policy structure.


title: On the computational equivalence of co-NP refutations of a matrix being a P-matrix

abstract: A P-matrix is a square matrix $X$ such that all principal submatrices of $X$ have positive determinant. Such matrices appear naturally in instances of the linear complementarity problem, where these are precisely the matrices for which the corresponding linear complementarity problem has a unique solution for any input vector. Testing whether or not a square matrix is a P-matrix is co-NP complete, so while it is possible to exhibit polynomially-sized witnesses for the fact that a matrix is not a P-matrix, it is believed that there is no efficient way to prove that a given matrix is a P-matrix. We will show that several well known witnesses for the fact that a matrix is not a P-matrix are computationally equivalent, so that we are able to convert between them in polynomial time, answering a question raised in arXiv:1811.03841 .


title: A global convergence theory for deep ReLU implicit networks via over-parameterization

abstract: Implicit deep learning has received increasing attention recently due to the fact that it generalizes the recursive prediction rules of many commonly used neural network architectures. Its prediction rule is provided implicitly based on the solution of an equilibrium equation. Although a line of recent empirical studies has demonstrated its superior performances, the theoretical understanding of implicit neural networks is limited. In general, the equilibrium equation may not be well-posed during the training. As a result, there is no guarantee that a vanilla (stochastic) gradient descent (SGD) training nonlinear implicit neural networks can converge. This paper fills the gap by analyzing the gradient flow of Rectified Linear Unit (ReLU) activated implicit neural networks. For an $m$-width implicit neural network with ReLU activation and $n$ training samples, we show that a randomly initialized gradient descent converges to a global minimum at a linear rate for the square loss function if the implicit neural network is \textit{over-parameterized}. It is worth noting that, unlike existing works on the convergence of (S)GD on finite-layer over-parameterized neural networks, our convergence results hold for implicit neural networks, where the number of layers is \textit{infinite}.


title: Influence of long range forces on the transition states and dynamics of NaCl ion-pair dissociation in water

abstract: We study NaCl ion-pair dissociation in a dilute aqueous solution using computer simulations both for the full system with long range Coulomb interactions and for a well chosen reference system with short range intermolecular interactions. Analyzing results using concepts from Local Molecular Field (LMF) theory and the recently proposed AI-based analysis tool "State predictive information bottleneck" (SPIB) we show that the system with short range interactions can accurately reproduce the transition rate for the dissociation process, the dynamics for moving between the underlying metastable states, and the transition state ensemble. Contributions from long range interactions can be largely neglected for these processes because long range forces from the direct interionic Coulomb interactions are almost completely canceled ($>90\%$) by those from solvent interactions over the length scale where the transition takes place. Thus for this important monovalent ion-pair system, short range forces alone are able to capture detailed consequences of the collective solvent motion, allowing the use of physically suggestive and computationally efficient short range models for the disassociation event. We believe that the framework here should be applicable to disentangling mechanisms for more complex processes such as multivalent ion disassociation, where previous work has suggested that long range contributions may be more important.

categories: physics.chem-ph cond-mat.soft cond-mat.stat-mech physics.comp-ph


title: Growth, characterization and Chern insulator state in MnBi$_2$Te$_4$ via the chemical vapor transport method

abstract: As the first intrinsic antiferromagnetic topological insulator, MnBi$_2$Te$_4$ has provided a platform to investigate the interplay of band topology and magnetism as well as the emergent phenomena arising from such an interplay. Here we report the chemical-vapor-transport (CVT) growth and characterization of MnBi$_2$Te$_4$, as well as the observation of the field-induced quantized Hall conductance in 6-layer devices. Through comparative studies between our CVT-grown and flux-grown MnBi$_2$Te$_4$ via magnetic, transport, scanning tunneling microscopy, and angle-resolved photoemission spectroscopy measurements, we find that CVT-grown MnBi$_2$Te$_4$ is marked with higher Mn occupancy on the Mn site, slightly higher Mn$_{\rm{Bi}}$ antisites, smaller carrier concentration and a Fermi level closer to the Dirac point. Furthermore, a 6-layer device made from the CVT-grown sample shows by far the highest mobility of 2500 cm$^2$V$\cdot$s in MnBi$_2$Te$_4$ devices with the quantized Hall conductance appearing at 1.8 K and 8 T. Our study provides a new route to obtain high-quality single crystals of MnBi$_2$Te$_4$ that are promising to make superior devices and realize emergent phenomena, such as the layer Hall effect and quantized anomalous hall effect, etc.


title: Peak algebras in combinatorial Hopf algebras

abstract: The peak algebra is originally introduced by Stembridge using enriched $P$-partitions. Using the character theory by Aguiar-Bergeron-Sottile, the peak algebra is also the image of $\Theta$, the universal morphism between certain combinatorial Hopf algebras. We extend the notion of peak algebras and theta maps to shuffle, tensor, and symmetric algebras. As examples, we study the peak algebras of symmetric functions in non-commuting variables and the graded associated Hopf algebra on permutations. We also introduce a new shuffle basis of quasi-symmetric functions that its elements are the eigenfunctions of $\Theta$. Using this new basis, we show that the peak algebra is the space spanned by the set of shuffle basis elements indexed by compositions whose all parts are odd.


title: Learned Robust PCA: A Scalable Deep Unfolding Approach for High-Dimensional Outlier Detection

abstract: Robust principal component analysis (RPCA) is a critical tool in modern machine learning, which detects outliers in the task of low-rank matrix reconstruction. In this paper, we propose a scalable and learnable non-convex approach for high-dimensional RPCA problems, which we call Learned Robust PCA (LRPCA). LRPCA is highly efficient, and its free parameters can be effectively learned to optimize via deep unfolding. Moreover, we extend deep unfolding from finite iterations to infinite iterations via a novel feedforward-recurrent-mixed neural network model. We establish the recovery guarantee of LRPCA under mild assumptions for RPCA. Numerical experiments show that LRPCA outperforms the state-of-the-art RPCA algorithms, such as ScaledGD and AltProj, on both synthetic datasets and real-world applications.


title: GM-Livox: An Integrated Framework for Large-Scale Map Construction with Multiple Non-repetitive Scanning LiDARs

abstract: With the ability of providing direct and accurate enough range measurements, light detection and ranging (LiDAR) is playing an essential role in localization and detection for autonomous vehicles. Since single LiDAR suffers from hardware failure and performance degradation intermittently, we present a multi-LiDAR integration scheme in this article. Our framework tightly couples multiple non-repetitive scanning LiDARs with inertial, encoder, and global navigation satellite system (GNSS) into pose estimation and simultaneous global map generation. Primarily, we formulate a precise synchronization strategy to integrate isolated sensors, and the extracted feature points from separate LiDARs are merged into a single sweep. The fused scans are introduced to compute the scan-matching correspondences, which can be further refined by additional real-time kinematic (RTK) measurements. Based thereupon, we construct a factor graph along with the inertial preintegration result, estimated ground constraints, and RTK data. For the purpose of maintaining a restricted number of poses for estimation, we deploy a keyframe based sliding-window optimization strategy in our system. The real-time performance is guaranteed with multi-threaded computation, and extensive experiments are conducted in challenging scenarios. Experimental results show that the utilization of multiple LiDARs boosts the system performance in both robustness and accuracy.


title: Learning with Algorithmic Supervision via Continuous Relaxations

abstract: The integration of algorithmic components into neural architectures has gained increased attention recently, as it allows training neural networks with new forms of supervision such as ordering constraints or silhouettes instead of using ground truth labels. Many approaches in the field focus on the continuous relaxation of a specific task and show promising results in this context. But the focus on single tasks also limits the applicability of the proposed concepts to a narrow range of applications. In this work, we build on those ideas to propose an approach that allows to integrate algorithms into end-to-end trainable neural network architectures based on a general approximation of discrete conditions. To this end, we relax these conditions in control structures such as conditional statements, loops, and indexing, so that resulting algorithms are smoothly differentiable. To obtain meaningful gradients, each relevant variable is perturbed via logistic distributions and the expectation value under this perturbation is approximated. We evaluate the proposed continuous relaxation model on four challenging tasks and show that it can keep up with relaxations specifically designed for each individual task.


title: Spectroscopic evidence of high temperature superconductivity in VSe2

abstract: High-resolution angle-resolved photoemission experiments reveal subtle modifications of the surface electronic structure of VSe2. Most remarkably, we show that superconductivity can be induced in VSe2 by the right selection of substrate and growth parameters. Evidence for the superconducting state comes from the simultaneous detection of spectral kink, quasiparticle peak, Fermi gap, and their evolution with the temperature. The observation of Bogoliubov-like back-bending bands at low temperatures, signaling electron-hole pairing, further supports the presence of superconductivity in this system. The photoemission experiment also provides evidence for a formation of a pseudogap state at high temperatures, characterized by the progressive quenching of the quasiparticle peak feature coexisting with a persistent gap at the Fermi level, a behavior reminiscent of high-Tc superconductors. We attributed the origin of superconductivity in these VSe2 films to a modified Fermi surface combined with the formation of van Hove singularity points with a binding energy corresponding to the chemical potential of the system. Although Tc cannot be accurately determined from photoemission data, observations based on the survival temperature of the quasiparticle peak suggest that Tc could be as high as 100 ? 5 K and substantially higher than previous reports for any transition metal dichalcogenides.


title: Efficient Evaluation of Exponential and Gaussian Functions on a Quantum Computer

abstract: The exponential and Gaussian functions are among the most fundamental and important operations, appearing ubiquitously throughout all areas of science, engineering, and mathematics. Whereas formally, it is well-known that any function may in principle be realized on a quantum computer, in practice present-day algorithms tend to be very expensive. In this work, we present algorithms for evaluating exponential and Gaussian functions efficiently on quantum computers. The implementations require a (generally) small number of multiplications, which represent the overall computational bottleneck. For a specific, realistic NISQ application, the Toffoli count of the exponential function is found to be reduced from 15,690 down to 912, when compared against a state-of-the art competing method by H\"aner and coworkers [arXiv:1805.12445], under the most favorable conditions for each method. For the corresponding Gaussian function comparison, the Toffoli count is reduced from 19,090 down to 704. Space requirements are also quite modest, to the extent that the aforementioned NISQ application can be implemented with as few as 71 logical qubits. More generally, the methods presented here could also be equally well applied in a fault-tolerant context, using error-corrected multiplications, etc.


title: Bandgap tuning in kerfed metastrips under extreme deformation

abstract: The process of kerfing enables planar structures with the ability to undergo dramatic out-of-plane deformation in response to static loads. Starting from flat and stiff sheets, kerfing allows for the formation of a wide variety of unconventional free-form shapes, making the process especially attractive for architectural applications. In this work, we investigate numerically and experimentally the bandgap behavior of densely cut kerfed strips. Our study reveals a rich landscape of bandgaps that is predominantly ascribable to the activation of resonant sub-units within the kerf unit cells. We also document how the extreme deformability of the strips under twisting and bending loads, enhanced by the meandering cut pattern, can serve as a powerful bandgap tuning mechanism.


title: Defocus Map Estimation and Deblurring from a Single Dual-Pixel Image

abstract: We present a method that takes as input a single dual-pixel image, and simultaneously estimates the image's defocus map -- the amount of defocus blur at each pixel -- and recovers an all-in-focus image. Our method is inspired from recent works that leverage the dual-pixel sensors available in many consumer cameras to assist with autofocus, and use them for recovery of defocus maps or all-in-focus images. These prior works have solved the two recovery problems independently of each other, and often require large labeled datasets for supervised training. By contrast, we show that it is beneficial to treat these two closely-connected problems simultaneously. To this end, we set up an optimization problem that, by carefully modeling the optics of dual-pixel images, jointly solves both problems. We use data captured with a consumer smartphone camera to demonstrate that, after a one-time calibration step, our approach improves upon prior works for both defocus map estimation and blur removal, despite being entirely unsupervised.


title: The mean square radius of the neutron distribution in the relativistic and non-relativistic mean field models

abstract: It is investigated why the root mean square radius of the point neutron distribution is smaller by about 0.1 fm in non-relativistic mean field models than in relativistic ones. The difference is shown to stem from the different values of the product of the effective mass and the strength of the one-body potential in the two frameworks. The values of those quantities are constrained by the Hugenholtz-Van Hove theorem. The neutron skin is not a simple function of the symmetry potential, but depends on the nucleon effective mass.


title: From black-bounce to traversable wormhole, and beyond

abstract: Key results from the literature pertaining to a class of nonsingular black hole mimickers are explored. The family of candidate spacetimes is for now labelled the `black-bounce' family, stemming from the original so-called `Simpson--Visser' spacetime in static spherical symmetry. All model geometries are analysed through the lens of standard general relativity, are globally free from curvature singularities, pass all weak-field observational tests, and smoothly interpolate between regular black holes and traversable wormholes. The discourse is segregated along geometrical lines, with candidate spacetimes each belonging to one of: static spherical symmetry, spherical symmetry with dynamics, and stationary axisymmetry.


title: Collective coordinates for the hybrid model

abstract: In the present work, we carry out the study of scattering solitons for the anti-kink/kink and kink/anti-kink configurations. Furthermore, we can observe the same effects as those described by D. Bazeia et al.. We apply the collective coordinate approximation method to describe both scattering configurations and verify that just as happens in the polynomial models $\phi^4$ and $\phi^6$, the method has its limitations regarding the initial scattering speeds. In such a way that, for certain initial speeds, the solution of collective coordinates agrees with the fullsimulation, and for other speeds, there is a discrepancy in the solutions obtained by these two methods. We also noticed that, considering the hybrid model, the null-vector problem persists for both configurations, and when trying to fix it, a singularity is created in moduli-space as well as in $\phi^4$.


title: Central Values of Degree Six L-functions: The Case of Hilbert Modular Forms

abstract: In this paper we give a formula for the central value of the completed $L$-function $L(s,Sym^{2} g\times f)$, where $f$ and $g$ are Hilbert newforms, by explicitly computing the local integrals appearing in the refined Gan-Gross-Prasad conjecture for $SL_{2}\times\tilde{SL_{2}}$. We also work out the rationality of this value in some special cases and give a conjecture for the general case.


title: Orientable triangulable manifolds are essentially quasigroups

abstract: We introduce an $n$-dimensional analogue of the construction of tessellated surfaces from finite groups first described by Herman and Pakianathan. Our construction is functorial and associates to each $n$-ary alternating quasigroup both a smooth, flat Riemannian $n$-manifold which we dub the open serenation of the quasigroup in question, as well as a topological $n$-manifold (the serenation of the quasigroup) which is a subspace of the metric completion of the open serenation. We prove that every connected orientable smooth manifold is serene, in the sense that each such manifold is a component of the serenation of some quasigroup. We prove some basic results about the variety of alternating $n$-quasigroups and note connections between our construction, Latin hypercubes, and Johnson graphs.


title: BotNet Detection on Social Media

abstract: As our reliance on social media platforms and web services increase day by day, exploiters view these platforms as an opportunity to manipulate our thoughts ad actions. These platforms have become an open playground for social bot accounts. Social bots not only learn human conversations, manners, and presence but also manipulate public opinion, act as scammers, manipulate stock markets, and so on. There has been evidence of bots manipulating people's opinions and thoughts which can be a great threat to democracy. Identification and prevention of such campaigns that release or create these bots have become critical. Our goal in this paper is to leverage web mining techniques to help detect fake bots on social media platforms such as Twitter, thereby mitigating the spread of disinformation.


title: The colored Jones polynomials as vortex partition functions

abstract: We construct 3D $\mathcal{N}=2$ abelian gauge theories on $\mathbb{S}^2 \times \mathbb{S}^1$ labeled by knot diagrams whose K-theoretic vortex partition functions, each of which is a building block of twisted indices, give the colored Jones polynomials of knots in $\mathbb{S}^3$. The colored Jones polynomials are obtained as the Wilson loop expectation values along knots in $SU(2)$ Chern-Simons gauge theories on $\mathbb{S}^3$, and then our construction provides an explicit correspondence between 3D $\mathcal{N}=2$ abelian gauge theories and 3D $SU(2)$ Chern-Simons gauge theories. We verify, in particular, the applicability of our constructions to a class of tangle diagrams of 2-bridge knots with certain specific twists.


title: Learned Construction Grammars Converge Across Registers Given Increased Exposure

abstract: This paper measures the impact of increased exposure on whether learned construction grammars converge onto shared representations when trained on data from different registers. Register influences the frequency of constructions, with some structures common in formal but not informal usage. We expect that a grammar induction algorithm exposed to different registers will acquire different constructions. To what degree does increased exposure lead to the convergence of register-specific grammars? The experiments in this paper simulate language learning in 12 languages (half Germanic and half Romance) with corpora representing three registers (Twitter, Wikipedia, Web). These simulations are repeated with increasing amounts of exposure, from 100k to 2 million words, to measure the impact of exposure on the convergence of grammars. The results show that increased exposure does lead to converging grammars across all languages. In addition, a shared core of register-universal constructions remains constant across increasing amounts of exposure.


title: Accurate and Generalizable Quantitative Scoring of Liver Steatosis from Ultrasound Images via Scalable Deep Learning

abstract: Background & Aims: Hepatic steatosis is a major cause of chronic liver disease. 2D ultrasound is the most widely used non-invasive tool for screening and monitoring, but associated diagnoses are highly subjective. We developed a scalable deep learning (DL) algorithm for quantitative scoring of liver steatosis from 2D ultrasound images. Approach & Results: Using retrospectively collected multi-view ultrasound data from 3,310 patients, 19,513 studies, and 228,075 images, we trained a DL algorithm to diagnose steatosis stages (healthy, mild, moderate, or severe) from ultrasound diagnoses. Performance was validated on two multi-scanner unblinded and blinded (initially to DL developer) histology-proven cohorts (147 and 112 patients) with histopathology fatty cell percentage diagnoses, and a subset with FibroScan diagnoses. We also quantified reliability across scanners and viewpoints. Results were evaluated using Bland-Altman and receiver operating characteristic (ROC) analysis. The DL algorithm demonstrates repeatable measurements with a moderate number of images (3 for each viewpoint) and high agreement across 3 premium ultrasound scanners. High diagnostic performance was observed across all viewpoints: area under the curves of the ROC to classify >=mild, >=moderate, =severe steatosis grades were 0.85, 0.90, and 0.93, respectively. The DL algorithm outperformed or performed at least comparably to FibroScan with statistically significant improvements for all levels on the unblinded histology-proven cohort, and for =severe steatosis on the blinded histology-proven cohort. Conclusions: The DL algorithm provides a reliable quantitative steatosis assessment across view and scanners on two multi-scanner cohorts. Diagnostic performance was high with comparable or better performance than FibroScan.


title: Are you doing what I say? On modalities alignment in ALFRED

abstract: ALFRED is a recently proposed benchmark that requires a model to complete tasks in simulated house environments specified by instructions in natural language. We hypothesize that key to success is accurately aligning the text modality with visual inputs. Motivated by this, we inspect how well existing models can align these modalities using our proposed intrinsic metric, boundary adherence score (BAS). The results show the previous models are indeed failing to perform proper alignment. To address this issue, we introduce approaches aimed at improving model alignment and demonstrate how improved alignment, improves end task performance.


title: Structural Study of Adsorbed Helium Films: New Approach with Synchrotron Radiation X-rays

abstract: A few atomic layers of helium adsorbed on graphite have been attracting much attention as one of the ideal quantum systems in two dimension. Although previous reports on neutron diffraction have shown fundamental structural information in these systems, there still remain many open questions. Here, we propose surface crystal truncation rod (CTR) scatterings using synchrotron radiation X-rays as a promising method to reveal surface and interface structures of helium films on graphite at temperatures below 2 K, based on the preliminary experimental results on a monolayer of He-4 on a thin graphite. Our estimation on heat generation by X-ray irradiations also suggests that CTR scatterings are applicable to even at system temperatures near 100 mK.


title: Why Lottery Ticket Wins? A Theoretical Perspective of Sample Complexity on Pruned Neural Networks

abstract: The \textit{lottery ticket hypothesis} (LTH) states that learning on a properly pruned network (the \textit{winning ticket}) improves test accuracy over the original unpruned network. Although LTH has been justified empirically in a broad range of deep neural network (DNN) involved applications like computer vision and natural language processing, the theoretical validation of the improved generalization of a winning ticket remains elusive. To the best of our knowledge, our work, for the first time, characterizes the performance of training a pruned neural network by analyzing the geometric structure of the objective function and the sample complexity to achieve zero generalization error. We show that the convex region near a desirable model with guaranteed generalization enlarges as the neural network model is pruned, indicating the structural importance of a winning ticket. Moreover, when the algorithm for training a pruned neural network is specified as an (accelerated) stochastic gradient descent algorithm, we theoretically show that the number of samples required for achieving zero generalization error is proportional to the number of the non-pruned weights in the hidden layer. With a fixed number of samples, training a pruned neural network enjoys a faster convergence rate to the desired model than training the original unpruned one, providing a formal justification of the improved generalization of the winning ticket. Our theoretical results are acquired from learning a pruned neural network of one hidden layer, while experimental results are further provided to justify the implications in pruning multi-layer neural networks.


title: NAS-Bench-360: Benchmarking Neural Architecture Search on Diverse Tasks

abstract: Most existing neural architecture search (NAS) benchmarks and algorithms prioritize well-studied tasks, e.g. image classification on CIFAR or ImageNet. This makes the performance of NAS approaches in more diverse areas poorly understood. In this paper, we present NAS-Bench-360, a benchmark suite to evaluate methods on domains beyond those traditionally studied in architecture search, and use it to address the following question: do state-of-the-art NAS methods perform well on diverse tasks? To construct the benchmark, we curate ten tasks spanning a diverse array of application domains, dataset sizes, problem dimensionalities, and learning objectives. Each task is carefully chosen to interoperate with modern CNN-based search methods while possibly being far-afield from its original development domain. To speed up and reduce the cost of NAS research, for two of the tasks we release the precomputed performance of 15,625 architectures comprising a standard CNN search space. Experimentally, we show the need for more robust NAS evaluation of the kind NAS-Bench-360 enables by showing that several modern NAS procedures perform inconsistently across the ten tasks, with many catastrophically poor results. We also demonstrate how NAS-Bench-360 and its associated precomputed results will enable future scientific discoveries by testing whether several recent hypotheses promoted in the NAS literature hold on diverse tasks. NAS-Bench-360 is hosted at https://nb360.ml.cmu.edu.


title: Adaptive Feedforward Reference Design for Active Vibration Rejection in Multi-Actuator Hard Disk Drives

abstract: In December 2017, Seagate unveiled the Multi Actuator Technology to double the data performance of the future generation hard disk drives (HDD). This technology will equip drives with two dual stage actuators (DSA) each comprising of a voice coil motor (VCM) actuator and a piezoelectric micro actuator (MA) operating on the same pivot point. Each DSA is responsible for controlling half of the drive's arms. As both the DSAs operate independently on the same pivot timber, the control forces and torques generated by one can affect the operation of the other and thereby worsening the performance drastically. In this paper, a robust adaptive feedforward controller is designed as an add-on controller to an existing stabilizing feedback controller to reject the disturbances transferred through the common pivot timber by shaping the references to the VCM actuator and the total output of the dual stage system.


title: Eigenvalues and cycles of consecutive lengths

abstract: As the counterpart of classical theorems on cycles of consecutive lengths due to Bondy and Bollob\'as in spectral graph theory, Nikiforov proposed the following open problem in 2008: What is the maximum $C$ such that for all positive $\varepsilon<C$ and sufficiently large $n$, every graph $G$ of order $n$ with spectral radius $\rho(G)>\sqrt{\lfloor\frac{n^2}{4}\rfloor}$ contains a cycle of length $\ell$ for each integer $\ell\in[3,(C-\varepsilon)n]$. We prove that $C\geq\frac{1}{4}$ by a novel method, improving the existing bounds. Besides several novel ideas, our proof technique is partly inspirited by the recent research on Ramsey numbers of star versus large even cycles due to Allen, {\L}uczak, Polcyn and Zhang, and with aid of a powerful spectral inequality. We also derive an Erd\H{o}s-Gallai-type edge number condition for even cycles, which may be of independent interest.


title: Predicting the Stereoselectivity of Chemical Transformations by Machine Learning

abstract: Stereoselective reactions (both chemical and enzymatic reactions) have been essential for origin of life, evolution, human biology and medicine. Since late 1960s, there have been numerous successes in the exciting new frontier of asymmetric catalysis. However, most industrial and academic asymmetric catalysis nowadays do follow the trial-and-error model, since the energetic difference for success or failure in asymmetric catalysis is incredibly small. Our current understanding about stereoselective reactions is mostly qualitative that stereoselectivity arises from differences in steric effects and electronic effects in multiple competing mechanistic pathways. Quantitatively understanding and modulating the stereoselectivity of for a given chemical reaction still remains extremely difficult. As a proof of principle, we herein present a novel machine learning technique, which combines a LASSO model and two Random Forest model via two Gaussian Mixture models, for quantitatively predicting stereoselectivity of chemical reactions. Compared to the recent ground-breaking approach [1], our approach is able to capture interactions between features and exploit complex data distributions, which are important for predicting stereoselectivity. Experimental results on a recently published dataset demonstrate that our approach significantly outperform [1]. The insight obtained from our results provide a solid foundation for further exploration of other synthetically valuable yet mechanistically intriguing stereoselective reactions.


title: Data-Driven Strictly Positive Real System Identification with prior System Knowledge

abstract: Strictly Positive Real (SPR) transfer functions arise in many areas of engineering like passivity theory in circuit analysis and adaptive control to name a few. In many physical systems, it is possible to conclude that the system is Positive Real (PR) or SPR but system identification algorithms might produce estimates which are not SPR. In this paper, an algorithm to approximate frequency response data with SPR transfer functions using Generalized Orthonormal Basis Functions (GOBFs) is presented. Prior knowledge of the system helps us to get approximate pole locations, which can then be used to construct GOBFs. Next, a convex optimization problem will be formulated to obtain an estimate of the SPR transfer function.


title: Fundamental Distinction between Intrinsic and Extrinsic Nonlinear Thermal Hall Effects

abstract: We theoretically investigated the fundamental distinction between intrinsic and extrinsic nonlinear thermal Hall effect in the presence of disorder at the second-order response to the temperature gradient in terms of the semi-classical Boltzmann equation. We found that, at low temperatures, the intrinsic contribution of the nonlinear thermal Hall conductivity is proportional to the square of temperature, whereas the extrinsic contributions (side-jump and skew-scattering) are independent of temperature. This distinct dependency on temperature provide a new approach to readily distinguish the intrinsic and extrinsic contributions. Specifically, we analysed the nonlinear thermal Hall effect for a tilted two-dimensional massive Dirac material. In particular, we showed that when the Fermi energy is located at the Dirac point, the signal is solely from the intrinsic mechanism; when the Fermi energy is higher, the extrinsic contributions are dominant, which are uncovered to be two to three orders of magnitude larger than the intrinsic contribution.