title: SomnNET: An SpO2 Based Deep Learning Network for Sleep Apnea Detection in Smartwatches

abstract: The abnormal pause or rate reduction in breathing is known as the sleep-apnea hypopnea syndrome and affects the quality of sleep of an individual. A novel method for the detection of sleep apnea events (pause in breathing) from peripheral oxygen saturation (SpO2) signals obtained from wearable devices is discussed in this paper. The paper details an apnea detection algorithm of a very high resolution on a per-second basis for which a 1-dimensional convolutional neural network -- which we termed SomnNET -- is developed. This network exhibits an accuracy of 97.08% and outperforms several lower resolution state-of-the-art apnea detection methods. The feasibility of model pruning and binarization to reduce the computational complexity is explored. The pruned network with 80% sparsity exhibited an accuracy of 89.75%, and the binarized network exhibited an accuracy of 68.22%. The performance of the proposed networks is compared against several state-of-the-art algorithms.


title: Proceedings of the Advanced Accelerator Concepts Seminar Series 2020

abstract: The Advanced Accelerator Concepts (AAC) Seminar Series 2020 (https://aacseminarseries.lbl.gov/), organized and hosted by the Lawrence Berkeley National Laboratory, consisted of nine weekly sessions, each one dedicated to a particular topic of interest within the framework of advanced accelerator concepts research. The Seminar Series was a fully-remote event that provided a forum for the advanced accelerator community. The AAC Seminar Series was held in lieu of the AAC 2020 Workshop (https://aac2020.lbl.gov/), originally planned for June 2020 and canceled due to the COVID-19 pandemic. Since its inception in 1982, the biennial AAC Workshop has become the principal US and international meeting for advanced particle accelerator research and development.


title: A Bayesian Surveillance Model to Track Variable Rainfall-Runoff Responses for Small Watersheds

abstract: Understanding dynamics of hydrological responses is essential in producing skillful runoff forecast. This can be quantitatively done by tracking changes in hydrology model parameters that represent physical characteristics. In this study, we implement a Bayesian estimation method in continuously estimating hydrology model parameters given observations of rainfall and runoff for small watersheds. The method is coupled with a conceptual hydrology model using a Gamma distribution-based Instantaneous Unit Hydrograph. The whole analytical framework is tested using synthetic data as well as observational data from the Fall Creek watershed. The results show that the Bayesian method can well track the hidden parameters that change inter-annually. Then the model is applied to examine temporal and spatial variability of the rainfall-runoff responses and we find 1) a systematic shift in the rainfall-runoff response for the Fall Creek watershed around 1943 and 2) a statistically significant relationship between rainfall-runoff responses and watershed sizes for selected NY watersheds. Our results demonstrate potential of the Bayesian estimation method as a rapid surveillance tool in monitoring and tracking changes of hydrological responses for small watersheds.


title: Productive Multitasking for Industrial Robots

abstract: The application of robotic solutions to small-batch production is challenging: economical constraints tend to dramatically limit the time for setting up new batches. Organizing robot tasks into modular software components, called skills, and allowing the assignment of multiple concurrent tasks to a single robot is potentially game-changing. However, due to cycle time constraints, it may be necessary for a skill to take over without waiting on another to terminate, and the available literature lacks a systematic approach in this case. In the present article, we fill the gap by (a) establishing the specifications of skills that can be sequenced with partial executions, (b) proposing an implementation based on the combination of finite-state machines and behavior trees, and (c) demonstrating the benefits of such skills through extensive trials in the environment of ARIAC (Agile Robotics for Industrial Automation Competition).


title: Bandwidth Allocation and Service Differentiation in D2D Wireless Networks

abstract: Inspired by a new feature in 5G NR called bandwidth part (BWP), this paper presents a bandwidth allocation (BA) model that allows one to adapt the bandwidth allocated to users depending on their data rate needs. Specifically, in adaptive BA, a wide bandwidth is divided into chunks of smaller bandwidths and the number of bandwidth chunks allocated to a user depends on its needs or type. Although BWP in 5G NR mandates allocation of a set of contiguous bandwidth chunks, our BA model also allows other assumptions on chunk allocation such as the allocation of any set of bandwidth chunks, as in, e.g., LTE resource allocation, where chunks are selected uniformly at random. The BA model studied here is probabilistic in that the user locations are assumed to form a realization of a Poisson point process and each user decides independently to be of a certain type with some probability. This model allows one to quantify spectrum sharing and service differentiation in this context, namely to predict what performance a user gets depending on its type as well as the overall performance. This is based on exact representations of key performance metrics for each user type, namely its success probability, the meta distribution of its signal-to-interference ratio, and its Shannon throughput. We show that, surprisingly, the higher traffic variability stemming from adaptive BA is beneficial: when comparing two networks using adaptive BA and having the same mean signal and the same mean interference powers, the network with higher traffic variability performs better for all these performance metrics. With respect to Shannon throughput, we observe that our BA model is roughly egalitarian per Hertz and leads to a linear service differentiation in aggregated throughput value.


title: Interpolating the Stochastic Heat and Wave Equations with Time-independent Noise: Solvability and Exact Asymptotics

abstract: In this article, we study a class of stochastic partial differential equations with fractional differential operators subject to some time-independent multiplicative Gaussian noise. We derive sharp conditions, under which a unique global $L^p(\Omega)$-solution exits for all $p\ge 2$. In this case, we derive exact moment asymptotics following the same strategy in a recent work by Balan et al [1]. In the case when there exits only a local solution, we determine the precise deterministic time, $T_2$, before which a unique $L^2(\Omega)$-solution exits, but after which the series corresponding to the $L^2(\Omega)$ moment of the solution blows up. By properly choosing the parameters, results in this paper interpolate the known results for both stochastic heat and wave equations.


title: Graded Galois Lattices and Closed Itemsets

abstract: The Galois lattice is a graphic method of representing knowledge structures. The first basic purpose in this paper is to introduce a new class of Galois lattices, called graded Galois lattices. As a direct result, one can obtain the notion of graded closed itemsets (sets of items), to extend the definition of closed itemsets. Our second important goal in this paper, is related to set a constructive method, computing the graded formal concepts and graded closed itemsets. We mean by a constructive method, a method that builds up a complete solution from scratch by sequentially adding components to a partial solution until the solution is complete. Besides of computational aspects, our methods in this paper are based on the strong results obtained by special mappings in the realm of domain theory. To reach the fertilized consequences and constructive algorithms, we need to push the study to the structures of Banach lattices.


title: Faster Exponential Algorithm for Permutation Pattern Matching

abstract: The Permutation Pattern Matching problem asks, given two permutations $\sigma$ on $n$ elements and $\pi$, whether $\sigma$ admits a subsequence with the same relative order as $\pi$ (or, in the counting version, how many such subsequences are there). This natural problem was shown by Bose, Buss and Lubiw [IPL 1998] to be NP-complete, and hence we should seek exact exponential time solutions. The asymptotically fastest such solution up to date, by Berendsohn, Kozma and Marx [IPEC 2019], works in $\mathcal{O}(1.6181^n)$ time. We design a simple and faster $\mathcal{O}(1.415^{n})$ time algorithm for both the detection and the counting version. We also prove that this is optimal among a certain natural class of algorithms.


title: Enabling Longitudinal Exploratory Analysis of Clinical COVID Data

abstract: As the COVID-19 pandemic continues to impact the world, data is being gathered and analyzed to better understand the disease. Recognizing the potential for visual analytics technologies to support exploratory analysis and hypothesis generation from longitudinal clinical data, a team of collaborators worked to apply existing event sequence visual analytics technologies to a longitudinal clinical data from a cohort of 998 patients with high rates of COVID-19 infection. This paper describes the initial steps toward this goal, including: (1) the data transformation and processing work required to prepare the data for visual analysis, (2) initial findings and observations, and (3) qualitative feedback and lessons learned which highlight key features as well as limitations to address in future work.


title: Inverse Sampling of Degenerate Datasets from a Linear Regression Line

abstract: When linear regression generates a relationship between a (dependent) scalar response and one or multiple independent variables, various datasets providing distinct graphical trends can develop resembling relationships based on the same statistical properties. Advanced statistical approaches, such as neural networks and machine learning methods, are of great necessity to process, characterize, and analyze these degenerate datasets. On the other hand, the accurate creation of purposedly degenerate datasets is essential to test new models in the research and education of applied statistics. In this light, the present study characterizes the famous Anscombe datasets and provides a general algorithm for creating multiple paired datasets of identical statistical properties.


title: Quantum Metric and Correlated States in Two-dimensional Systems

abstract: The recent realization of twisted, two-dimensional, bilayers exhibiting strongly correlated states has created a platform in which the relation between the properties of the electronic bands and the nature of the correlated states can be studied in unprecedented ways. The reason is that these systems allow extraordinary control of the electronic bands' properties, for example by varying the relative twist angle between the layers forming the system. In particular, in twisted bilayers the low energy bands can be tuned to be very flat and with a nontrivial quantum metric. This allows the quantitative and experimental exploration of the relation between the metric of Bloch quantum states and the properties of correlated states. In this work we first review the general connection between quantum metric and the properties of correlated states that break a continuous symmetry. We then discuss the specific case when the correlated state is a superfluid and show how the quantum metric is related to the superfluid stiffness. To exemplify such relation we show results for the case of superconductivity in magic angle twisted bilayer graphene. We conclude by discussing possible research directions to further elucidate the connection between quantum metric and correlated states' properties.


title: Quantendynamik offener Systeme, Trendbericht Theoretische Chemie 2020, Teil 3

abstract: Quantum mechanics describes the unitary time evolution of isolated systems. In reality, every quantum system interacts with its environment, leading to an irreversible loss of the phase relation. Path integral based methods provide a rigorous description of the process of decoherence. We provide an overview over recent numerical methods suited to describe the non-perturbative regime of system bath interaction. -- Quantenmechanik beschreibt die unit\"are Zeitentwicklung isolierter Systeme. In der Realit\"at interagiert jedes Quantensystem mit seiner Umwelt, was zu einem irreversiblen Verlust von Phasenbeziehungen f\"uhrt. Pfadintegralbasierte Methoden erlauben eine rigorose Beschreibung dieser Dekoh\"arenzprozesse.


title: On Approximate Nearest Neighbour Selection for Multi-Stage Dense Retrieval

abstract: Dense retrieval, which describes the use of contextualised language models such as BERT to identify documents from a collection by leveraging approximate nearest neighbour (ANN) techniques, has been increasing in popularity. Two families of approaches have emerged, depending on whether documents and queries are represented by single or multiple embeddings. ColBERT, the exemplar of the latter, uses an ANN index and approximate scores to identify a set of candidate documents for each query embedding, which are then re-ranked using accurate document representations. In this manner, a large number of documents can be retrieved for each query, hindering the efficiency of the approach. In this work, we investigate the use of ANN scores for ranking the candidate documents, in order to decrease the number of candidate documents being fully scored. Experiments conducted on the MSMARCO passage ranking corpus demonstrate that, by cutting of the candidate set by using the approximate scores to only 200 documents, we can still obtain an effective ranking without statistically significant differences in effectiveness, and resulting in a 2x speedup in efficiency.


title: Learning to discover: expressive Gaussian mixture models for multi-dimensional simulation and parameter inference in the physical sciences

abstract: We show that density models describing multiple observables with (i) hard boundaries and (ii) dependence on external parameters may be created using an auto-regressive Gaussian mixture model. The model is designed to capture how observable spectra are deformed by hypothesis variations, and is made more expressive by projecting data onto a configurable latent space. It may be used as a statistical model for scientific discovery in interpreting experimental observations, for example when constraining the parameters of a physical model or tuning simulation parameters according to calibration data. The model may also be sampled for use within a Monte Carlo simulation chain, or used to estimate likelihood ratios for event classification. The method is demonstrated on simulated high-energy particle physics data considering the anomalous electroweak production of a $Z$ boson in association with a dijet system at the Large Hadron Collider, and the accuracy of inference is tested using a realistic toy example. The developed methods are domain agnostic; they may be used within any field to perform simulation or inference where a dataset consisting of many real-valued observables has conditional dependence on external parameters.


title: ETA Prediction with Graph Neural Networks in Google Maps

abstract: Travel-time prediction constitutes a task of high importance in transportation networks, with web mapping services like Google Maps regularly serving vast quantities of travel time queries from users and enterprises alike. Further, such a task requires accounting for complex spatiotemporal interactions (modelling both the topological properties of the road network and anticipating events -- such as rush hours -- that may occur in the future). Hence, it is an ideal target for graph representation learning at scale. Here we present a graph neural network estimator for estimated time of arrival (ETA) which we have deployed in production at Google Maps. While our main architecture consists of standard GNN building blocks, we further detail the usage of training schedule methods such as MetaGradients in order to make our model robust and production-ready. We also provide prescriptive studies: ablating on various architectural decisions and training regimes, and qualitative analyses on real-world situations where our model provides a competitive edge. Our GNN proved powerful when deployed, significantly reducing negative ETA outcomes in several regions compared to the previous production baseline (40+% in cities like Sydney).


title: Heavy-tailed Streaming Statistical Estimation

abstract: We consider the task of heavy-tailed statistical estimation given streaming $p$-dimensional samples. This could also be viewed as stochastic optimization under heavy-tailed distributions, with an additional $O(p)$ space complexity constraint. We design a clipped stochastic gradient descent algorithm and provide an improved analysis, under a more nuanced condition on the noise of the stochastic gradients, which we show is critical when analyzing stochastic optimization problems arising from general statistical estimation problems. Our results guarantee convergence not just in expectation but with exponential concentration, and moreover does so using $O(1)$ batch size. We provide consequences of our results for mean estimation and linear regression. Finally, we provide empirical corroboration of our results and algorithms via synthetic experiments for mean estimation and linear regression.


title: Phonon induced rank-2 U(1) nematic liquid states

abstract: Fascinating new phases of matter can emerge from strong electron interactions in solids. In recent years, a new exotic class of many-body phases, described by generalized electromagnetism of symmetric rank-2 electric and magnetic fields and immobile charge excitations dubbed \textit{fractons}, has attracted wide attention. Besides having interesting properties in their own right, the models with generalized electromagnetism are also closely related to gapped fracton quantum orders, new phases of dipole-covering systems, as well as quantum information and quantum gravity. However, experimental realization of the rank-2 U(1) gauge theory is still absent and even known practical experimental routes are scarce. In this work, we propose a scheme of coupled optical phonons and nematic degrees of freedom, as well as several concrete experimental platforms for their realizations. We show that these systems can realize the electrostatics sector of the rank-2 U(1) gauge theory. A great advantage of the proposed scheme is that it requires only the basic ingredients of phonon and nematic physics, and hence may be applicable to a wide range of experimental realizations from liquid crystals to electron orbitals.


title: Chung-type law of the iterated logarithm and exact moduli of continuity for a class of anisotropic Gaussian random fields

abstract: We establish a Chung-type law of the iterated logarithm and the exact local and uniform moduli of continuity for a large class of anisotropic Gaussian random fields with a harmonizable-type integral representation and the property of strong local nondeterminism. Compared with the existing results in the literature, our results do not require the assumption of stationary increments and provide more precise upper and lower bounds for the limiting constants. The results are applicable to the solutions of a class of linear stochastic partial differential equations driven by a fractional-colored Gaussian noise, including the stochastic heat equation.


title: "Look! It's a Computer Program! It's an Algorithm! It's AI!": Does Terminology Affect Human Perceptions and Evaluations of Algorithmic Decision-Making Systems?

abstract: In the media, in policy-making, but also in research articles, algorithmic decision-making (ADM) systems are referred to as algorithms, artificial intelligence, and computer programs, amongst other terms. We hypothesize that such terminological differences can affect people's perceptions of properties of ADM systems, people's evaluations of systems in application contexts, and the replicability of research as findings may be influenced by terminological differences. In two studies (N = 397, N = 622), we show that terminology does indeed affect laypeople's perceptions of system properties (e.g., perceived complexity) and evaluations of systems (e.g., trust). Our findings highlight the need to be mindful when choosing terms to describe ADM systems, because terminology can have unintended consequences, and may impact the robustness and replicability of HCI research. Additionally, our findings indicate that terminology can be used strategically (e.g., in communication about ADM systems) to influence people's perceptions and evaluations of these systems.


title: Phase-space matrix representation of differential equations for obtaining the energy spectrum of model quantum systems

abstract: Employing the phase-space representation of second order ordinary differential equations we developed a method to find the eigenvalues and eigenfunctions of the 1-dimensional time independent Schr\"odinger equation for quantum model systems. The method presented simplifies some approaches shown in textbooks, based on asymptotic analyses of the time-independent Schr\"odinger equation, and power series methods with recurrence relations. In addition, the method presented here facilitates the understanding of the relationship between the ordinary differential equations of the mathematical physics and the time independent Schr\"odinger equation of physical models as the harmonic oscillator, the rigid rotor, the Hydrogen atom, and the Morse oscillator.


title: Excess and Deficiency of Extreme Multidimensional Random Fields

abstract: Probability distributions and densities are derived for the excess and deficiency of the intensity or instantaneous energy (quasi-static power) associated with a $p$-dimensional random vector field. Explicit expressions for the exact distributions are obtained for arbitrary threshold levels, together with simple approximate functions for relatively high or low thresholds. It is shown that precise expressions only require an expansion of order $p-1$ in the ratio of the excess height to the threshold level. Numerical simulations validate the analytical results.

categories: physics.data-an hep-ex hep-th stat.AP stat.ME


title: The Interplay Between Implicit Bias and Benign Overfitting in Two-Layer Linear Networks

abstract: The recent success of neural network models has shone light on a rather surprising statistical phenomenon: statistical models that perfectly fit noisy data can generalize well to unseen test data. Understanding this phenomenon of $\textit{benign overfitting}$ has attracted intense theoretical and empirical study. In this paper, we consider interpolating two-layer linear neural networks trained with gradient flow on the squared loss and derive bounds on the excess risk when the covariates satisfy sub-Gaussianity and anti-concentration properties, and the noise is independent and sub-Gaussian. By leveraging recent results that characterize the implicit bias of this estimator, our bounds emphasize the role of both the quality of the initialization as well as the properties of the data covariance matrix in achieving low excess risk.


title: Lie groups of Poisson diffeomorphisms

abstract: By considering suitable Poisson groupoids, we develop an approach to obtain Lie group structures on (subgroups of) the Poisson diffeomorphism groups of various classes of Poisson manifolds. As applications, we show that the Poisson diffeomorphism groups of (normal-crossing) log-symplectic, elliptic symplectic, scattering-symplectic and cosymplectic manifolds are regular infinite-dimensional Lie groups.


title: Linearization of Poisson groupoids

abstract: Motivated by a search for Lie group structures on groups of Poisson diffeomorphisms [24], we investigate linearizability of Poisson structures of Poisson groupoids around the unit section. After extending the Lagrangian neighbourhood theorem to the setting of cosymplectic Lie algebroids, we establish that dual integrations of triangular bialgebroids are always linearizable. Additionally, we show that the (non-dual) integration of a triangular Lie bialgebroid is linearizable whenever the $r$-matrix is of so-called cosymplectic type. The proof relies on the integration of a triangular Lie bialgebroid to a symplectic LA-groupoid, and in the process we define interesting new examples of double Lie algebroids and LA-groupoids. We also show that the product Poisson groupoid can only be linearizable when the Poisson structure on the unit space is regular.


title: Influence of Solvent Quality on Depletion Potentials in Colloid-Polymer Mixtures

abstract: As first explained by the classic Asakura-Oosawa (AO) model, effective attractive forces between colloidal particles induced by depletion of nonadsorbing polymers can drive demixing of colloid-polymer mixtures into colloid-rich and colloid-poor phases, with practical relevance for purification of water, stability of foods and pharmaceuticals, and macromolecular crowding in biological cells. By idealizing polymer coils as effective penetrable spheres, the AO model qualitatively captures the influence of polymer depletion on thermodynamic phase behavior of colloidal suspensions. In previous work, we extended the AO model to incorporate aspherical polymer conformations and showed that fluctuating shapes of random-walk coils can significantly modify depletion potentials [W. K. Lim and A. R. Denton, Soft Matter 12, 2247 (2016); J. Chem. Phys. 144, 024904 (2016)]. We further demonstrated that the shapes of polymers in crowded environments depend sensitively on solvent quality [W. J. Davis and A. R. Denton, J. Chem. Phys. 149, 124901 (2018)]. Here we apply Monte Carlo simulation to analyze the influence of solvent quality on depletion potentials in mixtures of hard sphere colloids and nonadsorbing polymer coils, modeled as ellipsoids whose principal radii fluctuate according to random-walk statistics. We consider both self-avoiding and non-self-avoiding random walks, corresponding to polymers in good and theta solvents, respectively. Our simulation results demonstrate that depletion of polymers of equal molecular weight induces much stronger attraction between colloids in good solvents than in theta solvents and confirm that depletion interactions are significantly influenced by aspherical polymer conformations.


title: Valence quark ratio in the proton

abstract: Beginning with precise data on the ratio of structure functions in deep inelastic scattering (DIS) from $^3$He and $^3$H, collected on the domain $0.19 \leq x_B \leq 0.83$, where $x_B$ is the Bjorken scaling variable, we employ a robust method for extrapolating such data to arrive at a model-independent result for the $x_B=1$ value of the ratio of neutron and proton structure functions. Combining this with information obtained in analyses of DIS from nuclei, corrected for target-structure dependence, we arrive at a prediction for the proton's valence-quark ratio: $\left. d_v/u_v \right|_{x_B\to 1} = 0.230 (57)$. Requiring consistency with this result presents a challenge to many descriptions of proton structure.


title: Super-resolution of spin configurations based on flow-based generative models

abstract: We present a super-resolution method for spin systems using a flow-based generative model that is a deep generative model with reversible neural network architecture. Starting from spin configurations on a two-dimensional square lattice, our model generates spin configurations of a larger lattice. As a flow-based generative model precisely estimates the distribution of the generated configurations, it can be combined with Monte Carlo simulation to generate large lattice configurations according to the Boltzmann distribution. Hence, the long-range correlation on a large configuration is reduced into the shorter one through the flow-based generative model. This alleviates the critical slowing down near the critical temperature. We demonstrated 8 times increased lattice size in the linear dimensions using our super-resolution scheme repeatedly. We numerically show that by performing simulations for $16\times 16$ configurations, our model can sample lattice configurations at $128\times 128$ on which the thermal average of physical quantities has good agreement with the one evaluated by the traditional Metropolis-Hasting Monte Carlo simulation.


title: Binary Signed-Digit Integers and the Stern Diatomic Sequence

abstract: Stern's diatomic sequence is a well-studied and simply defined sequence with many fascinating characteristics. The binary signed-digit representation of integers is an alternative representation of integers with much use in efficient computation, coding theory and cryptography. We link these two ideas here, showing that the number of $i$-bit binary signed-digit representations of an integer $n$ with $n<2^i$ is the $(2^i-n)^\text{th}$ element in Stern's diatomic sequence. This correspondence makes the vast range of results known for Stern's diatomic sequence available for consideration in the study of binary signed-digit integers.


title: A Class of Trees Having Near-Best Balance

abstract: Full binary trees naturally represent commutative non-associative products. There are many important examples of these products: finite-precision floating-point addition and NAND gates, among others. Balance in such a tree is highly desirable for efficiency in calculation. The best balance is attained with a divide-and-conquer approach. However, this may not be the optimal solution, since the success of many calculations is dependent on the grouping and ordering of the calculation, for reasons ranging from the avoidance of rounding error, to calculating with varying precision, to the placement of calculation within a heterogeneous system. We introduce a new class of computational trees having near-best balance in terms of the Colless index from mathematical phylogenetics. These trees are easily constructed from the binary decomposition of the number of terms in the problem. They also permit much more flexibility than the optimally balanced divide-and-conquer trees. This gives needed freedom in the grouping and ordering of calculation, and allows intelligent efficiency trade-offs.


title: Computing the viscous effect in early-time drop impact dynamics

abstract: The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. (2014b) show that the liquid-air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier-Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square root scaling relationship. We also predict the effect of interfacial tension at the liquid-gas interface on the drop impact, showing that it mediates the lift-off behavior.


title: Physics-informed neural networks for improving cerebral hemodynamics predictions

abstract: Determining brain hemodynamics plays a critical role in the diagnosis and treatment of various cerebrovascular diseases. In this work, we put forth a physics-informed deep learning framework that augments sparse clinical measurements with fast computational fluid dynamics (CFD) simulations to generate physically consistent and high spatiotemporal resolution of brain hemodynamic parameters. Transcranial Doppler (TCD) ultrasound is one of the most common techniques in the current clinical workflow that enables noninvasive and instantaneous evaluation of blood flow velocity within the cerebral arteries. However, it is spatially limited to only a handful of locations across the cerebrovasculature due to the constrained accessibility through the skull's acoustic windows. Our deep learning framework employs in-vivo real-time TCD velocity measurements at several locations in the brain and the baseline vessel cross-sectional areas acquired from 3D angiography images, and provides high-resolution maps of velocity, area, and pressure in the entire vasculature. We validated the predictions of our model against in-vivo velocity measurements obtained via 4D flow MRI scans. We then showcased the clinical significance of this technique in diagnosing the cerebral vasospasm (CVS) by successfully predicting the changes in vasospastic local vessel diameters based on corresponding sparse velocities measurements. The key finding here is that the combined effects of uncertainties in outlet boundary condition subscription and modeling physics deficiencies render the conventional purely physics-based computational models unsuccessful in recovering accurate brain hemodynamics. Nonetheless, fusing these models with clinical measurements through a data-driven approach ameliorates predictions of brain hemodynamic variables.

categories: physics.med-ph cs.LG physics.bio-ph physics.flu-dyn


title: Model Predictive Control approach to adaptive messaging intervention for physical activity

abstract: In this work, we have developed a framework for synthesizing data driven controllers for a class of uncertain switched systems arising in an application to physical activity interventions. In particular, we present an application of probabilistic model predictive control to design an efficient, tractable, and adaptive intervention using behavioral data sets i.e. physical activity behavior. The models of physical activity for each individual are provided for the design of controllers that maximize the probability of achieving a desired physical activity goal subject to intervention specifications. We have tailored the mixed-integer programming-based approach for evaluating the Model Predictive Control decision at each time step.


title: The use of the Born-Oppenheimer factorization in the phase-space representation of the time-independent Schr\"odinger equation for bilinearly coupled harmonic oscillators

abstract: A system of two bilinearly coupled harmonic oscillators has been solved analytically by using the Born-Oppenheimer (BO) product wavefunction ansatz and the phase-space bound trajectory approach [J. S. Molano et al., Chem. Phys. Lett. \textbf{76}(12), 138171 (2021)]. The bilinearly coupled oscillator system allows to obtain the analytical expression of the quantum system, facilitating comparison with the results of using the BO ansatz product. The analytical and BO wavefunctions are obtained as a product of parabolic cylinder functions. The arguments of the parabolic cylinder functions of the exact and the BO wavefunctions are related to the physical coordinates by linear transformations, represented by matrices $\mathsf{G}$ and $\mathsf{\tilde{G}}$ respectively. A $QU$ decomposition of the matrix $\mathsf{G}$ outputs an upper triangular matrix $\mathsf{U}$ that is closely related to the BO $\mathsf{\tilde{G}}$ matrix. The effect of the BO non-adiabatic coupling is analyzed on the stability of the phase-space equations of motion. The eigenvalues and eigenfunctions obtained by the Born-Oppenheimer phase-space trajectory approach show excellent agreement with the analytical solutions.


title: Improving Object Detection and Attribute Recognition by Feature Entanglement Reduction

abstract: We explore object detection with two attributes: color and material. The task aims to simultaneously detect objects and infer their color and material. A straight-forward approach is to add attribute heads at the very end of a usual object detection pipeline. However, we observe that the two goals are in conflict: Object detection should be attribute-independent and attributes be largely object-independent. Features computed by a standard detection network entangle the category and attribute features; we disentangle them by the use of a two-stream model where the category and attribute features are computed independently but the classification heads share Regions of Interest (RoIs). Compared with a traditional single-stream model, our model shows significant improvements over VG-20, a subset of Visual Genome, on both supervised and attribute transfer tasks.


title: Probing the Structure of Deuteron at Very Short Distances

abstract: We study the electrodisintegration of deuteron at quasi-elastic kinematics and high transferred momentum as a probe for the short distance structure in nuclei. In this reaction, an electron hits a nucleus of deuterium, which breaks up into a proton-neutron pair. We focus our attention on events where fast nucleons emerge, corresponding to nuclear configurations where the bound nucleons have a high relative momentum (exceeding 700 MeV/c). The present research is relevant to physical systems where high-density nuclear matter is present. This condition covers a wide range of physics, from neutron stars to nuclei stability and the repulsive nuclear core. The present work differs from previous studies in two crucial features. One is the definition of the deuteron wave function, which include terms of a relativistic origin that can be ordered based on their relative contribution to the deuteron's internal momentum distribution. These terms, related to the off-shell properties of the nucleon-nucleon bound-state, do not occur in non-relativistic quantum mechanics. However, they become increasingly important in describing configurations with a high nucleon-nucleon relative momentum. The second difference is that we account for the off-shell nature of the bound nucleon that enters on the definition of the (half-off-shell) electromagnetic current. We avoid many of the difficulties inherent to the relativistic nature of the processes involved by adopting a theoretical framework known as Light Front dynamics. Simultaneous simplifications in the definition of the relativistic wave function for the proton-neutron bound state and the treatment of the (half-off-shell) electromagnetic current for the bound nucleon are among the essential advantages resulting from the use of the Light Front dynamics.


title: Maneuver Identification Challenge

abstract: AI algorithms that identify maneuvers from trajectory data could play an important role in improving flight safety and pilot training. AI challenges allow diverse teams to work together to solve hard problems and are an effective tool for developing AI solutions. AI challenges are also a key driver of AI computational requirements. The Maneuver Identification Challenge hosted at maneuver-id.mit.edu provides thousands of trajectories collected from pilots practicing in flight simulators, descriptions of maneuvers, and examples of these maneuvers performed by experienced pilots. Each trajectory consists of positions, velocities, and aircraft orientations normalized to a common coordinate system. Construction of the data set required significant data architecture to transform flight simulator logs into AI ready data, which included using a supercomputer for deduplication and data conditioning. There are three proposed challenges. The first challenge is separating physically plausible (good) trajectories from unfeasible (bad) trajectories. Human labeled good and bad trajectories are provided to aid in this task. Subsequent challenges are to label trajectories with their intended maneuvers and to assess the quality of those maneuvers.


title: Studying strangeness and baryon production in small systems through $\Xi-$hadron correlations using the ALICE detector

abstract: These proceedings summarise recent measurements of angular correlations between the $\Xi$ baryon and identified hadrons in pp collisions at $\sqrt{s}=13\,\mathrm{TeV}$ using the ALICE detector. The results are compared with both string-based (PYTHIA8 with extensions) and core-corona (EPOS-LHC) models, to improve our understanding of strangeness and baryon production in small systems. The results favour baryon production through string junctions over diquark breaking, but the PYTHIA models fail at describing the relatively wide $\Xi-$strangeness jet peak, indicating stronger diffusion of strange quarks in data. On the other hand, EPOS-LHC is missing local conservation of quantum numbers, making it difficult to draw any conclusion about the core-corona model.


title: Generalized Real-World Super-Resolution through Adversarial Robustness

abstract: Real-world Super-Resolution (SR) has been traditionally tackled by first learning a specific degradation model that resembles the noise and corruption artifacts in low-resolution imagery. Thus, current methods lack generalization and lose their accuracy when tested on unseen types of corruption. In contrast to the traditional proposal, we present Robust Super-Resolution (RSR), a method that leverages the generalization capability of adversarial attacks to tackle real-world SR. Our novel framework poses a paradigm shift in the development of real-world SR methods. Instead of learning a dataset-specific degradation, we employ adversarial attacks to create difficult examples that target the model's weaknesses. Afterward, we use these adversarial examples during training to improve our model's capacity to process noisy inputs. We perform extensive experimentation on synthetic and real-world images and empirically demonstrate that our RSR method generalizes well across datasets without re-training for specific noise priors. By using a single robust model, we outperform state-of-the-art specialized methods on real-world benchmarks.


title: High energy spin excitations in the quantum spin liquid candidate Zn-barlowite probed by resonant inelastic x-ray scattering

abstract: A quantum spin liquid is a novel ground state that can support long-range entanglement between magnetic moments, resulting in exotic spin excitations involving fractionalized $S=\frac{1}{2}$ spinons. Here, we measure the excitations in single crystals of the spin liquid candidate Zn-barlowite using resonant inelastic X-ray scattering. By analyzing the incident polarization and temperature dependences, we deduce a clear magnetic scattering contribution forming a broad continuum that surprisingly extends up to $\sim$200 meV ($\sim$14$J$, where $J$ is the magnetic exchange). The excitation spectrum reveals that significant contributions arise from multiple pairs of spinons and/or antispinons at high energies.


title: Hardware-assisted Trusted Memory Disaggregation for Secure Far Memory

abstract: Memory disaggregation provides efficient memory utilization across network-connected systems. It allows a node to use part of memory in remote nodes in the same cluster. Recent studies have improved RDMA-based memory disaggregation systems, supporting lower latency and higher bandwidth than the prior generation of disaggregated memory. However, the current disaggregated memory systems manage remote memory only at coarse granularity due to the limitation of the access validation mechanism of RDMA. In such systems, to support fine-grained remote page allocation, the trustworthiness of all participating systems needs to be assumed, and thus a security breach in a node can propagate to the entire cluster. From the security perspective, the memory-providing node must protect its memory from memory-requesting nodes. On the other hand, the memory-requesting node requires the confidentiality and integrity protection of its memory contents even if they are stored in remote nodes. To address the weak isolation support in the current system, this study proposes a novel hardware-assisted memory disaggregation system. Based on the security features of FPGA, the logic in each per-node FPGA board provides a secure memory disaggregation engine. With its own networks, a set of FPGA-based engines form a trusted memory disaggregation system, which is isolated from the privileged software of each participating node. The secure memory disaggregation system allows fine-grained memory management in memory-providing nodes, while the access validation is guaranteed with the hardware-hardened mechanism. In addition, the proposed system hides the memory access patterns observable from remote nodes, supporting obliviousness. Our evaluation with FPGA implementation shows that such fine-grained secure disaggregated memory is feasible with comparable performance to the latest software-based techniques.


title: First Passage with Restart in Discrete Time: with applications to biased random walks on the half-line

abstract: In recent years, it has been well-established that adding a restart mechanism can alter the firstpassage statistics of a stochastic processes in useful and interesting ways. Though different mecha-nisms have been investigated, we derive a probability generating function for a discrete-time FirstPassage process Under Restart and use it to examine two examples, including a biased random walkon the non-negative integers.


title: Trade-off between deep learning for species identification and inference about predator-prey co-occurrence: Reproducible R workflow integrating models in computer vision and ecological statistics

abstract: Deep learning is used in computer vision problems with important applications in several scientific fields. In ecology for example, there is a growing interest in deep learning for automatizing repetitive analyses on large amounts of images, such as animal species identification. However, there are challenging issues toward the wide adoption of deep learning by the community of ecologists. First, there is a programming barrier as most algorithms are written in Python while most ecologists are versed in R. Second, recent applications of deep learning in ecology have focused on computational aspects and simple tasks without addressing the underlying ecological questions or carrying out the statistical data analysis to answer these questions. Here, we showcase a reproducible R workflow integrating both deep learning and statistical models using predator-prey relationships as a case study. We illustrate deep learning for the identification of animal species on images collected with camera traps, and quantify spatial co-occurrence using multispecies occupancy models. Despite average model classification performances, ecological inference was similar whether we analysed the ground truth dataset or the classified dataset. This result calls for further work on the trade-offs between time and resources allocated to train models with deep learning and our ability to properly address key ecological questions with biodiversity monitoring. We hope that our reproducible workflow will be useful to ecologists and applied statisticians. All material (source of the Rmarkdown notebook and auxiliary files) is available from https://github.com/oliviergimenez/computo-deeplearning-occupany-lynx.


title: Deep Reinforcement Learning in Computer Vision: A Comprehensive Survey

abstract: Deep reinforcement learning augments the reinforcement learning framework and utilizes the powerful representation of deep neural networks. Recent works have demonstrated the remarkable successes of deep reinforcement learning in various domains including finance, medicine, healthcare, video games, robotics, and computer vision. In this work, we provide a detailed review of recent and state-of-the-art research advances of deep reinforcement learning in computer vision. We start with comprehending the theories of deep learning, reinforcement learning, and deep reinforcement learning. We then propose a categorization of deep reinforcement learning methodologies and discuss their advantages and limitations. In particular, we divide deep reinforcement learning into seven main categories according to their applications in computer vision, i.e. (i)landmark localization (ii) object detection; (iii) object tracking; (iv) registration on both 2D image and 3D image volumetric data (v) image segmentation; (vi) videos analysis; and (vii) other applications. Each of these categories is further analyzed with reinforcement learning techniques, network design, and performance. Moreover, we provide a comprehensive analysis of the existing publicly available datasets and examine source code availability. Finally, we present some open issues and discuss future research directions on deep reinforcement learning in computer vision


title: Bayesian Estimation of the Hydroxyl Radical Diffusion Coefficient at Low Temperature and High Pressure from Atomistic Molecular Dynamics

abstract: The hydroxyl radical is the primary reactive oxygen species produced by the radiolysis of water, and is a significant source of radiation damage to living organisms. Mobility of the hydroxyl radical at low temperatures and/or high pressures is hence a potentially important factor in determining the challenges facing psychrophilic and/or barophilic organisms in high-radiation environments (e.g., ice-interface or undersea environments in which radiative heating is a potential heat and energy source). Here, we estimate the diffusion coefficient for the hydroxyl radical in aqueous solution, using a hierarchical Bayesian model based on atomistic molecular dynamics trajectories in TIP4P/2005 water over a range of temperatures and pressures.


title: On variational principles for polarization responses in electromechanical systems

abstract: Classical electrodynamics uses a dielectric constant to describe the polarization response of electromechanical systems to changes in an electric field. We generalize that description to include a wide variety of responses to changes in the electric field, as found in most systems and applications. Electromechanical systems can be found in many physical and biological applications, such as ion transport in membranes, batteries, and dielectric elastomers. We present a unified, thermodynamically consistent, variational framework for modeling electromechanical systems as they respond to changes in the electric field; that is to say, as they polarize. This framework is motivated and developed using the classical energetic variational approach (EnVarA). The coupling between the electric part and the chemo-mechanical parts of the system is described either by Lagrange multipliers or various energy relaxations. The classical polarization and its dielectrics and dielectric constants appear as outputs of this analysis. The Maxwell equations then become universal conservation laws of charge and current, conjoined to an electromechanical description of polarization. Polarization describes the entire electromechanical response to changes in the electric field and can sometimes be approximated as a dielectric constant or dielectric dispersion.


title: Learning Effective and Efficient Embedding via an Adaptively-Masked Twins-based Layer

abstract: Embedding learning for categorical features is crucial for the deep learning-based recommendation models (DLRMs). Each feature value is mapped to an embedding vector via an embedding learning process. Conventional methods configure a fixed and uniform embedding size to all feature values from the same feature field. However, such a configuration is not only sub-optimal for embedding learning but also memory costly. Existing methods that attempt to resolve these problems, either rule-based or neural architecture search (NAS)-based, need extensive efforts on the human design or network training. They are also not flexible in embedding size selection or in warm-start-based applications. In this paper, we propose a novel and effective embedding size selection scheme. Specifically, we design an Adaptively-Masked Twins-based Layer (AMTL) behind the standard embedding layer. AMTL generates a mask vector to mask the undesired dimensions for each embedding vector. The mask vector brings flexibility in selecting the dimensions and the proposed layer can be easily added to either untrained or trained DLRMs. Extensive experimental evaluations show that the proposed scheme outperforms competitive baselines on all the benchmark tasks, and is also memory-efficient, saving 60\% memory usage without compromising any performance metrics.


title: Bilateral Denoising Diffusion Models

abstract: Denoising diffusion probabilistic models (DDPMs) have emerged as competitive generative models yet brought challenges to efficient sampling. In this paper, we propose novel bilateral denoising diffusion models (BDDMs), which take significantly fewer steps to generate high-quality samples. From a bilateral modeling objective, BDDMs parameterize the forward and reverse processes with a score network and a scheduling network, respectively. We show that a new lower bound tighter than the standard evidence lower bound can be derived as a surrogate objective for training the two networks. In particular, BDDMs are efficient, simple-to-train, and capable of further improving any pre-trained DDPM by optimizing the inference noise schedules. Our experiments demonstrated that BDDMs can generate high-fidelity samples with as few as 3 sampling steps and produce comparable or even higher quality samples than DDPMs using 1000 steps with only 16 sampling steps (a 62x speedup).


title: Robust High-Resolution Video Matting with Temporal Guidance

abstract: We introduce a robust, real-time, high-resolution human video matting method that achieves new state-of-the-art performance. Our method is much lighter than previous approaches and can process 4K at 76 FPS and HD at 104 FPS on an Nvidia GTX 1080Ti GPU. Unlike most existing methods that perform video matting frame-by-frame as independent images, our method uses a recurrent architecture to exploit temporal information in videos and achieves significant improvements in temporal coherence and matting quality. Furthermore, we propose a novel training strategy that enforces our network on both matting and segmentation objectives. This significantly improves our model's robustness. Our method does not require any auxiliary inputs such as a trimap or a pre-captured background image, so it can be widely applied to existing human matting applications.


title: Ray tracing and finite element modeling of sound propagation in a compartment fire

abstract: A compartment fire (a fire in a room or building) creates temperature gradients and inhomogeneous time-varying temperature, density, and flow fields. This work compared experimental measurements of the room acoustic impulse/frequency response in a room with a fire to numerically modeled responses. The fire is modeled using Fire Dynamics Simulator (FDS). Acoustic modeling was performed using the temperature field computed by FDS. COMSOL Multiphysics was used for finite element acoustic modeling and Bellhop for ray-trace acoustics modeling. The results show that the fire causes wave-fronts to arrive earlier (due to the higher sound speed) and with more variation in the delay times (due to the sound speed perturbations). The frequency response shows that the modes are shifted up in frequency and high frequency (>2500 Hz) modes are significantly attenuated. Model results are compared with data and show good agreement in observed trends.


title: eROSITA calibration and performance verification phase: High-mass X-ray binaries in the Magellanic Clouds

abstract: During its performance verification phase, the soft X-ray instrument eROSITA aboard the Spektrum-Roentgen-Gamma(SRG) spacecraft observed large regions in the Magellanic Clouds, where almost 40 known high-mass X-ray binaries (HMXBs, including candidates) are located. We looked for new HMXBs in the eROSITA data, searched for pulsations in HMXB candidates and investigated the long-term behaviour of the full sample using archival X-ray and optical data. For sources sufficiently bright, a detailed spectral and temporal analysis of their eROSITA data was performed. A source detection analysis of the eROSITA images in different energy bands provided count rates and upper limits for the remaining sources. We report the discovery of a new Be/X-ray binary in the Large Magellanic Cloud. The transient SRGEt J052829.5-690345 was detected with a 0.2-8.0 keV luminosity of ~10^35 erg/s and exhibits a hard X-ray spectrum, typical for this class of HMXBs. The OGLE I-band light curve of the V~15.7 mag counterpart shows large variations by up to 0.75 mag, which occur quasi periodically with ~511 days. The eROSITA observations of the Small Magellanic Cloud covered 16 Be/X-ray binary pulsars, five of them were bright enough to accurately determine their current pulse period. The pulse periods for SXP 726 and SXP 1323 measured from eROSITA data are ~800 s and ~1006 s, respectively, far away from their discovery periods. Including archival XMM-Newton observations we update the spin-period history of the two long-period pulsars which show nearly linear trends in their period evolution since more than 15 years. The corresponding average spin-down rate for SXP 726 is 4.3 s/yr while SXP 1323 shows spin-up with a rate of -23.2 s/yr. We discuss the spin evolution of the two pulsars in the framework of quasi-spherical accretion.


title: A Systematic Study of Projection Biases in Weak Lensing Analysis

abstract: We present a systematic study of projection biases in the weak lensing analysis of the first year of data from the Dark Energy Survey (DES) experiment. In the analysis we used a $\Lambda$CDM model and three two-point correlation functions. We show that these biases are a consequence of projecting, or marginalizing, over parameters like $h$, $\Omega_b$, $n_s$ and $\Omega_\nu h^2$ that are both poorly constrained and correlated with the parameters of interest like $\Omega_m$, $\sigma_8$ and $S_8$. Covering the relevant parameter space we show that the projection biases are a function of where the true values of the poorly constrained parameters lie with respect to the parameter priors. For example, biases in the position of the posteriors can exceed the 1.5$\sigma$ level if the true values of $h$ and $n_s$ are close to the top of the prior's range and the true values of $\Omega_b$ and $\Omega_\nu h^2$ are close to the bottom of the range of their priors. We also show that in some cases the 1D credible intervals can be over-specified by as much as 30% and coverage can be as low as 27%. Finally we estimate these projection biases for the analysis of three and six years worth of DES data.


title: Towards merged-element transmons using silicon fins: the FinMET

abstract: A merged-element transmon (MET) device, based on silicon (Si) fins, is proposed and the first steps to form such a "FinMET" are demonstrated. This new application of fin technology capitalizes on the anisotropic etch of Si(111) relative to Si(110) to define atomically flat, high aspect ratio Si tunnel barriers with epitaxial superconductor contacts on the parallel side-wall surfaces. This process circumvents the challenges associated with the growth of low-loss insulating barriers on lattice matched superconductors. By implementing low-loss, intrinsic float-zone Si as the barrier material rather than commonly used, potentially lossy AlOx, the FinMET is expected to overcome problems with standard transmons by (1) reducing dielectric losses, (2) minimizing the formation of two-level system spectral features, (3) exhibiting greater control over barrier thickness and qubit frequency spread, especially when combined with commercial fin fabrication and atomic-layer digital etching; (4) potentially reducing the footprint by several orders of magnitude; and (5) allowing scalable fabrication. Here, as a first step to making such a device, the fabrication of Si fin capacitors on Si(110) substrates with shadow-deposited Al electrodes is demonstrated. These fin capacitors are then fabricated into lumped element resonator circuits and probed using low-temperature microwave measurements. Further thinning of silicon junctions towards the tunneling regime will enable the scalable fabrication of FinMET devices based on existing silicon technology, while simultaneously avoiding lossy amorphous dielectrics for the tunnel barriers.

categories: quant-ph cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph


title: BRST Quantization of General Relativity in Unimodular Gauge and Unimodular Gravity

abstract: Unimodular gravity (UG) is an important theory which may explain the smallness of the cosmological constant. To get insight into the covariant quantization of UG, we discuss the BRST quantization of General Relativity (GR) with a cosmological constant in the unimodular gauge. We develop a novel way to gauge fix the transverse diffeomorphism (TDiff) and then further to fulfill the unimodular gauge. This process requires the introduction of an additional pair of BRST doublets which decouple from the physical sector together with the other three pairs of BRST doublets for the TDiff. We show that the physical spectrum is the same as GR in the usual covariant gauge fixing. We then study a theory derived by making "Fourier transform" of GR in the unimodular gauge with respect to the cosmological constant as a candidate of "quantum UG." We clarify the difference from GR and point out problems in this theory.


title: Efficient On-Chip Communication for Parallel Graph-Analytics on Spatial Architectures

abstract: Large-scale graph processing has drawn great attention in recent years. Most of the modern-day datacenter workloads can be represented in the form of Graph Processing such as MapReduce etc. Consequently, a lot of designs for Domain-Specific Accelerators have been proposed for Graph Processing. Spatial Architectures have been promising in the execution of Graph Processing, where the graph is partitioned into several nodes and each node works in parallel. We conduct experiments to analyze the on-chip movement of data in graph processing on a Spatial Architecture. Based on the observations, we identify a data movement bottleneck, in the execution of such highly parallel processing accelerators. To mitigate the bottleneck we propose a novel power-law aware Graph Partitioning and Data Mapping scheme to reduce the communication latency by minimizing the hop counts on a scalable network-on-chip. The experimental results on popular graph algorithms show that our implementation makes the execution 2-5x faster and 2.7-4x energy-efficient by reducing the data movement time in comparison to a baseline implementation.


title: Inverse boundary value problems for polyharmonic operators with non-smooth coefficients

abstract: We consider inverse boundary value problems for polyharmonic operators and in particular, the problem of recovering the coefficients of terms up to order one. The main interest of our result is that it further relaxes the regularity required to establish uniqueness. The proof relies on an averaging technique introduced by Haberman and Tataru for the study of an inverse boundary value problem for a second order operator.


title: SOMTimeS: Self Organizing Maps for Time Series Clustering and its Application to Serious Illness Conversations

abstract: There is an increasing demand for scalable algorithms capable of clustering and analyzing large time series datasets. The Kohonen self-organizing map (SOM) is a type of unsupervised artificial neural network for visualizing and clustering complex data, reducing the dimensionality of data, and selecting influential features. Like all clustering methods, the SOM requires a measure of similarity between input data (in this work time series). Dynamic time warping (DTW) is one such measure, and a top performer given that it accommodates the distortions when aligning time series. Despite its use in clustering, DTW is limited in practice because it is quadratic in runtime complexity with the length of the time series data. To address this, we present a new DTW-based clustering method, called SOMTimeS (a Self-Organizing Map for TIME Series), that scales better and runs faster than other DTW-based clustering algorithms, and has similar performance accuracy. The computational performance of SOMTimeS stems from its ability to prune unnecessary DTW computations during the SOM's training phase. We also implemented a similar pruning strategy for K-means for comparison with one of the top performing clustering algorithms. We evaluated the pruning effectiveness, accuracy, execution time and scalability on 112 benchmark time series datasets from the University of California, Riverside classification archive. We showed that for similar accuracy, the speed-up achieved for SOMTimeS and K-means was 1.8x on average; however, rates varied between 1x and 18x depending on the dataset. SOMTimeS and K-means pruned 43% and 50% of the total DTW computations, respectively. We applied SOMtimeS to natural language conversation data collected as part of a large healthcare cohort study of patient-clinician serious illness conversations to demonstrate the algorithm's utility with complex, temporally sequenced phenomena.


title: Objective quantum fields, retrocausality and ontology

abstract: We compare different approaches to quantum ontology. In particular, we discuss an interpretation of quantum mechanics that we call objective quantum field theory (OQFT), which involves retrocausal fields. Here, objective implies the existence of fields independent of an observer, but not that the results of conjugate measurements are predetermined: the theory is contextual. The ideas and analyses of Einstein and Bohr through to more recent approaches to objective realism are discussed. We briefly describe measurement induced projections, the guided wave interpretation, many-universes, consistent histories, and modal theories. These earlier interpretations are compared with OQFT. We argue that this approach is compatible both with Bohr's quantum complementarity, and Einstein's objective realism.


title: Supercomputing Enabled Deployable Analytics for Disaster Response

abstract: First responders and other forward deployed essential workers can benefit from advanced analytics. Limited network access and software security requirements prevent the usage of standard cloud based microservice analytic platforms that are typically used in industry. One solution is to precompute a wide range of analytics as files that can be used with standard preinstalled software that does not require network access or additional software and can run on a wide range of legacy hardware. In response to the COVID-19 pandemic, this approach was tested for providing geo-spatial census data to allow quick analysis of demographic data for better responding to emergencies. These data were processed using the MIT SuperCloud to create several thousand Google Earth and Microsoft Excel files representative of many advanced analytics. The fast mapping of census data using Google Earth and Microsoft Excel has the potential to give emergency responders a powerful tool to improve emergency preparedness. Our approach displays relevant census data (total population, population under 15, population over 65, median age) per census block, sorted by county, through a Microsoft Excel spreadsheet (xlsx file) and Google Earth map (kml file). The spreadsheet interface includes features that allow users to convert between different longitude and latitude coordinate units. For the Google Earth files, a variety of absolute and relative colors maps of population density have been explored to provide an intuitive and meaningful interface. Using several hundred cores on the MIT SuperCloud, new analytics can be generated in a few minutes.


title: Gravitational Wave Distances in Horndeski Cosmology

abstract: Gravitational wave propagation encounters a spacetime friction from a running Planck mass in modified gravity, causing the luminosity distance to deviate from that in general relativity (or given by the photon luminosity distance to the source), thus making it a valuable cosmological probe. We present the exact expression for the cosmological distance deviation in Horndeski gravity including theories that have a $G_5$ term yet propagate at the speed of light. An especially simple result ensues for coupled Gauss-Bonnet gravity, which we use to show it does not give a viable cosmology. We also generalize such coupling, and review the important connection of gravitational wave cosmological distance deviations to growth of cosmic structure measured by redshift space distortions.


title: Subspace Differential Privacy

abstract: Many data applications have certain invariant constraints due to practical needs. Data curators who employ differential privacy need to respect such constraints on the sanitized data product as a primary utility requirement. Invariants challenge the formulation, implementation, and interpretation of privacy guarantees. We propose subspace differential privacy, to honestly characterize the dependence of the sanitized output on confidential aspects of the data. We discuss two design frameworks that convert well-known differentially private mechanisms, such as the Gaussian and the Laplace mechanisms, to subspace differentially private ones that respect the invariants specified by the curator. For linear queries, we discuss the design of near-optimal mechanisms that minimize the mean squared error. Subspace differentially private mechanisms rid the need for post-processing due to invariants, preserve transparency and statistical intelligibility of the output, and can be suitable for distributed implementation. We showcase the proposed mechanisms on the 2020 Census Disclosure Avoidance demonstration data, and a spatio-temporal dataset of mobile access point connections on a large university campus.


title: A note on the Sundaram--Stanley bijection (or, Viennot for up-down tableaux)

abstract: We give a direct proof of a result of Sundaram and Stanley: that the dimension of the space of invariant vectors in a $2k$-fold tensor product of the vector representation of $\mathfrak{sp}_{2n}$ equals the number of $(n+1)$-avoiding matchings of $2k$ points. This can be viewed as an extension of Schensted's theorem on longest decreasing subsequences. Our main tool is an extension of Viennot's geometric construction to the setting of up-down tableaux.


title: The microstructural dependence of ionic transport in bi-continuous nanoporous metal

abstract: Ionic transports in nanopores hold the key to unlocking the full potential of bi-continuous nanoporous (NP) metals as advanced electrodes in electrochemical devices. The precise control of the uniform NP metal structures also provides us a unique opportunity to understand how complex structures determine transports at nanoscales. For NP Au from the dealloying of a Ag-Au alloy, we can tune the pore size in the range of 13 nm to 2.4 microns and the porosity between 38% and 69% via isothermal coarsening. For NP Ag from the reduction-induced decomposition of AgCl, we can control additionally its structural hierarchy and pore orientation. We measure the effective ionic conductivities of 1 M NaClO4 through these NP metals as membranes, which range from 7% to 44% of that of a free solution, corresponding to calculated pore tortuosities between 2.7 and 1.3. The tortuosity of NP Au displays weak dependences on both the pore size and the porosity, consistent with the observed self-similarity in the coarsening, except for those of pores < 25 nm, which we consider deviating from the well-coarsened pore geometry. For NP Ag, the low tortuosity of the hierarchical structure can be explained with the Maxwell-Garnett equation and that of the oriented structure underlines the random orientation as the cause of slow transport in other NP metals. At last, we achieve high current densities of CO2 reduction with these two low-tortuosity NP Ags, demonstrating the significance of the structure-transport relationships for designing functional NP metals.


title: A New Interpolation Approach and Corresponding Instance-Based Learning

abstract: Starting from finding approximate value of a function, introduces the measure of approximation-degree between two numerical values, proposes the concepts of "strict approximation" and "strict approximation region", then, derives the corresponding one-dimensional interpolation methods and formulas, and then presents a calculation model called "sum-times-difference formula" for high-dimensional interpolation, thus develops a new interpolation approach, that is, ADB interpolation. ADB interpolation is applied to the interpolation of actual functions with satisfactory results. Viewed from principle and effect, the interpolation approach is of novel idea, and has the advantages of simple calculation, stable accuracy, facilitating parallel processing, very suiting for high-dimensional interpolation, and easy to be extended to the interpolation of vector valued functions. Applying the approach to instance-based learning, a new instance-based learning method, learning using ADB interpolation, is obtained. The learning method is of unique technique, which has also the advantages of definite mathematical basis, implicit distance weights, avoiding misclassification, high efficiency, and wide range of applications, as well as being interpretable, etc. In principle, this method is a kind of learning by analogy, which and the deep learning that belongs to inductive learning can complement each other, and for some problems, the two can even have an effect of "different approaches but equal results" in big data and cloud computing environment. Thus, the learning using ADB interpolation can also be regarded as a kind of "wide learning" that is dual to deep learning.


title: How well is angular momentum accretion modelled in semi-analytic galaxy formation models?

abstract: Gas cooling and accretion in haloes delivers mass and angular momentum onto galaxies. In this work, we investigate the accuracy of the modelling of this important process in several different semi-analytic (SA) galaxy formation models (GALFORM, L-GALAXIES and MORGANA) through comparisons with a hydrodynamical simulation performed with the moving-mesh code AREPO. Both SA models and the simulation were run without any feedback or metal enrichment, in order to focus on the cooling and accretion process. All of the SA models considered here assume that gas cools from a spherical halo. We found that the assumption that the gas conserves its angular momentum when moving from the virial radius, $r_{\rm vir}$, to the central region of the halo, $r\sim 0.1 r_{\rm vir}$, is approximately consistent with the results from our simulation, in which gas typically retains $70-80\%$ of its angular momentum during this process. We also found that, compared to the simulation, the MORGANA model tends to overestimate the mean specific angular momentum of cooled-down gas, the L-GALAXIES model also tends to overestimate this in low-redshift massive haloes, while the two older GALFORM models tend to underestimate the angular momentum. In general, the predictions of the new GALFORM cooling model developed by Hou et al. agree the best with the simulation.


title: Structures of secular resonances for inner test particles in hierarchical planetary systems

abstract: In this study, dynamics of secular resonances for inner test particles are investigated under the octupole-level approximation by taking non-perturbative approaches. In practice, webs of the major secular resonances are produced by identifying families of stable periodic orbits and the associated stable libration zones are obtained by analysing Poincar\'e surfaces of section. By taking different values of the factor $\epsilon$ ($\epsilon$ measures the contribution of octupole terms), the influences of the octupole-order terms upon the dynamical structures are evaluated. Under the condition of $\epsilon = 0$ (no octupole-order contribution), the dynamical model is totally integrable and there is only Kozai resonance arising in the phase space. When the factor $\epsilon$ is different from zero, the dynamical structure in the phase space becomes complicated due to varieties of secular resonances appearing. Numerical results further indicate that (a) distributions of libration centres and stable libration zones remain qualitatively similar with different values of $\epsilon$, (b) Kozai resonance disappears due to the chaotic motion in the low-eccentricity region, and (c) the chaotic area arising in the low-eccentricity region increases with the factor $\epsilon$. Secular resonances are the source of many important dynamical phenomena, such as chaos, orbit alignment and orbit flipping, and thus the results presented in this work could be useful to understand the secular dynamics for those high-eccentricity and/or high-inclination objects in hierarchical planetary systems.


title: Quantum Markov monogamy inequalities

abstract: Markovianity lies at the heart of communication problems. This in turn makes the information-theoretic characterization of Markov processes worthwhile. Data processing inequalities are ubiquitous in this sense, assigning necessary conditions for all Markov processes. We address here the problem of the information-theoretic analysis of constraints on Markov processes in the quantum regime. We show the existence of a novel class of quantum data processing inequalities called here quantum Markov monogamy inequalities. This new class of necessary conditions on quantum Markov processes is inspired by its counterpart for classical Markov processes, and thus providing a strong link between classical and quantum constraints on Markovianity. We go on to construct a family of multitime quantum Markov monogamy inequalities, based on the process tensor formalism and that exploits multitime correlations. We then show, by means of an explicit example, that the Markov monogamy inequalities can be stronger than the usual quantum data processing inequalities.


title: Lies My Teacher Told Me About Density Functional Theory: Seeing Through Them with the Hubbard Dimer

abstract: Most realistic calculations of moderately correlated materials begin with a ground-state density functional theory (DFT) calculation. While Kohn-Sham DFT is used in about 40,000 scientific papers each year, the fundamental underpinnings are not widely appreciated. In this chapter, we analyze the inherent characteristics of DFT in their simplest form, using the asymmetric Hubbard dimer as an illustrative model. We begin by working through the core tenets of DFT, explaining what the exact ground-state density functional yields and does not yield. Given the relative simplicity of the system, almost all properties of the exact exchange-correlation functional are readily visualized and plotted. Key concepts include the Kohn-Sham scheme, the behavior of the XC potential as correlations become very strong, the derivative discontinuity and the difference between KS gaps and true charge gaps, and how to extract optical excitations using time-dependent DFT. By the end of this text and accompanying exercises, the reader will improve their ability to both explain and visualize the concepts of DFT, as well as better understand where others may go wrong.


title: ChessMix: Spatial Context Data Augmentation for Remote Sensing Semantic Segmentation

abstract: Labeling semantic segmentation datasets is a costly and laborious process if compared with tasks like image classification and object detection. This is especially true for remote sensing applications that not only work with extremely high spatial resolution data but also commonly require the knowledge of experts of the area to perform the manual labeling. Data augmentation techniques help to improve deep learning models under the circumstance of few and imbalanced labeled samples. In this work, we propose a novel data augmentation method focused on exploring the spatial context of remote sensing semantic segmentation. This method, ChessMix, creates new synthetic images from the existing training set by mixing transformed mini-patches across the dataset in a chessboard-like grid. ChessMix prioritizes patches with more examples of the rarest classes to alleviate the imbalance problems. The results in three diverse well-known remote sensing datasets show that this is a promising approach that helps to improve the networks' performance, working especially well in datasets with few available data. The results also show that ChessMix is capable of improving the segmentation of objects with few labeled pixels when compared to the most common data augmentation methods widely used.


title: Factorizations in evaluation monoids of Laurent semirings

abstract: For a positive real number $\alpha$, let $\mathbb{N}_0[\alpha,\alpha^{-1}]$ be the semiring of all real numbers $f(\alpha)$ for $f(x)$ lying in $\mathbb{N}_0[x,x^{-1}]$, which is the semiring of all Laurent polynomials over the set of nonnegative integers $\mathbb{N}_0$. In this paper, we study various factorization properties of the additive structure of $\mathbb{N}_0[\alpha, \alpha^{-1}]$. We characterize when $\mathbb{N}_0[\alpha, \alpha^{-1}]$ is atomic. Then we characterize when $\mathbb{N}_0[\alpha, \alpha^{-1}]$ satisfies the ascending chain condition on principal ideals in terms of certain well-studied factorization properties. Finally, we characterize when $\mathbb{N}_0[\alpha, \alpha^{-1}]$ satisfies the unique factorization property and show that, when this is not the case, $\mathbb{N}_0[\alpha, \alpha^{-1}]$ has infinite elasticity.


title: Mars' atmospheric neon suggests volatile-rich primitive mantle

abstract: Martian atmospheric neon (Ne) has been detected by Viking and also found as trapped gas in Martian meteorites, though its abundance and isotopic composition have not been well determined. Because the timescale of Ne loss via atmospheric escape estimated from recent measurements with MAVEN is short (0.6--1 $\times$ 10$^8$ years), the abundance and isotope composition of Martian atmospheric Ne reflect recent atmospheric gas supply mostly from volcanic degassing. Thus, it can serve as a probe for the volatile content of the interior. Here we show that the tentatively-informed atmospheric Ne abundance suggests recent active volcanism and the mantle being richer in Ne than Earth's mantle today by more than a factor of 5--80. The estimated mantle Ne abundance requires efficient solar nebular gas capture or accretion of Ne-rich materials such as solar-wind-implanted dust in the planet formation stage, both of which provide important constraints on the abundance of other volatile elements in the interior and the accretion history of Mars. More precise determination of atmospheric Ne abundance and isotopic composition by in situ analysis or Mars sample return is crucial for distinguishing the possible origins of Ne.


title: Terahertz waves generation using the isomorphs of PPKTP crystal: a theoretical investigation

abstract: Highly efficient terahertz (THz) wave sources based on difference frequency generation (DFG) process in nonlinear optical crystals play an important role for the applications of THz wave. In order to find more novel nonlinear crystals, here we theoretically investigate the generation of THz wave using the isomorphs of periodically poled $\mathrm{KTiOPO_4}$ (PPKTP), including periodically poled RTP, KTA, RTA and CTA. By solving the cascaded difference frequency coupled wave equations, it is found that the intensities of the THz wave generated from the cascaded difference frequency processes are improved by 5.27, 2.87, 2.82, 3.03, and 2.76 times from the non-cascaded cases for KTP, RTP, KTA, RTA and CTA, respectively. The effects of the crystal absorption, the phase mismatch and the pump intensity are also analyzed in detail. This study might help to provide a stronger THz radiation source based on the nonlinear crystals.


title: TPH-YOLOv5: Improved YOLOv5 Based on Transformer Prediction Head for Object Detection on Drone-captured Scenarios

abstract: Object detection on drone-captured scenarios is a recent popular task. As drones always navigate in different altitudes, the object scale varies violently, which burdens the optimization of networks. Moreover, high-speed and low-altitude flight bring in the motion blur on the densely packed objects, which leads to great challenge of object distinction. To solve the two issues mentioned above, we propose TPH-YOLOv5. Based on YOLOv5, we add one more prediction head to detect different-scale objects. Then we replace the original prediction heads with Transformer Prediction Heads (TPH) to explore the prediction potential with self-attention mechanism. We also integrate convolutional block attention model (CBAM) to find attention region on scenarios with dense objects. To achieve more improvement of our proposed TPH-YOLOv5, we provide bags of useful strategies such as data augmentation, multiscale testing, multi-model integration and utilizing extra classifier. Extensive experiments on dataset VisDrone2021 show that TPH-YOLOv5 have good performance with impressive interpretability on drone-captured scenarios. On DET-test-challenge dataset, the AP result of TPH-YOLOv5 are 39.18%, which is better than previous SOTA method (DPNetV3) by 1.81%. On VisDrone Challenge 2021, TPHYOLOv5 wins 5th place and achieves well-matched results with 1st place model (AP 39.43%). Compared to baseline model (YOLOv5), TPH-YOLOv5 improves about 7%, which is encouraging and competitive.


title: Learning-based Predictive Beamforming for Integrated Sensing and Communication in Vehicular Networks

abstract: This paper investigates the integrated sensing and communication (ISAC) in vehicle-to-infrastructure (V2I) networks. To realize ISAC, an effective beamforming design is essential which however, highly depends on the availability of accurate channel tracking requiring large training overhead and computational complexity. Motivated by this, we adopt a deep learning (DL) approach to implicitly learn the features of historical channels and directly predict the beamforming matrix to be adopted for the next time slot to maximize the average achievable sum-rate of an ISAC system. The proposed method can bypass the need of explicit channel tracking process and reduce the signaling overhead significantly. To this end, a general sum-rate maximization problem with Cramer-Rao lower bounds (CRLBs)-based sensing constraints is first formulated for the considered ISAC system taking into account the multiple access interference. Then, by exploiting the penalty method, a versatile unsupervised DL-based predictive beamforming design framework is developed to address the formulated design problem. As a realization of the developed framework, a historical channels-based convolutional long short-term memory (LSTM) network (HCL-Net) is devised for predictive beamforming in the ISAC-based V2I network. Specifically, the convolution and LSTM modules are successively adopted in the proposed HCL-Net to exploit the spatial and temporal dependencies of communication channels to further improve the learning performance. Finally, simulation results show that the proposed predictive method not only guarantees the required sensing performance, but also achieves a satisfactory sum-rate that can approach the upper bound obtained by the genie-aided scheme with the perfect instantaneous channel state information available.


title: A new expansion of planetary disturbing function and applications to interior, co-orbital and exterior resonances with planets

abstract: In this study, a new expansion of planetary disturbing function is developed for describing the resonant dynamics of minor bodies with arbitrary inclinations and semimajor axis ratios. In practice, the disturbing function is expanded around circular orbits in the first step and then, in the second step, the resulting mutual interaction between circular orbits is expanded around a reference point. As usual, the resulting expansion is presented in the Fourier series form, where the force amplitudes are dependent on the semimajor axis, eccentricity and inclination and the harmonic arguments are linear combinations of the mean longitude, longitude of pericenter and longitude of ascending node of each mass. The resulting new expansion is valid for arbitrary inclinations and semimajor axis ratios. In the case of mean motion resonant configuration, the disturbing function can be easily averaged to produce the analytical expansion of resonant disturbing function. Based on the analytical expansion, the Hamiltonian model of mean motion resonances is formulated, and the resulting analytical developments are applied to Jupiter's inner and co-orbital resonances and Neptune's exterior resonances. Analytical expansion is validated by comparing the analytical results with the associated numerical outcomes.


title: The Phase Diagram of BPS Black Holes in AdS$_5$

abstract: Motivated by recent studies of supersymmetric black holes, we revisit the phase diagram of AdS black holes, whether BPS or not, with particular emphasis on the role of rotation. We develop BPS thermodynamics systematically and, in many explicit examples, we study its striking similarities with more familiar AdS black holes, as well as some differences. We highlight an important fugacity that preserves BPS saturation but is not captured by the supersymmetric index.


title: An improved search for the electric dipole moment of the $\tau$ lepton

abstract: We report a measurement of the electric dipole moment of the $\tau$ lepton ($d_\tau$) using an 833~fb$^{-1}$ data sample collected near the $\Upsilon(4S)$ resonance, with the Belle detector at the KEKB asymmetric-energy $e^+ e^-$ collider. Using an optimal observable method, we obtain the real and imaginary parts of $d_\tau$ as ${\rm Re}(d_\tau) = ( -0.62 \pm 0.63 ) \times 10^{-17} ~e{\rm cm}$ and ${\rm Im}(d_\tau) = ( -0.40 \pm 0.32 ) \times 10^{-17} ~e{\rm cm}$, respectively. These results are consistent with null electric dipole moment at the present level of experimental sensitivity and improve the sensitivity by about a factor of three.


title: Vision-Language Navigation: A Survey and Taxonomy

abstract: Vision-Language Navigation (VLN) tasks require an agent to follow human language instructions to navigate in previously unseen environments. This challenging field involving problems in natural language processing, computer vision, robotics, etc., has spawn many excellent works focusing on various VLN tasks. This paper provides a comprehensive survey and an insightful taxonomy of these tasks based on the different characteristics of language instructions in these tasks. Depending on whether the navigation instructions are given for once or multiple times, this paper divides the tasks into two categories, i.e., single-turn and multi-turn tasks. For single-turn tasks, we further subdivide them into goal-oriented and route-oriented based on whether the instructions designate a single goal location or specify a sequence of multiple locations. For multi-turn tasks, we subdivide them into passive and interactive tasks based on whether the agent is allowed to question the instruction or not. These tasks require different capabilities of the agent and entail various model designs. We identify progress made on the tasks and look into the limitations of existing VLN models and task settings. Finally, we discuss several open issues of VLN and point out some opportunities in the future, i.e., incorporating knowledge with VLN models and implementing them in the real physical world.


title: Identification of Peer Effects using Panel Data

abstract: We provide new identification results for panel data models with peer effects operating through unobserved individual heterogeneity. The results apply for general network structures governing peer interactions and allow for correlated effects. Identification hinges on a conditional mean restriction requiring exogenous mobility of individuals between groups over time. We apply our method to surgeon-hospital-year data to study take-up of keyhole surgery for cancer, finding a positive effect of the average individual heterogeneity of other surgeons practicing in the same hospital


title: A study of computational and conceptual complexities of compartment and agent based models

abstract: The ongoing COVID-19 pandemic highlights the essential role of mathematical models in understanding the spread of the virus along with a quantifiable and science-based prediction of the impact of various mitigation measures. Numerous types of models have been employed with various levels of success. This leads to the question of what kind of a mathematical model is most appropriate for a given situation. We consider two widely used types of models: equation-based models (such as standard compartmental epidemiological models) and agent-based models. We assess their performance by modeling the spread of COVID-19 on the Hawaiian island of Oahu under different scenarios. We show that when it comes to information crucial to decision making, both models produce very similar results. At the same time, the two types of models exhibit very different characteristics when considering their computational and conceptual complexity. Consequently, we conclude that choosing the model should be mostly guided by available computational and human resources.


title: Magnetoelastic anisotropy in Heusler-type Mn$_{2-\delta}$CoGa$_{1+\delta}$ films

abstract: Perpendicular magnetization is essential for high-density memory application using magnetic materials. High-spin polarization of conduction electrons is also required for realizing large electric signals from spin-dependent transport phenomena. Heusler alloy is a well-known material class showing the half-metallic electronic structure. However, its cubic lattice nature favors in-plane magnetization and thus minimizes the perpendicular magnetic anisotropy (PMA), in general. This study focuses on an inverse-type Heusler alloy, Mn$_{2-\delta}$CoGa$_{1+\delta}$ (MCG) with a small off-stoichiometry ($\delta$) , which is expected to be a half-metallic material. We observed relatively large uniaxial magnetocrystalline anisotropy constant ($K_\mathrm{u}$) of the order of 10$^5$ J/m$^3$ at room temperature in MCG films with a small tetragonal distortion of a few percent. A positive correlation was confirmed between the $c/a$ ratio of lattice constants and $K_\mathrm{u}$. Imaging of magnetic domains using Kerr microscopy clearly demonstrated a change in the domain patterns along with $K_\mathrm{u}$. X-ray magnetic circular dichroism (XMCD) was employed using synchrotron radiation soft x-ray beam to get insight into the origin for PMA. Negligible angular variation of orbital magnetic moment ($\Delta m_\mathrm{orb}$) evaluated using the XMCD spectra suggested a minor role of the so-called Bruno's term to $K_\mathrm{u}$. Our first principles calculation reasonably explained the small $\Delta m_\mathrm{orb}$ and the positive correlation between the $c/a$ ratio and $K_\mathrm{u}$. The origin of the magnetocrystalline anisotropy was discussed based on the second-order perturbation theory in terms of the spin--orbit coupling, claiming that the mixing of the occupied $\uparrow$- and the unoccupied $\downarrow$-spin states is responsible for the PMA of the MCG films.


title: Dynamo and the Adiabatic Invariant

abstract: The paper considers dynamo generated by a shallow fluid layer in a celestial body (planet or star). This dynamo is based on the extra invariant for interacting magnetic Rossby waves. The magnetohydrodynamics (MHD) is linearized on the background of strong toroidal magnetic field. The extra invariant is used to show that the background field is maintained.

categories: physics.plasm-ph astro-ph.EP nlin.AO physics.flu-dyn


title: Dynamic Structural Clustering on Graphs

abstract: Structural Clustering ($DynClu$) is one of the most popular graph clustering paradigms. In this paper, we consider $StrClu$ under two commonly adapted similarities, namely Jaccard similarity and cosine similarity on a dynamic graph, $G = \langle V, E\rangle$, subject to edge insertions and deletions (updates). The goal is to maintain certain information under updates, so that the $StrClu$ clustering result on~$G$ can be retrieved in $O(|V| + |E|)$ time, upon request. The state-of-the-art worst-case cost is $O(|V|)$ per update; we improve this update-time bound significantly with the $\rho$-approximate notion. Specifically, for a specified failure probability, $\delta^*$, and every sequence of $M$ updates (no need to know $M$'s value in advance), our algorithm, $DynELM$, achieves $O(\log^2 |V| + \log |V| \cdot \log \frac{M}{\delta^*})$ amortized cost for each update, at all times in linear space. Moreover, $DynELM$ provides a provable "sandwich" guarantee on the clustering quality at all times after \emph{each update} with probability at least $1 - \delta^*$. We further develop $DynELM$ into our ultimate algorithm, $DynStrClu$, which also supports cluster-group-by queries. Given $Q\subseteq V$, this puts the non-empty intersection of $Q$ and each $StrClu$ cluster into a distinct group. $DynStrClu$ not only achieves all the guarantees of $DynELM$, but also runs cluster-group-by queries in $O(|Q|\cdot \log |V|)$ time. We demonstrate the performance of our algorithms via extensive experiments, on 15 real datasets. Experimental results confirm that our algorithms are up to three orders of magnitude more efficient than state-of-the-art competitors, and still provide quality structural clustering results. Furthermore, we study the difference between the two similarities w.r.t. the quality of approximate clustering results.