title: Reconciling cosmic-ray transport theory with phenomenological models motivated by Milky-Way data

abstract: Phenomenological models of cosmic-ray (CR) transport in the Milky Way (MW) can reproduce a wide range of observations assuming that CRs scatter off of magnetic-field fluctuations with spectrum $\propto k^{-\delta}$ and $\delta \sim [1.4,1.67]$. We study the extent to which such models can be reconciled with current microphysical theories of CR transport, specifically self-confinement due to the streaming instability and/or extrinsic turbulence due to a cascade of MHD fast modes. We first review why it is that on their own neither theory is compatible with observations. We then highlight that CR transport is a strong function of local plasma conditions in the multi-phase interstellar medium (ISM), and may be diffusive due to turbulence in some regions and streaming due to self-confinement in others. A multi-phase combination of scattering mechanisms can in principle reproduce the main trends in the proton spectrum and the boron-to-carbon ratio (B/C). However, models with a combination of scattering by self-excited waves and fast-mode turbulence require significant fine-tuning due to fast-mode damping, unlike phenomenological models that assume undamped Kolmogorov turbulence. The assumption that fast modes follow a weak cascade is also not well justified theoretically, as the weak cascade is suppressed by wave steepening and weak-shock dissipation even in subsonic turbulence. These issues suggest that there may be a significant theoretical gap in our understanding of MHD turbulence. We discuss a few topics at the frontier of MHD turbulence theory that bear on this (possible) gap and that may be relevant for CR scattering.


title: Quantum nucleation of gas bubbles in liquid helium

abstract: Liquid helium under negative pressures represents a unique possibility for studying nucleation and growth dynamics of cavities at low temperatures down to absolute zero. We analyze the growth dynamics of cavities and determine the temperature behavior of quantum cavitation rate and the crossover temperature between the classical and quantum cavitation mechanisms. The energy dissipation processes, which differ in kind in normal fluid $^3$He and superfluid $^4$He, result in the different temperature behavior of cavitation rate in $^3$He and $^4$He below the thermal-quantum crossover temperature.


title: The case of an N-gon

abstract: We construct the indefinite theta series attached to N-gons in the symmetric space of an indefinite inner product space of signature (m-2,2) following the suggestions of section C in the recent paper of Alexandrov, Banerjee, Manschot, and Pioline. We prove the termwise absolute convergence of the holomorphic mock modular part of these series and also obtain an interpretation of the coefficients of this part as linking numbers. Thus we prove the convergence conjecture of ABMP provided none of the vectors in the collection CC={C_1,..., C_N} is a null vector. The use of linking numbers and a homotopy argument eliminates the need for an explicit parametrization of a surface S spanning the N-gon that was used in an essential way in our previous work. We indicate how our method could be carried over to a more general situation for signature (m-q,q) where higher homotopy groups are now involved. In the last section, we apply the method to the case of a dodecahedral cell in the symmetric space of a quadratic form of signature (m-3,3).


title: Design, Fabrication, Characterization and Reliability Study of CMOS-MEMS Lorentz-Force Magnetometers

abstract: This article presents several design techniques to fabricate micro-electro-mechanical systems (MEMS) using standard complementary metal-oxide semiconductor (CMOS) processes. They were applied to fabricate high yield CMOS-MEMS shielded Lorentz-force magnetometers (LFM). The multilayered metals and oxides of the back-end-of-line (BEOL), normally used for electronic routing, comprise the structural part of the MEMS. The most important fabrication challenges, modeling approaches and design solutions are discussed. Equations that predict the Q factor, sensitivity, Brownian noise and resonant frequency as a function of temperature, gas pressure and design parameters are presented and validated in characterization tests. A number of the fabricated magnetometers were packaged into Quad Flat No-leads (QFN) packages. We show this process can achieve yields above 95 % when the proper design techniques are adopted. Despite CMOS not being a process for MEMS manufacturing, estimated performance (sensitivity and noise level) is similar or superior to current commercial magnetometers and others built with MEMS processes. Additionally, typical offsets present in Lorentz-force magnetometers were prevented with a shielding electrode, whose efficiency is quantified. Finally, several reliability test results are presented, which demonstrate the robustness against high temperatures, magnetic fields and acceleration shocks.


title: The Point-to-Set Principle and the Dimensions of Hamel Bases

abstract: We prove that every real number in [0,1] is the Hausdorff dimension of a Hamel basis of the vector space of reals over the field of rationals. The logic of our proof is of particular interest. The statement of our theorem is classical; it does not involve the theory of computing. However, our proof makes essential use of algorithmic fractal dimension--a computability-theoretic construct--and the point-to-set principle of J. Lutz and N. Lutz (2018).


title: Quantifying nonlocality: how outperforming local quantum codes is expensive

abstract: Quantum low-density parity-check (LDPC) codes are a promising avenue to reduce the cost of constructing scalable quantum circuits. However, it is unclear how to implement these codes in practice. Seminal results of Bravyi & Terhal, and Bravyi, Poulin & Terhal have shown that quantum LDPC codes implemented through local interactions obey restrictions on their dimension $k$ and distance $d$. Here we address the complementary question of how many long-range interactions are required to implement a quantum LDPC code with parameters $k$ and $d$. In particular, in 2D we show that a quantum LDPC with distance $n^{1/2 + \epsilon}$ code requires $\Omega(n^{1/2 + \epsilon})$ interactions of length $\widetilde{\Omega}(n^{\epsilon})$. Further a code satisfying $k \propto n$ with distance $d \propto n^\alpha$ requires $\widetilde{\Omega}(n)$ interactions of length $\widetilde{\Omega}(n^{\alpha/2})$. Our results are derived using bounds on quantum codes from graph metrics. As an application of these results, we consider a model called a stacked architecture, which has previously been considered as a potential way to implement quantum LDPC codes. In this model, although most interactions are local, a few of them are allowed to be very long. We prove that limited long-range connectivity implies quantitative bounds on the distance and code dimension.


title: Anisotropic Flow Measurements of Identified Particles in the STAR Experiment

abstract: We report results of $v_1$ and $v_2$ for identified hadrons ($\pi^{\pm}, K^{\pm}, K_{S}^{0}, p, \phi$ and $\Lambda$) from Au+Au collisions at $\sqrt{s_{NN}}$ = 3 GeV and $v_2$ for $\pi^{\pm}, K^{\pm}, p$ and $\bar{p}$ at $\sqrt{s_{NN}}$ = 27 and 54.4 GeV using the STAR detector at RHIC. In Au+Au collisions at high energies, one finds that the values of $v_2$ are all positive and the number-of-constituent-quark (NCQ) scaling holds. On the other hand, results from collisions at 3 GeV, the midrapidity $v_2$ is negative for all hadrons and NCQ scaling is absent. In addition, the midrapidity $v_1$ slopes for all hadrons are found to be positive. Furthermore, the features of negative $v_2$ and positive $v_1$ slope at 3 GeV can be reproduced by calculations with baryonic mean-field potential in transport model. These results imply that in 3 GeV Au+Au collisions, the medium is characterized by baryonic interactions.


title: The dynamics of the TRAPPIST-1 system in the context of its formation

abstract: TRAPPIST-1 is an 0.09 $M_{\odot}$ star, which harbours a system of seven Earth-sized planets. Two main features stand out: (i) all planets have similar radii, masses, and compositions; and (ii) all planets are in resonance. Previous works have outlined a pebble-driven formation scenario where planets of similar composition form sequentially at the H$_2$O snowline (${\sim}0.1$ au for this low-mass star). It was hypothesized that the subsequent formation and migration led to the current resonant configuration. Here, we investigate whether the sequential planet formation model is indeed capable to produce the present-day resonant configuration, characterized by its two-body and three-body mean motion resonances structure. We carry out N-body simulations, accounting for type-I migration, stellar tidal damping, disc eccentricity-damping, and featuring a migration barrier located at the disc's inner edge. Due to the sequential migration, planets naturally form a chain of first-order resonances. But to explain the period ratios of the b/c/d-system, which are presently in higher-order resonances, we find that planets b and c must have marched across the migration barrier, into the gas-free cavity, before the disc has dispersed. We investigate both an early and late cavity infall scenario and find that the early infall model best matches the constraints, as well as being more probable. After the dispersal of the gaseous disc, stellar tidal torque also contributes towards a modest separation of the inner system. We outline how the insights obtained in this work can be applied to aid the understanding of other compact resonant planet systems.


title: Enhanced dispersion in intermittent multiphase flow

abstract: Transport in multiphase flow through porous media plays a central role in many biological, geological, and engineered systems. Here, we use numerical simulations of transport in immiscible two-phase flow to investigate dispersion in multiphase porous media flow. While dispersion in the main flow direction is similar to that of a single-phase flow and is governed by the porous media structure, we find that transverse dispersion exhibits fundamentally different dynamics. The repeated activation and deactivation of different flow pathways under the effect of capillary forces lead to intermittent flow patterns, strongly enhancing dispersion. We show that the transverse dispersion is controlled by the length scale of fluid clusters, and thus inversely proportional to the square root of the Bond number, the ratio of the force driving the flow and the surface tension. This result leads to a new scaling law relating multiphase flow dynamics to transport properties, opening a range of fundamental and engineering applications.


title: An Efficient and Scalable Collection of Fly-inspired Voting Units for Visual Place Recognition in Changing Environments

abstract: State-of-the-art visual place recognition performance is currently being achieved utilizing deep learning based approaches. Despite the recent efforts in designing lightweight convolutional neural network based models, these can still be too expensive for the most hardware restricted robot applications. Low-overhead VPR techniques would not only enable platforms equipped with low-end, cheap hardware but also reduce computation on more powerful systems, allowing these resources to be allocated for other navigation tasks. In this work, our goal is to provide an algorithm of extreme compactness and efficiency while achieving state-of-the-art robustness to appearance changes and small point-of-view variations. Our first contribution is DrosoNet, an exceptionally compact model inspired by the odor processing abilities of the fruit fly, Drosophyla melanogaster. Our second and main contribution is a voting mechanism that leverages multiple small and efficient classifiers to achieve more robust and consistent VPR compared to a single one. We use DrosoNet as the baseline classifier for the voting mechanism and evaluate our models on five benchmark datasets, assessing moderate to extreme appearance changes and small to moderate viewpoint variations. We then compare the proposed algorithms to state-of-the-art methods, both in terms of precision-recall AUC results and computational efficiency.


title: Complexity of magnetic-field turbulence at reconnection exhausts in the solar wind at 1 AU

abstract: Magnetic reconnection is a complex mechanism that converts magnetic energy into particle kinetic energy and plasma thermal energy in space and astrophysical plasmas. In addition, magnetic reconnection and turbulence appear to be intimately related in plasmas. We analyze the magnetic-field turbulence at the exhaust of four reconnection events detected in the solar wind using the Jensen-Shannon complexity-entropy index. The interplanetary magnetic field is decomposed into the LMN coordinates using the hybrid minimum variance technique. The first event is characterized by an extended exhaust period that allows us to obtain the scaling exponents of higher-order structure functions of magnetic-field fluctuations. By computing the complexity-entropy index we demonstrate that a higher degree of intermittency is related to lower entropy and higher complexity in the inertial subrange. We also compute the complexity-entropy index of three other reconnection exhaust events. For all four events, the $B_L$ component of the magnetic field displays a lower degree of entropy and higher degree of complexity than the $B_M$ and $B_N$ components. Our results show that coherent structures can be responsible for decreasing entropy and increasing complexity within reconnection exhausts in magnetic-field turbulence.


title: Theoretical investigation of decoherence channels in athermal phonon sensors

abstract: The creation and evolution of nonequilibrium phonons is central in applications ranging from cosmological particle searches to decoherence processes in qubits. However, the fundamental understanding of decoherence pathways for athermal phonon distributions in solid-state systems remains an open question. Using first-principles calculations, we investigate the primary decay channels of athermal phonons in two technologically relevant semiconductors -- Si and GaAs. We quantify the contributions of anharmonic, isotopic, and interfacial scattering in these materials. From this, we construct a model to estimate the thermal power in a readout scheme as a function of time. We discuss the implication of our results on noise limitations in current phonon sensor designs and strategies for improving coherence in next-generation phonon sensors.


title: Robust Folding of Elastic Origami

abstract: Self-folding origami, structures that are engineered flat to fold into targeted, three-dimensional shapes, have many potential engineering applications. Though significant effort in recent years has been devoted to designing fold patterns that can achieve a variety of target shapes, recent work has also made clear that many origami structures exhibit multiple folding pathways, with a proliferation of geometric folding pathways as the origami structure becomes complex. The competition between these pathways can lead to structures that are programmed for one shape, yet fold incorrectly. To disentangle the features that lead to misfolding, we introduce a model of self-folding origami that accounts for the finite stretching rigidity of the origami faces and allows the computation of energy landscapes that lead to misfolding. We find that, in addition to the geometrical features of the origami, the finite elasticity of the nearly-flat origami configurations regulates the proliferation of potential misfolded states through a series of saddle-node bifurcations. We apply our model to one of the most common origami motifs, the symmetric "bird's foot," a single vertex with four folds. We show that though even a small error in programmed fold angles induces metastability in rigid origami, elasticity allows one to tune resilience to misfolding. In a more complex design, the "Randlett flapping bird," which has thousands of potential competing states, we further show that the number of actual observed minima is strongly determined by the structure's elasticity. In general, we show that elastic origami with both stiffer folds and stiffer faces self-folds better.


title: A new threshold model based on tropical mathematics reveals network backbones

abstract: Maxmin-$\omega$ is a new threshold model, where each node in a network waits for the arrival of states from a fraction $\omega$ of neighborhood nodes before processing its own state, and subsequently transmitting it to downstream nodes. Repeating this sequence of events leads to periodic behavior. We show that maxmin-$\omega$ reduces to a smaller, simpler, system represented by tropical mathematics, which forges a useful link between the nodal state update times and circuits in the network. Thus, the behavior of the system can be analyzed directly from the smaller network structure and is computationally faster. We further show that these reduced networks: (i) are not unique; they are dependent on the initialisation time, (ii) tend towards periodic orbits of networks -- "attractor networks." In light of these features, we vary the initial condition and $\omega$ to obtain statistics on types of attractor networks. The results suggest the case $\omega =0.5$ to give the most stable system. We propose that the most prominent attractor networks of nodes and edges may be seen as a 'backbone' of the original network. We subsequently provide an algorithm to construct attractor networks, the main result being that they must necessarily contain circuits that are deemed critical and that can be identified without running the maxmin-$\omega$ system. Finally, we apply this work to the C. elegans neuron network, giving insight into its function and similar networks. This novel area of application for tropical mathematics adds to its theory, and provides a new way to identify important edges and circuits of a network.


title: Dynamical Spin-Orbit-Based Spin Transistor

abstract: Spin-orbit interaction (SOI) has been a key tool to steer and manipulate spin-dependent transport properties in two-dimensional electron gases. Here we demonstrate how spin currents can be created and efficiently read out in nano- or mesoscale conductors with time-dependent and spatially inhomogeneous Rashba SOI. Invoking an underlying non-Abelian SU(2) gauge structure we show how time-periodic spin-orbit fields give rise to spin electric forces and enable the generation of pure spin currents of the order of several hundred nano-Amperes. In a complementary way, by combining gauge transformations with "hidden" Onsager relations, we exploit spatially inhomogeneous Rashba SOI to convert spin currents (back) into charge currents. In combining both concepts, we devise a spin transistor that integrates efficient spin current generation, by employing dynamical SOI, with its experimentally feasible detection via conversion into charge signals. We derive general expressions for the respective spin- and charge conductance, covering large parameter regimes of SOI strength and driving frequencies, far beyond usual adiabatic approaches such as the frozen scattering matrix approximation. We check our analytical expressions and approximations with full numerical spin-dependent transport simulations and demonstrate that the predictions hold true in a wide range from low to high driving frequencies.


title: Scalable Fact-checking with Human-in-the-Loop

abstract: Researchers have been investigating automated solutions for fact-checking in a variety of fronts. However, current approaches often overlook the fact that the amount of information released every day is escalating, and a large amount of them overlap. Intending to accelerate fact-checking, we bridge this gap by grouping similar messages and summarizing them into aggregated claims. Specifically, we first clean a set of social media posts (e.g., tweets) and build a graph of all posts based on their semantics; Then, we perform two clustering methods to group the messages for further claim summarization. We evaluate the summaries both quantitatively with ROUGE scores and qualitatively with human evaluation. We also generate a graph of summaries to verify that there is no significant overlap among them. The results reduced 28,818 original messages to 700 summary claims, showing the potential to speed up the fact-checking process by organizing and selecting representative claims from massive disorganized and redundant messages.


title: Verification of Approximate Opacity via Barrier Certificates

abstract: This paper is motivated by the increasing security concerns of cyber-physical systems. Here, we develop a discretization-free verification scheme targeting an information-flow security property, called approximate initial-state opacity, for the class of discrete-time control systems. We propose notions of so-called augmented control barrier certificates in conjunction with specified regions of interest capturing the initial and secret sets of the system. Sufficient conditions for (the lack of) approximate initial-state opacity of discrete-time control systems are proposed based on the existence of the proposed barrier certificates. We further present an efficient computation method by casting the conditions for barrier certificates as sum-of-squares programming problems. The effectiveness of the proposed results is illustrated through two numerical examples.


title: Some non-standard biparametric Poincar\'e type inequalities through harmonic analysis

abstract: We show some non-standard Poincar\'e type estimates in the biparametric setting with appropriate weights. We will derive these results using variants from classical estimates exploiting the interplay between maximal functions and fractional integrals. We also provide a sharper result by using extrapolation techniques.


title: Charmonium in nuclear matter and nuclei

abstract: We present results for the $\eta_c$-nucleus bound state energies for various nuclei using an effective Lagrangians approach. The attractive potentials for the $\eta_c$ in the nuclear medium originate from the medium-modified intermediate $D D^{*}$ state in the $\eta_c$ self energy, using the local density approximation. Our results suggest that the $\eta_c$ should form bound states with all the nuclei considered


title: Cartesian Frames

abstract: We introduce a novel framework, the theory of Cartesian frames (CF), that gives powerful tools for manipulating sets of acts. The CF framework takes as its most fundamental building block that an agent can freely choose from a set of available actions. The framework uses the mathematics of Chu spaces to develop a calculus of those sets of actions, how those actions change at various levels of description, and how different agents' actions can combine when agents work in concert. We discuss how this framework might provide an illuminating perspective on issues in decision theory and formal epistemology.


title: Extinction time in growth models subject to geometric catastrophes

abstract: Recently, different dispersion strategies in population models subject to geometric catastrophes have been considered as strategies to improve the chance of po\-pu\-lation's survival. Such dispersion strategies have been contrasted with the strategy where there is no dispersion, comparing the probabilities of survival. In this article, we contrast survival strategies when extinction occurs almost surely, evaluating which strategy prolongs population's life span. Our results allow one to analyze what is the best strategy for survival based on parameters as the probability that each individual exposed to catastrophe survives, the growth rate of the colony, the type of dispersion and the spatial restrictions.


title: Dark and luminous mass components of Omega Centauri with stellar kinematics

abstract: We combine proper motion data from $Gaia$ EDR3 and HST with line-of-sight velocity data to study the stellar kinematics of the Omega Centauri globular cluster. Using a steady-state, axisymmetric dynamical model, we measure the distribution of both the dark and luminous mass components. Assuming both Gaussian and NFW mass profiles, depending on the dataset, we measure an integrated mass of $\lesssim 10^6$ M$_\odot$ within the Omega Centauri half-light radius for a dark component that is distinct from the luminous stellar component. For the HST and radial velocity data, models with a non-luminous mass component are strongly statistically preferred relative to a stellar mass-only model with a constant mass-to-light ratio. While a compact core of stellar remnants may account for a dynamical mass up to $\sim 5 \times 10^5$ M$_\odot$, they likely cannot explain the higher end of the range. This leaves open the possibility that this non-luminous dynamical mass component is comprised of non-baryonic dark matter. In comparison to the dark matter distributions around dwarf spheroidal galaxies, the Omega Centauri dark mass component is much more centrally concentrated. Interpreting the non-luminous mass distribution as particle dark matter, we use these results to obtain the J-factor, which sets the sensitivity to the annihilation cross section. For the datasets considered, the range of median J-factors is $\sim 10^{22} - 10^{24}$ GeV$^2$ cm$^{-5}$, which is larger than that obtained for any dwarf spheroidal galaxy.


title: A discontinuous Galerkin pressure correction scheme for the incompressible Navier-Stokes equations: stability and convergence

abstract: A discontinuous Galerkin pressure correction numerical method for solving the incompressible Navier-Stokes equations is formulated and analyzed. We prove unconditional stability of the propose scheme. Convergence of the discrete velocity is established by deriving a priori error estimates. Numerical results verify the convergence rates.


title: The first retrieval of a substellar subdwarf: A cloud-free SDSS J125637.13-022452.4

abstract: We present the first retrieval analysis of a substellar subdwarf, SDSS J125637.13-022452.4 (SDSS J1256-0224), using the Brewster retrieval code base. We find SDSS J1256-0224 is best fit by a cloud-free model with an ion (neutral H, H-, and electron) abundance corresponding to ion [Fe/H]=-1.5. However, this model is indistinguishable from a cloud-free model with ion [Fe/H]=-2.0 and a cloud-free model with ion Fe/H]=-1.5 assuming a subsolar carbon-to-oxygen ratio. We are able to constrain abundances for water, FeH, and CrH, with an inability to constrain any carbon-bearing species likely due to the low-metallicity of SDSS J1256-0224. We also present an updated spectral energy distribution (SED) and semi-empirical fundamental parameters. Our retrieval- and SED-based fundamental parameters agree with the Baraffe low-metallicity evolutionary models. From examining our "rejected" models (those with $\Delta$BIC>45), we find that we are able to retrieve gas abundances consistent with those of our best-fitting model. We find the cloud in these poorer fitting "cloudy" models is either pushed to the bottom of the atmosphere or made optically thin.


title: Facile equilibration of well-entangled semiflexible bead-spring polymer melts

abstract: The widely used double-bridging hybrid (DBH) method for equilibrating simulated entangled polymer melts [R. Auhl et al., J. Chem. Phys. v. 119, p. 12718, 2003] loses its effectiveness as chain stiffness increases into the semiflexible regime because the energy barriers associated with double-bridging Monte Carlo moves become prohibitively high. Here we overcome this issue by combining DBH with the use of core-softened pair potentials. This reduces the energy barriers substantially, allowing us to equilibrate melts with $N \simeq 40 N_e$ and chain stiffnesses all the way up to the isotropic-nematic transition using simulations of no more than 100 million timesteps. For semiflexible chains, our method is several times faster than standard DBH; we exploit this speedup to develop improved expressions for Kremer-Grest melts' chain-stiffness-dependent Kuhn length $\ell_K$ and entanglement length $N_e$.


title: A Benchmark Comparison of Visual Place Recognition Techniques for Resource-Constrained Embedded Platforms

abstract: Visual Place Recognition (VPR) has been a subject of significant research over the last 15 to 20 years. VPR is a fundamental task for autonomous navigation as it enables self-localization within an environment. Although robots are often equipped with resource-constrained hardware, the computational requirements of and effects on VPR techniques have received little attention. In this work, we present a hardware-focused benchmark evaluation of a number of state-of-the-art VPR techniques on public datasets. We consider popular single board computers, including ODroid, UP and Raspberry Pi 3, in addition to a commodity desktop and laptop for reference. We present our analysis based on several key metrics, including place-matching accuracy, image encoding time, descriptor matching time and memory needs. Key questions addressed include: (1) How does the performance accuracy of a VPR technique change with processor architecture? (2) How does power consumption vary for different VPR techniques and embedded platforms? (3) How much does descriptor size matter in comparison to today's embedded platforms' storage? (4) How does the performance of a high-end platform relate to an on-board low-end embedded platform for VPR? The extensive analysis and results in this work serve not only as a benchmark for the VPR community, but also provide useful insights for real-world adoption of VPR applications.


title: Rational approximations of irrational numbers

abstract: Given quantities $\Delta_1,\Delta_2,\dots\geqslant 0$, a fundamental problem in Diophantine approximation is to understand which irrational numbers $x$ have infinitely many reduced rational approximations $a/q$ such that $|x-a/q|<\Delta_q$. Depending on the choice of $\Delta_q$ and of $x$, this question may be very hard. However, Duffin and Schaeffer conjectured in 1941 that if we assume a "metric" point of view, the question is governed by a simple zero--one law: writing $\varphi$ for Euler's totient function, we either have $\sum_{q=1}^\infty \varphi(q)\Delta_q=\infty$ and then almost all irrational numbers (in the Lebesgue sense) are approximable, or $\sum_{q=1}^\infty\varphi(q)\Delta_q<\infty$ and almost no irrationals are approximable. We present the history of the Duffin--Schaeffer conjecture and the main ideas behind the recent work of Koukoulopoulos--Maynard that settled it.


title: Considerations on anomalous photon and Z-boson self-couplings from the Born-Infeld weak hypercharge action

abstract: We investigate the effects of the Born-Infeld action on the Abelian sector of the electroweak model. The consequence of this approach is the emergence of anomalous couplings in the neutral sector of the $Z$-gauge boson and photon. These new couplings consist of quartic interactions of the photon with the $Z$-particle, as for example, three-photon-and-one-$Z$ vertex. With that, we obtain the decay width of $Z \to 3\,\gamma$ from which we impose a bound on the Born-Infeld parameter. Other bounds are also obtained from the photon quartic couplings. Subsequently, we consider the presence of an external uniform magnetic field in connection with this Born-Infeld weak hypercharge model. The magnetic background field yields new kinematic effects, like the kinetic mixing between the photon and the $Z$-boson, and we obtain thereby the corresponding dispersion relations for the mixed photon-$Z$-particle system. Finally, we calculate the lowest-order modifications to the interaction energy for the anomalous coupling $3\gamma-Z$, within the framework of the gauge-invariant but path-dependent variables formalism. Our results show that the interaction energy contains a linear term leading to the confinement of static probe charges. With the help of the potential that comes out, interparticle forces are estimated.


title: Prospects of SETI by small size optical telescopes

abstract: In the present manuscript we consider the possibility of SETI by the small size (with diameters less than 1 m) optical telescopes. Calculations are performed for typical parameters of the mentioned type of telescopes. In particular, we show that the techno-signatures of Type-2.x and Type-3.x civilizations might be detected. It is demonstrated that it is possible to detect as hot megastructures (up to 4000 K) built around main sequence stars and pulsars as well as von Neumann extraterrestrial probes.


title: The Sharp Erd\H{o}s-Tur\'an Inequality

abstract: Erd\H{o}s and Tur\'an proved a classical inequality on the distribution of roots for a complex polynomial in 1950, depicting the fundamental interplay between the size of the coefficients of a polynomial and the distribution of its roots on the complex plane. Various results have been dedicated to improving the constant in this inequality, while the optimal constant remains open. In this paper, we give the optimal constant, i.e., prove the sharp Erd\H{o}s-Tur\'an inequality. To achieve this goal, we reformulate the inequality into an optimization problem, whose equilibriums coincide with a class of energy minimizers with the logarithmic interaction and external potentials. This allows us to study their properties by taking advantage of the recent development of energy minimization and potential theory, and to give explicit constructions via complex analysis. Finally the sharp Erd\H{o}s-Tur\'an inequality is obtained based on a thorough understanding of these equilibrium distributions.


title: SoK: Cryptographic Confidentiality of Data on Mobile Devices

abstract: Mobile devices have become an indispensable component of modern life. Their high storage capacity gives these devices the capability to store vast amounts of sensitive personal data, which makes them a high-value target: these devices are routinely stolen by criminals for data theft, and are increasingly viewed by law enforcement agencies as a valuable source of forensic data. Over the past several years, providers have deployed a number of advanced cryptographic features intended to protect data on mobile devices, even in the strong setting where an attacker has physical access to a device. Many of these techniques draw from the research literature, but have been adapted to this entirely new problem setting. This involves a number of novel challenges, which are incompletely addressed in the literature. In this work, we outline those challenges, and systematize the known approaches to securing user data against extraction attacks. Our work proposes a methodology that researchers can use to analyze cryptographic data confidentiality for mobile devices. We evaluate the existing literature for securing devices against data extraction adversaries with powerful capabilities including access to devices and to the cloud services they rely on. We then analyze existing mobile device confidentiality measures to identify research areas that have not received proper attention from the community and represent opportunities for future research.


title: Intermediate spin state and the B1-B2 transition in ferropericlase at tera-Pascal pressures

abstract: Ferropericlase (fp), (Mg$_{1-x}$Fe$_x$)O, the second most abundant mineral in the Earth's lower mantle, is expected to be an essential component of super-Earths' mantles. Here we present an ab initio investigation of the structure and magnetic ground state of fp up to $\sim$ 3 TPa with iron concentrations (x$_{Fe}$) varying from 0.03 to 0.12. Calculations were performed using LDA+U$_{sc}$ and PBE exchange-correlation functionals to elucidate the pressure range for which the Hubbard U (U$_{sc}$) is required. Similar to the end-members FeO and MgO, fp also undergoes a B1 to B2 phase transition that should be essential for modeling the structure and dynamics of super-Earths' mantles. This structural transition involves a simultaneous change in magnetic state from a low spin (LS) B1 phase with iron total spin S=0 to an intermediate spin (IS) B2 phase with S=1. This is a rare form of pressure/strain-induced magnetism produced by local cation coordination changes. Phonon calculations confirm the dynamical stability of the iron B2-IS state. Free energy calculations are then carried out including vibrational effects and electronic and magnetic entropy contributions. The phase diagram is then obtained for low concentration fp using a quasi-ideal solid solution model. For x$_{Fe}$ > 0.12 this approach is no longer valid. At ultra-high pressures, there is an IS to LS spin state change in Fe in the B2 phase, but the transition pressure depends sensitively on thermal electronic excitations and on x$_{Fe}$.


title: Subaru Hyper Suprime-Cam Survey of Cygnus OB2 Complex -- I: Introduction, Photometry and Source Catalog

abstract: Low mass star formation inside massive clusters is crucial to understand the effect of cluster environment on processes like circumstellar disk evolution, planet and brown dwarf formation. The young massive association of Cygnus OB2, with a strong feedback from massive stars, is an ideal target to study the effect of extreme environmental conditions on its extensive low-mass population. We aim to perform deep multi-wavelength studies to understand the role of stellar feedback on the IMF, brown dwarf fraction and circumstellar disk properties in the region. We introduce here, the deepest and widest optical photometry of 1.5$^\circ$ diameter region centred at Cygnus OB2 in r$_{2}$, i$_{2}$, z and Y-filters using Subaru Hyper Suprime-Cam (HSC). This work presents the data reduction, source catalog generation, data quality checks and preliminary results about the pre-main sequence sources. We obtain 713,529 sources in total, with detection down to $\sim$ 28 mag, 27 mag, 25.5 mag and 24.5 mag in r$_{2}$, i$_{2}$, z and Y-band respectively, which is $\sim$ 3 - 5 mag deeper than the existing Pan-STARRS and GTC/OSIRIS photometry. We confirm the presence of a distinct pre-main sequence branch by statistical field subtraction of the central 18$^\prime$ region. We find the median age of the region as $\sim$ 5 $\pm$ 2 Myrs with an average disk fraction of $\sim$ 9$\%$. At this age, combined with A$_V$ $\sim$ 6 - 8 mag, we detect sources down to a mass range $\sim$ 0.01 - 0.17 M$_\odot$. The deep HSC catalog will serve as the groundwork for further studies on this prominent active young cluster.


title: Alzheimers Dementia Detection using Acoustic & Linguistic features and Pre-Trained BERT

abstract: Alzheimers disease is a fatal progressive brain disorder that worsens with time. It is high time we have inexpensive and quick clinical diagnostic techniques for early detection and care. In previous studies, various Machine Learning techniques and Pre-trained Deep Learning models have been used in conjunction with the extraction of various acoustic and linguistic features. Our study focuses on three models for the classification task in the ADReSS (The Alzheimers Dementia Recognition through Spontaneous Speech) 2021 Challenge. We use the well-balanced dataset provided by the ADReSS Challenge for training and validating our models. Model 1 uses various acoustic features from the eGeMAPs feature-set, Model 2 uses various linguistic features that we generated from auto-generated transcripts and Model 3 uses the auto-generated transcripts directly to extract features using a Pre-trained BERT and TF-IDF. These models are described in detail in the models section.


title: SOCIALGYM: A Framework for Benchmarking Social Robot Navigation

abstract: Robots moving safely and in a socially compliant manner in dynamic human environments is an essential benchmark for long-term robot autonomy. However, it is not feasible to learn and benchmark social navigation behaviors entirely in the real world, as learning is data-intensive, and it is challenging to make safety guarantees during training. Therefore, simulation-based benchmarks that provide abstractions for social navigation are required. A framework for these benchmarks would need to support a wide variety of learning approaches, be extensible to the broad range of social navigation scenarios, and abstract away the perception problem to focus on social navigation explicitly. While there have been many proposed solutions, including high fidelity 3D simulators and grid world approximations, no existing solution satisfies all of the aforementioned properties for learning and evaluating social navigation behaviors. In this work, we propose SOCIALGYM, a lightweight 2D simulation environment for robot social navigation designed with extensibility in mind, and a benchmark scenario built on SOCIALGYM. Further, we present benchmark results that compare and contrast human-engineered and model-based learning approaches to a suite of off-the-shelf Learning from Demonstration (LfD) and Reinforcement Learning (RL) approaches applied to social robot navigation. These results demonstrate the data efficiency, task performance, social compliance, and environment transfer capabilities for each of the policies evaluated to provide a solid grounding for future social navigation research.


title: Chaos in self-gravitating many-body systems: Lyapunov time dependence of $N$ and the influence of general relativity

abstract: In self-gravitating $N$-body systems, small perturbations introduced at the start, or infinitesimal errors that are produced by the numerical integrator or are due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results that for relatively homogeneous systems, this rate of growth per crossing time increases with $N$ up to $N \sim 30$, but that for larger systems, the growth rate has a weaker scaling with $N$. For concentrated systems, however, the rate of exponential growth continues to scale with $N$. In relativistic self-gravitating systems, the rate of growth is almost independent of $N$. This effect, however, is only noticeable when the system's mean velocity approaches the speed of light to within three orders of magnitude. The chaotic behavior of systems with more than a dozen bodies for the usually adopted approximation of only solving the pairwise interactions in the Einstein-Infeld-Hoffmann equation of motion is qualitatively different than when the interaction terms (or cross terms) are taken into account. This result provides a strong motivation for follow-up studies on the microscopic effect of general relativity on orbital chaos, and on the influence of higher-order cross-terms in the Taylor-series expansion of the Einstein-Infeld-Hoffmann equations of motion.

categories: nlin.CD astro-ph.CO astro-ph.HE gr-qc physics.comp-ph


title: On a Generalization of Tupper's Formula for $m$ Colours and $n$ Dimensions

abstract: Tupper's formula $\frac{1}{2}<\bigg\lfloor \bmod \bigg(\lfloor \frac{y}{17}\rfloor 2^{-17\lfloor x \rfloor -\bmod (\lfloor y \rfloor,17)},2\bigg)\bigg\rfloor$ has an interesting property that for any monochrome image that can be represented by pixels in a two dimensional array of dimensions $106\times 17$, there exists a natural number $k$ such that the graph of the equation in the range $0\leq x <106$ and $k\leq y<k+17$, is that image. In this paper, we give a generalization for $m$ colours and $n$ dimensions. We give $m$ formulae consisting of $n$ free variables, with the property that, for any $n$ dimensional object of $m$ colours $C_1,\cdots, C_m$, that can be represented by hypervoxels(multidimensional analogue of pixel) in a $n$ dimensional array of dimensions $A_1\times \cdots \times A_n$, there exists a natural number $k$ such that, when the first formula is graphed using colour $C_1$, second formula is graphed using colour $C_2$,$\cdots$, $m$th formula is graphed using colour $C_m$ in the range $0\leq x_1<A_1$,$0\leq x_2<A_2,\cdots, 0\leq x_{n-1}<A_{n-1},k\leq x_n <k+A_n$, the union of all graphs is that $n$-dimensional object.


title: A multivariate approach to single-molecule thermopower and electric conductance measurements

abstract: We report a method using scanning tunnelling microscope single molecular break junction to simultaneously measure and correlate the single-molecule thermopower and electrical conductance. In contrast to previously reported approaches, it does not require custom-built electronics and takes advantage of a trace-by-trace calibration of the thermal offset at the Au/Au contact, thus greatly facilitating thermoelectric measurements at the single-molecule level. We report measurements of three molecules: 1,4-di(4-(ethynyl(phenylthioacetate)) benzene, 1,8-octanedithiol, and 4,4'-bipyridine, and determine single-molecule Seebeck coefficients of 12(3), 5(2), and -5(2) microV K-1, respectively. Furthermore, the method statistically correlates the Seebeck voltage offset, electrical conductance, and stretching displacement of the single-molecule junction, and allows for direct comparison with current-distance spectroscopy results obtained at constant bias.


title: Chiral Dirac Equation and Its Spacetime and CPT Symmetries

abstract: The Dirac equation with chiral symmetry is derived using the irreducible representations of the Poincar\'{e} group, the Lagrangian formalism, and a novel method of projection operators that takes as its starting point the minimal assumption of four linearly independent physical states. We thereby demonstrate the fundamental nature of this form of the Dirac equation. The resulting equation is then examined within the context of spacetime and CPT symmetries with a discussion of the implications for the general formulation of physical theories.


title: Report on the "The Future of the Shell" Panel at HotOS 2021

abstract: This document summarizes the challenges and possible research directions around the shell and its ecosystem, collected during and after the HotOS21 Panel on the future of the shell. The goal is to create a snapshot of what a number of researchers from various disciplines -- connected to the shell to varying degrees -- think about its future. We hope that this document will serve as a reference for future research on the shell and its ecosystem.


title: On the Fate of Interstellar Objects Captured by our Solar System

abstract: With the recent discoveries of interstellar objects `Oumuamua and Borisov traversing the solar system, understanding the dynamics of interstellar objects is more pressing than ever. These detections have highlighted the possibility that captured interstellar material could be trapped in our solar system. The first step in rigorously investigating this question is to calculate a capture cross section for interstellar objects as a function of hyperbolic excess velocity, which can be convolved with any velocity dispersion to compute a capture rate (Napier et. al. 2021). Although the cross section provides the first step toward calculating the mass of alien rocks residing in our solar system, we also need to know the lifetime of captured objects. We use an ensemble of N-body simulations to characterize a dynamical lifetime for captured interstellar objects and determines the fraction of surviving objects as a function of time (since capture). We also illuminate the primary effects driving their secular evolution. Finally, we use the resulting dynamical lifetime function to estimate the current inventory of captured interstellar material in the solar system. We find that capture from the field yields a steady state mass of only $\sim 10^{-13} M_{\oplus}$, whereas the mass remaining from capture events in the birth cluster is roughly $10^{-9} M_{\oplus}$.


title: Making Human-Like Trade-offs in Constrained Environments by Learning from Demonstrations

abstract: Many real-life scenarios require humans to make difficult trade-offs: do we always follow all the traffic rules or do we violate the speed limit in an emergency? These scenarios force us to evaluate the trade-off between collective norms and our own personal objectives. To create effective AI-human teams, we must equip AI agents with a model of how humans make trade-offs in complex, constrained environments. These agents will be able to mirror human behavior or to draw human attention to situations where decision making could be improved. To this end, we propose a novel inverse reinforcement learning (IRL) method for learning implicit hard and soft constraints from demonstrations, enabling agents to quickly adapt to new settings. In addition, learning soft constraints over states, actions, and state features allows agents to transfer this knowledge to new domains that share similar aspects. We then use the constraint learning method to implement a novel system architecture that leverages a cognitive model of human decision making, multi-alternative decision field theory (MDFT), to orchestrate competing objectives. We evaluate the resulting agent on trajectory length, number of violated constraints, and total reward, demonstrating that our agent architecture is both general and achieves strong performance. Thus we are able to capture and replicate human-like trade-offs from demonstrations in environments when constraints are not explicit.


title: Engineering infinite-range SU($n$) interactions with spin-orbit-coupled fermions in an optical lattice

abstract: We study multilevel fermions in an optical lattice described by the Hubbard model with on site SU($n$)-symmetric interactions. We show that in an appropriate parameter regime this system can be mapped onto a spin model with all-to-all SU($n$)-symmetric couplings. Raman pulses that address internal spin states modify the atomic dispersion relation and induce spin-orbit coupling, which can act as a synthetic inhomogeneous magnetic field that competes with the SU($n$) exchange interactions. We investigate the mean-field dynamical phase diagram of the resulting model as a function of $n$ and different initial configurations that are accessible with Raman pulses. Consistent with previous studies for $n=2$, we find that for some initial states the spin model exhibits two distinct dynamical phases that obey simple scaling relations with $n$. Moreover, for $n>2$ we find that dynamical behavior can be highly sensitive to initial intra-spin coherences. Our predictions are readily testable in current experiments with ultracold alkaline-earth(-like) atoms.


title: Exploring Decomposition for Table-based Fact Verification

abstract: Fact verification based on structured data is challenging as it requires models to understand both natural language and symbolic operations performed over tables. Although pre-trained language models have demonstrated a strong capability in verifying simple statements, they struggle with complex statements that involve multiple operations. In this paper, we improve fact verification by decomposing complex statements into simpler subproblems. Leveraging the programs synthesized by a weakly supervised semantic parser, we propose a program-guided approach to constructing a pseudo dataset for decomposition model training. The subproblems, together with their predicted answers, serve as the intermediate evidence to enhance our fact verification model. Experiments show that our proposed approach achieves the new state-of-the-art performance, an 82.7\% accuracy, on the TabFact benchmark.


title: Intel Optane DCPMM and Serverless Computing

abstract: This report describes 1) how we use Intel's Optane DCPMM in the memory Mode. We investigate the the scalability of applications on a single Optane machine, using Subgraph counting as memory-intensive graph problem. We test with various input graph and subtemplate sizes to determine its performance for different memory and CPU loads, as well as a comparison of performance on a single node Optane with a distributed set of nodes in a cluster using MPI. 2) We investigate the end-to-end execution delays in serverless computing and study concurrent function executions with cold starts. In future work, we will show that persistent memory machines may significantly improve concurrent function invocations in serverless computing including Amazon Lambda, Microsoft Azure Functions, Google Cloud Functions and IBM Cloud Functions (Apache OpenWhisk).


title: Galactic Cosmic Ray Acceleration with Steep Spectra

abstract: Galactic cosmic rays (CRs) are accelerated by astrophysical shocks, primarily supernova remnants (SNRs), via diffusive shock acceleration (DSA), an efficient mechanism that predicts power-law energy distributions of CRs. However, observations of both nonthermal SNR emission and Galactic CRs imply CR spectra that are steeper than the standard DSA prediction, $\propto E^{-2}$. Recent kinetic hybrid simulations suggest that such steep spectra may be the result of a "postcursor", or drift of CRs and magnetic structures with respect to the thermal plasma behind the shock. Using a semi-analytic model of non-linear DSA, we generalize this result to a wide range of astrophysical shocks. By accounting for the presence of a postcursor, we produce CR energy distributions that are substantially steeper than $E^{-2}$ and consistent with observations. Our formalism reproduces both modestly steep spectra of Galactic SNRs ($\propto E^{-2.2}$) and the very steep spectra of young radio supernovae ($\propto E^{-3}$).


title: An Investigation And Insight Into Terrorism In Nigeria

abstract: Terrorism is one of the most life-challenging threats facing humanity worldwide. The activities of terrorist organizations threaten peace, disrupts progress, and halt the development of any nation. Terrorist activities in Nigeria in the last decades have negatively affected economic growth and have drastically reduced the possibilities of foreign investments in Nigeria. In this paper, statistical and inferential insights are applied to the terrorist activities in Nigeria from 1970 to 2019. Using the Global Terrorism Database (GTD), insights are made on the occurrences of terrorist attacks, the localities of the target, and the successful and unsuccessful rates of such attacks. The Apriori algorithm is also used in this paper to draw hidden patterns from the GTD to aid in generating strong rules through database mining, resulting in relevant insights. This understanding of terrorist activities will provide security agencies with the needed information to be one step ahead of terrorists in making the right decisions targeted at curbing terrorism in Nigeria.


title: Social-Media Activity Forecasting with Exogenous Information Signals

abstract: Due to their widespread adoption, social media platforms present an ideal environment for studying and understanding social behavior, especially on information spread. Modeling social media activity has numerous practical implications such as supporting efforts to analyze strategic information operations, designing intervention techniques to mitigate disinformation, or delivering critical information during disaster relief operations. In this paper we propose a modeling technique that forecasts topic-specific daily volume of social media activities by using both exogenous signals, such as news or armed conflicts records, and endogenous data from the social media platform we model. Empirical evaluations with real datasets from two different platforms and two different contexts each composed of multiple interrelated topics demonstrate the effectiveness of our solution.


title: Spaceability on some classes of Banach spaces

abstract: In this paper, we study spaceability of subsets of generalized Orlicz and Lebesgue spaces associated to Banach function space. Also, we give some sufficient conditions for spaceability of subsets of a general Banach space which improves an important result on this topic. As an application, it is shown that the set of all bounded linear operators which are not positive semidefinite on a separable Hilbert space is spaceable.


title: Outbursts and stellar properties of the classical Be star HD 6226

abstract: The bright and understudied classical Be star HD 6226 has exhibited multiple outbursts in the last several years during which the star grew a viscous decretion disk. We analyze 659 optical spectra of the system collected from 2017-2020, along with a UV spectrum from the Hubble Space Telescope and high cadence photometry from both TESS and the KELT survey. We find that the star has a spectral type of B2.5IIIe, with a rotation rate of 74% of critical. The star is nearly pole-on with an inclination of $13.4$ degree. We confirm the spectroscopic pulsational properties previously reported, and report on three photometric oscillations from KELT photometry. The outbursting behavior is studied with equivalent width measurements of H$\alpha$ and H$\beta$, and the variations in both of these can be quantitatively explained with two frequencies through a Fourier analysis. One of the frequencies for the emission outbursts is equal to the difference between two photometric oscillations, linking these pulsation modes to the mass ejection mechanism for some outbursts. During the TESS observation time period of 2019 October 7 to 2019 November 2, the star was building a disk. With a large dataset of H$\alpha$ and H$\beta$ spectroscopy, we are able to determine the timescales of dissipation in both of these lines, similar to past work on Be stars that has been done with optical photometry. HD 6226 is an ideal target with which to study the Be disk-evolution given its apparent periodic nature, allowing for targeted observations with other facilities in the future.


title: Quantile-based fuzzy C-means clustering of multivariate time series: Robust techniques

abstract: Three robust methods for clustering multivariate time series from the point of view of generating processes are proposed. The procedures are robust versions of a fuzzy C-means model based on: (i) estimates of the quantile cross-spectral density and (ii) the classical principal component analysis. Robustness to the presence of outliers is achieved by using the so-called metric, noise and trimmed approaches. The metric approach incorporates in the objective function a distance measure aimed at neutralizing the effect of the outliers, the noise approach builds an artificial cluster expected to contain the outlying series and the trimmed approach eliminates the most atypical series in the dataset. All the proposed techniques inherit the nice properties of the quantile cross-spectral density, as being able to uncover general types of dependence. Results from a broad simulation study including multivariate linear, nonlinear and GARCH processes indicate that the algorithms are substantially effective in coping with the presence of outlying series (i.e., series exhibiting a dependence structure different from that of the majority), clearly poutperforming alternative procedures. The usefulness of the suggested methods is highlighted by means of two specific applications regarding financial and environmental series.


title: On physics-informed data-driven isotropic and anisotropic constitutive models through probabilistic machine learning and space-filling sampling

abstract: Data-driven constitutive modeling is an emerging field in computational solid mechanics with the prospect of significantly relieving the computational costs of hierarchical computational methods. Traditionally, these surrogates have been trained using datasets which map strain inputs to stress outputs directly. Data-driven constitutive models for elastic and inelastic materials have commonly been developed based on artificial neural networks (ANNs), which recently enabled the incorporation of physical laws in the construction of these models. However, ANNs do not offer convergence guarantees and are reliant on user-specified parameters. In contrast to ANNs, Gaussian process regression (GPR) is based on nonparametric modeling principles as well as on fundamental statistical knowledge and hence allows for strict convergence guarantees. GPR however has the major disadvantage that it scales poorly as datasets get large. In this work we present a physics-informed data-driven constitutive modeling approach for isostropic and anisotropic materials based on probabilistic machine learning that can be used in the big data context. The trained GPR surrogates are able to respect physical principles such as material frame indifference, material symmetry, thermodynamic consistency, stress-free undeformed configuration, and the local balance of angular momentum. Furthermore, this paper presents the first sampling approach that directly generates space-filling points in the invariant space corresponding to bounded domain of the gradient deformation tensor. Overall, the presented approach is tested on synthetic data from isotropic and anisotropic constitutive laws and shows surprising accuracy even far beyond the limits of the training domain, indicating that the resulting surrogates can efficiently generalize as they incorporate knowledge about the underlying physics.


title: From Steklov to Laplace: free boundary minimal surfaces with many boundary components

abstract: In the present paper, we study sharp isoperimetric inequalities for the first Steklov eigenvalue $\sigma_1$ on surfaces with fixed genus and large number $k$ of boundary components. We show that as $k\to \infty$ the free boundary minimal surfaces in the unit ball arising from the maximization of $\sigma_1$ converge to a closed minimal surface in the boundary sphere arising from the maximization of the first Laplace eigenvalue on the corresponding closed surface. For some genera, we prove that the corresponding areas converge at the optimal rate $\frac{\log k}{k}$. This result appears to provide the first examples of free boundary minimal surfaces in a compact domain converging to closed minimal surfaces in the boundary, suggesting new directions in the study of free boundary minimal surfaces, with many open questions proposed in the present paper. A similar phenomenon is observed for free boundary harmonic maps associated to conformally-constrained shape optimization problems.


title: A Tunnel Diode Oscillator for High Sensitivity Broad Temperature Range Susceptibility Measurements

abstract: We report the design and operation of a versatile tunnel diode oscillator for high sensitivity (50 ppb) measurements of the magnetic susceptibility of samples over a wide temperature range (1.7-100 K) that can be used with a simple liquid helium storage dewar. The design allows for the application of an electric field across the sample to search for magnetoelectric effects.


title: Using the Carnot cycle to determine changes of the phase transition temperature

abstract: The Clausius-Clapeyron relation and its analogs in other first-order phase transitions, such as type-I superconductors, are derived using very elementary methods, without appealing to the more advanced concepts of entropy or Gibbs free energy. The reasoning is based on Kelvin's formulation of the second law of thermodynamics, and should be accessible to high school students. After recalling some basic facts about the Carnot cycle, we present two very different systems that undergo discontinuous phase transitions (ice/water and normal/superconductor), and construct engines that exploit the properties of these systems to produce work. In each case, we show that if the transition temperature $T_tr$ were independent of other parameters, such as pressure or magnetic field, it would be possible to violate Kelvin's principle, i.e., to construct a perpetuum mobile of the second kind. Since the proposed cyclic processes can be realized reversibly in the limit of infinitesimal changes in temperature, their efficiencies must be equal to that of an ordinary Carnot cycle. We immediately obtain an equation of the form $dT /dX = f(T, X)$, which governs how the transition temperature changes with the parameter $X$.


title: Searching for New Physics in Rare $K$ and $B$ Decays without $|V_{cb}|$ and $|V_{ub}|$ Uncertainties

abstract: We reemphasize the strong dependence of the branching ratios $B(K^+\to\pi^+\nu\bar\nu)$ and $B(K_L\to\pi^0\nu\bar\nu)$ on $|V_{cb}|$ that is stronger than in rare $B$ decays, in particular for $K_L\to\pi^0\nu\bar\nu$. Thereby the persistent tension between inclusive and exclusive determinations of $|V_{cb}|$ weakens the power of these theoretically clean decays in the search for new physics (NP). We demonstrate how this uncertainty can be practically removed by considering within the SM suitable ratios of the two branching ratios between each other and with other observables like the branching ratios for $K_S\to\mu^+\mu^-$, $B_{s,d}\to\mu^+\mu^-$ and $B\to K(K^*)\nu\bar\nu$. We use as basic CKM parameters $V_{us}$, $|V_{cb}|$ and the angles $\beta$ and $\gamma$ in the unitarity triangle (UT). This avoids the use of the problematic $|V_{ub}|$. A ratio involving $B(K^+\to\pi^+\nu\bar\nu)$ and $B(B_s\to\mu^+\mu^-)$ while being $|V_{cb}|$-independent exhibits sizable dependence on the angle $\gamma$. It should be of interest for several experimental groups in the coming years. We point out that the $|V_{cb}|$-independent ratio of $B(B^+\to K^+\nu\bar\nu)$ and $B(B_s\to\mu^+\mu^-)$ from Belle II and LHCb signals a $1.8\sigma$ tension with its SM value. As a complementary test of the Standard Model, we propose to extract $|V_{cb}|$ from different observables as a function of $\beta$ and $\gamma$. We illustrate this with $\epsilon_K$, $\Delta M_d$ and $\Delta M_s$ finding tensions between these three determinations of $|V_{cb}|$ within the SM. From $\Delta M_s$ and $S_{\psi K_S}$ alone we find $|V_{cb}|=41.8(6)\times 10^{-3}$ and $|V_{ub}|=3.65(12)\times 10^{-3}$. We stress the importance of a precise measurement of $\gamma$. We obtain most precise SM predictions for considered branching ratios of rare K and B decays to date.


title: Collection of Propagating Electromagnetic Fields by Uncoated Probe

abstract: Understanding light-matter interaction at the nanoscale by observation of fine details of electromagnetic fields is achieved by bringing nanoscale probes into the nearfield of light sources, capturing information that is lost in the far field. Although metal coated probes are often used for nearfield microscopy, they strongly perturb the electromagnetic fields under study. Here, through experiment and simulation, we detail light collection by uncoated fiber probes, which minimize such perturbation. Second-harmonic light from intensely-irradiated sub-wavelength sub-surface features was imaged to avoid otherwise dominating fundamental light background, yielding clear nearfield details through a 50 nm aperture uncoated probe with ~23 nm optical resolution. Simulations shows how a metallic coating distorts optical nearfields and limits optical coupling into the probe in comparison to an uncoated probe.


title: Conditional Poisson Stochastic Beam Search

abstract: Beam search is the default decoding strategy for many sequence generation tasks in NLP. The set of approximate K-best items returned by the algorithm is a useful summary of the distribution for many applications; however, the candidates typically exhibit high overlap and may give a highly biased estimate for expectations under our model. These problems can be addressed by instead using stochastic decoding strategies. In this work, we propose a new method for turning beam search into a stochastic process: Conditional Poisson stochastic beam search. Rather than taking the maximizing set at each iteration, we sample K candidates without replacement according to the conditional Poisson sampling design. We view this as a more natural alternative to Kool et. al. 2019's stochastic beam search (SBS). Furthermore, we show how samples generated under the CPSBS design can be used to build consistent estimators and sample diverse sets from sequence models. In our experiments, we observe CPSBS produces lower variance and more efficient estimators than SBS, even showing improvements in high entropy settings.


title: Continuity Properties of Value Functions in Information Structures for Zero-Sum and General Games and Stochastic Teams

abstract: We study continuity properties of stochastic game problems with respect to various topologies on information structures, defined as probability measures characterizing a game. We will establish continuity properties of the value function under total variation, setwise, and weak convergence of information structures. Our analysis reveals that the value function for a bounded game is continuous under total variation convergence of information structures in both zero-sum games and team problems. Continuity may fail to hold under setwise or weak convergence of information structures, however, the value function exhibits upper semicontinuity properties under weak and setwise convergence of information structures for team problems, and upper or lower semicontinuity properties hold for zero-sum games when such convergence is through a Blackwell-garbled sequence of information structures. If the individual channels are independent, fixed, and satisfy a total variation continuity condition, then the value functions are continuous under weak convergence of priors. We finally show that value functions for players may not be continuous even under total variation convergence of information structures in general non-zero-sum games.


title: Matrix product states with large sites

abstract: We explore various ways to group orbitals into clusters in a matrix product state (MPS). We explain how a generic cluster MPS can often lead to an increase in computational cost and instead propose a special cluster structure, involving only the first and last orbitals/sites, with a wider scope for computational advantage. This structure is a natural formalism to describe correlated multireference (MR) theories. We demonstrate the flexibility and usefulness of this approach by implementing various uncontracted MR configuration interaction, perturbation and linearized coupled cluster theories using an MPS with large cluster sites. Applications to the nitrogen dimer, the chromium dimer, and benzene, including up to triple excitations in the external space, demonstrate the utility of an MPS with up to two large sites. We use our results to analyze the quality of different multireference approximations.

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


title: A Computational Approach for Checking Compliance with European View and Sunlight Exposure Criteria

abstract: The paper presents open-source computational workflows for assessing the "Exposure to sunlight" and "View out" criteria as defined in the European standard EN 17037 "Daylight in Buildings", issued by the European Committee for Standardization. In addition to these factors, the standard document also addresses daylight provision and protection from glare, both of which fall out of the scope of this paper. The purpose of the standard is stated as 'encouraging building designers to assess and ensure successfully daylit spaces'. The standard document proposes verification methods for performing such assessments, albeit without recommending a simulation procedure for computing the aforementioned criteria. The workflows proposed in this paper are arguably the first attempt to standardize these assessment methods using de-facto open-source standard technologies currently used in practice. The approach of this work is twofold: establish that the compliance check can be systematically performed on a 3D model by a novel simulation tool developed by the authors; and highlighting the additional assumptions that need to be implemented to build a robust and unambiguous tool within existing open-source frameworks.


title: Static analysis for coupled nonlinear Klein-Gordon equations with asymmetric parameter settings

abstract: Klein-Gordon equations describe the dynamics of waves/particles in sub-atomic scales. For a system of nonlinear Klein-Gordon equations, a systematic analysis of the time evolution for their spatially uniform solutions has been performed \cite{21takei}. In the study, the parameters (mass, wave propagation speed, and the force parameters) are chosen to be symmetric between the two single equations. Symmetric parameter settings are equivalent to assume the interacting two same particles. In this paper, for a system of nonlinear Klein-Gordon equations with asymmetric parameter settings, the time evolution for their spatially uniform solutions are studied. This is equivalent to assume the interacting two different particles. As a result, based on the high precision numerical scheme \cite{22takei}, the existence of divergent and bounded solutions that depend on parameter settings is revealed. The competition, coherence, and decoherence of different waves are shown to appear depending on the choice of asymmetrically-implemented parameter values.


title: Classification of (2+1)D invertible fermionic topological phases with symmetry

abstract: We provide a classification of invertible topological phases of interacting fermions with symmetry in two spatial dimensions for general fermionic symmetry groups $G_f$ and general values of the chiral central charge $c_-$. Here $G_f$ is a central extension of a bosonic symmetry group $G_b$ by fermion parity, $(-1)^F$, specified by a second cohomology class $[\omega_2] \in \mathcal{H}^2(G_b, \mathbb{Z}_2)$. Our approach proceeds by gauging fermion parity and classifying the resulting $G_b$ symmetry-enriched topological orders while keeping track of certain additional data and constraints. We perform this analysis through two perspectives, using $G$-crossed braided tensor categories and Spin$(2c_-)_1$ Chern-Simons theory coupled to a background $G$ gauge field. These results give a way to characterize and classify invertible fermionic topological phases in terms of a concrete set of data and consistency equations, which is more physically transparent and computationally simpler than the more abstract methods using cobordism theory and spectral sequences. Our results also generalize and provide a different approach to the recent classification of fermionic symmetry-protected topological phases by Wang and Gu, which have chiral central charge $c_- = 0$. We show how the 10-fold way classification of topological insulators and superconductors fits into our scheme, along with general non-perturbative constraints due to certain choices of $c_-$ and $G_f$. Mathematically, our results also suggest an explicit general parameterization of deformation classes of (2+1)D invertible topological quantum field theories with $G_f$ symmetry.

categories: cond-mat.str-el hep-th math-ph math.MP math.QA quant-ph


title: Can one predict a drop contact angle?

abstract: The study of wetting phenomena is of great interest due to the multifaceted technological applications of hydrophobic and hydrophilic surfaces. The theoretical approaches proposed by Wenzel and later by Cassie and Baxter to describe the behaviour of a droplet of water on a rough solid were extensively used and improved to characterize the apparent contact angle of a droplet. However, the equilibrium hypothesis implied in these models means that they are not always predictive of experimental contact angles due to strong metastabilities typically occurring on heterogeneous surfaces. A predictive scheme for contact angle is thus urgently needed both to characterise a surface by contact angle measurements and to design superhydrophobic and oleophobic surfaces with the desired properties, e.g., contact angle hysteresis. In this work a combination of Monte Carlo simulation and the string method is employed to calculate the free energy profile of a liquid droplet deposited on a pillared surface. For the analyzed surfaces, we show that there is only one minimum of the free energy that corresponds to the superhydrophobic wetting state while the wet state can present multiple minima. Furthermore, when the surface roughness decreases the amount of local minima observed in the free energy profile increases.


title: Security Analysis of Capsule Network Inference using Horizontal Collaboration

abstract: The traditional convolution neural networks (CNN) have several drawbacks like the Picasso effect and the loss of information by the pooling layer. The Capsule network (CapsNet) was proposed to address these challenges because its architecture can encode and preserve the spatial orientation of input images. Similar to traditional CNNs, CapsNet is also vulnerable to several malicious attacks, as studied by several researchers in the literature. However, most of these studies focus on single-device-based inference, but horizontally collaborative inference in state-of-the-art systems, like intelligent edge services in self-driving cars, voice controllable systems, and drones, nullify most of these analyses. Horizontal collaboration implies partitioning the trained CNN models or CNN tasks to multiple end devices or edge nodes. Therefore, it is imperative to examine the robustness of the CapsNet against malicious attacks when deployed in horizontally collaborative environments. Towards this, we examine the robustness of the CapsNet when subjected to noise-based inference attacks in a horizontal collaborative environment. In this analysis, we perturbed the feature maps of the different layers of four DNN models, i.e., CapsNet, Mini-VGG, LeNet, and an in-house designed CNN (ConvNet) with the same number of parameters as CapsNet, using two types of noised-based attacks, i.e., Gaussian Noise Attack and FGSM noise attack. The experimental results show that similar to the traditional CNNs, depending upon the access of the attacker to the DNN layer, the classification accuracy of the CapsNet drops significantly. For example, when Gaussian Noise Attack classification is performed at the DigitCap layer of the CapsNet, the maximum classification accuracy drop is approximately 97%.


title: Two-loop application of the Breitenlohner-Maison/'t Hooft-Veltman scheme with non-anticommuting $\gamma_5$: Full renormalization and symmetry-restoring counterterms in an abelian chiral gauge theory

abstract: We apply the BMHV scheme for non-anticommuting $\gamma_5$ to an abelian chiral gauge theory at the two-loop level. As our main result, we determine the full structure of symmetry-restoring counterterms up to the two-loop level. These counterterms turn out to have the same structure as at the one-loop level and a simple interpretation in terms of restoration of well-known Ward identities. In addition, we show that the ultraviolet divergences cannot be canceled completely by counterterms generated by field and parameter renormalization, and we determine needed UV divergent evanescent counterterms. The paper establishes the two-loop methodology based on the quantum action principle and direct computations of Slavnov-Taylor identity breakings. The same method will be applicable to nonabelian gauge theories.


title: Learning Predictive and Interpretable Timeseries Summaries from ICU Data

abstract: Machine learning models that utilize patient data across time (rather than just the most recent measurements) have increased performance for many risk stratification tasks in the intensive care unit. However, many of these models and their learned representations are complex and therefore difficult for clinicians to interpret, creating challenges for validation. Our work proposes a new procedure to learn summaries of clinical time-series that are both predictive and easily understood by humans. Specifically, our summaries consist of simple and intuitive functions of clinical data (e.g. falling mean arterial pressure). Our learned summaries outperform traditional interpretable model classes and achieve performance comparable to state-of-the-art deep learning models on an in-hospital mortality classification task.


title: Adapting conditional simulation using circulant embedding for irregularly spaced spatial data

abstract: Computing an ensemble of random fields using conditional simulation is an ideal method for retrieving accurate estimates of a field conditioned on available data and for quantifying the uncertainty of these realizations. Methods for generating random realizations, however, are computationally demanding, especially when the estimates are conditioned on numerous observed data and for large domains. In this article, a \textit{new}, \textit{approximate} conditional simulation approach is applied that builds on \textit{circulant embedding} (CE), a fast method for simulating stationary Gaussian processes. The standard CE is restricted to simulating stationary Gaussian processes (possibly anisotropic) on regularly spaced grids. In this work we explore two possible algorithms, namely local Kriging and nearest neighbor Kriging, that extend CE for irregularly spaced data points. We establish the accuracy of these methods to be suitable for practical inference and the speedup in computation allows for generating conditional fields close to an interactive time frame. The methods are motivated by the U.S. Geological Survey's software \textit{ShakeMap}, which provides near real-time maps of shaking intensity after the occurrence of a significant earthquake. An example for the 2019 event in Ridgecrest, California is used to illustrate our method.


title: Training Deep Spiking Auto-encoders without Bursting or Dying Neurons through Regularization

abstract: Spiking neural networks are a promising approach towards next-generation models of the brain in computational neuroscience. Moreover, compared to classic artificial neural networks, they could serve as an energy-efficient deployment of AI by enabling fast computation in specialized neuromorphic hardware. However, training deep spiking neural networks, especially in an unsupervised manner, is challenging and the performance of a spiking model is significantly hindered by dead or bursting neurons. Here, we apply end-to-end learning with membrane potential-based backpropagation to a spiking convolutional auto-encoder with multiple trainable layers of leaky integrate-and-fire neurons. We propose bio-inspired regularization methods to control the spike density in latent representations. In the experiments, we show that applying regularization on membrane potential and spiking output successfully avoids both dead and bursting neurons and significantly decreases the reconstruction error of the spiking auto-encoder. Training regularized networks on the MNIST dataset yields image reconstruction quality comparable to non-spiking baseline models (deterministic and variational auto-encoder) and indicates improvement upon earlier approaches. Importantly, we show that, unlike the variational auto-encoder, the spiking latent representations display structure associated with the image class.


title: DARSI: An Ontology for Facilitating the Development of Data Sharing and Use Agreements

abstract: The advantages of data sharing across organizations and disciplines are indisputable; although, sensitive and restricted data cannot be easily shared due to policies and legal matters. The research presented in this paper takes a step toward systematizing the sharing of sensitive and restricted research data by developing an ontology to frame and guide DSUA (Data Sharing and Usage Agreement) development. The paper provides background context, describes the ontology creation process, and introduces the Data Sharing Agreements for Restricted and Sensitive Information (DARSI) ontology. DARSI contains four top level classes, 20 sub-classes, 33 sub-categories, and 17 simple properties for categories applicable at various levels. The discussion provides further insight into the work accomplished, and the conclusion identifies next steps.


title: Cross-Modal Coherence for Text-to-Image Retrieval

abstract: Common image-text joint understanding techniques presume that images and the associated text can universally be characterized by a single implicit model. However, co-occurring images and text can be related in qualitatively different ways, and explicitly modeling it could improve the performance of current joint understanding models. In this paper, we train a Cross-Modal Coherence Modelfor text-to-image retrieval task. Our analysis shows that models trained with image--text coherence relations can retrieve images originally paired with target text more often than coherence-agnostic models. We also show via human evaluation that images retrieved by the proposed coherence-aware model are preferred over a coherence-agnostic baseline by a huge margin. Our findings provide insights into the ways that different modalities communicate and the role of coherence relations in capturing commonsense inferences in text and imagery.


title: Towards practical object detection for weed spraying in precision agriculture

abstract: The evolution of smaller, faster processors and cheaper digital storage mechanisms across the last 4-5 decades has vastly increased the opportunity to integrate intelligent technologies in a wide range of practical environments to address a broad spectrum of tasks. One exciting application domain for such technologies is precision agriculture, where the ability to integrate on-board machine vision with data-driven actuation means that farmers can make decisions about crop care and harvesting at the level of the individual plant rather than the whole field. This makes sense both economically and environmentally. However, the key driver for this capability is fast and robust machine vision -- typically driven by machine learning (ML) solutions and dependent on accurate modelling. One critical challenge is that the bulk of ML-based vision research considers only metrics that evaluate the accuracy of object detection and do not assess practical factors. This paper introduces three metrics that highlight different aspects relevant for real-world deployment of precision weeding and demonstrates their utility through experimental results.


title: Simulating the neutrino flux from the Spallation Neutron Source for the COHERENT experiment

abstract: The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is a pulsed source of neutrons and, as a byproduct of this operation, an intense source of pulsed neutrinos via stopped-pion decay. The COHERENT collaboration uses this source to investigate coherent elastic neutrino-nucleus scattering and other physics with a suite of detectors. This work includes a description of our Geant4 simulation of neutrino production at the SNS and the flux calculation which informs the COHERENT studies. We estimate the uncertainty of this calculation at about 10% based on validation against available low-energy pion production data.


title: Atlas of Urban Scaling Laws

abstract: Highly accurate estimates of the urban fractal dimension $D_f$ are obtained by implementing the Detrended Moving Average algorithm (DMA) on WorldView2 satellite high-resolution multi-spectral images covering the largest European cities. Higher fractal dimensions are systematically obtained for urban sectors (centrally located areas) than for suburban and peripheral areas, with $Df$ values ranging from $1.65$ to $1.90$ respectively.The exponents $\beta_s$ and $\beta_i$ of the scaling law $N^{\beta}$ with $N$ the population size, respectively for socio-economic and infrastructural variables, are evaluated for different urban and suburban sectors in terms of the fractal dimension $D_f$.Results confirm the range of empirical values reported in the literature. Urban scaling laws have been traditionally derived as if cities were zero-dimensional objects with the relevant feature related to a single homogeneous population value, thus neglecting the microscopic heterogeneity of the urban structure. Our findings allow one to go beyond this limit. High sensitive and repeatable satellite records yield robust local estimates of the Hurst and scaling exponents. Furthermore, the approach allows one to discriminate among different scaling theories, shedding light on the open issue of scaling phenomena, reconciling contradictory scientific perspectives and pave paths to the systematic adoption of the complex system science approach to urban landscape analysis.


title: Estimating the number of serial killers that were never caught

abstract: Many serial killers commit tens of murders. At the same time inter-murder intervals can be decades long. This suggests that some serial killers can die of an accident or a disease, having been never caught. We use the distribution of the killers by the number of murders, the distribution of the length of inter-murder intervals and USA life tables to estimate the number of the uncaught killers. The result is that in 20th century there were about seven of such killers. The most prolific of them likely committed over sixty murders.


title: On Bonus-Based Exploration Methods in the Arcade Learning Environment

abstract: Research on exploration in reinforcement learning, as applied to Atari 2600 game-playing, has emphasized tackling difficult exploration problems such as Montezuma's Revenge (Bellemare et al., 2016). Recently, bonus-based exploration methods, which explore by augmenting the environment reward, have reached above-human average performance on such domains. In this paper we reassess popular bonus-based exploration methods within a common evaluation framework. We combine Rainbow (Hessel et al., 2018) with different exploration bonuses and evaluate its performance on Montezuma's Revenge, Bellemare et al.'s set of hard of exploration games with sparse rewards, and the whole Atari 2600 suite. We find that while exploration bonuses lead to higher score on Montezuma's Revenge they do not provide meaningful gains over the simpler $\epsilon$-greedy scheme. In fact, we find that methods that perform best on that game often underperform $\epsilon$-greedy on easy exploration Atari 2600 games. We find that our conclusions remain valid even when hyperparameters are tuned for these easy-exploration games. Finally, we find that none of the methods surveyed benefit from additional training samples (1 billion frames, versus Rainbow's 200 million) on Bellemare et al.'s hard exploration games. Our results suggest that recent gains in Montezuma's Revenge may be better attributed to architecture change, rather than better exploration schemes; and that the real pace of progress in exploration research for Atari 2600 games may have been obfuscated by good results on a single domain.


title: Automated Feature-Topic Pairing: Aligning Semantic and Embedding Spaces in Spatial Representation Learning

abstract: Automated characterization of spatial data is a kind of critical geographical intelligence. As an emerging technique for characterization, Spatial Representation Learning (SRL) uses deep neural networks (DNNs) to learn non-linear embedded features of spatial data for characterization. However, SRL extracts features by internal layers of DNNs, and thus suffers from lacking semantic labels. Texts of spatial entities, on the other hand, provide semantic understanding of latent feature labels, but is insensible to deep SRL models. How can we teach a SRL model to discover appropriate topic labels in texts and pair learned features with the labels? This paper formulates a new problem: feature-topic pairing, and proposes a novel Particle Swarm Optimization (PSO) based deep learning framework. Specifically, we formulate the feature-topic pairing problem into an automated alignment task between 1) a latent embedding feature space and 2) a textual semantic topic space. We decompose the alignment of the two spaces into: 1) point-wise alignment, denoting the correlation between a topic distribution and an embedding vector; 2) pair-wise alignment, denoting the consistency between a feature-feature similarity matrix and a topic-topic similarity matrix. We design a PSO based solver to simultaneously select an optimal set of topics and learn corresponding features based on the selected topics. We develop a closed loop algorithm to iterate between 1) minimizing losses of representation reconstruction and feature-topic alignment and 2) searching the best topics. Finally, we present extensive experiments to demonstrate the enhanced performance of our method.


title: Point vortex dynamics in three-dimensional ageostrophic balanced flows

abstract: Geophysical turbulent flows, characterized by rapid rotation, quantified by small Rossby number, and stable stratification, often self-organize into a collection of coherent vortices, referred to as a vortex gas. The lowest order asymptotic expansion in Rossby number is quasigeostrophy which has purely horizontal velocities and cyclone-anticylone antisymmetry. Ageostrophic effects are important components of many geophysical flows and, as such, these phenomena are not well-modeled by quasigeostrophy. The next order correction in Rossby number, which includes ageostrohpic affects, is so-called balanced dynamics. Balanced dynamics includes ageostrophic vertical velocity and breaks the geostrophic cyclone-anticylone antisymmetry. Point vortex solutions are well known in 2d and quasigeostrophic dynamics and are useful for studying the vortex gas regime of geophysical turbulence. Here we find point vortex solutions in fully three dimensional continuously stratified $QG^{+1}$ dynamics, a particular formulation of balanced dynamics. Simulations of $QG^{+1}$ point vortices show several interesting features not captured by quasigeostrophic point vortices including significant vertical transport on long timescales. The ageostrophic component of $QG^{+1}$ point vortex dynamics renders them useful in modeling flows where quasigeostrophy filters out important physical processes.


title: Optically- and thermally-driven huge lattice orbital and spin angular momenta from spinning fullerenes

abstract: Lattice vibration in solids may carry angular momentum. But unlike the intrinsic spin of electrons, the lattice vibration is rarely rotational. To induce angular momentum, one needs to find a material that can accommodate a twisted normal mode, two orthogonal modes or excitation of magnons. If excitation is too strong, one may exceed the Lindemann limit, so the material melts. Therefore these methods are not ideal. Here, we theoretically propose a new route to phonon angular momentum in a molecular crystal $\rm C_{60}$. We find that a single laser pulse is able to inject a significant amount of angular momentum to $\rm C_{60}$, and the momentum transfer is helicity-dependent. Changing from right-circularly polarized light to left-circularly polarized light switches the direction of phonon angular momentum. On the ultrafast time scale, the orbital angular momentum change closely resembles the displacive excitation of coherent phonons, with a cosine-function dependence on time, different from the spin counterpart. Atomic displacements, even under strong laser excitation, remain far below the Lindemann criterion. Under thermal excitation, spinning $\rm C_{60}$ even at room temperature generates a huge angular momentum close to several hundred $\hbar$. Our finding opens the door to a large group of fullerenes, from $\rm C_{60}$, C$_{70}$ to their endohedral derivatives, where angular momentum can be generated through light or temperature. This paves the way to the phononic control electronic spin and harvesting thermal energy through phonon angular momentum.


title: Kinematics And Ellipsoidal Motion Of The Mid To Late M-Type Stars

abstract: Truly, kinematics in the Solar neighborhood has provided important information both for the structure and for the evolution of the Galaxy since the early 20th century. The relationship between Oort constants and the ratio of the velocity dispersion are important quantities in stellar kinematics. In the present paper, we calculated the kinematical parameters and the Oort constants of various samples of late to intermediate M-type stars. We calculated the velocity dispersion ({\sigma}_1,{\sigma}_2,{\sigma}_3 ) in units of km s-1 for the samples under study. The longitude of the vertex (l_2 ) having negative values with our analysis; i.e., the program I (538 stars), l_2 =-0_.^o 5410, program II (100 stars), l_2=-0_.^o 4937 and program III (60 stars), l_2=-0_.^o 9495. We calculate the Oort constants as A=14.69+-0.61 km s-1 kpc-1 and B=-16.70+-0.67 km s-1 kpc-1, and the rotational velocity V_o=257.38+-9.40 km s^(-1). A possible explanation for the overestimated values of the second Oort constants has been presented.


title: Weighted Low Rank Matrix Approximation and Acceleration

abstract: Low-rank matrix approximation is one of the central concepts in machine learning, with applications in dimension reduction, de-noising, multivariate statistical methodology, and many more. A recent extension to LRMA is called low-rank matrix completion (LRMC). It solves the LRMA problem when some observations are missing and is especially useful for recommender systems. In this paper, we consider an element-wise weighted generalization of LRMA. The resulting weighted low-rank matrix approximation technique therefore covers LRMC as a special case with binary weights. WLRMA has many applications. For example, it is an essential component of GLM optimization algorithms, where an exponential family is used to model the entries of a matrix, and the matrix of natural parameters admits a low-rank structure. We propose an algorithm for solving the weighted problem, as well as two acceleration techniques. Further, we develop a non-SVD modification of the proposed algorithm that is able to handle extremely high-dimensional data. We compare the performance of all the methods on a small simulation example as well as a real-data application.