title: Classical-to-Topological Transmission Line Couplers

abstract: Recent advances in topologically robust waveguiding for electromagnetic systems have presented opportunities for improving practical photonic and microwave devices. To bring this rich area of physics within the reach of application, it is critical for such systems to be interfaced with classical, continuous waveguiding and transmission line technology. This Letter presents a compact, highly efficient transition from a classical metallic transmission line to a topologically nontrivial line wave emulating the quantum spin Hall effect. A zero-gap antipodal slot line is used as the starting transmission line, which is then coupled to the topological metasurface via a field matching procedure. Additional modifications to the interface between the two structures to eliminate unwanted edge coupling improves transmission further. A simulated loss analysis isolates the effect of the transitions from the rest of the structure, showing a loss contribution of only 2.1% per classical-to-topological conversion. Using the transition, a quantitative characterization of the robustness of common topologically protected devices is presented. This design lays the foundation to integrate topologically robust metasurface transmission lines to traditional systems, opening the door to future uses of such structures in systems.


title: High-speed silicon microring modulator at 2-um waveband

abstract: We demonstrated a silicon integrated microring modulator working at the 2-um waveband with an L-shaped PN junction. 15-GHz 3-dB electro-optic bandwidth and <1 Vcm modulation efficiency for 45-Gbps NRZ-OOK signaling is achieved at 1960 nm.


title: Silicon microring synapses enable photonic deep learning beyond 9-bit precision

abstract: Deep neural networks (DNN) consist of layers of neurons interconnected by synaptic weights. A high bit-precision in weights is generally required to guarantee high accuracy in many applications. Minimizing error accumulation between layers is also essential when building large-scale networks. Recent demonstrations of photonic neural networks are limited in bit-precision due to crosstalk and the high sensitivity of optical components (e.g., resonators). Here, we experimentally demonstrate a record-high precision of 9 bits with a dithering control scheme for photonic synapses. We then numerically simulated the impact with increased synaptic precision on a wireless signal classification application. This work could help realize the potential of photonic neural networks for many practical, real-world tasks.

categories: physics.app-ph cs.SY eess.SY physics.optics


title: Distributional data analysis of accelerometer data from the NHANES database using nonparametric survey regression models

abstract: Accelerometers enable an objective measurement of physical activity levels among groups of individuals in free-living environments, providing high-resolution detail about physical activity changes at different time scales. Current approaches used in the literature for analyzing such data typically employ summary measures such as total inactivity time or compositional metrics. However, at the conceptual level, these methods have the potential disadvantage of discarding important information from recorded data when calculating these summaries and metrics since these typically depend on cut-offs related to exercise intensity zones chosen subjectively or even arbitrarily. Furthermore, much of the data collected in these studies follow complex survey designs. Then, using specific estimation strategies adapted to a particular sampling mechanism is mandatory. The aim of this paper is two-fold. First, a new functional representation of a distributional nature accelerometer data is introduced to build a complete individualized profile of each subject's physical activity levels. Second, we extend two nonparametric functional regression models, kernel smoothing and kernel ridge regression, to handle survey data and obtain reliable conclusions about the influence of physical activity in the different analyses performed in the complex sampling design NHANES cohort and so, show representation advantages.


title: Thermometric Calibration of the Ultrafast Relaxation Dynamics in Plasmonic Au Nanoparticles

abstract: The excitation of plasmonic nanoparticles by ultrashort laser pulses sets in motion a complex ultrafast relaxation process involving the gradual re-equilibration of the system's electron gas, lattice and environment. One of the major hurdles in studying these processes is the lack of direct measurements of the dynamic temperature evolution of the system subcomponents. We measured the dynamic optical response of ensembles of plasmonic Au nanoparticles following ultrashort-pulse excitation and we compared it with the corresponding static optical response as a function of the increasing temperature of the thermodynamic bath. Evaluating the two sets of data, the optical fingerprints of equilibrium or off-equilibrium responses could be clearly identified, allowing us to extract a dynamic thermometric calibration scale of the relaxation process, yielding the experimental ultrafast temperature evolution of the plasmonic particles as a function of time.


title: Tactile-RL for Insertion: Generalization to Objects of Unknown Geometry

abstract: Object insertion is a classic contact-rich manipulation task. The task remains challenging, especially when considering general objects of unknown geometry, which significantly limits the ability to understand the contact configuration between the object and the environment. We study the problem of aligning the object and environment with a tactile-based feedback insertion policy. The insertion process is modeled as an episodic policy that iterates between insertion attempts followed by pose corrections. We explore different mechanisms to learn such a policy based on Reinforcement Learning. The key contribution of this paper is to demonstrate that it is possible to learn a tactile insertion policy that generalizes across different object geometries, and an ablation study of the key design choices for the learning agent: 1) the type of learning scheme: supervised vs. reinforcement learning; 2) the type of learning schedule: unguided vs. curriculum learning; 3) the type of sensing modality: force/torque (F/T) vs. tactile; and 4) the type of tactile representation: tactile RGB vs. tactile flow. We show that the optimal configuration of the learning agent (RL + curriculum + tactile flow) exposed to 4 training objects yields an insertion policy that inserts 4 novel objects with over 85.0% success rate and within 3~4 attempts. Comparisons between F/T and tactile sensing, shows that while an F/T-based policy learns more efficiently, a tactile-based policy provides better generalization.


title: Quantum computing critical exponents

abstract: We show that the Variational Quantum-Classical Simulation algorithm admits a finite circuit depth scaling collapse when targeting the critical point of the transverse field Ising chain. The order parameter only collapses on one side of the transition due to a slowdown of the quantum algorithm when crossing the phase transition. In order to assess performance of the quantum algorithm and compute correlations in a system of up to 752 qubits, we use techniques from integrability to derive closed-form analytical expressions for expectation values with respect to the output of the quantum circuit. We also reduce a conjecture made by Ho and Hsieh about the exact preparation of the transverse field Ising ground state to a system of equations.


title: Revisiting Pulsar Velocities Using Gaia Data Release 2

abstract: Precise measurements of neutron star (NS) velocities provide critical clues to the supernova physics and evolution of binary systems. Based on Gaia Data Release 2 (DR2), we selected a sample of 24 young (<3 Myr) pulsars with precise parallax measurements and measured the velocity of their local standard of rest (LSR) and the velocity dispersion among their respective local stellar groups. The median velocity difference between thus calculated LSRs and the Galactic rotation model is ~7.6km/s, small compared to the typical velocity dispersion of ~27.5km/s. For pulsars off the Galactic plane, such differences grow significantly to as large as ~40 km/s. More importantly, the velocity dispersion of stars in the local group of low-velocity pulsars can be comparable to their transverse velocities, suggesting that the intrinsic velocities of NS progenitors should be taken into account when we consider their natal kicks and binary evolution. We also examined the double NS systems J0737- 3039A/B, and measured its transverse velocity to be 26(+18, -13) km/s assuming nearby Gaia sources being representative of its birth environment. This work demonstrated the feasibility and importance of using Gaia data to study the velocity of individual systems and velocity distribution of NSs.


title: Characterization of the dissipative mappings and their application to perturbations of dissipative-Hamiltonian systems

abstract: In this paper, we find necessary and sufficient conditions to identify pairs of matrices $X$ and $Y$ for which there exists $\Delta \in \mathbb C^{n,n}$ such that $\Delta+\Delta^*$ is positive semidefinite and $\Delta X=Y$. Such a $\Delta$ is called a dissipative mapping taking $X$ to $Y$. We also provide two different characterizations for the set of all dissipative mappings, and use them to characterize the unique dissipative mapping with minimal Frobenius norm. The minimal-norm dissipative mapping is then used to determine the distance to asymptotic instability for dissipative-Hamiltonian systems under general structure-preserving perturbations. We illustrate our results over some numerical examples and compare them with those of Mehl, Mehrmann and Sharma (Stability Radii for Linear Hamiltonian Systems with Dissipation Under Structure-Preserving Perturbations, SIAM J. Mat. Anal. Appl.\ 37 (4): 1625-1654, 2016).


title: Weighted homomorphisms on the $p$-analog of the Fourier-Stieltjes algebra induced by piecewise affine maps

abstract: In this paper, for $p\in(1,\infty)$ we study $p$-complete boundedness of weighted homomorphisms on the $p$-analog of the Fourier-Stieltjes algebras, $B_p(G)$, based on the $p$-operator space structure defined by the authors. Here, for a locally compact group $G$, the space $B_p(G)$ stands for Runde's definition of the $p$-analog of the Fourier-Stieltjes algebra and the implemented $p$-operator space structure is come from the duality between $B_p(G)$ and the algebra of universal $p$-pseudofunctions, $UPF_p(G)$. It is established that the homomorphism $\Phi_\alpha:B_p(G)\to B_p(H)$, defined by $\Phi(u)=u\circ\alpha$ on $Y$ and zero otherwise, is $p$-completely contractive when the continuous and proper map $\alpha :Y\subseteq H\to G$ is affine, and it is $p$-completely bounded whenever $\alpha$ is piecewise affine map. Moreover, we assume that $Y$ belongs to the coset ring generated by open and amenable subgroups of $H$. To obtain the result, by utilizing the properties of $QSL_p$-spaces and representations on them, the relation between $B_p(G/N)$ and a closed subalgebra of $B_p(G)$ is shown, where $N$ is a closed normal subgroup of $G$. Additionally, $p$-complete boundedness of several well-known maps on such algebras are obtained.


title: The XFaster Power Spectrum and Likelihood Estimator for the Analysis of Cosmic Microwave Background Maps

abstract: We present the XFaster analysis package. XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-$C_\ell$ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data, and also used as part of the Planck analysis. Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument. The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set.


title: Branch-and-cut algorithms for the covering salesman problem

abstract: The Covering Salesman Problem (CSP) is a generalization of the Traveling Salesman Problem in which the tour is not required to visit all vertices, as long as all vertices are covered by the tour. The objective of CSP is to find a minimum length Hamiltonian cycle over a subset of vertices that covers an undirected graph. In this paper, valid inequalities from the generalized traveling salesman problem are applied to the CSP in addition to new valid inequalities that explore distinct aspects of the problem. A branch-and-cut framework assembles exact and heuristic separation routines for integer and fractional CSP solutions. Computational experiments show that the proposed framework outperformed methodologies from literature with respect to optimality gaps. Moreover, optimal solutions were proven for several previously unsolved instances.


title: Probing the X-ray Soft Excess in the Narrow-Line Seyfert-1 (NLS1) Galaxy Tonantzintla S180

abstract: A detailed X-ray Soft Excess of the narrow-line Seyfert-1 galaxy Tonantzintla S180 (Ton s180) present. The study used four XMM-Newton observations for period from the year 2000 to the year 2016 taken from the XMM-Newton archive. Two different models were used to treat the X-ray Soft Excess for each data separately, the first model is power law component and the second one is two black body components. We found that, for all observations, the spectra fitted by using a power law component are poor and the output parameters much different from the previous studies, this model unsuccessfully reproduces the continuum shape of the spectrum, including the soft excess, of all observations (Nh = 3.62 +-1020 cm^{-2} and photon index for the hard and soft bands are 3.31 and 1.57). The study proved that it can rely upon the two black body components in treating X-ray Soft Excess of the narrow-line Seyfert-1 galaxy (Nh= 1.395 +- 10^20 cm^{-2} and The temperature of the two black-body components (KT) are 0.075 and 0.17 keV).


title: Direct laser writing for cardiac tissue engineering: a microfluidic heart on a chip with integrated transducers

abstract: We have designed and fabricated a microfluidic-based platform for sensing mechanical forces generated by cardiac microtissues in a highly-controlled microenvironment. Our fabrication approach combines Direct Laser Writing (DLW) lithography with soft lithography. At the center of our platform is a cylindrical volume, divided into two chambers by a cylindrical polydimethylsiloxane (PDMS) shell. Cells are seeded into the inner chamber from a top opening, and the microtissue assembles onto tailor-made attachment sites on the inner walls of the cylindrical shell. The outer chamber is electrically and fluidically isolated from the inner one by the cylindrical shell and is designed for actuation and sensing purposes. Externally applied pressure waves to the outer chamber deform parts of the cylindrical shell and thus allow us to exert time-dependent forces on the microtissue. Oscillatory forces generated by the microtissue similarly deform the cylindrical shell and change the volume of the outer chamber, resulting in measurable electrical conductance changes. We have used this platform to study the response of cardiac microtissues derived from human induced pluripotent stem cells (hiPSC) under prescribed mechanical loading and pacing.


title: Trends in eBusiness and eGovernment

abstract: The first chapter is a critical review and a case study in eBusiness, with special attention to the digital currencies resource and its possibilities. 2. chapter attempts to incorporate the UTAUT model with perceived risk theory to explore its impact on the intention to use m-government services. 3. chapter aims to assess the level of gender inclusivity in the municipal e-procurement processes in the City of Johannesburg as a case study. It uses a GAD approach. 4. chapter examines the impediments that derail the intensive uptake of eLearning programmes in a particular higher education institution. The study adopted an inductive research paradigm that followed a qualitative research strategy. Data were collected by means of one-on-one in-depth interviews from selected faculty members at a nominated institution of higher learning. 5. chapter investigated the role of KMS in enhancing the export performance of firms operating within the manufacturing sector in Zimbabwe. The study used a quantitative approach in which a survey questionnaire was distributed to 555 managers drawn from 185 manufacturing firms based in Harare. Data analyses involved the use of descriptive statistics, Spearman correlations and regression analysis. In the sixth chapter, a survey was undertaken on 131 SMEs from the Pelagonija region in order to determine the current level of SME digitalization within the region. It is aimed to compare with the EU average and to make conclusions on the impact of the SME digitalization on region GDP growth as well as revenues collection. The last chapter s purpose was to develop a measuring and modelling framework, an instrument of IBSQ for the South African banking sector. Snowball and convenience sampling, both non-probability techniques were used to recruit participants for the study. A total of 310 Internet banking customer responses were utilised in the analysis.


title: How Powerful are Performance Predictors in Neural Architecture Search?

abstract: Early methods in the rapidly developing field of neural architecture search (NAS) required fully training thousands of neural networks. To reduce this extreme computational cost, dozens of techniques have since been proposed to predict the final performance of neural architectures. Despite the success of such performance prediction methods, it is not well-understood how different families of techniques compare to one another, due to the lack of an agreed-upon evaluation metric and optimization for different constraints on the initialization time and query time. In this work, we give the first large-scale study of performance predictors by analyzing 31 techniques ranging from learning curve extrapolation, to weight-sharing, to supervised learning, to "zero-cost" proxies. We test a number of correlation- and rank-based performance measures in a variety of settings, as well as the ability of each technique to speed up predictor-based NAS frameworks. Our results act as recommendations for the best predictors to use in different settings, and we show that certain families of predictors can be combined to achieve even better predictive power, opening up promising research directions. Our code, featuring a library of 31 performance predictors, is available at https://github.com/automl/naslib.


title: Topological charge density waves at half-integer filling of a moir\'e superlattice

abstract: At partial filling of a flat band, strong electronic interactions may favor gapped states harboring emergent topology with quantized Hall conductivity. Emergent topological states have been found in partially filled Landau levels and Hofstadter bands; in both cases, a large magnetic field is required to engineer the underlying flat band. The recent observation of quantum anomalous Hall effects (QAH) in narrow band moir\'e systems has led to the theoretical prediction that such phases may be realized even at zero magnetic field. Here we report the experimental observation of insulators with Chern number $C=1$ in the zero magnetic field limit at $\nu=3/2$ and $7/2$ filling of the moir\'e superlattice unit cell in twisted monolayer-bilayer graphene (tMBG). Our observation of Chern insulators at half-integer values of $\nu$ suggests spontaneous doubling of the superlattice unit cell, in addition to spin- and valley-ferromagnetism. This is confirmed by Hartree-Fock calculations, which find a topological charge density wave ground state at half filling of the underlying $C=2$ band, in which the Berry curvature is evenly partitioned between occupied and unoccupied states. We find the translation symmetry breaking order parameter is evenly distributed across the entire folded superlattice Brillouin zone, suggesting that the system is in the flat band, strongly correlated limit. Our findings show that the interplay of quantum geometry and Coulomb interactions in moir\'e bands allows for topological phases at fractional superlattice filling that spontaneously break time-reversal symmetry, a prerequisite in pursuit of zero magnetic field phases harboring fractional statistics as elementary excitations or bound to lattice dislocations.


title: Some Exact Results in QCD-like Theories

abstract: I propose a controlled approximation to QCD-like theories with massless quarks by employing supersymmetric QCD perturbed by anomaly-mediated supersymmetry breaking. They have identical massless particle contents. Thanks to the ultraviolet-insensitivity of anomaly mediation, dynamics can be worked out exactly when $m \ll \Lambda$, where $m$ is the size of supersymmetry breaking and $\Lambda$ the dynamical scale of the gauge theory. I demonstrate that chiral symmetry is dynamically broken for $N_{f} \leq \frac{3}{2} N_{c}$ while the theories lead to non-trivial infrared fixed points for larger number of flavors. While there may be a phase transition as $m$ is increased beyond $\Lambda$, qualitative agreements with expectations in QCD are encouraging and suggest that two limits $m \ll \Lambda$ and $m \gg \Lambda$ may be in the same universality class.


title: Realizing topologically ordered states on a quantum processor

abstract: The discovery of topological order has revolutionized the understanding of quantum matter in modern physics and provided the theoretical foundation for many quantum error correcting codes. Realizing topologically ordered states has proven to be extremely challenging in both condensed matter and synthetic quantum systems. Here, we prepare the ground state of the toric code Hamiltonian using an efficient quantum circuit on a superconducting quantum processor. We measure a topological entanglement entropy near the expected value of $\ln2$, and simulate anyon interferometry to extract the braiding statistics of the emergent excitations. Furthermore, we investigate key aspects of the surface code, including logical state injection and the decay of the non-local order parameter. Our results demonstrate the potential for quantum processors to provide key insights into topological quantum matter and quantum error correction.


title: Survival times of supramassive neutron stars resulting from binary neutron star mergers

abstract: A binary neutron star (BNS) merger can lead to various outcomes, from indefinitely stable neutron stars, through supramassive (SMNS) or hypermassive (HMNS) neutron stars supported only temporarily against gravity, to black holes formed promptly after the merger. Up-to-date constraints on the BNS total mass and the neutron star equation of state suggest that a long-lived SMNS may form in $\sim 0.45-0.9$ of BNS mergers. A maximally rotating SMNS needs to lose $\sim 3-6\times 10^{52}$ erg of it's rotational energy before it collapses, on a fraction of the spin-down timescale. A SMNS formation imprints on the electromagnetic counterparts to the BNS merger. However, a comparison with observations reveals tensions. First, the distribution of collapse times is too wide and that of released energies too narrow (and the energy itself too large) to explain the observed distributions of internal X-ray plateaus, invoked as evidence for SMNS-powered energy injection. Secondly, the immense energy injection into the blastwave should lead to extremely bright radio transients which previous studies found to be inconsistent with deep radio observations of short gamma-ray bursts. Furthermore, we show that upcoming all-sky radio surveys will constrain the extracted energy distribution, independently of a GRB jet formation. Our results can be self-consistently understood, provided that most BNS merger remnants collapse shortly after formation (even if their masses are low enough to allow for SMNS formation). This naturally occurs if the remnant retains half or less of its initial energy by the time it enters solid body rotation.


title: How Robust are the Inferred Density and Metallicity of the Circumgalactic Medium?

abstract: Quantitative estimates of the basic properties of the circumgalactic medium(CGM), such as its density and metallicity, depend on the spectrum of incident UV background radiation. Models of UV background are known to have large variations, mainly because they are synthesized using poorly constrained parameters, which introduce uncertainty in the inferred properties of the CGM. Here, we quantify this uncertainty using a large set of new UV background models with physically motivated toy models of metal-enriched CGM. We find that, the inferred density and metallicity of low-density ($10^{-5}$ cm$^{-3}$) gas is uncertain by factors of 6.3 and 3.2, whereas high density ($10^{-3}$ cm$^{-3}$) gas by factors of 4 and 1.6, respectively. The variation in the shape of the UV background models is entirely responsible for such a variation in the metallicity while variation in the density arises from both normalization and shape of the UV background. Moreover, we find a harder (softer) UV background infers higher (lower) density and metallicity. We also study warm-hot gas at $T= 10^{5.5}$ K and find that metallicity is robustly estimated but the inferred density is uncertain by a factor of 3 to 5.4 for low to high-density gas. Such large uncertainties in density and metallicity may severely limit the studies of the CGM and demand better observational constraints on the input parameters used in synthesizing UV background.


title: The GOGREEN survey: the internal dynamics of clusters of galaxies at redshift 0.9-1.4

abstract: We aim to determine the mass, velocity anisotropy, and pseudo phase-space density profiles (M(r), beta(r), and Q(r), respectively) of clusters of galaxies at the highest redshifts investigated in detail so far. We combine the GOGREEN and GCLASS spectroscopic data-sets for 14 clusters with mass M200 > 10^14 Msolar at redshifts 0.9 < z < 1.4. We stack these 14 clusters into an ensemble cluster of 581 member galaxies with stellar mass > 10^9.5 M_solar. We use the MAMPOSSt method and the inversion of the Jeans equation technique to determine M(r) and beta(r). We then combine the results of the M(r) and beta(r) analysis to determine Q(r) for the ensemble cluster. The concentration c200 of the ensemble cluster M(r) is in excellent agreement with predictions from LambdaCDM cosmological numerical simulations, and with previous determinations for clusters of similar mass and at similar redshifts, obtained from gravitational lensing and X-ray data. We see no significant difference between the total mass density and either the galaxy number density distributions or the stellar mass distribution. Star-forming galaxies are spatially significantly less concentrated than quiescent galaxies. The orbits of cluster galaxies are isotropic near the center and more radial outside. Star-forming galaxies and galaxies of low stellar mass tend to move on more radially elongated orbits than quiescent galaxies and galaxies of high stellar mass. Q(r), determined either using the total mass or the number density profile, is very close to the power-law behavior predicted by numerical simulations. The internal dynamics of clusters at the highest redshift probed in detail so far are very similar to those of lower-redshift clusters, and in excellent agreement with predictions of numerical simulations. The clusters in our sample have already reached a high degree of dynamical relaxation. (Abridged)


title: An effective matrix model for dynamical end of the world branes in Jackiw-Teitelboim gravity

abstract: We study Jackiw-Teitelboim gravity with dynamical end of the world branes in asymptotically nearly AdS$_2$ spacetimes. We quantize this theory in Lorentz signature, and compute the Euclidean path integral summing over topologies including dynamical branes. The latter will be seen to exactly match with a modification of the SSS matrix model. The resolution of UV divergences in the gravitational instantons involving the branes will lead us to understand the matrix model interpretation of the Wilsonian effective theory perspective on the gravitational theory. We complete this modified SSS matrix model nonperturbatively by extending the integration contour of eigenvalues into the complex plane. Furthermore, we give a new interpretation of other phases in such matrix models. We derive an effective $W(\Phi)$ dilaton gravity, which exhibits similar physics semiclassically. In the limit of a large number of flavors of branes, the effective extremal entropy $S_{0,\text{eff}}$ has the form of counting the states of these branes.


title: Position-Dependent Correlation Function of Weak Lensing Convergence

abstract: We provide a systematic study of the position-dependent correlation function in weak lensing convergence maps and its relation to the squeezed limit of the three-point correlation function (3PCF) using state-of-the-art numerical simulations. We relate the position-dependent correlation function to its harmonic counterpart, i.e., the position-dependent power spectrum or equivalently the integrated bispectrum. We use a recently proposed improved fitting function, BiHalofit, for the bispectrum to compute the theoretical predictions as a function of source redshifts. In addition to low redshift results ($z_s=1.0-2.0$), we also provide results for maps inferred from lensing of the cosmic microwave background, i.e., $z_s=1100$. We include a {\em Euclid}-type realistic survey mask and noise. In agreement with the recent studies on the position-dependent power spectrum, we find that the results from simulations are consistent with the theoretical expectations when appropriate corrections are included. Performing a rough estimate, we find that the (S/N) for the detection of the position-dependent correlation function from {\em Euclid}-type mask with $f_{sky}=0.35$, can range between $6-12$ depending on the value of the intrinsic ellipticity distribution parameter $\sigma_{\epsilon} = 0.3-1.0$. For reconstructed $\kappa$ maps using an ideal CMB survey the (S/N) $\approx 1.8$. We also found that a $10\%$ deviation in $\sigma_8$ can be detected using IB for the optimistic case of $\sigma_\epsilon=0.3$ with a (S/N) $\approx 5$. The (S/N) for such detection in case of $\Omega_M$ is lower.


title: Bijections from Dyck and Motzkin meanders with catastrophes to pattern avoiding Dyck paths

abstract: In this note, we present constructive bijections from Dyck and Motzkin meanders with catastrophes to Dyck paths avoiding some patterns. As a byproduct, we deduce correspondences from Dyck and Motzkin excursions to restricted Dyck paths.


title: Bispectral dual Hahn polynomials with an arbitrary number of continuous parameters

abstract: We construct new examples of bispectral dual Hahn polynomials, i.e., orthogonal polynomials with respect to certain superposition of Christoffel and Geronimus transforms of the dual Hahn measure and which are also eigenfunctions of a higher order difference operator. The new examples have the novelty that they depend on an arbitrary number of continuous parameters. These are the first examples with this property constructed from the classical discrete families.


title: Scan Specific Artifact Reduction in K-space (SPARK) Neural Networks Synergize with Physics-based Reconstruction to Accelerate MRI

abstract: Purpose: To develop a scan-specific model that estimates and corrects k-space errors made when reconstructing accelerated Magnetic Resonance Imaging (MRI) data. Methods: Scan-Specific Artifact Reduction in k-space (SPARK) trains a convolutional-neural-network to estimate and correct k-space errors made by an input reconstruction technique by back-propagating from the mean-squared-error loss between an auto-calibration signal (ACS) and the input technique's reconstructed ACS. First, SPARK is applied to GRAPPA and demonstrates improved robustness over other scan-specific models, such as RAKI and residual-RAKI. Subsequent experiments demonstrate that SPARK synergizes with residual-RAKI to improve reconstruction performance. SPARK also improves reconstruction quality when applied to advanced acquisition and reconstruction techniques like 2D virtual coil (VC-) GRAPPA, 2D LORAKS, 3D GRAPPA without an integrated ACS region, and 2D/3D wave-encoded images. Results: SPARK yields 1.5x - 2x RMSE reduction when applied to GRAPPA and improves robustness to ACS size for various acceleration rates in comparison to other scan-specific techniques. When applied to advanced reconstruction techniques such as residual-RAKI, 2D VC-GRAPPA and LORAKS, SPARK achieves up to 20% RMSE improvement. SPARK with 3D GRAPPA also improves performance by ~2x and perceived image quality without a fully sampled ACS region. Finally, SPARK synergizes with non-cartesian 2D and 3D wave-encoding imaging by reducing RMSE between 20-25% and providing qualitative improvements. Conclusion: SPARK synergizes with physics-based acquisition and reconstruction techniques to improve accelerated MRI by training scan-specific models to estimate and correct reconstruction errors in k-space.


title: Proving Non-termination by Program Reversal

abstract: We present a new approach to proving non-termination of non-deterministic integer programs. Our technique is rather simple but efficient. It relies on a purely syntactic reversal of the program's transition system followed by a constraint-based invariant synthesis with constraints coming from both the original and the reversed transition system. The latter task is performed by a simple call to an off-the-shelf SMT-solver, which allows us to leverage the latest advances in SMT-solving. Moreover, our method offers a combination of features not present (as a whole) in previous approaches: it handles programs with non-determinism, provides relative completeness guarantees and supports programs with polynomial arithmetic. The experiments performed with our prototype tool RevTerm show that our approach, despite its simplicity and stronger theoretical guarantees, is at least on par with the state-of-the-art tools, often achieving a non-trivial improvement under a proper configuration of its parameters.


title: grASP: A Graph Based ASP-Solver and Justification System

abstract: Answer set programming (ASP) is a popular nonmonotonic-logic based paradigm for knowledge representation and solving combinatorial problems. Computing the answer set of an ASP program is NP-hard in general, and researchers have been investing significant effort to speed it up. The majority of current ASP solvers employ SAT solver-like technology to find these answer sets. As a result, justification for why a literal is in the answer set is hard to produce. There are dependency graph based approaches to find answer sets, but due to the representational limitations of dependency graphs, such approaches are limited. We propose a novel dependency graph-based approach for finding answer sets in which conjunction of goals is explicitly represented as a node which allows arbitrary answer set programs to be uniformly represented. Our representation preserves causal relationships allowing for justification for each literal in the answer set to be elegantly found. Performance results from an implementation are also reported. Our work paves the way for computing answer sets without grounding a program.


title: The Edges of Our Universe

abstract: This paper explores the fundamental causal limits on how much of the universe we can observe or affect. It distinguishes four principal regions: the affectable universe, the observable universe, the eventually observable universe, and the ultimately observable universe. It then shows how these (and other) causal limits set physical bounds on what spacefaring civilisations could achieve over the longterm future.


title: Energy-saving Cross-layer Optimization of Big Data Transfer Based on Historical Log Analysis

abstract: With the proliferation of data movement across the Internet, global data traffic per year has already exceeded the Zettabyte scale. The network infrastructure and end-systems facilitating the vast data movement consume an extensive amount of electricity, measured in terawatt-hours per year. This massive energy footprint costs the world economy billions of dollars partially due to energy consumed at the network end-systems. Although extensive research has been done on managing power consumption within the core networking infrastructure, there is little research on reducing the power consumption at the end-systems during active data transfers. This paper presents a novel cross-layer optimization framework, called Cross-LayerHLA, to minimize energy consumption at the end-systems by applying machine learning techniques to historical transfer logs and extracting the hidden relationships between different parameters affecting both the performance and resource utilization. It utilizes offline analysis to improve online learning and dynamic tuning of application-level and kernel-level parameters with minimal overhead. This approach minimizes end-system energy consumption and maximizes data transfer throughput. Our experimental results show that Cross-LayerHLA outperforms other state-of-the-art solutions in this area.


title: Learning Description Logic Ontologies. Five Approaches. Where Do They Stand?

abstract: The quest for acquiring a formal representation of the knowledge of a domain of interest has attracted researchers with various backgrounds into a diverse field called ontology learning. We highlight classical machine learning and data mining approaches that have been proposed for (semi-)automating the creation of description logic (DL) ontologies. These are based on association rule mining, formal concept analysis, inductive logic programming, computational learning theory, and neural networks. We provide an overview of each approach and how it has been adapted for dealing with DL ontologies. Finally, we discuss the benefits and limitations of each of them for learning DL ontologies.


title: Physics Informed Convex Artificial Neural Networks (PICANNs) for Optimal Transport based Density Estimation

abstract: Optimal Mass Transport (OMT) is a well studied problem with a variety of applications in a diverse set of fields ranging from Physics to Computer Vision and in particular Statistics and Data Science. Since the original formulation of Monge in 1781 significant theoretical progress been made on the existence, uniqueness and properties of the optimal transport maps. The actual numerical computation of the transport maps, particularly in high dimensions, remains a challenging problem. By Brenier's theorem, the continuous OMT problem can be reduced to that of solving a non-linear PDE of Monge-Ampere type whose solution is a convex function. In this paper, building on recent developments of input convex neural networks and physics informed neural networks for solving PDE's, we propose a Deep Learning approach to solve the continuous OMT problem. To demonstrate the versatility of our framework we focus on the ubiquitous density estimation and generative modeling tasks in statistics and machine learning. Finally as an example we show how our framework can be incorporated with an autoencoder to estimate an effective probabilistic generative model.


title: Solving Homotopy Domain Equations

abstract: In order to get $\lambda$-models with a rich structure of $\infty$-groupoid, which we call "homotopy $\lambda$-models", a general technique is described for solving domain equations on any cartesian closed $\infty$-category (c.c.i.) with enough points. Finally, the technique is applied in a particular c.c.i., where some examples of homotopy $\lambda$-models are given.


title: Two-Stage Gauss--Seidel Preconditioners and Smoothers for Krylov Solvers on a GPU cluster

abstract: Gauss-Seidel (GS) relaxation is often employed as a preconditioner for a Krylov solver or as a smoother for Algebraic Multigrid (AMG). However, the requisite sparse triangular solve is difficult to parallelize on many-core architectures such as graphics processing units (GPUs). In the present study, the performance of the traditional GS relaxation based on a triangular solve is compared with two-stage variants, replacing the direct triangular solve with a fixed number of inner Jacobi-Richardson (JR) iterations. When a small number of inner iterations is sufficient to maintain the Krylov convergence rate, the two-stage GS (GS2) often outperforms the traditional algorithm on many-core architectures. We also compare GS2 with JR. When they perform the same number of flops for SpMV (e.g. three JR sweeps compared to two GS sweeps with one inner JR sweep), the GS2 iterations, and the Krylov solver preconditioned with GS2, may converge faster than the JR iterations. Moreover, for some problems (e.g. elasticity), it was found that JR may diverge with a damping factor of one, whereas two-stage GS may improve the convergence with more inner iterations. Finally, to study the performance of the two-stage smoother and preconditioner for a practical problem, %(e.g. using tuned damping factors), these were applied to incompressible fluid flow simulations on GPUs.


title: Knowledge, Trust, Security and Covertness In Massively Distributed Social Platforms: An Epistemic Networks Approach

abstract: Social networks seeded, crystallized and structured by large-scale social platforms like Facebook, LinkedIn, Slack, Twitter, WhatsApp, Instagram and WeChat exhibit complicated epistemic structures and dynamics that arise from the nature of interactive knowledge that users/participants possess, share and shape. What users know about what other users know, about what other users know they know, and so forth, plays an important role in channeling user behavior on the platform as well as the structuring and dynamics of the social network the platform engenders. We use a novel approach to characterize knowledge states and flows in social networks, based on epistemic networks, or epinets, to study the relationship between the structure of platforms and the structure of social networks that form on them.


title: Beyond Short Clips: End-to-End Video-Level Learning with Collaborative Memories

abstract: The standard way of training video models entails sampling at each iteration a single clip from a video and optimizing the clip prediction with respect to the video-level label. We argue that a single clip may not have enough temporal coverage to exhibit the label to recognize, since video datasets are often weakly labeled with categorical information but without dense temporal annotations. Furthermore, optimizing the model over brief clips impedes its ability to learn long-term temporal dependencies. To overcome these limitations, we introduce a collaborative memory mechanism that encodes information across multiple sampled clips of a video at each training iteration. This enables the learning of long-range dependencies beyond a single clip. We explore different design choices for the collaborative memory to ease the optimization difficulties. Our proposed framework is end-to-end trainable and significantly improves the accuracy of video classification at a negligible computational overhead. Through extensive experiments, we demonstrate that our framework generalizes to different video architectures and tasks, outperforming the state of the art on both action recognition (e.g., Kinetics-400 & 700, Charades, Something-Something-V1) and action detection (e.g., AVA v2.1 & v2.2).


title: Synergies between Vera C. Rubin Observatory, Nancy Grace Roman Space Telescope, and Euclid Mission: Constraining Dark Energy with Type Ia Supernovae

abstract: We review the needs of the supernova community for improvements in survey coordination and data sharing that would significantly boost the constraints on dark energy using samples of Type Ia supernovae from the Vera C. Rubin Observatories, the \textit{Nancy Grace Roman Space Telescope}, and the \textit{Euclid} Mission. We discuss improvements to both statistical and systematic precision that the combination of observations from these experiments will enable. For example, coordination will result in improved photometric calibration, redshift measurements, as well as supernova distances. We also discuss what teams and plans should be put in place now to start preparing for these combined data sets. Specifically, we request coordinated efforts in field selection and survey operations, photometric calibration, spectroscopic follow-up, pixel-level processing, and computing. These efforts will benefit not only experiments with Type Ia supernovae, but all time-domain studies, and cosmology with multi-messenger astrophysics.


title: Exploring the structure of $^{29}$Ne

abstract: We apply a fully quantum mechanical Coulomb breakup theory under the aegis of post form finite-range distorted wave Born approximation to analyze the elastic Coulomb breakup of $^{29}$Ne on $^{208}$Pb at $244$\,MeV/u. We calculate several reaction observables to quantify its structural parameters. One-neutron removal cross-section is calculated to check the consistency of the ground state configuration of $^{29}$Ne with the available experimental data. A scrutiny of the parallel momentum distribution of the charged fragment reveals a full width at half maximum of $82$\,MeV/c, which is in good agreement with the experimental value and indicates a moderate halo for a nearly spherical $^{29}$Ne in the $^{28}$Ne$(0^+) \otimes 2p_{3/2}\nu$ ground state. The energy-angular distributions and average momentum of the charged fragment point to the absence of post-acceleration effects in the breakup process, a desirable result for the elastic breakup.


title: Site-dependent selection of atoms for homogeneous atom-cavity coupling

abstract: We demonstrate a method to obtain homogeneous atom-cavity coupling by selecting and keeping $^{87}$Rb atoms that are near maximally coupled to the cavity's standing-wave mode. We select atoms by imposing an AC Stark shift on the ground state hyperfine microwave transition frequency with light injected into the cavity. We then induce a spin flip with microwaves that are resonant for atoms that are near maximally coupled to the cavity mode of interest, after which, we use radiation pressure forces to remove from the cavity all the atoms in the initial spin state. Achieving greater homogeneity in the atom-cavity coupling will potentially enhance entanglement generation, intracavity driving of atomic transitions, cavity-optomechanics, and quantum simulations. This approach can easily be extended to other atomic species with microwave or optical transitions.


title: Monitoring of Particle Flux and LET Variations with Pulse Stretching Inverters

abstract: This work investigates the use of pulse stretching inverters for monitoring the variation of flux and Linear Energy Transfer (LET) of energetic particles. The basic particle detector consists of two cascaded pulse stretching (skew-sized) inverters designed in CMOS technology, and the required sensing area is obtained by connecting multiple two-inverter pulse stretching cells in parallel, and employing the required number of parallel arrays. The particle strikes are detected in terms of the Single Event Transients (SETs), and the detector provides the information on the SET count rate and SET pulse width variation, from which the particle flux and LET can be determined. The main advantage of the proposed solution is the possibility to sense the LET variations using purely digital processing logic. The SPICE simulations done on IHP 130 nm bulk CMOS technology have shown that the SET pulse width at the output of detector changes by 550 ps over the LET range from 1 to 100 MeVcm2mg-1. The proposed solution is intended to operate as an on-chip particle detector within the self-adaptive multiprocessing systems.


title: Universal features in panarthropod inter-limb coordination during forward walking

abstract: Terrestrial animals must often negotiate heterogeneous, varying environments. Accordingly, their locomotive strategies must adapt to a wide range of terrain, as well as to a range of speeds in order to accomplish different behavioral goals. Studies in \textit{Drosophila} have found that inter-leg coordination patterns (ICPs) vary smoothly with walking speed, rather than switching between distinct gaits as in vertebrates (e.g., horses transitioning between trotting and galloping). Such a continuum of stepping patterns implies that separate neural controllers are not necessary for each observed ICP. Furthermore, the spectrum of \textit{Drosophila} stepping patterns includes all canonical coordination patterns observed during forward walking in insects. This raises the exciting possibility that the controller in \textit{Drosophila} is common to all insects, and perhaps more generally to panarthropod walkers. Here, we survey and collate data on leg kinematics and inter-leg coordination relationships during forward walking in a range of arthropod species, as well as include data from a recent behavioral investigation into the tardigrade \textit{Hypsibius exemplaris}. Using this comparative dataset, we point to several functional and morphological features that are shared amongst panarthropods. The goal of the framework presented in this review is to emphasize the importance of comparative functional and morphological analyses in understanding the origins and diversification of walking in Panarthropoda.


title: Sharp bounds of Hankel determinants of second and third order for inverse functions of certain class of univalent functions

abstract: Let ${\mathcal A}$ be the class of functions that are analytic in the unit disc ${\mathbb D}$, normalized such that $f(z)=z+\sum_{n=2}^\infty a_nz^n$, and let class ${\mathcal U}(\lambda)$, $0<\lambda\le1$, consists of functions $f\in{\mathcal A}$, such that \[ \left |\left (\frac{z}{f(z)} \right )^{2}f'(z)-1\right | < \lambda\quad (z\in {\mathbb D}). \] In this paper we determine the sharp upper bounds for the Hankel determinants of second and third order for the inverse functions of functions from the class ${\mathcal U}(\lambda)$.


title: Demonstration of Shor encoding on a trapped-ion quantum computer

abstract: Fault-tolerant quantum error correction (QEC) is crucial for unlocking the true power of quantum computers. QEC codes use multiple physical qubits to encode a logical qubit, which is protected against errors at the physical qubit level. Here we use a trapped ion system to experimentally prepare $m$-qubit GHZ states and sample the measurement results to construct $m\times m$ logical states of the $[[m^2,1,m]]$ Shor code, up to $m=7$. The synthetic logical fidelity shows how deeper encoding can compensate for additional gate errors in state preparation for larger logical states. However, the optimal code size depends on the physical error rate and we find that $m=5$ has the best performance in our system. We further realize the direct logical encoding of the $[[9,1,3]]$ Shor code on nine qubits in a thirteen-ion chain for comparison, with $98.8(1)\%$ and $98.5(1)\%$ fidelity for state $\left\vert\pm\right\rangle_L$, respectively.


title: A Review of Particle Detectors for Space-Borne Self-Adaptive Fault-Tolerant Systems

abstract: The soft error rate (SER) of integrated circuits (ICs) operating in space environment may vary by several orders of magnitude due to the variable intensity of radiation exposure. To ensure the radiation hardness without compromising the system performance, it is necessary to implement the dynamic hardening mechanisms which can be activated under the critical radiation exposure. Such operating scenario requires the real-time detection of energetic particles responsible for the soft errors. Although numerous particle detection solutions have been reported, very few works address the on-chip particle detectors suited for the self-adaptive fault tolerant microprocessor systems for space missions. This work reviews the state-of-the-art particle detectors, with emphasis on two solutions for the self-adaptive systems: particle detector based on embedded SRAM and particle detector based on pulse stretching inverters.


title: Type Prediction Systems

abstract: Inferring semantic types for entity mentions within text documents is an important asset for many downstream NLP tasks, such as Semantic Role Labelling, Entity Disambiguation, Knowledge Base Question Answering, etc. Prior works have mostly focused on supervised solutions that generally operate on relatively small-to-medium-sized type systems. In this work, we describe two systems aimed at predicting type information for the following two tasks, namely, a TypeSuggest module, an unsupervised system designed to predict types for a set of user-entered query terms, and an Answer Type prediction module, that provides a solution for the task of determining the correct type of the answer expected to a given query. Our systems generalize to arbitrary type systems of any sizes, thereby making it a highly appealing solution to extract type information at any granularity.


title: Wet Extreme Mass Ratio Inspirals May Be More Common For Spaceborne Gravitational Wave Detection

abstract: Extreme Mass Ratio Inspirals (EMRIs) can be classified as dry EMRIs and wet EMRIs based on their formation mechanisms. Dry (or the "loss-cone") EMRIs, previously considered as the main EMRI sources for the Laser Interferometer Space Antenna, are primarily produced by multi-body scattering in the nuclear star cluster and gravitational capture. In this work, we highlight an alternative EMRI formation channel: (wet) EMRI formation assisted by the accretion flow around accreting galactic-center massive black holes (MBHs). In this channel, the accretion disk captures stellar-mass black holes that are intially moving on inclined orbits, and subsequently drives them to migrate towards the MBH - this process boosts the formation rate of EMRIs in such galaxies by orders of magnitude. Taking into account the fraction ($\mathcal O(10^{-2}-10^{-1})$) of active galactic nuclei where the MBHs are expected to be rapidly accreting, we forecast that wet EMRIs will contribute an important or even dominant fraction of all detectable EMRIs by spaceborne gravitational wave detectors.


title: Distance-duality in theories with a nonminimal coupling to gravity

abstract: We show that gravitational theories with a nonminimal coupling (NMC) to the matter fields lead to a violation of Etherington's distance-duality relation, which relates the luminosity and angular diameter distances. We derive constraints on power-law and exponential NMC models using existing measurements of type Ia supernovae and baryon acoustic oscillations throughout the redshift range $0<z<1.5$. These complement previous constrains derived from cosmic-microwave background radiation and big-bang nucleosynthesis data.


title: An efficient 146-line 3D sensitivity analysis code of stress-based topology optimization written in MATLAB

abstract: This paper presents an efficient and compact MATLAB code for three-dimensional stress-based sensitivity analysis. The 146 lines code includes the finite element analysis and p-norm stress sensitivity analysis based on the adjoint method. The 3D sensitivity analysis for p-norm global stress measure is derived and explained in detail accompanied by corresponding MATLAB code. The correctness of the analytical sensitivity is verified by comparison with finite difference approximation. The nonlinear optimization solver is chosen as the Method of moving asymptotes (MMA). Three typical volume-constrained stress minimization problems are presented to verify the effectiveness of sensitivity analysis code. The MATLAB code presented in this paper can be extended to resolve different stress related 3D topology optimization problems. The complete program for sensitivity analysis is given in the Appendix and is intended for educational purposes only.


title: Convergence of limit shapes for 2D near-critical first-passage percolation

abstract: We consider Bernoulli first-passage percolation on the triangular lattice in which sites have 0 and 1 passage times with probability $p$ and $1-p$, respectively. For each $p\in(0,p_c)$, let $\mathcal {B}(p)$ be the limit shape in the classical "shape theorem", and let $L(p)$ be the correlation length. We show that as $p\uparrow p_c$, the rescaled limit shape $L(p)^{-1}\mathcal {B}(p)$ converges to a Euclidean disk. This improves a result of Chayes et al. [J. Stat. Phys. 45 (1986) 933--951]. The proof relies on the scaling limit of near-critical percolation established by Garban et al. [J. Eur. Math. Soc. 20 (2018) 1195--1268], and uses the construction of the collection of continuum clusters in the scaling limit introduced by Camia et al. [Springer Proceedings in Mathematics \& Statistics, 299 (2019) 44--89].


title: Estimation Technique for a Contact Point Between two Materials in a Stationary Heat Transfer Problem

abstract: An inverse problem for a stationary heat transfer process is studied for a totally isolated bar on its lateral surface, of negligible diameter, made up of two consecutive sections of different, isotropic and homogeneous materials. At the left boundary, a Dirichlet type condition is imposed that represents a constant temperature source while a Robin type condition that models the heat dissipation by convection is considered at the right one. Many articles in the literature focus on thermal and stress analysis at the interface but no one is dedicated to the estimation of the contact point location between the two materials. In this work, it is assumed that the interface position is unknown. A technique to determine it from a unique noisy flow measurement at the right boundary is introduced. Necessary and sufficient conditions are derived in order to obtain the estimation of the interface point from a heat flux measured at the right boundary. Numerical solutions are obtained together with an expression for the estimation error. Moreover, an elasticity analysis is included to study the influence of data errors. The results show that our approach is useful for determining the location of the materials interface.


title: Optical trapping with structured light

abstract: Optical trapping describes the interaction between light and matter to manipulate micro-objects through momentum transfer. In the case of 3D trapping with a single beam, this is termed optical tweezers. Optical tweezers are a powerful and non-invasive tool for manipulating small objects, which have become indispensable in many fields, including physics, biology, soft condensed matter, amongst others. In the early days, optical trapping were typically used with a single Gaussian beam. In recent years, we have witnessed the rapid progress in the use of structured light beams with customized phase, amplitude and polarization in optical trapping. Unusual beam properties, such as phase singularities on-axis, propagation invariant nature, have opened up novel capabilities to the study of micromanipulation in liquid, air and vacuum. In this review, we summarize the recent advances in the field of optical trapping using structured light beams.


title: Modeling Censored Mobility Demand through Quantile Regression Neural Networks

abstract: Shared mobility services require accurate demand models for effective service planning. On the one hand, modeling the full probability distribution of demand is advantageous because the entire uncertainty structure preserves valuable information for decision-making. On the other hand, demand is often observed through the usage of the service itself, so that the observations are censored, as they are inherently limited by available supply. Since the 1980s, various works on Censored Quantile Regression models have performed well under such conditions. Further, in the last two decades, several papers have proposed to implement these models flexibly through Neural Networks. However, the models in current works estimate the quantiles individually, thus incurring a computational overhead and ignoring valuable relationships between the quantiles. We address this gap by extending current Censored Quantile Regression models to learn multiple quantiles at once and apply these to synthetic baseline datasets and datasets from two shared mobility providers in the Copenhagen metropolitan area in Denmark. The results show that our extended models yield fewer quantile crossings and less computational overhead without compromising model performance.


title: The Coronavirus is a Bioweapon: Analysing Coronavirus Fact-Checked Stories

abstract: The 2020 coronavirus pandemic has heightened the need to flag coronavirus-related misinformation, and fact-checking groups have taken to verifying misinformation on the Internet. We explore stories reported by fact-checking groups PolitiFact, Poynter and Snopes from January to June 2020, characterising them into six story clusters before then analyse time-series and story validity trends and the level of agreement across sites. We further break down the story clusters into more granular story types by proposing a unique automated method with a BERT classifier, which can be used to classify diverse story sources, in both fact-checked stories and tweets.


title: Linear differential operators with polynomial coefficients generating generalised Sylvester-Kac matrices

abstract: A method of generating differential operators is used to solve the spectral problem for a generalisation of the Sylvester-Kac matrix. As a by-product, we find a linear differential operator with polynomial coefficients of the first order that has a finite sequence of polynomial eigenfunctions generalising the operator considered by M. Kac. In addition, we explain spectral properties of two related tridiagonal matrices whose shape differ from our generalisation.


title: Uncertainty-Aware Annotation Protocol to Evaluate Deformable Registration Algorithms

abstract: Landmark correspondences are a widely used type of gold standard in image registration. However, the manual placement of corresponding points is subject to high inter-user variability in the chosen annotated locations and in the interpretation of visual ambiguities. In this paper, we introduce a principled strategy for the construction of a gold standard in deformable registration. Our framework: (i) iteratively suggests the most informative location to annotate next, taking into account its redundancy with previous annotations; (ii) extends traditional pointwise annotations by accounting for the spatial uncertainty of each annotation, which can either be directly specified by the user, or aggregated from pointwise annotations from multiple experts; and (iii) naturally provides a new strategy for the evaluation of deformable registration algorithms. Our approach is validated on four different registration tasks. The experimental results show the efficacy of suggesting annotations according to their informativeness, and an improved capacity to assess the quality of the outputs of registration algorithms. In addition, our approach yields, from sparse annotations only, a dense visualization of the errors made by a registration method. The source code of our approach supporting both 2D and 3D data is publicly available at https://github.com/LoicPeter/evaluation-deformable-registration.


title: Saturation bounds for smooth varieties

abstract: We prove bounds on the saturation degrees of homogeneous ideals (and their powers) defining smooth complex projective varieties. For example, we show that a classical statement due to Macualay for zero-dimensional complete intersection ideals holds for any smooth variety. For curves, we bound the saturation degree of powers in terms of the regularity.


title: Simulating near-field enhancement in transmission of airborne viruses with a quadrature-based model

abstract: Airborne viruses, such as influenza, tuberculosis, and SARS-CoV-2, are transmitted through virus-laden particles expelled when an infectious person sneezes, coughs, talks, or breathes. These virus-laden particles are more highly concentrated in the expiratory jet of an infectious person than in a well-mixed room, but this near-field enhancement in virion exposure has not been well quantified. Transmission of airborne viruses depends on factors that are inherently variable and, in many cases, poorly constrained, and quantifying this uncertainty requires large ensembles of model simulations that span the variability in input parameters. However, models that are well-suited to simulate the near-field evolution of respiratory particles are also computationally expensive, which limits the exploration of parametric uncertainty. In order to perform many simulations that span the wide variability in factors governing transmission, we developed the Quadrature-based model of Respiratory Aerosol and Droplets (QuaRAD). QuaRAD is an efficient framework for simulating the evolution of virus-laden particles after they are expelled from an infectious person, their deposition to the nasal cavity of a susceptible person, and the subsequent risk of initial infection. We simulated 10,000 scenarios to quantify the risk of initial infection by a particular virus, SARS-CoV-2. The predicted risk of infection was highly variable among scenarios and, in each scenario, was strongly enhanced near the infectious individual. In more than 50% of scenarios, the physical distancing needed to avoid near-field enhancements in airborne transmission was beyond the recommended safe distance of two meters (six feet) if the infectious person is not wearing a mask, though this distance defining the near-field extent was also highly variable among scenarios.


title: Monodromy Defects in Free Field Theories

abstract: We study co-dimension two monodromy defects in theories of conformally coupled scalars and free Dirac fermions in arbitrary $d$ dimensions. We characterise this family of conformal defects by computing the one-point functions of the stress-tensor and conserved current for Abelian flavour symmetries as well as two-point functions of the displacement operator. In the case of $d=4$, the normalisation of these correlation functions are related to defect Weyl anomaly coefficients, and thus provide crucial information about the defect conformal field theory. We provide explicit checks on the values of the defect central charges by calculating the universal part of the defect contribution to entanglement entropy, and further, we use our results to extract the universal part of the vacuum R\'enyi entropy. Moreover, we leverage the non-supersymmetric free field results to compute a novel defect Weyl anomaly coefficient in a $d=4$ theory of free $\mathcal{N}=2$ hypermultiplets. Including singular modes in the defect operator product expansion of fundamental fields, we identify notable relevant deformations in the singular defect theories and show that they trigger a renormalisation group flow towards an IR fixed point with the most regular defect OPE. We also study Gukov-Witten defects in free $d=4$ Maxwell theory and show that their central charges vanish.


title: Channel Estimation for Intelligent Reflecting Surface Assisted Wireless Communications

abstract: In this paper, the minimum mean square error (MMSE) channel estimation for intelligent reflecting surface (IRS) assisted wireless communication systems is investigated. In the considered setting, each row vector of the equivalent channel matrix from the base station (BS) to the users is shown to be Bessel $K$ distributed, and all these row vectors are independent of each other. By introducing a Gaussian scale mixture model, we obtain a closed-form expression for the MMSE estimate of the equivalent channel, and determine analytical upper and lower bounds on the mean square error. Using the central limit theorem, we conduct an asymptotic analysis of the MMSE estimate, and show that the upper bound on the mean square error of the MMSE estimate is equal to the asymptotic mean square error of the MMSE estimation when the number of reflecting elements at the IRS tends to infinity. Numerical simulations show that the gap between the upper and lower bounds are very small, and they almost overlap with each other at medium signal-to-noise ratio (SNR) levels and moderate number of elements at the IRS.


title: Sensing force and charge at the nanoscale with a single-molecule tether

abstract: Measuring the electrophoretic mobility of molecules is a powerful experimental approach for investigating biomolecular processes. A frequent challenge in the context of single-particle measurements is throughput, limiting the obtainable statistics. Here, we present a molecular force sensor and charge detector based on parallelised imaging and tracking of tethered double-stranded DNA functionalised with charged nanoparticles interacting with an externally applied electric field. Tracking the position of the tethered particle with simultaneous nanometre precision and microsecond temporal resolution allows us to detect and quantify electrophoretic forces down to the sub-piconewton scale. Furthermore, we demonstrate that this approach is capable of detecting changes to the particle charge state, as induced by the addition of charged biomolecules or changes to pH. Our approach provides an alternative route to studying structural and charge dynamics at the single-molecule level.


title: Obstruction Flat Rigidity of the CR 3-Sphere

abstract: On a bounded strictly pseudoconvex domain in $\mathbb{C}^n$, $n >1$, the smoothness of the Cheng-Yau solution to Fefferman's complex Monge-Amp\`ere equation up to the boundary is obstructed by a local curvature invariant of the boundary, the CR obstruction density $\mathcal{O}$. While local examples of obstruction flat CR manifolds are plentiful, the only known compact examples are the spherical CR manifolds. We consider the obstruction flatness problem for small deformations of the standard CR 3-sphere. That rigidity holds for the CR sphere was previously known (in all dimensions) for the case of embeddable CR structures, where it also holds at the infinitesimal level. In the 3-dimensional case, however, a CR structure need not be embeddable. While in the nonembeddable case we may no longer interpret the obstruction density $\mathcal{O}$ in terms of the boundary regularity of Fefferman's equation (or the logarithmic singularity of the Bergman kernel) the equation $\mathcal{O}\equiv 0$ is still of great interest, e.g., since it corresponds to the Bach flat equation of conformal gravity for the Fefferman space of the CR structure (a conformal Lorentzian 4-manifold). Unlike in the embeddable case, it turns out that in the nonembeddable case there is an infinite dimensional space of solutions to the linearized obstruction flatness equation on the standard CR 3-sphere and this space defines a natural complement to the tangent space of the embeddable deformations. In spite of this, we show that the CR 3-sphere does not admit nontrivial obstruction flat deformations, embeddable or nonembeddable.


title: High pCO$_2$ reduces sensitivity to CO$_2$ perturbations on temperate, Earth-like planets throughout most of habitable zone

abstract: The nearly logarithmic radiative impact of CO$_2$ means that planets near the outer edge of the liquid water habitable zone (HZ) require $\sim$10$^6$x more CO$_2$ to maintain temperatures conducive to standing liquid water on the planetary surface than their counterparts near the inner edge. This logarithmic radiative response also means that atmospheric CO$_2$ changes of a given mass will have smaller temperature effects on higher pCO$_2$ planets. Ocean pH is linked to atmospheric pCO$_2$ through seawater carbonate speciation and calcium carbonate dissolution/precipitation, and the response of pH to changes in pCO$_2$ also decreases at higher initial pCO$_2$. Here, we use idealized climate and ocean chemistry models to demonstrate that CO$_2$ perturbations large enough to cause catastrophic changes to surface temperature and ocean pH on low-pCO$_2$ planets in the innermost region of the HZ are likely to have much smaller effects on planets with higher pCO$_2$. Major bouts of extraterrestrial fossil fuel combustion or volcanic CO$_2$ outgassing on high-pCO$_2$ planets in the mid-to-outer HZ should have mild or negligible impacts on surface temperature and ocean pH. Owing to low pCO$_2$, Phanerozoic Earth's surface environment may be unusually volatile compared to similar planets receiving lower instellation.


title: Interlayer exciton mediated second harmonic generation in bilayer MoS2

abstract: Second harmonic generation (SHG) is a non-linear optical process, where two photons coherently combine into one photon of twice their energy. Efficient SHG occurs for crystals with broken inversion symmetry, such as transition metal dichalcogenide monolayers. Here we show tuning of non-linear optical processes in an inversion symmetric crystal. This tunability is based on the unique properties of bilayer MoS2, that shows strong optical oscillator strength for the intra- but also inter-layer exciton resonances. As we tune the SHG signal onto these resonances by varying the laser energy, the SHG amplitude is enhanced by several orders of magnitude. In the resonant case the bilayer SHG signal reaches amplitudes comparable to the off-resonant signal from a monolayer. In applied electric fields the interlayer exciton energies can be tuned due to their in-built electric dipole via the Stark effect. As a result the interlayer exciton degeneracy is lifted and the bilayer SHG response is further enhanced by an additional two orders of magnitude, well reproduced by our model calculations.


title: Liger for Next Generation Keck AO: Filter Wheel and Pupil Design

abstract: Liger is a next-generation near-infrared imager and integral field spectrograph (IFS) for the W.M. Keck Observatory designed to take advantage of the Keck All-Sky Precision Adaptive Optics (KAPA) upgrade. Liger will operate at spectral resolving powers between R$\sim$4,000 - 10,000 over a wavelength range of 0.8-2.4$\mu$m. Liger takes advantage of a sequential imager and spectrograph design that allows for simultaneous observations between the two channels using the same filter wheel and cold pupil stop. We present the design for the filter wheels and pupil mask and their location and tolerances in the optical design. The filter mechanism is a multi-wheel design drawing from the heritage of the current Keck/OSIRIS imager single wheel design. The Liger multi-wheel configuration is designed to allow future upgrades to the number and range of filters throughout the life of the instrument. The pupil mechanism is designed to be similarly upgradeable with the option to add multiple pupil mask options. A smaller wheel mechanism allows the user to select the desired pupil mask with open slots being designed in for future upgrade capabilities. An ideal pupil would match the shape of the image formed of the primary and would track its rotation. For different pupil shapes without tracking we model the additional exposure time needed to achieve the same signal to noise of an ideal pupil and determine that a set of fixed masks of different shapes provides a mechanically simpler system with little compromise in performance.


title: MetricNet: Towards Improved Modeling For Non-Intrusive Speech Quality Assessment

abstract: The objective speech quality assessment is usually conducted by comparing received speech signal with its clean reference, while human beings are capable of evaluating the speech quality without any reference, such as in the mean opinion score (MOS) tests. Non-intrusive speech quality assessment has attracted much attention recently due to the lack of access to clean reference signals for objective evaluations in real scenarios. In this paper, we propose a novel non-intrusive speech quality measurement model, MetricNet, which leverages label distribution learning and joint speech reconstruction learning to achieve significantly improved performance compared to the existing non-intrusive speech quality measurement models. We demonstrate that the proposed approach yields promisingly high correlation to the intrusive objective evaluation of speech quality on clean, noisy and processed speech data.


title: Liger for Next Generation Keck Adaptive Optics: Opto-Mechanical Dewar for Imaging Camera and Slicer

abstract: Liger is a next generation adaptive optics (AO) fed integral field spectrograph (IFS) and imager for the W. M. Keck Observatory. This new instrument is being designed to take advantage of the upgraded AO system provided by Keck All-Sky Precision Adaptive-optics (KAPA). Liger will provide higher spectral resolving power (R$\sim$4,000-10,000), wider wavelength coverage ($\sim$0.8-2.4 $\mu$m), and larger fields of view than any current IFS. We present the design and analysis for a custom-made dewar chamber for characterizing the Liger opto-mechanical system. This dewar chamber is designed to test and assemble the Liger imaging camera and slicer IFS components while being adaptable for future experiments. The vacuum chamber will operate below $10^{-5}$ Torr with a cold shield that will be kept below 90 K. The dewar test chamber will be mounted to an optical vibration isolation platform and further isolated from the cryogenic and vacuum systems with bellows. The cold head and vacuums will be mounted to a custom cart that will also house the electronics and computer that interface with the experiment. This test chamber will provide an efficient means of calibrating and characterizing the Liger instrument and performing future experiments.


title: Sampling discretization of the uniform norm

abstract: Discretization of the uniform norm of functions from a given finite dimensional subspace of continuous functions is studied. We pay special attention to the case of trigonometric polynomials with frequencies from an arbitrary finite set with fixed cardinality. We give two different proofs of the fact that for any $N$-dimensional subspace of the space of continuous functions it is sufficient to use $e^{CN}$ sample points for an accurate upper bound for the uniform norm. Previous known results show that one cannot improve on the exponential growth of the number of sampling points for a good discretization theorem in the uniform norm. Also, we prove a general result, which connects the upper bound on the number of sampling points in the discretization theorem for the uniform norm with the best $m$-term bilinear approximation of the Dirichlet kernel associated with the given subspace. We illustrate application of our technique on the example of trigonometric polynomials.


title: Surface-Normal Free-Space Beam Projection via Slow-Light Standing Wave Resonances in Extra-Large Near-Zero Index Grating Couplers

abstract: On-chip grating couplers directly connect photonic circuits to free-space light. The commonly used photonic gratings have been specialized for small areas, specific intensity profiles and non-vertical beam projection. This falls short of the precise and flexible wavefront control over large beam areas needed to empower emerging integrated miniaturized optical systems that leverage volumetric light matter interactions, including trapping, cooling, and interrogation of atoms, bio- and chemi- sensing and complex free-space interconnect. The large coupler size challenges general inverse design techniques, and solutions obtained by them are often difficult to physically understand and generalize. Here by posing the problem to a carefully constrained computational inverse design algorithm capable of large area structures, we discover a qualitatively new class of grating couplers. The numerically found solutions can be understood as coupling an incident photonic slab mode to a spatially extended slow-light (near-zero refractive index) region, backed by a Bragg reflector. The structure forms a spectrally broad standing wave resonance at the target wavelength, radiating vertically into free space. A reflection-less adiabatic transition critically couples the incident photonic mode to the resonance, and the numerically optimized lower cladding provides 70 $\%$ overall theoretical conversion efficiency. We have experimentally validated efficient surface normal collimated emission of $\approx$ 90 $\mu$m full width at half maximum Gaussian at the thermally tunable operating wavelength of $\approx$ 780 nm. The variable-mesh-deformation inverse design approach scales to extra-large photonic devices, while directly implementing the fabrication constraints. The deliberate choice of smooth parametrization resulted in a novel type of solution, which is both efficient and physically comprehensible.


title: Diverse Gaussian Noise Consistency Regularization for Robustness and Uncertainty Calibration

abstract: Deep neural networks achieve high prediction accuracy when the train and test distributions coincide. In practice though, various types of corruptions occur which deviate from this setup and cause severe performance degradations. Few methods have been proposed to address generalization in the presence of unforeseen domain shifts. In particular, digital noise corruptions arise commonly in practice during the image acquisition stage and present a significant challenge for current methods. In this paper, we propose a diverse Gaussian noise consistency regularization method for improving robustness of image classifiers under a variety of corruptions while still maintaining high clean accuracy. We derive bounds to motivate and understand the behavior of our Gaussian noise consistency regularization using a local loss landscape analysis. Our approach improves robustness against unforeseen noise corruptions by 4.2-18.4% over adversarial training and other strong diverse data augmentation baselines across several benchmarks. Furthermore, it improves robustness and uncertainty calibration by 3.7% and 5.5%, respectively, against all common corruptions (weather, digital, blur, noise) when combined with state-of-the-art diverse data augmentations.


title: Calibration of the Halpha Age-Activity relation for M dwarfs

abstract: In this work, we calibrate the relationship between Halpha emission and M dwarf ages. We compile a sample of 892 M dwarfs with Halpha equivalent width (HaEW) measurements from the literature that are either co-moving with a white dwarf of known age (21 stars) or in a known young association (871 stars). In this sample we identify 7 M dwarfs that are new candidate members of known associations. By dividing the stars into active and inactive categories according to their HaEW and spectral type (SpT), we find that the fraction of active dwarfs decreases with increasing age, and the form of the decline depends on SpT. Using the compiled sample of age-calibrators we find that HaEW and fractional Halpha luminosity (LHaLbol) decrease with increasing age. HaEW for SpT<M7 decreases gradually up until ~1Gyr. For older ages, we found only two early M dwarfs which are both inactive and seem to continue the gradual decrease. We also found 14 mid-type out of which 11 are inactive and present a significant decrease of HaEW, suggesting that the magnetic activity decreases rapidly after ~1Gyr. We fit LHaLbol versus age with a broken power-law and find an index of -0.11+0.02-0.01 for ages <~776Myr. The index becomes much steeper at older ages however a lack of field age-calibrators leaves this part of the relation far less constrained. Finally, from repeated independent measurements for the same stars we find that 94% of these has a level of HaEW variability <=5A at young ages (<1Gyr).


title: FBCNet: A Multi-view Convolutional Neural Network for Brain-Computer Interface

abstract: Lack of adequate training samples and noisy high-dimensional features are key challenges faced by Motor Imagery (MI) decoding algorithms for electroencephalogram (EEG) based Brain-Computer Interface (BCI). To address these challenges, inspired from neuro-physiological signatures of MI, this paper proposes a novel Filter-Bank Convolutional Network (FBCNet) for MI classification. FBCNet employs a multi-view data representation followed by spatial filtering to extract spectro-spatially discriminative features. This multistage approach enables efficient training of the network even when limited training data is available. More significantly, in FBCNet, we propose a novel Variance layer that effectively aggregates the EEG time-domain information. With this design, we compare FBCNet with state-of-the-art (SOTA) BCI algorithm on four MI datasets: The BCI competition IV dataset 2a (BCIC-IV-2a), the OpenBMI dataset, and two large datasets from chronic stroke patients. The results show that, by achieving 76.20% 4-class classification accuracy, FBCNet sets a new SOTA for BCIC-IV-2a dataset. On the other three datasets, FBCNet yields up to 8% higher binary classification accuracies. Additionally, using explainable AI techniques we present one of the first reports about the differences in discriminative EEG features between healthy subjects and stroke patients. Also, the FBCNet source code is available at https://github.com/ravikiran-mane/FBCNet.


title: Observation of two independent skyrmion phases in a chiral magnetic material

abstract: Magnetic materials can host skyrmions, which are topologically non-trivial spin textures. In chiral magnets with cubic lattice symmetry, all previously-observed skyrmion phases require thermal fluctuations to become thermodynamically stable in bulk materials, and therefore exist only at relatively high temperature, close to the helimagnetic transition temperature. Other stabilization mechanisms require a lowering of the cubic crystal symmetry. Here, we report the identification of a second skyrmion phase in Cu$_{2}$OSeO$_{3}$ at low temperature and in the presence of an applied magnetic field. The new skyrmion phase is thermodynamically disconnected from the well-known, nearly-isotropic, high-temperature phase, and exists, in contrast, when the external magnetic field is oriented along the $\langle100\rangle$ crystal axis only. Theoretical modelling provides evidence that the stabilization mechanism is given by well-known cubic anisotropy terms, and accounts for an additional observation of metastable helices tilted away from the applied field. The identification of two distinct skyrmion phases in the same material and the generic character of the underlying mechanism suggest a new avenue for the discovery, design, and manipulation of topological spin textures.


title: Reynolds number dependence of Lyapunov exponents of turbulence and fluid particles

abstract: The Navier-Stokes equations generate an infinite set of generalized Lyapunov exponents defined by different ways of measuring the distance between exponentially diverging perturbed and unperturbed solutions. This set is demonstrated to be similar, yet different, from the generalized Lyapunov exponent that provides moments of distance between two fluid particles below the Kolmogorov scale. We derive rigorous upper bounds on dimensionless Lyapunov exponent of the fluid particles that demonstrate the exponent's decay with Reynolds number $Re$ in accord with previous studies. In contrast, terms of cumulant series for exponents of the moments have power-law growth with $Re$. We demonstrate as an application that the growth of small fluctuations of magnetic field in ideal conducting turbulence is hyper-intermittent, being exponential in both time and Reynolds number. We resolve the existing contradiction between the theory, that predicts slow decrease of dimensionless Lyapunov exponent of turbulence with $Re$, and observations exhibiting quite fast growth. We demonstrate that it is highly plausible that a pointwise limit for the growth of small perturbations of the Navier-Stokes equations exists.

categories: physics.flu-dyn astro-ph.GA cond-mat.stat-mech nlin.CD


title: Revisiting Flat bands and localization

abstract: Flat bands imply lack of itinerancy due to some constraints that, in principle, results in anomalous behaviors with randomness. By a molecular orbital (MO) representation of the flat band systems, random MO models are introduced where the degeneracy due to the flat bands is preserved even with randomness. The zero modes of the chiral symmetric system with sublattice imbalance belong to the class. After explaining the generic flat band construction by MOs, several examples are discussed with numerical demonstration as sawtooth lattice in one dimension and hyper-Pyrochlore lattice in any $d$-dimensions that extends the Kagome ($d=2$) and Pyrochlore ($d=3$) lattices to general dimensions.


title: Numerical simulation of baffled circulating fluidized bed with Geldart B particles

abstract: As most classical fluidization technologies, Circulating Fluidized Beds (CFB) are widely employed in chemical and physical operations. However, CFB Reactors still have some disadvantages such as the non-uniform solid particles distribution and the backflow of the solid particles near the wall, which deteriorate the system mixing, and lead to low system conversion/heat and mass transfer rate. To improve the CFB performance, we propose in the current work to study the effect of the ring baffle configuration on the fluidization system hydrodynamics using a multifluid Eulerian CFD model, incorporating the Kinetic Theory of Granular Flow. In this approach, two different types of ring baffles, circular and trapezoidal, were considered and the corresponding results were compared. The results revealed that the incorporation of ring baffles improved the system mixing and reduced the backflow near the wall. Nevertheless, the shape of baffles had a limited impact on the system hydrodynamics.


title: Investigating the Reliability in Three RAID Storage Models and Effect of Ordering Replicas on Disks

abstract: One of the most important parts of cloud computing is storage devices, and Redundant Array of Independent Disks (RAID) systems are well known and frequently used storage devices. With the increasing production of data in cloud environments, we need high-reliable storage, given the importance of data. RAID system's reliability analysis is of particular significance in the area of cloud storage. Generally, data redundancy is used to create fault tolerance and increases the reliability of storage. This study has considered three examples of the simple RAID storage models and analyzed their reliability. All of which have a replication factor of two and have the same number of disks and reliabilities. The only difference is the model that they used for generating the redundancy. To compare these three models' reliability, we examined the degree of fault tolerance (FT) and calculated the models' reliability using the reliability block diagram (RBD). In this paper, the effect of redundancy's type and the blocks' arrangement on the system's reliability was investigated.


title: ICFP 2020 Post-Conference Report

abstract: This document describes the ICFP 2020 virtual conference, including the planning process and the criteria that informed its design, plus feedback from the post-conference survey. It is intended to provide a record of the event and give advice to future organizers of virtual conferences.


title: Long-Term Orbit Dynamics of Decommissioned Geostationary Satellites

abstract: In nominal mission scenarios, geostationary satellites perform end-of-life orbit maneuvers to reach suitable disposal orbits, where they do not interfere with operational satellites. This research investigates the long-term orbit evolution of decommissioned geostationary satellite under the assumption that the disposal maneuver does not occur and the orbit evolves with no control. The dynamical model accounts for all the relevant harmonics of the gravity field at the altitude of geostationary orbits, as well as solar radiation pressure and third-body perturbations caused by the Moon and the Sun. Orbit propagations are performed using two algorithms based on different equations of motion and numerical integration methods: (i) Gauss planetary equations for modified equinoctial elements with a Runge-Kutta numerical integration scheme based on 8-7th-order Dorman and Prince formulas; (ii) Cartesian state equations of motion in an Earth-fixed frame with a Runge-Kutta Fehlberg 7/8 integration scheme. The numerical results exhibit excellent agreement over integration times of decades. Some well-known phenomena emerge, such as the longitudinal drift due to the resonance between the orbital motion and Earth's rotation, attributable to the J22 term of the geopotential. In addition, the third-body perturbation due to Sun and Moon causes two major effects: (a) a precession of the orbital plane, and (b) complex longitudinal dynamics. This study proposes an analytical approach for the prediction of the precessional motion and shows its agreement with the orbit evolution obtained numerically. Moreover, long-term orbit propagations show that the above mentioned complex longitudinal dynamics persists over time scales of several decades. Frequent and unpredictable migrations toward different longitude regions occur, in contrast with the known effects due only to the J22 perturbation.


title: TreeToaster: Towards an IVM-Optimized Compiler

abstract: A compiler's optimizer operates over abstract syntax trees (ASTs), continuously applying rewrite rules to replace subtrees of the AST with more efficient ones. Especially on large source repositories, even simply finding opportunities for a rewrite can be expensive, as optimizer traverses the AST naively. In this paper, we leverage the need to repeatedly find rewrites, and explore options for making the search faster through indexing and incremental view maintenance (IVM). Concretely, we consider bolt-on approaches that make use of embedded IVM systems like DBToaster, as well as two new approaches: Label-indexing and TreeToaster, an AST-specialized form of IVM. We integrate these approaches into an existing just-in-time data structure compiler and show experimentally that TreeToaster can significantly improve performance with minimal memory overheads.


title: Evolution of Symbiosis in the Game of Life: Three Characteristics of Successful Symbiotes

abstract: In past work, we developed a computational model of the evolution of symbiotic entities (Model-S), based on Conway's Game of Life. In this article, we examine three trends that biologists have observed in the evolution of symbiotes. (1) Management: If one partner is able to control the symbiotic relation, this control can reduce conflict; thus, evolutionary selection favours symbiotes that have a manager. (2) Mutualism: Although partners in a symbiote often have conflicting needs, evolutionary selection favours symbiotes in which partners are better off together inside the symbiote than they would be as individuals outside of the symbiote. (3) Interaction: Repeated interaction among partners in symbiosis tends to promote increasing fitness due to evolutionary selection. We have added new components to Model-S that allow us to observe these three trends in runs of Model-S. The new components are analogous to the practice of staining cells in biology research, to reveal patterns that are not usually visible. When we measure the fitness of a symbiote by making it compete with other symbiotes, we find that fitter symbiotes have significantly more management, mutualism, and interaction than less fit symbiotes. These results confirm the trends observed in nature by biologists. Model-S allows biologists to study these evolutionary trends and other characteristics of symbiosis in ways that are not tractable with living organisms.


title: Curvature-driven feedback on aggregation-diffusion of proteins in lipid bilayers

abstract: Membrane bending is an extensively studied problem from both modeling and experimental perspectives because of the wide implications of curvature generation in cell biology. Many of the curvature generating aspects in membranes can be attributed to interactions between proteins and membranes. These interactions include protein diffusion and formation of aggregates due to protein-protein interactions in the plane of the membrane. Recently, we developed a model that couples the in-plane flow of lipids and diffusion of proteins with the out-of-plane bending of the membrane. Building on this work, here, we focus on the role of explicit aggregation of proteins on the surface of the membrane in the presence of membrane bending and diffusion. We develop a comprehensive framework that includes lipid flow, membrane bending, the entropy of protein distribution, along with an explicit aggregation potential and derive the governing equations for the coupled system. We compare this framework to the Cahn-Hillard formalism to predict the regimes in which the proteins form patterns on the membrane. We demonstrate the utility of this model using numerical simulations to predict how aggregation and diffusion, when coupled with curvature generation, can alter the landscape of membrane-protein interactions.


title: On random compact sets, equidecomposition, and domains of expansion in R^3

abstract: We study random compact subsets of R^3 which can be described as "random Menger sponges". We use those random sets to construct a pair of compact sets A and B in R^3 which are of the same positive measure, such that A can be covered by finitely many translates of B, B can be covered by finitely many translates of A, and yet A and B are not equidecomposable. Furthermore, we construct the first example of a compact subset of R^3 of positive measure which is not a domain of expansion. This answers a question of Adrian Ioana.