title: Topological quantum walks in momentum space with a Bose-Einstein condensate

abstract: We report the experimental implementation of discrete-time topological quantum walks of a Bose-Einstein condensate in momentum space. Introducing stroboscopic driving sequences to the generation of a momentum lattice, we show that the dynamics of atoms along the lattice is effectively governed by a periodically driven Su-Schrieffer-Heeger model, which is equivalent to a discrete-time topological quantum walk. We directly measure the underlying topological invariants through time-averaged mean chiral displacements, which are consistent with our experimental observation of topological phase transitions. We then observe interaction-induced localization in the quantum-walk dynamics, where atoms tend to populate a single momentum-lattice site under interactions that are non-local in momentum space. Our experiment opens up the avenue of investigating discrete-time topological quantum walks using cold atoms, where the many-body environment and tunable interactions offer exciting new possibilities.


title: Binary optimal linear codes from posets of the disjoint union of two chains

abstract: Recently, Chang and Hyun obtained some classes of binary optimal codes via simplicial complexes. In this letter, we utilize posets of the disjoint union of two chains to construct binary optimal linear codes.


title: Pruned Landmark Labeling Meets Vertex Centric Computation: A Surprisingly Happy Marriage!

abstract: In this paper, we study how the Pruned Landmark Labeling (PPL) algorithm can be parallelized in a scalable fashion, producing the same results as the sequential algorithm. More specifically, we parallelize using a Vertex-Centric (VC) computational model on a modern SIMD powered multicore architecture. We design a new VC-PLL algorithm that resolves the apparent mismatch between the inherent sequential dependence of the PLL algorithm and the Vertex- Centric (VC) computing model. Furthermore, we introduce a novel batch execution model for VC computation and the BVC-PLL algorithm to reduce the computational inefficiency in VC-PLL. Quite surprisingly, the theoretical analysis reveals that under a reasonable assumption, BVC-PLL has lower computational and memory access costs than PLL and indicates it may run faster than PLL as a sequential algorithm. We also demonstrate how BVC-PLL algorithm can be extended to handle directed graphs and weighted graphs and how it can utilize the hierarchical parallelism on a modern parallel computing architecture. Extensive experiments on real-world graphs not only show the sequential BVC-PLL can run more than two times faster than the original PLL, but also demonstrates its parallel efficiency and scalability.


title: Topological characterizations of an extended Su-Schrieffer-Heeger model

abstract: The Su-Schrieffer-Heeger (SSH) model perhaps is the easiest and the most basic model for topological excitations. Many variations and extensions of the SSH model have been proposed and explored to better understand both fundamental and novel aspects of topological physics. The SSH4 model has been proposed theoretically as an extended SSH model with higher dimension (the internal dimension changes from two to four). It has been proposed that the winding number in this system can be determined through a higher-dimensional extension of the mean chiral displacement measurement, however this has not yet been verified in experiment. Here we report the realization of this model with ultracold atoms in a momentum lattice. We verify the winding number through measurement of the mean chiral displacement in a system with higher internal dimension, we map out the topological phase transition in this system, and we confirm the topological edge state by observation of the quench dynamics when atoms are initially prepared at the system boundary.


title: Signatures of many-body localization and metastability by weak perturbation

abstract: Nonequilibrium dynamics in isolated quantum many-body systems displays a number of intriguing features, such as many-body localization (MBL) and prethermalization. Here we investigate a simple ladder system with disorder, in which various distinct dynamical features coexist and interplay. By exact diagonalization, we demonstrate that the system exhibits the signatures of an MBL-ergodic-MBL reentrant transition, metastability, and disorder-free MBL. We give an account of these properties by introducing a quasi-particle picture and interpreting the quasi-vacuum energy fluctuation as an effective disorder on the quasi-particle dynamics. It is speculated that the weak perturbation behavior is a finite-size effect, but its relaxation time scale increases with the system size.


title: Densely Residual Laplacian Super-Resolution

abstract: Super-Resolution convolutional neural networks have recently demonstrated high-quality restoration for single images. However, existing algorithms often require very deep architectures and long training times. Furthermore, current convolutional neural networks for super-resolution are unable to exploit features at multiple scales and weigh them equally, limiting their learning capability. In this exposition, we present a compact and accurate super-resolution algorithm namely, Densely Residual Laplacian Network (DRLN). The proposed network employs cascading residual on the residual structure to allow the flow of low-frequency information to focus on learning high and mid-level features. In addition, deep supervision is achieved via the densely concatenated residual blocks settings, which also helps in learning from high-level complex features. Moreover, we propose Laplacian attention to model the crucial features to learn the inter and intra-level dependencies between the feature maps. Furthermore, comprehensive quantitative and qualitative evaluations on low-resolution, noisy low-resolution, and real historical image benchmark datasets illustrate that our DRLN algorithm performs favorably against the state-of-the-art methods visually and accurately.


title: Crossover from charge density wave stabilized antiferromagnetism to superconductivity in Nd$_{1-x}$La$_x$NiC$_2$ compounds

abstract: The path from the charge density wave antiferromagnet NdNiC$_2$ to a noncentrosymmetric superconductor LaNiC$_2$ has been studied by gradual replacement of Nd by La ions. The evolution of physical properties has been explored by structural, magnetic, transport, magnetoresistance and specific heat measurements. With the substitution of La for Nd, the Peierls temperature is gradually suppressed, which falls within the BCS mean-field relation for chemical pressure with a critical concentration of $x_c$ = 0.38. As long as charge density wave is maintained, the antiferromagnetic ground state remains robust against doping and despite of a N\'eel temperature reduction shows a rapid and sharp magnetic transition. Once the CDW is completely suppressed, intermediate compounds of the Nd$_{1-x}$La$_x$NiC$_2$ series reveal symptoms of a gradual softening of the features associated with AFM transition and increase of the spin disorder. Immediately after the antiferromagnetic transition is depressed to zero temperature, the further incorporation of La ions results in the emergence of superconductivity. This crossover in the Nd$_{1-x}$La$_x$NiC$_2$ is discussed in the terms of the possible quantum critical point.


title: Evaluation of Abramowitz functions in the right half of the complex plane

abstract: A numerical scheme is developed for the evaluation of Abramowitz functions $J_n$ in the right half of the complex plane. For $n=-1,\, \ldots,\, 2$, the scheme utilizes series expansions for $|z|<1$ and asymptotic expansions for $|z|>R$ with $R$ determined by the required precision, and modified Laurent series expansions which are precomputed via a least squares procedure to approximate $J_n$ accurately and efficiently on each sub-region in the intermediate region $1\le |z| \le R$. For $n>2$, $J_n$ is evaluated via a recurrence relation. The scheme achieves nearly machine precision for $n=-1, \ldots, 2$, with the cost about four times of evaluating a complex exponential per function evaluation.


title: Direct Learning with Guarantees of the Difference DAG Between Structural Equation Models

abstract: Discovering cause-effect relationships between variables from observational data is a fundamental challenge in many scientific disciplines. However, in many situations it is desirable to directly estimate the change in causal relationships across two different conditions, e.g., estimating the change in genetic expression across healthy and diseased subjects can help isolate genetic factors behind the disease. This paper focuses on the problem of directly estimating the structural difference between two structural equation models (SEMs), having the same topological ordering, given two sets of samples drawn from the individual SEMs. We present an principled algorithm that can recover the difference SEM in $\mathcal{O}(d^2 \log p)$ samples, where $d$ is related to the number of edges in the difference SEM of $p$ nodes. We also study the fundamental limits and show that any method requires at least $\Omega(d' \log \frac{p}{d'})$ samples to learn difference SEMs with at most $d'$ parents per node. Finally, we validate our theoretical results with synthetic experiments and show that our method outperforms the state-of-the-art. Moreover, we show the usefulness of our method by using data from the medical domain.


title: Generation of Quantum Entanglement based on Electromagnetically Induced Transparency Media

abstract: Quantum entanglement is an essential ingredient for the absolute security of quantum communication. Generation of continuous-variable entanglement or two-mode squeezing between light fields based on the effect of electromagnetically induced transparency (EIT) has been systematically investigated in this work. Here, we propose a new scheme to enhance the degree of entanglement between probe and coupling fields of coherent-state light by introducing a two-photon detuning in the EIT system. This proposed scheme is more efficient than the conventional one, utilizing the dephasing rate of ground-state coherence, i.e., the decoherence rate to produce entanglement or two-mode squeezing which adds far more excess fluctuation or noise to the system. In addition, maximum degree of entanglement at a given optical depth can be achieved with a wide range of the coupling Rabi frequency and the two-photon detuning, showing our scheme is robust and flexible. It is also interesting to note that while EIT is the effect in the perturbation limit, i.e. the probe field being much weaker than the coupling field and treated as a perturbation, there exists an optimum ratio of the probe to coupling intensities to achieve the maximum entanglement. Our proposed scheme can advance the continuous-variable-based quantum technology and may lead to applications in quantum communication utilizing squeezed light.


title: Signatures of Self-Interacting dark matter on cluster density profile and subhalo distributions

abstract: Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($\sigma/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.


title: Ultra-Compact Coupling Structures for Heterogeneously Integrated Silicon Lasers

abstract: Due to the inherent in-direct bandgap nature of Silicon, heterogeneous integration of semiconductor lasers on Silicon on Insulator (SOI) is crucial for next-generation on-chip optical interconnects. Compact, high-efficient and high-tolerant couplers between III-V light source and silicon chips have been the challenge for photonic integrated circuit (PIC). Here, we redesign the taper adiabatic coupler with the total coupling length of only 4 {\mu}m, and propose another two novel slot coupler and bridge-SWG coupler with both coupling length of 7 {\mu}m, to heterogeneously integrate III-V lasers and silicon chips. We study theoretically the optical mode coupling process through the redesigned taper coupler, the final coupling results match well with the simulation in 3D-FDTD. The three compact couplers represent fundamental TE mode coupling efficiencies all over 90%, even 95.7% for bridge-SWG coupler, to the best of our knowledge, are also the shortest coupling structures (7 um). Moreover, these coupling structures also possess excellent fabrication tolerance.


title: Learning from Web Data with Self-Organizing Memory Module

abstract: Learning from web data has attracted lots of research interest in recent years. However, crawled web images usually have two types of noises, label noise and background noise, which induce extra difficulties in utilizing them effectively. Most existing methods either rely on human supervision or ignore the background noise. In this paper, we propose a novel method, which is capable of handling these two types of noises together, without the supervision of clean images in the training stage. Particularly, we formulate our method under the framework of multi-instance learning by grouping ROIs (i.e., images and their region proposals) from the same category into bags. ROIs in each bag are assigned with different weights based on the representative/discriminative scores of their nearest clusters, in which the clusters and their scores are obtained via our designed memory module. Our memory module could be naturally integrated with the classification module, leading to an end-to-end trainable system. Extensive experiments on four benchmark datasets demonstrate the effectiveness of our method.


title: A Neural-based Program Decompiler

abstract: Reverse engineering of binary executables is a critical problem in the computer security domain. On the one hand, malicious parties may recover interpretable source codes from the software products to gain commercial advantages. On the other hand, binary decompilation can be leveraged for code vulnerability analysis and malware detection. However, efficient binary decompilation is challenging. Conventional decompilers have the following major limitations: (i) they are only applicable to specific source-target language pair, hence incurs undesired development cost for new language tasks; (ii) their output high-level code cannot effectively preserve the correct functionality of the input binary; (iii) their output program does not capture the semantics of the input and the reversed program is hard to interpret. To address the above problems, we propose Coda, the first end-to-end neural-based framework for code decompilation. Coda decomposes the decompilation task into two key phases: First, Coda employs an instruction type-aware encoder and a tree decoder for generating an abstract syntax tree (AST) with attention feeding during the code sketch generation stage. Second, Coda then updates the code sketch using an iterative error correction machine guided by an ensembled neural error predictor. By finding a good approximate candidate and then fixing it towards perfect, Coda achieves superior performance compared to baseline approaches. We assess Coda's performance with extensive experiments on various benchmarks. Evaluation results show that Coda achieves an average of 82\% program recovery accuracy on unseen binary samples, where the state-of-the-art decompilers yield 0\% accuracy. Furthermore, Coda outperforms the sequence-to-sequence model with attention by a margin of 70\% program accuracy.


title: TensorNetwork on TensorFlow: Entanglement Renormalization for quantum critical lattice models

abstract: We use TensorNetwork [C. Roberts et al., arXiv: 1905.01330], a recently developed API for performing tensor network contractions using accelerated backends such as TensorFlow, to implement an optimization algorithm for the Multi-scale Entanglement Renormalization Ansatz (MERA). We use the MERA to approximate the ground state wave function of the infinite, one-dimensional transverse field Ising model at criticality, and extract conformal data from the optimized ansatz. Comparing run times of the optimization on CPUs vs. GPU, we report a very significant speed-up, up to a factor of 200, of the optimization algorithm when run on a GPU.


title: Magnon gravitomagnetoelectric effect in noncentrosymmetric antiferromagnetic insulators

abstract: We study the magnon contribution to the gravitomagnetoelectric (gravito-ME) effect, in which the magnetization is induced by a temperature gradient, in noncentrosymmetric antiferromagnetic insulators. This phenomenon is totally different from the ME effect, because the temperature gradient is coupled to magnons but an electric field is not. We derive a general formula of the gravito-ME susceptibility in terms of magnon wave functions and find that a difference in $g$ factors of magnetic ions is crucial. We also apply our formula to a specific model. Although the obtained gravito-ME susceptibility is small, we discuss several ways to enhance this phenomenon.


title: A sum involving the greatest-integer function

abstract: We determine properties of the set of values of $ [nx] - ([x]/1 + [2x]/2 + \cdots + [nx]/x) $ as $n$ and $x$ vary.


title: Differences in mechanical properties lead to anomalous phase separation in a model cell co-culture

abstract: During the morphogenesis of tissues and tumors, cells often interact with neighbors with different mechanical properties, but the understanding of its role is lacking. We use active Brownian Dynamics simulations to study a model co-culture consisting of two types of cells with the same size and self-propulsion speed, but different mechanical stiffness and cell-cell adhesion. As time evolves, the system phase separates out into clusters with distinct morphologies and transport properties for the two cell types. The density structure functions and the growth of cell clusters deviate from behavior characteristic of the phase separation in binary fluids. Our results capture emergent structure and motility previously observed in co-culture experiments and provide mechanistic insights into intercellular phase separation during development and disease.


title: Probing Quadratic Gravity with Binary Inspirals

abstract: In this paper, we study gravitational waves generated by binary systems within an extension of General Relativity which is described by the addition of quadratic in curvature tensor terms to the Einstein-Hilbert action. Treating quadratic gravity as an effective theory valid in the low energy/curvature regime, we argue that reliable calculations can be performed in the early inspiral phase, and furthermore, no flux of additional massive waves can be detected. We then compute the massive dipole (-1PN) leading corrections to the post-Newtonian (PN) expansion of the standard waveform. By confronting these theoretical calculations with available experimental data, we constrain both unknown parameters of quadratic gravity to be $0 \leq \gamma \, \lesssim 5.7\cdot 10^{76}$, and $-\frac{\gamma}{4} \leq \beta \, \lesssim - 4.2\cdot 10^{75}$.


title: A Concise Model for Multi-Criteria Chinese Word Segmentation with Transformer Encoder

abstract: Multi-criteria Chinese word segmentation (MCCWS) aims to exploit the relations among the multiple heterogeneous segmentation criteria and further improve the performance of each single criterion. Previous work usually regards MCCWS as different tasks, which are learned together under the multi-task learning framework. In this paper, we propose a concise but effective unified model for MCCWS, which is fully-shared for all the criteria. By leveraging the powerful ability of the Transformer encoder, the proposed unified model can segment Chinese text according to a unique criterion-token indicating the output criterion. Besides, the proposed unified model can segment both simplified and traditional Chinese and has an excellent transfer capability. Experiments on eight datasets with different criteria show that our model outperforms our single-criterion baseline model and other multi-criteria models. Source codes of this paper are available on Github https://github.com/acphile/MCCWS.


title: Synchronicity phenomenon in cluster patterns

abstract: It has been known that several objects such as cluster variables, coefficients, seeds, and $Y$-seeds in different cluster patterns with common exchange matrices share the same periodicity under mutations. We call it synchronicity phenomenon in cluster patterns. In this expository note we explain the mechanism of synchronicity based on several fundamental results on cluster algebra theory such as separation formulas, sign-coherence, Laurent positivity, duality, and detropicalization obtained by several authors. We also show that all synchronicity properties studied in this paper are naturally extended to cluster patterns of generalized cluster algebras, up to the Laurent positivity conjecture.


title: A global Torelli theorem for certain Calabi-Yau threefolds

abstract: We establish a global Torelli theorem for the complete family of Calabi-Yau threefolds arising from cyclic triple covers of $\mathbb P^3$ branched along stable hyperplane arrangements.


title: Analyzing GDPR Compliance Through the Lens of Privacy Policy

abstract: With the arrival of the European Union's General Data Protection Regulation (GDPR), several companies are making significant changes to their systems to achieve compliance. The changes range from modifying privacy policies to redesigning systems which process personal data. This work analyzes the privacy policies of large-scaled cloud services which seek to be GDPR compliant. The privacy policy is the main medium of information dissemination between the data controller and the users. We show that many services that claim compliance today do not have clear and concise privacy policies. We identify several points in the privacy policies which potentially indicate non-compliance; we term these GDPR vulnerabilities. We identify GDPR vulnerabilities in ten cloud services. Based on our analysis, we propose seven best practices for crafting GDPR privacy policies.


title: Supervised Contextual Embeddings for Transfer Learning in Natural Language Processing Tasks

abstract: Pre-trained word embeddings are the primary method for transfer learning in several Natural Language Processing (NLP) tasks. Recent works have focused on using unsupervised techniques such as language modeling to obtain these embeddings. In contrast, this work focuses on extracting representations from multiple pre-trained supervised models, which enriches word embeddings with task and domain specific knowledge. Experiments performed in cross-task, cross-domain and cross-lingual settings indicate that such supervised embeddings are helpful, especially in the low-resource setting, but the extent of gains is dependent on the nature of the task and domain. We make our code publicly available.


title: Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

abstract: We present results from a semicoherent search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the $\mathcal{J}$-statistic, and by analysing data from Advanced LIGO's second observing run. In the frequency range searched, from $60$ to $650\,\mathrm{Hz}$, we find no evidence of gravitational radiation. At $194.6\,\mathrm{Hz}$, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95\% confidence) of $h_0^{95\%} = 3.47 \times 10^{-25}$ when marginalising over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering.


title: Thermalization of non-stochastic Hamiltonian systems

abstract: Ability of dynamical systems to relax to equilibrium has been investigated since the invention of statistical mechanics, which establishes the connection between dynamics of many-body Hamiltonian systems and phenomenological thermodynamics. The key link in this connection is stochasticity, which translates the deterministic evolution of a dynamical system to its probabilistic exploration of the state space. To-date research focuses on determining the conditions of stochasticity for particular systems. Here we propose an alternative agenda and prove thermalization for non-stochastic Hamiltonian systems. This shows that thermalization happens in both stochastic and non-stochastic systems, reducing the need to rely on stochasticity in a "coarse-grained" analysis. The result is valid for an arbitrary classical Hamiltonian system and does not rely on the thermodynamic limit or the particular form of the interaction potential. It utilizes the property of adiabatic invariance, and reveals a deep relation between the structure of the microscopic Hamiltonian and macroscopic thermodynamics.


title: Deformations, renormgroup, symmetries, AdS/CFT

abstract: We consider the deformations of a supersymmetric quantum field theory by adding spacetime-dependent terms to the action. We propose to describe the renormalization of such deformations in terms of some cohomological invariants, a class of solutions of a Maurer-Cartan equation. We consider the strongly coupled limit of $N=4$ supersymmetric Yang-Mills theory. In the context of AdS/CFT correspondence, we explain what corresponds to our invariants in classical supergravity. There is a leg amputation procedure, which constructs a solution of the Maurer-Cartan equation from tree diagramms of SUGRA. We consider a particular example of the beta-deformation. It is known that the leading term of the beta-function is cubic in the parameter of the beta-deformation. We give a cohomological interpretation of this leading term. We conjecture that it is actually encoded in some simpler cohomology class, which is quadratic in the parameter of the beta-deformation.


title: Faster Distributed Deep Net Training: Computation and Communication Decoupled Stochastic Gradient Descent

abstract: With the increase in the amount of data and the expansion of model scale, distributed parallel training becomes an important and successful technique to address the optimization challenges. Nevertheless, although distributed stochastic gradient descent (SGD) algorithms can achieve a linear iteration speedup, they are limited significantly in practice by the communication cost, making it difficult to achieve a linear time speedup. In this paper, we propose a computation and communication decoupled stochastic gradient descent (CoCoD-SGD) algorithm to run computation and communication in parallel to reduce the communication cost. We prove that CoCoD-SGD has a linear iteration speedup with respect to the total computation capability of the hardware resources. In addition, it has a lower communication complexity and better time speedup comparing with traditional distributed SGD algorithms. Experiments on deep neural network training demonstrate the significant improvements of CoCoD-SGD: when training ResNet18 and VGG16 with 16 Geforce GTX 1080Ti GPUs, CoCoD-SGD is up to 2-3$\times$ faster than traditional synchronous SGD.


title: Steady-state squeezing and entanglement in a dissipatively coupled NOPO network

abstract: We investigate the steady-state photon-number squeezing and quantum entanglement in a network of nondegenerate optical parametric oscillators (NOPOs). We treat each NOPO with Shen's Raman laser model, whose lasing mode provides a photon-number-squeezed state. Two dissipatively coupled NOPOs satisfy Hillery-Zubairy's $HZ1$ entanglement criterion if they are pumped far above the threshold and the dissipative coupling is sufficiently larger than the NOPO cavity loss.


title: 4K-Memristor Analog-Grade Passive Crossbar Circuit

abstract: The superior density of passive analog-grade memristive crossbars may enable storing large synaptic weight matrices directly on specialized neuromorphic chips, thus avoiding costly off-chip communication. To ensure efficient use of such crossbars in neuromorphic computing circuits, variations of current-voltage characteristics of crosspoint devices must be substantially lower than those of memory cells with select transistors. Apparently, this requirement explains why there were so few demonstrations of neuromorphic system prototypes using passive crossbars. Here we report a 64x64 passive metal-oxide memristor crossbar circuit with ~99% device yield, based on a foundry-compatible fabrication process featuring etch-down patterning and low-temperature budget, conducive to vertical integration. The achieved ~26% variations of switching voltages of our devices were sufficient for programming 4K-pixel gray-scale patterns with an average tuning error smaller than 4%. The analog properties were further verified by experimentally demonstrating MNIST pattern classification with a fidelity close to the software-modeled limit for a network of this size, with an ~1% average error of import of ex-situ-calculated synaptic weights. We believe that our work is a significant improvement over the state-of-the-art passive crossbar memories in both complexity and analog properties.


title: On uniqueness of two meromorphic functions sharing a small function

abstract: In this paper, we have investigated the uniqueness problems of entire and meromorphic functions concerning differential polynomials sharing a small function. Our results radically extended and improved the results of Bhoosnurmath-Pujari and Harina - Anand not only by sharing small function instead of fixed point but also reducing the lower bound of $ n $. The authors Harina-Anand made plenty of mistakes in their paper. We have corrected all of them in a more convenient way. At last some open questions have been posed for further study in this direction.


title: Evolution of longevity, age at last birth and sexual conflict with grandmothering

abstract: We use a two-sex partial differential equation (PDE) model based on the Grandmother Hypothesis. We build on an earlier model by Kim et al. [27] by allowing for evolution in both longevity and age at last birth, and also assuming that post-fertile females support only their daughters' fertility. Similarly to Kim et al. [27], we find that only two locally stable equilibria exist: one corresponding to great ape-like longevities and the other corresponding to hunter-gatherer longevities. Our results show that grandmothering enables the transition between these two equilibria, without extending the end of fertility. Moreover, sensitivity analyses of the model show that male competition, arising from a skew in the mating sex ratio towards males, plays a significant role in determining whether the transition from great ape-like longevities to higher longevities is possible and the equilibrium value of the average adult lifespan. Whereas grandmothering effects have a significant impact on the equilibrium value of the average age at last birth and enable the transition to higher longevities, they have an insignificant impact on the equilibrium value of the average adult lifespan.


title: Why does women's fertility end in mid-life? Grandmothering and age at last birth

abstract: Great apes, the other living members of our hominid family, become decrepit before the age of forty and rarely outlive their fertile years. In contrast, women - even in high mortality hunter-gatherer populations - usually remain healthy and productive well beyond menopause. The grandmother hypothesis aims to account for the evolution of this distinctive feature of human life history. Our previous mathematical simulations of that hypothesis fixed the end of female fertility at the age of 45, based on the similarities among living hominids, and then modeled the evolution of human-like longevity from an ancestral state, like that of the great apes, due only to grandmother effects. A major modification here allows the age female fertility ends to vary as well, directly addressing a version of the question, influentially posed by GC Williams six decades ago: Why isn't menopause later in humans? Our model is an agent-based model (ABM) that accounts for the coevolution of both expected adult lifespan and end of female fertility as selection maximizes reproductive value. We find that grandmother effects not only drive the population from an equilibrium representing a great ape-like longevity to a new human-like longevity, they also maintain the observed termination of women's fertility before the age of 50.


title: Dynamics of groups of magnetically driven artificial microswimmers

abstract: Magnetically driven artificial microswimmers have the potential to revolutionize many biomedical technologies, such as minimally-invasive microsurgery, micro-particle manipulation, and localized drug delivery. However, many of these applications will require the controlled dynamics of teams of these micro-robots with minimal feedback. In this work, we study the motion and fluid dynamics produced by groups of artificial microswimmers driven by a torque induced through a uniform, rotating magnetic field. Through Stokesian dynamics simulations, we show that the swimmer motion produces a rotational velocity field in the plane orthogonal to the direction of the magnetic field's rotation, which causes two interacting swimmers to move in circular trajectories in this plane around a common center. The resulting over all motion is on a helical trajectory for the swimmers. We compare the highly rotational velocity field of the fluid to the velocity field generated by a rotlet, the point-torque singularity of Stokes flows, showing that this is a reasonable approximation on the time average. Finally, we study the motion of larger groups of swimmers and show that these groups tend to move coherently, especially when swimmer magnetizations are uniform. This coherence is achieved because the group center remains almost constant in the plane orthogonal to the net motion of the swimmers. The results in the paper will prove useful for controlling the ensemble dynamics of small collections of magnetic swimmers.


title: Adult sex ratio as an index for male strategy in primates

abstract: The adult sex ratio (ASR) is defined as the number of fertile males divided by the number of fertile females in a population. We build an ODE model with minimal age structure, in which males compete for paternities using either a multiple-mating or searching-then-guarding strategy, to investigate the value of ASR as an index for predicting which strategy males will adopt, with a focus in our investigation on the differences of strategy choice between chimpanzees Pan troglodytes and human hunter-gatherers Homo sapiens. Parameters in the model characterise aspects of life history and behaviour, and determine both dominant strategy and the ASR when the population is at or near equilibrium. Sensitivity analysis on the model parameters informs us that ASR is strongly influenced by parameters characterising life history, while dominant strategy is affected most strongly by the effectiveness of guarding (average length of time a guarded pair persists, and resistance to paternity theft) and moderately by some life history traits. For fixed effectiveness of guarding and other parameters, dominant strategy tends to change from multiple mating to guarding along a curve that aligns well with a contour of constant ASR, under variation of parameters such as longevity and age female fertility ends. This confirms the hypothesis that ASR may be a useful index for predicting the optimal male mating strategy, provided we have some limited information about ecology and behaviour.


title: A Taxonomy for Virtual and Augmented Reality in Education

abstract: In this paper, a taxonomy for VR/AR in education is presented that can help differentiate and categorise education experiences and provide indication as to why some applications of fail whereas others succeed. Examples will be presented to illustrate the taxonomy, including its use in developing and planning two current VR projects in our laboratory. The first project is a VR application for the training of Chemical Engineering students (and potentially industrial operators) on the use of a physical pilot plant facility. The second project involves the use of VR cinematography for enacting ethics scenarios (and thus ethical awareness and development) pertinent to engineering work situations.


title: Irregular vertex algebras

abstract: We introduce the notion of irregular vertex (operator) algebras. The irregular versions of fundamental properties, such as Goddard uniqueness theorem, associativity and operator product expansions are formulated and proved. We also give some elementary examples of irregular vertex operator algebras.


title: Generating Normal Networks via Leaf Insertion and Nearest Neighbor Interchange

abstract: Galled trees are studied as a recombination model in theoretic population genetics. This class of phylogenetic networks has been generalized to tree-child networks, normal networks and tree-based networks by relaxing a structural condition. Although these networks are simple, their topological structures have yet to be fully understood. It is well-known that all phylogenetic trees on $n$ taxa can be generated by the insertion of the $n$-th taxa to each edge of all the phylogenetic trees on $n-1$ taxa. We prove that all tree-child networks with $k$ reticulate nodes on $n$ taxa can be uniquely generated via three operations from all the tree-child networks with $k-1$ or $k$ reticulate nodes on $n-1$ taxa . An application of this result is found in counting tree-child networks and normal networks. In particular, a simple formula is given for the number of rooted phylogenetic networks with one reticulate node.


title: On the equational graphs over finite fields

abstract: In this paper, we generalize the notion of functional graph. Specifically, given an equation $E(X,Y) = 0$ with variables $X$ and $Y$ over a finite field $\mathbb{F}_q$ of odd characteristic, we define a digraph by choosing the elements in $\mathbb{F}_q$ as vertices and drawing an edge from $x$ to $y$ if and only if $E(x,y)=0$. We call this graph as equational graph. In this paper, we study the equational graphs when choosing $E(X,Y) = (Y^2 - f(X))(\lambda Y^2 - f(X))$ with $f(X)$ a polynomial over $\mathbb{F}_q$ and $\lambda$ a non-square element in $\mathbb{F}_q$. We show that if $f$ is a permutation polynomial over $\mathbb{F}_q$, then every connected component of the graph has a Hamiltonian cycle. Moreover, these Hamiltonian cycles can be used to construct balancing binary sequences. By making computations for permutation polynomials $f$ of low degree, it appears that almost all these graphs are strongly connected, and there are many Hamiltonian cycles in such a graph if it is connected.


title: Non-Fierz-Pauli bimetric theory from quadratic curvature gravity on Einstein manifolds

abstract: We show that, in four-dimensional spacetimes with an arbitrary Einstein metric, with and without a cosmological constant, perturbative dynamical degrees of freedom in generic quadratic-curvature gravity can be decoupled into massless and massive parts. The massive part has the structure identical to, modulo the over-all sign, the non-Fierz-Pauli-type massive gravity, and a further decomposition into the spin-2 and spin-0 sectors can be done. The equivalence at the level of equations of motion allows us to translate various observational bounds on the mass of extra fields into constraints on the coupling constants in quadratic curvature gravity. We find that the Weyl-squared term is confronting two apparently contradicting constraints on massive spin-2 fields from the inverse-square law experiments and observations of spinning black holes.


title: DP-LSSGD: A Stochastic Optimization Method to Lift the Utility in Privacy-Preserving ERM

abstract: Machine learning (ML) models trained by differentially private stochastic gradient descent (DP-SGD) have much lower utility than the non-private ones. To mitigate this degradation, we propose a DP Laplacian smoothing SGD (DP-LSSGD) to train ML models with differential privacy (DP) guarantees. At the core of DP-LSSGD is the Laplacian smoothing, which smooths out the Gaussian noise used in the Gaussian mechanism. Under the same amount of noise used in the Gaussian mechanism, DP-LSSGD attains the same DP guarantee, but in practice, DP-LSSGD makes training both convex and nonconvex ML models more stable and enables the trained models to generalize better. The proposed algorithm is simple to implement and the extra computational complexity and memory overhead compared with DP-SGD are negligible. DP-LSSGD is applicable to train a large variety of ML models, including DNNs. The code is available at \url{https://github.com/BaoWangMath/DP-LSSGD}.


title: Phonon-laser sensing in a hetero optomechanical crystal cavity

abstract: Micro- and nanomechanical resonators have emerged as promising platforms for sensing a broad range of physical properties such as mass, force, torque, magnetic field, and acceleration. The sensing performance relies critically on the motional mass, the mechanical frequency, and the linewidth of the mechanical resonator. Here, we demonstrate a hetero optomechanical crystal (OMC) cavity based on a silicon nanobeam structure. The cavity supports phonon lasing in a fundamental mechanical mode with a frequency of 5.91 GHz, an effective mass of 116 fg, and a mechanical linewidth narrowing from 3.3 MHz to 5.2 kHz, while the optomechanical coupling rate of is as high as 1.9 MHz. With this phonon laser, the on-chip sensing with a resolution of $\delta$$\lambda$/$\lambda$ = 1.0*10-8 can be attained, which is at least two orders of magnitude larger than that obtained with conventional silicon-based sensors. The use of a silicon-based hetero OMC cavity that harnesses phonon lasing could pave the way towards exciting, high-precision sensors that lend themselves to silicon monolithic integration and offer unprecedented sensitivity for broad physical sensing applications.


title: Holonomic Gates in Pseudo-Hermitian Quantum Systems

abstract: The time-dependent pseudo-Hermitian formulation of quantum mechanics allows to study open system dynamics in analogy to Hermitian quantum systems. In this setting, we show that the notion of holonomic quantum computation can equally be formulated for pseudo-Hermitian systems. Starting from a degenerate pseudo-Hermitian Hamiltonian we show that, in the adiabatic limit, a non-Abelian geometric phase emerges which realizes a pseudounitary quantum gate. We illustrate our findings by studying a pseudo-Hermitian gain/loss system which can be written in the form of a tripod Hamiltonian by using the biorthogonal representation. It is shown that this system allows for arbitrary pseudo-$\mathrm{U}(2)$ transformations acting on the dark subspace of the system.


title: The blow-up rate for a non-scaling invariant semilinear wave equations

abstract: We consider the semilinear wave equation $$\partial_t^2 u -\Delta u =f(u), \quad (x,t)\in \mathbb{R}^N\times [0,T),\qquad (1)$$ with $f(u)=|u|^{p-1}u\log^a (2+u^2)$, where $p>1$ and $a\in \mathbb{R}$. We show an upper bound for any blow-up solution of (1). Then, in the one space dimensional case, using this estimate and the logarithmic property, we prove that the exact blow-up rate of any singular solution of (1) is given by the ODE solution associated with $(1)$, namely $u'' =|u|^{p-1}u\log^a (2+u^2)$ Unlike the pure power case ($g(u)=|u|^{p-1}u$) the difficulties here are due to the fact that equation (1) is not scale invariant.


title: Delayed spin-up and persistent shift phenomena of Crab pulsar glitches: two sides of the same coin?

abstract: Pulsar glitches are sudden increase in their spin frequency, in most cases followed by the long timescale recovery process. As of this writing, about 546 glitches have been reported in 188 pulsars, the Crab pulsar is a special one with unique manifestations. This writing presents a statistic study on post-glitch observables of the Crab pulsar, especially the delayed spin-up in post-glitch phase and persistent shift in the slow-down rate of the star. By analyzing the radio data over 45 years, we find that two power law functions respectively fit the persistent shift and delayed spin-up timescales versus glitch size well, and we find a linear correlation between the persistent shift and delayed spin-up timescale from the consistency of the two fitting functions, probably indicating their same physical origin and may provide a new probe of interior physics of neutron stars.


title: Learning to Cope with Adversarial Attacks

abstract: The security of Deep Reinforcement Learning (Deep RL) algorithms deployed in real life applications are of a primary concern. In particular, the robustness of RL agents in cyber-physical systems against adversarial attacks are especially vital since the cost of a malevolent intrusions can be extremely high. Studies have shown Deep Neural Networks (DNN), which forms the core decision-making unit in most modern RL algorithms, are easily subjected to adversarial attacks. Hence, it is imperative that RL agents deployed in real-life applications have the capability to detect and mitigate adversarial attacks in an online fashion. An example of such a framework is the Meta-Learned Advantage Hierarchy (MLAH) agent that utilizes a meta-learning framework to learn policies robustly online. Since the mechanism of this framework are still not fully explored, we conducted multiple experiments to better understand the framework's capabilities and limitations. Our results shows that the MLAH agent exhibits interesting coping behaviors when subjected to different adversarial attacks to maintain a nominal reward. Additionally, the framework exhibits a hierarchical coping capability, based on the adaptability of the Master policy and sub-policies themselves. From empirical results, we also observed that as the interval of adversarial attacks increase, the MLAH agent can maintain a higher distribution of rewards, though at the cost of higher instabilities.


title: Effect of aggregation on adsorption phenomena

abstract: Adsorption at an attractive surface in a system with particles self-assembling into small clusters is studied by Molecular dynamics (MD) simulation. We assume Lennard-Jones plus repulsive Yukawa tail interactions, and focus on small densities. The relative increase of the temperature at the critical cluster concentration near the attractive surface (CCCS) shows a power-law dependence on the strength of the wall-particle attraction. At temperatures below the CCCS, the adsorbed layer consists of undeformed clusters if the wall-particle attraction is not too strong. Above the CCCS, or for strong attraction leading to flattening of the adsorbed aggregates, we obtain a monolayer that for strong or very strong attraction consists of flattened clusters or stripes respectively. The accumulated repulsion from the particles adsorbed at the wall leads to a repulsive barrier that slows down the adsorption process, and the accession time grows rapidly with the strength of the wall-particle attraction. Beyond the adsorbed layer of particles, a depletion region of a thickness comparable with the range of the repulsive tail of interactions occurs, and the density in this region decreases with increasing strength of the wall-particle attraction. At larger separations, the exponentially damped oscillations of density agree with theoretical predictions for self-assembling systems. Structural and thermal properties of the bulk are also determined. In particular, a new structural crossover associated with the maximum of the specific heat, and a double-peaked histogram of the cluster size distribution are observed.


title: Bias-Variance Trade-Off in Hierarchical Probabilistic Models Using Higher-Order Feature Interactions

abstract: Hierarchical probabilistic models are able to use a large number of parameters to create a model with a high representation power. However, it is well known that increasing the number of parameters also increases the complexity of the model which leads to a bias-variance trade-off. Although it is a classical problem, the bias-variance trade-off between hidden layers and higher-order interactions have not been well studied. In our study, we propose an efficient inference algorithm for the log-linear formulation of the higher-order Boltzmann machine using a combination of Gibbs sampling and annealed importance sampling. We then perform a bias-variance decomposition to study the differences in hidden layers and higher-order interactions. Our results have shown that using hidden layers and higher-order interactions have a comparable error with a similar order of magnitude and using higher-order interactions produce less variance for smaller sample size.


title: Background Subtraction using Adaptive Singular Value Decomposition

abstract: An important task when processing sensor data is to distinguish relevant from irrelevant data. This paper describes a method for an iterative singular value decomposition that maintains a model of the background via singular vectors spanning a subspace of the image space, thus providing a way to determine the amount of new information contained in an incoming frame. We update the singular vectors spanning the background space in a computationally efficient manner and provide the ability to perform block-wise updates, leading to a fast and robust adaptive SVD computation. The effects of those two properties and the success of the overall method to perform a state of the art background subtraction are shown in both qualitative and quantitative evaluations.


title: Photon trapping in static axially symmetric spacetime

abstract: Recently, several new characteristics have been introduced to describe null geodesic structure of strong gravitational field, such as photon regions, transversely trapping surfaces and some generalizations. They give an alternative and concise way to describe lensing and shadow features of compact objects with strong gravitational field without recurring to complete integration of the geodesic equations. Here we test this construction in the case of the Weyl metrics when geodesic equations are non-separable, and thus can not be integrated analytically, while the above characteristic surfaces and regions can be described in a closed form. We develop further our formalism for a class of static axially symmetric spacetimes introducing more detailed specification of transversely trapping surfaces in terms of their principal curvatures. Surprisingly, we find in the static case without spherical symmetry certain features, such as photon regions, previously known in the Kerr space. These photon regions can be regarded as photon spheres, "thickened" due to oblateness of the metric.


title: Pinpointing Performance Inefficiencies in Java

abstract: Many performance inefficiencies such as inappropriate choice of algorithms or data structures, developers' inattention to performance, and missed compiler optimizations show up as wasteful memory operations. Wasteful memory operations are those that produce/consume data to/from memory that may have been avoided. We present, JXPerf, a lightweight performance analysis tool for pinpointing wasteful memory operations in Java programs. Traditional byte-code instrumentation for such analysis (1) introduces prohibitive overheads and (2) misses inefficiencies in machine code generation. JXPerf overcomes both of these problems. JXPerf uses hardware performance monitoring units to sample memory locations accessed by a program and uses hardware debug registers to monitor subsequent accesses to the same memory. The result is a lightweight measurement at machine-code level with attribution of inefficiencies to their provenance: machine and source code within full calling contexts. JXPerf introduces only 7% runtime overhead and 7% memory overhead making it useful in production. Guided by JXPerf, we optimize several Java applications by improving code generation and choosing superior data structures and algorithms, which yield significant speedups.


title: Valuative dimension and monomial orders

abstract: The main result from this note provides a constructive characterization of the valuative dimension, which bears a strong analogy to Lombardi's constructive characterization of the Krull dimension. While Lombardi's characterization uses the lexicographic monomial order, ours uses the graded (reverse) lexicographic order or, in fact, any graded rational monomial order. Apart from this, the paper contains some related results and some examples which readers may find illuminating.


title: Lost in Translation: Loss and Decay of Linguistic Richness in Machine Translation

abstract: This work presents an empirical approach to quantifying the loss of lexical richness in Machine Translation (MT) systems compared to Human Translation (HT). Our experiments show how current MT systems indeed fail to render the lexical diversity of human generated or translated text. The inability of MT systems to generate diverse outputs and its tendency to exacerbate already frequent patterns while ignoring less frequent ones, might be the underlying cause for, among others, the currently heavily debated issues related to gender biased output. Can we indeed, aside from biased data, talk about an algorithm that exacerbates seen biases?


title: A neighbourhood theorem for submanifolds in generalized complex geometry

abstract: We study neighbourhoods of submanifolds in generalized complex geometry. Our first main result provides sufficient criteria for such a submanifold to admit a neighbourhood on which the generalized complex structure is B-field equivalent to a holomorphic Poisson structure. This is intimately tied with our second main result, which is a rigidity theorem for generalized complex deformations of holomorphic Poisson structures. Specifically, on a compact manifold with boundary we provide explicit conditions under which any generalized complex perturbation of a holomorphic Poisson structure is B-field equivalent to another holomorphic Poisson structure. The proofs of these results require two analytical tools: Hodge decompositions on almost complex manifolds with boundary, and the Nash-Moser algorithm. As a concrete application of these results, we show that on a four-dimensional generalized complex submanifold which is generically symplectic, a neighbourhood of the entire complex locus is B-field equivalent to a holomorphic Poisson structure. Furthermore, we use the neighbourhood theorem to develop the theory of blowing down submanifolds in generalized complex geometry.


title: Highly angular resolving beam separator based on total internal reflection

abstract: We present an optical element for the separation of superimposed beams which only differ in angle. The beams are angularly resolved and separated by total internal reflection at an air gap between two prisms. As a showcase application, we demonstrate the separation of superimposed beams of different diffraction orders directly behind acousto-optic modulators for an operating wavelength of 800nm. The wavelength as well as the component size can easily be adapted to meet the requirements of a wide variety of applications. The presented optical element allows to reduce the lengths of beam paths and thus to decrease laser system size and complexity.


title: Dynamical $t/U$ Expansion of the Doped Hubbard Model

abstract: We construct a new $U(1)$ slave spin representation for the single band Hubbard model in the large-$U$ limit. The mean-field theory in this representation is more amenable to describe both the spin-charge-separation physics of the Mott insulator at half-filling and the strange metal behavior at finite doping. By employing a dynamical Green's function theory for slave spins, we calculate the single-particle spectral function of electrons, and the result is comparable to that in dynamical mean field theories. We then formulate a dynamical $t/U$ expansion for the doped Hubbard model that reproduces the mean-field results at the lowest order of expansion. To the next order of expansion, it naturally yields an effective low-energy theory of a $t-J$ model for spinons self-consistently coupled to n $XXZ$ model for the slave spins. We show that the superexchange $J$ is renormalized by doping, in agreement with the Gutzwiller approximation. Surprisingly, we find a new ferromagnetic channel of exchange interactions which survives in the infinite $U$ limit, as a manifestation of the Nagaoka ferromagnetism.


title: Multiple Testing and Variable Selection along the path of the Least Angle Regression

abstract: We investigate multiple testing and variable selection using the Least Angle Regression (LARS) algorithm in high dimensions under the assumption of Gaussian noise. LARS is known to produce a piecewise affine solution path with change points referred to as the knots of the LARS path. The key to our results is an expression in closed form of the exact joint law of a $K$-tuple of knots conditional on the variables selected by LARS, namely the so-called post-selection joint law of the LARS knots. Numerical experiments demonstrate the perfect fit of our findings. This paper makes three main contributions. First, we build testing procedures on variables entering the model along the LARS path in the general design case when the noise level can be unknown. These testing procedures are referred to as the Generalized $t$-Spacing tests (GtSt) and we prove that they have an exact non-asymptotic level (i.e., the Type I error is exactly controlled). This extends work of (Taylor et al., 2014) where the spacing test works for consecutive knots and known variance. Second, we introduce a new exact multiple false negatives test after model selection in the general design case when the noise level may be unknown. We prove that this testing procedure has exact non-asymptotic level for general design and unknown noise level. Third, we give an exact control of the false discovery rate under orthogonal design assumption. Monte Carlo simulations and a real data experiment are provided to illustrate our results in this case. Of independent interest, we introduce an equivalent formulation of the LARS algorithm based on a recursive function.

categories: math.ST cs.IT math.IT stat.ME stat.ML stat.TH


title: Access Balancing in Storage Systems by Labeling Partial Steiner Systems

abstract: Storage architectures ranging from minimum bandwidth regenerating encoded distributed storage systems to declustered-parity RAIDs can be designed using dense partial Steiner systems in order to support fast reads, writes, and recovery of failed storage units. In order to ensure good performance, popularities of the data items should be taken into account and the frequencies of accesses to the storage units made as uniform as possible. A proposed combinatorial model ranks items by popularity and assigns data items to elements in a dense partial Steiner system so that the sums of ranks of the elements in each block are as equal as possible. By developing necessary conditions in terms of independent sets, we demonstrate that certain Steiner systems must have a much larger difference between the largest and smallest block sums than is dictated by an elementary lower bound. In contrast, we also show that certain dense partial $S(t, t+1, v)$ designs can be labeled to realize the elementary lower bound. Furthermore, we prove that for every admissible order $v$, there is a Steiner triple system $(S(2, 3, v))$ whose largest difference in block sums is within an additive constant of the lower bound.


title: Recursion scheme for the largest $\beta$-Wishart-Laguerre eigenvalue and Landauer conductance in quantum transport

abstract: The largest eigenvalue distribution of the Wishart-Laguerre ensemble, indexed by Dyson parameter $\beta$ and Laguerre parameter $a$, is fundamental in multivariate statistics and finds applications in diverse areas. Based on a generalization of the Selberg integral, we provide an effective recursion scheme to compute this distribution explicitly in both the original model, and a fixed-trace variant, for $a,\beta$ non-negative integers and finite matrix size. For $\beta = 2$ this circumvents known symbolic evaluation based on determinants which become impractical for large dimensions. Our exact results have immediate applications in the areas of multiple channel communication and bipartite entanglement. Moreover, we are also led to the exact solution of a long standing problem of finding a general result for Landauer conductance distribution in a chaotic mesoscopic cavity with two ideal leads. Thus far, exact closed-form results for this were available only in the Fourier-Laplace space or could be obtained on a case-by-case basis.

categories: math-ph cond-mat.mes-hall cond-mat.stat-mech math.MP stat.AP stat.CO


title: A linear method for camera pair self-calibration and multi-view reconstruction with geometrically verified correspondences

abstract: We examine 3D reconstruction of architectural scenes in unordered sets of uncalibrated images. We introduce a linear method to self-calibrate and find the metric reconstruction of a camera pair. We assume unknown and different focal lengths but otherwise known internal camera parameters and a known projective reconstruction of the camera pair. We recover two possible camera configurations in space and use the Cheirality condition, that all 3D scene points are in front of both cameras, to disambiguate the solution. We show in two Theorems, first that the two solutions are in mirror positions and then the relations between their viewing directions. Our new method performs on par (median rotation error $\Delta R = 3.49^{\circ}$) with the standard approach of Kruppa equations ($\Delta R = 3.77^{\circ}$) for self-calibration and 5-Point algorithm for calibrated metric reconstruction of a camera pair. We reject erroneous image correspondences by introducing a method to examine whether point correspondences appear in the same order along $x, y$ image axes in image pairs. We evaluate this method by its precision and recall and show that it improves the robustness of point matches in architectural and general scenes. Finally, we integrate all the introduced methods to a 3D reconstruction pipeline. We utilize the numerous camera pair metric recontructions using rotation-averaging algorithms and a novel method to average focal length estimates.


title: n-dimensional PDM non-linear oscillators: Linearizability and Euler-Lagrange or Newtonian invariance

abstract: We argue that, under multidimensional position-dependent mass (PDM) settings, the Euler-Lagrange textbook invariance falls short and turned out to be vividly incomplete and/or insecure for a set of PDM-Lagrangians. We show that the transition from Euler-Lagrange component presentation to Newtonian vector presentation is necessary and vital to guarantee invariance. The totality of the Newtonian vector equations of motion is shown to be more comprehensive and instructive than the Euler-Lagrange component equations of motion (they do not run into conflict with each other though). We have successfully used the Newtonian invariance amendment, along with some nonlocal space-time point transformation recipe, to extract exact solutions for a set of n-dimensional nonlinear PDM-oscillators. They are, Mathews-Lakshmanan type-I PDM-oscillators, power-law type-I PDM-oscillators, the Mathews-Lakshmanan type-II PDM-oscillators, the power-law type-II PDM-oscillators, and some nonlinear shifted Mathews-Lakshmanan type-I PDM-oscillators.


title: Large scale Lasso with windowed active set for convolutional spike sorting

abstract: Spike sorting is a fundamental preprocessing step in neuroscience that is central to access simultaneous but distinct neuronal activities and therefore to better understand the animal or even human brain. But numerical complexity limits studies that require processing large scale datasets in terms of number of electrodes, neurons, spikes and length of the recorded signals. We propose in this work a novel active set algorithm aimed at solving the Lasso for a classical convolutional model. Our algorithm can be implemented efficiently on parallel architecture and has a linear complexity w.r.t. the temporal dimensionality which ensures scaling and will open the door to online spike sorting. We provide theoretical results about the complexity of the algorithm and illustrate it in numerical experiments along with results about the accuracy of the spike recovery and robustness to the regularization parameter.

categories: stat.ML cs.LG eess.SP math.OC stat.AP stat.CO


title: Metasurface integrated Vertical Cavity Surface Emitting Lasers for programmable directional lasing emissions

abstract: Vertical cavity surface emitting lasers (VCSELs) have made indispensable contributions to the development of modern optoelectronic technologies. However, arbitrary beam shaping of VCSELs within a compact system still remains inaccessible till now. The emerging ultra-thin flat optical structures, namely metasurfaces, offer a powerful technique to manipulate electromagnetic fields with subwavelength spatial resolution. Here, we show that the monolithic integration of dielectric metasurfaces with VCSELs enables remarkable arbitrary control of the laser beam profiles, including self-collimation, Bessel and Vortex lasers, with high efficiency. Such wafer-level integration of metasurface through VCSELs-compatible technology simplifies the assembling process and preserves the high performance of the VCSELs. We envision that our approach can be implemented in various wide-field applications, such as optical fibre communications, laser printing, smartphones, optical sensing, face recognition, directional displays and ultra-compact light detection and ranging (LiDAR).


title: Representations of simple Jordan superalgebras

abstract: This paper completes description of categories of representations of finite-dimensional simple unital Jordan superalgebras over algebraically closed field of characteristic zero.


title: Identification of Time-varying in situ Signals in Quantum Circuits

abstract: The identification of time-varying \textit{in situ} signals is crucial for characterizing the dynamics of quantum processes occurring in highly isolated environments. Under certain circumstances, they can be identified from time-resolved measurements via Ramsey interferometry experiments, but only with very special probe systems can the signals be explicitly read out, and a theoretical analysis is lacking on whether the measurement data are sufficient for unambiguous identification. In this paper, we formulate this problem as the invertibility of the underlying quantum input-output system, and derive the algebraic identifiability criterion and the algorithm for numerically identifying the signals. The criterion and algorithm can be applied to both closed and open quantum systems, and their effectiveness is demonstrated by numerical examples.


title: Ground-state cooling of an magnomechanical resonator induced by magnetic damping

abstract: Quantum manipulation of mechanical resonators has been widely applied in fundamental physics and quantum information processing. Among them, cooling the mechanical system to its quantum ground state is regarded as a key step. In this work, we propose a scheme which one can realize ground-state cooling of resonator in a cavity magnomechanical system. The system consists of a microwave cavity and a small ferromagnetic sphere, in which phonon-magnon coupling and cavity photon-magnon coupling can be achieved via magnetostrictive interaction and magnetic dipole interaction, respectively. After adiabatically eliminating the cavity mode, an effective Hamiltonian which consists of magnon and mechanical modes is obtained. Within experimentally feasible parameters, we demonstrate that the ground-state cooling of the magnomechanical resonator can be achieved by extra magnetic damping. Unlike optomechanical cooling, magnomechanical interaction is utilized to realize the cooling of resonators. We further illustrate the ground-state cooling can be effectively controlled by the external magnetic field.


title: Sample Efficient Learning of Path Following and Obstacle Avoidance Behavior for Quadrotors

abstract: In this paper we propose an algorithm for the training of neural network control policies for quadrotors. The learned control policy computes control commands directly from sensor inputs and is hence computationally efficient. An imitation learning algorithm produces a policy that reproduces the behavior of a path following control algorithm with collision avoidance. Due to the generalization ability of neural networks, the resulting policy performs local collision avoidance of unseen obstacles while following a global reference path. The algorithm uses a time-free model predictive path-following controller as a supervisor. The controller generates demonstrations by following few example paths. This enables an easy to implement learning algorithm that is robust to errors of the model used in the model predictive controller. The policy is trained on the real quadrotor, which requires collision-free exploration around the example path. An adapted version of the supervisor is used to enable exploration. Thus, the policy can be trained from a relatively small number of examples on the real quadrotor, making the training sample efficient.


title: Analytical design of flat-top transmission filters composed of several resonant structures

abstract: Resonant properties of composite structures consisting of several identical resonant diffractive structures (e.g. multilayer thin-film structures or guided-mode resonance gratings) separated by phase-shift layers are investigated theoretically. Using the scattering matrix formalism, we analytically demonstrate that, at properly chosen thicknesses of the phase-shift layers, the composite structures comprising two or four resonant diffractive structures with a Lorentzian transmittance profile optically implement the Butterworth filters of the order two or three, respectively, and enable achieving flat-top transmission spectra with steep slopes and low sidebands. In addition, we show that the composite structures consisting of three or four second-order Butterworth filters can accurately approximate the fourth- or fifth-order Butterworth filters, respectively. The presented theoretical results are confirmed by rigorous numerical simulations of composite structures consisting of the so-called W-structures (simple three-layer resonant structures comprising a high-index core layer and two low-index cladding layers in a high-index dielectric environment). The simulation results confirm the formation of flat-top transmittance peaks, the shape of which fully agrees with the derived theoretical description. Moreover, we demonstrate an exceptionally simple mechanism of controlling the transmittance peak width, which consists in changing the thicknesses of the cladding layers of the initial W-structure and enables generating flat-top transmission peaks with a significantly subnanometer size.


title: Clearly Discriminate the Continuum Band and Exciton State of the Hybrid Lead Bromide Perovskite

abstract: Electronic states of the hybrid perovskite enable their promising applications as distinctive optoelectronic materials. The understanding of their electronic structures and charge characters remains highly controversial. The electronic mechanism such as reabsorption, Urbach tail and indirect band for interpreting dual-peak emissions is one of the controversial focuses. Herein, we report that through heterojunction enhanced exciton dissociation and global tracing of multiple radiative electronic states across wide temperature regions, we have succeeded in directly observing free carrier emissions from the hybrid lead bromide perovskite and clearly discriminating the direct continuum band and exciton states. The widely-concerned dual-peak emissions are clarified to be excitonic, arising from two types of exciton states of the perovskite. These excitons possess giant binding energies and superior phase stability compared to conventional inorganic semiconductors, providing important implications for exploiting the excitonic mechanism for realizing novel optoelectronic applications.


title: Tutorial: Complexity analysis of Singular Value Decomposition and its variants

abstract: We compared the regular Singular Value Decomposition (SVD), truncated SVD, Krylov method and Randomized PCA, in terms of time and space complexity. It is well-known that Krylov method and Randomized PCA only performs well when k << n, i.e. the number of eigenpair needed is far less than that of matrix size. We compared them for calculating all the eigenpairs. We also discussed the relationship between Principal Component Analysis and SVD.


title: Safe Contextual Bayesian Optimization for Sustainable Room Temperature PID Control Tuning

abstract: We tune one of the most common heating, ventilation, and air conditioning (HVAC) control loops, namely the temperature control of a room. For economical and environmental reasons, it is of prime importance to optimize the performance of this system. Buildings account from 20 to 40% of a country energy consumption, and almost 50% of it comes from HVAC systems. Scenario projections predict a 30% decrease in heating consumption by 2050 due to efficiency increase. Advanced control techniques can improve performance; however, the proportional-integral-derivative (PID) control is typically used due to its simplicity and overall performance. We use Safe Contextual Bayesian Optimization to optimize the PID parameters without human intervention. We reduce costs by 32% compared to the current PID controller setting while assuring safety and comfort to people in the room. The results of this work have an immediate impact on the room control loop performances and its related commissioning costs. Furthermore, this successful attempt paves the way for further use at different levels of HVAC systems, with promising energy, operational, and commissioning costs savings, and it is a practical demonstration of the positive effects that Artificial Intelligence can have on environmental sustainability.


title: ARMIN: Towards a More Efficient and Light-weight Recurrent Memory Network

abstract: In recent years, memory-augmented neural networks(MANNs) have shown promising power to enhance the memory ability of neural networks for sequential processing tasks. However, previous MANNs suffer from complex memory addressing mechanism, making them relatively hard to train and causing computational overheads. Moreover, many of them reuse the classical RNN structure such as LSTM for memory processing, causing inefficient exploitations of memory information. In this paper, we introduce a novel MANN, the Auto-addressing and Recurrent Memory Integrating Network (ARMIN) to address these issues. The ARMIN only utilizes hidden state ht for automatic memory addressing, and uses a novel RNN cell for refined integration of memory information. Empirical results on a variety of experiments demonstrate that the ARMIN is more light-weight and efficient compared to existing memory networks. Moreover, we demonstrate that the ARMIN can achieve much lower computational overhead than vanilla LSTM while keeping similar performances. Codes are available on github.com/zoharli/armin.


title: Non-user Inclusive Design for Maintaining Harmony of Real-Virtual Human Interaction in Augmented Reality

abstract: Augmented reality enables the illusion of contents such as objects and humans in the virtual world co-existing with users in the real world. However, non-users who are not aware of the presence of the virtual world and dynamically move nearby might either cause a conflict by directly breaking into space where a user is talking to a Virtual Human (VH), or be troubled when try to avoid disturbing the user. To maintain harmony and keep both the user's and non-users' comfort, we propose a method that controls the VH to adjust its own position to avoid such potential conflict. The difficulty to address this problem is that the agent must avoid potential conflict in a natural way to keep the user away from feeling unnatural. Our idea is to endow the VH with three capabilities: anticipating non-users walking around, understanding how to establish and maintain proper formation to adapt to the environment, and planning to avoid conflicts by shifting formation in advance. We develop a non-user inclusive spatial formation model that realizes natural arrangement shift corresponding to the environment based on theoretical sources from literature. We implemented our proposed model into a VH behavior planning system to achieve natural conflict avoidance. Evaluation experiments showed that it successfully reduces potential conflicts caused by non-users.


title: Uncovering the Semantics of Wikipedia Categories

abstract: The Wikipedia category graph serves as the taxonomic backbone for large-scale knowledge graphs like YAGO or Probase, and has been used extensively for tasks like entity disambiguation or semantic similarity estimation. Wikipedia's categories are a rich source of taxonomic as well as non-taxonomic information. The category 'German science fiction writers', for example, encodes the type of its resources (Writer), as well as their nationality (German) and genre (Science Fiction). Several approaches in the literature make use of fractions of this encoded information without exploiting its full potential. In this paper, we introduce an approach for the discovery of category axioms that uses information from the category network, category instances, and their lexicalisations. With DBpedia as background knowledge, we discover 703k axioms covering 502k of Wikipedia's categories and populate the DBpedia knowledge graph with additional 4.4M relation assertions and 3.3M type assertions at more than 87% and 90% precision, respectively.


title: Bifractal nature of chromosome contact maps

abstract: Modern biological techniques such as Hi-C permit to measure probabilities that different chromosomal regions are close in space. These probabilities can be visualised as matrices called contact maps. In this paper, we introduce a multifractal analysis of chromosomal contact maps. Our analysis reveals that Hi-C maps are bifractal, i.e. complex geometrical objects characterized by two distinct fractal dimensions. To rationalize this observation, we introduce a model that describes chromosomes as a hierarchical set of nested domains and we solve it exactly. The predicted multifractal spectrum is in excellent quantitative agreement with experimental data. Moreover, we show that our theory yields to a more robust estimation of the scaling exponent of the contact probability than existing methods. By applying this method to experimental data, we detect subtle conformational changes among chromosomes during differentiation of human stem cells.


title: Efficient Neural Interaction Function Search for Collaborative Filtering

abstract: In collaborative filtering (CF), interaction function (IFC) plays the important role of capturing interactions among items and users. The most popular IFC is the inner product, which has been successfully used in low-rank matrix factorization. However, interactions in real-world applications can be highly complex. Thus, other operations (such as plus and concatenation), which may potentially offer better performance, have been proposed. Nevertheless, it is still hard for existing IFCs to have consistently good performance across different application scenarios. Motivated by the recent success of automated machine learning (AutoML), we propose in this paper the search for simple neural interaction functions (SIF) in CF. By examining and generalizing existing CF approaches, an expressive SIF search space is designed and represented as a structured multi-layer perceptron. We propose an one-shot search algorithm that simultaneously updates both the architecture and learning parameters. Experimental results demonstrate that the proposed method can be much more efficient than popular AutoML approaches, can obtain much better prediction performance than state-of-the-art CF approaches, and can discover distinct IFCs for different data sets and tasks


title: Throughput Scaling of Covert Communication over Wireless Adhoc Networks

abstract: We consider the problem of covert communication over wireless adhoc networks in which (roughly) $n$ legitimate nodes (LNs) and $n^{\kappa}$ for $0<\kappa<1$ non-communicating warden nodes (WNs) are randomly distributed in a square of unit area. Each legitimate source wants to communicate with its intended destination node while ensuring that every WN is unable to detect the presence of the communication. In this scenario, we study the throughput scaling law. Due the covert communication constraint, the transmit powers are necessarily limited. Under this condition, we introduce a preservation region around each WN. This region serves to prevent transmission from the LNs and to increase the transmit power of the LNs outside the preservation regions. For the achievability results, multi-hop (MH), hierarchical cooperation (HC), and hybrid HC-MH schemes are utilized with some appropriate modifications. In the proposed MH and hybrid schemes, because the preservation regions may impede communication along direct data paths, the data paths are suitably modified by taking a detour around each preservation region. To avoid the concentration of detours resulting extra relaying burdens, we distribute the detours evenly over a wide region. In the proposed HC scheme, we control the symbol power and the scheduling of distributed multiple-input multiple-output transmission. We also present matching upper bounds on the throughput scaling under the assumption that every active LN consumes the same average transmit power over the time period in which the WNs observe the channel outputs.


title: A study of a nonlocal problem with Robin boundary conditions arising from MEMS technology

abstract: In the current work we study a nonlocal parabolic problem with Robin boundary conditions. The problem arises from the study of an idealized electrically actuated MEMS (Micro-Electro-Mechanical System) device. Initially we study the steady-state problem in one dimension and we present some numerical results regarding its bifurcation diagram. Next the $N-$dimensional nonlocal problem steady-state problem is investigated and a Pohozaev-type identity is first obtained which then facilitates derivation of an estimate of the pull-in voltage for radially symmetric domains. Later the time-dependent problem is investigated and global-in-time as well as quenching results are obtained again for general and radially symmetric domains. We close the current work with a numerical investigation of the presented nonlocal model via an adaptive numerical method. Various numerical experiments are presented verifying the obtained analytical results.


title: Feature ranking of active region source properties in solar flare forecasting and the uncompromised stochasticity of flare occurrence

abstract: Solar flares originate from magnetically active regions but not all solar active regions give rise to a flare. Therefore, the challenge of solar flare prediction benefits by an intelligent computational analysis of physics-based properties extracted from active region observables, most commonly line-of-sight or vector magnetograms of the active-region photosphere. For the purpose of flare forecasting, this study utilizes an unprecedented 171 flare-predictive active region properties, mainly inferred by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) in the course of the European Union Horizon 2020 FLARECAST project. Using two different supervised machine learning methods that allow feature ranking as a function of predictive capability, we show that: i) an objective training and testing process is paramount for the performance of every supervised machine learning method; ii) most properties include overlapping information and are therefore highly redundant for flare prediction; iii) solar flare prediction is still - and will likely remain - a predominantly probabilistic challenge.


title: Robustness Against Transactional Causal Consistency

abstract: Distributed storage systems and databases are widely used by various types of applications. Transactional access to these storage systems is an important abstraction allowing application programmers to consider blocks of actions (i.e., transactions) as executing atomically. For performance reasons, the consistency models implemented by modern databases are weaker than the standard serializability model, which corresponds to the atomicity abstraction of transactions executing over a sequentially consistent memory. Causal consistency for instance is one such model that is widely used in practice. In this paper, we investigate application-specific relationships between several variations of causal consistency and we address the issue of verifying automatically if a given transactional program is robust against causal consistency, i.e., all its behaviors when executed over an arbitrary causally consistent database are serializable. We show that programs without write-write races have the same set of behaviors under all these variations, and we show that checking robustness is polynomial time reducible to a state reachability problem in transactional programs over a sequentially consistent shared memory. A surprising corollary of the latter result is that causal consistency variations which admit incomparable sets of behaviors admit comparable sets of robust programs. This reduction also opens the door to leveraging existing methods and tools for the verification of concurrent programs (assuming sequential consistency) for reasoning about programs running over causally consistent databases. Furthermore, it allows to establish that the problem of checking robustness is decidable when the programs executed at different sites are finite-state.


title: Spin Logic Devices via Electric Field Controlled Magnetization Reversal by Spin-Orbit Torque

abstract: We describe a spin logic device with controllable magnetization switching of perpendicularly magnetized ferromagnet / heavy metal structures on a ferroelectric (1-x)[Pb(Mg1/3Nb2/3)O3]-x[PbTiO3] (PMN-PT) substrate using current-induced spin-orbit torque. The devices were operated without an external magnetic field and controlled by voltages as low as 10 V applied across the PMN-PT substrate, which is much lower compared to previous reports (500 V). The deterministic switching with smaller voltage was realized from the virgin state of the PMN-PT. Ferroelectric simulation shows the unsaturated minor loop exhibits obvious asymmetries in the polarizations. Larger polarization can be induced from the initial ferroelectric state, while it is difficult for opposite polarization. The XNOR, AND, NAND and NOT logic functions were demonstrated by the deterministic magnetization switching from the interaction between the spin-orbit torque and electric field at the PMN-PT/Pt interface. The nonvolatile spin logic scheme in this work is simple, scalable, programmable, which are favorable in the logic-in-memory design with low energy consumption.


title: Existence of quantum symmetries for graphs on up to seven vertices: a computer based approach

abstract: The symmetries of a finite graph are described by its automorphism group; in the setting of Woronowicz's quantum groups, a notion of a quantum automorphism group has been defined by Banica capturing the quantum symmetries of the graph. In general, there are more quantum symmetries than symmetries and it is a non-trivial task to determine when this is the case for a given graph: The question is whether or not the algebra associated to the quantum automorphism group is commutative. We use Gr\"obner base computations in order to tackle this problem; the implementation uses GAP and the SINGULAR package LETTERPLACE. We determine the existence of quantum symmetries for all connected, undirected graphs without multiple edges and without self-edges, for up to seven vertices. As an outcome, we infer within our regime that a classical automorphism group of order one or two is an obstruction for the existence of quantum symmetries.


title: Category-Theoretic Foundations of "STCLang: State Thread Composition as a Foundation for Monadic Dataflow Parallelism"

abstract: This manuscript gives a category-theoretic foundation to the composition of State Threads as a Foundation for Monadic Dataflow Parallelism. It serves as a supplementary formalization of the concepts introduced in the Article "STCLang: State Thread Composition as a Foundation for Monadic Dataflow Parallelism", as published in Proceedings of the 12th ACM SIGPLAN International Symposium on Haskell (Haskell'19).


title: Maximal Entanglement: Applications in Quantum Information and Particle Physics

abstract: In this PhD thesis, several aspects regarding maximal entanglement are analyzed. In the first chapter, Bell Inequalities are analyzed from an operational perspective as well as novel Bell inequalities are obtained together with their optimal settings for a maximal violation. Multipartite figures of merit, in particular, the hyperdeterminant, are the subject of the second chapter. They are applied to detect quantum phase transitions in several spin models. The third chapter focuses on the simulation of the XY model in a quantum computer. The quantum circuit obtained is tested in three current quantum devices. Quantum computers must be able to generate and hold highly entangled states in order to show a quantum advantage. This statement is analyzed in chapter four, where quantum circuits for maximally entangled states are presented. Finally, chapter five analyzes how maximal entanglement is generated at its fundamental level. Maximal entanglement constraints the QED vertex interaction and predicts a value for the weak mixing angle close to the experimental value.


title: Formulating causal questions and principled statistical answers

abstract: Although review papers on causal inference methods are now available, there is a lack of introductory overviews on what they can render and on the guiding criteria for choosing one particular method. This tutorial gives an overview in situations where an exposure of interest is set at a chosen baseline (`point exposure') and the target outcome arises at a later time point. We first phrase relevant causal questions and make a case for being specific about the possible exposure levels involved and the populations for which the question is relevant. Using the potential outcomes framework, we describe principled definitions of causal effects and of estimation approaches classified according to whether they invoke the no unmeasured confounding assumption (including outcome regression and propensity score-based methods) or an instrumental variable with added assumptions. We discuss challenges and potential pitfalls and illustrate application using a `simulation learner', that mimics the effect of various breastfeeding interventions on a child's later development. This involves a typical simulation component with generated exposure, covariate, and outcome data that mimic those from an observational or randomised intervention study. The simulation learner further generates various (linked) exposure types with a set of possible values per observation unit, from which observed as well as potential outcome data are generated. It thus provides true values of several causal effects. R code for data generation and analysis is available on www.ofcaus.org, where SAS and Stata code for analysis is also provided.


title: A sufficient criterion to determine planar self-Cheeger sets

abstract: We show a sufficient criterion to determine if a planar set $\Omega$ is a minimizer of the prescribed curvature functional among all of its subsets. As a special case, we derive a sufficient criterion to determine if $\Omega$ is a self-Cheeger set, i.e. if it minimizes the ratio $P(E)/|E|$ among all of its subsets. Specifically, if a Jordan domain $\Omega$ possesses the interior disk property of radius $|\Omega|/P(\Omega)$, then it is a self-Cheeger set; if it possesses the strict interior disk property then it is a minimal Cheeger set, i.e. the unique minimizer. As a side effect we provide a way to build self-Cheeger sets.


title: Borexino and General Neutrino Interactions

abstract: We derive constraints on all possible general neutrino-electron interactions (scalar, vector, pseudoscalar, axialvector and tensor) using the recent real time Borexino event rate measurements of $pp$, $pep$ and $^{7}Be$ solar neutrinos. The limits improve several previous ones from TEXONO and CHARM-II for incoming electron and muon neutrinos, and are the first ones for the tau flavor. Future improvements by next-generation solar neutrino experiments are also studied. The limits extend the physics reach of solar neutrino measurements to TeV-scale physics. Finally, the different properties of the new interactions for Dirac and Majorana neutrinos are discussed.


title: Sturmian ground states in classical lattice-gas models

abstract: We construct for the first time examples of non-frustrated, two-body, infinite-range, one-dimensional classical lattice-gas models without periodic ground-state configurations. Ground-state configurations of our models are Sturmian sequences defined by irrational rotations on the circle. We present minimal sets of forbidden patterns which define Sturmian sequences in a unique way. Our interactions assign positive energies to forbidden patterns and are equal to zero otherwise. We illustrate our construction by the well-known example of the Fibonacci sequences.

categories: math-ph cond-mat.stat-mech math.DS math.MP


title: Stimulated spin noise in an activated crystal

abstract: In the spin noise spectroscopy, the magnetic susceptibility spectrum is known to be provided by the spin-system untouched by any external perturbation, or, better to say, disturbed only by its thermal bath. We propose a new version of the spin noise spectroscopy, with the detected magnetization (Faraday-rotation) noise being stimulated by an external fluctuating magnetic field with a quasi-white spectrum. Experimental study of the stimulated spin noise performed on a $\text{BaF}_2 : \text{U}^{3+}$ crystal in a longitudinal magnetic field has revealed specific features of this approach and allowed us to identify the Van-Vleck and population-related contributions to the AC susceptibility of the system and to discover unusual magnetic-field dependence of the longitudinal spin relaxation rate in low magnetic fields. It is shown that spectra of the stimulated and spontaneous spin noise, being both closely related to the spin-system magnetic susceptibility, are still essentially different. Distinctions between the two types of the spin-noise spectra and two approaches to the spin noise spectroscopy are discussed.


title: Diffusion-influenced reactions on non-spherical partially absorbing axisymmetric surfaces

abstract: The calculation of the diffusion-controlled reaction rate for partially absorbing, non-spherical boundaries presents a formidable problem of broad relevance. In this paper we take the reference case of a spherical boundary and work out a perturbative approach to get a simple analytical formula for the first-order correction to the diffusive flux onto a non-spherical partially absorbing surface of revolution. To assess the range of validity of this formula, we derive exact and approximate expressions for the reaction rate in the case of partially absorbing prolate and oblate spheroids. We also present numerical solutions by a finite-element method that extend the validity analysis beyond spheroidal shapes. Our perturbative solution provides a handy way to quantify the effect of non-sphericity on the rate of capture in the general case of partial surface reactivity.


title: On the conditional distribution of the mean of the two closest among a set of three observations

abstract: Chemical analyses of raw materials are often repeated in duplicate or triplicate. The assay values obtained are then combined using a predetermined formula to obtain an estimate of the true value of the material of interest. When duplicate observations are obtained, their average typically serves as an estimate of the true value. On the other hand, the "best of three" method involves taking three measurements and using the average of the two closest ones as estimate of the true value. In this paper, we consider another method which potentially involves three measurements. Initially two measurements are obtained and if their difference is sufficiently small, their average is taken as estimate of the true value. However, if the difference is too large then a third independent measurement is obtained. The estimator is then defined as the average between the third observation and the one among the first two which is closest to it. Our focus in the paper is the conditional distribution of the estimate in cases where the initial difference is too large. We find that the conditional distributions are markedly different under the assumption of a normal distribution and a Laplace distribution.