title: Black widow evolution: magnetic braking by an ablated wind

abstract: Black widows are close binary systems in which a millisecond pulsar is orbited by a companion a few per cent the mass of the sun. It has been suggested that the pulsar's rotationally powered $\gamma$-ray luminosity gradually evaporates the companion, eventually leaving behind an isolated millisecond pulsar. The evaporation efficiency is determined by the temperature $T_{\rm ch}\propto F^{2/3}$ to which the outflow is heated by the flux $F$ on a dynamical time-scale. Evaporation is most efficient for companions that fill their Roche lobes. In this case, the outflow is dominated by a cap around the L1 point with an angle $\theta_g\sim (T_{\rm ch}/T_g)^{1/2}$, and the evaporation time is $t_{\rm evap}=0.46(T_{\rm ch}/T_g)^{-2}\textrm{ Gyr}$, where $T_g>T_{\rm ch}$ is the companion's virial temperature. We apply our model to the observed black widow population, which has increased substantially over the last decade, considering each system's orbital period, companion mass, and pulsar spin-down power. While the original (Fruchter et al. 1988) black widow evaporates its companion on a few Gyr time-scale, direct evaporation on its own is too weak to explain the overall population. We propose instead that the evaporative wind couples to the companion's magnetic field, removes angular momentum from the binary, and maintains stable Roche-lobe overflow. While a stronger wind carries more mass, it also reduces the Alfv\'en radius, making this indirect magnetic braking mechanism less dependent on the flux $t_{\rm mag}\propto t_{\rm evap}^{1/3}$. This reduces the scatter in evolution times of observed systems, thus better explaining the combined black widow and isolated millisecond pulsar populations.


title: Why binary interaction does not necessarily dominate the formation of Wolf-Rayet stars at low metallicity

abstract: Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z.


title: Lieb Schultz Mattis-Type Theorems and Other Non-perturbative Results for Strongly Correlated Systems with Conserved Dipole Moments

abstract: Non-perturbative constraints on many body physics--such as the famous Lieb-Schultz-Mattis theorem--are valuable tools for studying strongly correlated systems. To this end, we present a number of non-perturbative results that constrain the low-energy physics of systems having conserved dipole moments. We find that for these systems, a unique translationally invariant gapped ground state is only possible if the polarization of the system is integer. Furthermore, if a lattice system also has $U(1)$ subsystem charge conservation symmetry, a unique gapped ground state is only possible if the particle filling along these subsystems is integer. We also apply these methods to spin systems, and determine criteria for the existence of a new type of magnetic response plateau in the presence of a non-uniform magnetic field. Finally, we formulate a version of Luttinger's theorem for 1D systems consisting of dipoles.


title: LoCuSS: exploring the connection between local environment, star formation and dust mass in Abell 1758

abstract: We explore the connection between dust and star formation, in the context of environmental effects on galaxy evolution. In particular, we exploit the susceptibility of dust to external processes to assess the influence of dense environment on star-forming galaxies. We have selected cluster Abell 1758 from the Local Cluster Substructure Survey (LoCuSS). Its complex dynamical state is an ideal test-bench to track dust removal and destruction in galaxies due to merger and accretion shocks. We present a systematic panchromatic study (from 0.15 $\rm \mu$m with GALEX to 500 $\rm \mu$m with Herschel) of spectroscopically confirmed star-forming cluster galaxies at intermediate redshift. We observe that the main subclusters (A1758N and A1758S) belong to two separate large-scale structures, with no overlapping galaxy members. Star-forming cluster members are distributed preferentially outside cluster central regions, and are not grouped in substructures. Rather, these galaxies are being funneled towards the main subclusters along separate accretion filaments. Additionally, we present the first study of dust-to-stellar (DTS) mass ratio used as indicator for local environmental influence on galaxy evolution. Star-forming cluster members show lower mean values (32% at 2.4$\rm \sigma$) of DTS mass ratio and lower levels of infrared emission from birth clouds with respect to coeval star-forming field galaxies. This picture is consistent with the majority of star-forming cluster members infalling in isolation. Upon accretion, star-formation is observed to decrease and warm dust is destroyed due to heating from the intracluster medium radiation, ram-pressure stripping and merger shocks.


title: A panchromatic spatially-resolved analysis of nearby galaxies -- I. Sub-kpc scale Main Sequence in grand-design spirals

abstract: We analyse the spatially resolved relation between stellar mass (M$_{\star}$) and star formation rate (SFR) in disk galaxies (i.e. the Main Sequence, MS). The studied sample includes eight nearby face-on grand-design spirals, e.g. the descendant of high-redshift, rotationally-supported star-forming galaxies. We exploit photometric information over 23 bands, from the UV to the far-IR, from the publicly available DustPedia database to build spatially resolved maps of stellar mass and star formation rates on sub-galactic scales of 0.5-1.5 kpc, by performing a spectral energy distribution fitting procedure that accounts for both the observed and the obscured star formation processes, over a wide range of internal galaxy environments (bulges, spiral arms, outskirts). With more than 30 thousands physical cells, we have derived a definition of the local spatially resolved MS per unit area for disks, $\log(\Sigma_{SFR})$=0.82log$(\Sigma_{*})$-8.69. This is consistent with the bulk of recent results based on optical IFU, using the H$\alpha$ line emission as a SFR tracer. Our work extends the analysis at lower sensitivities in both M$_{\star}$ and SFR surface densities, up to a factor $\sim$ 10. The self consistency of the MS relation over different spatial scales, from sub-galactic to galactic, as well as with a rescaled correlation obtained for high redshift galaxies, clearly proves its universality.


title: Cradles of the first stars: self-shielding, halo masses, and multiplicity

abstract: The formation of Population III (Pop III) stars is a critical step in the evolution of the early universe. To understand how these stars affected their metal-enriched descendants, the details of how, why and where Pop III formation takes place needs to be determined. One of the processes that is assumed to greatly affect the formation of Pop III stars is the presence of a Lyman-Werner (LW) radiation background, that destroys H$_2$, a necessary coolant in the creation of Pop III stars. Self-shielding can alleviate the effect the LW background has on the H$_2$ within haloes. In this work, we perform a cosmological simulation to study the birthplaces of Pop III stars, using the adaptive mesh refinement code Enzo. We investigate the distribution of host halo masses and its relationship to the LW background intensity. Compared to previous work, haloes form Pop III stars at much lower masses, up to a factor of a few, due to the inclusion of H$_2$ self-shielding. We see no relationship between the LW intensity and host halo mass. Most haloes form multiple Pop III stars, with a median number of four, up to a maximum of 16, at the instance of Pop III formation. Our results suggest that Pop III star formation may be less affected by LW radiation feedback than previously thought and that Pop III multiple systems are common.


title: New Selection Rules from Angular Momentum Conservation

abstract: We derive the generalized partial wave expansion for $M \rightarrow N$ scattering amplitude in terms of spinor helicity variables. The basis amplitudes of the expansion with definite angular momentum $j$ consist of the Poincare Clebsch-Gordan coefficients, while $j$ constrains the UV physics that could generate the corresponding operators at tree level. Moreover, we obtain a series of selection rules that restrict the anomalous dimension matrix of effective operators and the way how effective operators contribute to some $2 \rightarrow N$ amplitudes at the loop level.


title: Wind morphology around cool evolved stars in binaries: the case of slowly accelerating oxygen-rich outflows

abstract: The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars.


title: Two delays in white dwarf evolution revealed by Gaia

abstract: By comparing two age indicators of high-mass white dwarfs derived from Gaia data, two discoveries have been made recently: one is the existence of a cooling anomaly that produces the Q branch structure on the Hertzsprung--Russell diagram, the other is the existence of double-white-dwarf merger products. The former poses a challenge for white dwarf cooling models, and the latter has implications on binary evolution and type-Ia supernovae.


title: On the Replicability of Combining Word Embeddings and Retrieval Models

abstract: We replicate recent experiments attempting to demonstrate an attractive hypothesis about the use of the Fisher kernel framework and mixture models for aggregating word embeddings towards document representations and the use of these representations in document classification, clustering, and retrieval. Specifically, the hypothesis was that the use of a mixture model of von Mises-Fisher (VMF) distributions instead of Gaussian distributions would be beneficial because of the focus on cosine distances of both VMF and the vector space model traditionally used in information retrieval. Previous experiments had validated this hypothesis. Our replication was not able to validate it, despite a large parameter scan space.


title: Comparing the orthogonal and unitary functor calculi

abstract: The orthogonal and unitary calculi give a method to study functors from the category of real or complex inner product spaces to the category of based topological spaces. We construct functors between the calculi from the complexification-realification adjunction between real and complex inner product spaces. These allow for movement between the versions of calculi, and comparisons between the Taylor towers produced by both calculi. We show that when the inputted orthogonal functor is weakly polynomial, the Taylor tower of the functor restricted through realification and the restricted Taylor tower of the functor agree up to weak equivalence. We further lift the homotopy level comparison of the towers to a commutative diagram of Quillen functors relating the model categories for orthogonal calculus and the model categories for unitary calculus.


title: Quantifying Computational Advantage of Grover's Algorithm with the Trace Speed

abstract: Despite intensive research, the physical origin of the speed-up offered by quantum algorithms remains mysterious. No general physical quantity, like, for instance, entanglement, can be singled out as the essential useful resource. Here we report a close connection between the trace speed and the quantum speed-up in Grover's search algorithm implemented with pure and pseudo-pure states. For a noiseless algorithm, we find a one-to-one correspondence between the quantum speed-up and the polarization of the pseudo-pure state, which can be connected to a wide class of quantum statistical speeds. For time-dependent partial depolarization and for interrupted Grover searches, the speed-up is specifically bounded by the maximal trace speed that occurs during the algorithm operations. Our results quantify the quantum speed-up with a physical resource that is experimentally measurable and related to multipartite entanglement and quantum coherence.


title: Are universal "anomalous" properties of glasses at low temperatures truly universal?

abstract: The specific heat Cp and other properties of glasses (ranging from amorphous solids to disordered crystals) at low temperatures, are well known to be markedly different from those in fully-ordered crystals. For decades, this qualitative, and even quantitative, universal behavior of glasses has been thoroughly studied. However, a clear understanding of its origin and microscopic nature, needless to say a closed theory, is still lacking. To shed light on this matter, I review the situation in this work, mainly by compiling and discussing measured low-temperature Cp data of many glasses and disordered crystals, as well as highlighting a few exceptions to that "universality rule". Thus, one can see that, in contrast to other low-temperature properties of glasses, the magnitude of the "glassy" Cp excess at low temperature is far from being universal. Even worse, some molecular crystals without a clear sign of disorder exhibit linear coefficients in Cp larger than those found in many amorphous solids, whereas a few of the latter show negligible values.


title: Deep Residual Dense U-Net for Resolution Enhancement in Accelerated MRI Acquisition

abstract: Typical Magnetic Resonance Imaging (MRI) scan may take 20 to 60 minutes. Reducing MRI scan time is beneficial for both patient experience and cost considerations. Accelerated MRI scan may be achieved by acquiring less amount of k-space data (down-sampling in the k-space). However, this leads to lower resolution and aliasing artifacts for the reconstructed images. There are many existing approaches for attempting to reconstruct high-quality images from down-sampled k-space data, with varying complexity and performance. In recent years, deep-learning approaches have been proposed for this task, and promising results have been reported. Still, the problem remains challenging especially because of the high fidelity requirement in most medical applications employing reconstructed MRI images. In this work, we propose a deep-learning approach, aiming at reconstructing high-quality images from accelerated MRI acquisition. Specifically, we use Convolutional Neural Network (CNN) to learn the differences between the aliased images and the original images, employing a U-Net-like architecture. Further, a micro-architecture termed Residual Dense Block (RDB) is introduced for learning a better feature representation than the plain U-Net. Considering the peculiarity of the down-sampled k-space data, we introduce a new term to the loss function in learning, which effectively employs the given k-space data during training to provide additional regularization on the update of the network weights. To evaluate the proposed approach, we compare it with other state-of-the-art methods. In both visual inspection and evaluation using standard metrics, the proposed approach is able to deliver improved performance, demonstrating its potential for providing an effective solution.


title: On the Computational Viability of Quantum Optimization for PMU Placement

abstract: Using optimal phasor measurement unit placement as a prototypical problem, we assess the computational viability of the current generation D-Wave Systems 2000Q quantum annealer for power systems design problems. We reformulate minimum dominating set for the annealer hardware, solve the reformulation for a standard set of IEEE test systems, and benchmark solution quality and time to solution against the CPLEX Optimizer and simulated annealing. For some problem instances the 2000Q outpaces CPLEX. For instances where the 2000Q underperforms with respect to CPLEX and simulated annealing, we suggest hardware improvements for the next generation of quantum annealers.


title: A Blockchain-Encryption-Based Approach to Protect Fog Federations from Rogue Nodes

abstract: People have used cloud computing approach to store their data remotely. As auspicious as this approach is, it brings forth many challenges: from data security to time latency issues with data computation as well as delivery to end users. Fog computing has emerged as an extension for cloud computing to bring data processing and storage close to end-users; however, it minimizes the time latency issue but still suffers from data security challenges. For instance, when a fog node providing services to end users is compromised, the users' data security can be violated. Thus, this paper proposes a secure and fine-grained data access control scheme by integrating the CP-ABE algorithm and blockchain concept to prevent fog nodes from violating end users' data security in a situation where a compromised fog node is being ousted. We also classify the fog nodes into fog federations, based on their attributes such as services and locations, to minimize the time latency and communication overhead between fog nodes and cloud server. Further, the exploitation and integration of the blockchain concept and the CP-ABE algorithm enables fog nodes in the same fog federation to perform the authorization process in a distributed manner. In addition, to solve time latency and communication overhead problems, we equip every fog node with an off-chain database to store most frequently accessed data files for specific time, and with an on-chain access control policies table (On-chain Files Tracking Table) which must be protected from being tampered by malicious (rogue) fog nodes. Therefore, blockchain plays a vital role here as it is tamper-proof by nature. We demonstrate our scheme's efficiency and feasibility by designing algorithms and conducting a security analysis. The provided analysis shows that the proposed scheme is efficient and feasible in ousting malicious (rogue) fog nodes.


title: Regularization of the cosmological sector of loop quantum gravity with bosonic matter and the related problems with general covariance of quantum corrections

abstract: This article concerns the problems regarding different lattice regularization techniques for the matter fields of Hamiltonian constraints defined in the framework of loop quantum gravity. The analysis is formulated in the phase space-reduced cosmological model of the hypothetical theory of canonical quantum general relativity. This article explains why a different than links-related lattice smearing of fields leads to a local violation of general covariance. This happens by assuming, for instance, the nodes-related smearing. Therefore, this problem occurs in the case of any polymerlike scalar field quantization method by breaking the background independence of the semiclassical predictions. In consequence, the diffeomorphism symmetry that depends on a links distribution is broken locally at the level of generally relativistic corrections. Moreover, by using the phase space-reduced gauge fixing technique to analyze this issue, the results are general and they concern any coupling with the links-regularized gravitational degrees of freedom in loop quantum gravity. Therefore, they lead to the following no-go conclusion. Any lattice smearing of matter, not defined by using the geometrical distribution specified by the links-fluxes duality, violates the general principle of relativity.


title: On dark atoms, massive dark photons and millicharged sub-components

abstract: We present a simple model of two dark matter species with opposite millicharge that can form electrically neutral bound states via the exchange of a massive dark photon. If bound state formation is suppressed at low temperatures, a sub-dominant fraction of millicharged particles remains at late times, which can give rise to interesting features in the 21 cm absorption profile at cosmic dawn. The dominant neutral component, on the other hand, can have dipole interactions with ordinary matter, leading to non-standard signals in direct detection experiments. We identify the parameter regions predicting a percent-level ionisation fraction and study constraints from laboratory searches for dark matter scattering and dark photon decays.


title: Broadband variability and correlation study of 3C 279 during flare of 2017-2018

abstract: A multiwavelength temporal and spectral analysis of flares of 3C 279 during November 2017--July 2018 are presented in this work. Three bright gamma-ray flares were observed simultaneously in X-ray and Optical/UV along with a prolonged quiescent state. A "harder-when-brighter" trend is observed in both gamma-rays and X-rays during the flaring period. The gamma-ray light curve for all the flares are binned in one-day time bins and a day scale variability is observed. Variability time constrains the size and location of the emission region to 2.1$\times$10$^{16}$ cm and 4.4$\times$10$^{17}$ cm, respectively. The fractional variability reveals that the source is more than 100\% variable in gamma-rays and it decreases towards the lower energy. A cross-correlation study of the emission from different wavebands is done using the \textit{DCF} method, which shows a strong correlation between them without any time lags. The zero time lag between different wavebands suggest their co-spatial origin. This is the first time 3C 279 has shown a strong correlation between gamma-rays and X-rays emission with zero time lag. A single zone emission model was adopted to model the multiwavelength SEDs by using the publicly available code GAMERA. The study reveals that a higher jet power in electrons is required to explain the gamma-ray flux during the flaring state, as much as, ten times of that required for the quiescent state. However, more jet power in magnetic field has been observed during the quiescent state compared to the flaring state.


title: Scalar correlation functions for a double-well potential in de Sitter space

abstract: We use the spectral representation of the stochastic Starobinsky-Yokoyama approach to compute correlation functions in de Sitter space for a scalar field with a symmetric or asymmetric double-well potential. The terms in the spectral expansion are determined by the eigenvalues and eigenfunctions of the time-independent Fokker-Planck differential operator, and we solve them numerically. The long-distance asymptotic behaviour is given by the lowest state in the spectrum, but we demonstrate that the magnitude of the coeffients of different terms can be very different, and the correlator can be dominated by different terms at different distances. This can give rise to potentially observable cosmological signatures. In many cases the dominant states in the expansion do not correspond to small fluctuations around a minimum of the potential and are therefore not visible in perturbation theory. We discuss the physical interpretation these states, which can be present even when the potential has only one minimum.


title: Gaussian Null Coordinates for Rotating Charged Black Holes and Conserved Charges

abstract: Motivated by the study of conserved Aretakis charges for a scalar field on the horizon of an extremal black hole, we construct the metrics for certain classes of four-dimensional and five-dimensional extremal rotating black holes in Gaussian null coordinates. We obtain these as expansions in powers of the radial coordinate, up to sufficient order to be able to compute the Aretakis charges. The metrics we consider are for 4-charge black holes in four-dimensional STU supergravity (including the Kerr-Newman black hole in the equal-charge case) and the general 3-charge black holes in five-dimensional STU supergravity. We also investigate the circumstances under which the Aretakis charges of an extremal black hole can be mapped by conformal inversion of the metric into Newman-Penrose charges at null infinity. We show that while this works for four-dimensional static black holes, a simple radial inversion fails in rotating cases because a necessary conformal symmetry of the massless scalar equation breaks down. We also discuss that a massless scalar field in dimensions higher than four does not have any conserved Newman-Penrose charge, even in a static asymptotically flat spacetime.


title: Blaschke-Singular-Outer factorization of free non-commutative functions

abstract: By classical results of Herglotz and F. Riesz, any bounded analytic function in the complex unit disk has a unique inner-outer factorization. Here, a bounded analytic function is called \emph{inner} or \emph{outer} if multiplication by this function defines an isometry or has dense range, respectively, as a linear operator on the Hardy Space, $H^2$, of analytic functions in the complex unit disk with square-summable Taylor series. This factorization can be further refined; any inner function $\theta$ decomposes uniquely as the product of a \emph{Blaschke inner} function and a \emph{singular inner} function, where the Blaschke inner contains all the vanishing information of $\theta$, and the singular inner factor has no zeroes in the unit disk. We prove an exact analog of this factorization in the context of the full Fock space, identified as the \emph{Non-commutative Hardy Space} of analytic functions defined in a certain multi-variable non-commutative open unit disk.


title: Pulsation among TESS A and B stars and the Maia variables

abstract: Classification of over 50000 TESS stars in sectors 1-18 has resulted in the detection of 766 pulsating main sequence B stars as well as over 5000 $\delta$ Scuti, 2300 $\gamma$ Doradus and 114 roAp candidates. It is found that a distinct group of $\beta$ Cephei stars, defined as B stars with frequencies greater than 2.5 d$^{-1}$ and predicted by models, does not exist. Instead, high frequencies are to be found over the whole B-star range, eventually merging with $\delta$ Scuti stars. The cool B stars pulsating in high frequencies are the Maia variables. It is shown that Maia variables are not rapidly-rotating $\beta$ Cephei or SPB stars. In the region where $\beta$ Cephei variables are found, the proportion of pulsating stars is larger and amplitudes are higher and a considerable fraction pulsate in a single mode and low rotation rate. There is no distinct SPB instability region. Stars pulsating sole in low frequencies are found among all B stars. At most, only one-third of B stars appear to pulsate. These results, as well as the fact that a large fraction of A and B stars show rotational modulation, indicate a need for a revision of current ideas regarding stars with radiative envelopes.


title: Working mechanism and behavior of Collison nebulizer

abstract: The Collison nebulizer (as well as its variations) has been widely used for generating fine aerosol droplets of a few microns from liquids of viscosity up to 1000 centipoise. It was originally developed for producing medical aerosols in inhalation therapy, and now has become an important component as pneumatic atomizer in the Aerosol Jet(R) direct-write system for additive manufacturing. Qualitative descriptions of its working mechanism were given in the literature as an expanding high-speed gas jet creates a negative pressure to syphon liquid into the jet stream, where the liquid is subsequently blown into sheets, filaments, and eventually droplets. But quantitative analysis and in-depth understanding have been lacking until rather recently. In this work, we present a logical description of the working mechanism of Collison nebulizer based on OpenFOAM(R) CFD analysis of compressible jet flow in the jet expansion channel. The positive-feedback liquid aspiration mechanism becomes clear by examining the CFD results as the jet expansion channel geometry is varied. As a consequence, the output mist density can be rather insensitive to the liquid viscosity, which is illustrated by a set of experiments with the Collison-type pneumatic atomizer in an Aerosol Jet(R) direct-write system. Thus, an intrinsic self-regulation mechanism is elegantly incorporated in the Collison nebulizer design. As intuitively expected, experimental data also supports the notion of the existence of an upper limit for liquid holdup in the limited space of jet expansion channel; therefore, the output mist density cannot increase indefinitely by increasing the atomization gas flow rate.


title: Are exoplanetesimals differentiated?

abstract: Metals observed in the atmospheres of white dwarfs suggest that many have recently accreted planetary bodies. In some cases, the compositions observed suggest the accretion of material dominantly from the core (or the mantle) of a differentiated planetary body. Collisions between differentiated exoplanetesimals produce such fragments. In this work, we take advantage of the large numbers of white dwarfs where at least one siderophile (core-loving) and one lithophile (rock-loving) species have been detected to assess how commonly exoplanetesimals differentiate. We utilise N-body simulations that track the fate of core and mantle material during the collisional evolution of planetary systems to show that most remnants of differentiated planetesimals retain core fractions similar to their parents, whilst some are extremely core-rich or mantle-rich. Comparison with the white dwarf data for calcium and iron indicates that the data are consistent with a model in which $66^{+4}_{-6}\%$ have accreted the remnants of differentiated planetesimals, whilst $31^{+5}_{-5}\%$ have Ca/Fe abundances altered by the effects of heating (although the former can be as high as $100\%$, if heating is ignored). These conclusions assume pollution by a single body and that collisional evolution retains similar features across diverse planetary systems. These results imply that both collisions and differentiation are key processes in exoplanetary systems. We highlight the need for a larger sample of polluted white dwarfs with precisely determined metal abundances to better understand the process of differentiation in exoplanetary systems.


title: The shape of a seed bank tree

abstract: We derive the asymptotic behavior of the total, active and inactive branch lengths of the seed bank coalescent, when the size of the initial sample grows to infinity. Those random variables have important applications for populations evolving under some seed bank effects, such as plants and bacteria, and for some cases of structured populations like metapopulations. The proof relies on the study of the tree at a stopping time corresponding to the first time that a deactivated lineage reactivates. We also give conditional sampling formulas for the random partition and we study the system at the time of the first deactivation of a lineage. All these results provide a good picture of the different regimes and behaviors of the block-counting process of the seed bank coalescent.


title: Hysteresis and Stabillity

abstract: Hysteresis can be defined from a dynamical systems perspective with respect to equilibrium points. Consequently, hysteresis naturally lends itself as a topic to illustrate and extend concepts in a dynamical systems course. A number of examples exhibiting hysteresis, most motivated by applications, are presented. Although the examples can be used to construct student exercises, specific questions are listed in an appendix. A brief extension on hysteresis in partial differential equations is also included.


title: Does GW190425 require an alternative formation pathway than a fast-merging channel?

abstract: The LIGO/Virgo Scientific Collaboration (LVC) recently announced the detection of a compact object binary merger, GW190425, with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, and individual component masses in the range of about 1.1 to 2.5 $M_{\odot}$. If the constituent compact objects are neutron stars, then the total mass is five standard deviations higher than the mean of $2.66\pm 0.12$ M$_{\odot}$ for Galactic binary neutron stars (BNSs). LVC suggests such massive BNS systems are born from a fast-merging channel, and therefore, their non-detection in the Galaxy to be due to a selection effect. However, we are unable to reconcile the inferred formation efficiency from the reported merger rate, $\mathcal{R}_{\rm GW190425}=460^{+1050}_{-390}$ yr$^{-1}$ Gpc$^{-3}$, with predictions from our own study for fast-merging BNS systems. Moreover, the comparable merger rates of GW190425, and GW170817 are possibly in tension with our results for two reasons: (i) more massive systems are expected to have a lower formation rate, and (ii) fast merging channels should constitute $\lesssim 10\%$ of the total BNS systems if case BB unstable mass transfer is permitted to take place as a formation pathway. We argue that to account for the high merger rate of GW190425 as a BNS system requires: (i) revisiting our understanding of NS formation in supernova explosions, or (ii): that more massive NSs need to be preferentially born with either very weak or very high magnetic fields so that they would be undetectable in the radio surveys. Perhaps massive NSs detected in NS-white dwarf binaries are our clues to the formation path of GW190425 systems.


title: Detection of VHE gamma-ray transients with monitoring facilities

abstract: The observation of Very High Energy gamma rays (VHE, E > 100 GeV) led us to the identification of extremely energetic processes and particle acceleration sites both within our Galaxy and beyond. We expect that VHE facilities, like CTA, will explore these sources with an unprecedented level of detail. However, the transient and unpredictable nature of many important processes requires the development of proper monitoring strategies, to observe them. With this study, we estimate the properties of VHE transients that can be effectively detected by monitoring facilities. We use data collected by the Fermi-LAT instrument, during its monitoring campaign, to select events that are likely associated with VHE emission. We use this sample to estimate the frequency, the luminosity and the time-scales of different transients, focusing on blazar flares and Gamma Ray Bursts (GRBs). We discuss how the balance between Field of View, sensitivity and duty cycle of an observatory affects the likelihood to detect transients that occur at the inferred rates and we conclude describing the contribution that current and near-future monitoring facilities can bring to the identification and study of VHE transient emission.


title: CHIPKIT: An agile, reusable open-source framework for rapid test chip development

abstract: The current trend for domain-specific architectures (DSAs) has led to renewed interest in research test chips to demonstrate new specialized hardware. Tape-outs also offer huge pedagogical value garnered from real hands-on exposure to the whole system stack. However, successful tape-outs demand hard-earned experience, and the design process is time consuming and fraught with challenges. Therefore, custom chips have remained the preserve of a small number of research groups, typically focused on circuit design research. This paper describes the CHIPKIT framework. We describe a reusable SoC subsystem which provides basic IO, an on-chip programmable host, memory and peripherals. This subsystem can be readily extended with new IP blocks to generate custom test chips. We also present an agile RTL development flow, including a code generation tool calledVGEN. Finally, we outline best practices for full-chip validation across the entire design cycle.


title: Fast Generation of Big Random Binary Trees

abstract: random_tree() is a linear time and space C++ implementation able to create trees of up to a billion nodes for genetic programming and genetic improvement experiments. A 3.60GHz CPU can generate more than 18 million random nodes for GP program trees per second.


title: Application of Convolutional Neural Networks to Identify Stellar Feedback Bubbles in CO Emission

abstract: We adopt the deep learning method CASI (Convolutional Approach to Shell Identification) and extend it to 3D (CASI-3D) to identify signatures of stellar feedback in molecular line spectra, such as 13CO. We adopt magneto-hydrodynamics simulations that study the impact of stellar winds in a turbulent molecular cloud as an input to generate synthetic observations. We apply the 3D radiation transfer code radmc-3d to model 13CO (J=1-0) line emission from the simulated clouds. We train two CASI-3d models: ME1 predicts only the position of feedback, while MF predicts the fraction of the mass coming from feedback in each voxel. We adopt 75% of the synthetic observations as the training set and assess the accuracy of the two models with the remaining data. We demonstrate that model ME1 identifies bubbles in simulated data with 95% accuracy, and model MF predicts the bubble mass within 4% of the true value. We use bubbles previously visually identified in Taurus in 13CO to validate the models and show both perform well on the highest confidence bubbles. We apply our two models on the full 98 deg2 FCRAO 13CO survey of the Taurus cloud. Models ME1 and MF predict feedback gas mass of 2894 M and 302 M, respectively. When including a correction factor for missing energy due to the limited velocity range of the 13CO data cube, model ME1 predicts feedback kinetic energies of 4.0*1e46 ergs and 1.5*1e47 ergs with/without subtracting the cloud velocity gradient. Model MF predicts feedback kinetic energy of 9.6*1e45 ergs and 2.8*1e46 ergs with/without subtracting the cloud velocity gradient. Model ME1 predicts bubble locations and properties consistent with previous visually identified bubbles. However, model MF demonstrates that feedback properties computed based on visual identifications are significantly overestimated due to line of sight confusion and contamination from background and foreground gas.


title: Human mortality at extreme age

abstract: We use a combination of extreme value theory, survival analysis and computer-intensive methods to analyze the mortality of Italian and French semi-supercentenarians for whom there are validated records. After accounting for the effects of the sampling frame, there appears to be a constant rate of mortality beyond age 108 years and no difference between countries and cohorts. These findings are consistent with previous work based on the International Database on Longevity and suggest that any physical upper bound for humans is so large that it is unlikely to be approached. There is no evidence of differences in survival between women and men after age 108 in the Italian data and the International Database on Longevity; however survival is lower for men in the French data.


title: CATVI: Conditional and Adaptively Truncated Variational Inference for Hierarchical Bayesian Nonparametric Models

abstract: Current variational inference methods for hierarchical Bayesian nonparametric models can neither characterize the correlation structure among latent variables due to the mean-field setting, nor infer the true posterior dimension because of the universal truncation. To overcome these limitations, we propose the conditional and adaptively truncated variational inference method (CATVI) by maximizing the nonparametric evidence lower bound and integrating Monte Carlo into the variational inference framework. CATVI enjoys several advantages over traditional methods, including a smaller divergence between variational and true posteriors, reduced risk of underfitting or overfitting, and improved prediction accuracy. Empirical studies on three large datasets reveal that CATVI applied in Bayesian nonparametric topic models substantially outperforms competing models, providing lower perplexity and clearer topic-words clustering.


title: Examining the Effects of Emotional Valence and Arousal on Takeover Performance in Conditionally Automated Driving

abstract: In conditionally automated driving, drivers have difficulty in takeover transitions as they become increasingly decoupled from the operational level of driving. Factors influencing takeover performance, such as takeover lead time and the engagement of non-driving related tasks, have been studied in the past. However, despite the important role emotions play in human-machine interaction and in manual driving, little is known about how emotions influence drivers takeover performance. This study, therefore, examined the effects of emotional valence and arousal on drivers takeover timeliness and quality in conditionally automated driving. We conducted a driving simulation experiment with 32 participants. Movie clips were played for emotion induction. Participants with different levels of emotional valence and arousal were required to take over control from automated driving, and their takeover time and quality were analyzed. Results indicate that positive valence led to better takeover quality in the form of a smaller maximum resulting acceleration and a smaller maximum resulting jerk. However, high arousal did not yield an advantage in takeover time. This study contributes to the literature by demonstrating how emotional valence and arousal affect takeover performance. The benefits of positive emotions carry over from manual driving to conditionally automated driving while the benefits of arousal do not.


title: A NECTAr-based upgrade for the Cherenkov cameras of the H.E.S.S. 12-meter telescopes

abstract: The High Energy Stereoscopic System (H.E.S.S.) is one of the three arrays of imaging atmospheric Cherenkov telescopes (IACTs) currently in operation. It is composed of four 12-meter telescopes and a 28-meter one, and is sensitive to gamma rays in the energy range ~30 GeV - 100 TeV. The cameras of the 12-m telescopes recently underwent a substantial upgrade, with the goal of improving their performance and robustness. The upgrade involved replacing all camera components except for the photomultiplier tubes (PMTs). This meant developing new hardware for the trigger, readout, power, cooling and mechanical systems, and new software for camera control and data acquisition. Several novel technologies were employed in the cameras: the readout is built around the new NECTAr digitizer chip, developed for the next generation of IACTs; the camera electronics is fully controlled and read out via Ethernet using a combination of FPGA and embedded ARM computers; the software uses modern libraries such as Apache Thrift, ZMQ and Protocol buffers. This work describes in detail the design and the performance of the upgraded cameras.


title: Stable homotopy groups of spheres: From dimension 0 to 90

abstract: Using techniques in motivic homotopy theory, especially the theorem of Gheorghe, the second and the third author on the isomorphism between motivic Adams spectral sequence for $C{\tau}$ and the algebraic Novikov spectral sequence for $BP_*$, we compute the classical and motivic stable homotopy groups of spheres from dimension 0 to 90, except for some carefully enumerated uncertainties.


title: Unoriented Virtual Khovanov Homology

abstract: The Jones polynomial and Khovanov homology of a classical link are invariants that depend upon an initial choice of orientation for the link. In this paper, we give a Khovanov homology theory for unoriented virtual links. The graded Euler characteristic of this homology is proportional to a similarly-defined unoriented Jones polynomial for virtual links, which is a new invariant in the category of non-classical virtual links. The unoriented Jones polynomial continues to satisfy an important property of the usual one: for classical or even virtual links, the unoriented Jones polynomial evaluated at one is two to the power of the number of components of the link. As part of extending the main results of this paper to non-classical virtual links, a new framework for computing integral Khovanov homology is described that can be efficiently and effectively implemented on a computer. We define an unoriented Lee homology theory for virtual links based upon the unoriented version of Khovanov homology.


title: Donald Duck Holiday Game: A numerical analysis of a Game of the Goose role-playing variant

abstract: The 1996 Donald Duck Holiday Game is a role-playing variant of the historical Game of the Goose, involving characters with unique attributes, event squares, and random event cards. The objective of the game is to reach the camping before any other player does. We develop a Monte Carlo simulation model that automatically plays the game and enables analyzing its key characteristics. We assess the game on various metrics relevant to each playability. Numerical analysis shows that, on average, the game takes between 69 and 123 rounds to complete, depending on the number of players. However, durations over one hour (translated to human play time) occur over 25% of the games, which might reduce the quality of the gaming experience. Furthermore, we show that two characters are about 30% likely to win than the other three, primarily due to being exposed to fewer random events. We argue that the richer narrative of role-playing games may extend the duration for which the game remains enjoyable, such that the metrics cannot directly be compared to those of the traditional Game-of-the-Goose. Based on our analysis, we provide several suggestions to improve the game balance with only slight modifications. In a broader sense, we demonstrate that a basic Monte Carlo simulation suffices to analyze Game-of-the-Goose role-playing variants, verify how they score on criteria that contribute to an enjoyable game, and detect possible anomalies.


title: Realizing split-pulse x-ray photon correlation spectroscopy to measure ultrafast dynamics in complex matter

abstract: Split-pulse x-ray photon correlation spectroscopy has been proposed as one of the unique capabilities made possible with the x-ray free electron lasers. It enables characterization of atomic scale structural dynamics that dictates the macroscopic properties of various disordered material systems. Central to the experimental concept are x-ray optics that are capable of splitting individual coherent femtosecond x-ray pulse into two distinct pulses, introduce an adjustable time delay between them, and then recombine the two pulses at the sample position such that they generate two coherent scattering patterns in rapid succession. Recent developments in such optics showed that, while true 'amplitude splitting' optics at hard x-ray wavelengths remains a technical challenge, wavefront and wavelength splitting are both feasible, able to deliver two micron sized focused beams to the sample with sufficient relative stability. Here, we however show that the conventional approach to speckle visibility spectroscopy using these beam splitting techniques can be problematic, even leading to a decoupling of speckle visibility and material dynamics. In response, we discuss the details of the experimental approaches and data analysis protocols for addressing issues caused by subtle beam dissimilarities for both wavefront and wavelength splitting setups. We also show that in some scattering geometries, the Q-space mismatch can be resolved by using two beams of slightly different incidence angle and slightly different wavelengths at the same time. Instead of measuring the visibility of weak speckle patterns, the time correlation in sample structure is encoded in the 'side band' of the spatial autocorrelation of the summed speckle patterns, and can be retrieved straightforwardly from the experimental data. We demonstrate this with a numerical simulation.


title: Statistical Inference of the Value Function for Reinforcement Learning in Infinite Horizon Settings

abstract: Reinforcement learning is a general technique that allows an agent to learn an optimal policy and interact with an environment in sequential decision making problems. The goodness of a policy is measured by its value function starting from some initial state. The focus of this paper is to construct confidence intervals (CIs) for a policy's value in infinite horizon settings where the number of decision points diverges to infinity. We propose to model the action-value state function (Q-function) associated with a policy based on series/sieve method to derive its confidence interval. When the target policy depends on the observed data as well, we propose a SequentiAl Value Evaluation (SAVE) method to recursively update the estimated policy and its value estimator. As long as either the number of trajectories or the number of decision points diverges to infinity, we show that the proposed CI achieves nominal coverage even in cases where the optimal policy is not unique. Simulation studies are conducted to back up our theoretical findings. We apply the proposed method to a dataset from mobile health studies and find that reinforcement learning algorithms could help improve patient's health status. A Python implementation of the proposed procedure is available at https://github.com/shengzhang37/SAVE.


title: Sub-femtometer scale color charge correlations in the proton

abstract: Color charge correlations in the proton at moderately small $x\sim 0.1$ are extracted from its light-cone wave function. The charge fluctuations are far from Gaussian and they exhibit interesting dependence on impact parameter and on the relative transverse momentum (or distance) of the gluon probes. We provide initial conditions for small-$x$ Balitsky-Kovchegov evolution of the dipole scattering amplitude with impact parameter and $\hat r \cdot \hat b$ dependence, and with non-zero $C$-odd component due to three-gluon exchange. Lastly, we compute the (forward) Weizsaecker-Williams gluon distributions, including the distribution of linearly polarized gluons, up to fourth order in $A^+$. The correction due to the quartic correlator provides a transverse momentum scale, $q > 0.5$ GeV, for nearly maximal polarization.


title: Packing and covering balls in graphs excluding a minor

abstract: We prove that for every integer $t\ge 1$ there exists a constant $c_t$ such that for every $K_t$-minor-free graph $G$, and every set $S$ of balls in $G$, the minimum size of a set of vertices of $G$ intersecting all the balls of $S$ is at most $c_t$ times the maximum number of vertex-disjoint balls in $S$. This was conjectured by Chepoi, Estellon, and Vax\`es in 2007 in the special case of planar graphs and of balls having the same radius.


title: Cosets of monodromies and quantum representations

abstract: We use geometric methods to show that given any $3$-manifold $M$, and $g$ a sufficiently large integer, the mapping class group $\mathrm{Mod}(\Sigma_{g,1})$ contains a coset of an abelian subgroup of rank $\lfloor \frac{g}{2}\rfloor,$ consisting of pseudo-Anosov monodromies of open-book decompositions in $M.$ We prove a similar result for rank two free cosets of $\mathrm{Mod}(\Sigma_{g,1}).$ These results have applications to a conjecture of Andersen, Masbaum and Ueno about quantum representations of surface mapping class groups. For surfaces with boundary, and large enough genus, we construct cosets of abelian and free subgroups of their mapping class groups consisting of elements that satisfy the conjecture. The mapping tori of these elements are fibered 3-manifolds that satisfy a weak form of the Turaev-Viro invariants volume conjecture.


title: Heteroglossia: In-Situ Story Ideation with the Crowd

abstract: Ideation is essential for creative writing. Many authors struggle to come up with ideas throughout the writing process, yet modern writing tools fail to provide on-the-spot assistance for writers when they get stuck. This paper introduces Heteroglossia, an add-on for Google Docs that allows writers to elicit story ideas from the online crowd using their text editors. Writers can share snippets of their working drafts and ask the crowd to provide follow-up story ideas based on it. Heteroglossia employs a strategy called "role play", where each worker is assigned a fictional character in a story and asked to brainstorm plot ideas from that character's perspective. Our deployment with two experienced story writers shows that Heteroglossia is easy to use and can generate interesting ideas. Heteroglossia allows us to gain insight into how future technologies can be developed to support ideation in creative writing


title: Security Vetting Process of Smart-home Assistant Applications: A First Look and Case Studies

abstract: The popularity of smart-home assistant systems such as Amazon Alexa and Google Home leads to a booming third-party application market (over 70,000 applications across the two stores). While existing works have revealed security issues in these systems, it is not well understood how to help application developers to enforce security requirements. In this paper, we perform a preliminary case study to examine the security vetting mechanisms adopted by Amazon Alexa and Google Home app stores. With a focus on the authentication mechanisms between Alexa/Google cloud and third-party application servers (i.e. endpoints), we show the current security vetting is insufficient as developer mistakes can not be effectively detected and notified. A weak authentication would allow attackers to spoof the cloud to insert/retrieve data into/from the application endpoints. We validate the attack through ethical proof-of-concept experiments. To confirm vulnerable applications have indeed passed the security vetting and entered the markets, we develop a heuristic-based searching method. We find 219 real-world Alexa endpoints that carry the vulnerability, many of which are related to critical applications that control smart home devices and electronic cars. We have notified Amazon and Google about our findings and offered our suggestions to mitigate the issue.


title: Search for an excited lepton that decays via a contact interaction to a lepton and two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: Results are presented from a search for events containing an excited lepton (electron or muon) produced in association with an ordinary lepton of the same flavor and decaying to a lepton and two hadronic jets. Both the production and the decay of the excited leptons are assumed to occur via a contact interaction with a characteristic energy scale $\Lambda$. The branching fraction for the decay mode under study increases with the mass of the excited lepton and is the most sensitive channel for very heavy excited leptons. The analysis uses a sample of proton-proton collisions collected by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 77.4 fb$^{-1}$. The four-body invariant mass of the two lepton plus two jet system is used as the primary discriminating variable. No significant excess of events beyond the expectation for standard model processes is observed. Assuming that $\Lambda$ is equal to the mass of the excited leptons, excited electrons and muons with masses below 5.6 and 5.7 TeV, respectively, are excluded at 95% confidence level. These are the best limits to date.


title: $A$-Numerical radius orthogonality and parallelism of semi-Hilbertian space operators and their applications

abstract: In this paper, we aim to introduce and characterize the concept of numerical radius orthogonality of operators on a complex Hilbert space $\mathcal{H}$ which are bounded with respect to the semi-norm induced by a positive operator $A$ on $\mathcal{H}$. Moreover, a characterization of the $A$-numerical radius parallelism for $A$-rank one operators is proved. As applications of the obtained results, we obtain some $\mathbb{A}$-numerical radius inequalities of operator matrices where $\mathbb{A}$ is the operator diagonal matrix with diagonal entries are positive operator $A$. Some other related results are also investigated.


title: Quasi-Limiting Behavior of Drifted Brownian Motion

abstract: A Quasi-Stationary Distribution (QSD)for a Markov process with an almost surely hit absorbing state is a time-invariant initial distribution for the process conditioned on not being absorbed by any given time. An initial distribution for the process is in the domain of attraction of some QSD $\nu$ if the distribution of the process a time $t$, conditioned not to be absorbed by time $t$ converges to $\nu$. In this work study mostly Brownian motion with constant drift on the half line $[0,\infty)$ absorbed at $0$. Previous work by Martinez et al. identifies all QSDs and provides a nearly complete characterization for their domain of attraction. Specifically, it was shown that if the distribution a well-defined exponential tail (including the case of lighter than any exponential tail), then it is in the domain of attraction of a QSD determined by the exponent. In this work we 1. Obtain a new approach to existing results, explaining the direct relation between a QSD and an initial distribution in its domain of attraction. 2. Study the behavior under a wide class of initial distributions whose tail is heavier than exponential, and obtain no-trivial limits under appropriate scaling.


title: An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation

abstract: In this work, we propose and develop an arbitrary-order adaptive discontinuous Petrov-Galerkin (DPG) method for the nonlinear Grad-Shafranov equation. An ultraweak formulation of the DPG scheme for the equation is given based on a minimal residual method. The DPG scheme has the advantage of providing more accurate gradients compared to conventional finite element methods, which is desired for numerical solutions to the Grad-Shafranov equation. The numerical scheme is augmented with an adaptive mesh refinement approach, and a criterion based on the residual norm in the minimal residual method is developed to achieve dynamic refinement. Nonlinear solvers for the resulting system are explored and a Picard iteration with Anderson acceleration is found to be efficient to solve the system. Finally, the proposed algorithm is implemented in parallel on MFEM using a domain-decomposition approach, and our implementation is general, supporting arbitrary order of accuracy and general meshes. Numerical results are presented to demonstrate the efficiency and accuracy of the proposed algorithm.


title: Live Exploration with Mobile Robots in a Dynamic Ring, Revisited

abstract: The graph exploration problem requires a group of mobile robots, initially placed arbitrarily on the nodes of a graph, to work collaboratively to explore the graph such that each node is eventually visited by at least one robot. One important requirement of exploration is the {\em termination} condition, i.e., the robots must know that exploration is completed. The problem of live exploration of a dynamic ring using mobile robots was recently introduced in [Di Luna et al., ICDCS 2016]. In it, they proposed multiple algorithms to solve exploration in fully synchronous and semi-synchronous settings with various guarantees when $2$ robots were involved. They also provided guarantees that with certain assumptions, exploration of the ring using two robots was impossible. An important question left open was how the presence of $3$ robots would affect the results. In this paper, we try to settle this question in a fully synchronous setting and also show how to extend our results to a semi-synchronous setting. In particular, we present algorithms for exploration with explicit termination using $3$ robots in conjunction with either (i) unique IDs of the robots and edge crossing detection capability (i.e., two robots moving in opposite directions through an edge in the same round can detect each other), or (ii) access to randomness. The time complexity of our deterministic algorithm is asymptotically optimal. We also provide complementary impossibility results showing that there does not exist any explicit termination algorithm for $2$ robots. The theoretical analysis and comprehensive simulations of our algorithm show the effectiveness and efficiency of the algorithm in dynamic rings. We also present an algorithm to achieve exploration with partial termination using $3$ robots in the semi-synchronous setting.


title: Topology-Aware Cooperative Data Protection in Blockchain-Based Decentralized Storage Networks

abstract: From currency to cloud storage systems, the continuous rise of the blockchain technology is moving various information systems towards decentralization. Blockchain-based decentralized storage networks (DSNs) offer significantly higher privacy and lower costs to customers compared with centralized cloud storage associated with specific vendors. Coding is required in order to retrieve data stored on failing components. While coding solutions for centralized storage have been intensely studied, topology-aware coding for heterogeneous DSNs have not yet been discussed. In this paper, we propose a joint coding scheme where each node receives extra protection through the cooperation with nodes in its neighborhood in a heterogeneous DSN with any given topology. As an extension of, which also subsumes, our prior work on coding for centralized cloud storage, our proposed construction preserves desirable properties such as scalability and flexibility in networks with varying topologies.


title: Multi-Robot Formation Control Using Reinforcement Learning

abstract: In this paper, we present a machine learning approach to move a group of robots in a formation. We model the problem as a multi-agent reinforcement learning problem. Our aim is to design a control policy for maintaining a desired formation among a number of agents (robots) while moving towards a desired goal. This is achieved by training our agents to track two agents of the group and maintain the formation with respect to those agents. We consider all agents to be homogeneous and model them as unicycle [1]. In contrast to the leader-follower approach, where each agent has an independent goal, our approach aims to train the agents to be cooperative and work towards the common goal. Our motivation to use this method is to make a fully decentralized multi-agent formation system and scalable for a number of agents.


title: On Demand Solid Texture Synthesis Using Deep 3D Networks

abstract: This paper describes a novel approach for on demand volumetric texture synthesis based on a deep learning framework that allows for the generation of high quality 3D data at interactive rates. Based on a few example images of textures, a generative network is trained to synthesize coherent portions of solid textures of arbitrary sizes that reproduce the visual characteristics of the examples along some directions. To cope with memory limitations and computation complexity that are inherent to both high resolution and 3D processing on the GPU, only 2D textures referred to as "slices" are generated during the training stage. These synthetic textures are compared to exemplar images via a perceptual loss function based on a pre-trained deep network. The proposed network is very light (less than 100k parameters), therefore it only requires sustainable training (i.e. few hours) and is capable of very fast generation (around a second for $256^3$ voxels) on a single GPU. Integrated with a spatially seeded PRNG the proposed generator network directly returns an RGB value given a set of 3D coordinates. The synthesized volumes have good visual results that are at least equivalent to the state-of-the-art patch based approaches. They are naturally seamlessly tileable and can be fully generated in parallel.


title: Rethinking Curriculum Learning with Incremental Labels and Adaptive Compensation

abstract: Like humans, deep networks have been shown to learn better when samples are organized and introduced in a meaningful order or curriculum. Conventional curriculum learning schemes introduce samples in their order of difficulty. This forces models to begin learning from a subset of the available data while adding the external overhead of evaluating the difficulty of samples. In this work, we propose Learning with Incremental Labels and Adaptive Compensation (LILAC), a two-phase method that incrementally increases the number of unique output labels rather than the difficulty of samples while consistently using the entire dataset throughout training. In the first phase, Incremental Label Introduction, we partition data into mutually exclusive subsets, one that contains a subset of the ground-truth labels and another that contains the remaining data attached to a pseudo-label. Throughout the training process, we recursively reveal unseen ground-truth labels in fixed increments until all the labels are known to the model. In the second phase, Adaptive Compensation, we optimize the loss function using altered target vectors of previously misclassified samples. The target vectors of such samples are modified to a smoother distribution to help models learn better. On evaluating across three standard image benchmarks, CIFAR-10, CIFAR-100, and STL-10, we show that LILAC outperforms all comparable baselines. Further, we detail the importance of pacing the introduction of new labels to a model as well as the impact of using a smooth target vector.


title: Simulate forest trees by integrating L-system and 3D CAD files

abstract: In this article, we propose a new approach for simulating trees, including their branches, sub-branches, and leaves. This approach combines the theory of biological development, mathematical models, and computer graphics, producing simulated trees and forest with full geometry. Specifically, we adopt the Lindenmayer process to simulate the branching pattern of trees and modify the available measurements and dimensions of 3D CAD developed object files to create natural looking sub-branches and leaves. Randomization has been added to the placement of all branches, sub branches and leaves. To simulate a forest, we adopt Inhomogeneous Poisson process to generate random locations of trees. Our approach can be used to create complex structured 3D virtual environment for the purpose of testing new sensors and training robotic algorithms. We look forward to applying this approach to test biosonar sensors that mimick bats' fly in the simulated environment.


title: Dynamics of laterally-coupled pairs of spin-VCSELs

abstract: A newly-developed normal mode model of laser dynamics in a generalised array of waveguides is applied to extend the spin-flip model (SFM) to pairs of evanescently-coupled spin-VCSELS. The effect of high birefringence is explored, revealing new dynamics and regions of bistability. It is shown that optical switching of the polarisation states of the lasers may be controlled through the optical pump and that, under certain conditions, the polarisation of one laser may be switched by controlling the intensity and polarisation in the other.


title: The high-energy environment and atmospheric escape of the mini-Neptune K2-18 b

abstract: K2-18 b is a transiting mini-Neptune that orbits a nearby (38 pc) cool M3 dwarf and is located inside its region of temperate irradiation. We report on the search for hydrogen escape from the atmosphere K2-18 b using Lyman-$\alpha$ transit spectroscopy with the Space Telescope Imaging Spectrograph (STIS) instrument installed on the Hubble Space Telescope (HST). We analyzed the time-series of the fluxes of the stellar Lyman-$\alpha$ emission of K2-18 in both its blue- and redshifted wings. We found that the average blueshifted emission of K2-18 decreases by $67\% \pm 18\%$ during the transit of the planet compared to the pre-transit emission, tentatively indicating the presence of H atoms escaping vigorously and being blown away by radiation pressure. This interpretation is not definitive because it relies on one partial transit. Based on the reconstructed Lyman-$\alpha$ emission of K2-18, we estimate an EUV irradiation between $10^1-10^2$ erg s$^{-1}$ cm$^{-2}$ and a total escape rate in the order of $10^8$ g s$^{-1}$. The inferred escape rate suggests that the planet will lose only a small fraction (< 1%) of its mass and retain its volatile-rich atmosphere during its lifetime. More observations are needed to rule out stellar variability effects, confirm the in-transit absorption and better assess the atmospheric escape and high-energy environment of K2-18 b.


title: Street-level Travel-time Estimation via Aggregated Uber Data

abstract: Estimating temporal patterns in travel times along road segments in urban settings is of central importance to traffic engineers and city planners. In this work, we propose a methodology to leverage coarse-grained and aggregated travel time data to estimate the street-level travel times of a given metropolitan area. Our main focus is to estimate travel times along the arterial road segments where relevant data are often unavailable. The central idea of our approach is to leverage easy-to-obtain, aggregated data sets with broad spatial coverage, such as the data published by Uber Movement, as the fabric over which other expensive, fine-grained datasets, such as loop counter and probe data, can be overlaid. Our proposed methodology uses a graph representation of the road network and combines several techniques such as graph-based routing, trip sampling, graph sparsification, and least-squares optimization to estimate the street-level travel times. Using sampled trips and weighted shortest-path routing, we iteratively solve constrained least-squares problems to obtain the travel time estimates. We demonstrate our method on the Los Angeles metropolitan-area street network, where aggregated travel time data is available for trips between traffic analysis zones. Additionally, we present techniques to scale our approach via a novel graph pseudo-sparsification technique.


title: MOBSTER: Identifying Candidate Magnetic O Stars through Rotational Modulation of TESS Photometry

abstract: Being relatively rare, the properties of magnetic O stars are not fully understood. To date fewer than a dozen of these stars have been confirmed, making any inference of their global properties uncertain due to small number statistics. To better understand these objects it is necessary to increase the known sample. The MOBSTER collaboration aims to do this by identifying candidate magnetic O, B, and A stars from the identification of rotational modulation in high-precision photometry from the Transiting Exoplanet Survey Satellite (TESS). Here we discuss the collaboration's efforts to detect rotational modulation in TESS targets to identify candidate magnetic O stars for future spectropolarimetric observations.


title: Fixed Points of the Set-Based Bellman Operator

abstract: Motivated by uncertain parameters encountered in Markov decision processes (MDPs), we study the effect of parameter uncertainty on Bellman operator-based methods. Specifically, we consider a family of MDPs where the cost parameters are from a given compact set. We then define a Bellman operator acting on an input set of value functions to produce a new set of value functions as the output under all possible variations in the cost parameters. Finally we prove the existence of a fixed point of this set-based Bellman operator by showing that it is a contractive operator on a complete metric space.


title: Understanding and mitigating gradient pathologies in physics-informed neural networks

abstract: The widespread use of neural networks across different scientific domains often involves constraining them to satisfy certain symmetries, conservation laws, or other domain knowledge. Such constraints are often imposed as soft penalties during model training and effectively act as domain-specific regularizers of the empirical risk loss. Physics-informed neural networks is an example of this philosophy in which the outputs of deep neural networks are constrained to approximately satisfy a given set of partial differential equations. In this work we review recent advances in scientific machine learning with a specific focus on the effectiveness of physics-informed neural networks in predicting outcomes of physical systems and discovering hidden physics from noisy data. We will also identify and analyze a fundamental mode of failure of such approaches that is related to numerical stiffness leading to unbalanced back-propagated gradients during model training. To address this limitation we present a learning rate annealing algorithm that utilizes gradient statistics during model training to balance the interplay between different terms in composite loss functions. We also propose a novel neural network architecture that is more resilient to such gradient pathologies. Taken together, our developments provide new insights into the training of constrained neural networks and consistently improve the predictive accuracy of physics-informed neural networks by a factor of 50-100x across a range of problems in computational physics. All code and data accompanying this manuscript are publicly available at \url{https://github.com/PredictiveIntelligenceLab/GradientPathologiesPINNs}.


title: Deep convolutional neural networks for multi-planar lung nodule detection: improvement in small nodule identification

abstract: Objective: In clinical practice, small lung nodules can be easily overlooked by radiologists. The paper aims to provide an efficient and accurate detection system for small lung nodules while keeping good performance for large nodules. Methods: We propose a multi-planar detection system using convolutional neural networks. The 2-D convolutional neural network model, U-net++, was trained by axial, coronal, and sagittal slices for the candidate detection task. All possible nodule candidates from the three different planes are combined. For false positive reduction, we apply 3-D multi-scale dense convolutional neural networks to efficiently remove false positive candidates. We use the public LIDC-IDRI dataset which includes 888 CT scans with 1186 nodules annotated by four radiologists. Results: After ten-fold cross-validation, our proposed system achieves a sensitivity of 94.2% with 1.0 false positive/scan and a sensitivity of 96.0% with 2.0 false positives/scan. Although it is difficult to detect small nodules (i.e. < 6 mm), our designed CAD system reaches a sensitivity of 93.4% (95.0%) of these small nodules at an overall false positive rate of 1.0 (2.0) false positives/scan. At the nodule candidate detection stage, results show that a multi-planar method is capable to detect more nodules compared to using a single plane. Conclusion: Our approach achieves good performance not only for small nodules, but also for large lesions on this dataset. This demonstrates the effectiveness and efficiency of our developed CAD system for lung nodule detection. Significance: The proposed system could provide support for radiologists on early detection of lung cancer.


title: Measured velocity spectra and neutron densities of the PF2 ultracold-neutron beam ports at the Institut Laue--Langevin

abstract: Ultracold neutrons (UCNs) are a useful tool for fundamental physics experiments. They can be used to probe the lifetime of free neutrons, search for new CP violating processes and exotic interactions beyond the Standard Model, perform Ramsey spectroscopy, and carry out neutron-optical interference experiments. All of these experiments require high neutron count rates for good statistics. For optimal exploitation of experimental beam time, these experiments need to be prepared and, at times, even simulated in advance. To this end, it is crucial to know the velocity-dependent UCN flux at each beam position. Knowing the absolute neutron flux also allows for an absolute calibration of previously gathered data. Using the same time-of-fight experimental setup, we have measured the differential neutron flux of three out of the four UCN beam ports at the PF2 instrument at Institut Laue--Langevin, Grenoble. These beam ports are commonly used for UCN flux experiments and proof-of-principle tests.


title: High Frequency Microwave Emission of a Tri-layer Magnetic Tunnel Junction in the Absence of External Bias-Magnetic Field

abstract: We perform an experimental study of DC current induced microwave emission in a magnetic tunnel junction (MTJ) consisting of three active magnetic layers. For this tri-layer structure, in addition to a conventional bilayer orthogonal MTJ containing a perpendicular free layer and an in-plane fixed layer, a second perpendicular layer has been introduced. We found that the microwave emission frequency induced by spin-transfer torque (STT) can reach as high as 6 GHz in the absence of any applied magnetic field. Moreover, microwave emission is observed for both current polarizations where a redshift is seen with increase in magnitude of current. We discuss spin-dynamics of the observed bi-directional high-frequency emission and the physical origin of the red-shift. The distinct microwave emission properties exhibited in this tri-layer MTJ structure could potentially be useful for future applications in nanoscale spintronics devices such as microwave communication and neuromorphic computing.


title: Geometric Motion Planning for Affine Control Systems with Indefinite Boundary Conditions and Free Terminal Time

abstract: The problem of motion planning for affine control systems consists of designing control inputs that drive a system from a well-defined initial to final states in a desired amount of time. For control systems with drift, however, understanding which final states are reachable in a given time, or reciprocally the amount of time needed to reach a final state, is often the most difficult part of the problem. We address this issue in this paper and introduce a new method to solve motion planning problems for affine control systems, where the motion desired can have indefinite boundary conditions and the time required to perform the motion is free. The method extends on our earlier work on motion planning for systems without drift. A canonical example of parallel parking of a unicycle with constant linear velocity is provided in this paper to demonstrate our algorithm.


title: Visual Storytelling via Predicting Anchor Word Embeddings in the Stories

abstract: We propose a learning model for the task of visual storytelling. The main idea is to predict anchor word embeddings from the images and use the embeddings and the image features jointly to generate narrative sentences. We use the embeddings of randomly sampled nouns from the groundtruth stories as the target anchor word embeddings to learn the predictor. To narrate a sequence of images, we use the predicted anchor word embeddings and the image features as the joint input to a seq2seq model. As opposed to state-of-the-art methods, the proposed model is simple in design, easy to optimize, and attains the best results in most automatic evaluation metrics. In human evaluation, the method also outperforms competing methods.


title: Neutron imaging of liquid-liquid systems containing paramagnetic salt solutions

abstract: The method of neutron imaging was adopted to map the concentration evolution of aqueous paramagnetic Gd(NO3)3 solutions. Magnetic manipulation of the paramagnetic liquid within a miscible nonmagnetic liquid is possible by countering density-difference driven convection. The formation of salt fingers caused by double-diffusive convection in a liquid-liquid system of Gd(NO3)3 and Y(NO3)3 solutions can be prevented by the magnetic field gradient force.

categories: physics.app-ph cond-mat.other physics.chem-ph physics.flu-dyn


title: Hardware Implementation of Neural Self-Interference Cancellation

abstract: In-band full-duplex systems can transmit and receive information simultaneously on the same frequency band. However, due to the strong self-interference caused by the transmitter to its own receiver, the use of non-linear digital self-interference cancellation is essential. In this work, we describe a hardware architecture for a neural network-based non-linear self-interference (SI) canceller and we compare it with our own hardware implementation of a conventional polynomial based SI canceller. In particular, we present implementation results for a shallow and a deep neural network SI canceller as well as for a polynomial SI canceller. Our results show that the deep neural network canceller achieves a hardware efficiency of up to $312.8$ Msamples/s/mm$^2$ and an energy efficiency of up to $0.9$ nJ/sample, which is $2.1\times$ and $2\times$ better than the polynomial SI canceller, respectively. These results show that NN-based methods applied to communications are not only useful from a performance perspective, but can also be a very effective means to reduce the implementation complexity.


title: Chance Constrained Covariance Control for Linear Stochastic Systems With Output Feedback

abstract: We consider the problem of steering, via output feedback, the state distribution of a discrete-time, linear stochastic system from an initial Gaussian distribution to a terminal Gaussian distribution with prescribed mean and maximum covariance, subject to probabilistic path constraints on the state. The filtered state is obtained via a Kalman filter, and the problem is formulated as a deterministic convex program in terms of the distribution of the filtered state. We observe that, in the presence of constraints on the state covariance, and in contrast to classical Linear Quadratic Gaussian (LQG) control, the optimal feedback control depends on both the process noise and the observation model. The effectiveness of the proposed approach is verified using a numerical example.


title: Private Computation with Individual and Joint Privacy

abstract: This paper considers the problem of single-server Private Computation (PC) in the presence of Side Information (SI). In this problem, there is a server that stores $K$ i.i.d. messages, and a user who has a subset of $M$ uncoded messages or a coded linear combination of them as side information, where the identities of these messages are unknown to the server. The user wants to privately compute (via downloading information from the server) a linear combination of a subset of $D$ other messages, where the identities of these messages must be kept private individually or jointly. For each setting, we define the capacity as the supremum of all achievable download rates. We characterize the capacity of both PC with coded and uncoded SI when individual privacy is required, for all $K, M, D$. Our results indicate that both settings have the same capacity. In addition, we establish a non-trivial lower bound on the capacity of PC with coded SI when joint privacy is required, for a range of parameters $K, M, D$. This lower bound is the same as the lower bound we previously established on the capacity of PC with uncoded SI when joint privacy is required.


title: Irrigable Measures for Weighted Irrigation Plans

abstract: A model of irrigation network, where lower branches must be thicker in order to support the weight of the higher ones, was recently introduced in BSUN2 [4]. This leads to a countable family of ODEs, describing the thickness of every branch, solved by backward induction. The present paper determines what kind of measures can be irrigated with a finite weighted cost. Indeed, the boundedness of the cost depends on the dimension of the support of the irrigated measure, and also on the asymptotic properties of the ODE which determines the thickness of branches.


title: Learning Transformation-Aware Embeddings for Image Forensics

abstract: A dramatic rise in the flow of manipulated image content on the Internet has led to an aggressive response from the media forensics research community. New efforts have incorporated increased usage of techniques from computer vision and machine learning to detect and profile the space of image manipulations. This paper addresses Image Provenance Analysis, which aims at discovering relationships among different manipulated image versions that share content. One of the main sub-problems for provenance analysis that has not yet been addressed directly is the edit ordering of images that share full content or are near-duplicates. The existing large networks that generate image descriptors for tasks such as object recognition may not encode the subtle differences between these image covariates. This paper introduces a novel deep learning-based approach to provide a plausible ordering to images that have been generated from a single image through transformations. Our approach learns transformation-aware descriptors using weak supervision via composited transformations and a rank-based quadruplet loss. To establish the efficacy of the proposed approach, comparisons with state-of-the-art handcrafted and deep learning-based descriptors, and image matching approaches are made. Further experimentation validates the proposed approach in the context of image provenance analysis.


title: Tidal Disruption Flares from Stars on Marginally Bound and Unbound Orbits

abstract: We study the mass fallback rate of tidally disrupted stars on marginally bound and unbound orbits around a supermassive black hole (SMBH) by performing three-dimensional smoothed particle hydrodynamic (SPH) simulations with three key parameters: the star is modeled by a polytrope with two different indexes ($n=1.5$ and $3$). The stellar orbital properties are characterized by five orbital eccentricities ranging from $e=0.98$ to $1.02$ and five different penetration factors ranging from $\beta=1$ to $3$, where $\beta$ represents the ratio of the tidal disruption to pericenter distance radii. We derive analytic formulae for the mass fallback rate as a function of the stellar density profile, orbital eccentricity, and penetration factor. Moreover, two critical eccentricities to classify tidal disruption events (TDEs) into five different types: eccentric ($e<e_{\rm crit,1}$), marginally eccentric ($e_{\rm crit,1}\lesssim{e}<1$), purely parabolic ($e=1$), marginally hyperbolic ($1<e<e_{\rm crit,2}$), and hyperbolic ($e\gtrsim{e_{\rm crit,2}}$) TDEs, are reevaluated as $e_{\rm crit,1}=1-2q^{-1/3}\beta^{k-1}$ and $e_{\rm crit,2}=1+2q^{-1/3}\beta^{k-1}$, where $q$ is the ratio of the SMBH to stellar masses and $0<k\lesssim2$. We find the asymptotic slope of the mass fallback rate varies with the TDE type. The asymptotic slope approaches $-5/3$ for parabolic TDEs, is steeper for marginally eccentric TDEs, and is flatter for marginally hyperbolic TDEs. The mass fallback rates of the marginally eccentric TDEs are much larger than the parabolic TDE case, while those of marginally hyperbolic TDEs are much smaller. Marginally unbound TDEs could be an origin of a very low density gas disk around a dormant SMBH.


title: Solution sets of systems of equations over finite lattices and semilattices

abstract: Solution sets of systems of homogeneous linear equations over fields are characterized as being subspaces, i.e., sets that are closed under linear combinations. Our goal is to characterize solution sets of systems of equations over arbitrary finite algebras by a similar closure condition. We show that solution sets are always closed under the centralizer of the clone of term operations of the given algebra; moreover, the centralizer is the only clone that could characterize solution sets. If every centralizer-closed set is the set of all solutions of a system of equations over a finite algebra, then we say that the algebra has Property (SDC). Our main result is the description of finite lattices and semilattices with Property (SDC): we prove that a finite lattice has Property (SDC) if and only if it is a Boolean lattice, and a finite semilattice has Property (SDC) if and only if it is distributive.


title: ExoMol molecular line lists -- XXXVII: spectra of acetylene

abstract: A new ro-vibrational line list for the ground electronic state of the main isotopologue of acetylene, $^{12}$C$_2$H$_2$, is computed as part of the ExoMol project. The aCeTY line list covers the transition wavenumbers up to 10,000 cm$^{-1}$ ($ \lambda >1$ $\mu$m), with lower and upper energy levels up to 12,000 cm$^{-1}$ and 22,000 cm$^{-1}$ considered, respectively. The calculations are performed up to a maximum value for the vibrational angular momentum, $K_{\rm max}=L_{\rm max}$ = 16, and maximum rotational angular momentum, $J$ = 99. Higher values of $J$ were not within the specified wavenumber window. The aCeTY line list is considered to be complete up to 2200 K, making it suitable for use in characterising high-temperature exoplanet or cool stellar atmospheres. Einstein-A coefficients, which can directly be used to calculate intensities at a particular temperature, are computed for 4.3 billion (4,347,381,911) transitions between 5 million (5,160,803) energy levels. We make comparisons against other available data for $^{12}$C$_2$H$_2$, and demonstrate this to be the most complete line list available. The line list is available in electronic form from the online CDS and ExoMol databases.


title: Finite strain constitutive modeling for shape memory alloys considering transformation-induced plasticity and two-way shape memory effect

abstract: This work presents a three-dimensional constitutive model for shape memory alloys considering the TRansformation-Induced Plasticity (TRIP) as well as the Two-Way Shape Memory Effect (TWSME) through a large deformation framework. The presented logarithmic strain based model is able to capture the large strains and rotations exhibited by SMAs under general thermomechanical cycling. By using the martensitic volume fraction, transformation strain, internal stress, and TRIP strain tensors as internal state variables, the model is capable to capture the stress-dependent TRIP generation when SMAs are subjected to a multiaxial stress state, as well as the TWSME for thermomechanically trained SMAs under load-free conditions. A detailed implementation procedure of the proposed model is presented through a user-defined material subroutine within a finite element framework allowing for solving different Boundary Value Problems (BVPs). Comprehensive instruction on calibrating the model parameters as well as the derivation of continuum tangent stiffness matrix are also provided. In the end, the simulated cyclic pseudoelastic and actuation responses by the presented model for a wide range of SMA material systems under both uniaxial and multiaxial stress states are compared against experimental results to validate the proposed modeling capabilities.


title: Deep Learning Stereo Vision at the edge

abstract: We present an overview of the methodology used to build a new stereo vision solution that is suitable for System on Chip. This new solution was developed to bring computer vision capability to embedded devices that live in a power constrained environment. The solution is constructured as a hybrid between classical Stereo Vision techniques and deep learning approaches. The stereoscopic module is composed of two separate modules: one that accelerates the neural network we trained and one that accelerates the front-end part. The system is completely passive and does not require any structured light to obtain very compelling accuracy. With respect to the previous Stereo Vision solutions offered by the industries we offer a major improvement is robustness to noise. This is mainly possible due to the deep learning part of the chosen architecture. We submitted our result to Middlebury dataset challenge. It currently ranks as the best System on Chip solution. The system has been developed for low latency applications which require better than real time performance on high definition videos.


title: Two-body weak decays of doubly charmed baryons

abstract: The hadronic two-body weak decays of the doubly charmed baryons $\Xi_{cc}^{++}, \Xi_{cc}^+$ and $\Omega_{cc}^+$ are studied in this work. To estimate the nonfactorizable contributions, we work in the pole model for the $P$-wave amplitudes and current algebra for $S$-wave ones. For the $\Xi_{cc}^{++}\to \Xi_c^+\pi^+$ mode, we find a large destructive interference between factorizable and nonfactorizable contributions for both $S$- and $P$-wave amplitudes. Our prediction of $\sim 0.70\%$ for its branching fraction is smaller than the earlier estimates in which nonfactorizable effects were not considered, but agrees nicely with the result based on an entirely different approach, namely, the covariant confined quark model. On the contrary, a large constructive interference was found in the $P$-wave amplitude by Dhir and Sharma, leading to a branching fraction of order $(7-16)\%$. Using the current results for the absolute branching fractions of $(\Lambda_c^+,\Xi_c^+)\to p K^-\pi^+$ and the LHCb measurement of $\Xi_{cc}^{++}\to\Xi_c^+\pi^+$ relative to $\Xi_{cc}^{++}\to\Lambda_c^+ K^- \pi^+\pi^+$, we obtain $\B(\Xi_{cc}^{++}\to\Xi_c^+\pi^+)_{\rm expt}\approx (1.83\pm1.01)\%$ after employing the latest prediction of $\B(\Xi_{cc}^{++}\to\Sigma_c^{++}\overline{K}^{*0})$. Our prediction of $\mathcal{B}(\Xi_{cc}^{++}\to\Xi_c^+\pi^+)\approx 0.7\%$ is thus consistent with the experimental value but in the lower end. It is important to pin down the branching fraction of this mode in future study. Factorizable and nonfactorizable $S$-wave amplitudes interfere constructively in $\Xi_{cc}^+\to\Xi_c^0\pi^+$. Its large branching fraction of order 4\% may enable experimentalists to search for the $\Xi_{cc}^+$ through this mode. That is, the $\Xi_{cc}^+$ is reconstructed through the $\Xi_{cc}^+\to\Xi_c^0\pi^+$ followed by the decay chain $\Xi_c^0\to \Xi^-\pi^+\to p\pi^-\pi^-\pi^+$.


title: On the Pilot Contamination Attack in Multi-Cell Multiuser Massive MIMO Networks

abstract: In this paper, we analyze pilot contamination (PC) attacks on a multi-cell massive multiple-input multiple-output (MIMO) network with correlated pilots. We obtain correlated pilots using a user capacity-achieving pilot sequence design. This design relies on an algorithm which designs correlated pilot sequences based on signal-to-interference-plus-noise ratio (SINR) requirements for all the legitimate users. The pilot design is capable of achieving the SINR requirements for all users even in the presence of PC. However, this design has some intrinsic limitations and vulnerabilities, such as a known pilot sequence and the non-zero cross-correlation among different pilot sequences. We reveal that such vulnerabilities may be exploited by an active attacker to increase PC in the network. Motivated by this, we analyze the correlated pilot design for vulnerabilities that can be exploited by an active attacker. Based on this analysis, we develop an effective active attack strategy in the massive MIMO network with correlated pilot sequences. Our examinations reveal that the user capacity region of the network is significantly reduced in the presence of the active attack. Importantly, the SINR requirements for the worst-affected users may not be satisfied even with an infinite number of antennas at the base station.