title: Three Efficient, Low-Complexity Algorithms for Automatic Color Trapping

abstract: Color separations (most often cyan, magenta, yellow, and black) are commonly used in printing to reproduce multi-color images. For mechanical reasons, these color separations are generally not perfectly aligned with respect to each other when they are rendered by their respective imaging stations. This phenomenon, called color plane misregistration, causes gap and halo artifacts in the printed image. Color trapping is an image processing technique that aims to reduce these artifacts by modifying the susceptible edge boundaries to create small, unnoticeable overlaps between the color planes. We propose three low-complexity algorithms for automatic color trapping which hide the effects of small color plane mis-registrations. Our algorithms are designed for software or embedded firmware implementation. The trapping method they follow is based on a hardware-friendly technique proposed by J. Trask (JTHBCT03) which is too computationally expensive for software or firmware implementation. The first two algorithms are based on the use of look-up tables (LUTs). The first LUT-based algorithm corrects all registration errors of one pixel in extent and reduces several cases of misregistration errors of two pixels in extent using only 727 Kbytes of storage space. This algorithm is particularly attractive for implementation in the embedded firmware of low-cost formatter-based printers. The second LUT-based algorithm corrects all types of misregistration errors of up to two pixels in extent using 3.7 Mbytes of storage space. The third algorithm is a hybrid one that combines LUTs and feature extraction to minimize the storage requirements (724 Kbytes) while still correcting all misregistration errors of up to two pixels in extent. This algorithm is suitable for both embedded firmware implementation on low-cost formatter-based printers and software implementation on host-based printers.


title: Angle Feedback for NOMA Transmission in mmWave Drone Networks

abstract: In this paper, we consider an unmanned aerial vehicle (UAV) based wireless network using non-orthogonal multiple access (NOMA) transmission in millimeter-wave frequencies to deliver broadband data in a spectrally efficient fashion at hotspot scenarios. The necessity for the NOMA transmitter to gather information on user channel quality becomes a major drawback in practical deployments. We therefore consider various limited feedback schemes for NOMA transmission, to relieve the complexity of tracking and feeding back the full channel state information (CSI) of the users. In particular, through beamforming we allow NOMA to exploit the space domain, and hence the user angle emerges as a promising (yet novel) limited feedback scheme. We show that as the user region for NOMA transmission gets wider, the users become more distinctive at the transmitter side with respect to their angles, making user angle feedback a better alternative than distance feedback in such scenarios. We rigorously derive and analyze the outage sum rate performance for NOMA transmission considering various user ordering strategies involving full CSI, angle, and distance feedback schemes. Our analytical results for NOMA outage sum rates using those feedback schemes match closely with simulations, and provide useful insights on properly choosing a limited feedback scheme for different deployment geometries and operating configurations.


title: On the uniqueness of the maximum parsimony tree for data with up to two substitutions: an extension of the classic Buneman theorem in phylogenetics

abstract: One of the main aims of phylogenetics is the reconstruction of the correct evolutionary tree when data concerning the underlying species set are given. These data typically come in the form of DNA, RNA or protein alignments, which consist of various characters (also often referred to as sites). Often, however, tree reconstruction methods based on criteria like maximum parsimony may fail to provide a unique tree for a given dataset, or, even worse, reconstruct the `wrong' tree (i.e. a tree that differs from the one that generated the data). On the other hand it has long been known that if the alignment consists of all the characters that correspond to edges of a particular tree, i.e. they all require exactly $k=1$ substitution to be realized on that tree, then this tree will be recovered by maximum parsimony methods. This is based on Buneman's theorem in mathematical phylogenetics. It is the goal of the present manuscript to extend this classic result as follows: We prove that if an alignment consists of all characters that require exactly $k=2$ substitutions on a particular tree, this tree will always be the unique maximum parsimony tree (and we also show that this can be generalized to characters which require at most $k=2$ substitutions). In particular, this also proves a conjecture based on a recently published observation by Goloboff et al. affirmatively for the special case of $k=2$.


title: Millimeter-wave Extended NYUSIM Channel Model for Spatial Consistency

abstract: Commonly used drop-based channel models cannot satisfy the requirements of spatial consistency for millimeter-wave (mmWave) channel modeling where transient motion or closely-spaced users need to be considered. A channel model having \textit{spatial consistency} can capture the smooth variations of channels, when a user moves, or when multiple users are close to each other in a local area within, say, 10 m in an outdoor scenario. Spatial consistency is needed to support the testing of beamforming and beam tracking for massive multiple-input and multiple-output (MIMO) and multi-user MIMO in fifth-generation (5G) mmWave mobile networks. This paper presents a channel model extension and an associated implementation of spatial consistency in the NYUSIM channel simulation platform. Along with a mathematical model, we use measurements where the user moved along a street and turned at a corner over a path length of 75 m in order to derive realistic values of several key parameters such as correlation distance and the rate of cluster birth and death, that are shown to provide spatial consistency for NYUSIM in an urban microcell street canyon scenario.


title: Smoothed Hinge Loss and $\ell^{1}$ Support Vector Machines

abstract: A new algorithm is presented for solving the soft-margin Support Vector Machine (SVM) optimization problem with an $\ell^{1}$ penalty. This algorithm is designed to require a modest number of passes over the data, which is an important measure of its cost for very large data sets. The algorithm uses smoothing for the hinge-loss function, and an active set approach for the $\ell^{1}$ penalty.


title: A multi-wavelength view of the neutron star binary 1FGL J1417.7--4402: A progenitor to canonical millisecond pulsars

abstract: The Fermi $\gamma$-ray source 1FGL J1417.7--4407 (J1417) is a compact X-ray binary with a neutron star primary and a red giant companion in a $\sim$5.4 day orbit. This initial conclusion, based on optical and X-ray data, was confirmed when a 2.66 ms radio pulsar was found at the same location (and with the same orbital properties) as the optical/X-ray source. However, these initial studies found conflicting evidence about the accretion state and other properties of the binary. We present new optical, radio, and X-ray observations of J1417 that allow us to better understand this unusual system. We show that one of the main pieces of evidence previously put forward for an accretion disk---the complex morphology of the persistent H$\alpha$ emission line---can be better explained by the presence of a strong, magnetically driven stellar wind from the secondary and its interaction with the pulsar wind. The radio spectral index derived from VLA/ATCA observations is broadly consistent with that expected from a millisecond pulsar, further disfavoring an accretion disk scenario. X-ray observations show evidence for a double-peaked orbital light curve, similar to that observed in some redback millisecond pulsar binaries and likely due to an intrabinary shock. Refined optical light curve fitting gives a distance of 3.1$\pm$0.6 kpc, confirmed by a Gaia DR2 parallax measurement. At this distance the X-ray luminosity of J1417 is (1.0$^{+0.4}_{-0.3}$) $\times 10^{33}$ erg s$^{-1}$, which is more luminous than all known redback systems in the rotational-powered pulsar state, perhaps due to the wind from the giant companion. The unusual phenomenology of this system and its differing evolutionary path from redback millisecond pulsar binaries points to a new eclipsing pulsar "spider" subclass that is a possible progenitor of normal field millisecond pulsar binaries.


title: A Tutorial on Clique Problems in Communications and Signal Processing

abstract: Since its first use by Euler on the problem of the seven bridges of K\"onigsberg, graph theory has shown excellent abilities in solving and unveiling the properties of multiple discrete optimization problems. The study of the structure of some integer programs reveals equivalence with graph theory problems making a large body of the literature readily available for solving and characterizing the complexity of these problems. This tutorial presents a framework for utilizing a particular graph theory problem, known as the clique problem, for solving communications and signal processing problems. In particular, the paper aims to illustrate the structural properties of integer programs that can be formulated as clique problems through multiple examples in communications and signal processing. To that end, the first part of the tutorial provides various optimal and heuristic solutions for the maximum clique, maximum weight clique, and $k$-clique problems. The tutorial, further, illustrates the use of the clique formulation through numerous contemporary examples in communications and signal processing, mainly in maximum access for non-orthogonal multiple access networks, throughput maximization using index and instantly decodable network coding, collision-free radio frequency identification networks, and resource allocation in cloud-radio access networks. Finally, the tutorial sheds light on the recent advances of such applications, and provides technical insights on ways of dealing with mixed discrete-continuous optimization problems.


title: Toward a Measurement of the Transverse Peculiar Velocity of Galaxy Pairs

abstract: The transverse peculiar velocities caused by the mass distribution of large-scale structure provide a test of the theoretical matter power spectrum and the cosmological parameters that contribute to its shape. Typically, the matter density distribution of the nearby Universe is measured through redshift or line-of-sight peculiar velocity surveys. However, both methods require model-dependent distance measures to place the galaxies or to differentiate peculiar velocity from the Hubble expansion. In this paper, we use the correlated proper motions of galaxy pairs from the VLBA Extragalactic Proper Motion Catalog to place limits on the transverse peculiar velocity of galaxy pairs with comoving separations <1500 Mpc without a reliance on precise distance measurements. The relative proper motions of galaxy pairs across the line of sight can be directly translated into relative peculiar velocities because no proper motion will occur in a homogeneous expansion. We place a 3 sigma limit on the relative proper motion of pairs with comoving separations < 100 Mpc of -17.4 microas/yr < thetadot / sin theta < 19.8 microas/yr. We also confirm that large-separation objects (> 200 Mpc) are consistent with pure Hubble expansion to within ~ 5.3 microas/yr (1 sigma). Finally, we predict that Gaia end-of-mission proper motions will be able to significantly detect the mass distribution of large-scale structure on length scales < 25 Mpc. This future detection will allow a test of the shape of the theoretical mass power spectrum without a reliance on precise distance measurements.


title: Aiming to Know You Better Perhaps Makes Me a More Engaging Dialogue Partner

abstract: There have been several attempts to define a plausible motivation for a chit-chat dialogue agent that can lead to engaging conversations. In this work, we explore a new direction where the agent specifically focuses on discovering information about its interlocutor. We formalize this approach by defining a quantitative metric. We propose an algorithm for the agent to maximize it. We validate the idea with human evaluation where our system outperforms various baselines. We demonstrate that the metric indeed correlates with the human judgments of engagingness.


title: Non-asymptotic bounds for sampling algorithms without log-concavity

abstract: Discrete time analogues of ergodic stochastic differential equations (SDEs) are one of the most popular and flexible tools for sampling high-dimensional probability measures. Non-asymptotic analysis in the $L^2$ Wasserstein distance of sampling algorithms based on Euler discretisations of SDEs has been recently developed by several authors for log-concave probability distributions. In this work we replace the log-concavity assumption with a log-concavity at infinity condition. We provide novel $L^2$ convergence rates for Euler schemes, expressed explicitly in terms of problem parameters. From there we derive non-asymptotic bounds on the distance between the laws induced by Euler schemes and the invariant laws of SDEs, both for schemes with standard and with randomised (inaccurate) drifts. We also obtain bounds for the hierarchy of discretisation, which enables us to deploy a multi-level Monte Carlo estimator. Our proof relies on a novel construction of a coupling for the Markov chains that can be used to control both the $L^1$ and $L^2$ Wasserstein distances simultaneously. Finally, we provide a weak convergence analysis that covers both the standard and the randomised (inaccurate) drift case. In particular, we reveal that the variance of the randomised drift does not influence the rate of weak convergence of the Euler scheme to the SDE.


title: Convergence in High Probability of the Quantum Diffusion in a Random Band Matrix Model

abstract: We consider Hermitian random band matrices $H$ in $d \geq 1 $ dimensions. The matrix elements $H_{xy},$ indexed by $x, y \in \Lambda \subset \mathbb{Z}^d,$ are independent, uniformly distributed random variable if $|x-y| $ is less than the band width $W,$ and zero otherwise. We update the previous results of the converge of quantum diffusion in a random band matrix model from convergence of the expectation to convergence in high probability. The result is uniformly in the size $|\Lambda| $ of the matrix.


title: Limit order books, diffusion approximations and reflected SPDEs: from microscopic to macroscopic models

abstract: Motivated by a zero-intelligence approach, the aim of this paper is to connect the microscopic (discrete price and volume), mesoscopic (discrete price and continuous volume) and macroscopic (continuous price and volume) frameworks for the modelling of limit order books, with a view to providing a natural probabilistic description of their behaviour in a high to ultra high-frequency setting. Starting with a microscopic framework, we first examine the limiting behaviour of the order book process when order arrival and cancellation rates are sent to infinity and when volumes are considered to be of infinitesimal size. We then consider the transition between this mesoscopic model and a macroscopic model for the limit order book, obtained by letting the tick size tend to zero. The macroscopic limit can then be described using reflected SPDEs which typically arise in stochastic interface models. We then use financial data to discuss a possible calibration procedure for the model and illustrate numerically how it can reproduce observed behaviour of prices. This could then be used as a market simulator for short-term price prediction or for testing optimal execution strategies.


title: Delayed Stellar Mass Assembly in the Low Surface Brightness Dwarf Galaxy KDG215

abstract: We present HI spectral line and optical broadband images of the nearby low surface brightness dwarf galaxy KDG215. The HI images, acquired with the Karl G. Jansky Very Large Array (VLA), reveal a dispersion dominated ISM with only weak signatures of coherent rotation. The HI gas reaches a peak mass surface density of 6 M$_{\odot}$ pc$^{-2}$ at the location of the peak surface brightness in the optical and the UV. Although KDG215 is gas-rich, the H$\alpha$ non-detection implies a very low current massive star formation rate. In order to investigate the recent evolution of this system, we have derived the recent and lifetime star formation histories from archival Hubble Space Telescope images. The recent star formation history shows a peak star formation rate $\sim$1 Gyr ago, followed by a decreasing star formation rate to the present day quiescent state. The cumulative star formation history indicates that a significant fraction of the stellar mass assembly in KDG215 has occurred within the last 1.25 Gyr. KDG215 is one of only a few known galaxies which demonstrates such a delayed star formation history. While the ancient stellar population (predominantly red giants) is prominent, the look-back time by which 50% of the mass of all stars ever formed had been created is among the youngest of any known galaxy.


title: Quantized Magnetic Flux and the Magneto-halon Effect in a Critical Superconductor

abstract: Employing the standard worldline-vortex mapping, we conclude that at the critical temperature, superconductors demonstrate the magneto-halon effect with respect to the quantized net magnetic flux generated by a solenoid inserted into the system. The effect is a direct counterpart of the recently revealed halon effect in terms of the quantized particle charge of a static impurity in the two-dimensional U(1) quantum-critical environment. The flux-loop model (a.k.a. frozen lattice superconductor) with a quasi-solenoid perturbation proves to be the model of choice for qualitative and quantitative description of both the halon and the magneto-halon effects.


title: Improving Super-Resolution Methods via Incremental Residual Learning

abstract: Recently, Convolutional Neural Networks (CNNs) have shown promising performance in super-resolution (SR). However, these methods operate primarily on Low Resolution (LR) inputs for memory efficiency but this limits, as we demonstrate, their ability to (i) model high frequency information; and (ii) smoothly translate from LR to High Resolution (HR) space. To this end, we propose a novel Incremental Residual Learning (IRL) framework to address these mentioned issues. In IRL, first we select a typical SR pre-trained network as a master branch. Next we sequentially train and add residual branches to the main branch, where each residual branch is learned to model accumulated residuals of all previous branches. We plug state of the art methods in IRL framework and demonstrate consistent performance improvement on public benchmark datasets to set a new state of the art for SR at only approximately 20% increase in training time.


title: ALMA Thermal Observations of Europa

abstract: We present four daytime thermal images of Europa taken with the Atacama Large Millimeter Array. Together, these images comprise the first spatially resolved thermal dataset with complete coverage of Europa's surface. The resulting brightness temperatures correspond to a frequency of 233 GHz (1.3 mm) and a typical linear resolution of roughly 200 km. At this resolution, the images capture spatially localized thermal variations on the scale of geologic and compositional units. We use a global thermal model of Europa to simulate the ALMA observations in order to investigate the thermal structure visible in the data. Comparisons between the data and model images suggest that the large-scale daytime thermal structure on Europa largely results from bolometric albedo variations across the surface. Using bolometric albedos extrapolated from Voyager measurements, a homogenous model reproduces these patterns well, but localized discrepancies exist. These discrepancies can be largely explained by spatial inhomogeneity of the surface thermal properties. Thus, we use the four ALMA images to create maps of the surface thermal inertia and emissivity at our ALMA wavelength. From these maps, we identify a region of either particularly high thermal inertia or low emissivity near 90 degrees West and 23 degrees North, which appears anomalously cold in two of our images.


title: Clusters of cavity solitons bounded by conical radiation

abstract: We introduce a new class of self-sustained states, which may exist as single solitons or form multisoliton clusters, in driven passive cylindrical microresonators. Remarkably, such states are stabilized by the radiation they emit, which strongly breaks spatial symmetry and leads to the appearance of long polychromatic conical tails. The latter induce long-range soliton interactions that make possible the formation of clusters, which can be stable if their spatial arrangement is non-collinear with the soliton rotation direction in the microcavity. The clusters are intrinsically two-dimensional and, also, spatially rich. The mechanism behind the formation of the clusters is explained using soliton clustering theory. Our results bring fundamental understanding of a new class of multidimensional cavity solitons and may lead to the development of monolithic multi-soliton sources.


title: $C^{1,\alpha}$-subelliptic regularity on SU(3) and compact, semi-simple Lie groups

abstract: Let the vector fields $X_1, ... , X_{6}$ form an orthonormal basis of ${\mathcal H}$, the orthogonal complement of a Cartan subalgebra (of dimension $2$) in SU(3). We prove that weak solutions $u$ to the degenerate subelliptic $p$-Laplacian $$ \Delta_{\mathcal{H},{p}} u(x)=\sum_{i=1}^{6} X_i^{*}\left(|\nabla_{\hspace{-0.1cm} {\mathcal H}} u|^{p-2}X_{i}u \right) =0,$$ have H\"older continuous horizontal derivatives $\nabla_{\hspace{-0.1cm}{\mathcal H}} u=(X_1u, \ldots, X_{6}u)$ for $p\ge 2$. We also prove that a similar result holds for all compact connected semisimple Lie groups.


title: Carbon Vacancy Formation in Binary Transition Metal Carbides from Density Functional Theory

abstract: We have investigated the trends in the formation of carbon vacancies in binary transition metal (TM) carbides using density functional calculations for two common TM carbide crystal structures, namely the B1 (rock-salt) and WC structure types. The TM are taken from group IV (Ti, Zr, Hf), V (V, Nb, Ta) and VI (Cr, Mo, W) of the periodic table, as well as Sc from group III. For B1-structured TM carbides, the general trend is that it is easier for C vacancies to be formed as the number of valence electrons in the system increases. The exception is ScC where C vacancies are rather easy to form. It is also clear that the formation of C vacancies depends on the growth conditions: For TM-rich conditions, B1-structured carbides will always favour C vacancy formation. For C-rich conditions it is energetically favourable to form C vacancies in ScC, VC, NbC, CrC, MoC and WC. Experimentally large C vacancy concentrations are found in B1-structured TM carbides and our calculations are in line with these observations. In fact, TiC, ZrC, HfC and TaC are the only B1-structured TM carbides that do not spontaneously favour C vacancies and then only for C-rich conditions. C vacancies may still be present in these systems, however, due to the high temperatures used when growing TM carbides. In the case of WC-structured TM carbides, C vacancy formation is not energetically favourable irrespective of the growth conditions, which is the reason why this structure is only found close to a one-to-one metal-to-carbon ratio. ...


title: Observation of Brillouin optomechanical strong coupling with an 11 GHz mechanical mode

abstract: Achieving cavity-optomechanical strong coupling with high-frequency phonons provides a rich avenue for quantum technology development including quantum state-transfer, memory, and transduction, as well as enabling several fundamental studies of macroscopic phononic degrees-of-freedom. Reaching such coupling with GHz mechanical modes however has proved challenging, with a prominent hindrance being material- and surface-induced-optical absorption in many materials. Here, we circumvent these challenges and report the observation of optomechanical strong coupling to a high frequency (11 GHz) mechanical mode of a fused-silica whispering-gallery microresonator via the electrostrictive Brillouin interaction. Using an optical heterodyne detection scheme, the anti-Stokes light backscattered from the resonator is measured and normal-mode splitting and an avoided crossing are observed in the recorded spectra, providing unambiguous signatures of strong coupling. The optomechanical coupling rate reaches values as high as $G/2\pi = 39 \ \text{MHz}$ through the use of an auxiliary pump resonance, where the coupling dominates both the optical ($\kappa/2\pi = 3 \ \text{MHz}$) and the mechanical ($\gamma_\text{m}/2\pi = 21 \ \text{MHz}$) amplitude decay rates. Our findings provide a promising new approach for optical quantum control using light and sound.


title: Galois group action and Jordan decomposition of characters of finite reductive groups with connected center

abstract: Let $\mathbf{G}$ be a connected reductive group with connected center defined over $\mathbb{F}_q$, with Frobenius morphism F. Given an irreducible complex character $\chi$ of $\mathbf{G}^F$ with its Jordan decomposition, and a Galois automorphism $\sigma \in \mathrm{Gal}(\overline{\mathbb{Q}}/\mathbb{Q})$, we give the Jordan decomposition of the image ${^\sigma \chi}$ of $\chi$ under the action of $\sigma$ on its character values.


title: Satellite Positioning with Large Constellations

abstract: Modern global navigation satellite system receivers can access signals from several satellite constellations (including GPS, GLONASS, Galileo, BeiDou). Once these constellations are all fully operational, a typical receiver can expect to have on the order of 40-50 satellites in view. Motivated by that observation, this paper presents an asymptotic analysis of positioning algorithms in the large-constellation regime. We determine the exact asymptotic behavior for both pseudo-range and carrier-phase positioning. One interesting insight from our analysis is that the standard carrier-phase positioning approach based on resolving the carrier-phase integer ambiguities fails for large satellite constellations. Instead, we adopt a Bayesian approach, in which the ambiguities are treated as noise terms and not explicitly estimated.


title: Language Identification in Code-Mixed Data using Multichannel Neural Networks and Context Capture

abstract: An accurate language identification tool is an absolute necessity for building complex NLP systems to be used on code-mixed data. Lot of work has been recently done on the same, but there's still room for improvement. Inspired from the recent advancements in neural network architectures for computer vision tasks, we have implemented multichannel neural networks combining CNN and LSTM for word level language identification of code-mixed data. Combining this with a Bi-LSTM-CRF context capture module, accuracies of 93.28% and 93.32% is achieved on our two testing sets.


title: Phonon-induced multi-color correlations in hBN single-photon emitters

abstract: Color centers in hexagonal boron nitride have shown enormous promise as single-photon sources, but a clear understanding of electron-phonon interaction dynamics is critical to their development for quantum communications or quantum simulations. We demonstrate photon antibunching in the filtered auto- and cross-correlations $g^{(2)}_{lm}(\tau)$ between zero-, one- and two-phonon replicas of defect luminescence. Moreover, we combine autocorrelation measurements with a violation of the Cauchy-Schwarz inequality in the filtered cross-correlation measurements to distinguish a low quantum-efficiency defect from phonon replicas of a bright defect. With no background correction, we observe single photon purity of $g^{(2)}(0)=0.20$ in a phonon replica and cross-spectral correlations of $g^{(2)}_{lm}(0)=0.18$ between a phonon replica and the zero phonon line. These results illustrate a coherent interface between visible photons and mid-infrared phonons and provide a clear path toward control of photon-phonon entanglement in 2D materials.


title: Exploring a Unified Attention-Based Pooling Framework for Speaker Verification

abstract: The pooling layer is an essential component in the neural network based speaker verification. Most of the current networks in speaker verification use average pooling to derive the utterance-level speaker representations. Average pooling takes every frame as equally important, which is suboptimal since the speaker-discriminant power is different between speech segments. In this paper, we present a unified attention-based pooling framework and combine it with the multi-head attention. Experiments on the Fisher and NIST SRE 2010 dataset show that involving outputs from lower layers to compute the attention weights can outperform average pooling and achieve better results than vanilla attention method. The multi-head attention further improves the performance.


title: Efficient Propagation of Uncertainties in Manufacturing Supply Chains: Time Buckets, L-leap and Multilevel Monte Carlo

abstract: Uncertainty propagation of large scale discrete supply chains can be prohibitive when a large number of events occur during the simulated period and discrete event simulations (DES) are costly. We present a time bucket method to approximate and accelerate the DES of supply chains. Its stochastic version, which we call the L(logistic)-leap method, can be viewed as an extension of the leap methods, e.g., tau-leap, D-leap, developed in the chemical engineering community for the acceleration of stochastic DES of chemical reactions. The L-leap method instantaneously updates the system state vector at discrete time points and the production rates and policies of a supply chain are assumed to be stationary during each time bucket. We propose to use Multilevel Monte Carlo (MLMC) to efficiently propagate the uncertainties in a supply chain network, where the levels are naturally defined by the sizes of the time buckets of the simulations. We demonstrate the efficiency and accuracy of our methods using four numerical examples derived from a real world manufacturing material flow. In these examples, our multilevel L-leap approach can be faster than the standard Monte Carlo (MC) method by one or two orders of magnitudes without compromising the accuracy.


title: Heterogeneous and rate-dependent streptavidin-biotin unbinding revealed by high-speed force spectroscopy and atomistic simulations

abstract: Receptor-ligand interactions are essential for biological function and their binding strength is commonly explained in terms of static lock-and-key models based on molecular complementarity. However, detailed information of the full unbinding pathway is often lacking due, in part, to the static nature of atomic structures and ensemble averaging inherent to bulk biophysics approaches. Here we combine molecular dynamics and high-speed force spectroscopy on the streptavidin-biotin complex to determine the binding strength and unbinding pathways over the widest dynamic range. Experiment and simulation show excellent agreement at overlapping velocities and provided evidence of the unbinding mechanisms. During unbinding, biotin crosses multiple energy barriers and visits various intermediate states far from the binding pocket while streptavidin undergoes transient induced fits, all varying with loading rate. This multistate process slows down the transition to the unbound state and favors rebinding, thus explaining the long lifetime of the complex. We provide an atomistic, dynamic picture of the unbinding process, replacing a simple two-state picture with one that involves many routes to the lock and rate-dependent induced-fit motions for intermediates, which might be relevant for other receptor-ligand bonds.


title: Multi-Higgs doublet models: the Higgs-fermion couplings and their sum rules

abstract: This is the second of a series of papers that explores the physical parameterization, sum rules and unitarity bounds arising from a non-minimal scalar sector of the Standard Model (SM) that consists of N Higgs doublets. In this paper, we focus on the structure and implication of the Yukawa interactions that couple the N scalar doublets to the SM fermions. We employ the charged Higgs basis, which is defined as the basis of scalar fields such that the neutral scalar field vacuum expectation value resides entirely in one of the N scalar doublet fields, and the charged components of the remaining N-1 scalar doublet fields are the physical (mass-eigenstate) charged Higgs fields. Based on the structure of the Yukawa Lagrangian of the model (and as a consequence of tree-level unitarity), one may deduce numerous sum rules, several of which have not appeared previously in the literature. These sum rules can be used to uncover intimate relations between the structure of the Higgs-fermion couplings and the scalar/gauge couplings. In particular, we show that the approximate alignment limit, in which the W+W- and ZZ couplings to the observed Higgs boson are approximately SM-like, imposes significant constraints on the Higgs-fermion couplings.


title: Strong Jump Inversion

abstract: We say that a structure $\mathcal{A}$ admits \emph{strong jump inversion} provided that for every oracle $X$, if $X'$ computes $D(\mathcal{C})'$ for some $\mathcal{C}\cong\mathcal{A}$, then $X$ computes $D(\mathcal{B})$ for some $\mathcal{B}\cong\mathcal{A}$. Jockusch and Soare \cite{JS} showed that there are low linear orderings without computable copies, but Downey and Jockusch \cite{DJ} showed that every Boolean algebra admits strong jump inversion. More recently, D.\ Marker and R.\ Miller \cite{MM} have shown that all countable models of $DCF_0$ (the theory of differentially closed fields of characteristic $0$) admit strong jump inversion. We establish a general result with sufficient conditions for a structure $\mathcal{A}$ to admit strong jump inversion. Our conditions involve an enumeration of $B_1$-types, where these are made up of formulas that are Boolean combinations of existential formulas. Our general result applies to some familiar kinds of structures, including some classes of linear orderings and trees. We do not get the result of Downey and Jockusch for arbitrary Boolean algebras, but we do get a result for Boolean algebras with no $1$-atom, with some extra information on the complexity of the isomorphism. Our general result gives the result of Marker and Miller. In order to apply our general result, we produce a computable enumeration of the types realized in models of $DCF_0$. This also yields the fact that the saturated model of $DCF_0$ has a decidable copy.


title: Generalized micropolar continualization of 1D beam lattices

abstract: The enhanced continualization approach proposed in this paper is aimed to overcome some drawbacks observed in the homogenization of beam lattices. To this end an enhanced homogenization technique is proposed and formulated to obtain consistent micropolar continuum models of the beam lattices and able to simulate with good approximation the boundary layer effects and the Floquet-Bloch spectrum of the Lagrangian model. The continualization technique here proposed is based on a transformation of the difference equation of motion of the discrete system via a proper down-scaling law into a pseudo-differential problem; a further McLaurin approximation is applied to obtain a higher order differential problem. The formulation is carried out for simple one-dimensional beam lattices that are, nevertheless, characterized by a rather wide variety of static and dynamic behaviors: the rod lattice, the beam lattice with node rotations and a 1D beam lattice model with generalized displacements. Higher order models may be obtained which are characterized by differential problems involving non-local inertia terms together with spatial high gradient terms. Moreover, the homogenized models obtained by the proposed enhanced continualization technique turn out to be energetically consistent and provide a good simulation of both the static response and of the acoustic spectrum of the original discrete models. The proposed homogenization procedure is first presented for the simple case of monoatomic axial chains. The beam lattice with node rotations and displacement prevented exhibits, in the static regime, decaying oscillations of the nodal rotation in the boundary layer which is well simulated by the homogenized model obtained by the proposed approach.


title: Medium polarization effects in 3SD1 spin-triplet pairing

abstract: Stimulated by the still puzzling competition between spin-singlet and spin-triplet pairing in nuclei, the 3SD1 neutron-proton pairing is investigated in the framework of BCS theory of nuclear matter. The medium polarization effects are included in the single particle spectrum and also in the pairing interaction starting from the G-matrix, calculated in the Brueckner-Hartree-Fock approximation. The vertex corrections due to spin and isospin collective excitations of the medium are determined from the Bethe-Salpeter equation in the RPA limit, taking into account the tensor correlations. It is found that the self-energy corrections confine the superfluid state to very low-density, while remarkably quenching the magnitude of the energy gap, while the induced interaction has an attractive effect. The interplay between spin-singlet and spin-triplet pairing is discussed in nuclear matter as well as in finite nuclei.


title: Statistical inference and feasibility determination: a nonasymptotic approach

abstract: We develop non-asymptotically justified methods for hypothesis testing about the $p-$dimensional coefficients $\theta^{*}$ in (possibly nonlinear) regression models. Given a function $h:\,\mathbb{R}^{p}\mapsto\mathbb{R}^{m}$, we consider the null hypothesis $H_{0}:\,h(\theta^{*})\in\Omega$ against the alternative hypothesis $H_{1}:\,h(\theta^{*})\notin\Omega$, where $\Omega$ is a nonempty closed subset of $\mathbb{R}^{m}$ and $h$ can be nonlinear in $\theta^{*}$. Our (nonasymptotic) control on the Type I and Type II errors holds for fixed $n$ and does not rely on well-behaved estimation error or prediction error; in particular, when the number of restrictions in $H_{0}$ is large relative to $p-n$, we show it is possible to bypass the sparsity assumption on $\theta^{*}$ (for both Type I and Type II error control), regularization on the estimates of $\theta^{*}$, and other inherent challenges in an inverse problem. We also demonstrate an interesting link between our framework and Farkas' lemma (in math programming) under uncertainty, which points to some potential applications of our method outside traditional hypothesis testing.


title: Fixation Probabilities in Weakly Compressible Fluid Flows

abstract: Competition between biological species in marine environments is affected by the motion of the surrounding fluid. An effective 2D compressibility can arise, for example, from the convergence and divergence of water masses at the depth at which passively traveling photosynthetic organisms are restricted to live. In this report, we seek to quantitatively study genetics under flow. To this end, we couple an off-lattice agent-based simulation of two populations in 1D to a weakly compressible velocity field--first a sine wave and then a shell model of turbulence. We find for both cases that even in a regime where the overall population structure is approximately unaltered, the flow can significantly diminish the effect of a selective advantage on fixation probabilities. We understand this effect in terms of the enhanced survival of organisms born at sources in the flow and the influence of Fisher genetic waves.


title: Integration in power-bounded $T$-convex valued fields

abstract: This is the second installment of a series of papers aimed at developing a theory of Hrushovski-Kazhdan style motivic integration for certain types of nonarchimedean $o$-minimal fields, namely power-bounded $T$-convex valued fields, and closely related structures. The main result in the first installment is a canonical isomorphism between the Grothendieck rings of certain categories of definable sets, which is understood as a universal additive invariant or a generalized Euler characteristic because the categories do not carry volume forms. Here we introduce two types of volume forms into each of the relevant categories, one takes values in the value group and the other in the finer RV-sort. The resulting isomorphisms respect Jacobian transformations --- that is, the change of variables formula holds --- and hence are regarded as motivic integrals. As in the classical theory of integration, one is often led to consider locally constant functions with bounded support in various situations. For the space of such functions, the construction may be fine-tuned so as to become more amenable to applications. The modifications are nevertheless substantial and constitute the bulk of the technical work.


title: Global classical solutions to the continuous coagulation equation with collisional breakage

abstract: Existence and uniqueness of mass-conserving classical solutions to the continuous coagulation equation with collisional breakage are investigated for an unbounded class of collision kernels and a particular case of the distribution function. The distribution function may have a possibility to attain singularity at the origin. The proof of the existence result relies on the compactness method. Moreover, a uniqueness result is shown. In addition, it is observed that the uniqueness solution is mass-conserving.


title: Unstable entropy of partially hyperbolic diffeomorphisms along non-compact subsets

abstract: Given a partially hyperbolic diffeomorphism $f:M \rightarrow M$ defined on a compact Riemannian manifold $M$, in this paper we define the concept of unstable topological entropy of $f$ on a set $Y \subset M$ not necessarily compact and we extend a theorem of R. Bowen proving that, for an ergodic $f$-invariant measure $\mu$, the unstable measure theoretical entropy of $f$ is upper bounded by the unstable topological entropy of $f$ on any set of full $\mu$-measure. We also define a notion of unstable topological entropy of $f$ using a Hausdorff dimension like characterization.


title: A finitely presented ${E}_{\infty}$-prop II: cellular context

abstract: We construct, using finitely many generating cell and relations, props in the category of CW-complexes with the property that their associated operads are models for the $E_\infty$-operad. We use one of these to construct a cellular $E_\infty$-bialgebra structure on the interval and derive from it a natural cellular $E_\infty$-coalgebra structure on the geometric realization of a simplicial set which, passing to cellular chains, recovers up to signs the Barratt-Eccles and Surjection coalgebra structures introduced by Berger-Fresse and McClure-Smith. We use another prop, a quotient of the first, to relate our constructions to earlier work of Kaufmann and prove a conjecture of his. This is the second of two papers in a series, the first investigates analogue constructions in the category of chain complexes.


title: Zeros of polynomials with four-term recurrence and linear coefficients

abstract: This paper investigates the zero distribution of a sequence of polynomials $\left\{ P_{m}(z)\right\} _{m=0}^{\infty}$ generated by the reciprocal of $1+ct+B(z)t^{2}+A(z)t^{3}$ where $c\in\mathbb{R}$ and $A(z)$, $B(z)$ are real linear polynomials. We study necessary and sufficient conditions for the reality of the zeros of $P_{m}(z)$. Under these conditions, we find an explicit interval containing these zeros, whose union forms a dense subset of this interval.


title: Unifying fast scrambling, thermalization and entanglement through the measurement of FOTOCs in the Dicke model

abstract: Scrambling of quantum information is the process by which information initially stored in the local degrees of freedom of a quantum many-body system spreads over its many-body degrees of freedom, becoming inaccessible to local probes and thus apparently lost. Scrambling and entanglement are key concepts reconciling seemingly unrelated behaviors including thermalization of isolated quantum systems and information loss in black holes, and have revolutionized our understanding of non-equilibrium phenomena. Here, we demonstrate that a family of fidelity out-of-time-order correlators (FOTOCs), recently measured in a trapped-ion quantum simulator via time reversal of the many-body dynamics followed by a fidelity measurement, can serve as a unifying diagnostic tool that elucidates the intrinsic connection between fast scrambling, volume law entanglement, ergodicity, quantum chaos, and the associated butterfly effect in the semiclassical dynamics of the system. We demonstrate the utility of FOTOCs by computing them in the Dicke model, an iconic model in quantum optics, recently implemented in atomic and trapped-ion setups. This model describes the coupling of a large spin to an oscillator and features rich behaviors, including a quantum phase transition and chaos. Here, we show that FOTOCs provide a direct measure of the spin-phonon Renyi entropy and quantum thermalization. Moreover, we connect the FOTOCs to the variance of simple operators, allowing us to observe fast scrambling in the parameter regime where the system's classical trajectories are chaotic, and to explicitly relate the quantum and classical Lyapunov exponents in a truly quantum many-body system. Our results open a path for the experimental use of FOTOCs to quantify fast scrambling, determine bounds on quantum information processing and to identify possible candidates of black hole analogs in controllable quantum systems.


title: Search for lepton-flavour-violating decays of Higgs-like bosons

abstract: A search is presented for a Higgs-like boson with mass in the range 45 to 195 GeV/$c^2$ decaying into a muon and a tau lepton. The dataset consists of proton-proton interactions at a centre-of-mass energy of 8 TeV, collected by the LHCb experiment, corresponding to an integrated luminosity of 2 fb$^{-1}$. The tau leptons are reconstructed in both leptonic and hadronic decay channels. An upper limit on the production cross-section multiplied by the branching fraction at 95% confidence level is set and ranges from 22 pb for a boson mass of 45 GeV/$c^2$ to 4 pb for a mass of 195 GeV/$c^2$.


title: Particle Acceleration and Heating in a Turbulent Solar Corona

abstract: Turbulence, magnetic reconnection, and shocks can be present in explosively unstable plasmas, forming a new electromagnetic environment, which we call here turbulent reconnection, and where spontaneous formation of current sheets takes place. We will show that the heating and the acceleration of particles is the result of the synergy of stochastic (second order Fermi) and systematic (first order Fermi) acceleration inside fully developed turbulence. The solar atmosphere is magnetically coupled to a turbulent driver (the convection zone), therefore the appearance of turbulent reconnection in the solar atmosphere is externally driven. Turbulent reconnection, once it is established in the solar corona, drives the coronal heating and particle acceleration.


title: Nontangential estimates on layer potentials and the Neumann problem for higher order elliptic equations

abstract: We solve the Neumann problem, with nontangential estimates, for higher order divergence form elliptic operators with variable $t$-independent coefficients. Our results are accompanied by nontangential estimates on higher order layer potentials.


title: Geometrical effects on mobility

abstract: In this paper we analyze the effect of randomly deleting streets of a synthetic city on the statistics of displacements. Our city is constituted initially by a set of streets that form a regular tessellation of the euclidean plane. Therefore we will have three types of cities, formed by squares, triangles or hexagons. We studied the complementary cumulative distribution function for displacements (CCDF). For the whole set of streets the CCDF is a stretched exponential, and as streets are deleted this function becomes a linear function and then two clear different exponentials. This behavior is qualitatively the same for all the tessellations. Most of this functions has been reported in the literature when studying the displacements of individuals based on cell data trajectories and GPS information. However, in the light of this work, the appearance of different functions for displacements CCDF can be attributed to the connectivity of the underlying street network. It is remarkably that for some proportion of streets we got a linear function for such function, and as far as we know this behavior has not been reported nor considered. Therefore, it is advisable to analyze experimental in the light of connectivity of the street network to make correlations with the present work.


title: Low-Complexity Reconfigurable MIMO for Millimeter Wave Communications

abstract: The performance of millimeter wave (mmWave) multiple-input multiple-output (MIMO) systems is limited by the sparse nature of propagation channels and the restricted number of radio frequency (RF) chains at transceivers. The introduction of reconfigurable antennas offers an additional degree of freedom on designing mmWave MIMO systems. This paper provides a theoretical framework for studying the mmWave MIMO with reconfigurable antennas. Based on the virtual channel model, we present an architecture of reconfigurable mmWave MIMO with beamspace hybrid analog-digital beamformers and reconfigurable antennas at both the transmitter and the receiver. We show that employing reconfigurable antennas can provide throughput gain for the mmWave MIMO. We derive the expression for the average throughput gain of using reconfigurable antennas in the system, and further derive the expression for the outage throughput gain for the scenarios where the channels are (quasi) static. Moreover, we propose a low-complexity algorithm for reconfiguration state selection and beam selection. Our numerical results verify the derived expressions for the throughput gains and demonstrate the near-optimal throughput performance of the proposed low-complexity algorithm.


title: Multinomial Models with Linear Inequality Constraints: Overview and Improvements of Computational Methods for Bayesian Inference

abstract: Many psychological theories can be operationalized as linear inequality constraints on the parameters of multinomial distributions (e.g., discrete choice analysis). These constraints can be described in two equivalent ways: Either as the solution set to a system of linear inequalities or as the convex hull of a set of extremal points (vertices). For both representations, we describe a general Gibbs sampler for drawing posterior samples in order to carry out Bayesian analyses. We also summarize alternative sampling methods for estimating Bayes factors for these model representations using the encompassing Bayes factor method. We introduce the R package multinomineq, which provides an easily-accessible interface to a computationally efficient implementation of these techniques.


title: Orbital- and $k_z$-selective hybridisation of Se 4p and Ti 3d states in the charge density wave phase of TiSe$_2$

abstract: We revisit the enduring problem of the $2\times{}2\times{}2$ charge density wave (CDW) order in TiSe$_2$, utilising photon energy-dependent angle-resolved photoemission spectroscopy to probe the full three-dimensional high- and low-temperature electronic structure. Our measurements demonstrate how a mismatch of dimensionality between the 3D conduction bands and the quasi-2D valence bands in this system leads to a hybridisation that is strongly $k_z$-dependent. While such a momentum-selective coupling can provide the energy gain required to form the CDW, we show how additional "passenger" states remain, which couple only weakly to the CDW and thus dominate the low-energy physics in the ordered phase of TiSe$_2$.


title: Mathematical modelling indicates that lower activity of the haemostatic system in neonates is primarily due to lower prothrombin concentration

abstract: Haemostasis is governed by a highly complex system of interacting proteins. Due to the central role of thrombin, thrombin generation and specifically the thrombin generation curve (TGC) is commonly used as an indicator of haemostatic activity. Functional characteristics of the haemostatic system in neonates and children are significantly different compared with adults; at the same time plasma levels of haemostatic proteins vary considerably with age. However, relating one to the other has been difficult, both due to significant inter-individual differences for individuals of similar age and the complexity of the biochemical reactions underlying haemostasis. Mathematical modelling has been very successful at representing the biochemistry of blood clotting. In this study we address the challenge of large inter-individual variability by parameterising the Hockin-Mann model with data from individual patients, across different age groups from neonates to adults. Calculating TGCs for each patient of a specific age group provides us with insight into the variability of haemostatic activity across that age group. From our model we observe that two commonly used metrics for haemostatic activity are significantly lower in neonates than in older patients. Because both metrics are strongly determined by prothrombin and prothrombin levels are considerably lower in neonates we conclude that decreased haemostatic activity in neonates is due to lower prothrombin availability.


title: Polarization rotation in Bi$_{\mathbf{4}}$Ti$_{\mathbf{3}}$O$_{\mathbf{12}}$ by isovalent doping at the fluorite sublattice

abstract: Bismuth titanate, Bi$_4$Ti$_3$O$_{12}$ (BiT), is a complex layered ferroelectric material that is composed of three perovskite-like units and one fluorite-like unit stacked alternatively along the $c$-direction. The ground state crystal structure is monoclinic with the spontaneous polarization (~50 $\mu$C/cm$^{2}$) along the in-plane $b$-direction. BiT typically grows along the $c$-direction in thin film form and having the polarization vector aligned with the growth orientation can be beneficial for several potential device applications. It is well known that judicious doping of ferroelectrics is an effective method in adjusting the magnitude and the orientation of the spontaneous polarization. Here, we show using first-principles density functional theory and a detailed phonon analysis that Bi atoms in the fluorite-like layers have significantly more impact on the magnitude and orientation of the spontaneous polarization vector as compared to the perovskite-like layer. The low energy hard phonon modes are characterized by fluorite-like layers experiencing transverse displacements and large changes in Born effective charges on Bi atoms. Thus, the breaking of symmetry caused by doping of Bi sites within the fluorite-like layer leads to the formation of uncancelled permanent dipole moments along the $c$-direction. This provides an opportunity for doping the Bi site in the fluorite-like layer. Isovalent dopants P, As, and Sb were studied. P is found to be most effective in the reorientation of the spontaneous polarization. It leads to a three-fold enhancement of the $c$-component of polarization and to a commensurate rotation of the spontaneous polarization vector by 36.2$^{\circ}$ towards the $c$-direction.


title: A compact $G_2$-calibrated manifold with first Betti number $b_1=1$

abstract: We construct a compact formal 7-manifold with a closed $G_2$-structure and with first Betti number $b_1=1$, which does not admit any torsion-free $G_2$-structure, that is, it does not admit any $G_2$-structure such that the holonomy group of the associated metric is a subgroup of $G_2$. We also construct associative calibrated (hence volume-minimizing) 3-tori with respect to this closed $G_2$-structure and, for each of those 3-tori, we show a 3-dimensional family of non-trivial associative deformations. We also construct a fibration of our 7-manifold over $S^2\times S^1$ with generic fiber a (non-calibrated) coassociative 4-torus and some singular fibers.


title: Data Analysis and Phenomenological Cosmology

abstract: In the era of precision cosmology, even percentage level effects are significant on cosmological observables. The recent tension between the local and global values of $H_0$ is much more significant than this, and any possible solution might rely on us going beyond the standard $\Lambda$CDM cosmological model. For much smaller, yet potentially significant effects, spatial curvature from averaging and cosmological backreaction on observational predictions could play a role. This is especially true with the higher precision of new observational data and improved statistical techniques. In this paper, we discuss the observational viability of a class of physically motivated cosmologies which can be parametrized by a phenomenological two-scale backreaction model with decoupled spatial curvature parameters and two Hubble scales. Using the latest JLA Supernovae data together with some of the latest BAO data, we perform a Bayesian model selection analysis and find that the phenomenological models are not favoured over the standard $\Lambda$CDM cosmological model. Although there is still a preference for non-zero and unequal dynamic and geometric spatial curvatures, there is little evidence for differing Hubble scales within these phenomenological template models.


title: Control Distance and Energy Scaling of Complex Networks

abstract: It has recently been shown that the average energy required to control a subset of nodes in a complex network scales exponentially with the cardinality of the subset. While the mean scales exponentially, the variance of the control energy over different subsets of nodes is large and has as of yet not been explained. Here, we provide an explanation of the large variance as a result of both the length of the path that connects control inputs to the target nodes and the redundancy of paths of shortest length. Our first result provides an exact upper bound of the control energy as a function of path length between driver node and target node along an infinite path graph. We also show that the energy estimation is still very accurate even when finite size effects are taken into account. Our second result refines the upper bound that takes into account not only the length of the path, but also the redundancy of paths. We show that it improves the upper bound approximation by an order of magnitude or more. Finally, we lay out the foundations for a more accurate estimation of the control energy for the multi-target and multi-driver problem.


title: Lectures on Polyfolds and Symplectic Field Theory

abstract: This is a lecture note prepared for the SFT 9 workshop in Augsburg, Germany. The text describes a polyfold approach to the construction of symplectic field theory and focuses on the perturbation and transversality theory.


title: Experimental quantum tomography assisted by multiply symmetric states in higher dimensions

abstract: High-dimensional quantum information processing has become a mature field of research with several different approaches being adopted for the encoding of $D$-dimensional quantum systems. Such progress has fueled the search of reliable quantum tomographic methods aiming for the characterization of these systems, being most of these methods specifically designed for a given scenario. Here, we report on a new tomographic method based on multiply symmetric states and on experimental investigations to study its performance in higher dimensions. Unlike other methods, it is guaranteed to exist in any dimension and provides a significant reduction in the number of measurement outcomes when compared to standard quantum tomography. Furthermore, in the case of odd dimensions, the method requires the least possible number of measurement outcomes. In our experiment we adopt the technique where high-dimensional quantum states are encoded using the linear transverse momentum of single photons and are controlled by spatial light modulators. Our results show that fidelities of $0.984\pm0.009$ with ensemble sizes of only $1.5\times10^5$ photons in dimension $D=15$ can be obtained in typical laboratory conditions, thus showing its practicability in higher dimensions.


title: Revealing evolutionary constraints on proteins through sequence analysis

abstract: Statistical analysis of alignments of large numbers of protein sequences has revealed "sectors" of collectively coevolving amino acids in several protein families. Here, we show that selection acting on any functional property of a protein, represented by an additive trait, can give rise to such a sector. As an illustration of a selected trait, we consider the elastic energy of an important conformational change within an elastic network model, and we show that selection acting on this energy leads to correlations among residues. For this concrete example and more generally, we demonstrate that the main signature of functional sectors lies in the small-eigenvalue modes of the covariance matrix of the selected sequences. However, secondary signatures of these functional sectors also exist in the extensively-studied large-eigenvalue modes. Our simple, general model leads us to propose a principled method to identify functional sectors, along with the magnitudes of mutational effects, from sequence data. We further demonstrate the robustness of these functional sectors to various forms of selection, and the robustness of our approach to the identification of multiple selected traits.


title: Faster PET Reconstruction with Non-Smooth Priors by Randomization and Preconditioning

abstract: Uncompressed clinical data from modern positron emission tomography (PET) scanners are very large, exceeding 350 million data points (projection bins). The last decades have seen tremendous advancements in mathematical imaging tools many of which lead to non-smooth (i.e. non-differentiable) optimization problems which are much harder to solve than smooth optimization problems. Most of these tools have not been translated to clinical PET data, as the state-of-the-art algorithms for non-smooth problems do not scale well to large data. In this work, inspired by big data machine learning applications, we use advanced randomized optimization algorithms to solve the PET reconstruction problem for a very large class of non-smooth priors which includes for example total variation, total generalized variation, directional total variation and various different physical constraints. The proposed algorithm randomly uses subsets of the data and only updates the variables associated with these. While this idea often leads to divergent algorithms, we show that the proposed algorithm does indeed converge for any proper subset selection. Numerically, we show on real PET data (FDG and florbetapir) from a Siemens Biograph mMR that about ten projections and backprojections are sufficient to solve the MAP optimisation problem related to many popular non-smooth priors; thus showing that the proposed algorithm is fast enough to bring these models into routine clinical practice.


title: MobilityMirror: Bias-Adjusted Transportation Datasets

abstract: We describe customized synthetic datasets for publishing mobility data. Private companies are providing new transportation modalities, and their data is of high value for integrative transportation research, policy enforcement, and public accountability. However, these companies are disincentivized from sharing data not only to protect the privacy of individuals (drivers and/or passengers), but also to protect their own competitive advantage. Moreover, demographic biases arising from how the services are delivered may be amplified if released data is used in other contexts. We describe a model and algorithm for releasing origin-destination histograms that removes selected biases in the data using causality-based methods. We compute the origin-destination histogram of the original dataset then adjust the counts to remove undesirable causal relationships that can lead to discrimination or violate contractual obligations with data owners. We evaluate the utility of the algorithm on real data from a dockless bike share program in Seattle and taxi data in New York, and show that these adjusted transportation datasets can retain utility while removing bias in the underlying data.


title: Validity and Failure of the Boltzmann Weight

abstract: The dynamics and thermostatistics of a classical inertial XY model, characterized by long-range interactions, are investigated on $d$-dimensional lattices ($d=1,2,$ and 3), through molecular dynamics. The interactions between rotators decay with the distance $r_{ij}$ like~$1/r_{ij}^{\alpha}$ ($\alpha \geq 0$), where $\alpha\to\infty$ and $\alpha=0$ respectively correspond to the nearest-neighbor and infinite-range interactions. We verify that the momenta probability distributions are Maxwellians in the short-range regime, whereas $q$-Gaussians emerge in the long-range regime. Moreover, in this latter regime, the individual energy probability distributions are characterized by long tails, corresponding to $q$-exponential functions. The present investigation strongly indicates that, in the long-range regime, central properties fall out of the scope of Boltzmann-Gibbs statistical mechanics, depending on $d$ and $\alpha$ through the ratio $\alpha/d$.


title: Coulomb Log for Conductivity of Dense Plasmas

abstract: The Coulomb log (log {\Lambda}) approximation is widely used to approximate electron transport coefficients in dense plasmas. It is a classical approximation to the momentum transport cross section. The accuracy of this approximation for electrical conductivity in dense plasmas is assessed by comparing to fully quantum mechanical calculations for realistic scattering potentials. It is found that the classical approximation is accurate to +/-10% when log {\Lambda} > 3, irrespective of the plasma species. The thermodynamic regime (density and temperature) for which log {\Lambda} > 3 corresponds to does, however, strongly depend on the material. For increasing Z, log {\Lambda} is greater than 3 for increasingly high temperatures and lower densities.


title: Assessing sources of uncertainty in length-based estimates of body growth in populations of fishes and macroinvertebrates with bootstrapped ELEFAN

abstract: The determination of rates of body growth is the first step in many aquatic population studies and fisheries stock assessments. ELEFAN (Electronic LEngth Frequency ANalysis) is a widely used method to fit a growth curve to length-frequency distribution (LFD) data. However, up to now, it was not possible to assess its accuracy or the uncertainty inherent of this method, or to obtain confidence intervals for growth parameters within an unconstrained search space. In this study, experiments were conducted to assess the precision and accuracy of bootstrapped and single-fit ELEFAN-based curve fitting methods, using synthetic LFDs with known input parameters and a real data set of Abra alba shell lengths. The comparison of several types of bootstrap experiments and their outputs (95% confidence intervals and confidence contour plots) provided a first glimpse into the accuracy of modern ELEFAN-based fit methods. The main components of uncertainty (precision and reproducibility of fit algorithms, seed effects, sample size and matrix information content) could be assessed from partial bootstraps. Uncertainty was mainly determined by LFD matrix size, total number of non-zero bins and the sampling of large-sized individuals. A new pseudo-Rsquared index for the goodness-of-fit of VBGF models to LFD data is proposed. For a large, perfect synthetic data set, pseudo-RsquaredPhi was very high (88 to 100%), indicating an excellent fit of the VBGF model. The small Abra alba data set showed a low pseudo-RsquaredPhi, from to 54% to 68%, indicating the need for more samples and a larger LFD data matrix. New, robust, bootstrap-based methods for curve fitting are presented and discussed. This study demonstrates a promising new path for length-based analyses of growth and mortality in natural populations, which are the basis for a new suite of methods that are included in the new fishboot package.


title: Polar Convolution

abstract: The Moreau envelope is one of the key convexity-preserving functional operations in convex analysis, and it is central to the development and analysis of many approaches for convex optimization. This paper develops the theory for an analogous convolution operation, called the polar envelope, specialized to gauge functions. Many important properties of the Moreau envelope and the proximal map are mirrored by the polar envelope and its corresponding proximal map. These properties include smoothness of the envelope function, uniqueness and continuity of the proximal map, which play important roles in duality and in the construction of algorithms for gauge optimization. A suite of tools with which to build algorithms for this family of optimization problems is thus established.


title: The Planar Modular Partition Monoid

abstract: The primary contribution of this thesis is to introduce and examine the planar modular partition monoid for parameters $m, k \in \mathbb{Z}_{>0}$, which has simultaneously and independently generated interest from other researchers as outlined within. Our collective understanding of related monoids, in particular the Jones, Brauer and partition monoids, along with the algebras they generate, has heavily influenced the direction of research by a significant number of mathematicians and physicists. Examples include Schur-Weyl type dualities in representation theory along with Potts, ice-type and Andrew-Baxter-Forrester models from statistical mechanics, giving strong motivation for the planar modular partition monoid to be examined. The original results contained within this thesis relating to the planar modular partition monoid are: the establishment of generators; recurrence relations for the cardinality of the monoid; recurrence relations for the cardinality of Green's $\mathcal{R}$, $\mathcal{L}$ and $\mathcal{D}$ relations; and a conjecture on relations that appear to present the planar modular partition monoid when $m=2$. For diagram semigroups that are closed under vertical reflections, characterisations of Green's $\mathcal{R}$, $\mathcal{L}$ and $\mathcal{H}$ relations have previously been established using the upper and lower patterns of bipartitions. We give a characterisation of Green's $\mathcal{D}$ relation with a similar flavour for diagram semigroups that are closed under vertical reflections. $<paragraph\ missing\ due\ to\ arXiv\ character\ limit\ on\ abstracts\ -\ see\ thesis\ file>$ A further contribution of this thesis is a reasonably comprehensive exposition of the fundamentals of diagram semigroups, which have traditionally been approached from the representation theory side and have since blossomed into a thriving area of research in their own right.


title: Orientable 4-dimensional Poincar\'e Complexes have Reducible Spivak Fibrations

abstract: We show that the Spivak normal fibration of an orientable 4-dimensional Poincar\'e complex has a vector bundle reduction.


title: Five-body choreography on the algebraic lemniscate is a potential motion

abstract: In a remarkable paper of 2003 by Fujiwara et al. \cite{Fujiwara2003}, a figure-eight three-body choreography on the algebraic lemniscate by Bernoulli was discovered. Such a choreography was found to be driven by the action of a pairwise potential $ V(r_{ij}) $, depending only on the mutual relative distances $r_{ij}, i,j=1,2,3$. In the present Letter we show that two different choreographies of five bodies on the same algebraic lemniscate exist and correspond to solutions of ten coupled Newton equations of motion with a pairwise interaction potential. For each choreography the explicit form of the potential is found and ten constants of motion are presented, thus, it is superintegrable.


title: Calculus on Dual Real Numbers

abstract: We present the basic theory of calculus on dual real numbers, and prove the counterpart of the ordinary fundamental theorem of calculus in the context of dual real numbers.


title: Meager Sets, Games and Singular Cardinals

abstract: We show that a statement concerning the existence of winning strategies of limited memory in an infinite two-person topological game is equivalent to a weak version of the Singular Cardinals Hypothesis.


title: Implications of loop-top origin for microwave, hard X-ray, and low-energy gamma-ray emissions from behind the limb flares

abstract: The Fermi gamma-ray Space Telescope (Fermi) has detected hard X-ray (HXR) and gamma-ray photons from three flares, which according to \stereo occurred in active regions behind the limb of the Sun as delineated by near Earth instruments. For two of these flares \r has provided HXR images with sources located just above the limb, presumably from the loop top (LT) region of a relatively large loop. Fermi-Gamma-ray Burst Monitor has detected HXRs and gamma-rays, and RSTN has detected microwaves emissions with similar light curves. This paper presents a quantitative analysis of these multi-wavelength observations assuming that HXRs and microwaves are produced by electrons accelerated at the LT source, with emphasize on the importance of the proper treatment of escape of the particles from the acceleration-source region and the trans-relativistic nature of the analysis. The observed spectra are used to determine the magnetic field and relativistic electron spectra. It is found that a simple power-law in momentum (with cut off above a few 100 MeV) agrees with all observations, but in energy space a broken power law spectrum (steepening at rest mass energy) may be required. It is also shown that the production of the $>100$ MeV photons detected by The Fermi-Large Area Telescope at the LT source would require more energy compared to photospheric emission. These energies are smaller than that required for electrons, so that the possibility that all the emissions originate in the LT cannot be ruled out on energetic grounds. However, the differences in the light curves and emission centroids of HXRs and $>100$ MeV gamma-rays favour a different source for the latter.


title: Persistence and time periodic positive solutions of doubly nonlocal Fisher-KPP equations in time periodic and space heterogeneous media

abstract: In this paper, we investigate the asymptotic dynamics of Fisher-KPP equations with nonlocal dispersal operator and nonlocal reaction term in time periodic and space heterogeneous media. We first show the global existence and boundedness of nonnegative solutions, and then obtain some sufficient conditions ensuring the uniform persistence. In particular, we study the existence, uniqueness and global stability of positive time periodic solutions under several different conditions.


title: Signature-based Non-orthogonal Multiple Access (S-NOMA) for Massive Machine-Type Communications in 5G

abstract: The problem of providing massive connectivity in Internet-of-Things (IoT) with a limited number of available resources motivates the non-orthogonal multiple access (NOMA) solutions. In this article, we provide a comprehensive review of the signature-based NOMA (S-NOMA) schemes as potential candidates for IoT. The signature in S-NOMA represents the way the data stream of an active device is spread over available resources in a non-orthogonal manner. It can be designed based on device-specific codebook structures, delay patterns, spreading sequences, interleaving patterns, and scrambling sequences. Additionally, we present the detection algorithms employed to decode each device's data from non-orthogonally superimposed signals at the receiver. The bit error rate of different S-NOMA schemes is simulated in impulsive noise environments, which can be important in machine-type communications. Simulation results show that the performance of the S-NOMA schemes degrades under such conditions. Finally, research challenges in S-NOMA oriented IoT are presented.


title: The Spatially Resolved Dust-to-Metals Ratio in M101

abstract: The dust-to-metals ratio describes the fraction of the heavy elements contained in dust grains, and its variation provides key insights into the life cycle of dust. We measure the dust-to-metals ratio in M101, a nearby galaxy with a radial metallicity (Z) gradient spanning $\sim$1 dex. We fit the dust spectral energy distribution from 100 to 500 $\mu m$ with five variants of the modified blackbody dust emission model in which we vary the temperature distribution and how emissivity depends on wavelength. Among them, the model with a single temperature blackbody modified by a broken power-law emissivity gives the statistically best fit and physically most plausible results. Using these results, we show that the dust-to-gas ratio is proportional to $\rm Z^{1.7}$. This implies that the dust-to-metals ratio is not constant in M101, but decreases as a function of radius, equivalent to a lower fraction of metals trapped in dust at low metallicity (large radius). The dust-to-metals ratio in M101 remains at or above what would be predicted by the minimum depletion level of metals observed in the Milky Way. Our current knowledge of metallicity-dependent CO-to-H$_2$ conversion factor suggests that variations in the conversion factor cannot be responsible for the dust-to-metals ratio trends we observe. This change of dust-to-metals ratio is significantly correlated with molecular hydrogen fraction, which suggests that the accretion of gas phase metals onto existing dust grains could be a mechanism contributing to a variable dust-to-metals ratio.


title: Residual orbit of BC bump for CSNS/RCS

abstract: In a beam cycle, there is a chicane bump (BC bump) in the injection region of China Spallation Neutron Source (CSNS). It is a core part of the injection system and the important guarantee that the Linac beam injecting into the rapid cycling synchrotron (RCS). During the beam commissioning, the residual orbit of BC bump, i.e. the closedorbit distortion (COD) due to the imperfect BC bump, was an important problem which can affect the orbit correction and beam loss control of the RCS. In this paper, a method to measure the residual orbit of BC bump was studied and given. By using the three auxiliary powers of BC magnets, a method to correct the residual orbit of BC bump was developed. The numerical simulation results showed that the measurement and correction methods had worked well and they can be applied to the beam commissioning of CSNS/RCS. The measurement results in the machine study showed that the residual orbit of BC bump is not large and can be corrected by the three auxiliary powers of BC magnets.


title: Deciding the status of controversial phonemes using frequency distributions; an application to semiconsonants in Spanish

abstract: Exploiting the fact that natural languages are complex systems, the present exploratory article proposes a direct method based on frequency distributions that may be useful when making a decision on the status of problematic phonemes, an open problem in linguistics. The main notion is that natural languages, which can be considered from a complex outlook as information processing machines, and which somehow manage to set appropriate levels of redundancy, already "made the choice" whether a linguistic unit is a phoneme or not, and this would be reflected in a greater smoothness in a frequency versus rank graph. For the particular case we chose to study, we conclude that it is reasonable to consider the Spanish semiconsonant /w/ as a separate phoneme from its vowel counterpart /u/, on the one hand, and possibly also the semiconsonant /j/ as a separate phoneme from its vowel counterpart /i/, on the other. As language has been so central a topic in the study of complexity, this discussion grants us, in addition, an opportunity to gain insight into emerging properties in the broader complex systems debate.


title: Three-terminal normal-superconductor junction as thermal transistor

abstract: We propose a thermal transistor based on a three-terminal normal-superconductor (NS) junction with superconductor terminal acting as the base. The emergence of heat amplification is due to the negative differential thermal conductance (NDTC) effect for the NS diode in which the normal side maintains a higher temperature. The temperature dependent superconducting energy gap is responsible for the NDTC. By controlling quantum dot levels and their coupling strengths to the terminals, a huge heat amplification factor can be achieved. The setup offers an alternative tuning scheme of heat amplification factor and may find use in cryogenic applications.


title: On Deep Neural Networks for Detecting Heart Disease

abstract: Heart disease is the leading cause of death, and experts estimate that approximately half of all heart attacks and strokes occur in people who have not been flagged as "at risk." Thus, there is an urgent need to improve the accuracy of heart disease diagnosis. To this end, we investigate the potential of using data analysis, and in particular the design and use of deep neural networks (DNNs) for detecting heart disease based on routine clinical data. Our main contribution is the design, evaluation, and optimization of DNN architectures of increasing depth for heart disease diagnosis. This work led to the discovery of a novel five layer DNN architecture - named Heart Evaluation for Algorithmic Risk-reduction and Optimization Five (HEARO-5) -- that yields best prediction accuracy. HEARO-5's design employs regularization optimization and automatically deals with missing data and/or data outliers. To evaluate and tune the architectures we use k-way cross-validation as well as Matthews correlation coefficient (MCC) to measure the quality of our classifications. The study is performed on the publicly available Cleveland dataset of medical information, and we are making our developments open source, to further facilitate openness and research on the use of DNNs in medicine. The HEARO-5 architecture, yielding 99% accuracy and 0.98 MCC, significantly outperforms currently published research in the area.


title: Statistical Neurodynamics of Deep Networks: Geometry of Signal Spaces

abstract: Statistical neurodynamics studies macroscopic behaviors of randomly connected neural networks. We consider a deep layered feedforward network where input signals are processed layer by layer. The manifold of input signals is embedded in a higher dimensional manifold of the next layer as a curved submanifold, provided the number of neurons is larger than that of inputs. We show geometrical features of the embedded manifold, proving that the manifold enlarges or shrinks locally isotropically so that it is always embedded conformally. We study the curvature of the embedded manifold. The scalar curvature converges to a constant or diverges to infinity slowly. The distance between two signals also changes, converging eventually to a stable fixed value, provided both the number of neurons in a layer and the number of layers tend to infinity. This causes a problem, since when we consider a curve in the input space, it is mapped as a continuous curve of fractal nature, but our theory contradictorily suggests that the curve eventually converges to a discrete set of equally spaced points. In reality, the numbers of neurons and layers are finite and thus, it is expected that the finite size effect causes the discrepancies between our theory and reality. We need to further study the discrepancies to understand their implications on information processing.


title: Comparison of hard X-ray spectra of luminous radio galaxies and their radio-quiet counterparts

abstract: We study the differences in X-ray properties of luminous radio galaxies and their radio-quiet counterparts. In order to avoid biases associated with the dependence of X-ray properties on the black hole mass and Eddington ratio, the radio-loud and radio-quiet objects are selected to cover similar ranges of these parameters. Our studies are based on the X-ray data from the Swift/BAT catalogue. We confirm previous results that radio galaxies are on average X-ray-louder than radio-quiet AGNs, but find that their spectral slopes are very similar. This suggests that in radio-loud and radio-quiet AGNs the hard X-rays are produced in the same region and by the same mechanism. We argue that this region can be associated with the hot, geometrically thick, central portion of the accretion flow, where production of hard X-rays is likely to be dominated by Comptonization of the optical/UV radiation of the truncated 'cold' accretion disc by hot electrons. The larger X-ray luminosities in radio-loud AGNs may result from larger radiative efficiencies of the innermost portions of the accretion flows around faster rotating black holes.


title: Renormalization group analysis of phase transitions in the two dimensional Majorana-Hubbard model

abstract: A lattice of interacting Majorana modes can occur in a superconducting film on a topological insulator in a magnetic field. The phase diagram as a function of interaction strength for the square lattice was analyzed recently using a combination of mean field theory and renormalization group methods, and was found to include second order phase transitions. One of these corresponds to spontaneous breaking of an emergent U(1) symmetry, for attractive interactions. Despite the fact that the U(1) symmetry is not exact, this transition was claimed to be in a supersymmetric universality class when time reversal symmetry is present and in the conventional XY universality class otherwise. Another second order transition was predicted for repulsive interactions with time reversal symmetry to be in the same universality class as the transition occurring in the Gross-Neveu model, despite the fact that the U(1) symmetry is not exact in the Majorana model. We analyze these phase transitions using a modified $\epsilon$-expansion, confirming the previous conclusions.


title: Fisher Information and Natural Gradient Learning of Random Deep Networks

abstract: A deep neural network is a hierarchical nonlinear model transforming input signals to output signals. Its input-output relation is considered to be stochastic, being described for a given input by a parameterized conditional probability distribution of outputs. The space of parameters consisting of weights and biases is a Riemannian manifold, where the metric is defined by the Fisher information matrix. The natural gradient method uses the steepest descent direction in a Riemannian manifold, so it is effective in learning, avoiding plateaus. It requires inversion of the Fisher information matrix, however, which is practically impossible when the matrix has a huge number of dimensions. Many methods for approximating the natural gradient have therefore been introduced. The present paper uses statistical neurodynamical method to reveal the properties of the Fisher information matrix in a net of random connections under the mean field approximation. We prove that the Fisher information matrix is unit-wise block diagonal supplemented by small order terms of off-block-diagonal elements, which provides a justification for the quasi-diagonal natural gradient method by Y. Ollivier. A unitwise block-diagonal Fisher metrix reduces to the tensor product of the Fisher information matrices of single units. We further prove that the Fisher information matrix of a single unit has a simple reduced form, a sum of a diagonal matrix and a rank 2 matrix of weight-bias correlations. We obtain the inverse of Fisher information explicitly. We then have an explicit form of the natural gradient, without relying on the numerical matrix inversion, which drastically speeds up stochastic gradient learning.


title: Optimizing the MIMO Cellular Downlink: Multiplexing, Diversity, or Interference Nulling?

abstract: A base-station (BS) equipped with multiple antennas can use its spatial dimensions in three different ways: (1) to serve multiple users, thereby achieving a multiplexing gain, (2) to provide spatial diversity in order to improve user rates and (3) to null interference in neighboring cells. This paper answers the following question: What is the optimal balance between these three competing benefits? We answer this question in the context of the downlink of a cellular network, where multi-antenna BSs serve multiple single-antenna users using zero-forcing beamforming with equal power assignment, while nulling interference at a subset of out-of-cell users. Any remaining spatial dimensions provide transmit diversity for the scheduled users. Utilizing tools from stochastic geometry, we show that, surprisingly, to maximize the per-BS ergodic sum rate, with an optimal allocation of spatial resources, interference nulling does not provide a tangible benefit. The strategy of avoiding inter-cell interference nulling, reserving some fraction of spatial resources for multiplexing and using the rest to provide diversity, is already close-to-optimal in terms of the sum-rate. However, interference nulling does bring significant benefit to cell-edge users, particularly when adopting a range-adaptive nulling strategy where the size of the cooperating BS cluster is increased for cell-edge users.


title: Selective versions of $\theta$-density

abstract: In [8] the authors initiate the study of selective versions of the notion of $\theta$-separability in non-regular spaces. In this paper we continue this investigation by establishing connections between the familiar cardinal numbers arising in the set theory of the real line, and game-theoretic assertions regarding $\theta$-separability.


title: Optical TEMPEST

abstract: Research on optical TEMPEST has moved forward since 2002 when the first pair of papers on the subject emerged independently and from widely separated locations in the world within a week of each other. Since that time, vulnerabilities have evolved along with systems, and several new threat vectors have consequently appeared. Although the supply chain ecosystem of Ethernet has reduced the vulnerability of billions of devices through use of standardised PHY solutions, other recent trends including the Internet of Things (IoT) in both industrial settings and the general population, High Frequency Trading (HFT) in the financial sector, the European General Data Protection Regulation (GDPR), and inexpensive drones have made it relevant again for consideration in the design of new products for privacy. One of the general principles of security is that vulnerabilities, once fixed, sometimes do not stay that way.


title: Intersections, circuits, and colorability of line segments

abstract: We derive sharp upper and lower bounds on the number of intersection points and closed regions that can occur in sets of line segments with certain structure, in terms of the number of segments. We consider sets of segments whose underlying planar graphs are Halin graphs, cactus graphs, maximal planar graphs, and triangle-free planar graphs, as well as randomly produced segment sets. We also apply these results to a variant of the Erd\H{o}s-Faber-Lov\'asz (EFL) Conjecture stating that the intersection points of $m$ segments can be colored with $m$ colors so that no segment contains points with the same color. We investigate an optimization problem related to the EFL Conjecture for line segments, determine its complexity, and provide some computational approaches.


title: Hamiltonian engineering for adiabatic quantum computation: Lessons from shortcuts to adiabaticity

abstract: We discuss applications of shortcuts to adiabaticity (STA) to adiabatic quantum computation. After reviewing the fundamental properties and the present status of STA from the author's personal point of view, we apply the method to the adiabatic algorithm of the Grover's problem. We discuss two possible implementations of STA for the adiabatic quantum computations: the method of quantum adiabatic brachistochrone and the Lewis-Riesenfeld invariant-based inverse engineering.


title: A Search for Extended Radio Sources in 1.3 sr of the VLA Sky Survey (VLASS)

abstract: We report preliminary results of a visual inspection of ~4300 deg$^2$ covered by 4414 images of the 3-GHz VLA Sky Survey (VLASS, epoch 1.1) in search of extended radio structures. Over 7600 positions were registered, and for a subset of 270 of the most promising candidates their host objects were searched in optical and infrared surveys, as well as in other public databases. Among these 270, we found 9 new giant radio galaxies with a projected radio extent exceeding 1 Mpc, of which two are likely quasars at a redshift of z$\sim$2, and another one had hitherto been misidentified with a low-redshift host. Apart from various sources with unusual radio morphologies, we found numerous remnant-type double radio sources with very diffuse lobes in VLASS, both with and without radio cores. We confirm that the high angular resolution and observing frequency of VLASS is well suited for the optical identification of radio sources, and despite its limited sensitivity to extended emission, allows to reveal large radio galaxies based on morphological features unrecognizable in lower-resolution surveys.


title: Propagation Estimates for Two-cluster Scattering Channels of N-body Schr\"odinger Operators

abstract: In this paper we prove propagation estimates for two-cluster scattering channels of N-body Schr\"odinger operators. These estimates are based on the estimate similar to Mourre's commutator estimate and the method of Skibsted. We also obtain propagation estimates with better indices using projections onto almost invariant subspaces close to two-cluster scattering channels. As an application of these estimates we obtain the resolvent estimate for two-cluster scattering channels and microlocal propagation estimates in three-body problems without projections. Our method clearly illustrates evolution of the solutions of the Schr\"odinger equation.


title: Quantum metrology out of equilibrium

abstract: We address open quantum systems out-of-equilibrium as effective quantum probes for the characterisation of their environment. We discuss estimation schemes for parameters driving a de-phasing evolution of the probe and then focus on qubits, establishing a relationship between the quantum Fisher information and the residual coherence of the probe. Finally, we apply our results to the characterisation of the ohmicity parameter of a bosonic environment.


title: Efficient sparse semismooth Newton methods for the clustered lasso problem

abstract: We focus on solving the clustered lasso problem, which is a least squares problem with the $\ell_1$-type penalties imposed on both the coefficients and their pairwise differences to learn the group structure of the regression parameters. Here we first reformulate the clustered lasso regularizer as a weighted ordered-lasso regularizer, which is essential in reducing the computational cost from $O(n^2)$ to $O(n\log (n))$. We then propose an inexact semismooth Newton augmented Lagrangian ({\sc Ssnal}) algorithm to solve the clustered lasso problem or its dual via this equivalent formulation, depending on whether the sample size is larger than the dimension of the features. An essential component of the {\sc Ssnal} algorithm is the computation of the generalized Jacobian of the proximal mapping of the clustered lasso regularizer. Based on the new formulation, we derive an efficient procedure for its computation. Comprehensive results on the global convergence and local linear convergence of the {\sc Ssnal} algorithm are established. For the purpose of exposition and comparison, we also summarize/design several first-order methods that can be used to solve the problem under consideration, but with the key improvement from the new formulation of the clustered lasso regularizer. As a demonstration of the applicability of our algorithms, numerical experiments on the clustered lasso problem are performed. The experiments show that the {\sc Ssnal} algorithm substantially outperforms the best alternative algorithm for the clustered lasso problem.


title: Can 3D Pose be Learned from 2D Projections Alone?

abstract: 3D pose estimation from a single image is a challenging task in computer vision. We present a weakly supervised approach to estimate 3D pose points, given only 2D pose landmarks. Our method does not require correspondences between 2D and 3D points to build explicit 3D priors. We utilize an adversarial framework to impose a prior on the 3D structure, learned solely from their random 2D projections. Given a set of 2D pose landmarks, the generator network hypothesizes their depths to obtain a 3D skeleton. We propose a novel Random Projection layer, which randomly projects the generated 3D skeleton and sends the resulting 2D pose to the discriminator. The discriminator improves by discriminating between the generated poses and pose samples from a real distribution of 2D poses. Training does not require correspondence between the 2D inputs to either the generator or the discriminator. We apply our approach to the task of 3D human pose estimation. Results on Human3.6M dataset demonstrates that our approach outperforms many previous supervised and weakly supervised approaches.