title: Diverse Linguistic Features for Assessing Reading Difficulty of Educational Filipino Texts

abstract: In order to ensure quality and effective learning, fluency, and comprehension, the proper identification of the difficulty levels of reading materials should be observed. In this paper, we describe the development of automatic machine learning-based readability assessment models for educational Filipino texts using the most diverse set of linguistic features for the language. Results show that using a Random Forest model obtained a high performance of 62.7% in terms of accuracy, and 66.1% when using the optimal combination of feature sets consisting of traditional and syllable pattern-based predictors.


title: The Inelastic Market Hypothesis: A Microstructural Interpretation

abstract: We attempt to reconcile Gabaix and Koijen's (GK) recent Inelastic Market Hypothesis (IMH) with the order-driven view of markets that emerged within the microstructure literature in the past 20 years. We review the most salient empirical facts and arguments that give credence to the idea that market price fluctuations are mostly due to order flow, whether informed or non-informed. We show that the Latent Liquidity Theory of price impact makes a precise prediction for GK's multiplier $M$, which measures by how many dollars, on average, the market value of a company goes up if one buys one dollar worth of its stocks. Our central result is that $M$ is of order unity, as found by GK, and increases with the volatility of the stock and decreases with the fraction of the market cap. traded daily. We discuss several empirical results suggesting that the lion's share of volatility is due to trading activity. We argue that the IMH holds for all asset classes, beyond the case of stock markets considered by GK.


title: Assignment of a Synthetic Population for Activity-based Modelling employing Publicly Available Data

abstract: Agent based modelling has acquired the spotlight in the transportation domain both in scientific literature and in real life applications, thanks to its capability to deal with the ever-growing complexity of transportation systems, including future disrupting mobility technologies and services such as automated driving, Mobility as a Service, and micromobility. Different software emerged, dedicated to the simulation of disaggregate travel demand framing individual choices based on the profile of each agent. Still, the actual research work exploiting these models is scarce and the professionals with the knowledge to use them are few. This may be ascribed to the large amount of needed input data or to a lack of commercial solutions and of research production detailing the process leading to the actual simulations. In this paper, a methodology to spatially assign a synthetic population by exploiting publicly available aggregate data is presented and implemented on a case study. In doing so, the paper provides a systematic approach for a quick and efficient treatment of the data needed for activity-based demand generation. Finally, the obtained dataset, representing a synthetic population of the city of Tallinn, Estonia, and its spatial assignment, is described so that it may be exploited by fellow researchers, since both the tools needed for spatial assignment and the resulting dataset are made available as open source.


title: Mean field game equations with underlying jump-diffusion process

abstract: We consider a couple of integrodifferential PDEs arising from a stochastic Markovian control problem subjected to initial-terminal conditions. These equations correspond to the MFG system for a controlled jump-diffusion process. We prove that for a specific choice of the control function the expectation of the jump-diffusion process can be found explicitly. The study is an extension of similar results known for the pure diffusion process. As an example, we show how this can be applied to the problem of investors evaluating the trend of an asset when choosing an optimal portfolio.


title: Bipartite Graft III: General Case

abstract: This paper is a sequel of our previous paper (N. Kita: Bipartite graft {II}: Cathedral decomposition for combs. arXiv preprint arXiv:2101.06678, 2021). In our previous paper, a graft analogue of the Dulmage-Mendelsohn decomposition has been introduced for comb bipartite grafts. In this paper, we prove how this result can be extended for general bipartite grafts.


title: HR-Crime: Human-Related Anomaly Detection in Surveillance Videos

abstract: The automatic detection of anomalies captured by surveillance settings is essential for speeding the otherwise laborious approach. To date, UCF-Crime is the largest available dataset for automatic visual analysis of anomalies and consists of real-world crime scenes of various categories. In this paper, we introduce HR-Crime, a subset of the UCF-Crime dataset suitable for human-related anomaly detection tasks. We rely on state-of-the-art techniques to build the feature extraction pipeline for human-related anomaly detection. Furthermore, we present the baseline anomaly detection analysis on the HR-Crime. HR-Crime as well as the developed feature extraction pipeline and the extracted features will be publicly available for further research in the field.


title: Fourier orthogonal series on a paraboloid

abstract: We study orthogonal structures and Fourier orthogonal series on the surface of a paraboloid $\mathbb{V}_0^{d+1} = \{(x,t): \|x\| = \sqrt{t}, \, x \in \mathbb{R}^d, \, 0\le t<1\}$. The reproducing kernels of the orthogonal polynomials with respect to $t^\beta(1-t)^\gamma$ on $\mathbb{V}_0^{d+1}$ are related to the reproducing kernels of the Jacobi polynomials on the parabolic domain $\{(x_1,x_2): x_1^2 \le x_2 \le 1\}$ in $\mathbb{R}^2$. This connection serves as an essential tool for our study of the Fourier orthogonal series on the surface of the paraboloid, which allow us, in particular, to study the convergence of the Ces\`aro means on the surface. Analogous results are also established for the solid paraboloid bounded by $\mathbb{V}_0^{d+1}$ and the hyperplane $t=1$.


title: Threshold scattering for the focusing NLS with a repulsive Dirac delta potential

abstract: We establish the scattering of solutions to the focusing mass supercritical nonlinear Schr\"odinger equation with a repulsive Dirac delta potential \[ i\partial_{t}u+\partial^{2}_{x}u+\gamma\delta(x)u+|u|^{p-1}u=0, \quad (t,x)\in {\mathbb R}\times{\mathbb R}, \] at the mass-energy threshold, namely, when $E_{\gamma}(u_{0})[M(u_{0})]^{\sigma}=E_{0}(Q)[M(Q)]^{\sigma}$ where $u_{0}\in H^{1}({\mathbb R})$ is the initial data, $Q$ is the ground state of the free NLS on the real line ${\mathbb R}$, $E_{\gamma}$ is the energy, $M$ is the mass and $\sigma=(p+3)/(p-5)$. We also prove failure of the uniform space-time bounds at the mass-energy threshold.


title: SyDog: A Synthetic Dog Dataset for Improved 2D Pose Estimation

abstract: Estimating the pose of animals can facilitate the understanding of animal motion which is fundamental in disciplines such as biomechanics, neuroscience, ethology, robotics and the entertainment industry. Human pose estimation models have achieved high performance due to the huge amount of training data available. Achieving the same results for animal pose estimation is challenging due to the lack of animal pose datasets. To address this problem we introduce SyDog: a synthetic dataset of dogs containing ground truth pose and bounding box coordinates which was generated using the game engine, Unity. We demonstrate that pose estimation models trained on SyDog achieve better performance than models trained purely on real data and significantly reduce the need for the labour intensive labelling of images. We release the SyDog dataset as a training and evaluation benchmark for research in animal motion.


title: Bayesian analysis of the prevalence bias: learning and predicting from imbalanced data

abstract: Datasets are rarely a realistic approximation of the target population. Say, prevalence is misrepresented, image quality is above clinical standards, etc. This mismatch is known as sampling bias. Sampling biases are a major hindrance for machine learning models. They cause significant gaps between model performance in the lab and in the real world. Our work is a solution to prevalence bias. Prevalence bias is the discrepancy between the prevalence of a pathology and its sampling rate in the training dataset, introduced upon collecting data or due to the practioner rebalancing the training batches. This paper lays the theoretical and computational framework for training models, and for prediction, in the presence of prevalence bias. Concretely a bias-corrected loss function, as well as bias-corrected predictive rules, are derived under the principles of Bayesian risk minimization. The loss exhibits a direct connection to the information gain. It offers a principled alternative to heuristic training losses and complements test-time procedures based on selecting an operating point from summary curves. It integrates seamlessly in the current paradigm of (deep) learning using stochastic backpropagation and naturally with Bayesian models.


title: Modulating the temporal dynamics of nonlinear ultrafast plasmon resonances

abstract: Spatio-temporal control of ultrafast plasmon resonances has gained research interest in recent years because of their tremendous implications in nonlinear optics and ultrafast quantum technology. In particular, the lifetime of ultrashort plasmon oscillations has become a debatable subject in recent experimental and theoretical studies to fulfill the future challenges concerning their effective employment in the vast applications of the plasmonic industry. Here, we examined the temporal properties of nonlinear plasmonic modes in metal nanostructures by interacting them with quantum objects in the weak coupling regime in order to distinguish it from the fundamental plasmonic mode. First of all, we present an analytical description of nonlinear ultrafast dynamics of localized surface plasmon resonances when the second harmonic plasmon mode interacts with long-lived dark mode or quantum emitter. Later, the coupled plasmonic system is realized in two different ways to control the lifetime of second harmonic mode by coupling, i) driven mode to dark mode (or long lifetime quantum emitter) ii)itself to dark mode (or long lifetime quantum emitter). The driven-dissipative dynamics are solved through a numerical technique governing the spatial and temporal changes in the second harmonic plasmonic response supported by AuNP. Finally, the lifetime enhancement of nonlinear plasmon mode is manifested by performing FDTD simulations for a nonlinear plasmonic system of Au nanoparticles coupled with a long lifetime quantum emitter.


title: Ph\"anomen-Signal-Modell: Formalismus, Graph und Anwendung

abstract: If we consider information as the basis of action, it may be of interest to examine the flow and acquisition of information between the actors in traffic. The central question is, which signals an automaton has to receive, decode or send in road traffic in order to act safely and in a conform manner to valid standards. The phenomenon-signal-model is a method to structure the problem, to analyze and to describe this very signal flow. Explaining the basics, structure and application of this method is the aim of this paper. -- Betrachtet man Information als Grundlage des Handelns, so wird es interessant sein, Fluss und Erfassung von Information zwischen den Akteuren des Verkehrsgeschehens zu untersuchen. Die zentrale Frage ist, welche Signale ein Automat im Stra{\ss}enverkehr empfangen, decodieren oder senden muss, um konform zu geltenden Ma{\ss}st\"aben und sicher zu agieren. Das Ph\"anomen-Signal-Modell ist eine Methode, das Problemfeld zu strukturieren, eben diesen Signalfluss zu analysieren und zu beschreiben. Der vorliegende Aufsatz erkl\"art Grundlagen, Aufbau und Anwendung dieser Methode.


title: Sphaleron bound in some nonstandard cosmology scenarios

abstract: In the first scenario, we revise the upper bound of the cutoff scale of the dimension-six potential Higgs operator required for a successful electroweak baryogenesis in the case of a modified expansion history caused by the existence of a non-interacting scalar field at the time phase transition happens. The upper bound 860 GeV of the cutoff scale in the conventional case can be improved to 1 TeV in the modified expansion case under certain conditions. In the second scenario, we consider the Randall-Sundrum type II model. We show that the lower bound of the five-dimensional Planck scale in this model, which is determined from the validity of Newtonian gravitational potential at small distance, turns out to be crucial in eliminating this model as a viable candidate to satisfy the sphaleron bound; however, again modifying the expansion history by including a non-interacting scalar field at the electroweak scale can then make this model satisfy the sphaleron bound with a certain parameter space.


title: Oscillating scalar dissipating in a medium

abstract: We study how oscillations of a scalar field condensate are damped due to dissipative effects in a thermal medium. Our starting point is a non-linear and non-local condensate equation of motion descending from a 2PI-resummed effective action derived in the Schwinger-Keldysh formalism appropriate for non-equilibrium quantum field theory. We solve this non-local equation by means of multiple-scale perturbation theory appropriate for time-dependent systems, obtaining approximate analytic solutions valid for very long times. The non-linear effects lead to power-law damping of oscillations, that at late times transition to exponentially damped ones characteristic for linear systems. These solutions describe the evolution very well, as we demonstrate numerically in a number of examples. We then approximate the non-local equation of motion by a Markovianised one, resolving the ambiguities appearing in the process, and solve it utilizing the same methods to find the very same leading approximate solution. This comparison justifies the use of Markovian equations at leading order. The standard time-dependent perturbation theory in comparison is not capable of describing the non-linear condensate evolution beyond the early time regime of negligible damping. The macroscopic evolution of the condensate is interpreted in terms of microphysical particle processes. Our results have implications for the quantitative description of the decay of cosmological scalar fields in the early Universe, and may also be applied to other physical systems.


title: Possible Realization of Optical Quadratic and Dirac Points in Woodpile Photonic Crystals

abstract: The simulation of fermionic relativistic physics, e.g., Dirac and Weyl physics, has led to the discovery of many unprecedented phenomena in photonics, of which the optical-frequency realization is, however, still challenging. Here, surprisingly, we discover that the woodpile photonic crystals commonly used for optical frequency applications host exotic fermion-like relativistic degeneracies: a Dirac nodal line and a fourfold quadratic point, as protected by the nonsymmorphic crystalline symmetry. Deforming the woodpile photonic crystal leads to the emergence of type-II Dirac points from the fourfold quadratic point. Such type-II Dirac points can be detected by its anomalous refraction property which is manifested as a giant birefringence in a slab setup. Our findings provide a promising route towards 3D optical Dirac physics in all-dielectric photonic crystals.


title: An Intelligent Energy Management Framework for Hybrid-Electric Propulsion Systems Using Deep Reinforcement Learning

abstract: Hybrid-electric propulsion systems powered by clean energy derived from renewable sources offer a promising approach to decarbonise the world's transportation systems. Effective energy management systems are critical for such systems to achieve optimised operational performance. However, developing an intelligent energy management system for applications such as ships operating in a highly stochastic environment and requiring concurrent control over multiple power sources presents challenges. This article proposes an intelligent energy management framework for hybrid-electric propulsion systems using deep reinforcement learning. In the proposed framework, a Twin-Delayed Deep Deterministic Policy Gradient agent is trained using an extensive volume of historical load profiles to generate a generic energy management strategy. The strategy, i.e. the core of the energy management system, can concurrently control multiple power sources in continuous state and action spaces. The proposed framework is applied to a coastal ferry model with multiple fuel cell clusters and a battery, achieving near-optimal cost performance when applied to novel future voyages.


title: BoA-PTA, A Bayesian Optimization Accelerated Error-Free SPICE Solver

abstract: One of the greatest challenges in IC design is the repeated executions of computationally expensive SPICE simulations, particularly when highly complex chip testing/verification is involved. Recently, pseudo transient analysis (PTA) has shown to be one of the most promising continuation SPICE solver. However, the PTA efficiency is highly influenced by the inserted pseudo-parameters. In this work, we proposed BoA-PTA, a Bayesian optimization accelerated PTA that can substantially accelerate simulations and improve convergence performance without introducing extra errors. Furthermore, our method does not require any pre-computation data or offline training. The acceleration framework can either be implemented to speed up ongoing repeated simulations immediately or to improve new simulations of completely different circuits. BoA-PTA is equipped with cutting-edge machine learning techniques, e.g., deep learning, Gaussian process, Bayesian optimization, non-stationary monotonic transformation, and variational inference via parameterization. We assess BoA-PTA in 43 benchmark circuits against other SOTA SPICE solvers and demonstrate an average 2.3x (maximum 3.5x) speed-up over the original CEPTA.


title: Analytic characterization of monotone Hopf-harmonics

abstract: We study solutions of the inner-variational equation associated with the Dirichlet energy in the plane, given homeomorphic Sobolev boundary data. We prove that such a solution is monotone if and only if its Jacobian determinant does not change sign. These solutions, called monotone Hopf-harmonics, are a natural alternative to harmonic homeomorphisms. Examining the topological behavior of a solution (not a priori monotone) on the trajectories of Hopf quadratic differentials plays a sizable role in our arguments.


title: How much pre-training is enough to discover a good subnetwork?

abstract: Neural network pruning is useful for discovering efficient, high-performing subnetworks within pre-trained, dense network architectures. More often than not, it involves a three-step process -- pre-training, pruning, and re-training -- that is computationally expensive, as the dense model must be fully pre-trained. While previous work has revealed through experiments the relationship between the amount of pre-training and the performance of the pruned network, a theoretical characterization of such dependency is still missing. Aiming to mathematically analyze the amount of dense network pre-training needed for a pruned network to perform well, we discover a simple theoretical bound in the number of gradient descent pre-training iterations on a two-layer, fully-connected network, beyond which pruning via greedy forward selection [61] yields a subnetwork that achieves good training error. Interestingly, this threshold is shown to be logarithmically dependent upon the size of the dataset, meaning that experiments with larger datasets require more pre-training for subnetworks obtained via pruning to perform well. Lastly, we empirically validate our theoretical results on a multi-layer perceptron trained on MNIST.


title: Pseudo-symmetric pairs for Kac-Moody algebras

abstract: Lie algebra involutions and their fixed-point subalgebras give rise to symmetric spaces and real forms of complex Lie algebras, and are well-studied in the context of symmetrizable Kac-Moody algebras. In this paper we study a generalization. Namely, we introduce the concept of a pseudo-involution, an automorphism which is only required to act involutively on a stable Cartan subalgebra, and the concept of a pseudo-fixed-point subalgebra, a natural substitute for the fixed-point subalgebra. In the symmetrizable Kac-Moody setting, we give a comprehensive discussion of pseudo-involutions of the second kind, the associated pseudo-fixed-point subalgebras, restricted root systems and Weyl groups, in terms of generalizations of Satake diagrams.


title: ECLARE: Extreme Classification with Label Graph Correlations

abstract: Deep extreme classification (XC) seeks to train deep architectures that can tag a data point with its most relevant subset of labels from an extremely large label set. The core utility of XC comes from predicting labels that are rarely seen during training. Such rare labels hold the key to personalized recommendations that can delight and surprise a user. However, the large number of rare labels and small amount of training data per rare label offer significant statistical and computational challenges. State-of-the-art deep XC methods attempt to remedy this by incorporating textual descriptions of labels but do not adequately address the problem. This paper presents ECLARE, a scalable deep learning architecture that incorporates not only label text, but also label correlations, to offer accurate real-time predictions within a few milliseconds. Core contributions of ECLARE include a frugal architecture and scalable techniques to train deep models along with label correlation graphs at the scale of millions of labels. In particular, ECLARE offers predictions that are 2 to 14% more accurate on both publicly available benchmark datasets as well as proprietary datasets for a related products recommendation task sourced from the Bing search engine. Code for ECLARE is available at https://github.com/Extreme-classification/ECLARE.


title: Speech2AffectiveGestures: Synthesizing Co-Speech Gestures with Generative Adversarial Affective Expression Learning

abstract: We present a generative adversarial network to synthesize 3D pose sequences of co-speech upper-body gestures with appropriate affective expressions. Our network consists of two components: a generator to synthesize gestures from a joint embedding space of features encoded from the input speech and the seed poses, and a discriminator to distinguish between the synthesized pose sequences and real 3D pose sequences. We leverage the Mel-frequency cepstral coefficients and the text transcript computed from the input speech in separate encoders in our generator to learn the desired sentiments and the associated affective cues. We design an affective encoder using multi-scale spatial-temporal graph convolutions to transform 3D pose sequences into latent, pose-based affective features. We use our affective encoder in both our generator, where it learns affective features from the seed poses to guide the gesture synthesis, and our discriminator, where it enforces the synthesized gestures to contain the appropriate affective expressions. We perform extensive evaluations on two benchmark datasets for gesture synthesis from the speech, the TED Gesture Dataset and the GENEA Challenge 2020 Dataset. Compared to the best baselines, we improve the mean absolute joint error by 10--33%, the mean acceleration difference by 8--58%, and the Fr\'echet Gesture Distance by 21--34%. We also conduct a user study and observe that compared to the best current baselines, around 15.28% of participants indicated our synthesized gestures appear more plausible, and around 16.32% of participants felt the gestures had more appropriate affective expressions aligned with the speech.


title: Functional estimation in log-concave location families

abstract: Let $\{P_{\theta}:\theta \in {\mathbb R}^d\}$ be a log-concave location family with $P_{\theta}(dx)=e^{-V(x-\theta)}dx,$ where $V:{\mathbb R}^d\mapsto {\mathbb R}$ is a known convex function and let $X_1,\dots, X_n$ be i.i.d. r.v. sampled from distribution $P_{\theta}$ with an unknown location parameter $\theta.$ The goal is to estimate the value $f(\theta)$ of a smooth functional $f:{\mathbb R}^d\mapsto {\mathbb R}$ based on observations $X_1,\dots, X_n.$ In the case when $V$ is sufficiently smooth and $f$ is a functional from a ball in a H\"older space $C^s,$ we develop estimators of $f(\theta)$ with minimax optimal error rates measured by the $L_2({\mathbb P}_{\theta})$-distance as well as by more general Orlicz norm distances. Moreover, we show that if $d\leq n^{\alpha}$ and $s>\frac{1}{1-\alpha},$ then the resulting estimators are asymptotically efficient in H\'ajek-LeCam sense with the convergence rate $\sqrt{n}.$ This generalizes earlier results on estimation of smooth functionals in Gaussian shift models. The estimators have the form $f_k(\hat \theta),$ where $\hat \theta$ is the maximum likelihood estimator and $f_k: {\mathbb R}^d\mapsto {\mathbb R}$ (with $k$ depending on $s$) are functionals defined in terms of $f$ and designed to provide a higher order bias reduction in functional estimation problem. The method of bias reduction is based on iterative parametric bootstrap and it has been successfully used before in the case of Gaussian models.


title: Time evolution of a mean-field generalized contact process

abstract: We investigate the macroscopic time evolution and stationary states of a mean field generalized contact process in $\mathbb{R}^d$. The model is described by a coupled set of nonlinear integral-differential equations. It was inspired by a model of neurons with discrete voltages evolving by a stochastic integrate and fire mechanism. We obtain a complete solution in the spatially uniform case and partial solutions in the general case. The system has one or more fixed points and also traveling wave solutions.

categories: math.PR cond-mat.dis-nn cond-mat.stat-mech nlin.AO q-bio.PE


title: Localization-enhanced dissipation in a generalized Aubry-Andr\'{e}-Harper model coupled with Ohmic baths

abstract: In this work, the exact dynamics of excitation in the generalized Aubry-Andr\'{e}-Harper model coupled with an Ohmic-type environment is discussed by evaluating the survival probability and inverse participation ratio of the state of system. In contrast to the common belief that localization will preserve the information of the initial state in the system against dissipation into the environment, our study found that strong localization can enhance the dissipation of quantum information instead. By a thorough examination of the dynamics, we show that the coherent transition between the energy state of system is crucial for understanding this unusual behavior. Under this circumstance, the coupling induced energy exchange between the system and its environment can induce the periodic population of excitation on the states of system. As a result, the stable or localization-enhanced decaying of excitation can be observed, dependent on the energy difference between the states of system. This point is verified in further by checking the varying of dynamics of excitation in the system when the coupling between the system and environment is more strong.


title: MOG cosmology without dark matter and the cosmological constant

abstract: In this work, we investigate the MOdified Gravity (MOG) theory for dynamics of the universe and compare the results with the $\Lambda$CDM cosmology. We study the background cosmological properties of the MOG model and structure formation at the linear perturbation level. We compare the two models with the currently available cosmological data by using statistical Bayesian analyses. After obtaining updated constraints on the free parameters, we use some methods of model selection to assist in choosing the more consistent model such as the reduced chi-squared ($\chi^2_{\rm red}$) and a number of the basic information criteria such as the Akaike Information Criterion (AIC), the Bayes factor or Bayesian Information Criterion (BIC), and Deviance Information Criterion (DIC). MOG model appears to be consistent with the $\Lambda$CDM model by the results of $\chi^2_{\rm red}$ and DIC for an overall statistical analysis using the background data and the linear growth of structure formation.


title: Supercritical accretion of stellar-mass compact objects in active galactic nuclei

abstract: Accretion disks of active galactic nuclei (AGN) have been proposed as promising sites for producing both (stellar-mass) compact object mergers and extreme mass ratio inspirals. Along with the disk-assisted migration/evolution process, ambient gas materials inevitably accrete onto the compact objects. The description of this process is subject to significant theoretical uncertainties in previous studies. It was commonly assumed that either an Eddington accretion rate or a Bondi accretion rate (or any rate in between) takes place, although these two rates can differ from each other by several orders of magnitude. As a result, the mass and spin evolution of compact objects within AGN disks are essentially unknown. In this work, we construct a relativistic supercritical inflow-outflow model for black hole (BH) accretion. We show that the radiation efficiency of the supercritical accretion of a stellar-mass BH (sBH) is generally too low to explain the proposed electromagnetic counterpart of GW190521. Applying this model to sBHs embedded in AGN disks, we find that, although the gas inflow rates at Bondi radii of these sBHs are in general highly super-Eddington, a large fraction of inflowing gas eventually escapes as outflows so that only a small fraction accretes onto the sBH, resulting in mildly super-Eddington BH absorption in most cases. We also implement this inflow-outflow model to study the evolution of neutron stars (NS) and white dwarfs (WD) in AGN disks, taking into account corrections from star sizes and star magnetic fields. It turns out to be difficult for WDs to grow to the Chandrasekhar limit via accretion because WDs are spun up more efficiently to reach the shedding limit before the Chandrasekhar limit. For NSs the accretion-induced collapse is possible if NS magnetic fields are sufficiently strong, keeping the NS in a slow rotation state during accretion.


title: RLTutor: Reinforcement Learning Based Adaptive Tutoring System by Modeling Virtual Student with Fewer Interactions

abstract: A major challenge in the field of education is providing review schedules that present learned items at appropriate intervals to each student so that memory is retained over time. In recent years, attempts have been made to formulate item reviews as sequential decision-making problems to realize adaptive instruction based on the knowledge state of students. It has been reported previously that reinforcement learning can help realize mathematical models of students learning strategies to maintain a high memory rate. However, optimization using reinforcement learning requires a large number of interactions, and thus it cannot be applied directly to actual students. In this study, we propose a framework for optimizing teaching strategies by constructing a virtual model of the student while minimizing the interaction with the actual teaching target. In addition, we conducted an experiment considering actual instructions using the mathematical model and confirmed that the model performance is comparable to that of conventional teaching methods. Our framework can directly substitute mathematical models used in experiments with human students, and our results can serve as a buffer between theoretical instructional optimization and practical applications in e-learning systems.


title: Quantum Representation Theory and Manin matrices I: finite-dimensional case

abstract: We construct Quantum Representation Theory which describes quantum analogue of representations in frame of "non-commutative linear geometry" developed by Manin. To do it we generalise the internal hom-functor to the case of adjunction with a parameter and construct a general approach to representations of a monoid in a symmetric monoidal category with a parameter subcategory. Quantum Representation Theory is obtained by application of this approach to a monoidal category of some class of graded algebras with Manin product, where the parameter subcategory consists of connected finitely generated quadratic algebras. We formulate this theory in the language of Manin matrices and obtain quantum analogues of direct sum and tensor product of representations. Finally, we give some examples of quantum representations.


title: Opinion Prediction with User Fingerprinting

abstract: Opinion prediction is an emerging research area with diverse real-world applications, such as market research and situational awareness. We identify two lines of approaches to the problem of opinion prediction. One uses topic-based sentiment analysis with time-series modeling, while the other uses static embedding of text. The latter approaches seek user-specific solutions by generating user fingerprints. Such approaches are useful in predicting user's reactions to unseen content. In this work, we propose a novel dynamic fingerprinting method that leverages contextual embedding of user's comments conditioned on relevant user's reading history. We integrate BERT variants with a recurrent neural network to generate predictions. The results show up to 13\% improvement in micro F1-score compared to previous approaches. Experimental results show novel insights that were previously unknown such as better predictions for an increase in dynamic history length, the impact of the nature of the article on performance, thereby laying the foundation for further research.


title: Writhing and hockling instabilities in twisted elastic fibers

abstract: The buckling and twisting of slender, elastic fibers is a deep and well-studied field. A slender elastic rod that is twisted with respect to a fixed end will spontaneously form a loop, or hockle, to relieve the torsional stress that builds. Further twisting results in the formation of plectonemes -- a helical excursion in the fiber that extends with additional twisting. Here we use an idealized, micron-scale experiment to investigate the energy stored, and subsequently released, by hockles and plectonemes as they are pulled apart, in analogy with force spectroscopy studies of DNA and protein folding. Hysteresis loops in the snapping and unsnapping inform the stored energy in the twisted fiber structures.


title: Multivariate $\alpha$-normal distributions

abstract: The Weibull distribution can be obtained using a power transformation from the standard exponential distribution. In this article, we will consider a symmetrized power transformation of a random variable with the standard normal distribution. We will call its distribution the $\alpha$-{\it normal (Gaussian) distribution}. We examine properties of this distribution in detail. We calculate moments and consider the moment problem of $\alpha$-normal distribution. We derive the formula of its differential entropy and (exponential) Orlicz norm. % of $\alpha$-normal random variables. Moreover, we define the joint distribution function of the multivariate $\alpha$-normal distribution as a meta-Gaussian distribution with $\alpha$-normal marginals. We consider also the limiting distribution as $\alpha$ tends to infinity.


title: Towards explainable artificial intelligence (XAI) for early anticipation of traffic accidents

abstract: Traffic accident anticipation is a vital function of Automated Driving Systems (ADSs) for providing a safety-guaranteed driving experience. An accident anticipation model aims to predict accidents promptly and accurately before they occur. Existing Artificial Intelligence (AI) models of accident anticipation lack a human-interpretable explanation of their decision-making. Although these models perform well, they remain a black-box to the ADS users, thus difficult to get their trust. To this end, this paper presents a Gated Recurrent Unit (GRU) network that learns spatio-temporal relational features for the early anticipation of traffic accidents from dashcam video data. A post-hoc attention mechanism named Grad-CAM is integrated into the network to generate saliency maps as the visual explanation of the accident anticipation decision. An eye tracker captures human eye fixation points for generating human attention maps. The explainability of network-generated saliency maps is evaluated in comparison to human attention maps. Qualitative and quantitative results on a public crash dataset confirm that the proposed explainable network can anticipate an accident on average 4.57 seconds before it occurs, with 94.02% average precision. In further, various post-hoc attention-based XAI methods are evaluated and compared. It confirms that the Grad-CAM chosen by this study can generate high-quality, human-interpretable saliency maps (with 1.23 Normalized Scanpath Saliency) for explaining the crash anticipation decision. Importantly, results confirm that the proposed AI model, with a human-inspired design, can outperform humans in the accident anticipation.


title: Self Context and Shape Prior for Sensorless Freehand 3D Ultrasound Reconstruction

abstract: 3D ultrasound (US) is widely used for its rich diagnostic information. However, it is criticized for its limited field of view. 3D freehand US reconstruction is promising in addressing the problem by providing broad range and freeform scan. The existing deep learning based methods only focus on the basic cases of skill sequences, and the model relies on the training data heavily. The sequences in real clinical practice are a mix of diverse skills and have complex scanning paths. Besides, deep models should adapt themselves to the testing cases with prior knowledge for better robustness, rather than only fit to the training cases. In this paper, we propose a novel approach to sensorless freehand 3D US reconstruction considering the complex skill sequences. Our contribution is three-fold. First, we advance a novel online learning framework by designing a differentiable reconstruction algorithm. It realizes an end-to-end optimization from section sequences to the reconstructed volume. Second, a self-supervised learning method is developed to explore the context information that reconstructed by the testing data itself, promoting the perception of the model. Third, inspired by the effectiveness of shape prior, we also introduce adversarial training to strengthen the learning of anatomical shape prior in the reconstructed volume. By mining the context and structural cues of the testing data, our online learning methods can drive the model to handle complex skill sequences. Experimental results on developmental dysplasia of the hip US and fetal US datasets show that, our proposed method can outperform the start-of-the-art methods regarding the shift errors and path similarities.


title: Demonstration of Decentralized, Physics-Driven Learning

abstract: In typical artificial neural networks, neurons adjust according to global calculations of a central processor, but in the brain neurons and synapses self-adjust based on local information. Contrastive learning algorithms have recently been proposed to train physical systems, such as fluidic, mechanical, or electrical networks, to perform machine learning tasks from local evolution rules. However, to date such systems have only been implemented in silico due to the engineering challenge of creating elements that autonomously evolve based on their own response to two sets of global boundary conditions. Here we introduce and implement a physics-driven contrastive learning scheme for a network of variable resistors, using circuitry to locally compare the response of two identical networks subjected to the two different sets of boundary conditions. Using this innovation, our system effectively trains itself, optimizing its resistance values without use of a central processor or external information storage. Once the system is trained for a specified allostery, regression, or classification task, the task is subsequently performed rapidly and automatically by the physical imperative to minimize power dissipation in response to the given voltage inputs. We demonstrate that, unlike typical computers, such learning systems are robust to extreme damage (and thus manufacturing defects) due to their decentralized learning. Our twin-network approach is therefore readily scalable to extremely large or nonlinear networks where its distributed nature will be an enormous advantage; a laboratory network of only 500 edges will already outpace its in silico counterpart.


title: Risk Averse Bayesian Reward Learning for Autonomous Navigation from Human Demonstration

abstract: Traditional imitation learning provides a set of methods and algorithms to learn a reward function or policy from expert demonstrations. Learning from demonstration has been shown to be advantageous for navigation tasks as it allows for machine learning non-experts to quickly provide information needed to learn complex traversal behaviors. However, a minimal set of demonstrations is unlikely to capture all relevant information needed to achieve the desired behavior in every possible future operational environment. Due to distributional shift among environments, a robot may encounter features that were rarely or never observed during training for which the appropriate reward value is uncertain, leading to undesired outcomes. This paper proposes a Bayesian technique which quantifies uncertainty over the weights of a linear reward function given a dataset of minimal human demonstrations to operate safely in dynamic environments. This uncertainty is quantified and incorporated into a risk averse set of weights used to generate cost maps for planning. Experiments in a 3-D environment with a simulated robot show that our proposed algorithm enables a robot to avoid dangerous terrain completely in two out of three test scenarios and accumulates a lower amount of risk than related approaches in all scenarios without requiring any additional demonstrations.


title: Holographic thermal correlators revisited

abstract: We study 2-point correlation functions for scalar operators in position space through holography including bulk cubic couplings as well as higher curvature couplings to the square of the Weyl tensor. We focus on scalar operators with large conformal dimensions. This allows us to use the geodesic approximation for propagators. In addition to the leading order contribution, captured by geodesics anchored at the insertion points of the operators on the boundary and probing the bulk geometry thoroughly studied in the literature, the first correction is given by a Witten diagram involving both the bulk cubic coupling and the higher curvature couplings. As a result, this correction is proportional to the VEV of a neutral operator $O_k$ and thus probes the interior of the black hole exactly as in the case studied by Grinberg and Maldacena [13]. The form of the correction matches the general expectations in CFT and allows to identify the contributions of $T^nO_k$ (being $T^n$ the general contraction of n energy-momentum tensors) to the 2-point function. This correction is actually the leading term for off-diagonal correlators (i.e. correlators for operators of different conformal dimension), which can then be computed holographically in this way.


title: Odd-even staggering and shell effects of charge radii for nuclei with even $Z$ from $36$ to $38$ and from $52$ to $62$

abstract: A unified theoretical model reproducing charge radii of known atomic nuclei plays an essential role in making extrapolations for unknown nuclei. Recently developed new ansatz which phenomenologically takes into account the neutron-proton short-range correlations ($np$-SRCs) can describe the discontinuity properties and odd-even staggering (OES) effect of charge radii along isotopic chains remarkably well. In this work, we further review the modified root-mean-square (rms) charge radii formula in the framework of relativistic mean field (RMF) theory. The charge radii are calculated along various isotopic chains that include the nuclei featuring the $N=50$ and $82$ magic shells. Our results suggest that RMF with and without considering a correction term give an almost similar trend of nuclear size for some isotopic chains with open proton shell, especially the abrupt increases across the strong neutron closed shells and the OES behaviors. This reflects that the $np$-SRCs have almost no influence for some nuclei due to the strong coupling between different levels around Fermi surface. The weakening OES behavior of nuclear charge radii is observed generally at completely filled neutron shells and this may be proposed as a signature of magic indicator.


title: A Psychologically Informed Part-of-Speech Analysis of Depression in Social Media

abstract: In this work, we provide an extensive part-of-speech analysis of the discourse of social media users with depression. Research in psychology revealed that depressed users tend to be self-focused, more preoccupied with themselves and ruminate more about their lives and emotions. Our work aims to make use of large-scale datasets and computational methods for a quantitative exploration of discourse. We use the publicly available depression dataset from the Early Risk Prediction on the Internet Workshop (eRisk) 2018 and extract part-of-speech features and several indices based on them. Our results reveal statistically significant differences between the depressed and non-depressed individuals confirming findings from the existing psychology literature. Our work provides insights regarding the way in which depressed individuals are expressing themselves on social media platforms, allowing for better-informed computational models to help monitor and prevent mental illnesses.


title: Vector fields and differential forms on the orbit space of a proper action

abstract: In this paper we study differential forms and vector fields on the orbit space of a proper action of a Lie group on a smooth manifold, defining them as multilinear maps on the generators of infinitesimal diffeomorphisms, respectively. This yields an intrinsic view of vector fields and differential forms on the orbit space.


title: Enhanced Regular Corecursion for Data Streams

abstract: We propose a simple calculus for processing data streams (infinite flows of data series), represented by finite sets of equations built on stream operators. Furthermore, functions defining streams are regularly corecursive, that is, cyclic calls are detected, avoiding non-termination as happens with ordinary recursion in the call-by-value evaluation strategy. As we illustrate by several examples, the combination of such two mechanisms provides a good compromise between expressive power and decidability. Notably, we provide an algorithm to check that the stream returned by a function call is represented by a well-formed set of equations which actually admits a unique solution, hence access to an arbitrary element of the returned stream will never diverge.


title: Parameter-Robust Preconditioning for Oseen Iteration Applied to Stationary and Instationary Navier--Stokes Control

abstract: We derive novel, fast, and parameter-robust preconditioned iterative methods for steady and time-dependent Navier--Stokes control problems. Our approach may be applied to time-dependent problems which are discretized using backward Euler or Crank--Nicolson, and is also a valuable candidate for Stokes control problems discretized using Crank--Nicolson. The key ingredients of the solver are a saddle-point type approximation for the linear systems, an inner iteration for the $(1,1)$-block accelerated by a preconditioner for convection--diffusion control, and an approximation to the Schur complement based on a potent commutator argument applied to an appropriate block matrix. A range of numerical experiments validate the effectiveness of our new approach.


title: Concept Study of a Radio Array Embedded in a Deep Gen2-like Optical Array

abstract: The IceCube Neutrino Observatory has discovered a diffuse astrophysical flux up to 10 PeV and is now planning a large extension with IceCube-Gen2, including an optical array and a large radio array at shallow depth [1]. Neutrino searches for energies >100 PeV are best done with such shallow radio detectors like the Askaryan Radio Array (ARA) or similar (buried as deep as 200 meters below the surface) as they are cheaper to deploy. This poster explores the potential of opportunistically burying radio antennas within the planned IceCube-Gen2 detector volume (between 1350 meters and 2600 meters below the surface). A hybrid detection of events in optical and radio could substantially improve the uncertainty of neutrino cascade direction as radio signals do not scatter in ice. We show the first results of simulating neutrinos from an astrophysical and a cosmogenic flux interacting with 9760 ARA-style vertically polarized radio antennas distributed evenly across 122 strings.


title: The Instrument of the Imaging X-ray Polarimetry Explorer

abstract: While X-ray Spectroscopy, Timing and Imaging have improved verymuch since 1962, when the first astronomical non-solar source was discovered, especially with the launch of Newton/X-ray Multi-Mirror Mission, Rossi/X-ray Timing Explorer and Chandra/Advanced X-ray Astrophysics Facility, the progress of X-ray polarimetry has been meager. This is in part due to the lack of sensitive polarization detectors, in part due to the fate of approved missions and in part because the celestial X-ray sources appeared less polarized than expected. Only one positive measurement has been available until now. Indeed the eight Orbiting Solar Observatory measured the polarization of the Crab nebula in the 70s. The advent of techniques of microelectronics allowed for designing a detector based on the photoelectric effect in gas in an energy range where the optics are efficient in focusing X-rays. Herewe describe the Instrument, which is the major contribution of the Italian collaboration to the SmallExplorer mission called IXPE, the Imaging X-ray Polarimetry Explorer, which will be flown in late 2021. The instrument, is composed of three Detector Units, based on this technique, and a Detector Service Unit. Three Mirror Modules provided by Marshall Space Flight Center focus X-rays onto the detectors. In the following we will show the technological choices, their scientific motivation and the results from the calibration of the Instrument. IXPE will perform imaging, timing and energy resolved polarimetry in the 2-8 keV energy band opening this window of X-ray astronomy to tens of celestial sources of almost all classes.


title: Planning of Power Grasps Using Infinite Program Under Complementary Constraints

abstract: We propose an optimization-based approach to plan power grasps. Central to our method is a reformulation of grasp planning as an infinite program under complementary constraints (IPCC), which allows contacts to happen between arbitrary pairs of points on the object and the robot gripper. We show that IPCC can be reduced to a conventional finite-dimensional nonlinear program (NLP) using a kernel-integral relaxation. Moreover, the values and Jacobian matrices of the kernel-integral can be evaluated efficiently using a modified Fast Multipole Method (FMM). We further guarantee that the planned grasps are collision-free using primal barrier penalties. We demonstrate the effectiveness, robustness, and efficiency of our grasp planner on a row of challenging 3D objects and high-DOF grippers, such as Barrett Hand and Shadow Hand, where our method achieves superior grasp qualities over competitors.


title: Chiral discrimination by polarization singularities of a metal sphere

abstract: The practical applications of chiral discrimination are usually limited by the weak chiral response of enantiomers and the high complexity of detection methods. Here, we propose to use the C lines (i.e., lines of polarization singularities) emerged in light scattering by a metal sphere to detect small chiral particles. Using full-wave numerical simulations and analytical multipole expansions, we determined the absorption dissymmetry of deep-subwavelength helices (i.e., chiral particles) at different positions on the C lines. We find that the absorption dissymmetry can be much larger than that induced by circularly polarized plane wave excitation, and it can be attributed to the asymmetric absorption of the induced electric and magnetic dipoles. We also show that the absorption dissymmetry strongly depends on the anisotropy of the helices. The results may find applications in optical manipulations, optical sensing, and chiral quantum optics.


title: A paucity problem associated with a shifted integer analogue of the divisor function

abstract: Suppose that $\theta$ is irrational. Then almost all elements $\nu\in {\mathbb Z}[\theta]$ that may be written as a $k$-fold product of the shifted integers $n+\theta$ $(n\in {\mathbb N})$ are thus represented essentially uniquely.


title: Distribution of twin primes in repeating sequences of prime factors

abstract: We present a deterministic relationship between relative primes and twin primes in successively larger sequences of the natural numbers. This enables setting a finite lower limit on the occurrence of actual twin primes in an unbounded selection of domains.


title: Proof that Half-Harmonic Oscillators become Full-Harmonic Oscillators after the "Wall Slides Away"

abstract: Normally, the half-harmonic oscillator is active when $x>0$ and absent when $x<0$. From a canonical quantization perspective, this leads to odd eigenfunctions being present while even eigenfunctions are absent. In that case, only the usual odd eigenfunctions will appear if the wall slides to negative infinity. However, if an affine quantization is used, sliding the wall away shows that all the odd and even eigenfunctions are encountered, exactly like any full-harmonic oscillator. We provide numerical support for this.


title: A Hybrid Ensemble Feature Selection Design for Candidate Biomarkers Discovery from Transcriptome Profiles

abstract: The discovery of disease biomarkers from gene expression data has been greatly advanced by feature selection (FS) methods, especially using ensemble FS (EFS) strategies with perturbation at the data level (i.e., homogeneous, Hom-EFS) or method level (i.e., heterogeneous, Het-EFS). Here we proposed a Hybrid EFS (Hyb-EFS) design that explores both types of perturbation to improve the stability and the predictive power of candidate biomarkers. With this, Hyb-EFS aims to disrupt associations of good performance with a single dataset, single algorithm, or a specific combination of both, which is particularly interesting for better reproducibility of genomic biomarkers. We investigated the adequacy of our approach for microarray data related to four types of cancer, carrying out an extensive comparison with other ensemble and single FS approaches. Five FS methods were used in our experiments: Wx, Symmetrical Uncertainty (SU), Gain Ratio (GR), Characteristic Direction (GeoDE), and ReliefF. We observed that the Hyb-EFS and Het-EFS approaches attenuated the large performance variation observed for most single FS and Hom-EFS across distinct datasets. Also, the Hyb-EFS improved upon the stability of the Het-EFS within our domain. Comparing the Hyb-EFS and Het-EFS composed of the top-performing selectors (Wx, GR, and SU), our hybrid approach surpassed the equivalent heterogeneous design and the best Hom-EFS (Hom-Wx). Interestingly, the rankings produced by our Hyb-EFS reached greater biological plausibility, with a notably high enrichment for cancer-related genes and pathways. Thus, our experiments suggest the potential of the proposed Hybrid EFS design in discovering candidate biomarkers from microarray data. Finally, we provide an open-source framework to support similar analyses in other domains, both as a user-friendly application and a plain Python package.


title: Modeling and Design of IRS-Assisted Multi-Link FSO Systems

abstract: In this paper, we investigate the modeling and design of intelligent reflecting surface (IRS)-assisted optical communication systems which are deployed to relax the line-of-sight (LOS) requirement in multi-link free space optical (FSO) systems. The FSO laser beams incident on the optical IRSs have a Gaussian power intensity profile and a nonlinear phase profile, whereas the plane waves in radio frequency (RF) systems have a uniform power intensity profile and a linear phase profile. Given these substantial differences, the results available for IRS-assisted RF systems are not applicable to IRS-assisted FSO systems. Therefore, we develop a new analytical channel model for point-to-point IRS-assisted FSO systems based on the Huygens-Fresnel principle. Our analytical model captures the impact of the size, position, and orientation of the IRS as well as its phase shift profile on the end-to-end channel. To allow the sharing of the optical IRS by multiple FSO links, we propose three different protocols, namely the time division (TD), IRS-division (IRSD), and IRS homogenization (IRSH) protocols. The proposed protocols address the specific characteristics of FSO systems including the non-uniformity and possible misalignment of the laser beams. Furthermore, to compare the proposed IRS sharing protocols, we analyze the bit error rate (BER) and the outage probability of IRS-assisted multi-link FSO systems in the presence of inter-link interference. Our simulation results validate the accuracy of the proposed analytical channel model for IRS-assisted FSO systems and confirm that this model is applicable for both large and intermediate IRS-receiver lens distances. Moreover, in the absence of misalignment errors, the IRSD protocol outperforms the other protocols, whereas in the presence of misalignment errors, the IRSH protocol performs significantly better than the IRSD protocol.


title: Design Simulation of Czerny-Turner Configuration-based Raman Spectrometer using physical optics propagation algorithm

abstract: We report the design simulation of the Raman spectrometer using optical system design software Zemax. The design is based on the Czerny-turner configuration which includes an optical system consisting of an entrance slit, two concave mirrors, reflecting type diffraction grating, and an image detector. The system's modelling approach is suggested by introducing the corresponding relationship between detector pixels and wavelength, linear CCD receiving surface length, and image surface dimension. Simulations have been carried out using the POP (Physical Optics Propagation) algorithm. Spot diagram, relative illumination, irradiance plot, Modulation Transfer Function (MTF), geometric, and encircled energy simulated for designing the Raman spectrometer. The simulation results of the Raman spectrometer's using a 527 nm wavelength laser as an excitation light source are presented. The present optical system is designed in sequential mode and a Raman spectrum observed in the range of 530 to 630 nm. The analysis shows that the system's image efficiency is higher, predicting that it is possible to build an efficient and cost-effective Raman spectrometer for optical diagnostics.


title: Inverse Reinforcement Learning for Strategy Identification

abstract: In adversarial environments, one side could gain an advantage by identifying the opponent's strategy. For example, in combat games, if an opponents strategy is identified as overly aggressive, one could lay a trap that exploits the opponent's aggressive nature. However, an opponent's strategy is not always apparent and may need to be estimated from observations of their actions. This paper proposes to use inverse reinforcement learning (IRL) to identify strategies in adversarial environments. Specifically, the contributions of this work are 1) the demonstration of this concept on gaming combat data generated from three pre-defined strategies and 2) the framework for using IRL to achieve strategy identification. The numerical experiments demonstrate that the recovered rewards can be identified using a variety of techniques. In this paper, the recovered reward are visually displayed, clustered using unsupervised learning, and classified using a supervised learner.


title: Scaling Theory of 3D Magnetic Reconnection Spreading

abstract: We develop a first-principles scaling theory of the spreading of three-dimensional (3D) magnetic reconnection of finite extent in the out of plane direction. This theory addresses systems with or without an out of plane (guide) magnetic field, and with or without Hall physics. The theory reproduces known spreading speeds and directions with and without guide fields, unifying previous knowledge in a single theory. New results include: (1) Reconnection spreads in a particular direction if an x-line is induced at the interface between reconnecting and non-reconnecting regions, which is controlled by the out of plane gradient of the electric field in the outflow direction. (2) The spreading mechanism for anti-parallel collisionless reconnection is convection, as is known, but for guide field reconnection it is magnetic field bending. We confirm the theory using 3D two-fluid and resistive-magnetohydrodynamics simulations. (3) The theory explains why anti-parallel reconnection in resistive-magnetohydrodynamics does not spread. (4) The simulation domain aspect ratio, associated with the free magnetic energy, influences whether reconnection spreads or convects with a fixed x-line length. (5) We perform a simulation initiating anti-parallel collisionless reconnection with a pressure pulse instead of a magnetic perturbation, finding spreading is unchanged rather than spreading at the magnetosonic speed as previously suggested. The results provide a theoretical framework for understanding spreading beyond systems studied here, and are important for applications including two-ribbon solar flares and reconnection in Earth's magnetosphere.


title: Fair Representation Learning using Interpolation Enabled Disentanglement

abstract: With the growing interest in the machine learning community to solve real-world problems, it has become crucial to uncover the hidden reasoning behind their decisions by focusing on the fairness and auditing the predictions made by these black-box models. In this paper, we propose a novel method to address two key issues: (a) Can we simultaneously learn fair disentangled representations while ensuring the utility of the learned representation for downstream tasks, and (b)Can we provide theoretical insights into when the proposed approach will be both fair and accurate. To address the former, we propose the method FRIED, Fair Representation learning using Interpolation Enabled Disentanglement. In our architecture, by imposing a critic-based adversarial framework, we enforce the interpolated points in the latent space to be more realistic. This helps in capturing the data manifold effectively and enhances the utility of the learned representation for downstream prediction tasks. We address the latter question by developing a theory on fairness-accuracy trade-offs using classifier-based conditional mutual information estimation. We demonstrate the effectiveness of FRIED on datasets of different modalities - tabular, text, and image datasets. We observe that the representations learned by FRIED are overall fairer in comparison to existing baselines and also accurate for downstream prediction tasks. Additionally, we evaluate FRIED on a real-world healthcare claims dataset where we conduct an expert aided model auditing study providing useful insights into opioid ad-diction patterns.


title: Deep Predictive Learning of Carotid Stenosis Severity

abstract: Carotid artery stenosis is the narrowing of carotid arteries, which supplies blood to the neck and head. In this work, we train a model to predict the severity of the stenosis blockage based on SRUC criteria variables and other patient information. We implement classic machine learning methods, decision trees and random forests, used in a previous experiment. In addition, we improve the accuracy through the use of the state-of-art Augmented Neural ODE deep learning method. Through systematic and theory-rooted analysis, we examine different parameters to achieve an accuracy of about 77%. These results show the strong potential in applying recently developing deep learning methods, while simultaneously suggesting that the current data provided by the SRUC criteria may be insufficient to predict stenosis severity at a high performance level.


title: A Multifunctional Sub-10nm Transistor

abstract: Nano-electronic integrated circuit technology is exclusively based on MOSFET transistor due to its scalability down to the nanometer range. On the other hand, Bipolar Junction Transistor (BJT), which provides unmatched analog characteristics and frequency response, cannot be scaled to nanometer regime without the loss of transistor action. Here a versatile nanoscale transistor is introduced that provides identical BJT behavior and expands its capabilities. The new transistor uses CMOS fabrication technology and creates BJT emitter, base, and collector via electric fields. By allowing carrier modulation during operation, its current gain can be changed at least by five orders of magnitude. This property introduces novel adaptive, variable gain, and programmable analog modules into existing electronic circuit design and manufacturing. A NOT gate version of this device with the critical dimension of 7 nm operates at 730 GHz, and its three-stage ring oscillator exhibits a frequency of 240 GHz. With proper gate biasing, it can also operate as a nanoscale MOSFET, easily alleviating short-channel effects.


title: Filling the G_ap_s: Multivariate Time Series Imputation by Graph Neural Networks

abstract: Dealing with missing values and incomplete time series is a labor-intensive, tedious, inevitable task when handling data coming from real-world applications. Effective spatio-temporal representations would allow imputation methods to reconstruct missing temporal data by exploiting information coming from sensors at different locations. However, standard methods fall short in capturing the nonlinear time and space dependencies existing within networks of interconnected sensors and do not take full advantage of the available - and often strong - relational information. Notably, most state-of-the-art imputation methods based on deep learning do not explicitly model relational aspects and, in any case, do not exploit processing frameworks able to adequately represent structured spatio-temporal data. Conversely, graph neural networks have recently surged in popularity as both expressive and scalable tools for processing sequential data with relational inductive biases. In this work, we present the first assessment of graph neural networks in the context of multivariate time series imputation. In particular, we introduce a novel graph neural network architecture, named GRIN, which aims at reconstructing missing data in the different channels of a multivariate time series by learning spatio-temporal representations through message passing. Empirical results show that our model outperforms state-of-the-art methods in the imputation task on relevant real-world benchmarks with mean absolute error improvements often higher than 20%.


title: Coupled and uncoupled sign-changing spikes of singularly perturbed elliptic systems

abstract: We study the existence and asymptotic behavior of solutions having positive and sign-changing components to the singularly perturbed system of elliptic equations \begin{equation*} \begin{cases} -\varepsilon^2\Delta u_i+u_i=\mu_i|u_i|^{p-2}u_i + \sum\limits_{\substack{j=1 \\ j \not=i}}^\ell\lambda_{ij}\beta_{ij}|u_j|^{\alpha_{ij}}|u_i|^{\beta_{ij} -2}u_i,\\ u_i \in H^1_0(\Omega), \quad u_i\neq 0, \qquad i=1,\ldots,\ell, \end{cases} \end{equation*} in a bounded domain $\Omega$ in $\mathbb{R}^N$, with $N\geq 4$, $\varepsilon>0$, $\mu_i>0$, $\lambda_{ij}=\lambda_{ji}<0$, $\alpha_{ij}, \beta_{ij}>1$, $\alpha_{ij}=\beta_{ji}$, $\alpha_{ij} + \beta_{ij} = p\in (2,2^*)$, and $2^{*}:=\frac{2N}{N-2}$. If $\Omega$ is the unit ball we obtain solutions with a prescribed combination of positive and nonradial sign-changing components exhibiting two different types of asymptotic behavior as $\varepsilon\to 0$: solutions whose limit profile is a rescaling of a solution with positive and nonradial sign-changing components of the limit system \begin{equation*} \begin{cases} -\Delta u_i+u_i=\mu_i|u_i|^{p-2}u_i + \sum\limits_{\substack{j=1 \\ j \not=i}}^\ell\lambda_{ij}\beta_{ij}|u_j|^{\alpha_{ij}}|u_i|^{\beta_{ij} -2}u_i,\\ u_i \in H^1(\mathbb{R}^N), \quad u_i\neq 0, \qquad i=1,\ldots,\ell, \end{cases} \end{equation*} and solutions whose limit profile is a solution of the uncoupled system, i.e., after rescaling and translation, the limit profile of the $i$-th component is a positive or a nonradial sign-changing solution to the equation $$-\Delta u+u=\mu_i|u|^{p-2}u,\qquad u \in H^1(\mathbb{R}^N), \qquad u\neq 0.$$


title: Universal bounds for fixed point iterations via optimal transport metrics

abstract: We present a self-contained analysis of a particular family of metrics over the set of non-negative integers. We show that these metrics, which are defined through a nested sequence of optimal transport problems, provide tight estimates for general Krasnosel'skii-Mann fixed point iterations for non-expansive maps. We also describe some of their very special properties, including their monotonicity and the so-called "convex quadrangle inequality" that yields a greedy algorithm to compute them efficiently.


title: Improving Grasp Stability with Rotation Measurement from Tactile Sensing

abstract: Rotational displacement about the grasping point is a common grasp failure when an object is grasped at a location away from its center of gravity. Tactile sensors with soft surfaces, such as GelSight sensors, can detect the rotation patterns on the contacting surfaces when the object rotates. In this work, we propose a model-based algorithm that detects those rotational patterns and measures rotational displacement using the GelSight sensor. We also integrate the rotation detection feedback into a closed-loop regrasping framework, which detects the rotational failure of grasp in an early stage and drives the robot to a stable grasp pose. We validate our proposed rotation detection algorithm and grasp-regrasp system on self-collected dataset and online experiments to show how our approach accurately detects the rotation and increases grasp stability.


title: Conditional Bures Metric for Domain Adaptation

abstract: As a vital problem in classification-oriented transfer, unsupervised domain adaptation (UDA) has attracted widespread attention in recent years. Previous UDA methods assume the marginal distributions of different domains are shifted while ignoring the discriminant information in the label distributions. This leads to classification performance degeneration in real applications. In this work, we focus on the conditional distribution shift problem which is of great concern to current conditional invariant models. We aim to seek a kernel covariance embedding for conditional distribution which remains yet unexplored. Theoretically, we propose the Conditional Kernel Bures (CKB) metric for characterizing conditional distribution discrepancy, and derive an empirical estimation for the CKB metric without introducing the implicit kernel feature map. It provides an interpretable approach to understand the knowledge transfer mechanism. The established consistency theory of the empirical estimation provides a theoretical guarantee for convergence. A conditional distribution matching network is proposed to learn the conditional invariant and discriminative features for UDA. Extensive experiments and analysis show the superiority of our proposed model.


title: Practical Adoption of Cloud Computing in Power Systems- Drivers, Challenges, Guidance, and Real-world Use Cases

abstract: Motivated by The Federal Energy Regulatory Commission's (FERC) recent direction and ever-growing interest in cloud adoption by power utilities, a Task Force was established to assist power system practitioners with secure, reliable and cost-effective adoption of cloud technology to meet various business needs. This paper summarizes the business drivers, challenges, guidance, and best practices for cloud adoption in power systems from the Task Force's perspective, after extensive review and deliberation by its members, including grid operators, utility companies, software vendors, and cloud providers. The paper begins by enumerating various business drivers for cloud adoption in the power industry. It follows with the discussion of the challenges and risks of migrating power grid utility workloads to the cloud. Next, for each corresponding challenge or risk, the paper provides appropriate guidance. Notably, the guidance is directed toward power industry professionals who are considering cloud solutions and are yet hesitant about the practical execution. Finally, to tie all the sections together, the paper documents various real-world use cases of cloud technology in the power system domain, which both the power industry practitioners and software vendors can look toward to design and select their own future cloud solutions. We hope that the information in this paper will serve as helpful guidance for the development of NERC guidelines and standards relevant to cloud adoption in the industry.


title: High-precision mapping of diamond crystal strain using quantum interferometry

abstract: Crystal strain variation imposes significant limitations on many quantum sensing and information applications for solid-state defect qubits in diamond. Thus, precision measurement and control of diamond crystal strain is a key challenge. Here, we report diamond strain measurements with a unique set of capabilities, including micron-scale spatial resolution, millimeter-scale field-of-view, and a two order-of-magnitude improvement in volume-normalized sensitivity over previous work [1], reaching $5(2) \times 10^{-8}/\sqrt{\rm{Hz}\cdot\rm{\mu m}^3}$ (with spin-strain coupling coefficients representing the dominant systematic uncertainty). We use strain-sensitive spin-state interferometry on ensembles of nitrogen vacancy (NV) color centers in single-crystal CVD bulk diamond with low strain gradients. This quantum interferometry technique provides insensitivity to magnetic-field inhomogeneity from the electronic and nuclear spin bath, thereby enabling long NV ensemble electronic spin dephasing times and enhanced strain sensitivity. We demonstrate the strain-sensitive measurement protocol first on a scanning confocal laser microscope, providing quantitative measurement of sensitivity as well as three-dimensional strain mapping; and second on a wide-field imaging quantum diamond microscope (QDM). Our strain microscopy technique enables fast, sensitive characterization for diamond material engineering and nanofabrication; as well as diamond-based sensing of strains applied externally, as in diamond anvil cells or embedded diamond stress sensors, or internally, as by crystal damage due to particle-induced nuclear recoils.

categories: quant-ph cond-mat.mtrl-sci hep-ex physics.app-ph physics.ins-det


title: QCD factorization for hadronic quarkonium production at high $p_T$

abstract: Heavy quarkonium production at high transverse momentum ($p_T$) in hadronic collisions is explored in the QCD factorization approach. We find that the leading power in the $1/p_T$ expansion is responsible for high $p_T$ regime, while the next-to-leading power contribution is necessary for the low $p_T$ region. We present the first numerical analysis of the scale evolution of coupled twist-2 and twist-4 fragmentation functions (FFs) for heavy quarkonium production and demonstrate that the QCD factorization approach is capable of describing the $p_T$ spectrum of hadronic $J/\psi$ production at the LHC.


title: A graphical multi-fidelity Gaussian process model, with application to emulation of heavy-ion collisions

abstract: With advances in scientific computing and mathematical modeling, complex scientific phenomena such as galaxy formations and rocket propulsion can now be reliably simulated. Such simulations can however be very time-intensive, requiring millions of CPU hours to perform. One solution is multi-fidelity emulation, which uses data of different fidelities to train an efficient predictive model which emulates the expensive simulator. For complex scientific problems and with careful elicitation from scientists, such multi-fidelity data may often be linked by a directed acyclic graph (DAG) representing its scientific model dependencies. We thus propose a new Graphical Multi-fidelity Gaussian Process (GMGP) model, which embeds this DAG structure (capturing scientific dependencies) within a Gaussian process framework. We show that the GMGP has desirable modeling traits via two Markov properties, and admits a scalable algorithm for recursive computation of the posterior mean and variance along at each depth level of the DAG. We also present a novel experimental design methodology over the DAG given an experimental budget, and propose a nonlinear extension of the GMGP via deep Gaussian processes. The advantages of the GMGP are then demonstrated via a suite of numerical experiments and an application to emulation of heavy-ion collisions, which can be used to study the conditions of matter in the Universe shortly after the Big Bang. The proposed model has broader uses in data fusion applications with graphical structure, which we further discuss.


title: Quasiperiodicity and blowup in integrable subsystems of nonconservative nonlinear Schr\"odinger equations

abstract: In this paper, we study the dynamics of a class of nonlinear Schr\"odinger equation $ i u_t = \triangle u + u^p $ for $ x \in \mathbb{T}^d$. We prove that the PDE is integrable on the space of non-negative Fourier coefficients, in particular that each Fourier coefficient of a solution can be explicitly solved by quadrature. Within this subspace we demonstrate a large class of (quasi)periodic solutions all with the same frequency, as well as solutions which blowup in finite time in the $L^2$ norm.


title: Human Evaluation of Creative NLG Systems: An Interdisciplinary Survey on Recent Papers

abstract: We survey human evaluation in papers presenting work on creative natural language generation that have been published in INLG 2020 and ICCC 2020. The most typical human evaluation method is a scaled survey, typically on a 5 point scale, while many other less common methods exist. The most commonly evaluated parameters are meaning, syntactic correctness, novelty, relevance and emotional value, among many others. Our guidelines for future evaluation include clearly defining the goal of the generative system, asking questions as concrete as possible, testing the evaluation setup, using multiple different evaluation setups, reporting the entire evaluation process and potential biases clearly, and finally analyzing the evaluation results in a more profound way than merely reporting the most typical statistics.


title: Motion Planning for Variable Topology Trusses: Reconfiguration and Locomotion

abstract: Truss robots are highly redundant parallel robotic systems that can be applied in a variety of scenarios. The variable topology truss (VTT) is a class of modular truss robots. As self-reconfigurable modular robots, a VTT is composed of many edge modules that can be rearranged into various structures depending on the task. These robots change their shape by not only controlling joint positions as with fixed morphology robots, but also reconfiguring the connectivity between truss members in order to change their topology. The motion planning problem for VTT robots is difficult due to their varying morphology, high dimensionality, the high likelihood for self-collision, and complex motion constraints. In this paper, a new motion planning framework to dramatically alter the structure of a VTT is presented. It can also be used to solve locomotion tasks that are much more efficient compared with previous work. Several test scenarios are used to show its effectiveness. Supplementary materials are available at https://www.modlabupenn.org/vtt-motion-planning/.


title: Observation of complex functional cortical patterns in brain cognition

abstract: Synaptic plasticity and neuron cross-talk are some of the important key mechanisms underlying formation of dynamic clusters of active neurons. The essence of this study is to model and decipher the mechanism of emergence of a task-specific functional connectivity among a population of neurons. We have used our proposed neuron activity pattern (NAP) model to define an assorted range of interactions and simulated across a 2D lattice with two-state (firing/non-firing) neurons. We observed the emergence of scale-free, hierarchically organized ordered patterns of active neurons, that we call as functional cortical patterns (FCPs), only at the near-critical phase transition. We have done extensive topological characterization of FCPs and tested its congruency with the functional brain networks obtained from the empirical electrophysiological data of a visual stimulus task. Our results of network theory attributes, fractal analysis and functional cartography have confirmed the transitioning effect while generation of FCPs. This gives us strong implication that the functional neuronal system supports far-from-equilibrium dynamics upon receiving a stimulus. Our investigation for the long, critical and short-ranged interactions has lead to an interpretation that near-critical transition phase is a spectrum of critical temperatures (coupling strengths) being affected by the distance between the neurons. The range of interaction among neurons plays an important role in defining the functional severity of the neuronal circuitry. We propose that interaction range as a function of coupling strength drives connectivity in an augmenting fashion as we observed an increase in coupling strength from short to long-range interactions, among the population of neurons. This gives us insight towards the intensity of cognitive behaviour, which is attained by multiple stimulus attempts.


title: Scrutinizing $GW$-based methods using the Hubbard dimer

abstract: Using the simple (symmetric) Hubbard dimer, we analyze some important features of the $GW$ approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot $GW$ method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard $GW$ approximation and find, in particular, that i) some neutral excitation energies become complex when the electron-electron interaction $U$ increases, which can be traced back to the approximate nature of the $GW$ quasiparticle energies; ii) the BSE formalism yields accurate correlation energies over a wide range of $U$ when the trace (or plasmon) formula is employed; iii) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); iv) the trace formula has the correct behavior for weak (\ie, small $U$) interaction, unlike the ACFDT expression.

categories: physics.chem-ph cond-mat.mtrl-sci cond-mat.str-el physics.comp-ph


title: Realizing normal group-velocity dispersion in free space via angular dispersion

abstract: It has long been thought that normal group-velocity dispersion (GVD) cannot be produced in free space via angular dispersion. Indeed, conventional diffractive or dispersive components such as gratings or prisms produce only anomalous GVD. We identify the conditions that must be fulfilled by the angular dispersion introduced into a plane-wave pulse to yield normal GVD. We then utilize a pulsed-beam shaper capable of introducing arbitrary angular-dispersion profiles to symmetrically produce normal and anomalous GVD in free space, which are realized here on the same footing for the first time.


title: Adaptive interface-Mesh un-Refinement (AiMuR) based Sharp-Interface Level-Set-Method for Two-Phase Flow

abstract: Adaptive interface-Mesh un-Refinement (AiMuR) based Sharp-Interface Level-Set-Method (SI-LSM) is proposed for both uniform and non-uniform Cartesian-Grid. The AiMuR involves interface location based dynamic un-refinement (with merging of the four control volumes) of the Cartesian grid away from the interface. The un-refinement is proposed for the interface solver only. A detailed numerical methodology is presented for the AiMuR and ghost-fluid method based SI-LSM. Advantage of the novel as compared to the traditional SI-LSM is demonstrated with a detailed qualitative as well as quantitative performance study, involving the SI-LSMs on both coarse grid and fine grid, for three sufficiently different two-phase flow problems: dam break, breakup of a liquid jet and drop coalescence. A superior performance of AiMuR based SI-LSM is demonstrated - the AiMuR on a coarser non-uniform grid ($NU_{c}^{AiMuR}$) is almost as accurate as the traditional SI-LSM on a uniform fine grid ($U_{f}$) and takes a computational time almost same as that by the traditional SI-LSM on a uniform coarse grid ($U_{c}$). The AMuR is different from the existing Adaptive Mesh Refinement (AMR) as the former involves only mesh un-refinement while the later involves both refinement and un-refinement of the mesh. Moreover, the proposed computational development is significant since the present adaptive un-refinement strategy is much simpler to implement as compared to that for the commonly used adaptive refinement strategies. The proposed numerical development can be extended to various other multi-physics, multi-disciplinary and multi-scale problems involving interfaces.


title: Iterative Deliberation via Metric Aggregation

abstract: We investigate an iterative deliberation process for an agent community wishing to make a joint decision. We develop a general model consisting of a community of n agents, each with their initial ideal point in some metric space (X, d), such that in each iteration of the iterative deliberation process, all agents move slightly closer to the current winner, according to some voting rule R. For several natural metric spaces and suitable voting rules for them, we identify conditions under which such an iterative deliberation process is guaranteed to converge.


title: Exotic vortices and twisted holomorphic maps

abstract: We consider the exotic vortex equations on compact Riemann surfaces. These generalise the well-known Jackiw-Pi and Ambj{\o}rn-Olesen vortex equations and arise as equations for Bogomolny-Prasad-Sommerfield-like configurations in nonrelativistic Chern-Simons-matter theories. We show that (exotic) vortex solutions in $U(N)$ gauge theories with a single fundamental flavour on compact Riemann surfaces can be constructed from gauged holomorphic maps into complex space forms of dimension $N$. We conjecture, with some evidence, that all such solutions arise this way. This leads us to insights regarding the moduli theory of exotic vortices, including the identification of interesting selection rules for nonAbelian exotic vortices on the sphere. We also consider the classical and quantum dynamics of exotic vortices. Using recent localisation results, we give the Witten index of the quantum mechanics describing low-temperature Abelian exotic vortex dynamics.


title: Chest ImaGenome Dataset for Clinical Reasoning

abstract: Despite the progress in automatic detection of radiologic findings from chest X-ray (CXR) images in recent years, a quantitative evaluation of the explainability of these models is hampered by the lack of locally labeled datasets for different findings. With the exception of a few expert-labeled small-scale datasets for specific findings, such as pneumonia and pneumothorax, most of the CXR deep learning models to date are trained on global "weak" labels extracted from text reports, or trained via a joint image and unstructured text learning strategy. Inspired by the Visual Genome effort in the computer vision community, we constructed the first Chest ImaGenome dataset with a scene graph data structure to describe $242,072$ images. Local annotations are automatically produced using a joint rule-based natural language processing (NLP) and atlas-based bounding box detection pipeline. Through a radiologist constructed CXR ontology, the annotations for each CXR are connected as an anatomy-centered scene graph, useful for image-level reasoning and multimodal fusion applications. Overall, we provide: i) $1,256$ combinations of relation annotations between $29$ CXR anatomical locations (objects with bounding box coordinates) and their attributes, structured as a scene graph per image, ii) over $670,000$ localized comparison relations (for improved, worsened, or no change) between the anatomical locations across sequential exams, as well as ii) a manually annotated gold standard scene graph dataset from $500$ unique patients.


title: A non-renormalizable neutrino mass model with $\mathbf{S}_{3}\otimes \mathbf{Z}_{2}$ symmetry

abstract: The lepton sector is studied within a flavored non-renormalizable model where the $\mathbf{S}_{3}\otimes \mathbf{Z}_{2}$ flavor symmetry drives the Yukawa couplings. In this framework, the effective neutrino mass, that comes from the type II see-saw mechanism, as well as the charged lepton mass matrices are hierarchical and possessed (under a benchmark in the charged sector) a kind of Fritzsch textures that achieve to accommodate the mixing angles in good agreement with the last experimental data. The model predicts a CP-violating phase consistent with the experimental data and the branching ratio for the lepton flavor violation process, $\mu\rightarrow e\gamma$, is well below the current bound.


title: Statistical learning method for predicting density-matrix based electron dynamics

abstract: We develop a statistical method to learn a molecular Hamiltonian matrix from a time-series of electron density matrices. We extend our previous method to larger molecular systems by incorporating physical properties to reduce dimensionality, while also exploiting regularization techniques like ridge regression for addressing multicollinearity. With the learned Hamiltonian we can solve the Time-Dependent Hartree-Fock (TDHF) equation to propagate the electron density in time, and predict its dynamics for field-free and field-on scenarios. We observe close quantitative agreement between the predicted dynamics and ground truth for both field-off trajectories similar to the training data, and field-on trajectories outside of the training data.


title: Less is more: balancing noise reduction and data retention in fMRI with data-driven scrubbing

abstract: Artifacts in functional MRI (fMRI) data cause deviations from common distributional assumptions, introduce spatial and temporal outliers, and reduce the signal-to-noise ratio of the data -- all of which can have negative consequences for downstream statistical analysis. Scrubbing is a technique for excluding fMRI volumes thought to be contaminated by artifacts and generally comes in two flavors. Motion scrubbing based on subject head motion-derived measures is popular but suffers from a number of drawbacks, especially high rates of censoring of individual volumes and entire subjects. Alternatively, data-driven scrubbing methods like DVARS are based on observed noise in the processed fMRI timeseries and may avoid some of these issues. Here we propose "projection scrubbing", a novel data-driven scrubbing method based on a statistical outlier detection framework and strategic dimension reduction, including independent component analysis (ICA), to isolate artifactual variation. We undertake a comprehensive comparison of motion scrubbing with data-driven projection scrubbing and DVARS. We argue that an appropriate metric for the success of scrubbing is maximal data retention subject to reasonable performance on typical benchmarks of functional connectivity. We find that stringent motion scrubbing yields worsened validity, worsened reliability, and produced small improvements to fingerprinting. Meanwhile, data-driven scrubbing methods tend to yield greater improvements to fingerprinting while not generally worsening validity or reliability. Importantly, however, data-driven scrubbing excludes a fraction of the number of volumes or entire sessions compared to motion scrubbing. The ability of data-driven fMRI scrubbing to improve data retention without negatively impacting the quality of downstream analysis has major implications for sample sizes in population neuroscience research.


title: StudyMe: A New Mobile App for User-Centric N-of-1 Trials

abstract: N-of-1 trials are multi-crossover self-experiments that allow individuals to systematically evaluate the effect of interventions on their personal health goals. Although several tools for N-of-1 trials exist, none support non-experts in conducting their own user-centric trials. In this study we present StudyMe, an open-source mobile application that is freely available from https://play.google.com/store/apps/details?id=health.studyu.me and offers users flexibility and guidance in configuring every component of their trials. We also present research that informed the development of StudyMe. Through an initial survey with 272 participants, we learned that individuals are interested in a variety of personal health aspects and have unique ideas on how to improve them. In an iterative, user-centered development process with intermediate user tests we developed StudyMe that also features an educational part to communicate N-of-1 trial concepts. A final empirical evaluation of StudyMe showed that all participants were able to create their own trials successfully using StudyMe and the app achieved a very good usability rating. Our findings suggest that StudyMe provides a significant step towards enabling individuals to apply a systematic science-oriented approach to personalize health-related interventions and behavior modifications in their everyday lives.


title: Accurate full configuration interaction correlation energy estimates for five- and six-membered rings

abstract: Following our recent work on the benzene molecule [\href{https://doi.org/10.1063/5.0027617}{J.~Chem.~Phys.~\textbf{153}, 176101 (2020)}], itself motivated by the blind challenge of Eriksen \textit{et al.} [\href{https://doi.org/10.1021/acs.jpclett.0c02621}{J.~Phys.~Chem.~Lett.~\textbf{11}, 8922 (2020)}] on the same system, we report accurate full configuration interaction (FCI) frozen-core correlation energy estimates for twelve five- and six-membered ring molecules in the standard correlation-consistent double-$\zeta$ Dunning basis set (cc-pVDZ). Our FCI correlation energy estimates, with estimated error smaller than 1 millihartree, are based on energetically optimized-orbital selected configuration interaction (SCI) calculations performed with the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) algorithm. Having at our disposal these accurate reference energies, the respective performance and convergence properties of several popular and widely-used families of single-reference quantum chemistry methods are investigated. In particular, we study the convergence properties of i) the M{\o}ller-Plesset perturbation series up to fifth-order (MP2, MP3, MP4, and MP5), ii) the iterative approximate coupled-cluster series CC2, CC3, and CC4, and iii) the coupled-cluster series CCSD, CCSDT, and CCSDTQ. The performance of the ground-state gold standard CCSD(T) as well as the completely renormalized CC model, CR-CC(2,3), are also investigated.

categories: physics.chem-ph cond-mat.str-el nucl-th physics.comp-ph quant-ph


title: Evolution of Saturn's north polar color and cloud structure between 2012 and 2017 inferred from Cassini VIMS and ISS observations

abstract: Cassini/ISS imagery and Cassini/VIMS spectral imaging observations from 0.35 to 5.12 microns show that between 2012 and 2017 the region poleward of the Saturn's northern hexagon changed from dark blue/green to a moderately brighter gold color, except for the inner eye region (88.2 deg - 90 deg N), which remained relatively unchanged. These and even more dramatic near-IR changes can be reproduced by an aerosol model of four compact layers consisting of a stratospheric haze at an effective pressure near 50 mbar, a deeper haze of putative diphosphine particles typically near 300 mbar, an ammonia cloud layer with a base pressure between 0.4 bar and 1.3 bar, and a deeper cloud of a possible mix of NH4SH and water ice particles within the 2.7 to 4.5 bar region. Our analysis of the background clouds between the discrete features shows that between 2013 and 2016 the effective pressures of most layers changed very little, except for the ammonia ice layer, which decreased from about 1 bar to 0.4 bar near the edge of the eye, but increased to 1 bar inside the eye. Inside the hexagon there were large increases in optical depth, by up to a factor of 10 near the eye for the putative diphosphine layer and by a factor of four over most of the hexagon interior. Inside the eye, aerosol optical depths were very low, suggesting downwelling motions. The high contrast between eye and surroundings in 2016 was due to substantial increases in optical depths outside the eye. The color change from blue/green to gold inside most of the hexagon region can be explained in our model almost entirely by changes in the stratospheric haze, which increased between 2013 and 2016 by a factor of four in optical depth and by almost a factor of three in the short-wavelength peak imaginary index.