title: XNAT-PIC: Extending XNAT to Preclinical Imaging Centers

abstract: Molecular imaging generates large volumes of heterogeneous biomedical imagery with an impelling need of guidelines for handling image data. Although several successful solutions have been implemented for human epidemiologic studies, few and limited approaches have been proposed for animal population studies. Preclinical imaging research deals with a variety of machinery yielding tons of raw data but the current practices to store and distribute image data are inadequate. Therefore, standard tools for the analysis of large image datasets need to be established. In this paper, we present an extension of XNAT for Preclinical Imaging Centers (XNAT-PIC). XNAT is a worldwide used, open-source platform for securely hosting, sharing, and processing of clinical imaging studies. Despite its success, neither tools for importing large, multimodal preclinical image datasets nor pipelines for processing whole imaging studies are yet available in XNAT. In order to overcome these limitations, we have developed several tools to expand the XNAT core functionalities for supporting preclinical imaging facilities. Our aim is to streamline the management and exchange of image data within the preclinical imaging community, thereby enhancing the reproducibility of the results of image processing and promoting open science practices.


title: Origin of the Moon

abstract: The Earth-Moon system is unusual in several respects. The Moon is roughly 1/4 the radius of the Earth - a larger satellite-to-planet size ratio than all known satellites other than Pluto's Charon. The Moon has a tiny core, perhaps with only ~1% of its mass, in contrast to Earth whose core contains nearly 30% of its mass. The Earth-Moon system has a high total angular momentum, implying a rapidly spinning Earth when the Moon formed. In addition, the early Moon was hot and at least partially molten with a deep magma ocean. Identification of a model for lunar origin that can satisfactorily explain all of these features has been the focus of decades of research.


title: On Boundary Conditions in the sub-mesh interaction of the Particle-Particle-Particle-Mesh Algorithm

abstract: The Particle-Particle-Particle-Mesh algorithm elegantly extends the standard Particle-In-Cell scheme by direct summation of interaction that happens over distances below or around mesh size. Generally, this allows for a more accurate description of Coulomb interactions and improves precision in the prediction of key observables. Nevertheless, most implementations neglect electrostatic boundary conditions for the short-ranged interaction that are directly summed. In this paper a variational description of the Particle-Particle-Particle-Mesh algorithm will be developed for the first time and subsequently used to derive temporally and spatially discrete equations of motion. We show that the error committed by neglecting boundary conditions on the short scale is directly tied to the discretization error induced by the computational grid.


title: Towards Optimal Heterogeneity in Lattice Structures

abstract: We present a multi-phase design parameterization to obtain optimized heterogeneous lattice structures. The 3D domain is discretized into a cubical grid wherein each cube has eight distinct unit cell types or phases. When all phases are present, the domain resembles a densely connected ground structure. The cross-section area of beam segments in lattice units, modelled using Timoshenko beam theory, are the design variables. All beam segments in a particular lattice phase have the same area of cross-section to keep the number of design variables low. The optimization problem is formulated for stiff structures and is solved using the optimality criteria algorithm. We present a case study to show the superiority of topology-optimized heterogeneous structures over uniform lattices of a single phase. In order to interpret the phase composition, we perform four basic load tests on single phases, namely, tension or compression, shear, bending, and torsion. The phases are ranked based on the stiffness corresponding to individual loading conditions. The results show that in the optimized structure, the local internal load configuration drives the selection of phases. We also note that micropolar elasticity captures the bulk behaviour of heterogeneous lattice structures, and helps not only to interpret the optimality of phases but also to improve the computational efficiency of the proposed optimization technique.


title: A Flexible Rolling Regression Framework for Time-Varying SIRD models: Application to COVID-19

abstract: The present paper introduces a data-driven framework for describing the time-varying nature of an SIRD model in the context of COVID-19. By embedding a rolling regression in a mixed integer bilevel nonlinear programming problem, our aim is to provide the research community with a model that reproduces accurately the observed changes in the number of infected, recovered, and death cases, while providing information about the time dependency of the parameters that govern the SIRD model. We propose this optimization model and a genetic algorithm to tackle its solution. Moreover, we test this algorithm with 2020 COVID-19 data from the state of Minnesota and found that our results are consistent both qualitatively and quantitatively, thus proving that the framework proposed is an effective an flexible tool to describe the dynamics of a pandemic.


title: A VLA Census of the Galactic H II Region Population

abstract: The Milky Way contains thousands of H II region candidates identified by their characteristic mid-infrared morphology, but lacking detections of ionized gas tracers such as radio continuum or radio recombination line emission. These targets thus remain unconfirmed as H II regions. With only $\sim$2500 confirmed H II regions in the Milky Way, Galactic surveys are deficient by several thousand nebulae when compared to external galaxies with similar star formation rates. Using sensitive 9 GHz radio continuum observations with the Karl G. Jansky Very Large Array (VLA), we explore a sample of H II region candidates in order to set observational limits on the actual total population of Galactic H II regions. We target all infrared-identified "radio quiet" sources from the WISE Catalog of Galactic H II regions between $245^{\circ}\geq\ell\geq90^{\circ}$ with infrared diameters less than 80$^{\prime\prime}$. We detect radio continuum emission from 50% of the targeted H II region candidates, providing strong evidence that most of the radio quiet candidates are bona fide HII regions. We measure the peak and integrated radio flux densities and compare the inferred Lyman continuum fluxes using models of OB-stars. We conclude that stars of approximately spectral type B2 and earlier are able to create H II regions with similar infrared and radio continuum morphologies as the more luminous H II regions created by O-stars. From our 50% detection rate of "radio quiet" sources, we set a lower limit of $\sim$7000 for the H II region population of the Galaxy. Thus the vast majority of the Milky Way's H II regions remain to be discovered.


title: Obstacle Avoidance onboard MAVs using a FMCW RADAR

abstract: Micro Air Vehicles (MAVs) are increasingly being used for complex or hazardous tasks in enclosed and cluttered environments such as surveillance or search and rescue. With this comes the necessity for sensors that can operate in poor visibility conditions to facilitate with navigation and avoidance of objects or people. Radar sensors in particular can provide more robust sensing of the environment when traditional sensors such as cameras fail in the presence of dust, fog or smoke. While extensively used in autonomous driving, miniature FMCW radars on MAVs have been relatively unexplored. This study aims to investigate to what extent this sensor is of use in these environments by employing traditional signal processing such as multi-target tracking and velocity obstacles. The viability of the solution is evaluated with an implementation on board a MAV by running trial tests in an indoor environment containing obstacles and by comparison with a human pilot, demonstrating the potential for the sensor to provide a more robust sense and avoid function in fully autonomous MAVs.


title: Cross-Gradient Aggregation for Decentralized Learning from Non-IID data

abstract: Decentralized learning enables a group of collaborative agents to learn models using a distributed dataset without the need for a central parameter server. Recently, decentralized learning algorithms have demonstrated state-of-the-art results on benchmark data sets, comparable with centralized algorithms. However, the key assumption to achieve competitive performance is that the data is independently and identically distributed (IID) among the agents which, in real-life applications, is often not applicable. Inspired by ideas from continual learning, we propose Cross-Gradient Aggregation (CGA), a novel decentralized learning algorithm where (i) each agent aggregates cross-gradient information, i.e., derivatives of its model with respect to its neighbors' datasets, and (ii) updates its model using a projected gradient based on quadratic programming (QP). We theoretically analyze the convergence characteristics of CGA and demonstrate its efficiency on non-IID data distributions sampled from the MNIST and CIFAR-10 datasets. Our empirical comparisons show superior learning performance of CGA over existing state-of-the-art decentralized learning algorithms, as well as maintaining the improved performance under information compression to reduce peer-to-peer communication overhead. The code is available here on GitHub.


title: Convergence and Inequality in Research Globalization

abstract: The catch-up effect and the Matthew effect offer opposing characterizations of globalization: the former predicts an eventual convergence as the poor can grow faster than the rich due to free exchanges of complementary resources, while the latter, a deepening inequality between the rich and the poor. To understand these effects on the globalization of research, we conduct an in-depth study based on scholarly and patent publications covering STEM research from 218 countries/regions over the past four decades, covering more than 55 million scholarly articles and 1.7 billion citations. Unique to this investigation is the simultaneous examination of both the research output and its impact in the same data set, using a novel machine learning based measure, called saliency, to mitigate the intrinsic biases in quantifying the research impact. The results show that the two effects are in fact co-occurring: there are clear indications of convergence among the high income and upper middle income countries across the STEM fields, but a widening gap is developing that segregates the lower middle and low income regions from the higher income regions. Furthermore, the rate of convergence varies notably among the STEM sub-fields, with the highly strategic area of Artificial Intelligence (AI) sandwiched between fields such as Medicine and Materials Science that occupy the opposite ends of the spectrum. The data support the argument that a leading explanation of the Matthew effect, namely, the preferential attachment theory, can actually foster the catch-up effect when organizations from lower income countries forge substantial research collaborations with those already dominant. The data resoundingly show such collaborations benefit all parties involved, and a case of role reversal can be seen in the Materials Science field where the most advanced signs of convergence are observed.


title: A note on the generalized-hypergeometric solutions of general and single-confluent Heun equations

abstract: We review the series solutions of the general and single-confluent Heun equations in terms of powers, ordinary-hypergeometric and confluent-hypergeometric functions. The conditions under which the expansions reduce to finite sums as well as the cases when the coefficients of power-series expansions obey two-term recurrence relations are discussed. Infinitely many cases for which a solution of a general or single-confluent Heun equation is written through a single generalized hypergeometric function are indicated. It is shown that the parameters of this generalized hypergeometric function are the roots of a certain polynomial equation, which we explicitly present for any given order.


title: Nucleon axial-vector and pseudoscalar form factors, and PCAC relations

abstract: We use a continuum quark+diquark approach to the nucleon bound-state problem in relativistic quantum field theory to deliver parameter-free predictions for the nucleon axial and induced pseudoscalar form factors, $G_A$ and $G_P$, and unify them with the pseudoscalar form factor $G_5$ or, equivalently, the pion-nucleon form factor $G_{\pi NN}$. We explain how partial conservation of the axial-vector current and the associated Goldberger-Treiman relation are satisfied once all necessary couplings of the external current to the building blocks of the nucleon are constructed consistently; in particular, we fully resolve the seagull couplings to the diquark-quark vertices associated with the axial-vector and pseudoscalar currents. Among the results we describe, the following are worth highlighting. A dipole form factor defined by an axial charge $g_A=G_A(0)=1.25(3)$ and a mass-scale $M_A = 1.23(3) m_N$, where $m_N$ is the nucleon mass, can accurately describe the pointwise behaviour of $G_A$. Concerning $G_P$, we obtain the pseudoscalar charge $g_p^\ast = 8.80(23)$, and find that the pion pole dominance approach delivers a reliable estimate of the directly computed result. Our computed value of the pion-nucleon coupling constant, $g_{\pi NN}/m_N =14.02(33)/{\rm GeV}$ is consistent with recent precision determinations.


title: The Topological Origin of the Peierls-Nabarro Barrier

abstract: Crystals and other condensed matter systems described by density waves often exhibit dislocations. Here we show, by considering the topology of the ground state manifolds (GSMs) of such systems, that dislocations in the density phase field always split into disclinations, and that the disclinations themselves are constrained to sit at particular points in the GSM. Consequently, the topology of the GSM forbids zero-energy dislocation glide, giving rise to a Peirels-Nabarro barrier.

categories: cond-mat.soft cond-mat.mtrl-sci math-ph math.MP


title: A Game-Theoretic Analysis of Cross-Ledger Swaps with Packetized Payments

abstract: We propose a game-theoretic framework to study the outcomes of packetized payments, a cross-ledger transaction protocol, with strategic and possibly malicious agents. We derive the transaction failure rate and demonstrate that without disciplinary mechanisms, packetized payments are likely to be incomplete. Our analysis suggests that collateral deposits can prevent malicious agents from taking advantage of the protocol. We further infer that the deposit amount should depend on the underlying asset price volatility or that it should be dynamically adjusted as the price changes.


title: Evaporation of four-dimensional dynamical black holes sourced by the quantum trace anomaly

abstract: We study the evaporation of a four-dimensional spherically symmetric black hole formed in a gravitational collapse. We analyze the back-reaction of a massless quantum scalar field conformally coupled to the scalar curvature by means of the semiclassical Einstein equations. We show that the evaporation is linked to an ingoing negative energy flux at the dynamical horizon and that this flux is induced by the quantum matter trace anomaly outside the black hole horizon whenever a suitable averaged energy condition is satisfied. For illustrative purposes, we evaluate the negative ingoing flux and the corresponding rate of evaporation in the case of a null radiating star described by the Vaidya spacetime.


title: The Heavy Quark Expansion for the Charm Quark

abstract: We show that one can re-arrange the Heavy Quark Expansion for inclusive weak decays of charmed hadrons in such a way that the resulting expansion is an expansion in $\Lambda_{\rm QCD} / m_c$ and $\alpha_s (m_c)$ with order-one coefficients. Unlike in the case of the bottom quark, the leading term includes not only the contribution of the free-quark decay, but also a tower of terms related to matrix elements of four quark operators.


title: Observer Path Planning for Maximum Information

abstract: This paper is concerned with finding an optimal path for an observer, or sensor, moving at a constant speed, which is to estimate the position of a stationary target, using only bearing angle measurements. The generated path is optimal in the sense that, along the path, information, and thus the efficiency of a potential estimator employed, is maximized. In other words, an observer path is deemed optimal if it maximizes information so that the location of the target is estimated with smallest uncertainty, in some sense. We formulate this problem as an optimal control problem maximizing the determinant of the Fisher information matrix, which is one of the possible measures of information. We derive analytical results for optimality using the Maximum Principle. We carry out numerical experiments and discuss the multiple (locally) optimal solutions obtained. We check graphically that the necessary conditions of optimality are verified by the numerical solutions. Finally we provide a comprehensive list of possible extensions for future work.


title: Capelin: Data-Driven Capacity Procurement for Cloud Datacenters using Portfolios of Scenarios -- Extended Technical Report

abstract: Cloud datacenters provide a backbone to our digital society. Inaccurate capacity procurement for cloud datacenters can lead to significant performance degradation, denser targets for failure, and unsustainable energy consumption. Although this activity is core to improving cloud infrastructure, relatively few comprehensive approaches and support tools exist for mid-tier operators, leaving many planners with merely rule-of-thumb judgement. We derive requirements from a unique survey of experts in charge of diverse datacenters in several countries. We propose Capelin, a data-driven, scenario-based capacity planning system for mid-tier cloud datacenters. Capelin introduces the notion of portfolios of scenarios, which it leverages in its probing for alternative capacity-plans. At the core of the system, a trace-based, discrete-event simulator enables the exploration of different possible topologies, with support for scaling the volume, variety, and velocity of resources, and for horizontal (scale-out) and vertical (scale-up) scaling. Capelin compares alternative topologies and for each gives detailed quantitative operational information, which could facilitate human decisions of capacity planning. We implement and open-source Capelin, and show through comprehensive trace-based experiments it can aid practitioners. The results give evidence that reasonable choices can be worse by a factor of 1.5-2.0 than the best, in terms of performance degradation or energy consumption.


title: Spam Prevention Using zk-SNARKs for Anonymous Peer-to-Peer Content Sharing Systems

abstract: Decentralized unpermissioned peer-to-peer networks are inherently vulnerable to spam when they allow arbitrary participants to submit content to a common public index or registry; preventing this is difficult due to the absence of a central arbitrator who can act as a gate-keeper. For this reason indexing of new content even in otherwise decentralized networks (e.g. Bittorrent with DHT, IPFS) has generally been left to centralized services such as torrent sites. Decentralized methods for spam prevention, such as Web of Trust, already exist[1][2] but they require submitters to assume pseudonymous identities and establish trust over time. In this paper we present a method of spam prevention that works under the assumption that the participants are fully anonymous and do not want different submissions of theirs to be linked to each other. By spam we do not specifically mean unsolicited advertising; rather it is the practice of adding a large amount of content in a short time, saturating the capacity of the network, and causing denial of service to peers. The purpose of our solution is to prevent users from saturating the system, and can be described as rate-limiting. The system should be censorship resistant: it should not be possible for submissions to be excluded because of the content itself, or for users to be excluded based on what they submit. We first discuss a solution based on a single, centralized rate-limiter that is censorship resistant and anonymous, then we extend this to a fully decentralized blockchain-based system and present methods to make it economical and scalable.


title: Variance Reduced Training with Stratified Sampling for Forecasting Models

abstract: In large-scale time series forecasting, one often encounters the situation where the temporal patterns of time series, while drifting over time, differ from one another in the same dataset. In this paper, we provably show under such heterogeneity, training a forecasting model with commonly used stochastic optimizers (e.g. SGD) potentially suffers large variance on gradient estimation, and thus incurs long-time training. We show that this issue can be efficiently alleviated via stratification, which allows the optimizer to sample from pre-grouped time series strata. For better trading-off gradient variance and computation complexity, we further propose SCott (Stochastic Stratified Control Variate Gradient Descent), a variance reduced SGD-style optimizer that utilizes stratified sampling via control variate. In theory, we provide the convergence guarantee of SCott on smooth non-convex objectives. Empirically, we evaluate SCott and other baseline optimizers on both synthetic and real-world time series forecasting problems, and demonstrate SCott converges faster with respect to both iterations and wall clock time.


title: A 3D Printing Hexacopter: Design and Demonstration

abstract: 3D printing using robots has garnered significant interest in manufacturing and construction in recent years. A robot's versatility paired with the design freedom of 3D printing offers promising opportunities for how parts and structures are built in the future. However, 3D printed objects are still limited in size and location due to a lack of vertical mobility of ground robots. These limitations severely restrict the potential of the 3D printing process. To overcome these limitations, we develop a hexacopter testbed that can print via fused deposition modeling during flight. We discuss the design of this testbed and develop a simple control strategy for initial print tests. By successfully performing these initial print tests, we demonstrate the feasibility of this approach and lay the groundwork for printing 3D parts and structures with drones.


title: Poisson Yang-Baxter equations and $\mathcal{O}$-operators of Poisson superalgebras

abstract: We investigate connections between $\mathcal {O}$-operators of Poisson superalgebras and skew-symmetric solutions of the Poisson Yang-Baxter equation (PYBE). We prove that a skew-symmetric solution of the PYBE on a Poisson superalgebra can be interpreted as an $\mathcal {O}$-operator associated to the co-regular representation. We show that this connection can be enhanced with symplectic forms when considering non-degenerate skew-symmetric solutions. We also show that $\mathcal {O}$-operators associated to a general representation could give skew-symmetric solutions of the PYBE in certain semi-direct product of Poisson superalgebras.


title: Confined Vortex Surface and Irreversibility. 1. Properties of Exact solution

abstract: We revise the steady vortex surface theory following the recent finding of asymmetric vortex sheets (AM,2021). These surfaces avoid the Kelvin-Helmholtz instability by adjusting their discontinuity and shape. The vorticity collapses to the sheet only in an exceptional case considered long ago by Burgers and Townsend, where it decays as a Gaussian on both sides of the sheet. In generic asymmetric vortex sheets (Shariff,2021), vorticity leaks to one side or another, making such sheets inadequate for vortex sheet statistics and anomalous dissipation. We conjecture that the vorticity in a turbulent flow collapses on a special kind of surface (confined vortex surface, or CVS), satisfying some equations involving the tangent components of the local strain tensor. The most important qualitative observation is that the inequality needed for this solution's stability breaks the Euler dynamics' time reversibility. We interpret this as dynamic irreversibility. We have also represented the enstrophy as a surface integral, conserved in the Navier-Stokes equation in the turbulent limit, with vortex stretching and viscous diffusion terms exactly canceling each other on the CVS surfaces. We have studied the CVS equations for the cylindrical vortex surface for an arbitrary constant background strain with two different eigenvalues. This equation reduces to a particular version of the stationary Birkhoff-Rott equation for the 2D flow with an extra nonanalytic term. We study some general properties of this equation and reduce its solution to a fixed point of a map on a sphere, guaranteed to exist by the Brouwer theorem.


title: Geospatial Perspective Reprojections for Ground-Based Sky Imaging System

abstract: Sky imaging systems use lenses to acquire images concentrating light beams in a sensor. The light beams received by the sky imager have an elevation angle with respect to the device normal. Thus, the pixels in the image contain information from different areas of the sky within the imaging system field of view. The area of the field of view contained in the pixels increases as the elevation angle of the incident light beams decreases. When the sky imager is mounted on a solar tracker, the light beam's angle of incidence in a pixel varies over time. This investigation formulates and compares two geospatial reprojections that transform the original euclidean frame of the imager plane to the geospatial atmosphere cross-section where the sky imager field of view intersects the cloud layer. One assumes that an object (i.e., cloud) moving in the troposphere is sufficiently far so the Earth's surface is approximated \emph{flat}. The other transformation takes into account the curvature of the Earth in the portion of the atmosphere (i.e., voxel) that is recorded. The results show that the differences between the dimensions calculated by both geospatial transformations are in the order of magnitude of kilometers when the Sun's elevation angle is below $30^\circ$.


title: Eigenvalues of singular measures and Connes noncommutative integration

abstract: For a singular measure $\mu$, Ahlfors regular of order $\alpha>0,$ with compact support in $\mathbb{R}^{\mathbf{N}}$ and a pseudodifferential operator $\mathbf{A}$ of order $-l=-\mathbf{N}/2$ we consider the compact operator $\mathbf{T}(P,\mathbf{A}) = \mathbf{A}^*P\mathbf{A}.$ Here $P$ is the signed measure, $P=V\mu$ with density $V$ belonging to the Orlicz class $L^{\Psi,\mu}$ with $\Psi(t)=(t+1)\log(t+1)-t.$ Using eigenvalue estimates for such operators, obtained in \texttt{arXiv:2011.14877}, we establish eigenvalue asymptotics of $\mathbf{T}(P,\mathbf{A})$ for a class of measures, including the ones supported on uniformly rectifiable sets. These results lead to the measurability in the sense of A.Connes of operators $\mathbf{T}(P,\mathbf{A})$ and a formula for the singular trace of these operators, producing a noncommutative version of integral with respect to singular measure.


title: OPTIMADE, an API for exchanging materials data

abstract: The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API through worked examples on each of the public materials databases that support the full API specification.


title: Non-perturbative methods for NN singular interactions

abstract: Chiral Effective Field Theory ($\chi$EFT) has been extensively used to study the $NN$ interaction during the last three decades. In Effective Field Theories (EFTs) the renormalization is performed order by order including the necessary counter terms. Due to the strong character of the $NN$ interaction a non-perturbative resummation is needed. In this work we will review some of the methods proposed to completely remove cutoff dependencies. The methods covered are renormalization with boundary conditions, renormalization with one counter term in momentum space (or equivalently substractive renormalization) and the exact $N/D$ method. The equivalence between the methods up to one renormalization condition will be checked showing results in the $NN$ system. The exact $N/D$ method allows to go beyond the others, and using a toy model it is shown how it can renormalize singular repulsive interactions.


title: Wold Decomposition on Odometer Semigroups

abstract: We establish a Wold-type decomposition for isometric and isometric Nica-covariant representations of the odometer semigroup. These generalize the Wold-type decomposition for commuting pairs of isometries due to Popovici and for pairs of doubly commuting isometries due to S{\l}oci\'nski.


title: Sibyl: Understanding and Addressing the Usability Challenges of Machine Learning In High-Stakes Decision Making

abstract: Machine learning (ML) is being applied to a diverse and ever-growing set of domains. In many cases, domain experts - who often have no expertise in ML or data science - are asked to use ML predictions to make high-stakes decisions. Multiple ML usability challenges can appear as result, such as lack of user trust in the model, inability to reconcile human-ML disagreement, and ethical concerns about oversimplification of complex problems to a single algorithm output. In this paper, we investigate the ML usability challenges that present in the domain of child welfare screening through a series of collaborations with child welfare screeners. Following the iterative design process between the ML scientists, visualization researchers, and domain experts (child screeners), we first identified four key ML challenges and honed in on one promising explainable ML technique to address them (local factor contributions). Then we implemented and evaluated our visual analytics tool, Sibyl, to increase the interpretability and interactivity of local factor contributions. The effectiveness of our tool is demonstrated by two formal user studies with 12 non-expert participants and 13 expert participants respectively. Valuable feedback was collected, from which we composed a list of design implications as a useful guideline for researchers who aim to develop an interpretable and interactive visualization tool for ML prediction models deployed for child welfare screeners and other similar domain experts.


title: Naturally Light Dirac Neutrinos from SU(6)

abstract: A known mechanism for obtaining naturally light Dirac neutrinos is implemented in the context of $SU(6) \to SU(5) \times U(1)_N$.


title: Coherent control and distinguishability of quantum channels via PBS-diagrams

abstract: We introduce a graphical language for coherent control of general quantum channels inspired by practical quantum optical setups involving polarising beam splitters (PBS). As standard completely positive trace preserving maps are known not to be appropriate to represent coherently controlled quantum channels, we propose to instead use purified channels, an extension of Stinespring's dilation. We characterise the observational equivalence of purified channels in various coherent-control contexts, paving the way towards a faithful representation of quantum channels under coherent control.


title: Sequential Place Learning: Heuristic-Free High-Performance Long-Term Place Recognition

abstract: Sequential matching using hand-crafted heuristics has been standard practice in route-based place recognition for enhancing pairwise similarity results for nearly a decade. However, precision-recall performance of these algorithms dramatically degrades when searching on short temporal window (TW) lengths, while demanding high compute and storage costs on large robotic datasets for autonomous navigation research. Here, influenced by biological systems that robustly navigate spacetime scales even without vision, we develop a joint visual and positional representation learning technique, via a sequential process, and design a learning-based CNN+LSTM architecture, trainable via backpropagation through time, for viewpoint- and appearance-invariant place recognition. Our approach, Sequential Place Learning (SPL), is based on a CNN function that visually encodes an environment from a single traversal, thus reducing storage capacity, while an LSTM temporally fuses each visual embedding with corresponding positional data -- obtained from any source of motion estimation -- for direct sequential inference. Contrary to classical two-stage pipelines, e.g., match-then-temporally-filter, our network directly eliminates false-positive rates while jointly learning sequence matching from a single monocular image sequence, even using short TWs. Hence, we demonstrate that our model outperforms 15 classical methods while setting new state-of-the-art performance standards on 4 challenging benchmark datasets, where one of them can be considered solved with recall rates of 100% at 100% precision, correctly matching all places under extreme sunlight-darkness changes. In addition, we show that SPL can be up to 70x faster to deploy than classical methods on a 729 km route comprising 35,768 consecutive frames. Extensive experiments demonstrate the... Baseline code available at https://github.com/mchancan/deepseqslam


title: Instability-induced patterning of a jelling jet

abstract: When a thin stream of aqueous sodium alginate is extruded into a reacting calcium chloride bath, it polymerizes into a soft elastic tube that spontaneously forms helical coils due to the ambient fluid drag. We quantify the onset of this drag-induced instability and its nonlinear evolution using experiments, and explain the results using a combination of scaling, theory and simulations. By co-extruding a second (internal) liquid within the aqueous sodium alginate jet and varying the rates of co-extrusion of the two liquids, as well as the diameter of the jet, we show that we can tune the local composition of the composite filament and the nature of the ensuing instabilities to create soft filaments of variable relative buoyancy, shape and mechanical properties. All together, by harnessing the fundamental varicose (jetting) and sinuous (buckling) instabilities associated with the extrusion of a jelling filament, we show that it is possible to print complex three-dimensional filamentous structures in the ambient fluid.


title: A continuous-state cellular automata algorithm for global optimization

abstract: Cellular automata are capable of developing complex behaviors based on simple local interactions between their elements. Some of these characteristics have been used to propose and improve meta-heuristics for global optimization; however, the properties offered by the evolution rules in cellular automata have not yet been used directly in optimization tasks. Inspired by the complexity that various evolution rules of cellular automata can offer, the continuous-state cellular automata algorithm (CCAA) is proposed. In this way, the CCAA takes advantage of different evolution rules to maintain a balance that maximizes the exploration and exploitation properties in each iteration. The efficiency of the CCAA is proven with 33 test problems widely used in the literature, 4 engineering applications that were also used in recent literature, and the design of adaptive infinite-impulse response (IIR) filters, testing 10 full-order IIR reference functions. The numerical results prove its competitiveness in comparison with state-of-the-art algorithms. The source codes of the CCAA are publicly available at https://github.com/juanseck/CCAA.git


title: The surface brightness - colour relations based on eclipsing binary stars and calibrated with Gaia EDR3

abstract: The surface brightness -- colour relation (SBCR) is a basic tool in establishing precise and accurate distances within the Local Group. Detached eclipsing binary stars with accurately determined radii and trigonometric parallaxes allow for a calibration of the SBCRs with unprecedented accuracy. We analysed four nearby eclipsing binary stars containing late F-type main sequence components: AL Ari, AL Dor, FM Leo and BN Scl. We determined very precise spectroscopic orbits and combined them with high precision ground- and space-based photometry. We derived the astrophysical parameters of their components with mean errors of 0.1% for mass and 0.4% for radius. We combined those four systems with another 24 nearby eclipsing binaries with accurately known radii from the literature for which $Gaia$ EDR3 parallaxes are available, in order to derive the SBCRs. The resulting SBCRs cover stellar spectral types from B9 V to G7 V. For calibrations we used Johnson optical $B$ and $V$, $Gaia$ $G_{\rm BP}$ and $G$ and 2MASS $JHK$ bands. The most precise relations are calibrated using the infrared $K$ band and allow to predict angular diameters of A-, F-, and G-type dwarf and subgiant stars with a precision of 1%.


title: Effects of overlapping sources on cosmic shear estimation: Statistical sensitivity and pixel-noise bias

abstract: In Stage-IV imaging surveys, a significant amount of the cosmologically useful information is due to sources whose images overlap with those of other sources on the sky. The cosmic shear signal is primarily encoded in the estimated shapes of observed galaxies and thus directly impacted by overlaps. We introduce a framework based on the Fisher formalism to analyze effects of overlapping sources (blending) on the estimation of cosmic shear. For the Rubin Observatory Legacy Survey of Space and Time (LSST), we present the expected loss in statistical sensitivity for the ten-year survey due to blending. We find that for approximately 62% of galaxies that are likely to be detected in full-depth LSST images, at least 1% of the flux in their pixels is from overlapping sources. We also find that the statistical correlations between measures of overlapping galaxies and, to a much lesser extent the higher shot noise level due to their presence, decrease the effective number density of galaxies, $N_{eff}$, by $\sim$18%. We calculate an upper limit on $N_{eff}$ of 39.4 galaxies per arcmin$^2$ in $r$ band. We study the impact of varying stellar density on $N_{eff}$ and illustrate the diminishing returns of extending the survey into lower Galactic latitudes. We extend the Fisher formalism to predict the increase in pixel-noise bias due to blending for maximum-likelihood (ML) shape estimators. We find that noise bias is sensitive to the particular shape estimator and measure of ensemble-average shape that is used, and properties of the galaxy that include redshift-dependent quantities such as size and luminosity.


title: DP-InstaHide: Provably Defusing Poisoning and Backdoor Attacks with Differentially Private Data Augmentations

abstract: Data poisoning and backdoor attacks manipulate training data to induce security breaches in a victim model. These attacks can be provably deflected using differentially private (DP) training methods, although this comes with a sharp decrease in model performance. The InstaHide method has recently been proposed as an alternative to DP training that leverages supposed privacy properties of the mixup augmentation, although without rigorous guarantees. In this work, we show that strong data augmentations, such as mixup and random additive noise, nullify poison attacks while enduring only a small accuracy trade-off. To explain these finding, we propose a training method, DP-InstaHide, which combines the mixup regularizer with additive noise. A rigorous analysis of DP-InstaHide shows that mixup does indeed have privacy advantages, and that training with k-way mixup provably yields at least k times stronger DP guarantees than a naive DP mechanism. Because mixup (as opposed to noise) is beneficial to model performance, DP-InstaHide provides a mechanism for achieving stronger empirical performance against poisoning attacks than other known DP methods.


title: Hyperconnected Megacity Parcel Logistic: Joint Parcel Routing and Containerized Consolidation

abstract: In high-speed hyperconnected parcel logistics, intermediate hubs play a critical role in reaching economies of scale through consolidating disperse flows of goods. However, resorting small-size parcels at every hub is resource-intensive and can increase hubs' workload and parcels' total travel time. Such resorting can be reduced by smart containerized consolidation, encapsulating together parcels sharing service level and a subsequent destination. In this study, we introduce an optimization model enabling to assess the impact of such consolidation in reducing the expected total pickup-to-delivery times over a parcel logistic network and its consequences on urban logistic operations sustainability. We provide empirical results for a synthetic urban environment, contrasting consolidation performance over different operational and tactical capabilities, network configurations and demand patterns.


title: Understanding Species Abundance Distributions in Complex Ecosystems of Interacting Species

abstract: Niche and neutral theory are two prevailing, yet much debated, ideas in ecology proposed to explain the patterns of biodiversity. Whereas niche theory emphasizes selective differences between species and interspecific interactions in shaping the community, neutral theory supposes functional equivalence between species and points to stochasticity as the primary driver of ecological dynamics. In this work, we draw a bridge between these two opposing theories. Starting from a Lotka-Volterra (LV) model with demographic noise and random symmetric interactions, we analytically derive the stationary population statistics and species abundance distribution (SAD). Using these results, we demonstrate that the model can exhibit three classes of SADs commonly found in niche and neutral theories and found conditions that allow an ecosystem to transition between these various regimes. Thus, we reconcile how neutral-like statistics may arise from a diverse community with niche differentiation.


title: Computing Sum of Sources over a Classical-Quantum MAC

abstract: We consider the problem of communicating a general bivariate function of two classical sources observed at the encoders of a classical-quantum multiple access channel. Building on the techniques developed for the case of a classical channel, we propose and analyze a coding scheme based on coset codes. The proposed technique enables the decoder recover the desired function without recovering the sources themselves. We derive a new set of sufficient conditions that are weaker than the current known for identified examples. This work is based on a new ensemble of coset codes that are proven to achieve the capacity of a classical-quantum point-to-point channel.


title: Uncertainty guided semi-supervised segmentation of retinal layers in OCT images

abstract: Deep convolutional neural networks have shown outstanding performance in medical image segmentation tasks. The usual problem when training supervised deep learning methods is the lack of labeled data which is time-consuming and costly to obtain. In this paper, we propose a novel uncertainty-guided semi-supervised learning based on a student-teacher approach for training the segmentation network using limited labeled samples and a large number of unlabeled images. First, a teacher segmentation model is trained from the labeled samples using Bayesian deep learning. The trained model is used to generate soft segmentation labels and uncertainty maps for the unlabeled set. The student model is then updated using the softly segmented samples and the corresponding pixel-wise confidence of the segmentation quality estimated from the uncertainty of the teacher model using a newly designed loss function. Experimental results on a retinal layer segmentation task show that the proposed method improves the segmentation performance in comparison to the fully supervised approach and is on par with the expert annotator. The proposed semi-supervised segmentation framework is a key contribution and applicable for biomedical image segmentation across various imaging modalities where access to annotated medical images is challenging


title: Minimax Model Learning

abstract: We present a novel off-policy loss function for learning a transition model in model-based reinforcement learning. Notably, our loss is derived from the off-policy policy evaluation objective with an emphasis on correcting distribution shift. Compared to previous model-based techniques, our approach allows for greater robustness under model misspecification or distribution shift induced by learning/evaluating policies that are distinct from the data-generating policy. We provide a theoretical analysis and show empirical improvements over existing model-based off-policy evaluation methods. We provide further analysis showing our loss can be used for off-policy optimization (OPO) and demonstrate its integration with more recent improvements in OPO.


title: Information-Centric Networking in Wireless Environments: Security Risks and Challenges

abstract: Information-Centric Networking (ICN) has emerged as a paradigm to cope with the lack of built-in security primitives and efficient mechanisms for content distribution of today's Internet. However, deploying ICN in a wireless environment poses a different set of challenges compared to a wired environment, especially when it comes to security. In this paper, we present the security issues that may arise and the attacks that may occur from different points of view when ICN is deployed in wireless environments. The discussed attacks may target both applications and the ICN network itself by exploiting elements of the ICN architecture, such as content names and in-network content caches. Furthermore, we discuss potential solutions to the presented issues and countermeasures to the presented attacks. Finally, we identify future research opportunities and directions.


title: A numerical method for computing the Jordan Canonical Form

abstract: The Jordan Canonical Form of a matrix is highly sensitive to perturbations, and its numerical computation remains a formidable challenge. This paper presents a regularization theory that establishes a well-posed least squares problem of finding the nearest staircase decomposition in the matrix bundle of the highest codimension. A two-staged algorithm is developed for computing the numerical Jordan Canonical Form. At the first stage, the method calculates the Jordan structure of the matrix and an initial approximation to the multiple eigenvalues. The staircase decomposition is then constructed by an iterative algorithm at the second stage. As a result, the numerical Jordan Canonical decomposition along with multiple eigenvalues can be computed with high accuracy even if the underlying matrix is perturbed.


title: Deep J-Sense: Accelerated MRI Reconstruction via Unrolled Alternating Optimization

abstract: Accelerated multi-coil magnetic resonance imaging reconstruction has seen a substantial recent improvement combining compressed sensing with deep learning. However, most of these methods rely on estimates of the coil sensitivity profiles, or on calibration data for estimating model parameters. Prior work has shown that these methods degrade in performance when the quality of these estimators are poor or when the scan parameters differ from the training conditions. Here we introduce Deep J-Sense as a deep learning approach that builds on unrolled alternating minimization and increases robustness: our algorithm refines both the magnetization (image) kernel and the coil sensitivity maps. Experimental results on a subset of the knee fastMRI dataset show that this increases reconstruction performance and provides a significant degree of robustness to varying acceleration factors and calibration region sizes.


title: Dissipative preparation of W states in trapped ion systems

abstract: We present protocols for dissipative entanglement of three trapped-ion qubits and discuss a scheme that uses sympathetic cooling as the dissipation mechanism. This scheme relies on tailored destructive interference to generate any one of six entangled W states in a three-ion qubit space. Using a beryllium-magnesium ion crystal as an example system, we theoretically investigate the protocol's performance and the effects of likely error sources, including thermal secular motion of the ion crystal, calibration imperfections, and spontaneous photon scattering. We estimate that a fidelity of $\sim$ 98 % may be achieved in typical trapped ion experiments with $\sim$ 1 ms interaction time. These protocols avoid timescale hierarchies for faster preparation of entangled states.


title: Robust Estimation of Loss Models for Lognormal Insurance Payment Severity Data

abstract: The primary objective of this scholarly work is to develop two estimation procedures - maximum likelihood estimator (MLE) and method of trimmed moments (MTM) - for the mean and variance of lognormal insurance payment severity data sets affected by different loss control mechanism, for example, truncation (due to deductibles), censoring (due to policy limits), and scaling (due to coinsurance proportions), in insurance and financial industries. Maximum likelihood estimating equations for both payment-per-payment and payment-per-loss data sets are derived which can be solved readily by any existing iterative numerical methods. The asymptotic distributions of those estimators are established via Fisher information matrices. Further, with a goal of balancing efficiency and robustness and to remove point masses at certain data points, we develop a dynamic MTM estimation procedures for lognormal claim severity models for the above-mentioned transformed data scenarios. The asymptotic distributional properties and the comparison with the corresponding MLEs of those MTM estimators are established along with extensive simulation studies. Purely for illustrative purpose, numerical examples for 1500 US indemnity losses are provided which illustrate the practical performance of the established results in this paper.


title: Building a testable shear viscosity across the QCD phase diagram

abstract: Current experiments at the Relativistic Heavy Ion Collider (RHIC) are probing finite baryon densities where the shear viscosity to enthalpy ratio $\eta T/w$ of the Quark Gluon Plasma remains unknown. We use the Hadron Resonance Gas (HRG) model with the most up-to-date hadron list to calculate $\eta T/w$ at low temperatures and at finite baryon densities $\rho_B$. We then match $\eta T/w$ to a QCD-based shear viscosity calculation within the deconfined phase to create a table across $\left\{T,\mu_B\right\}$ for different cross-over and critical point scenarios at a specified location. We find that these new $\eta T/w(T,\mu_B)$ values would require initial conditions at significantly larger $\rho_B$, compared to ideal hydrodynamic trajectories, in order to reach the same freeze-out point.


title: Preliminaries on the Accurate Estimation of the Hurst Exponent Using Time Series

abstract: This article explores the required amount of time series points from a high-speed computer network to accurately estimate the Hurst exponent. The methodology consists in designing an experiment using estimators that are applied to time series addresses resulting from the capture of high-speed network traffic, followed by addressing the minimum amount of point required to obtain in accurate estimates of the Hurst exponent. The methodology addresses the exhaustive analysis of the Hurst exponent considering bias behaviour, standard deviation, and Mean Squared Error using fractional Gaussian noise signals with stationary increases. Our results show that the Whittle estimator successfully estimates the Hurst exponent in series with few points. Based on the results obtained, a minimum length for the time series is empirically proposed. Finally, to validate the results, the methodology is applied to real traffic captures in a high-speed computer network.


title: Anticyclotomic $\large\mu$-invariants of residually reducible Galois Representations

abstract: Let $E$ be an elliptic curve over an imaginary quadratic field $K$, and $p$ be an odd prime such that the residual representation $E[p]$ is reducible. The $\mu$-invariant of the fine Selmer group of $E$ over the anticyclotomic $\mathbb{Z}_p$-extension of $K$ is studied. We do not impose the Heegner hypothesis on $E$, thus allowing certain primes of bad reduction to decompose infinitely in the anticyclotomic $\mathbb{Z}_p$-extension. It is shown that the fine $\mu$-invariant vanishes if certain explicit conditions are satisfied. Further, a partial converse is proven.


title: Pseudo-labeling for Scalable 3D Object Detection

abstract: To safely deploy autonomous vehicles, onboard perception systems must work reliably at high accuracy across a diverse set of environments and geographies. One of the most common techniques to improve the efficacy of such systems in new domains involves collecting large labeled datasets, but such datasets can be extremely costly to obtain, especially if each new deployment geography requires additional data with expensive 3D bounding box annotations. We demonstrate that pseudo-labeling for 3D object detection is an effective way to exploit less expensive and more widely available unlabeled data, and can lead to performance gains across various architectures, data augmentation strategies, and sizes of the labeled dataset. Overall, we show that better teacher models lead to better student models, and that we can distill expensive teachers into efficient, simple students. Specifically, we demonstrate that pseudo-label-trained student models can outperform supervised models trained on 3-10 times the amount of labeled examples. Using PointPillars [24], a two-year-old architecture, as our student model, we are able to achieve state of the art accuracy simply by leveraging large quantities of pseudo-labeled data. Lastly, we show that these student models generalize better than supervised models to a new domain in which we only have unlabeled data, making pseudo-label training an effective form of unsupervised domain adaptation.


title: Spatial asymptotics of Green's function and applications

abstract: We study the spatial asymptotics of Green's function for the 1d Schrodinger operator with operator-valued decaying potential. The bounds on the entropy of the spectral measures are obtained. They are used to establish the presence of a.c. spectrum


title: Canonical currents and heights for K3 surfaces

abstract: We construct canonical positive currents and heights on the boundary of the ample cone of a K3 surface. These are equivariant for the automorphism group and fit together into a continuous family, defined over an enlarged boundary of the ample cone. Along the way, we construct preferred representatives for certain height functions and currents on elliptically fibered surfaces.


title: An Alternative Practice of Tropical Convolution to Traditional Convolutional Neural Networks

abstract: Convolutional neural networks (CNNs) have been used in many machine learning fields. In practical applications, the computational cost of convolutional neural networks is often high with the deepening of the network and the growth of data volume, mostly due to a large amount of multiplication operations of floating-point numbers in convolution operations. To reduce the amount of multiplications, we propose a new type of CNNs called Tropical Convolutional Neural Networks (TCNNs) which are built on tropical convolutions in which the multiplications and additions in conventional convolutional layers are replaced by additions and min/max operations respectively. In addition, since tropical convolution operators are essentially nonlinear operators, we expect TCNNs to have higher nonlinear fitting ability than conventional CNNs. In the experiments, we test and analyze several different architectures of TCNNs for image classification tasks in comparison with similar-sized conventional CNNs. The results show that TCNN can achieve higher expressive power than ordinary convolutional layers on the MNIST and CIFAR10 image data set. In different noise environments, there are wins and losses in the robustness of TCNN and ordinary CNNs.


title: Dynamical dark energy models in the light of Gravitational-Wave Transient Catalogues

abstract: The study of current gravitational waves (GW) catalogues provide an interesting model independent way to understand further the nature of dark energy. In this work, we present an update of the constrains related to dynamical dark energy parametrisations using recent Gravitational-Wave Transient catalogues (GWTC-1 and GWTC-2) along with Type Ia supernova (SNeIa) and Cosmic Chronometers (CC) catalogues. According to our Bayesian results using the full SNeIa+CC+GW database, the $\Lambda$CDM model shows a strong preference against two dark energy parameterisation known as Barboza-Alcaniz (BA) and the Low Correlation (LC) models. Also, we obtain a very strong preference against the Chevallier-Polarski-Linder (CPL) model. Furthermore, we generated a mock GW catalogue and estimate that we require approximately 1000 standard sirens to have a constrain of $H_0$ within 1\% relative error, quantity that is out of reach of current standard sirens candidates in GWTC-1 and GWTC-2 catalogues.


title: Structural Twinning-induced Insulating Phase in CrN (111) Films

abstract: Electronic states of a correlated material can be effectively modified by structural variations delivered from a single-crystal substrate. In this letter, we show that the CrN films grown on MgO (001) substrates have a (001) orientation, whereas the CrN films on {\alpha}-Al2O3 (0001) substrates are oriented along (111) direction parallel to the surface normal. Transport properties of CrN films are remarkably different depending on crystallographic orientations. The critical thickness for the metal-insulator transition (MIT) in CrN 111 films is significantly larger than that of CrN 001 films. In contrast to CrN 001 films without apparent defects, scanning transmission electron microscopy results reveal that CrN 111 films exhibit strain-induced structural defects, e. g. the periodic horizontal twinning domains, resulting in an increased electron scattering facilitating an insulating state. Understanding the key parameters that determine the electronic properties of ultrathin conductive layers is highly desirable for future technological applications.


title: Design of an Affordable Prosthetic Arm Equipped with Deep Learning Vision-Based Manipulation

abstract: Many amputees throughout the world are left with limited options to personally own a prosthetic arm due to the expensive cost, mechanical system complexity, and lack of availability. The three main control methods of prosthetic hands are: (1) body-powered control, (2) extrinsic mechanical control, and (3) myoelectric control. These methods can perform well under a controlled situation but will often break down in clinical and everyday use due to poor robustness, weak adaptability, long-term training, and heavy mental burden during use. This paper lays the complete outline of the design process of an affordable and easily accessible novel prosthetic arm that reduces the cost of prosthetics from $10,000 to $700 on average. The 3D printed prosthetic arm is equipped with a depth camera and closed-loop off-policy deep learning algorithm to help form grasps to the object in view. Current work in reinforcement learning masters only individual skills and is heavily focused on parallel jaw grippers for in-hand manipulation. In order to create generalization, which better performs real-world manipulation, the focus is specifically on using the general framework of Markov Decision Process (MDP) through scalable learning with off-policy algorithms such as deep deterministic policy gradient (DDPG) and to study this question in the context of grasping a prosthetic arm. We were able to achieve a 78% grasp success rate on previously unseen objects and generalize across multiple objects for manipulation tasks. This work will make prosthetics cheaper, easier to use and accessible globally for amputees. Future work includes applying similar approaches to other medical assistive devices where a human is interacting with a machine to complete a task.


title: Near-Room Temperature Ferromagnetic Insulating State in Highly Distorted LaCoO2.5 with CoO5 Square Pyramids

abstract: Dedicated control of oxygen vacancies is an important route to functionalizing complex oxide films. It is well-known that tensile strain significantly lowers the oxygen vacancy formation energy, whereas compressive strain plays a minor role. Thus, atomically reconstruction by extracting oxygen from a compressive-strained film is challenging. Here we report an unexpected LaCoO2.5 phase with a zigzag-like oxygen vacancy ordering through annealing a compressive-strained LaCoO3 in vacuum. The synergetic tilt and distortion of CoO5 square pyramids with large La and Co shifts are quantified using scanning transmission electron microscopy. The large in-plane expansion of CoO5 square pyramids weaken the crystal-field splitting and facilitated the ordered high-spin state of Co2+, which produces an insulating ferromagnetic state with a Curie temperature of ~284 K and a saturation magnetization of ~0.25 {\mu}B/Co. These results demonstrate that extracting targeted oxygen from a compressive-strained oxide provides an opportunity for creating unexpected crystal structures and novel functionalities.


title: Dimensional Control of Octahedral Tilt in SrRuO3 via Infinite-layered Oxides

abstract: Manipulation of octahedral distortion at atomic length scale is an effective means to tune the physical ground states of functional oxides. Previous work demonstrates that epitaxial strain and film thickness are variable parameters to modify the octahedral rotation and tilt. However, selective control of bonding geometry by structural propagation from adjacent layers is rarely studied. Here we propose a new route to tune the ferromagnetic response in SrRuO3 (SRO) ultrathin layers by oxygen coordination of adjacent SrCuO2 (SCO) layers. The infinite-layered CuO2 in SCO exhibits a structural transformation from "planar-type" to "chain-type" as reducing film thickness. These two orientations dramatically modify the polyhedral connectivity at the interface, thus altering the octahedral distortion of SRO. The local structural variation changes the spin state of Ru and hybridization strength between Ru 4d and O 2p orbitals, leading to a significant change in the magnetoresistance and anomalous Hall resistivity of SRO layers. These findings could launch further investigations into adaptive control of magnetoelectric properties in quantum oxide heterostructures using oxygen coordination.


title: Experimental Mathematics Approach to Gauss Diagrams Realizability

abstract: A Gauss diagram (or, more generally, a chord diagram) consists of a circle and some chords inside it. Gauss diagrams are a well-established tool in the study of topology of knots and of planar and spherical curves. Not every Gauss diagram corresponds to a knot (or an immersed curve); if it does, it is called realizable. A classical question of computational topology asked by Gauss himself is which chords diagrams are realizable. An answer was first discovered in the 1930s by Dehn, and since then many efficient algorithms for checking realizability of Gauss diagrams have been developed. Recent studies in Grinblat-Lopatkin (2018,2020) and Biryukov (2019) formulated especially simple conditions related to realizability which are expressible in terms of parity of chords intersections. The simple form of these conditions opens an opportunity for experimental investigation of Gauss diagrams using constraint satisfaction and related techniques. In this paper we report on our experiments with Gauss diagrams of small sizes (up to 11 chords) using implementations of these conditions and other algorithms in logic programming language Prolog. In particular, we found a series of counterexamples showing that that realizability criteria established by Grinblat and Lopatkin (2018,2020) and Biryukov (2019) are not completely correct.


title: Accurate calculation of excitonic signatures in the absorption spectrum of BiSBr using semiconductor Bloch equations

abstract: In order to realize the significant potential of optical materials such as metal halides, computational techniques which give accurate optical properties are needed, which can work hand-in-hand with experiments to generate high efficiency devices. In this work a computationally efficient technique based on semiconductor Bloch equations (SBEs) is developed and applied to the material BiSBr. This approach gives excellent agreement with the experimental optical gap, and also agrees closely with the excitonic stabilisation energy and the absorption spectrum computed using the far more computationally demanding \textit{ab initio} Bethe-Salpeter approach. The SBE method is a good candidate for theoretical spectroscopy on large- or low dimensional systems which are too computationally expensive for an \textit{ab initio} treatment.


title: Estimation of attenuation of gravitational waves by Bose-Einstein condensate dark matter halos using Bogoliubov-de Gennes theory

abstract: We consider a gravitational plane wave passing through a galactic dark matter halo composed of weakly self-interacting, self-gravitating, Bose-Einstein condensate of ultralight particles. Treating the gravitational wave as a time dependent perturbation, we study energy transfer between the gravitational wave and the Bose-Einstein condensate by applying linear response theory to a non-uniform condensate described by the Bogoliubov-de Gennes theory, and compute the fractional loss in gravitational wave energy. We apply our results to investigate the extent to which this loss effects the estimation of the distance between the gravitational wave source and the earth. We show that the effect is negligible.


title: Design and Analysis of Delayed Bit-Interleaved Coded Modulation with LDPC Codes

abstract: This paper investigates the design and performance of delayed bit-interleaved coded modulation (DBICM) with low-density parity-check (LDPC) codes. For Gray labeled square $M$-ary quadrature amplitude modulation (QAM) constellations, we investigate the optimal delay scheme with the largest spectrum efficiency of DBICM for a fixed maximum number of delayed time slots and a given signal-to-noise ratio. When analyzing the capacity of DBICM, we find two important properties: the capacity improvement due to delayed coded bits being mapped to the real and imaginary parts of the transmitted symbols are independent of each other; a pair of delay schemes with delayed coded bits having identical bit-channel capacity lead to equivalent DBICM capacity. Using these two properties, we efficiently optimize the delay scheme for any uniform Gray-QAM systems. Furthermore, these two properties enable efficient LDPC code designs regarding unequal error protection via bit-channel type classifications. Moreover, we use protograph-based extrinsic information transfer charts to jointly optimize degree distributions and channel assignments of LDPC codes and propose a constrained progressive edge growth like algorithm to jointly construct LDPC codes and bit-interleavers for DBICM, taking distinctive bit-channel's capacity into account. Simulation results demonstrate that the designed LDPC coded DBICM systems significantly outperform LDPC coded BICM systems.


title: A vesicle microrheometer for high-throughput viscosity measurements of lipid and polymer membranes

abstract: Viscosity is a key property of cell membranes that controls mobility of embedded proteins and membrane remodeling. Measuring it is challenging because existing approaches involve complex experimental designs and/or models, and the applicability of some is limited to specific systems and membrane compositions. As a result there is scarcity of data and the reported values for membrane viscosity vary by orders of magnitude for the same system. Here, we show how viscosity of bilayer membranes can be obtained from the transient deformation of giant unilamellar vesicles. The approach enables a non-invasive, probe-independent and high-throughput measurement of the viscosity of bilayers made of lipids or polymers with a wide range of compositions and phase state. Pure lipid and single-phase mixed bilayers are found to behave as Newtonian fluids with strain-rate independent viscosity, while phase-separated and diblock-copolymers systems exhibit shear-thinning in the explored range of strain rates 1-2000 $s^{-1}$. The results also reveal that electrically polarized bilayers can be significantly more viscous than charge-neutral bilayers. These findings suggest that biomembrane viscosity is a dynamic property that can be actively modulated not only by composition but also by membrane polarization, e.g., as in action potentials.


title: Surrogate-assisted cooperative signal optimization for large-scale traffic networks

abstract: Reasonable setting of traffic signals can be very helpful in alleviating congestion in urban traffic networks. Meta-heuristic optimization algorithms have proved themselves to be able to find high-quality signal timing plans. However, they generally suffer from performance deterioration when solving large-scale traffic signal optimization problems due to the huge search space and limited computational budget. Directing against this issue, this study proposes a surrogate-assisted cooperative signal optimization (SCSO) method. Different from existing methods that directly deal with the entire traffic network, SCSO first decomposes it into a set of tractable sub-networks, and then achieves signal setting by cooperatively optimizing these sub-networks with a surrogate-assisted optimizer. The decomposition operation significantly narrows the search space of the whole traffic network, and the surrogate-assisted optimizer greatly lowers the computational burden by reducing the number of expensive traffic simulations. By taking Newman fast algorithm, radial basis function and a modified estimation of distribution algorithm as decomposer, surrogate model and optimizer, respectively, this study develops a concrete SCSO algorithm. To evaluate its effectiveness and efficiency, a large-scale traffic network involving crossroads and T-junctions is generated based on a real traffic network. Comparison with several existing meta-heuristic algorithms specially designed for traffic signal optimization demonstrates the superiority of SCSO in reducing the average delay time of vehicles.


title: Linear and Nonlinear Viscoelasticity of Concentrated Thermoresponsive Microgel Suspensions

abstract: This is an integrated experimental and theoretical study of the dynamics and rheology of self-crosslinked, slightly charged, temperature responsive soft Poly(N-isopropylacrylamide) (pNIPAM) microgels over a wide range of concentration and temperature spanning the sharp change in particle size and intermolecular interactions across the lower critical solution temperature (LCST). Dramatic, non-monotonic changes in viscoelasticity are observed with temperature, with distinctive concentration dependences in the dense fluid, glassy, and soft-jammed states. Motivated by our experimental observations, we formulate a minimalistic model for the size dependence of a single microgel particle and the change of interparticle interaction from purely repulsive to attractive upon heating. Using microscopic equilibrium and time-dependent statistical mechanical theories, theoretical predictions are quantitatively compared with experimental measurements of the shear modulus. Good agreement is found for the nonmonotonic temperature behavior that originates as a consequence of the competition between reduced microgel packing fraction and increasing interpar-ticle attractions. Testable predictions are made for nonlinear rheological properties such as the yield stress and strain. To the best of our knowledge, this is the first attempt to quantitatively understand in a unified manner the viscoelasticity of dense, temperature-responsive microgel suspensions spanning a wide range of temperatures and concentrations.


title: Quantum circuits with many photons on a programmable nanophotonic chip

abstract: Growing interest in quantum computing for practical applications has led to a surge in the availability of programmable machines for executing quantum algorithms. Present day photonic quantum computers have been limited either to non-deterministic operation, low photon numbers and rates, or fixed random gate sequences. Here we introduce a full-stack hardware-software system for executing many-photon quantum circuits using integrated nanophotonics: a programmable chip, operating at room temperature and interfaced with a fully automated control system. It enables remote users to execute quantum algorithms requiring up to eight modes of strongly squeezed vacuum initialized as two-mode squeezed states in single temporal modes, a fully general and programmable four-mode interferometer, and genuine photon number-resolving readout on all outputs. Multi-photon detection events with photon numbers and rates exceeding any previous quantum optical demonstration on a programmable device are made possible by strong squeezing and high sampling rates. We verify the non-classicality of the device output, and use the platform to carry out proof-of-principle demonstrations of three quantum algorithms: Gaussian boson sampling, molecular vibronic spectra, and graph similarity.


title: Privacy-Preserving Decentralized Multi-Agent Cooperative Optimization -- Paradigm Design and Privacy Analysis

abstract: Large-scale multi-agent cooperative control problems have materially enjoyed the scalability, adaptivity, and flexibility of decentralized optimization. However, due to the mandatory iterative communications between the agents and the system operator, the decentralized architecture is vulnerable to malicious attacks and privacy breach. Current research on addressing privacy preservation of both agents and the system operator in cooperative decentralized optimization with strongly coupled objective functions and constraints is still primitive. To fill in the gaps, this paper proposes a novel privacy-preserving decentralized optimization paradigm based on Paillier cryptosystem. The proposed paradigm achieves ideal correctness and security, as well as resists attacks from a range of adversaries. The efficacy and efficiency of the proposed approach are verified via numerical simulations and a real-world physical platform.


title: Robust Place Recognition using an Imaging Lidar

abstract: We propose a methodology for robust, real-time place recognition using an imaging lidar, which yields image-quality high-resolution 3D point clouds. Utilizing the intensity readings of an imaging lidar, we project the point cloud and obtain an intensity image. ORB feature descriptors are extracted from the image and encoded into a bag-of-words vector. The vector, used to identify the point cloud, is inserted into a database that is maintained by DBoW for fast place recognition queries. The returned candidate is further validated by matching visual feature descriptors. To reject matching outliers, we apply PnP, which minimizes the reprojection error of visual features' positions in Euclidean space with their correspondences in 2D image space, using RANSAC. Combining the advantages from both camera and lidar-based place recognition approaches, our method is truly rotation-invariant, and can tackle reverse revisiting and upside down revisiting. The proposed method is evaluated on datasets gathered from a variety of platforms over different scales and environments. Our implementation and datasets are available at https://git.io/image-lidar


title: Should multilevel methods for discontinuous Galerkin discretizations use discontinuous interpolation operators?

abstract: Multi-level preconditioners for Discontinuous Galerkin (DG) discretizations are widely used to solve elliptic equations, and a main ingredient of such solvers is the interpolation operator to transfer information from the coarse to the fine grid. Classical interpolation operators give continuous interpolated values, but since DG solutions are naturally discontinuous, one might wonder if one should not use discontinuous interpolation operators for DG discretizations. We consider a discontinuous interpolation operator with a parameter that controls the discontinuity, and determine the optimal choice for the discontinuity, leading to the fastest solver for a specific 1D symmetric interior penalty DG discretization model problem. We show in addition that our optimization delivers a perfectly clustered spectrum with a high geometric multiplicity, which is very advantageous for a Krylov solver using the method as its preconditioner. Finally, we show the applicability of the optimal choice to higher dimensions.


title: The Color and Brightness of the F-Corona Inferred from the 2019 July 2 Total Solar Eclipse

abstract: Total solar eclipses (TSEs) provide a unique opportunity to quantify the properties of the K-corona (electrons), F-corona (dust) and E-corona (ions) continuously from the solar surface out to a few solar radii. We apply a novel inversion method to separate emission from the K- and F-corona continua using unpolarized total brightness (tB) observations from five 0.5 nm bandpasses acquired during the 2019 July 2 TSE between 529.5 nm and 788.4 nm. The wavelength dependence relative to the photosphere (i.e., color) of the F-corona itself is used to infer the tB of the K- and F-corona for each line-of-sight. We compare our K-corona emission results with the Mauna Loa Solar Observatory (MLSO) K-Cor polarized brightness (pB) observations from the day of the eclipse, and the forward modeled K-corona intensity from the Predictive Science Inc. (PSI) Magnetohydrodynamic (MHD) model prediction. Our results are generally consistent with previous work and match both the MLSO data and PSI-MHD predictions quite well, supporting the validity of our approach and of the PSI-MHD model. However, we find that the tB of the F-corona is higher than expected in the low corona, perhaps indicating that the F-corona is slightly polarized -- challenging the common assumption that the F-corona is entirely unpolarized.


title: A Computational Design and Evaluation Tool for 3D Structures with Planar Surfaces

abstract: Three dimensional (3D) structures composed of planar surfaces can be build out of accessible materials using easier fabrication technique with shorter fabrication time. To better design 3D structures with planar surfaces, realistic models are required to understand and evaluate mechanical behaviors. Existing design tools are either effort-consuming (e.g. finite element analysis) or bounded by assumptions (e.g. numerical solutions). In this project, We have built a computational design tool that is (1) capable of rapidly and inexpensively evaluating planar surfaces in 3D structures, with sufficient computational efficiency and accuracy; (2) applicable to complex boundary conditions and loading conditions, both isotropic materials and orthotropic materials; and (3) suitable for rapid accommodation when design parameters need to be adjusted. We demonstrate the efficiency and necessity of this design tool by evaluating a glass table as well as a wood bookcase, and iteratively designing an origami gripper to satisfy performance requirements. This design tool gives non-expert users as well as engineers a simple and effective modus operandi in structural design.


title: Rank bias for elliptic curves mod $p$

abstract: We conjecture that, for a fixed prime $p$, rational elliptic curves with higher rank tend to have more points mod $p$. We show that there is an analogous bias for modular forms with respect to root numbers, and conjecture that the order of the rank bias for elliptic curves is greater than that of the root number bias for modular forms.


title: An inexact proximal point method for variational inequality on Hadamard manifolds

abstract: In this paper we present an inexact proximal point method for variational inequality problem on Hadamard manifolds and study its convergence properties. The proposed algorithm is inexact in two sense. First, each proximal subproblem is approximated by using the enlargement of the vector field in consideration and then the next iterated is obtained by solving this subproblem allowing a suitable error tolerance. As an application, we obtain an inexact proximal point method for constrained optimization problems, equilibrium problems and nonlinear optimization problems on Hadamard manifolds.


title: On the Hubble constant tension in the SNe Ia Pantheon sample

abstract: The Hubble constant ($H_0$) tension between Type Ia Supernovae (SNe Ia) and Planck measurements ranges from 4 to 6 $\sigma$. To investigate this tension, we estimate $H_{0}$ in the $\Lambda$CDM and $w_{0}w_{a}$CDM models by dividing the Pantheon sample, the largest compilation of SNe Ia, into 3, 4, 20 and 40 bins. We fit the extracted $H_{0}$ values with a function mimicking the redshift evolution: $g(z)={H_0}(z)=\tilde{H}_0/(1+z)^\alpha$, where $\alpha$ indicates an evolutionary parameter and $\tilde{H}_0=H_0$ at $z=0$. We set the absolute magnitude of SNe Ia so that $H_0=73.5\,\, \textrm{km s}^{-1}\,\textrm{Mpc}^{-1}$, and we fix fiducial values for $\Omega_{0m}^{\Lambda CDM}=0.298$ and $\Omega_{0m}^{w_{0}w_{a}CDM}=0.308$. We find that $H_0$ evolves with redshift, showing a slowly decreasing trend, with $\alpha$ coefficients consistent with zero only from 1.2 to 2.0 $\sigma$. Although the $\alpha$ coefficients are compatible with 0 in 3 $\sigma$, this however may affect cosmological results. We measure locally a variation of $H_0(z=0)-H_0(z=1)=0.4\, \textrm{km s}^{-1}\,\textrm{Mpc}^{-1}$ in 3 and 4 bins. Extrapolating ${H_0}(z)$ to $z=1100$, the redshift of the last scattering surface, we obtain values of $H_0$ compatible in 1 $\sigma$ with Planck measurements independently of cosmological models and number of bins we investigated. Thus, we have reduced the $H_0$ tension from $54\%$ to $72\%$ for the $\Lambda$CDM and $w_{0}w_{a}$CDM models, respectively. If the decreasing trend of $H_0(z)$ is real, it could be due to astrophysical selection effects or to modified gravity.


title: Arc plasma ablation of quartz crystals

abstract: Spherical quartz stones of around 1 cm in diameter have been exposed to anodic arc discharges in a helium atmosphere at 300 Torr. The arc current flowing between the graphite electrodes was set either in continuous DC mode (30-150 A) or in pulsed mode at 2 Hz (220 A peak). The ablation rate in each sample was systematically measured after several seconds of arc plasma treatment. Optical emission spectroscopy (OES) diagnostics and 2-D fluid simulations of the arc discharge have shed light on the heat flux transport and the heating mechanisms of the quartz crystals. A linear correlation is found between the absorbed power density and the resulting rate of penetration, which yields a maximal value of 15 cm/h for approximately 150 W/cm2. The linear fit on the slope provides a specific energy of 40 kJ/cm3. The incident energy flux onto the sample surface promoted a phase transition from crystalline to glassy silica, as characterized via Raman spectroscopy. This study points out the strong potential of arc plasma technology for geothermal drilling applications.


title: Central Moment Lattice Boltzmann Method on a Rectangular Lattice

abstract: Simulating inhomogeneous flows with different characteristic scales in different coordinate directions using the collide-and-stream based lattice Boltzmann methods (LBM) can be accomplished efficiently using rectangular lattice grids. We develop and investigate a new rectangular central moment LBM based on non-orthogonal moment basis (referred to as RC-LBM). The equilibria to which the central moments relax under collision in this approach are obtained from matching with those corresponding to the continuous Maxwell distribution. A Chapman-Enskog analysis is performed to derive the correction terms to the second order moment equilibria involving the grid aspect ratio and velocity gradients that restores the isotropy of the viscous stress tensor and eliminates the non-Galilean invariant cubic velocity terms of the resulting hydrodynamical equations. A special case of this rectangular formulation involving the raw moments (referred to as the RNR-LBM) is also constructed. The resulting schemes represent a considerable simplification, especially for the transformation matrices and isotropy corrections, and improvement over the existing MRT-LB schemes on rectangular lattice grids that use orthogonal moment basis. Numerical validation study of both the RC-LBM and RNR-LBM for a variety of benchmark flow problems are performed that show good accuracy at various grid aspect ratios. The ability of our proposed schemes to simulate flows using relatively lower grid aspect ratios than considered in prior rectangular LB approaches is demonstrated. Furthermore, simulations reveal the superior stability characteristics of the RC-LBM over RNR-LBM in handling shear flows at lower viscosities and/or higher characteristic velocities. In addition, computational advantages of using our rectangular LB formulation in lieu of that based on the square lattice is shown.


title: Interaction between a Waveguide Mode and an Absorptive Medium

abstract: Absorption spectroscopy studies the absorption of electromagnetic wave in a material sample under test. The interacting electromagnetic wave can be propagating in free space or in a waveguide. The waveguide-based absorption spectroscopy method has several advantages compared to the free space setup. The effect of the waveguide cross section on the interaction between the waveguide mode and the sample can be expressed by a factor, called interaction factor. In this article, a new formulation for the interaction factor is derived. It is shown that this factor is inversely proportional to the energy velocity of the waveguide mode.


title: Deblurring Processor for Motion-Blurred Faces Based on Generative Adversarial Networks

abstract: Low-quality face image restoration is a popular research direction in today's computer vision field. It can be used as a pre-work for tasks such as face detection and face recognition. At present, there is a lot of work to solve the problem of low-quality faces under various environmental conditions. This paper mainly focuses on the restoration of motion-blurred faces. In increasingly abundant mobile scenes, the fast recovery of motion-blurred faces can bring highly effective speed improvements in tasks such as face matching. In order to achieve this goal, a deblurring method for motion-blurred facial image signals based on generative adversarial networks(GANs) is proposed. It uses an end-to-end method to train a sharp image generator, i.e., a processor for motion-blurred facial images. This paper introduce the processing progress of motion-blurred images, the development and changes of GANs and some basic concepts. After that, it give the details of network structure and training optimization design of the image processor. Then we conducted a motion blur image generation experiment on some general facial data set, and used the pairs of blurred and sharp face image data to perform the training and testing experiments of the processor GAN, and gave some visual displays. Finally, MTCNN is used to detect the faces of the image generated by the deblurring processor, and compare it with the result of the blurred image. From the results, the processing effect of the deblurring processor on the motion-blurred picture has a significant improvement both in terms of intuition and evaluation indicators of face detection.


title: Actinide crystallization and fission reactions in cooling white dwarf stars

abstract: The first solids that form as a cooling white dwarf (WD) starts to crystallize are expected to be greatly enriched in actinides. This is because the melting points of WD matter scale as $Z^{5/3}$ and actinides have the largest charge $Z$. We estimate that the solids may be so enriched in actinides that they could support a fission chain reaction. This reaction could ignite carbon burning and lead to the explosion of an isolated WD in a thermonuclear supernova (SN Ia). Our mechanism could potentially explain SN Ia with sub-Chandrasekhar ejecta masses and short delay times.


title: A high occurrence of nuclear star clusters in faint Coma galaxies, and the roles of mass and environment

abstract: We use deep high resolution \textit{HST/ACS} imaging of two fields in the core of the Coma cluster to investigate the occurrence of nuclear star clusters (NSCs) in quiescent dwarf galaxies as faint as $M_{I} = -10$ mag. We employ a hierarchical Bayesian logistic regression framework to model the faint end of the nucleation fraction ($f_{n}$) as a function of both galaxy luminosity and environment. We find that $f_n$ is remarkably high in Coma: at $M_{I} \approx -13$ mag half of the cluster dwarfs still host prominent NSCs. Comparison with dwarf systems in nearby clusters and groups shows that, within the uncertainties, the rate at which the probability of nucleation varies with galaxy luminosity is nearly universal. On the other hand, the fraction of nucleated galaxies at fixed luminosity does exhibit an environmental dependence. More massive environments feature higher nucleation fractions and fainter values of the half-nucleation luminosity, which roughly scales with host halo virial mass as $L_{I,f_{n50}} \propto \mathcal{M}_{200}^{-0.2}$. Our results reinforce the role of galaxy luminosity/mass as a major driver of the efficiency of NSC formation and also indicate a clear secondary dependence on the environment, hence paving the way to more refined theoretical models.


title: Periodic Updates for Constrained OCO with Application to Large-Scale Multi-Antenna Systems

abstract: In many dynamic systems, decisions on system operation are updated over time, and the decision maker requires an online learning approach to optimize its strategy in response to the changing environment. When the loss and constraint functions are convex, this belongs to the general family of online convex optimization (OCO). In existing OCO works, the environment is assumed to vary in a time-slotted fashion, while the decisions are updated at each time slot. However, many wireless communication systems permit only periodic decision updates, i.e., each decision is fixed over multiple time slots, while the environment changes between the decision epochs. The standard OCO model is inadequate for these systems. Therefore, in this work, we consider periodic decision updates for OCO. We aim to minimize the accumulation of time-varying convex loss functions, subject to both short-term and long-term constraints. Information about the loss functions within the current update period may be incomplete and is revealed to the decision maker only after the decision is made. We propose an efficient algorithm, termed Periodic Queueing and Gradient Aggregation (PQGA), which employs novel periodic queues together with possibly multi-step aggregated gradient descent to update the decisions over time. We derive upper bounds on the dynamic regret, static regret, and constraint violation of PQGA. As an example application, we study the performance of PQGA in a large-scale multi-antenna system shared by multiple wireless service providers. Simulation results show that PQGA converges fast and substantially outperforms the known best alternative.


title: Minimization of p-Laplacian via the Finite Element Method in MATLAB

abstract: Minimization of energy functionals is based on a discretization by the finite element method and optimization by the trust-region method. A key tool is a local evaluation of the approximated gradients together with sparsity of the resulting Hessian matrix. We describe a vectorized MATLAB implementation of the p-Laplace problem in one and two space-dimensions, however it is easily applicable to other energy formulations.


title: Differentiable Neural Architecture Learning for Efficient Neural Network Design

abstract: Automated neural network design has received ever-increasing attention with the evolution of deep convolutional neural networks (CNNs), especially involving their deployment on embedded and mobile platforms. One of the biggest problems that neural architecture search (NAS) confronts is that a large number of candidate neural architectures are required to train, using, for instance, reinforcement learning and evolutionary optimisation algorithms, at a vast computation cost. Even recent differentiable neural architecture search (DNAS) samples a small number of candidate neural architectures based on the probability distribution of learned architecture parameters to select the final neural architecture. To address this computational complexity issue, we introduce a novel \emph{architecture parameterisation} based on scaled sigmoid function, and propose a general \emph{Differentiable Neural Architecture Learning} (DNAL) method to optimize the neural architecture without the need to evaluate candidate neural networks. Specifically, for stochastic supernets as well as conventional CNNs, we build a new channel-wise module layer with the architecture components controlled by a scaled sigmoid function. We train these neural network models from scratch. The network optimization is decoupled into the weight optimization and the architecture optimization. We address the non-convex optimization problem of neural architecture by the continuous scaled sigmoid method with convergence guarantees. Extensive experiments demonstrate our DNAL method delivers superior performance in terms of neural architecture search cost. The optimal networks learned by DNAL surpass those produced by the state-of-the-art methods on the benchmark CIFAR-10 and ImageNet-1K dataset in accuracy, model size and computational complexity.


title: Exoplanet Statistics and Theoretical Implications

abstract: In the last few years, significant advances have been made in understanding the distributions of exoplanet populations and the architecture of planetary systems. We review the recent progress of planet statistics, with a focus on the inner <~ 1 AU region of the planetary system that has been fairly thoroughly surveyed by the Kepler mission. We also discuss the theoretical implications of these statistical results for planet formation and dynamical evolution.


title: Identifying the Coronal Source Regions of Solar Wind Streams from Total Solar Eclipse Observations and in situ Measurements Extending Over a Solar Cycle

abstract: This Letter capitalizes on a unique set of total solar eclipse observations, acquired between 2006 and 2020, in white light, Fe XI 789.2 nm ($\rm T_{fexi}$ = $1.2 \pm 0.1$ MK) and Fe XIV 530.3 nm ($\rm T_{fexiv}$ = $ 1.8 \pm 0.1$ MK) emission, complemented by in situ Fe charge state and proton speed measurements from ACE/SWEPAM-SWICS, to identify the source regions of different solar wind streams. The eclipse observations reveal the ubiquity of open structures, invariably associated with Fe XI emission from $\rm Fe^{10+}$, hence a constant electron temperature, $\rm T_{c}$ = $\rm T_{fexi}$, in the expanding corona. The in situ Fe charge states are found to cluster around $\rm Fe^{10+}$, independently of the 300 to 700 km $\rm s^{-1}$ stream speeds, referred to as the continual solar wind. $\rm Fe^{10+}$ thus yields the fiducial link between the continual solar wind and its $\rm T_{fexi}$ sources at the Sun. While the spatial distribution of Fe XIV emission, from $\rm Fe^{13+}$, associated with streamers, changes throughout the solar cycle, the sporadic appearance of charge states $> \rm Fe^{11+}$, in situ, exhibits no cycle dependence regardless of speed. These latter streams are conjectured to be released from hot coronal plasmas at temperatures $\ge \rm T_{fexiv}$ within the bulge of streamers and from active regions, driven by the dynamic behavior of prominences magnetically linked to them. The discovery of continual streams of slow, intermediate and fast solar wind, characterized by the same $\rm T_{fexi}$ in the expanding corona, places new constraints on the physical processes shaping the solar wind.


title: A low phase noise microwave source for high performance CPT Rb atomic clock

abstract: Phase noise of the frequency synthesizer is one of the main limitations to the short-term stability of microwave atomic clocks. In this work, we demonstrated a low-noise, simple-architecture microwave frequency synthesizer for a coherent population trapping (CPT) clock. The synthesizer is mainly composed of a 100 MHz oven controlled crystal oscillator (OCXO), a microwave comb generator and a direct digital synthesizer (DDS). The absolute phase noises of 3.417 GHz signal are measured to be -55 dBc/Hz, -81 dBc/Hz, -111 dBc/Hz and -134 dBc/Hz, respectively, for 1 Hz, 10 Hz, 100 Hz and 1 kHz offset frequencies, which shows only 1 dB deterioration at the second harmonic of the modulation frequency of the atomic clock. The estimated frequency stability of intermodulation effect is 4.7*10^{-14} at 1s averaging time, which is about half order of magnitude lower than that of the state-of-the-art CPT Rb clock. Our work offers an alternative microwave synthesizer for high-performance CPT Rb atomic clock.